Identification and characterisation of a cytochrome P450 gene and processed pseudogene from an arachnid: the cattle tick, Boophilus microplus

We isolated and sequenced the first known cytochrome P450 gene and pseudogene from an arachnid, the cattle tick, Boophilus microplus. Bath the gene and pseudogene belong to the family CYP4, but a new subfamily, CYP4W, had to be created for these genes because they are substantially different to other CYP4 genes. The gene, CPP4W1, has greatest homology with CYP4C1 from a cockroach, Blaberus discoidalis. The predicted molecular weight of the protein encoded by CYP4W1 (63 KDa) is greater than that of the other CYP4 genes. The pseudogene, CYP4W1P, is probably a processed pseudogene derived from the functional gene CYP4W1. This is only the third CYP processed pseudogene to be identified. The pseudogene is 98% identical to the functional gene, CYP4W1, therefore we hypothesise that this pseudogene evolved recently from the functional gene. The CYP4 genes from arthropods have diverged from each other more than those of mammals; consequently the phylogeny of the arthropod genes could not be resolved. (C) 1999 Elsevier Science Ltd. All rights reserved.

To be or not to be a pseudogene: a molecular epidemiological approach to the mclx genes and its impact in tuberculosis

Tuberculosis presents a myriad of symptoms, progression routes and propagation patterns not yet fully understood. Whereas for a long time research has focused solely on the patient immunity and overall susceptibility, it is nowadays widely accepted that the genetic diversity of its causative agent, Mycobacterium tuberculosis, plays a key role in this dynamic. This study focuses on a particular family of genes, the mclxs (Mycobacterium cyclase/LuxR-like genes), which codify for a particular and nearly mycobacterial-exclusive combination of protein domains. mclxs genes were found to be pseudogenized by frameshift-causing insertion(s)/deletion(s) in a considerable number of M. tuberculosis complex strains and clinical isolates. To discern the functional implications of the pseudogenization, we have analysed the pattern of frameshift-causing mutations in a group of M. tuberculosis isolates while taking into account their microbial-, patient- and disease-related traits. Our logistic regression-based analyses have revealed disparate effects associated with the transcriptional inactivation of two mclx genes. In fact, mclx2 (Rv1358) pseudogenization appears to be primarily driven by the microbial phylogenetic background, being mainly related to the Euro-American (EAm) lineage; on the other hand, mclx3 (Rv2488c) presents a higher tendency for pseudogenization among isolates from patients born on the Western Pacific area, and from isolates causing extra-pulmonary infections. These results contribute to the overall knowledge on the biology of M. tuberculosis infection, whereas at the same time launch the necessary basis for the functional assessment of these so far overlooked genes.

Identification of Mycobacterium tuberculosis complex based on amplification and sequencing of the oxyR pseudogene from stored Ziehl-Neelsen-stained sputum smears in Brazil

A cross-sectional analysis of stored Ziehl-Neelsen (ZN)-stained sputum smear slides (SSS) obtained from two public tuberculosis referral laboratories located in Juiz de Fora, Minas Gerais, was carried out to distinguish Mycobacterium bovis from other members of the Mycobacterium tuberculosis complex (MTC). A two-step approach was used to distinguish M. bovis from other members of MTC: (i) oxyR pseudogene amplification to detect MTC and, subsequently, (ii) allele-specific sequencing based on the polymorphism at position 285 of this gene. The oxyR pseudogene was successfully amplified in 100 of 177 (56.5%) SSS available from 99 individuals. No molecular profile of M. bovis was found. Multivariate analysis indicated that acid-fast bacilli (AFB) results and the source laboratory were associated (p < 0.05) with oxyR pseudogene amplification. SSS that were AFB++ SSS showed more oxyR pseudogene amplification than those with AFB0, possibly due to the amount of DNA. One of the two source laboratories presented a greater chance of oxyR pseudogene amplification, suggesting that differences in sputum conservation between laboratories could have influenced the preservation of DNA. This study provides evidence that stored ZN-SSS can be used for the molecular detection of MTC.

False phylogenies on wood mice due to cryptic cytochrome-b pseudogene.

The phylogeny and phylogeography of the Old World wood mice (subgenus Sylvaemus, genus Apodemus, Muridae) are well-documented. Nevertheless, the distributions of species, such as A. fulvipectus and A. ponticus remain dubious, as well as their phylogenetic relationships with A. sylvaticus. We analysed samples of Apodemus spp. across Europe using the mitochondrial cytochrome-b gene (cyt-b) and compared the DNA and amino-acid compositions of previously published sequences. The main result stemming from this study is the presence of a well-differentiated lineage of Sylvaemus including samples of various species (A. sylvaticus, A. fulvipectus, A. ponticus) from distant locations, which were revealed to be nuclear copies of the mitochondrial cyt-b. The presence of this cryptic pseudogene in published sequences is supported by different pathways. This has led to important errors in previous molecular trees and hence to partial misinterpretations in the phylogeny of Apodemus.

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.

