Evolutionary origin of regulatory regions of retrogenes in Drosophila

Background

Retrogenes are processed copies of other genes. This duplication mechanism produces a copy of the parental gene that should not contain introns, and usually does not contain cis-regulatory regions. Here, we computationally address the evolutionary origin of promoter and other cis-regulatory regions in retrogenes using a total of 94 Drosophila retroposition events we recently identified. Previous tissue expression data has revealed that a large fraction of these retrogenes are specifically and/or highly expressed in adult testes of Drosophila.

Results

In this work, we infer that retrogenes do not generally carry regulatory regions from aberrant upstream or normal transcripts of their parental genes, and that expression patterns of neighboring genes are not consistently shared by retrogenes. Additionally, transposable elements do not appear to substantially provide regulatory regions to retrogenes. Interestingly, we find that there is an excess of retrogenes in male testis neighborhoods that is not explained by insertional biases of the retroelement machinery used for retroposition.

Conclusion

We conclude that retrogenes' regulatory regions mostly do not represent a random set of existing regulatory regions. On the contrary, our conclusion is that selection is likely to have played an important role in the persistence of autosomal testis biased retrogenes. Selection in favor of retrogenes inserted in male testis neighborhoods and at the sequence level to produce testis expression is postulated to have occurred.

In vivo analysis of developmentally and evolutionarily dynamic protein-DNA interactions regulating transcription of the Pgk2 gene during mammalian spermatogenesis

Transcription of the testis-specific Pgk2 gene is selectively activated in primary spermatocytes to provide a source of phosphoglycerate kinase that is critical to normal motility and fertility of mammalian spermatozoa. We examined dynamic changes in protein-DNA interactions at the Pgk2 gene promoter during murine spermatogenesis in vivo by performing genomic footprinting and chromatin immunoprecipitation assays with enriched populations of murine spermatogenic cells at stages prior to, during, and following transcription of this gene. We found that genes encoding the testis-specific homeodomain factor PBX4 and its coactivator, PREP1, are expressed in patterns that mirror expression of the Pgk2 gene and that these factors become bound to the Pgk2 enhancer in cells in which this gene is actively expressed. We therefore suggest that these factors, along with CREM and SP3, direct stage- and cell type-specific transcription of the Pgk2 gene during spermatogenesis. We propose that binding of PBX4, plus its coactivator PREP1, is a rate-limiting step leading to the initiation of tissue-specific transcription of the Pgk2 gene. This study provides insight into the developmentally dynamic establishment of tissue-specific protein-DNA interactions in vivo. It also allows us to speculate about the events that led to tissue-specific regulation of the Pgk2 gene during mammalian evolution.

Epigenetic regulation of testis-specific gene expression

Cellular differentiation is mediated by differential gene expression. The cells of the testis are no exception. Indeed, recent studies based on microarray and expressed sequence tag analyses have revealed dynamic changes in gene expression during spermatogenesis. The autosomal phosphoglycerate kinase gene Pgk2 is an example of a gene that is tightly regulated during spermatogenesis, with transcription initiating in primary spermatocytes and ceasing in postmeiotic spermatids. Our studies show that this tissue-specific transcription is regulated at multiple levels, including binding of ubiquitous and testis-specific transcription factors, reconfiguration of nucleosomes, decondensation of chromatin, modifications of specific residues in histones, and demethylation of DNA. Our ability to isolate relatively pure populations of different spermatogenic cell types has allowed us to determine the order of these events that lead up to initiation of transcription of the Pgk2 gene. By examining similar parameters in transgenic mice carrying various portions of Pgk2 regulatory sequence, we have localized specific signals responsible for each of these regulatory events. Together these results demonstrate that activation of testis-specific transcription of the Pgk2 gene involves a multifaceted, ordered cascade of epigenetic events that potentiates this locus in preparation for the initiation of transcription.

Extensive gene traffic on the mammalian X chromosome

Mammalian sex chromosomes have undergone profound changes since evolving from ancestral autosomes. By examining retroposed genes in the human and mouse genomes, we demonstrate that, during evolution, the mammalian X chromosome has generated and recruited a disproportionately high number of functional retroposed genes, whereas the autosomes experienced lower gene turnover. Most autosomal copies originating from X-linked genes exhibited testis-biased expression. Such export is incompatible with mutational bias and is likely driven by natural selection to attain male germline function. However, the excess recruitment is consistent with a combination of both natural selection and mutational bias.

