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

A large-scale phylogeny-guided analysis of pseudogenes in Pseudomonas aeruginosa bacterium

Pseudogenes, once considered "junk DNA" based on the incorrect assumption that the absence of full coding potential means a complete lack of functionality, have recently become a subject of significant interest in the scientific community. Concurrently, it is widely assumed that bacterial genomes are compact and have a high density of coding genes with little room for non-coding genes, including pseudogenes. A key aspect of genome annotation is the correct identification of genes and the distinction between coding genes and pseudogenes, as it directly impacts functional and comparative genomics studies. In this study, we analyzed the genomic data of 4,699 strains of the bacterium Pseudomonas aeruginosa (P. aeruginosa) as they exhibit high variability in the number of annotated pseudogenes. In particular, we looked for correlations between the number of pseudogenes and other genomic and meta-features of the strains. We identified clusters of orthologous genes and pseudogenes and compared cluster size distributions and length homogeneity within clusters. We then mapped and examined orthology relationships between genes and pseudogenes. Additionally, we generated a phylogenetic tree of the strains and found that phylogenetically related strains are more homogeneous in the number of pseudogenes and share a significant amount of pseudogenes. Finally, we delved into clusters of orthologous genes and pseudogenes and quantified their phylogenetic neighborhood, classifying pseudogenes into evolutionary preserved pseudogenes, mis-annotated pseudogenes, or pseudogenes formed by failed horizontal transfer events. This in-depth study provides important insights that can be incorporated into pseudogene annotation pipelines in the future. IMPORTANCE Accurate annotation of genes and pseudogenes is vital for comparative genomics analysis. Recent studies have shown that bacterial pseudogenes have an important role in regulatory processes and can provide insight into the evolutionary history of homologous genes or the genome as a whole. Due to pseudogenes' nature as non-functional genes, there is no commonly accepted definition of a pseudogene, which poses difficulties in verifying the annotation through experimental methods and resolving discrepancies among different annotation techniques. Our study introduces an in-depth analysis of annotated genes and pseudogenes and insights that can be incorporated into improved pseudogene annotation pipelines in the future.

Pseudogenes and their potential functions in hematopoiesis

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

Gauging the trends of pseudogenes in plants

Pseudogenes, the debilitated parts of ancient genes, were previously scrapped off as junk or discarded genes with no functional significance. Pseudogenes have come under scrutiny for their functionality, since recent studies have unveiled their importance in the regulation of their corresponding parent genes and various biological mechanisms. Despite the enormous occurrence of pseudogenes in plants, the lack of experimental validation has contributed toward their unresolved roles in gene regulation. Contrarily, most of the studies associated with gene regulation have been mainly reported for humans, mice, and other mammalian genomes. Consequently, in order to present a cumulative report on plant-based pseudogenes research, an attempt has been made to assemble multiple studies presenting the pseudogene classification, the prediction and the determination of comparative accuracies of various computational pipelines, and recent trends in analyzing their biological functions, and regulatory mechanisms. This review represents the classical, as well as the recent advances on pseudogene identification and their potential roles in transcriptional regulation, which could possibly invigorate the quality of genome annotation, evolutionary analysis, and complexity surrounding the regulatory pathways in plants. Thus, when the ambiguous boundary girdling the pseudogenes eventually recedes on account of their explicit orchestration role, research in flora would no longer saunter compared to that on fauna.

The pseudogene problem and RT-qPCR data normalization; SYMPK: a suitable reference gene for papillary thyroid carcinoma

In RT-qPCR, accuracy requires multiple levels of standardization, but results could be obfuscated by human errors and technical limitations. Data normalization against suitable reference genes is critical, yet their observed expression can be confounded by pseudogenes. Eight reference genes were selected based on literature review and analysis of papillary thyroid carcinoma (PTC) microarray data. RNA extraction and cDNA synthesis were followed by RT-qPCR amplification in triplicate with exon-junction or intron-spanning primers. Several statistical analyses were applied using Microsoft Excel, NormFinder, and BestKeeper. In normal tissues, the least correlation of variation (CqCV%) and the lowest maximum fold change (MFC) were respectively recorded for PYCR1 and SYMPK. In PTC tissues, SYMPK had the lowest CqCV% (5.16%) and MFC (1.17). According to NormFinder, the best reference combination was SYMPK and ACTB (stability value = 0.209). BestKeeper suggested SYMPK as the best reference in both normal (r = 0.969) and PTC tissues (r = 0.958). SYMPK is suggested as the best reference gene for overcoming the pseudogene problem in RT-qPCR data normalization, with a stability value of 0.319.

dreamBase: DNA modification, RNA regulation and protein binding of expressed pseudogenes in human health and disease