The GENCODE pseudogene resource.

BACKGROUND: Pseudogenes have long been considered as nonfunctional genomic sequences. However, recent evidence suggests that many of them might have some form of biological activity, and the possibility of functionality has increased interest in their accurate annotation and integration with functional genomics data. RESULTS: As part of the GENCODE annotation of the human genome, we present the first genome-wide pseudogene assignment for protein-coding genes, based on both large-scale manual annotation and in silico pipelines. A key aspect of this coupled approach is that it allows us to identify pseudogenes in an unbiased fashion as well as untangle complex events through manual evaluation. We integrate the pseudogene annotations with the extensive ENCODE functional genomics information. In particular, we determine the expression level, transcription-factor and RNA polymerase II binding, and chromatin marks associated with each pseudogene. Based on their distribution, we develop simple statistical models for each type of activity, which we validate with large-scale RT-PCR-Seq experiments. Finally, we compare our pseudogenes with conservation and variation data from primate alignments and the 1000 Genomes project, producing lists of pseudogenes potentially under selection. CONCLUSIONS: At one extreme, some pseudogenes possess conventional characteristics of functionality; these may represent genes that have recently died. On the other hand, we find interesting patterns of partial activity, which may suggest that dead genes are being resurrected as functioning non-coding RNAs. The activity data of each pseudogene are stored in an associated resource, psiDR, which will be useful for the initial identification of potentially functional pseudogenes.

Presence of the RHD pseudogene and the hybrid RHD-CE-Ds gene in Brazilians with the D-negative phenotype

The molecular basis for RHD pseudogene or RHDpsi is a 37-bp insertion in exon 4 of RHD. This insertion, found in two-thirds of D-negative Africans, appears to introduce a stop codon at position 210. The hybrid RHD-CE-Ds, where the 3' end of exon 3 and exons 4 to 8 are derived from RHCE, is associated with the VS+V- phenotype, and leads to a D-negative phenotype in people of African origin. We determined whether Brazilian blood donors of heterogeneous ethnic origin had RHDpsi and RHD-CE-Ds. DNA from 206 blood donors were tested for RHDpsi by a multiplex PCR that detects RHD, RHDpsi and the C and c alleles of RHCE. The RHD genotype was determined by comparison of size of amplified products associated with the RHD gene in both intron 4 and exon 10/3'-UTR. VS was determined by amplification of exon 5 of RHCE, and sequencing of PCR products was used to analyze C733G (Leu245Val). Twenty-two (11%) of the 206 D-negative Brazilians studied had the RHDpsi, 5 (2%) had the RHD-CE-Ds hybrid gene associated with the VS+V- phenotype, and 179 (87%) entirely lacked RHD. As expected, RHD was deleted in all the 50 individuals of Caucasian descent. Among the 156 individuals of African descent, 22 (14%) had inactive RHD and 3% had the RHD-CE-Ds hybrid gene. These data confirm that the inclusion of two different multiplex PCR for RHD is essential to test the D-negative Brazilian population in order to avoid false-positive typing of polytransfused patients and fetuses.

Pseudogene rescue: an adaptive mechanism of codon reassignment

Alterations to the genetic code – codon reassignments – have occurred many times in life’s history, despite the fact that genomes are coadapted to their genetic codes and therefore alterations are likely to be maladaptive. A potential mechanism for adaptive codon reassignment, which could trigger either a temporary period of codon ambiguity or a permanent genetic code change, is the reactivation of a pseudogene by a nonsense suppressor mutant transfer RNA. I examine the population genetics of each stage of this process and find that pseudogene rescue is plausible and also readily explains some features of extant variability in genetic codes.

Phylogenetic relationships among Brazilian howler monkeys, genus Alouatta (Platyrrhini, Atelidae), based on γ¹-globin pseudogene sequences

