Not so pseudo: the evolutionary history of protein phosphatase 1 regulatory subunit 2 and related pseudogenes
© Korrodi-Gregório et al.; licensee BioMed Central Ltd. 2013
Received: 15 May 2013
Accepted: 29 October 2013
Published: 6 November 2013
Pseudogenes are traditionally considered “dead” genes, therefore lacking biological functions. This view has however been challenged during the last decade. This is the case of the Protein phosphatase 1 regulatory subunit 2 (PPP1R2) or inhibitor-2 gene family, for which several incomplete copies exist scattered throughout the genome.
In this study, the pseudogenization process of PPP1R2 was analyzed. Ten PPP1R2-related pseudogenes (PPP1R2P1-P10), highly similar to PPP1R2, were retrieved from the human genome assembly present in the databases. The phylogenetic analysis of mammalian PPP1R2 and related pseudogenes suggested that PPP1R2P7 and PPP1R2P9 retroposons appeared before the great mammalian radiation, while the remaining pseudogenes are primate-specific and retroposed at different times during Primate evolution. Although considered inactive, four of these pseudogenes seem to be transcribed and possibly possess biological functions. Given the role of PPP1R2 in sperm motility, the presence of these proteins was assessed in human sperm, and two PPP1R2-related proteins were detected, PPP1R2P3 and PPP1R2P9. Signatures of negative and positive selection were also detected in PPP1R2P9, further suggesting a role as a functional protein.
The results show that contrary to initial observations PPP1R2-related pseudogenes are not simple bystanders of the evolutionary process but may rather be at the origin of genes with novel functions.
KeywordsPPP1 PPP1R2 Retroposons Pseudogenization Evolution Vertebrates
In the past, pseudogenes were generally regarded as functionally inert, due to the presence of several disabling features that prevent their expression (e.g. premature stop codons, frameshift mutations, no promoter regions, etc.), and therefore their evolution has been considered to be neutral . However, this view has been challenged by new evidences, which demonstrate that certain pseudogenes are functionally active [1, 2]. The GENCODE, a sub-project of the ENCODE (ENCyclopedia Of DNA Elements from the National Human Genome Research Institute, NHGRI), has estimated the number of pseudogenes in the human genome to be near 14,000 . From these, ~6% were identified has potentially transcribed by computational models and almost half of them validated by RT-PCR-Seq techniques . Indeed, pseudogenes can be functional at the DNA, RNA or protein levels and have a function related or independent of the parental gene . At the DNA level, pseudogenes can regulate other genes by pseudogene insertion in the non-coding or coding region of the target gene and regulate the parental counterpart gene by gene conversion, homologous recombination and through regulatory sequences . Concerning the RNA level, pseudogene RNAs can compete with the parental mRNA for miRNAs, RNA binding proteins and/or translational machinery binding, as well as, functioning as siRNAs and thereby inhibiting the parental gene expression. Pseudogenes can also function in unrelated genes as long non-coding RNAs, by encoding miRNA precursors or even compete for miRNAs . At the protein level, pseudogenic proteins can have the same activity of the parental protein but function in different tissues, subcellular localization and/or pathophysiological conditions [5–11]. Pseudogenic proteins with altered functions might also affect the activity of the parental ones . If a pseudogene mRNA is translated to a functional pseudogenic protein, this gene is often called a retrogene . Pseudogenes can also produce truncated proteins that can function as antigenic peptides in the surface of the cells to stimulate the immune system against the malignant cells . Pseudogenes have already been associated with several pathological conditions such as cancer , diabetes  and neurodegenerative diseases .
One promising model to understand the functional relevance of pseudogenization is the protein phosphatase 1 regulatory subunit 2 (PPP1R2). This protein, also known as inhibitor-2 (I2), was one of the first regulatory subunits identified as an inhibitor and binding partner of the Ser/Thr phosphoprotein phosphatase 1 (PPP1). PPP1R2 forms a stable complex with PPP1 catalytic subunit (PPP1C) blocking the active site and inhibiting it potently, being the reactivation triggered by phosphorylation [16–21]. The PPP1C/PPP1R2 complex has been implied in several processes such as cardiac function [22–24], mitosis and meiosis [25–30], tubulin acetylation and neuronal cell survival . Also, it has been previously shown that a PPP1CC2/PPP1R2-like complex is important in the acquisition of sperm motility [32, 33].
The PPP1R2 gene is conserved throughout all eukaryotes, from yeast to humans, with homologues found even in plants [34, 35]. In the human genome, as observed for other ancient PPP1 inhibitors such as PPP1R8 (NIPP1) and PPP1R11 (I3), several sequences have been identified that are highly similar to PPP1R2 . For PPP1R2, nine loci were found that present hallmark features of processed pseudogenes. These related sequences were collectively named PPP1R2 pseudogenes and were numbered from 1 to 9 (PPP1R2P1-P9) . These pseudogenes are found scattered in the genome due to random retrotransposition phenomena that consist on the reverse transcription of cellular RNAs and random insertion into the nuclear genome [36, 37]. Past studies identified four PPP1R2 pseudogenes at the messenger RNA level using high throughput techniques. PPP1R2P1 and PPP1R2P2 were discovered in human [38, 39], PPP1R2P3 in human and crab-eating macaque (Macaca fascicularis) and PPP1R2P9 (also called I4) was found in human and mouse (Mus musculus) [40–43].
In this work we performed an exhaustive search for PPP1R2 pseudogenes in publicly available mammalian genome databases in order to infer their evolutionary history. In the collected pseudogenes, an assay for detection of the proteins was conducted. Our results show that evolution and pseudogenization of PPP1R2 gene may be correlated with the formation of new genes and the gain of new specific functions.
Results and discussion
A total of 119 sequences were retrieved from the NCBI and Ensembl databases by blasting against the human PPP1R2 mRNA sequence. Ten pseudogenes were obtained from human sequences, increasing by one the previous number reported in the literature . All pseudogenes obtained are intronless and with a truncated 5’UTR meaning that are processed pseudogenes. The parental human PPP1R2 CDS (618 bp) covers 17% of the entire mRNA (3475 bp); even the pseudogenes with the lowest coverage contain the parental CDS, with the exception of PPP1R2P7 that only comprises part of the 3’UTR.
In order to increase the reliability of the alignment for the phylogenetic reconstruction, we selected sequences with >85% coverage and >60% similarity with the human PPP1R2 CDS. By doing this, 81 sequences were included in the tree that represented all the pseudogenes with the exception of PPP1R2P7 (Additional file 1: Table S1). The unused sequences encompassed pseudogenic fragments and sequences where identity with PPP1R2 was detected mostly outside the CDS (e.g. PPP1R2P7) or presented truncated CDS (e.g. some PPP1R2P8 and PPP1R2P9).
Evidences for functionality of PPP1R2-related pseudogenes
Features such as the existence of transcriptional related data, presence of regulatory elements, mRNA stability (e.g. UTRs, polyA signals), translation initiator sequence and complete ORFs (no truncations or disabling mutations) are indicators of the putative functionality of genes. A search for such features was conducted in order to verify the potential functionality of the PPP1R2 pseudogenes.
The Gene Expression Omnibus (GEO, NCBI) and Gene Expression Atlas (GXA, Ensembl) public repositories contain expression data for PPP1R2P1. The presence of promoters, enhancers and other regulatory elements could be an explanation for PPP1R2P1 transcriptional related data (174 GEO and 2 GXA), although basal transcription should not be set aside. Concerning the mRNA stability, only part of the 5’UTR (238 bp), due to the low processivity of the reverse transcriptase, and part of the 3’UTR (506 bp) are present. Therefore, the stability might be compromised although a polyA signal (ATTAAA) is present near the 3’UTR terminus (position 1361, Figure 3). Regarding the translation, the Kozak sequence, important for translation initiation, is present in the parental gene and is conserved in PPP1R2P1. Altogether, these results suggest that at least in humans, PPP1R2P1 is expressed and might be functionally relevant. Although we cannot set aside the low quality of some of the assembled genomes, in other primates the ORF of PPP1R2P1 has frameshift disruptions that introduce premature stop codons, indicating that in these species might not produce a putative functional protein, or if so the protein might be truncated.
The sequence of PPP1R2P3 is complete, without any frameshifts or element repeats disruptions (Figure 3). The sequence was truncated at the 5’UTR, as expected due to the low processivity of the reverse transcriptase, and in the 3’UTR it lost two of the four polyA signals that may lead to a short ~1500-1600 nt message. We have previously found, by a yeast two hybrid screening of human testis cDNA, using as bait PPP1CC1, one clone assigned to PPP1R2P3 [43, 45]. A search for PPP1R2P3 ESTs in databases revealed that this is one of the most represented PPP1R2 pseudogenes (72 GEO, 313 GXA), being highly detected in testis (14 ESTs in Unigene). Together with our previous data, this strongly suggests that this pseudogene is transcribed. Two independent reports using mass spectrometry have also assigned peptides to PPP1R2P3 [46, 47]. However, these peptides share the sequence with both parental gene and PPP1R2P3, being most probably misassigned. Nonetheless, we have shown recently by mass spectrometry the presence of PPP1R2P3 in human sperm samples .
The PPP1R2P9 sequences retrieved have not been disrupted, at least in primates (Figure 6). However, the 5’UTR of the parental gene is absent and the 3’UTR is truncated (671 bp in humans). At the 3’UTR there is a single polyA signal at nucleotide position 1088, according to the human sequence, which suggests that a shorter message is produced. Sequence repeats, deletions, unknown and known sequence insertions were only found in the PPP1R2-like sequences (Figure 6). The only exceptions are in mouse and rat where the 3’UTR was deleted in the parental PPP1R2P9 (Figure 6).
This pseudogene is the one with more transcriptional related data (1086 GEO, 128 GXA) and has many ESTs in testis (9 ESTs in Unigene) like PPP1R2P3. PPP1R2P9 was originally found in cDNA libraries of human germ cell tumors, binding to PPP1C directly and in heat stable extracts inhibits this phosphatase potently with an IC50 of 0.2nM . Also, we have recently identified PPP1R2P9 as an interacting partner of PPP1CA by yeast-two hybrid in human brain . This suggests that silent regulatory areas are present in the region were PPP1R2P9 was retroposed and that during the evolution PPP1R2P9 might have retained or gained the capacity to be transcribed. In spite of this, there is no data suggesting the translation of PPP1R2P9. Considering the ORF, all species show a continuous ORF with no or small truncations at the C-terminus (e.g. in mouse and rat), with the exception of pig where no protein translation was obtained from the ORF.
Evidences of non-coding nature of PPP1R2-related pseudogenes
Considering the other pseudogenes sequences, many insertions in PPP1R2P8 lead to a completely disrupted ORF and missense mutations in PPP1R2P4/P10 and PPP1R2P6 lead to premature stop codons (Figures 4 and 5). Also, since these pseudogenes showed low coverage to the parental gene, most of the 3’UTR is missing and so, regulatory elements such as polyadenylation signals that are important for the transcript cleavage and stability are absent. This indicates that no transcription or translation should be expected from these pseudogenes, which corroborates with the fact that no expression was found in ESTs and high-throughput databases with the exception of PPP1R2P4. Considering the pseudogenes with the highest coverage in relation to the parental gene, PPP1R2P2 and PPP1R2P5, no ORF disruptions were found but many missense mutations were found in all species analyzed that lead to premature stop codons (Figure 4). All four polyadenylation signals present in the parental PPP1R2 mRNA are conserved in PPP1R2P2. Although protein expression is unlikely, PPP1R2P2 message was found by qPCR in human testis but not in peripheral blood leukocytes . Also, two experiments from ArrayExpress, report the up/down regulation of this pseudogene in prostate adenocarcinoma and in a prostate transcriptomic study performed in a Caucasian population . These results might be artifacts or could be due to other PPP1R2 pseudogenes/parental gene since this pseudogene is located in chromosome 21 that has low density (~232 genes, only surpassed by the Y chromosome with 130 genes), and as expected, the processed pseudogene density is also low, 34 , making the transcription highly unlikely.
Detection of PPP1R2-related proteins
PPP1R2P9 presence in human sperm
Protein size( aa)
Range ( start - end )
Number of spectra
This analysis was based on the fact that the molecular weight of these PPP1R2-related proteins should be similar to the parental one (PPP1R2-23.0kDa), being therefore present in the same region where the band was extracted to mass spectrometry analysis. The antibody used to immunoprecipitate PPP1R2-related proteins was raised against a peptide containing amino acid residues 134–147 from the mouse PPP1R2 sequence (Figure 7). This antibody was used previously to detect PPP1R2 [48, 58]. In the 14-residue region, PPP1R2P1 and PPP1R2P9 have two and three substitutions respectively, when comparing to PPP1R2 sequence. We predicted that using this antibody, we were also able to detect the other PPP1R2-related proteins. Mass spectrometry data identified 23 MSMS (tandem mass spectrometry) spectra corresponding to 8 different peptides matching unequivocally to PPP1R2P9 (Figure 7 and Table 1) and 3 MSMS spectra corresponding to one peptide matching unequivocally to PPP1R2P3 .
The sequence coverage obtained for PPP1R2P9 was 36.5% and the mascot score levels were 623.41 (in addition, spectra were manually evaluated). This is the first time that PPP1R2P9 protein is detected, being clearly recovered from human ejaculated sperm. Additionally, these results also indicate that native PPP1R2-related proteins are indeed heat stable and migrate at the same position as the parental PPP1R2. Lastly, no peptides were recovered for PPP1R2P1 using this method, which might suggest the absence of the protein, at least from the sperm cells.
Signatures of selection
Pseudogenes have been regarded as being derived from functional-encoding genomic DNA sequences that have accumulated disabling mutations (frameshifts and premature stop codons) that make them non-coding protein genes. This lack of function predicts that pseudogenes are not under selective pressures and thus evolve neutrally (reviewed in ). Nevertheless, this view keeps being challenged by the accumulation of examples of transcribed pseudogenes with several acknowledged functions (e.g. regulation of the expression of paralogous genes through the generation of small-interfering RNAs) . Signatures of selection, in addition to sequence conservation, have been considered as obvious indicators of the functional importance of pseudogenes .
Here, by using six ML methods, signatures of both positive and negative selection were detected in the PPP1R2P9 pseudogene, as well as in the parental gene PPP1R2 (Additional file 2: Table S2). Signatures of negative selection were far more evident than those of positive selection, for both genes. Four methods, REL, FEL, SLAC and FUBAR, showed sites negatively selected, with most being detected by more than one method. Signatures of positive selection were principally detected by FEL and MEME methods. The codons 92 and 120, for PPP1R2, and the codons 6, 208 and 211, for PPP1R2P9, were detected by at least two separate methods. No detection was obtained for the PAML method. It is known that sperm-expressed genes present in chromosome X tend to be positively selected when compared with X-linked non-sperm genes and with sperm-expressed autosomal genes [61, 62]. This evolutionary pressure is due to their hemizigous expression in males that will favor advantageous mutations and remove any deleterious one. PPP1R2P9 is not evolving neutrally and may thus be expressed, further supporting a functional role for this pseudogene.
Retropositions from the PPP1R2 gene are ancient, prior to the great radiation of the mammals, as supported by the presence of PPP1R2P9 and PPP1R2P7 in the different groups of mammals. All the other pseudogenes found in humans are primate-specific and were retroposed at different times during the evolution of this group. For instance, PPP1R2P3 exists only in the members of the Hominoidea family, whereas PPP1R2P8, the most distinct, is present in all groups and was retroposed ~42.6-65.2 Mya. This reveals that retropositions have occurred in waves and in a unique way similar to the Alu repeats explosion that occurred ~40-50 Mya, after the divergence of simian ancestors from the prosimians (lemurs and lorises). The recent pseudogene duplication in humans, PPP1R2P4, and in chimpanzee, PPP1R2P1, suggests that evolution of pseudogenes is still an active process.
As suggested by the presence of an uninterrupted ORF, ESTs and polyA signals, PPP1R2P9 (along with PPP1R2P1 and PPP1R2P3) appears to be transcribed. Moreover, the finding of positive and negative selection signatures suggests that it could be functionally relevant. Indeed, we confirmed that two PPP1R2-related proteins are translated in human sperm (PPP1R2P3 and PPP1R2P9), and are heat stable in their native form . The importance of these PPP1R2-related proteins in physiological conditions, such as spermatogenesis and sperm physiology, should be assessed in future studies. Besides this, PPP1R2P1, PPP1R2P3 and PPP1R2P9 were found to be associated with pathological conditions [15, 38, 40, 42, 63]. Thus, assessing their ratios may be considered as a diagnostic tool in the future.
Furthermore, it has been shown that pseudogenes can regulate their parental counterparts at the RNA level either by siRNA or by competition for positive and negative stabilizing factors and miRNAs . Although PPP1R2P2, PPP1R2P4 and PPP1R2P10 translation is very unlikely, their expression is documented and so, it is feasible these pseudogenes could regulate the parental PPP1R2 message levels and therefore its function.
These observations indicate that PPP1R2 pseudogenes have possible biological functions rather than acting as non-functional relics as initially believed. Their evolution process might be in part related with the formation of new genes and the gain of new specific functions. Therefore, their designation as pseudogenes should be reevaluated.
The human PPP1R2 mRNA sequence (GenBank accession number NM_006241.4) was used to detect orthologs and pseudogene-related sequences by performing a BLAST search on GenBank, from National Center for Biotechnology Information (NCBI, http://BLAST.ncbi.nlm.nih.gov/) and Ensembl (http://www.ensembl.org/Multi/blastview) databases against all available mammalian reference genomic sequences. Only sequences with more than 60% of sequence similarity and with query coverage of more than 35% were recovered. Genomic sequences flanking the retrieved sequences were also manually inspected for missing parts, especially at the 3’UTR.
Evolutionary tree reconstruction and divergence times
The retrieved sequences (Additional file 1: Table S1) were visually inspected and aligned using ClustalW implemented in BioEdit 188.8.131.52 . For phylogenetic reconstruction, and to improve accuracy, only sequences encompassing >85% coverage of the human PPP1R2 CDS (nucleotide positions 377–994 of the mRNA sequence) and with >60% of sequence similarity were included in the alignment. In order to determine the phylogenetic relationships between the PPP1R2 gene and related pseudogenes, the best-fit model of nucleotide substitution was first assessed using the program jModelTest v0.1.1  under the Akaike Information Criterion (AIC). A maximum likelihood (ML) phylogeny was inferred using the software GARLI v1.0  by indicating the best nucleotide substitution model. No starting topology was defined and the program was set to run until no significant topology improvement (as defined by the default settings) was found after 1000000 generations. Five independent runs were performed to check the consistency of the estimates. The support of each node was assessed using 1000 bootstrap replicates. For each bootstrap replicate, the number of generations was set at 100000, above the generation where the last topological improvements were found for each of the five independent replicates. A 50% majority-rule consensus tree of the 1000 bootstrap replicates was created using PAUP* . The support values at each node of the consensus tree were added to the best tree found by GARLI.
Pseudogene classification and conserved linkage
Sequences obtained from the BLAST queries were analyzed in terms of presence of intronic regions, polyA traits (PolyApred, http://www.imtech.res.in/raghava/polyapred/), truncation of the 5’UTR and chromosomal location. Chromosomal locations were obtained from the GenBank database (Additional file 1: Table S1). Pseudogenes located in the same chromosome and nearby and/or with intronic regions were classified as duplicated pseudogenes. Pseudogenes that were located in different chromosomes and had polyA traits, truncation of the 5’UTR and no introns were classified as processed pseudogenes. Furthermore, genes flanking each human PPP1R2 pseudogene and conserved among mammals were selected. Conserved linkage, meaning conservation of synteny and also conservation of the gene order, was then searched for in order to provide insights regarding their orthology.
Distance to closest and repeated regions
The distance of each pseudogene to the closest neighboring gene, not taking into account the presence of nearby pseudogenes, was calculated. Repeated sequences were detected by submitting each pseudogene sequence to the program RepeatMasker from Institute for Systems Biology, Seattle, Washington, USA (http://www.repeatmasker.org/).
Signatures of natural selection
Coding sequences evolving neutrally present a ratio (ω) of non-synonymous (dN) over synonymous substitutions (dS) that do not significantly deviate from one. An excess of non-synonymous substitutions over synonymous substitutions (dN > dS) might indicate positive selection, suggesting that the replacement might be advantageous, while negative selection results from the scarcity of non-synonymous substitutions (dN < dS), indicating that a particular mutation most likely is deleterious and is being removed from the gene pool. Pseudogenes are considered to evolve neutrally (reviewed in ).
Maximum-likelihood codon-based tests were used to test for statistically significant signatures of selection in PPP1R2 and related-pseudogenes. Nevertheless, only PPP1R2P9 sequences were analyzed since at least 10 sequences are required to robustly detect signatures of selection . Signatures of positive and negative selection were searched for in Datamonkey webserver (http://www.datamonkey.org) that uses the HyPhy package . The best-fitting nucleotide model (GTR + G) was determined using the automated tool provided by Datamonkey. Five models were used: single likelihood ancestor counting (SLAC), fixed-effect likelihood (FEL), random effect likelihood (REL), fast unbiased bayesian approximation (FUBAR) and mixed effects model of evolution (MEME). SLAC is based on the reconstruction of the ancestral sequences and the counts of dS and dN at each codon position of the phylogeny. FEL estimates the ratio of dN/dS on a site-by-site basis, without assuming an a priori distribution across sites while REL fits a distribution of rates across sites and then infers the substitution rate for individual sites. FUBAR detects selection much faster than the other methods and to leverage Bayesian MCMC to robustly account for parameter estimation errors. Finally, MEME is capable of identifying instances of both episodic and pervasive positive selection at the level of an individual site. Sites with P values <0.1 for SLAC, FEL and MEME, posterior probability of >0.9 for FUBAR, and Bayes Factor >50 for REL were considered as being under selection. CODEML (PAML version 4, ) was also used to detect positive selection by comparing a null model and a model that allows positive selection (M1 vs. M2 and M7 vs. M8). The contrasting models were compared by computing twice the difference in the natural logs of the likelihoods (2ΔlnL). In the site-specific models that allow the ratio ω to vary among codons, we performed Likelihood Ratio Tests (LRTs) with 2 degrees of freedom to compare the following models (NS sites): M1 (nearly neutral evolution ω0 = 0, ω1 = 1) with M2 (neutral and positive selection: ω 0 = 0, ω 1 = 1, ω 2 > 1) and M7 (beta-distributed negative selection: 0 # ω # 1) with M8 (beta-distributed negative selection and positive selection: 0 # ω1 # 1, ω2 >1) [2, 73]. Only amino acids identified in M8 by using the Bayes Empirical Bayes (BEB) approach and with posterior probability >95% were considered as evolving under positive selection. For the initial working topology, ML trees were constructed using MEGA5  with substitution nucleotide models determined by the software: TN93 + I and partial deletion (95% cut-off) for PPP1R2P9 and K2 + G with G = 4 and partial deletion (95% cut-off) for PPP1R2.
Sperm extracts and immunoprecipitation
Since testis is one of the organs where most pseudogenes are expressed  and spermatozoa are the final product of spermatogenesis, the presence of some of the studied pseudogenes was tested in human sperm. Ejaculated sperm was collected from healthy donors by masturbation into an appropriate sterile container. Spermograms were performed by experienced technicians and only samples with normal parameters were used . Informed consents were signed allowing samples to be used for scientific purposes. The study was conducted in accordance with the guidelines of the “Helsinki Declaration”. In brief, sperm was lysed in 1 × RIPA buffer (radioimmunoprecipitation buffer, Millipore Iberica S.A.U., Madrid, Spain) supplemented with protease inhibitors (10 mM benzamidine, 1.5 μM aprotinin, 5 μM pepstatin A, 2 μM leupeptin, 1 mM PMSF), sonicated 3 × 10 sec and centrifuged at 16000 g for 20 min, at 4°C. RIPA supernatant sperm extract was immunoprecipitated using Dynabeads® Protein G (Life Technologies S.A., Madrid, Spain) and 1 μg of rabbit anti-PPP1R2 (against a mouse PPP1R2 peptide, amino acids 134–147) with standard direct immunoprecipitation procedure . Also, an independent RIPA supernatant sperm extract was prepared, boiled in a water bath for 30 min, chilled on ice for 2 min and centrifuged at 16000 g for 20 min, 4°C to obtain a heat stable extract.
For mass spectrometry analysis, the immunoprecipitate and the heat stable extract were resolved by 10% SDS-PAGE along with purified positive controls. Gels were stained with Coomassie blue colloidal (Sigma-Aldrich Química, S.A., Sintra, Portugal) using standard procedures . Bands were then excised from the gel using commercial PPP1R2 band as control and destained. An overnight digestion with trypsin (Promega, Madison, Wisconsin, USA) was performed and resulting peptides were extracted and prepared for mass spectrometry analysis using an Orbitrap Velos mass spectrometer as described elsewhere . Subsequent generated data were imported to ProteinScapeTM (Bruker Daltonik GmbH, Bremen, Germany, ) and analyzed using MASCOT (version 2.2.0, Matrix Science, London, UK, ) search algorithm. Proteins were considered to be identified if the Mascot score (ProteinScapeTM) was higher than 65.
We thank Srinivasan Vijayaraghavan (Anatomy and Cell Biology Department, Kent State University, Kent, Ohio, USA) for kindly providing the rabbit anti-PPP1R2 antibody. The authors have declared that no conflict of interest exists.
This work was supported by Programa Operacional Potencial Humano (POPH) – Quadro de Referência Estratégico Nacional (QREN) funds from the European Social Fund and Portuguese Ministério da Educação e Ciência through the Portuguese Foundation for Science and Technology-FCT–(SFRH/BD/41751/2007 to L.K.G., SFRH/BPD/73512/2010 to J.A., SFRH/BPD/43264/2008 to J.M.F. and SFRH/BPD/27021/2006 to P.J.E.). This work was also supported by FEDER (Fundo Europeu de Desenvolvimento Regional) funds through the Programa Operacional Factores de Competitividade (COMPETE program) and Portuguese national funds through the FCT projects PTDC/BIA-BEC/103158/2008 and PTDC/CVT/108490/2008. Also, the project “Genomics Applied to Genetic Resources,” co-financed by North Portugal Regional Operational Program 2007/2013 (ON.2–O Novo Norte), under the National Strategic Reference Framework, through the European Regional Development Fund, supported this work.
- Balakirev ES, Ayala FJ: Pseudogenes: are they “junk” or functional DNA?. Annu Rev Genet. 2003, 37 (1): 123-151. 10.1146/annurev.genet.37.040103.103949.PubMedView ArticleGoogle Scholar
- Zheng D, Gerstein MB: The ambiguous boundary between genes and pseudogenes: the dead rise up, or do they?. Trends Genet. 2007, 23 (5): 219-224. 10.1016/j.tig.2007.03.003.PubMedView ArticleGoogle Scholar
- Pei B, Sisu C, Frankish A, Howald C, Habegger L, Mu XJ, Harte R, Balasubramanian S, Tanzer A, Diekhans M, et al: The GENCODE pseudogene resource. Genome Biol. 2012, 13 (9): R51-10.1186/gb-2012-13-9-r51.PubMed CentralPubMedView ArticleGoogle Scholar
- Poliseno L: Pseudogenes: newly discovered players in human cancer. Sci Signal. 2012, 5 (242): re5-PubMedView ArticleGoogle Scholar
- McEntee G, Minguzzi S, O'Brien K, Ben Larbi N, Loscher C, O'Fagain C, Parle-McDermott A: The former annotated human pseudogene dihydrofolate reductase-like 1 (DHFRL1) is expressed and functional. Proc Natl Acad Sci U S A. 2011, 108 (37): 15157-15162. 10.1073/pnas.1103605108.PubMed CentralPubMedView ArticleGoogle Scholar
- McCarrey JR, Thomas K: Human testis-specific PGK gene lacks introns and possesses characteristics of a processed gene. Nature. 1987, 326 (6112): 501-505. 10.1038/326501a0.PubMedView ArticleGoogle Scholar
- Pai HV, Kommaddi RP, Chinta SJ, Mori T, Boyd MR, Ravindranath V: A frameshift mutation and alternate splicing in human brain generate a functional form of the pseudogene cytochrome P4502D7 that demethylates codeine to morphine. J Biol Chem. 2004, 279 (26): 27383-27389. 10.1074/jbc.M402337200.PubMedView ArticleGoogle Scholar
- Zhang J, Wang X, Li M, Han J, Chen B, Wang B, Dai J: NANOGP8 is a retrogene expressed in cancers. FEBS J. 2006, 273 (8): 1723-1730. 10.1111/j.1742-4658.2006.05186.x.PubMedView ArticleGoogle Scholar
- Sun C, Orozco O, Olson DL, Choi E, Garber E, Tizard R, Szak S, Sanicola M, Carulli JP: CRIPTO3, a presumed pseudogene, is expressed in cancer. Biochem Biophys Res Commun. 2008, 377 (1): 215-220. 10.1016/j.bbrc.2008.09.113.PubMedView ArticleGoogle Scholar
- Takahashi K, Mitsui K, Yamanaka S: Role of ERas in promoting tumour-like properties in mouse embryonic stem cells. Nature. 2003, 423 (6939): 541-545. 10.1038/nature01646.PubMedView ArticleGoogle Scholar
- Rohozinski J, Edwards CL, Anderson ML: Does expression of the retrogene UTP14c in the ovary pre-dispose women to ovarian cancer?. Med Hypotheses. 2012, 78 (4): 446-449. 10.1016/j.mehy.2011.12.016.PubMedView ArticleGoogle Scholar
- Zou M, Baitei EY, Alzahrani AS, Al-Mohanna F, Farid NR, Meyer B, Shi Y: Oncogenic activation of MAP kinase by BRAF pseudogene in thyroid tumors. Neoplasia. 2009, 11 (1): 57-65.PubMed CentralPubMedView ArticleGoogle Scholar
- Kaessmann H, Vinckenbosch N, Long M: RNA-based gene duplication: mechanistic and evolutionary insights. Nat Rev Genet. 2009, 10 (1): 19-31.PubMed CentralPubMedView ArticleGoogle Scholar
- Chiefari E, Iiritano S, Paonessa F, Le Pera I, Arcidiacono B, Filocamo M, Foti D, Liebhaber SA, Brunetti A: Pseudogene-mediated posttranscriptional silencing of HMGA1 can result in insulin resistance and type 2 diabetes. Nat Commun. 2010, 1: 40-PubMedView ArticleGoogle Scholar
- Costa V, Esposito R, Aprile M, Ciccodicola A: Non-coding RNA and pseudogenes in neurodegenerative diseases: “the (un)usual suspects”. Front Genet. 2012, 3: 231-PubMed CentralPubMedView ArticleGoogle Scholar
- DePaoli-Roach AA: Synergistic phosphorylation and activation of ATP-Mg-dependent phosphoprotein phosphatase by F A/GSK-3 and casein kinase II (PC0.7). J Biol Chem. 1984, 259 (19): 12144-12152.PubMedGoogle Scholar
- Aitken A, Holmes CFB, Campbell DG, Resink TJ, Cohen P, Leung CTW, Williams DH: Amino acid sequence at the site on protein phosphatase inhibitor-2, phosphorylated by glycogen synthase kinase-3. Biochim Biophys Acta Protein Struct Mol Enzymol. 1984, 790 (3): 288-291. 10.1016/0167-4838(84)90034-7.View ArticleGoogle Scholar
- Holmes CF, Kuret J, Chisholm AA, Cohen P: Identification of the sites on rabbit skeletal muscle protein phosphatase inhibitor-2 phosphorylated by casein kinase-II. Biochim Biophys Acta. 1986, 870 (3): 408-416. 10.1016/0167-4838(86)90248-7.PubMedView ArticleGoogle Scholar
- Wang QM, Guan K-L, Roach PJ, DePaoli-Roach AA: Phosphorylation and activation of the ATP-Mg-dependent protein phosphatase by the mitogen-activated protein kinase. J Biol Chem. 1995, 270 (31): 18352-18358. 10.1074/jbc.270.31.18352.PubMedView ArticleGoogle Scholar
- Puntoni F, Villamoruzzi E: Phosphorylation of the inhibitor-2 of protein phosphatase-1 by cdc2-cyclin B and GSK3. Biochem Biophys Res Commun. 1995, 207 (2): 732-739. 10.1006/bbrc.1995.1248.PubMedView ArticleGoogle Scholar
- Agarwal-Mawal A, Paudel HK: Neuronal Cdc2-like protein kinase (Cdk5/p25) is associated with protein phosphatase 1 and phosphorylates inhibitor-2. J Biol Chem. 2001, 276 (26): 23712-23718. 10.1074/jbc.M010002200.PubMedView ArticleGoogle Scholar
- Kirchhefer U, Baba HA, Bokník P, Breeden KM, Mavila N, Brüchert N, Justus I, Matus M, Schmitz W, DePaoli-Roach AA, et al: Enhanced cardiac function in mice overexpressing protein phosphatase inhibitor-2. Cardiovasc Res. 2005, 68 (1): 98-108. 10.1016/j.cardiores.2005.05.019.PubMedView ArticleGoogle Scholar
- Yamada M, Ikeda Y, Yano M, Yoshimura K, Nishino S, Aoyama H, Wang L, Aoki H, Matsuzaki M: Inhibition of protein phosphatase 1 by inhibitor-2 gene delivery ameliorates heart failure progression in genetic cardiomyopathy. FASEB J. 2006, 20 (8): 1197-1199. 10.1096/fj.05-5299fje.PubMedView ArticleGoogle Scholar
- Brüchert N, Mavila N, Boknik P, Baba HA, Fabritz L, Gergs U, Kirchhefer U, Kirchhof P, Matus M, Schmitz W, et al: Inhibitor-2 prevents protein phosphatase 1-induced cardiac hypertrophy and mortality. Am J Physiol Heart Circ Physiol. 2008, 295 (4): H1539-H1546. 10.1152/ajpheart.00515.2008.PubMed CentralPubMedView ArticleGoogle Scholar
- Eto M, Bock R, Brautigan DL, Linden DJ: Cerebellar long-term synaptic depression requires PKC-mediated activation of CPI-17, a myosin/moesin phosphatase inhibitor. Neuron. 2002, 36 (6): 1145-1158. 10.1016/S0896-6273(02)01107-8.PubMedView ArticleGoogle Scholar
- Leach C, Shenolikar S, Brautigan DL: Phosphorylation of phosphatase inhibitor-2 at centrosomes during mitosis. J Biol Chem. 2003, 278 (28): 26015-26020. 10.1074/jbc.M300782200.PubMedView ArticleGoogle Scholar
- Satinover DL, Leach CA, Stukenberg PT, Brautigan DL: Activation of aurora-A kinase by protein phosphatase inhibitor-2, a bifunctional signaling protein. Proc Natl Acad Sci U S A. 2004, 101 (23): 8625-8630. 10.1073/pnas.0402966101.PubMed CentralPubMedView ArticleGoogle Scholar
- Li M, Stukenberg PT, Brautigan DL: Binding of phosphatase inhibitor-2 to prolyl isomerase pin1 modifies specificity for mitotic phosphoproteins. Biochemistry. 2007, 47 (1): 292-300.PubMedView ArticleGoogle Scholar
- Wang W, Stukenberg PT, Brautigan DL: Phosphatase inhibitor-2 balances protein phosphatase 1 and aurora B kinase for chromosome segregation and cytokinesis in human retinal epithelial cells. Mol Biol Cell. 2008, 19 (11): 4852-4862. 10.1091/mbc.E08-05-0460.PubMed CentralPubMedView ArticleGoogle Scholar
- Khandani A, Mohtashami M, Camirand A: Inhibitor-2 induced M-phase arrest in Xenopus cycling egg extracts is dependent on MAPK activation. Cell Mol Biol Lett. 2011, 16 (4): 669-688. 10.2478/s11658-011-0030-z.PubMedView ArticleGoogle Scholar
- Zambrano CA, Egaña JT, Núñez MT, Maccioni RB, González-Billault C: Oxidative stress promotes τ dephosphorylation in neuronal cells: the roles of cdk5 and PP1. Free Radic Biol Med. 2004, 36 (11): 1393-1402. 10.1016/j.freeradbiomed.2004.03.007.PubMedView ArticleGoogle Scholar
- Vijayaraghavan S, Stephens DT, Trautman K, Smith GD, Khatra B, da Cruz e Silva EF, Greengard P: Sperm motility development in the epididymis is associated with decreased glycogen synthase kinase-3 and protein phosphatase 1 activity. Biol Reprod. 1996, 54 (3): 709-718. 10.1095/biolreprod54.3.709.PubMedView ArticleGoogle Scholar
- Smith GD, Wolf DP, Trautman KC, da Cruz e Silva EF, Greengard P, Vijayaraghavan S: Primate sperm contain protein phosphatase 1, a biochemical mediator of motility. Biol Reprod. 1996, 54 (3): 719-727. 10.1095/biolreprod54.3.719.PubMedView ArticleGoogle Scholar
- Ceulemans H, Stalmans W, Bollen M: Regulator-driven functional diversification of protein phosphatase-1 in eukaryotic evolution. Bioessays. 2002, 24 (4): 371-381. 10.1002/bies.10069.PubMedView ArticleGoogle Scholar
- Li M, Satinover DL, Brautigan DL: Phosphorylation and functions of inhibitor-2 family of proteins. Biochemistry. 2007, 46 (9): 2380-2389. 10.1021/bi602369m.PubMedView ArticleGoogle Scholar
- Maestre J, Tchenio T, Dhellin O, Heidmann T: mRNA retroposition in human cells: processed pseudogene formation. EMBO J. 1995, 14 (24): 6333-6338.PubMed CentralPubMedGoogle Scholar
- Weiner AM, Deininger PL, Efstratiadis A: Nonviral retroposons: genes, pseudogenes, and transposable elements generated by the reverse flow of genetic information. Annu Rev Biochem. 1986, 55 (1): 631-661. 10.1146/annurev.bi.55.070186.003215.PubMedView ArticleGoogle Scholar
- Wu I, Moses MA: Cloning of a cDNA encoding an isoform of human protein phosphatase inhibitor 2 from vascularized breast tumor. DNA Seq. 2001, 11 (6): 515-518.PubMedGoogle Scholar
- Yamada Y, Watanabe H, Miura F, Soejima H, Uchiyama M, Iwasaka T, Mukai T, Sakaki Y, Ito T: A comprehensive analysis of allelic methylation status of CpG islands on human chromosome 21q. Genome Res. 2004, 14 (2): 247-266. 10.1101/gr.1351604.PubMed CentralPubMedView ArticleGoogle Scholar
- Shirato H, Shima H, Sakashita G, Nakano T, Ito M, Lee EYC, Kikuchi K: Identification and characterization of a novel protein inhibitor of type 1 protein phosphatase. Biochemistry. 2000, 39 (45): 13848-13855. 10.1021/bi001326n.PubMedView ArticleGoogle Scholar
- Gerhard DS, Wagner L, Feingold EA, Shenmen CM, Grouse LH, Schuler G, Klein SL, Old S, Rasooly R, Good P, et al: The status, quality, and expansion of the NIH full-length cDNA project: the mammalian gene collection (MGC). Genome Res. 2004, 14 (10B): 2121-2127.PubMedView ArticleGoogle Scholar
- Eilon T, Barash I: Distinct gene-expression profiles characterize mammary tumors developed in transgenic mice expressing constitutively active and C-terminally truncated variants of STAT5. BMC Genomics. 2009, 10 (1): 231-10.1186/1471-2164-10-231.PubMed CentralPubMedView ArticleGoogle Scholar
- Fardilha M, Esteves SL, Korrodi-Gregório L, Vintém AP, Domingues SC, Rebelo S, Morrice N, Cohen PT, da Cruz e Silva OA, da Cruz e Silva EF: Identification of the human testis protein phosphatase 1 interactome. Biochem Pharmacol. 2011, 82 (10): 1403-1415. 10.1016/j.bcp.2011.02.018.PubMedView ArticleGoogle Scholar
- Murphy WJ, Eizirik E, O'Brien SJ, Madsen O, Scally M, Douady CJ, Teeling E, Ryder OA, Stanhope MJ, de Jong WW, et al: Resolution of the early placental mammal radiation using Bayesian phylogenetics. Science. 2001, 294 (5550): 2348-2351. 10.1126/science.1067179.PubMedView ArticleGoogle Scholar
- Fardilha M, Esteves SLC, Korrodi-Gregório L, Pelech S, da Cruz e Silva OAB, da Cruz e Silva E: Protein phosphatase 1 complexes modulate sperm motility and present novel targets for male infertility. Mol Hum Reprod. 2011, 17 (8): 466-477. 10.1093/molehr/gar004.PubMedView ArticleGoogle Scholar
- Mayya V, Han DK: Phosphoproteomics by mass spectrometry: insights, implications, applications and limitations. Expert Rev Proteomics. 2009, 6 (6): 605-618. 10.1586/epr.09.84.PubMed CentralPubMedView ArticleGoogle Scholar
- Gauci S, Helbig AO, Slijper M, Krijgsveld J, Heck AJR, Mohammed S: Lys-N and trypsin cover complementary parts of the phosphoproteome in a refined SCX-based approach. Anal Chem. 2009, 81 (11): 4493-4501. 10.1021/ac9004309.PubMedView ArticleGoogle Scholar
- Korrodi-Gregorio L, Ferreira M, Vintem AP, Wu W, Muller T, Marcus K, Vijayaraghavan S, Brautigan DL, da Cruz e Silva OA, Fardilha M, et al: Identification and characterization of two distinct PPP1R2 isoforms in human spermatozoa. BMC Cell Biol. 2013, 14: 15-10.1186/1471-2121-14-15.PubMed CentralPubMedView ArticleGoogle Scholar
- Esteves SL, Domingues SC, da Cruz e Silva OA, Fardilha M, da Cruz e Silva EF: Protein phosphatase 1alpha interacting proteins in the human brain. OMICS. 2012, 16 (1–2): 3-17.PubMed CentralPubMedView ArticleGoogle Scholar
- Montgomery SB, Sammeth M, Gutierrez-Arcelus M, Lach RP, Ingle C, Nisbett J, Guigo R, Dermitzakis ET: Transcriptome genetics using second generation sequencing in a Caucasian population. Nature. 2010, 464 (7289): 773-777. 10.1038/nature08903.PubMedView ArticleGoogle Scholar
- Ohshima K, Hattori M, Yada T, Gojobori T, Sakaki Y, Okada N: Whole-genome screening indicates a possible burst of formation of processed pseudogenes and Alu repeats by particular L1 subfamilies in ancestral primates. Genome Biol. 2003, 4 (11): R74-10.1186/gb-2003-4-11-r74.PubMed CentralPubMedView ArticleGoogle Scholar
- Kleene KC, Mulligan E, Steiger D, Donohue K, Mastrangelo M-A: The mouse gene encoding the testis-specific isoform of Poly(A) binding protein (Pabp2) is an expressed retroposon: intimations that gene expression in spermatogenic cells facilitates the creation of new genes. J Mol Evol. 1998, 47 (3): 275-281. 10.1007/PL00006385.PubMedView ArticleGoogle Scholar
- Marques AC, Dupanloup I, Vinckenbosch N, Reymond A, Kaessmann H: Emergence of young human genes after a burst of retroposition in primates. PLoS Biol. 2005, 3 (11): e357-10.1371/journal.pbio.0030357.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang C-J, Lin W-Y, Chang C-M, Choo K-B: Transcription of the rat testis-specific Rtdpoz-T1 and -T2 retrogenes during embryo development: co-transcription and frequent exonisation of transposable element sequences. BMC Mol Biol. 2009, 10 (1): 74-10.1186/1471-2199-10-74.PubMed CentralPubMedView ArticleGoogle Scholar
- Hurley TD, Yang J, Zhang L, Goodwin KD, Zou Q, Cortese M, Dunker AK, DePaoli-Roach AA: Structural basis for regulation of protein phosphatase 1 by inhibitor-2. J Biol Chem. 2007, 282 (39): 28874-28883. 10.1074/jbc.M703472200.PubMedView ArticleGoogle Scholar
- Yang J, Hurley TD, DePaoli-Roach AA: Interaction of inhibitor-2 with the catalytic subunit of type 1 protein phosphatase. J Biol Chem. 2000, 275 (30): 22635-22644. 10.1074/jbc.M003082200.PubMedView ArticleGoogle Scholar
- Varmuza S, Jurisicova A, Okano K, Hudson J, Boekelheide K, Shipp EB: Spermiogenesis is impaired in mice bearing a targeted mutation in the protein phosphatase 1cgamma gene. Dev Biol. 1999, 205 (1): 98-110. 10.1006/dbio.1998.9100.PubMedView ArticleGoogle Scholar
- Chakrabarti R, Cheng L, Puri P, Soler D, Vijayaraghavan S: Protein phosphatase PP1[gamma]2 in sperm morphogenesis and epididymal initiation of sperm motility. Asian J Androl. 2007, 9 (4): 445-452. 10.1111/j.1745-7262.2007.00307.x.PubMedView ArticleGoogle Scholar
- Vinckenbosch N, Dupanloup I, Kaessmann H: Evolutionary fate of retroposed gene copies in the human genome. Proc Natl Acad Sci U S A. 2006, 103 (9): 3220-3225. 10.1073/pnas.0511307103.PubMed CentralPubMedView ArticleGoogle Scholar
- Khachane A, Harrison P: Assessing the genomic evidence for conserved transcribed pseudogenes under selection. BMC Genomics. 2009, 10 (1): 435-10.1186/1471-2164-10-435.PubMed CentralPubMedView ArticleGoogle Scholar
- Torgerson DG, Singh RS: Enhanced adaptive evolution of sperm-expressed genes on the mammalian X chromosome. Heredity. 2006, 96 (1): 39-44.PubMedGoogle Scholar
- Dorus S, Wasbrough ER, Busby J, Wilkin EC, Karr TL: Sperm proteomics reveals intensified selection on mouse sperm membrane and acrosome genes. Mol Biol Evol. 2010, 27 (6): 1235-1246. 10.1093/molbev/msq007.PubMed CentralPubMedView ArticleGoogle Scholar
- Pandita A, Balasubramaniam A, Perrin R, Shannon P, Guha A: Malignant and benign ganglioglioma: a pathological and molecular study. Neuro Oncol. 2007, 9 (2): 124-134. 10.1215/15228517-2006-029.PubMed CentralPubMedView ArticleGoogle Scholar
- Pink RC, Wicks K, Caley DP, Punch EK, Jacobs L, Carter DR: Pseudogenes: pseudo-functional or key regulators in health and disease?. RNA. 2011, 17 (5): 792-798. 10.1261/rna.2658311.PubMed CentralPubMedView ArticleGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008, 25 (7): 1253-1256. 10.1093/molbev/msn083.PubMedView ArticleGoogle Scholar
- Zwickl DJ: Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion, Dissertation. 2006, Austin, TX: University of TexasGoogle Scholar
- Swofford DL: PAUP*. phylogenetic analysis using parsimony (*and other methods) version 4. 2003, Sunderland, MA: Sinauer AssociatesGoogle Scholar
- Hedges SB, Dudley J, Kumar S: TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics. 2006, 22 (23): 2971-2972. 10.1093/bioinformatics/btl505.PubMedView ArticleGoogle Scholar
- Poon AF, Frost SD, Pond SL: Detecting signatures of selection from DNA sequences using Datamonkey. Methods Mol Biol. 2009, 537: 163-183. 10.1007/978-1-59745-251-9_8.PubMedView ArticleGoogle Scholar
- Pond SL, Frost SD: Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics. 2005, 21 (10): 2531-2533. 10.1093/bioinformatics/bti320.PubMedView ArticleGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
- Yang Z, Nielsen R, Goldman N, Pedersen AM: Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000, 155 (1): 431-449.PubMed CentralPubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.PubMed CentralPubMedView ArticleGoogle Scholar
- Zheng D, Frankish A, Baertsch R, Kapranov P, Reymond A, Choo SW, Lu Y, Denoeud F, Antonarakis SE, Snyder M, et al: Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution. Genome Res. 2007, 17 (6): 839-851. 10.1101/gr.5586307.PubMed CentralPubMedView ArticleGoogle Scholar
- World Health Organization: WHO Laboratory for examination of human semen and sperm-cervical mucus interaction. Collection and examination of human semen. 1999, Cambridge: Cambridge University Press, 4-30.Google Scholar
- Thiele H, Glandorf J, Hufnagel P: Bioinformatics strategies in life sciences: from data processing and data warehousing to biological knowledge extraction. J Integr Bioinform. 2010, 7 (1): 141-PubMedGoogle Scholar
- Perkins DN, Pappin DJC, Creasy DM, Cottrell JS: Probability-based protein identification by searching sequence databases using mass spectrometry data. Electrophoresis. 1999, 20 (18): 3551-3567. 10.1002/(SICI)1522-2683(19991201)20:18<3551::AID-ELPS3551>3.0.CO;2-2.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.