The origin of new genes: glimpses from the young and old

Genome data have revealed great variation in the numbers of genes in different organisms, which indicates that there is a fundamental process of genome evolution: the origin of new genes. However, there has been little opportunity to explore how genes with new functions originate and evolve. The study of ancient genes has highlighted the antiquity and general importance of some mechanisms of gene origination, and recent observations of young genes at early stages in their evolution have unveiled unexpected molecular and evolutionary processes.

Retroposon compensatory mechanism hypothesis not supported: Zfa knockout mice are fertile

It is hypothesized that autosomal retroposons compensate for the loss of their inactivated essential X-chromosome progenitors during spermatogenesis. Here we test this Retroposon Compensatory Mechanism (RCM) hypothesis using the Zfy gene family. The mouse autosomal retroposon Zfa is expressed in testes at the same developmental time points at which Zfx levels decline, which correspond to the time of male sex chromosome inactivation, suggesting that Zfa may compensate for the loss of Zfx during spermatogenesis. We examined the effect of Zfa-targeted mutagenesis on spermatogenesis in three genetically distinct mouse strains. Surprisingly, Zfa knockout mice showed no detectable fertility, sperm count, or testes morphology defects. We therefore conclude that Zfa is not an essential gene for spermatogenesis and fertility. This surprising finding now challenges the RCM hypothesis at least for the Zfy gene family. It also forces us to reevaluate the original data underpinning the RCM hypothesis for this family and to propose alternative hypotheses.

A third isoform of cytochrome c oxidase subunit VIII is present in mammals

The terminal enzyme of the mitochondrial respiratory chain, cytochrome c oxidase (COX), contains three mitochondrial and ten nuclear encoded subunits in mammals. Three of the nuclear subunits (VIa, VIIa, and VIII) have muscle and non-muscle-specific isoforms, subunit IV contains a lung-specific isoform, and subunit VIb contains a testes-specific isoform. For subunit VIII, the smallest nuclear encoded COX polypeptide, we have now found a third gene (COX 8-3), which has been identified in human, lemur, rat, and mouse, suggesting that it is present in a broad range of Eutherian mammals. Sequence similarity and gene structure support the homology of COX8-3 to the other subunit VIII isoforms, indicating that all three are the product of gene duplications. COX VIII-3 protein is mitochondrially-targeted, as shown by a fluorescent COX VIII3/DsRed fusion protein. Both the mitochondrial targeting and its sequence conservation suggest that COXVIII-3 functions as part of the COX holoenzyme and could have a tissue-specific role, as is the case for the other two isoforms. Questions remain about where COX8-3 is predominantly expressed. However, detection of full-length cDNAs, lower levels of sequence divergence at the first and second codon positions compared to the third, and a conserved gene structure indicate that COX VIII-3 is an expressed gene whose origin dates to at least 91 million years ago.

Retroposed New Genes Out of the X in Drosophila

New genes that originated by various molecular mechanisms are an essential component in understanding the evolution of genetic systems. We investigated the pattern of origin of the genes created by retroposition in Drosophila. We surveyed the wholeDrosophila melanogaster genome for such new retrogenes and experimentally analyzed their functionality and evolutionary process. These retrogenes, functional as revealed by the analysis of expression, substitution, and population genetics, show a surprisingly asymmetric pattern in their origin. There is a significant excess of retrogenes that originate from the X chromosome and retropose to autosomes; new genes retroposed from autosomes are scarce. Further, we found that most of these X-derived autosomal retrogenes had evolved a testis expression pattern. These observations may be explained by natural selection favoring those new retrogenes that moved to autosomes and avoided the spermatogenesis X inactivation, and suggest the important role of genome position for the origin of new genes.

[The sequence data from this study have been submitted to GenBank under accession nos. AY150701–AY150797. The following individuals kindly provided reagents, samples, or unpublished information as indicated in the paper: M.-L. Wu, F. Lemeunier, and P. Gibert.]

Analysis of the cDNA and encoded protein of the human testis-specific PGK-2 gene