Os guaribas, do gênero Alouatta, que são os primatas do Novo Mundo com maior distribuição geográfica, têm sido colocados em três grupos de espécies: o grupo Alouatta palliata da América central, e os grupos sulamericanos Alouatta seniculus e Alouatta caraya. Este último é monotípico, mas o grupo A. seniculus inclui pelo menos três espécies (A. seniculus, A. belzebul e A. fusca). Neste estudo, foram seqüenciados aproximadamente 600 pares de base do pseudogene globina g1 nas quatro espécies brasileiras (A. seniculus, A. belzebul, A. fusca e A. caraya). Os métodos de máxima parcimônia e máxima verossimilhança produziram árvores filogenéticas com o mesmo arranjo: {A. caraya [A. seniculus (A. fusca, A. belzebul)]}. A árvore mais parcimoniosa apresentou valores de bootstrap maiores de 82% para todos os agrupamentos, e valores de força de ligação de pelo menos 2, apoiando o agrupamento irmão de A. fusca e A. belzebul. O estudo também confirmou a presença em A. fusca do elemento de inserção Alu, com 150 pares de base, e uma deleção de 1,8 kb no pseudogene globina g1 já conhecidos nas demais espécies de guaribas. A classificação cladística baseada em dados moleculares é congruente com as de estudos morfológicos, com um isolamento claro do grupo monoespecífico A. caraya em relação ao grupo A. seniculus.

8p23 beta-defensin copy number determination by single-locus pseudogene-based paralog ratio tests risk bias due to low-frequency sequence variations.

BACKGROUND The copy number variation (CNV) in beta-defensin genes (DEFB) on human chromosome 8p23 has been proposed to contribute to the phenotypic differences in inflammatory diseases. However, determination of exact DEFB CN is a major challenge in association studies. Quantitative real-time PCR (qPCR), paralog ratio tests (PRT) and multiplex ligation-dependent probe amplification (MLPA) have been extensively used to determine DEFB CN in different laboratories, but inter-method inconsistencies were observed frequently. In this study we asked which one is superior among the three methods for DEFB CN determination. RESULTS We developed a clustering approach for MLPA and PRT to statistically correlate data from a single experiment. Then we compared qPCR, a newly designed PRT and MLPA for DEFB CN determination in 285 DNA samples. We found MLPA had the best convergence and clustering results of the raw data and the highest call rate. In addition, the concordance rates between MLPA or PRT and qPCR (32.12% and 37.99%, respectively) were unacceptably low with underestimated CN by qPCR. Concordance rate between MLPA and PRT (90.52%) was high but PRT systematically underestimated CN by one in a subset of samples. In these samples a sequence variant which caused complete PCR dropout of the respective DEFB cluster copies was found in one primer binding site of one of the targeted paralogous pseudogenes. CONCLUSION MLPA is superior to PRT and even more to qPCR for DEFB CN determination. Although the applied PRT provides in most cases reliable results, such a test is particularly sensitive to low-frequency sequence variations preferably accumulating in loci like pseudogenes which are most likely not under selective pressure. In the light of the superior performance of multiplex assays, the drawbacks of such single PRTs could be overcome by combining more test markers.

Mouse trp2, the homologue of the human trpc2 pseudogene, encodes mTrp2, a store depletion-activated capacitative Ca2+ entry channel

Capacitative Ca2+ entry (CCE) is Ca2+ entering after stimulation of inositol 1,4,5-trisphosphate (IP3) formation and initiation of Ca2+ store depletion. One hallmark of CCE is that it can also be triggered merely by store depletion, as occurs after inhibition of internal Ca2+ pumps with thapsigargin. Evidence has accumulated in support of a role of transient receptor potential (Trp) proteins as structural subunits of a class of Ca2+-permeable cation channels activated by agonists that stimulate IP3 formation—very likely through a direct interaction between the IP3 receptor and a Trp subunit of the Ca2+ entry channel. The role of Trp’s in Ca2+ entry triggered by store depletion alone is less clear. Only a few of the cloned Trp’s appear to enhance this type of Ca2+ entry, and when they do, the effect requires special conditions to be observed, which native CCE does not. Here we report the full-length cDNA of mouse trp2, the homologue of the human trp2 pseudogene. Mouse Trp2 is shown to be readily activated not only after stimulation with an agonist but also by store depletion in the absence of an agonist. In contrast to other Trp proteins, Trp2-mediated Ca2+ entry activated by store depletion is seen under the same conditions that reveal endogenous store depletion-activated Ca2+ entry, i.e., classical CCE. The findings support the general hypothesis that Trp proteins are subunits of store- and receptor-operated Ca2+ channels.

Digging for dead genes: an analysis of the characteristics of the pseudogene population in the Caenorhabditis elegans genome