Although thousands of pseudogenes have been annotated in the human genome, their transcriptional regulation, expression profiles and functional mechanisms are largely unknown. In this study, we developed dreamBase (http://rna.sysu.edu.cn/dreamBase) to facilitate the investigation of DNA modification, RNA regulation and protein binding of potential expressed pseudogenes from multidimensional high-throughput sequencing data. Based on ∼5500 ChIP-seq and DNase-seq datasets, we identified genome-wide binding profiles of various transcription-associated factors around pseudogene loci. By integrating ∼18 000 RNA-seq data, we analysed the expression profiles of pseudogenes and explored their co-expression patterns with their parent genes in 32 cancers and 31 normal tissues. By combining microRNA binding sites, we demonstrated complex post-transcriptional regulation networks involving 275 microRNAs and 1201 pseudogenes. We generated ceRNA networks to illustrate the crosstalk between pseudogenes and their parent genes through competitive binding of microRNAs. In addition, we studied transcriptome-wide interactions between RNA binding proteins (RBPs) and pseudogenes based on 458 CLIP-seq datasets. In conjunction with epitranscriptome sequencing data, we also mapped 1039 RNA modification sites onto 635 pseudogenes. This database will provide insights into the transcriptional regulation, expression, functions and mechanisms of pseudogenes as well as their roles in biological processes and diseases.

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.

Evolutionary direction of processed pseudogenes

While some pseudogenes have been reported to play important roles in gene regulation, little is known about the possible relationship between pseudogene functions and evolutionary process of pseudogenes, or about the forces responsible for the pseudogene evolution. In this study, we characterized human processed pseudogenes in terms of evolutionary dynamics. Our results show that pseudogenes tend to evolve toward: lower GC content, strong dinucleotide bias, reduced abundance of transcription factor binding motifs and short palindromes, and decreased ability to form nucleosomes. We explored possible evolutionary forces that shaped the evolution pattern of pseudogenes, and concluded that mutations in pseudogenes are likely determined, at least partially, by neighbor-dependent mutational bias and recombination-associated selection.

H2A/K pseudogene mutation may promote cell proliferation

Little attention has been paid to the histone H2A/K pseudogene. Results from our laboratory showed that 7 of 10 kidney cancer patients carried a mutant H2A/K pseudogene; therefore, we were interested in determining the relationship between mutant H2A/K and cell proliferation. We used shotgun and label-free proteomics methods to study whether mutant H2A/K lncRNAs affected cell proliferation. Quantitative proteomic analysis indicated that the expression of mutant H2A/K lncRNAs resulted in the upregulation of many oncogenes, which promoted cell proliferation. Further interaction analyses revealed that a proliferating cell nuclear antigen (PCNA)-protein interaction network, with PCNA in the center, contributes to cell proliferation in cells expressing the mutant H2A/K lncRNAs. Western blotting confirmed the critical upregulation of PCNA by mutant H2A/K lncRNA expression. Finally, the promotion of cell proliferation by mutant H2A/K lncRNAs (C290T, C228A and A45G) was confirmed using cell proliferation assays. Although we did not determine the exact mechanism by which the oncogenes were upregulated by the mutant H2A/K lncRNAs, we confirmed that the mutant H2A/K lncRNAs promoted cell proliferation by upregulating PCNA and other oncogenes. The hypothesis that cell proliferation is promoted by the mutant H2A/K lncRNAs was supported by the protein expression and cell proliferation assay results. Therefore, mutant H2A/K lncRNAs may be a new factor in renal carcinogenesis.

Gene-pseudogene evolution: a probabilistic approach

Over the last decade, methods have been developed for the reconstruction of gene trees that take into account the species tree. Many of these methods have been based on the probabilistic duplication-loss model, which describes how a gene-tree evolves over a species-tree with respect to duplication and losses, as well as extension of this model, e.g., the DLRS (Duplication, Loss, Rate and Sequence evolution) model that also includes sequence evolution under relaxed molecular clock. A disjoint, almost as recent, and very important line of research has been focused on non protein-coding, but yet, functional DNA. For instance, DNA sequences being pseudogenes in the sense that they are not translated, may still be transcribed and the thereby produced RNA may be functional.

We extend the DLRS model by including pseudogenization events and devise an MCMC framework for analyzing extended gene families consisting of genes and pseudogenes with respect to this model, i.e., reconstructing gene-trees and identifying pseudogenization events in the reconstructed gene-trees. By applying the MCMC framework to biologically realistic synthetic data, we show that gene-trees as well as pseudogenization points can be inferred well. We also apply our MCMC framework to extended gene families belonging to the Olfactory Receptor and Zinc Finger superfamilies. The analysis indicate that both these super families contains very old pseudogenes, perhaps so old that it is reasonable to suspect that some are functional. In our analysis, the sub families of the Olfactory Receptors contains only lineage specific pseudogenes, while the sub families of the Zinc Fingers contains pseudogene lineages common to several species.

Nucleosome Organization around Pseudogenes in the Human Genome

Pseudogene, disabled copy of functional gene, plays a subtle role in gene expression and genome evolution. The first step in deciphering RNA-level regulation of pseudogenes is to understand their transcriptional activity. So far, there has been no report on possible roles of nucleosome organization in pseudogene transcription. In this paper, we investigated the effect of nucleosome positioning on pseudogene transcription. For transcribed pseudogenes, the experimental nucleosome occupancy shows a prominent depletion at the regions both upstream of pseudogene start positions and downstream of pseudogene end positions. Intriguingly, the same depletion is also observed for nontranscribed pseudogenes, which is unexpected since nucleosome depletion in those regions is thought to be unnecessary in light of the nontranscriptional property of those pseudogenes. The sequence-dependent prediction of nucleosome occupancy shows a consistent pattern with the experimental data-based analysis. Our results indicate that nucleosome positioning may play important roles in both the transcription initiation and termination of pseudogenes.

Construction and evaluation of a fluorescence-based live attenuated Escherichia coli delivery system for generating oral vaccine candidate

Enter toxigenic Escherichia coli (ETEC) is a major pathogen of swine industry that can have a substantial impact on morbidity and mortality. Therefore, it is necessary to develop effective vaccines for the prevention of ETEC infection. Live attenuated bacteria delivery system are effective tools for mucosal immunization. The purpose of this study was to construct a novel delivery system that can present the LTR192G-STb fusion protein as oral vaccine candidate. Firstly, the PRPL-mKate2 fluorescent cassette was inserted into the genome (yaiT pseudogene) of an attenuated E. coli by homologous recombination methods to construct the delivery system O142(yaiT::PRPL-mKate2). Secondly, the oral vaccine O142(yaiT:: LT192-STb) (ER-B) was derived for replacing the PRPL-mKate2 by LT192-STb fusion gene, and then it was tested for its feasibility as oral vaccine candidate. Subsequently, BALB/c mice were orogastrically immunized with ER-B. Results showed that mice orally immunized with ER-B produced high levels of specific IgA and IgG antibodies. The induced antibodies demonstrated neutralizing effects to enter toxins LT and STb. In addition, results of cellular immune responses showed that stimulation index values of immunized mice were significantly higher than the control group (P < 0.05) and with a marked shift towards Th 2 immunity. These data indicated that the recombinant E. coli ER-B could be a valuable candidate of future vaccines against ETEC infection.

Making novel proteins from pseudogenes

MOTIVATION

Recently, we made synthetic proteins from non-coding DNA of Escherichia coli. Encouraged by this, we asked: can we artificially express pseudogenes into novel and functional proteins? What kind of structures would be generated? Would these proteins be stable? How would the organism respond to the artificial reactivation of pseudogenes?

RESULTS

To answer these questions, we studied 16 full-length protein equivalents of pseudogenes. The sequence-based predictions indicated interesting molecular and cellular functional roles for pseudogene-derived proteins. Most of the proteins were predicted to be involved in the amino acid biosynthesis, energy metabolism, purines and pyrimidine biosynthesis, central intermediary metabolism, transport and binding. Interestingly, many of the pseudogene-derived proteins were predicted to be enzymes. Furthermore, proteins showed strong evidence of stable tertiary structures. The prediction scores for structure, function and stability were found to be favorable in most of the cases.

IMPACT

To our best knowledge, this is the first such report that predicts the possibility of making functional and stable proteins from pseudogenes. In future, it would be interesting to experimentally synthesize and validate these predictions.

Comparative analysis of pseudogenes across three phyla

Pseudogenes are degraded fossil copies of genes. Here, we report a comparison of pseudogenes spanning three phyla, leveraging the completed annotations of the human, worm, and fly genomes, which we make available as an online resource. We find that pseudogenes are lineage specific, much more so than protein-coding genes, reflecting the different remodeling processes marking each organism's genome evolution. The majority of human pseudogenes are processed, resulting from a retrotranspositional burst at the dawn of the primate lineage. This burst can be seen in the largely uniform distribution of pseudogenes across the genome, their preservation in areas with low recombination rates, and their preponderance in highly expressed gene families. In contrast, worm and fly pseudogenes tell a story of numerous duplication events. In worm, these duplications have been preserved through selective sweeps, so we see a large number of pseudogenes associated with highly duplicated families such as chemoreceptors. However, in fly, the large effective population size and high deletion rate resulted in a depletion of the pseudogene complement. Despite large variations between these species, we also find notable similarities. Overall, we identify a broad spectrum of biochemical activity for pseudogenes, with the majority in each organism exhibiting varying degrees of partial activity. In particular, we identify a consistent amount of transcription (∼15%) across all species, suggesting a uniform degradation process. Also, we see a uniform decay of pseudogene promoter activity relative to their coding counterparts and identify a number of pseudogenes with conserved upstream sequences and activity, hinting at potential regulatory roles.

Trichinella pseudospiralis vs. T. spiralis thymidylate synthase gene structure and T. pseudospiralis thymidylate synthase retrogene sequence

Background

Thymidylate synthase is a housekeeping gene, designated ancient due to its role in DNA synthesis and ubiquitous phyletic distribution. The genomic sequences were characterized coding for thymidylate synthase in two species of the genus Trichinella, an encapsulating T. spiralis and a non-encapsulating T. pseudospiralis.

Methods

Based on the sequence of parasitic nematode Trichinella spiralis thymidylate synthase cDNA, PCR techniques were employed.

Results

Each of the respective gene structures encompassed 6 exons and 5 introns located in conserved sites. Comparison with the corresponding gene structures of other eukaryotic species revealed lack of common introns that would be shared among selected fungi, nematodes, mammals and plants. The two deduced amino acid sequences were 96% identical. In addition to the thymidylate synthase gene, the intron-less retrocopy, i.e. a processed pseudogene, with sequence identical to the T. spiralis gene coding region, was found to be present within the T. pseudospiralis genome. This pseudogene, instead of the gene, was confirmed by RT-PCR to be expressed in the parasite muscle larvae.

Conclusions

Intron load, as well as distribution of exon and intron phases in thymidylate synthase genes from various sources, point against the theory of gene assembly by the primordial exon shuffling and support the theory of evolutionary late intron insertion into spliceosomal genes. Thymidylate synthase pseudogene expressed in T. pseudospiralis muscle larvae is designated a retrogene.

Systems biology approach to integrative comparative genomics

As more and more systems biology approaches are used to investigate the different types of biological macromolecules, increasing numbers of whole genomic studies are now available for a large array of organisms. Whether it is genomics, transcriptomics, proteomics, interactomics or metabolomics, the full complement of genomic information on all different levels can be juxtaposed between different organisms to reveal similarities or differences, and even to provide consensus models. At the intersection of comparative genomics and systems biology lies great possibility for discovery, analysis and prediction. This paper explores this nexus and the relationship from four general levels: DNA, RNA, protein and extragenomic. For each level, we provide an overview of the methods, discuss the potential challenges and survey the current research. Finally, we suggest some organizing principles and make proposals for new areas that will be important for future research.

Computational methods of identification of pseudogenes based on functionality: entropy and GC content

Spectral entropy and GC content analyses reveal comprehensive structural features of DNA sequences. To illustrate the significance of these features, we analyze the β-esterase gene cluster, including the Est-6 gene and the ψEst-6 putative pseudogene, in seven species of the Drosophila melanogaster subgroup. The spectral entropies show distinctly lower structural ordering for ψEst-6 than for Est-6 in all species studied. However, entropy accumulation is not a completely random process for either gene and it shows to be nucleotide dependent. Furthermore, GC content in synonymous positions is uniformly higher in Est-6 than in ψEst-6, in agreement with the reduced GC content generally observed in pseudogenes and nonfunctional sequences. The observed differences in entropy and GC content reflect an evolutionary shift associated with the process of pseudogenization and subsequent functional divergence of ψEst-6 and Est-6 after the duplication event. The data obtained show the relevance and significance of entropy and GC content analyses for pseudogene identification and for the comparative study of gene-pseudogene evolution.

Evolutionarily dynamic roles of a PUF RNA-binding protein in the somatic development of Caenorhabditis briggsae

Gene duplication and divergence has emerged as an important aspect of developmental evolution. The genomes of Caenorhabditis nematodes encode an ancient family of PUF RNA-binding proteins. Most have been implicated in germline development, and are often redundant with paralogs of the same sub-family. An exception is Cbr-puf-2 (one of three Caenorhabditis briggsae PUF-2 sub-family paralogs), which is required for development past the second larval stage. Here, we provide a detailed functional characterization of Cbr-puf-2. The larval arrest of Cbr-puf-2 mutant animals is caused by inefficient breakdown of bacterial food, which leads to starvation. Cbr-puf-2 is required for the normal grinding cycle of the muscular terminal bulb during early larval stages, and is transiently expressed in this tissue. In addition, rescue of larval arrest reveals that Cbr-puf-2 also promotes normal vulval development. It is expressed in the anchor cell (which induces vulval fate) and vulval muscles, but not in the vulva precursor cells (VPCs) themselves. This contrasts with the VPC-autonomous repression of vulval development described for the Caenorhabditis elegans homologs fbf-1/2. These different roles for PUF proteins occur even as the vulva and pharynx maintain highly conserved anatomies across Caenorhabditis, indicating pervasive developmental system drift (DSD). Because Cbr-PUF-2 shares RNA-binding specificity with its paralogs and with C. elegans FBF, we suggest that functional novelty of RNA-binding proteins evolves through changes in the site of their expression, perhaps in concert with cis-regulatory evolution in target mRNAs.

Human Transcriptome and Chromatin Modifications: An ENCODE Perspective

A decade-long project, led by several international research groups, called the Encyclopedia of DNA Elements (ENCODE), recently released an unprecedented amount of data. The ambitious project covers transcriptome, cistrome, epigenome, and interactome data from more than 1,600 sets of experiments in human. To make use of this valuable resource, it is important to understand the information it represents and the techniques that were used to generate these data. In this review, we introduce the data that ENCODE generated, summarize the observations from the data analysis, and revisit a computational approach that ENCODE used to predict gene expression, with a focus on the human transcriptome and its association with chromatin modifications.

Pseudogenes: pseudo or real functional elements?

Pseudogenes are genomic remnants of ancient protein-coding genes which have lost their coding potentials through evolution. Although broadly existed, pseudogenes used to be considered as junk or relics of genomes which have not drawn enough attentions of biologists until recent years. With the broad applications of high-throughput experimental techniques, growing lines of evidence have strongly suggested that some pseudogenes possess special functions, including regulating parental gene expression and participating in the regulation of many biological processes. In this review, we summarize some basic features of pseudogenes and their functions in regulating development and diseases. All of these observations indicate that pseudogenes are not purely dead fossils of genomes, but warrant further exploration in their distribution, expression regulation and functions. A new nomenclature is desirable for the currently called 'pseudogenes' to better describe their functions.

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.

The origin and evolution of new genes

New genes are a major source of genetic innovation in genomes. However, until recently, understanding how new genes originate and how they evolve was hampered by the lack of appropriate genetic datasets. The advent of the genomic era brought about a revolution in the amount of data available to study new genes. For the first time, decades-old theoretical principles could be tested empirically and novel and unexpected avenues of research opened up. This chapter explores how genomic data can and is being used to study both the origin and evolution of new genes and the surprising discoveries made thus far.

Pseudogenes: pseudo-functional or key regulators in health and disease?

Pseudogenes have long been labeled as "junk" DNA, failed copies of genes that arise during the evolution of genomes. However, recent results are challenging this moniker; indeed, some pseudogenes appear to harbor the potential to regulate their protein-coding cousins. Far from being silent relics, many pseudogenes are transcribed into RNA, some exhibiting a tissue-specific pattern of activation. Pseudogene transcripts can be processed into short interfering RNAs that regulate coding genes through the RNAi pathway. In another remarkable discovery, it has been shown that pseudogenes are capable of regulating tumor suppressors and oncogenes by acting as microRNA decoys. The finding that pseudogenes are often deregulated during cancer progression warrants further investigation into the true extent of pseudogene function. In this review, we describe the ways in which pseudogenes exert their effect on coding genes and explore the role of pseudogenes in the increasingly complex web of noncoding RNA that contributes to normal cellular regulation.

Two complementary recessive genes in duplicated segments control etiolation in rice

The main objective of this study was to identify the genes causing etiolation in a rice mutant, the thylakoids of which were scattered. Three populations were employed to map the genes for etiolation using bulked segregant analysis. Genetic analysis confirmed that etiolation was controlled by two recessive genes, et11 and et12, which were fine mapped to an approximately 147-kb region and an approximately 209-kb region on the short arms of chromosomes 11 and 12, respectively. Both regions were within the duplicated segments on chromosomes 11 and 12. They possessed a highly similar sequence of 38 kb at the locations of a pair of duplicated genes with protein sequences very similar to that of HCF152 in Arabidopsis that are required for the processing of chloroplast RNA. These genes are likely the candidates for et11 and et12. Expression profiling was used to compare the expression patterns of paralogs in the duplicated segments. Expression profiling indicated that the duplicated segments had been undergone concerted evolution, and a large number of the paralogs within the duplicated segments were functionally redundant like et11 and et12.

Processed pseudogenes are located preferentially in regions of low recombination rates in the human genome

The aim of this article is to demonstrate possible recombination-associated evolutionary forces affecting the genomic distribution of processed pseudogenes. The relationship between recombination rate and the distribution of processed pseudogenes is analysed in the human genome. The results show that processed pseudogenes preferentially accumulate in regions of low recombination rates and this correlation cannot be explained by indirect relationships with GC content and gene density. Several explanatory models for the observation are discussed. A model of selection against ectopic recombination is tested based on the difference in distribution pattern between two classes of processed pseudogenes, which differ in the possibility of stimulating ectopic recombination. Our results indicate that the correlation between processed pseudogene density and recombination rate is probably results, in part, from the selection against ectopic recombination between closely located homologous processed pseudogenes. We also found a length effect in processed pseudogene distribution, namely long processed pseudogenes are located more preferentially in regions of low recombination rates than short ones.

The ABC transporter gene family of Daphnia pulex

Background

The large gene superfamily of ABC (ATP-binding cassette) transporters encodes membrane proteins involved in trafficking processes across biological membranes and further essential cell biological functions. ABC transporters are evolutionary ancient and involved in the biochemical defence against toxicants. We report here a genome-wide survey of ABC proteins of Daphnia pulex, providing for the first time information on ABC proteins in crustacea, a primarily aquatic arthropod subphylum of high ecological and economical importance.

Results

We identified 64 ABC proteins in the Daphnia genome, which possesses members of all current ABC subfamilies A to H. To unravel phylogenetic relationships, ABC proteins of Daphnia were compared to those from yeast, worm, fruit fly and human. A high conservation of Daphnia of ABC transporters was observed for proteins involved in fundamental cellular processes, including the mitochondrial half transporters of the ABCB subfamily, which function in iron metabolism and transport of Fe/S protein precursors, and the members of subfamilies ABCD, ABCE and ABCF, which have roles in very long chain fatty acid transport, initiation of gene transcription and protein translation, respectively. A number of Daphnia proteins showed one-to-one orthologous relationships to Drosophila ABC proteins including the sulfonyl urea receptor (SUR), the ecdysone transporter ET23, and the eye pigment precursor transporter scarlet. As the fruit fly, Daphnia lacked homologues to the TAP protein, which plays a role in antigene processing, and the cystic fibrosis transmembrane conductance regulator (CFTR), which functions as a chloride channel. Daphnia showed two proteins homologous to MDR (multidrug resistance) P-glycoproteins (ABCB subfamily) and six proteins homologous to MRPs (multidrug resistance-associated proteins) (ABCC subfamily). However, lineage specific gene duplications in the ABCB and ABCC subfamilies complicated the inference of function. A particularly high number of gene duplications were observed in the ABCG and ABCH subfamilies, which have 23 and 15 members, respectively.

Conclusion

The in silico characterisation of ABC transporters in the Daphnia pulex genome revealed that the complement of ABC transporters is as complex in crustaceans as that other metazoans. Not surprisingly, among currently available genomes, Daphnia ABC transporters most closely resemble those of the fruit fly, another arthropod.

Evolution of the Caenorhabditis elegans genome

A fundamental problem in genome biology is to elucidate the evolutionary forces responsible for generating nonrandom patterns of genome organization. As the first metazoan to benefit from full-genome sequencing, Caenorhabditis elegans has been at the forefront of research in this area. Studies of genomic patterns, and their evolutionary underpinnings, continue to be augmented by the recent push to obtain additional full-genome sequences of related Caenorhabditis taxa. In the near future, we expect to see major advances with the onset of whole-genome resequencing of multiple wild individuals of the same species. In this review, we synthesize many of the important insights to date in our understanding of genome organization and function that derive from the evolutionary principles made explicit by theoretical population genetics and molecular evolution and highlight fertile areas for future research on unanswered questions in C. elegans genome evolution. We call attention to the need for C. elegans researchers to generate and critically assess nonadaptive hypotheses for genomic and developmental patterns, in addition to adaptive scenarios. We also emphasize the potential importance of evolution in the gonochoristic (female and male) ancestors of the androdioecious (hermaphrodite and male) C. elegans as the source for many of its genomic and developmental patterns.

Comparative genomics and molecular dynamics of DNA repeats in eukaryotes

Repeated elements can be widely abundant in eukaryotic genomes, composing more than 50% of the human genome, for example. It is possible to classify repeated sequences into two large families, "tandem repeats" and "dispersed repeats." Each of these two families can be itself divided into subfamilies. Dispersed repeats contain transposons, tRNA genes, and gene paralogues, whereas tandem repeats contain gene tandems, ribosomal DNA repeat arrays, and satellite DNA, itself subdivided into satellites, minisatellites, and microsatellites. Remarkably, the molecular mechanisms that create and propagate dispersed and tandem repeats are specific to each class and usually do not overlap. In the present review, we have chosen in the first section to describe the nature and distribution of dispersed and tandem repeats in eukaryotic genomes in the light of complete (or nearly complete) available genome sequences. In the second part, we focus on the molecular mechanisms responsible for the fast evolution of two specific classes of tandem repeats: minisatellites and microsatellites. Given that a growing number of human neurological disorders involve the expansion of a particular class of microsatellites, called trinucleotide repeats, a large part of the recent experimental work on microsatellites has focused on these particular repeats, and thus we also review the current knowledge in this area. Finally, we propose a unified definition for mini- and microsatellites that takes into account their biological properties and try to point out new directions that should be explored in a near future on our road to understanding the genetics of repeated sequences.

Retrosequence formation restructures the yeast genome

Retrosequences generated by reverse transcription of mRNA transcripts have a substantial influence on gene expression patterns, generation of novel gene functions, and genome organization. The Ty1 retrotransposon is a major source of RT activity in the yeast, Saccharomyces cerevisiae, and Ty1 retromobility is greatly elevated in strains lacking telomerase. We report that Ty1-dependent formation of retrosequences derived from single-copy gene transcripts is progressively elevated as yeast cells senesce in the absence of telomerase. Retrosequences are frequently fused to Ty1 sequences, and occasionally to sequences from other mRNA transcripts, forming chimeric pseudogenes. Efficient retrosequence formation requires the homologous recombination gene RAD52. Selection for retrosequence formation is correlated with a high frequency of chromosome rearrangements in telomerase-negative yeast. Ty1-associated retrosequences were present at the breakpoint junctions of four chromosomes analyzed in detail. Our results support a role for reverse transcripts in promoting chromosome rearrangements.

Gene decay in archaea

The gene-dense chromosomes of archaea and bacteria were long thought to be devoid of pseudogenes, but with the massive increase in available genome sequences, whole genome comparisons between closely related species have identified mutations that have rendered numerous genes inactive. Comparative analyses of sequenced archaeal genomes revealed numerous pseudogenes, which can constitute up to 8.6% of the annotated coding sequences in some genomes. The largest proportion of pseudogenes is created by gene truncations, followed by frameshift mutations. Within archaeal genomes, large numbers of pseudogenes contain more than one inactivating mutation, suggesting that pseudogenes are deleted from the genome more slowly in archaea than in bacteria. Although archaea seem to retain pseudogenes longer than do bacteria, most archaeal genomes have unique repertoires of pseudogenes.

FGF: a web tool for Fishing Gene Family in a whole genome database

Gene duplication is an important process in evolution. The availability of genome sequences of a number of organisms has made it possible to conduct comprehensive searches for duplicated genes enabling informative studies of their evolution. We have established the FGF (Fishing Gene Family) program to efficiently search for and identify gene families. The FGF output displays the results as visual phylogenetic trees including information on gene structure, chromosome position, duplication fate and selective pressure. It is particularly useful to identify pseudogenes and detect changes in gene structure. FGF is freely available on a web server at http://fgf.genomics.org.cn/

Origins and impact of constraints in evolution of gene families

Recent investigations of high-throughput genomic and phenomic data have uncovered a variety of significant but relatively weak correlations between a gene's functional and evolutionary characteristics. In particular, essential genes and genes with paralogs have a slight propensity to evolve more slowly than nonessential genes and singletons, respectively. However, given the weakness and multiplicity of these associations, their biological relevance remains uncertain. Here, we show that existence of an essential paralog can be used as a specific and strong gauge of selection. We partition gene families in several genomes into two classes: those that include at least one essential gene (E-families) and those without essential genes (N-families). We find that weaker purifying selection causes N-families to evolve in a more dynamic regime with higher rates both of duplicate fixation and pseudogenization. Because genes in E-families are subject to significantly stronger purifying selection than those in N-families, they survive longer and exhibit greater sequence divergence. Longer average survival time also allows for divergence of upstream regulatory regions, resulting in change of transcriptional context among paralogs in E-families. These findings are compatible with differential division of ancestral functions (subfunctionalization) or emergence of novel functions (neofunctionalization) being the prevalent modes of evolution of paralogs in E-families as opposed to pseudogenization (nonfunctionalization), which is the typical fate of paralogs in N-families. Unlike other characteristics of genes, such as essentiality, existence of paralogs, or expression level, membership in an E-family or an N-family strongly correlates with the level of selection and appears to be a major determinant of a gene's evolutionary fate.

Processed pseudogenes are more abundant in human and mouse X chromosomes than in autosomes

Two different hypotheses have been proposed to explain the observation that some genomes contain more processed pseudogenes than others. One predicts that processed pseudogene abundance is inversely proportional to the substrate specificity of the reverse transcriptase that generates processed pseudogenes. The other predicts that the amount of processed pseudogenes found in genomes is proportional to the length of oogenesis. Here, we test the oogenesis hypothesis by analyzing the data from 6 studies that described the number of pseudogenes on different chromosomes of the human and/or mouse genomes. Our results show a significant overabundance of processed pseudogenes in the X chromosomes and a significant underrepresentation of processed pseudogenes in the Y chromosome of the human genome. These observations support the hypothesis that the number of processed pseudogenes is proportional to the length of oogenesis.

Processed pseudogenes, processed genes, and spontaneous mutations in the Arabidopsis genome

We identified 411 processed sequences in the Arabidopsis thaliana genome based on the fact that they have lost their intron(s) and have a length that is at least 95% of the length of the gene that gave rise to them. These sequences were generated by 230 different genes and clearly originated from retrotranspositons events because most of them (91%) have a poly(A)-tail. They are composed of 376 sequences with frame shifts and/or premature stop codons (processed pseudogenes) and 35 sequences without disablements (processed genes). Eleven of these processed genes are likely functional retrotransposed genes because they have low Ka/Ks ratios and high Ks values, and their sequences match numerous Arabidopsis ESTs. Processed sequences are mostly randomly distributed in the Arabidopsis genome and their rate of accumulation has steadily been decreasing since it peaked some 50 MYA. In contrast with the situation observed in mammals, the processed sequences found in the Arabidopsis genome originate from genes with high copy numbers and not from highly expressed genes. The patterns of spontaneous mutations in Arabidopsis are slightly different than those of mammals but are similar to those observed in Drosophila. This suggests that methylated cytosine deamination is less frequent in Arabidopsis than in mammals.

Positive and negative selection in the beta-esterase gene cluster of the Drosophila melanogaster subgroup

We examine the pattern of molecular evolution of the beta-esterase gene cluster, including the Est-6 and psiEst-6 genes, in eight species of the Drosophila melanogaster subgroup. Using maximum likelihood estimates of nonsynonymous/synonymous rate ratios, we show that the majority of Est-6 sites evolves under strong (48% of sites) or moderate (50% of sites) negative selection and a minority of sites (1.5%) is under significant positive selection. Est-6 sites likely to be under positive selection are associated with increased intraspecific variability. One positively selected site is responsible for the EST-6 F/S allozyme polymorphism; the same site is responsible for the EST-6 functional divergence between species of the melanogaster subgroup. For psiEst-6 83.7% sites evolve under negative selection, 16% sites evolve neutrally, and 0.3% sites are under positive selection. The positively selected sites of psiEst-6 are located at the beginning and at the end of the gene, where there is reduced divergence between D. melanogaster and D. simulans; these regions of psiEst-6 could be involved in regulation or some other function. Branch-site-specific analysis shows that the evolution of the melanogaster subgroup underwent episodic positive selection. Collating the present data with previous results for the beta-esterase genes, we propose that positive and negative selection are involved in a complex relationship that may be typical of the divergence of duplicate genes as one or both duplicates evolve a new function.

Reverse transcription in genome evolution

As reverse transcription is predominantly performed by retrotransposable elements, the process is often entirely associated with the propagation of these elements. However, as a unique tool for transmitting information from the dynamic RNA to the more inert DNA, reverse transcription has been instrumental in shaping extant genomes. This review aims at presenting the diversity by which reverse transcription has influenced modern genomic structures.

mRNA surveillance of expressed pseudogenes in C. elegans

Messenger RNAs (mRNAs) that contain premature translation termination codons (PTCs) are targeted for rapid degradation in all eukaryotes tested. The mechanisms of nonsense-mediated mRNA decay (NMD) have been described in considerable detail, but the biological roles of NMD in wild-type organisms are poorly understood. mRNAs of wild-type organisms known to be degraded by NMD ("natural targets" of NMD) include by-products of regulated alternative splicing, out-of-frame mRNAs derived from unproductive gene rearrangements, cytoplasmic pre-mRNAs, endogenous retroviral and transposon RNAs, and mRNAs having upstream open reading frames or other unusual sequence features. NMD may function to eliminate aberrant PTC-containing mRNAs in order to protect cells from expression of potentially deleterious truncated proteins. Pseudogenes are nonfunctional genes or gene fragments that accumulate mutations through genetic drift. Such mutations will often introduce shifts of reading frame and/or PTCs, and mRNAs of expressed pseudogenes may thus be substrates of NMD. We demonstrate that mRNAs expressed from C. elegans pseudogenes are degraded by NMD and discuss possible implications for both mRNA surveillance and protein evolution. We describe an expressed pseudogene that encodes a small nucleolar RNA (snoRNA) within an intron and suggest this represents an evolutionary intermediate between snoRNA-encoding host genes that do or do not encode proteins.

Statistical alignment of retropseudogenes and their functional paralogs

We describe a model for the sequence evolution of a processed pseudogene and its paralog from a common protein-coding ancestor. The model accounts for substitutions, insertions, and deletions and combines nucleotide- and codon-level mutation models. We give a dynamic programming method for calculating the likelihood of homology between two sequences in the model and describe the accompanying alignment algorithm. We also describe how ancestral codons can be computed when the same gene produced multiple pseudogene homologs. We apply our methods to the evolution of human cytochrome c.

The innate allergenicity of helminth parasites

Helminth parasites are well known to induce an immune response in their hosts characterised by elevated IgE, peripheral blood or local tissue eosinophilia, and in some cases, intestinal mastocytosis. This immunological response has a strong T-helper 2 (Th2) cytokine bias and is reminiscent of the immunological constellation found in allergic diseases. However, the molecular forces driving the Th2 response to helminth parasites are still not understood. By using the human hookworm parasite Necator americanus as an example, the authors of the current article propose that in the course of its life cycle, this parasite becomes innately allergenic through the secretion of a molecular array designed to promote tissue migration and homing, feeding and survival against immunological attack. This complex array comprises proteases, lectins and other classes of molecules. Subsequent immunological and physiological events seemingly protect the host from both the allergic sequelae of exposure to environmental allergens and, moreover, from the parasite itself.

Computational identification of 69 retroposons in Arabidopsis

Retroposition is a shot-gun strategy of the genome to achieve evolutionary diversities by mixing and matching coding sequences with novel regulatory elements. We have identified 69 retroposons in the Arabidopsis (Arabidopsis thaliana) genome by a computational approach. Most of them were derivatives of mature mRNAs, and 20 genes contained relics of the reverse transcription process, such as truncations, deletions, and extra sequence additions. Of them, 22 are processed pseudogenes, and 52 genes are likely to be actively transcribed, especially in tissues from apical meristems (roots and flowers). Functional compositions of these retroposon parental genes imply that not the mRNA itself but its expression in gamete cells defines a suitable template for retroposition. The presence/absence patterns of retroposons can be used as cladistic markers for biogeographic research. Effects of human and the Mediterranean Pleistocene refugia in Arabidopsis biogeographic distributions were revealed based on two recent retroposons (At1g61410 and At5g52090). An evolutionary rate of new gene creation by retroposition was calculated as 0.6 genes per million years. Retroposons can also be used as molecular fossils of the parental gene expressions in ancient time. Extensions of 3' untranslated regions for those expressed parental genes are revealed as a possible trend of plant transcriptome evolution. In addition, we reported the first plant functional chimeric gene that adapts to intercompartmental transport by capturing two additional exons after retroposition.

Identification of a novel protein D3UPCA from Halobacterium salinarum and prediction of its function

Halobacterium salinarum is an extremely halophilic archaea, which is able to live in highly saline environments. In a recent study, several halophilic archaea were found to have the ability to biodegrade organic hydrocarbon pollutants, but protein information regarding hydrocarbon degradation and tolerance in halophilic archaea has been relatively rare. In this study, the protein expression profile of H. salinarum cultured under different diesel concentrations (0, 2 and 4%) was investigated by two-dimensional gel electrophoresis. Proteins which increased their expression levels in diesel media were identified using matrix-assisted laser desorption/ionization-time of flight and electrospray ionization-tandem mass spectrometry analysis. Among these, a protein spot (named D3UPCA) which was up-regulated about nine-fold and found to have COG3388, an uncharacterized protein conserved in archaea, was selected in order to further characterize its functions. The D3UPCA coding gene (named d3upca) was cloned and expressed in Escherichia coli, and purified by the glutathione-S-transferase-fusion method. The function of the protein was estimated using various bioinformatics tools and was predicted to be related to the regulation of transcription and/or translation of genes needed to tolerate stresses associated within the presence of diesel oil.

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.

Expression analysis of ABC transporters reveals differential functions of tandemly duplicated genes in Caenorhabditis elegans

We have previously identified 60 predicted ABC transporter genes in the Caenorhabditis elegans genome and classified them into eight groups. As an initial step towards understanding how these putative ABC genes work in worms, we generated promoter-fluorescent protein fusions for the entire family to address when and where these genes are turned on in vivo. Both Aequoria green fluorescent protein (GFP) and Discosoma red fluorescent protein (RFP) were used as reporters in our transgenic assay. Observable expression is more frequently seen from fusions to genes in subfamilies B, C, D and E than those in subfamilies A and G. Sixteen worm ABC genes are found in tandem duplications, forming two four-gene clusters and four two-gene clusters. Fifteen out of the 16 duplicated gene promoters drove different or partially overlapping expression patterns, suggesting active functions for these duplicated genes. Furthermore, our results suggest that an internal promoter can cause differential expression of genes within an operon. Finally, our observations suggest that it is possible for coding sequences to function as a regulatory region for a neighbouring gene.

Bovine ?-defensins: Identification and characterization of novel bovine ?-defensin genes and their expression in mammarygland tissue

beta-Defensin genes code for multifunctional peptides with a broad-range antimicrobial activity. In this project we hypothesized that beta-defensin genes may be candidate genes for resistance to mastitis. In this article we describe the identification and genomic characterization of eight bovine beta-defensin genes, including six novel defensin genes and two pseudogenes. Expression in the bovine mammary gland of one of the novel genes, DEFB401, has been demonstrated, as well as the expression of LAP, TAP, DEFB1, BNBD3, BNBD9, and BNBD12. For genomic characterization, 20 BACs from two different bovine BAC libraries (RZPD numbers 750 and 754) were isolated by PCR screening with beta-defensin consensus primers derived from published sequences. PCR products from BACs generated with consensus primers have been subcloned and sequenced, revealing a total of 16 genes and two pseudogenes. Six novel beta-defensin genes share the typical exon-intron structure and are highly homologous to published bovine beta-defensin genes. They are named DEFB401- DEFB405 and LAP-like, and two novel pseudogenes are named EBD-P and EBD-P2. Analysis of mammary gland tissue-derived cDNA from nine cows with different clinical findings demonstrated the expression of several beta-defensin genes mentioned above. First results indicate that the lactational status of the cow presumably has no influence on gene expression. Competent knowledge of antimicrobial activity of beta-defensins from literature, the abundance of beta-defensin mRNA in the bovine mammary gland, and the inducibility of some genes give first evidence that beta-defensins may play a role in local host defense during udder infections.

Retention of enzyme gene duplicates by subfunctionalization

Duplication-degeneration-complementation (DDC) describes a process by which evolving duplicates of a pleiotropic ancestral gene divide up the multiple functions of the ancestor between them (i.e. subfunctionalize), and this ultimately frustrates the rate of pseudogene formation. Focusing explicitly on enzyme-like pleiotropic function, we model DDC driven by sequence divergence between duplicates. The model incorporates an idealized sequence-function mapping in which enzyme-substrate binding affinity is related to hydrophobic versus polar (HP) amino-acid composition of tertiary structure about the binding pocket. In this sense, a transparent coupling between physical-chemical function of an enzyme and sequence evolution is presented.

Comprehensive analysis of pseudogenes in prokaryotes: widespread gene decay and failure of putative horizontally transferred genes

A comprehensive analysis of the occurrence of pseudogenes in a diverse selection of 64 prokaryote genomes identified around 7,000 candidate pseudogenes. A large fraction of prokaryote pseudogenes seems to have arisen from failed horizontal-transfer events.

Background

Pseudogenes often manifest themselves as disabled copies of known genes. In prokaryotes, it was generally believed (with a few well-known exceptions) that they were rare.

Results

We have carried out a comprehensive analysis of the occurrence of pseudogenes in a diverse selection of 64 prokaryote genomes. Overall, we find a total of around 7,000 candidate pseudogenes. Moreover, in all the genomes surveyed, pseudogenes occur in at least 1 to 5% of all gene-like sequences, with some genomes having considerably higher occurrence. Although many large populations of pseudogenes arise from large, diverse protein families (for example, the ABC transporters), notable numbers of pseudogenes are associated with specific families that do not occur that widely. These include the cytochrome P450 and PPE families (PF00067 and PF00823) and others that have a direct role in DNA transposition.

Conclusions

We find suggestive evidence that a large fraction of prokaryote pseudogenes arose from failed horizontal transfer events. In particular, we find that pseudogenes are more than twice as likely as genes to have anomalous codon usage associated with horizontal transfer. Moreover, we found a significant difference in the number of horizontally transferred pseudogenes in pathogenic and non-pathogenic strains of Escherichia coli.

A new role for expressed pseudogenes as ncRNA: regulation of mRNA stability of its homologous coding gene

We have earlier generated a mutant mouse in a course of making a transgenic line that exhibited interesting heterozygote phenotypes, which exhibited failure to thrive, severe bone deformities, and polycystic kidneys. This mutant mouse provided a clue to uncover a unique role of expressed pseudogenes. In this mutant the transgene was integrated into the vicinity of the expressing pseudogene of Makorin1 called Makorin1-p1. This insertion reduced transcription of the Makorin1-p1, resulting in destabilization of the Makorin1 mRNA in trans via a cis-acting RNA decay element within the 5' region of Makorin1 that is homologous between Makorin1 and Makorin1-p1. These findings demonstrate a novel and specific regulatory role of an expressed pseudogene as well as functional significance for noncoding RNAs. Next, we developed an original algorithm to determine how many pseudogenes are expressed. Based on our examination 2-3% of human processed pseudogenes are expressed using the most strict criteria. Interestingly, the mouse has a much smaller proportion of expressed pseudogenes (0.5-1%). Pseudogenes are functionally less constrained, and have accumulated more mutations than translated genes. If they have some functions in gene regulation, this property would allow more rapid functional diversification than protein-coding genes. In addition, some genetic phenomena that exhibit incomplete penetrance might be attributed to "mutation" or "variation" of pseudogenes.