Because of their unique function, germ cells require unique gene products. Thus, although the glycolytic enzyme phosphoglycerate kinase (PGK) is required in all metabolically active cell types, there are two functional PGK genes in the mammalian genome, one, PGK-1, that is X-linked and ubiquitously expressed in all somatic tissues, and a second, PGK-2, that is autosomal and expressed only in spermatogenic cells. Expression of the PGK-2 gene may function solely to compensate for repressed expression of the PGK-1 gene due to X-chromosome inactivation in spermatocytes. Alternative y, the PGK-2 gene could encode an isozyme with unique characteristics that are beneficial to spermatozoa. We have isolated a cDNA of the human PGK-2 gene and used this as probe to demonstrate that transcription of this gene in spermatocytes and spermatids coincides with a period of repressed transcription of the X-linked PGK-1 gene during spermatogenesis in the human testis. We have also analyzed the amino acid sequence and protein characteristics of the PGK-2 isozyme deduced from this cDNA and compared them with that of the human PGK-1 isozyme to show that known structural and functional motifs are conserved in both proteins. Finally, we have examined the distribution of the PGK-1 and PGK-2 isozymes during spermatogenesis in the mouse to show that while the PGK-2 protein does not appear to possess any unique intracellular localization signal, it is more stable in vivo than the PGK-1 protein.

The promoter of the Poly(A) binding protein 2 (Pabp2) retroposon is derived from the 5'-untranslated region of the Pabp1 progenitor gene

The mouse Pabp2 retroposon encodes an isoform of poly(A) binding protein that is expressed in meiotic and early haploid spermatogenic cells. In the present study, we have determined the transcription start site of the Pabp2 gene to clarify the source of its promoter, a prerequisite for expression of retroposons and preservation of their function by natural selection. The 5' end of the mouse Pabp2 retroposon exhibits extensive similarity to the entire 5' UTR of the human PABP1 mRNA, but there is no similarity upstream of the transcription start site of the human PABP1 mRNA, indicating that the Pabp2 gene lacks 5' flanking sequences of the parental PABP1 gene. Oligonucleotide-directed RNase H cleavage and 5' rapid amplification of cDNA ends both indicate that the transcription start site of the mouse Pabp2 gene is located approximately 330 bases downstream of the capsite of the PABP1 mRNA, indicating that the Pabp2 promoter is derived from the PABP1 5' UTR.

Testis-specific murine centrin, Cetn1: genomic characterization and evidence for retroposition of a gene encoding a centrosome protein

Centrin is a centrosome component in species from yeast to humans. Here, the mouse centrin 1 gene (Cetn1) is analyzed with respect to its genomic structure, chromosome localization, tissue-specific expression, and phylogenetic relationship to the other mouse centrin genes and their human orthologs. Cetn1 is an intronless gene located on chromosome 18A2 that encodes a 172-amino-acid protein with a predicted molecular mass of 19,696 Da (pI 4.61) and all of the structural features common to centrin. Cetn1 possesses the sequence features of an expressed retroposon: the gene lacks introns, the open reading frame is not interrupted by stop codons, and the coding region is flanked by a pair of direct repeats. Reverse transcriptase-polymerase chain reaction and Northern blot analysis demonstrate that Cetn1 expression is limited exclusively to the testis in adult male mice. Cetn1 expression is first seen in the neonatal testis at 14 days postpartum, reaching adult levels by day 17. These observations provide new insight into the regulation, function, and evolutionary history of centrin in higher eukaryotes.

Multiple elements influence transcriptional regulation from the human testis-specific PGK2 promoter in transgenic mice

The PGK2 gene is expressed in a strictly tissue-specific manner in meiotic spermatocytes and postmeiotic spermatids during spermatogenesis in eutherian mammals. Previous results indicate that this is regulated at the transcriptional level by core promoter sequences that bind ubiquitous transcription factors and by sequences in a 40-base pair (bp) upstream enhancer region (E1/E4) that bind tissue-specific transcription factors. Transgenic mice carrying different PGK2 promoter sequences linked to the chloramphenicol acetyltransferase (CAT) reporter gene, one containing only the 40-bp E1/E4 enhancer sequence plus the core promoter and two containing 515 bp of PGK2 promoter but with either the E1/E4 enhancer region or the Sp1-binding site in the core promoter disrupted by in vitro mutagenesis, all showed levels of expression reduced to less than half that of the wild-type 515 PGK2/CAT transgene. These results indicate that multiple factor-binding regions normally regulate initiation of transcription from the PGK2 promoter. The single disruption of any one of these binding activities reduced, but did not abolish, transgene expression. This is consistent with an "enhanceosome"-like function in this promoter involving multiple bound activator proteins that interact in a combinatorial manner to synergistically promote testis-specific transcription.