Pseudogenes are non-functioning copies of genes in genomic DNA, which may either result from reverse transcription from an mRNA transcript (processed pseudogenes) or from gene duplication and subsequent disablement (non-processed pseudogenes). As pseudogenes are apparently ‘dead’, they usually have a variety of obvious disablements (e.g., insertions, deletions, frameshifts and truncations) relative to their functioning homologs. We have derived an initial estimate of the size, distribution and characteristics of the pseudogene population in the Caenorhabditis elegans genome, performing a survey in ‘molecular archaeology’. Corresponding to the 18 576 annotated proteins in the worm (i.e., in Wormpep18), we have found an estimated total of 2168 pseudogenes, about one for every eight genes. Few of these appear to be processed. Details of our pseudogene assignments are available from http://bioinfo.mbb.yale.edu/genome/worm/pseudogene. The population of pseudogenes differs significantly from that of genes in a number of respects: (i) pseudogenes are distributed unevenly across the genome relative to genes, with a disproportionate number on chromosome IV; (ii) the density of pseudogenes is higher on the arms of the chromosomes; (iii) the amino acid composition of pseudogenes is midway between that of genes and (translations of) random intergenic DNA, with enrichment of Phe, Ile, Leu and Lys, and depletion of Asp, Ala, Glu and Gly relative to the worm proteome; and (iv) the most common protein folds and families differ somewhat between genes and pseudogenes—whereas the most common fold found in the worm proteome is the immunoglobulin fold and the most common ‘pseudofold’ is the C-type lectin. In addition, the size of a gene family bears little overall relationship to the size of its corresponding pseudogene complement, indicating a highly dynamic genome. There are in fact a number of families associated with large populations of pseudogenes. For example, one family of seven-transmembrane receptors (represented by gene B0334.7) has one pseudogene for every four genes, and another uncharacterized family (represented by gene B0403.1) is approximately two-thirds pseudogenic. Furthermore, over a hundred apparent pseudogenic fragments do not have any obvious homologs in the worm.

Human steroidogenic acute regulatory protein: functional activity in COS-1 cells, tissue-specific expression, and mapping of the structural gene to 8p11.2 and a pseudogene to chromosome 13.

Steroidogenic acute regulatory protein (StAR) appears to mediate the rapid increase in pregnenolone synthesis stimulated by tropic hormones. cDNAs encoding StAR were isolated from a human adrenal cortex library. Human StAR, coexpressed in COS-1 cells with cytochrome P450scc and adrenodoxin, increased pregnenolone synthesis > 4-fold. A major StAR transcript of 1.6 kb and less abundant transcripts of 4.4 and 7.5 kb were detected in ovary and testis. Kidney had a lower amount of the 1.6-kb message. StAR mRNA was not detected in other tissues including placenta. Treatment of granulosa cells with 8-bromo-adenosine 3',5'-cyclic monophosphate for 24 hr increased StAR mRNA 3-fold or more. The structural gene encoding StAR was mapped using somatic cell hybrid mapping panels to chromosome 8p. Fluorescence in situ hybridization placed the StAR locus in the region 8p11.2. A StAR pseudogene was mapped to chromosome 13. We conclude that StAR expression is restricted to tissues that carry out mitochondrial sterol oxidations subject to acute regulation by cAMP and that StAR mRNA levels are regulated by cAMP.

A new confirmatory Neisseria gonorrhoeae real-time PCR assay targeting the porA pseudogene

The Roche Cobas Amplicor system is widely used for the detection of Neisseria gonorrhoeae but is known to cross react with some commensal Neisseria spp. Therefore, a confirmatory test is required. The most common target for confirmatory tests is the cppB gene of N. gonorrhoeae. However, the cppB gene is also present in other Neisseria spp. and is absent in some N. gonorrhoeae isolates. As a result, laboratories targeting this gene run the risk of obtaining both false-positive and false-negative results. In the study presented here, a newly developed N. gonorrhoeae LightCycler assay (NGpapLC) targeting the N. gonorrhoeae porA pseudogene was tested. The NGpapLC assay was used to test 282 clinical samples, and the results were compared to those obtained using a testing algorithm combining the Cobas Amplicor System (Roche Diagnostics, Sydney, Australia) and an in-house LightCycler assay targeting the cppB gene (cppB-LC). In addition, the specificity of the NGpapLC assay was investigated by testing a broad panel of bacteria including isolates of several Neisseria spp. The NGpapLC assay proved to have comparable clinical sensitivity to the cppB-LC assay. In addition; testing of the bacterial panel showed the NGpapLC assay to be highly specific for N. gonorrhoeae DNA. The results of this study show the NGpapLC assay is a suitable alternative to the cppB-LC assay for confirmation of N. gonorrhoeae-positive results obtained with Cobas Amplicor.