Genome of the marsupial Monodelphis domestica reveals innovation in non-coding sequences

The evolution of the DLK1-DIO3 imprinted domain in mammals

A comprehensive, domain-wide comparative analysis of genomic imprinting between mammals that imprint and those that do not can provide valuable information about how and why imprinting evolved. The imprinting status, DNA methylation, and genomic landscape of the Dlk1-Dio3 cluster were determined in eutherian, metatherian, and prototherian mammals including tammar wallaby and platypus. Imprinting across the whole domain evolved after the divergence of eutherian from marsupial mammals and in eutherians is under strong purifying selection. The marsupial locus at 1.6 megabases, is double that of eutherians due to the accumulation of LINE repeats. Comparative sequence analysis of the domain in seven vertebrates determined evolutionary conserved regions common to particular sub-groups and to all vertebrates. The emergence of Dlk1-Dio3 imprinting in eutherians has occurred on the maternally inherited chromosome and is associated with region-specific resistance to expansion by repetitive elements and the local introduction of noncoding transcripts including microRNAs and C/D small nucleolar RNAs. A recent mammal-specific retrotransposition event led to the formation of a completely new gene only in the eutherian domain, which may have driven imprinting at the cluster.

Mammals have two copies of each gene in their somatic cells, and most of these gene pairs are regulated and expressed simultaneously. A fraction of mammalian genes, however, is subject to imprinting—a chemical modification that marks a gene according to its parental origin, so that one parent's copy is expressed while the other parent's copy is silenced. How and why this process evolved is the subject of much speculation. Here we have shown that all the genes in one genomic region, Dlk1-Dio3, which are imprinted in placental mammals such as mouse and human, are not imprinted in marsupial (wallaby) or monotreme (platypus) mammals. This is in contrast to a small number of other imprinted genes that are imprinted in marsupials and other therian mammals and indicates that imprinting arose at each genomic domain at different stages of mammalian evolution. We have compared the sequence of the Dlk1-Dio3 region between seven vertebrate species and identified sequences that are differentially represented in mammals that imprint compared to those that do not. Our data indicate that once imprinted gene regulation is acquired in a domain, it becomes evolutionarily constrained to remain unchanged.

A comparative analysis of genomic imprinting between mammals that imprint and those that don't has provided insights into how and why imprinting evolved.

A novel mode of enhancer evolution: the Tal1 stem cell enhancer recruited a MIR element to specifically boost its activity

Altered cis-regulation is thought to underpin much of metazoan evolution, yet the underlying mechanisms remain largely obscure. The stem cell leukemia TAL1 (also known as SCL) transcription factor is essential for the normal development of blood stem cells and we have previously shown that the Tal1 +19 enhancer directs expression to hematopoietic stem cells, hematopoietic progenitors, and to endothelium. Here we demonstrate that an adjacent region 1 kb upstream (+18 element) is in an open chromatin configuration and carries active histone marks but does not function as an enhancer in transgenic mice. Instead, it boosts activity of the +19 enhancer both in stable transfection assays and during differentiation of embryonic stem (ES) cells carrying single-copy reporter constructs targeted to the Hprt locus. The +18 element contains a mammalian interspersed repeat (MIR) which is essential for the +18 function and which was transposed to the Tal1 locus approximately 160 million years ago at the time of the mammalian/marsupial branchpoint. Our data demonstrate a previously unrecognized mechanism whereby enhancer activity is modulated by a transposon exerting a "booster" function which would go undetected by conventional transgenic approaches.

The status of dosage compensation in the multiple X chromosomes of the platypus

Dosage compensation has been thought to be a ubiquitous property of sex chromosomes that are represented differently in males and females. The expression of most X-borne genes is equalized between XX females and XY males in therian mammals (marsupials and “placentals”) by inactivating one X chromosome in female somatic cells. However, compensation seems not to be strictly required to equalize the expression of most Z-borne genes between ZZ male and ZW female birds. Whether dosage compensation operates in the third mammal lineage, the egg-laying monotremes, is of considerable interest, since the platypus has a complex sex chromosome system in which five X and five Y chromosomes share considerable genetic homology with the chicken ZW sex chromosome pair, but not with therian XY chromosomes. The assignment of genes to four platypus X chromosomes allowed us to examine X dosage compensation in this unique species. Quantitative PCR showed a range of compensation, but SNP analysis of several X-borne genes showed that both alleles are transcribed in a heterozygous female. Transcription of 14 BACs representing 19 X-borne genes was examined by RNA-FISH in female and male fibroblasts. An autosomal control gene was expressed from both alleles in nearly all nuclei, and four pseudoautosomal BACs were usually expressed from both alleles in male as well as female nuclei, showing that their Y loci are active. However, nine X-specific BACs were usually transcribed from only one allele. This suggests that while some genes on the platypus X are not dosage compensated, other genes do show some form of compensation via stochastic transcriptional inhibition, perhaps representing an ancestral system that evolved to be more tightly controlled in placental mammals such as human and mouse.

Dosage compensation equalizes the expression of genes found on sex chromosomes so that they are equally expressed in females and males. In placental and marsupial mammals, this is accomplished by silencing one of the two X chromosomes in female cells. In birds, dosage compensation seems not to be strictly required to balance the expression of most genes on the Z chromosome between ZZ males and ZW females. Whether dosage compensation exists in the third group of mammals, the egg-laying monotremes, is of considerable interest, particularly since the platypus has five different X and five different Y chromosomes. As part of the platypus genome project, genes have now been assigned to four of the five X chromosomes. We have shown that there is some evidence for dosage compensation, but it is variable between genes. Most interesting are our results showing that there is a difference in the probability of expression for X-specific genes, with about 50% of female cells having two active copies of an X gene while the remainder have only one. This means that, although the platypus has the variable compensation characteristic of birds, it also has some level of inactivation, which is characteristic of dosage compensation in other mammals.

Delineating slowly and rapidly evolving fractions of the Drosophila genome

Evolutionary conservation is an important indicator of function and a major component of bioinformatic methods to identify non-protein-coding genes. We present a new Bayesian method for segmenting pairwise alignments of eukaryotic genomes while simultaneously classifying segments into slowly and rapidly evolving fractions. We also describe an information criterion similar to the Akaike Information Criterion (AIC) for determining the number of classes. Working with pairwise alignments enables detection of differences in conservation patterns among closely related species. We analyzed three whole-genome and three partial-genome pairwise alignments among eight Drosophila species. Three distinct classes of conservation level were detected. Sequences comprising the most slowly evolving component were consistent across a range of species pairs, and constituted approximately 62-66% of the D. melanogaster genome. Almost all (>90%) of the aligned protein-coding sequence is in this fraction, suggesting much of it (comprising the majority of the Drosophila genome, including approximately 56% of non-protein-coding sequences) is functional. The size and content of the most rapidly evolving component was species dependent, and varied from 1.6% to 4.8%. This fraction is also enriched for protein-coding sequence (while containing significant amounts of non-protein-coding sequence), suggesting it is under positive selection. We also classified segments according to conservation and GC content simultaneously. This analysis identified numerous sub-classes of those identified on the basis of conservation alone, but was nevertheless consistent with that classification. Software, data, and results available at www.maths.qut.edu.au/-keithj/. Genomic segments comprising the conservation classes available in BED format.

Cone visual pigments of monotremes: Filling the phylogenetic gap

We have determined the sequence and genomic organization of the genes encoding the cone visual pigment of the platypus (Ornithorhynchus anatinus) and the echidna (Tachyglossus aculeatus), and inferred their spectral properties and evolutionary pathways. We prepared platypus and echidna retinal RNA and used primers of the middle-wave-sensitive (MWS), long-wave-sensitive (LWS), and short-wave sensitive (SWS1) pigments corresponding to coding sequences that are highly conserved among mammals; to PCR amplify the corresponding pigment sequences. Amplification from the retinal RNA revealed the expression of LWS pigment mRNA that is homologous in sequence and spectral properties to the primate LWS visual pigments. However, we were unable to amplify the mammalian SWS1 pigment from these two species, indicating this gene was lost prior to the echidna-platypus divergence (∼21 MYA). Subsequently, when the platypus genome sequence became available, we found an LWS pigment gene in a conserved genomic arrangement that resembles the primate pigment, but, surprisingly we found an adjacent (∼20 kb) SWS2 pigment gene within this conserved genomic arrangement. We obtained the same result after sequencing the echidna genes. The encoded SWS2 pigment is predicted to have a wavelength of maximal absorption of about 440 nm, and is paralogous to SWS pigments typically found in reptiles, birds, and fish but not in mammals. This study suggests the locus control region (LCR) has played an important role in the conservation of photo receptor gene arrays and the control of their spatial and temporal expression in the retina in all mammals. In conclusion, a duplication event of an ancestral cone visual pigment gene, followed by sequence divergence and selection gave rise to the LWS and SWS2 visual pigments. So far, the echidna and platypus are the only mammals that share the gene structure of the LWS-SWS2 pigment gene complex with reptiles, birds and fishes.

On the extent and origins of genic novelty in the phylum Nematoda

BACKGROUND

The phylum Nematoda is biologically diverse, including parasites of plants and animals as well as free-living taxa. Underpinning this diversity will be commensurate diversity in expressed genes, including gene sets associated specifically with evolution of parasitism.

METHODS AND FINDINGS

Here we have analyzed the extensive expressed sequence tag data (available for 37 nematode species, most of which are parasites) and define over 120,000 distinct putative genes from which we have derived robust protein translations. Combined with the complete proteomes of Caenorhabditis elegans and Caenorhabditis briggsae, these proteins have been grouped into 65,000 protein families that in turn contain 40,000 distinct protein domains. We have mapped the occurrence of domains and families across the Nematoda and compared the nematode data to that available for other phyla. Gene loss is common, and in particular we identify nearly 5,000 genes that may have been lost from the lineage leading to the model nematode C. elegans. We find a preponderance of novelty, including 56,000 nematode-restricted protein families and 26,000 nematode-restricted domains. Mapping of the latest time-of-origin of these new families and domains across the nematode phylogeny revealed ongoing evolution of novelty. A number of genes from parasitic species had signatures of horizontal transfer from their host organisms, and parasitic species had a greater proportion of novel, secreted proteins than did free-living ones.

CONCLUSIONS

These classes of genes may underpin parasitic phenotypes, and thus may be targets for development of effective control measures.

SCPP gene evolution and the dental mineralization continuum

Many genes critical to vertebrate skeletal mineralization are members of the secretory calcium-binding phosphoprotein (SCPP) gene family, which has evolved by gene duplication from a single ancestral gene. In humans, mutations in some of these SCPP genes have been associated with various diseases related to dentin or enamel hypoplasia. Recently, systematic searches for SCPP genes of various species have allowed us to investigate the history of phylogenetically variable dental tissues as a whole. One important conclusion is that not all disease-associated SCPP genes are present in tetrapods, and teleost fish probably have none, even in toothed species, having acquired their complement of SCPP genes through an independent duplication history. Here, we review comparative analyses of mineralized dental tissues, with particular emphasis on the use of SCPPs, within and between tetrapods and teleosts. Current knowledge suggests a close relationship among bone, dentin, teleost fish enameloid (enamel-like hard tissue), and tetrapod enamel. These tissues thus form a mineralized-tissue continuum. Contemporary dental tissues have evolved from an ancestral continuum through lineage-specific modifications.

Assignment of isochores for all completely sequenced vertebrate genomes using a consensus

A new consensus isochore assignment method and a database of isochore maps for all completely sequenced vertebrate genomes are presented.

We show that although the currently available isochore mapping methods agree on the isochore classification of about two-thirds of the human DNA, they produce significantly different results with regard to the location of isochore boundaries and isochore length distribution. We present a new consensus isochore assignment method based on majority voting and provide IsoBase, a comprehensive on-line database of isochore maps for all completely sequenced vertebrate genomes.

Conservation of the H19 noncoding RNA and H19-IGF2 imprinting mechanism in therians

Comparisons between eutherians and marsupials suggest limited conservation of the molecular mechanisms that control genomic imprinting in mammals. We have studied the evolution of the imprinted IGF2-H19 locus in therians. Although marsupial orthologs of protein-coding exons were easily identified, the use of evolutionarily conserved regions and low-stringency Bl2seq comparisons was required to delineate a candidate H19 noncoding RNA sequence. The therian H19 orthologs show miR-675 and exon structure conservation, suggesting functional selection on both features. Transcription start site sequences and poly(A) signals are also conserved. As in eutherians, marsupial H19 is maternally expressed and paternal methylation upstream of the gene originates in the male germline, encompasses a CTCF insulator, and spreads somatically into the H19 gene. The conservation in all therians of the mechanism controlling imprinting of the IGF2-H19 locus suggests a sequential model of imprinting evolution.

Optimization of volumetric computed tomography for skeletal analysis of model genetic organisms

Forward and reverse genetics now allow researchers to understand embryonic and postnatal gene function in a broad range of species. Although some genetic mutations cause obvious morphological change, other mutations can be more subtle and, without adequate observation and quantification, might be overlooked. For the increasing number of genetic model organisms examined by the growing field of phenomics, standardized but sensitive methods for quantitative analysis need to be incorporated into routine practice to effectively acquire and analyze ever-increasing quantities of phenotypic data. In this study, we present platform-independent parameters for the use of microscopic x-ray computed tomography (microCT) for phenotyping species-specific skeletal morphology of a variety of different genetic model organisms. We show that microCT is suitable for phenotypic characterization for prenatal and postnatal specimens across multiple species.

Pervasive positive selection on duplicated and nonduplicated vertebrate protein coding genes

A stringent branch-site codon model was used to detect positive selection in vertebrate evolution. We show that the test is robust to the large evolutionary distances involved. Positive selection was detected in 77% of 884 genes studied. Most positive selection concerns a few sites on a single branch of the phylogenetic tree: Between 0.9% and 4.7% of sites are affected by positive selection depending on the branches. No functional category was overrepresented among genes under positive selection. Surprisingly, whole genome duplication had no effect on the prevalence of positive selection, whether the fish-specific genome duplication or the two rounds at the origin of vertebrates. Thus positive selection has not been limited to a few gene classes, or to specific evolutionary events such as duplication, but has been pervasive during vertebrate evolution.

Chromosomal gene movements reflect the recent origin and biology of therian sex chromosomes

Mammalian sex chromosomes stem from ancestral autosomes and have substantially differentiated. It was shown that X-linked genes have generated duplicate intronless gene copies (retrogenes) on autosomes due to this differentiation. However, the precise driving forces for this out-of-X gene "movement" and its evolutionary onset are not known. Based on expression analyses of male germ-cell populations, we here substantiate and extend the hypothesis that autosomal retrogenes functionally compensate for the silencing of their X-linked housekeeping parental genes during, but also after, male meiotic sex chromosome inactivation (MSCI). Thus, sexually antagonistic forces have not played a major role for the selective fixation of X-derived gene copies in mammals. Our dating analyses reveal that although retrogenes were produced ever since the common mammalian ancestor, selectively driven retrogene export from the X only started later, on the placental mammal (eutherian) and marsupial (metatherian) lineages, respectively. Together, these observations suggest that chromosome-wide MSCI emerged close to the eutherian-marsupial split approximately 180 million years ago. Given that MSCI probably reflects the spread of the recombination barrier between the X and Y, crucial for their differentiation, our data imply that these chromosomes became more widely differentiated only late in the therian ancestor, well after the divergence of the monotreme lineage. Thus, our study also provides strong independent support for the recent notion that our sex chromosomes emerged, not in the common ancestor of all mammals, but rather in the therian ancestor, and therefore are much younger than previously thought.

Regulation of mammalian gene expression by retroelements and non-coding tandem repeats

Genomes of higher eukaryotes contain abundant non-coding repeated sequences whose overall biological impact is unclear. They comprise two categories. The first consists of retrotransposon-derived elements. These are three major families of retroelements (LINEs, SINEs and LTRs). SINEs are clustered in gene-rich regions and are found in promoters of genes while LINEs are concentrated in gene-poor regions and are depleted from promoters. The second class consists of non-coding tandem repeats (satellite DNAs and TTAGGG arrays), which are associated with mammalian centromeres, heterochromatin and telomeres. Terminal TTAGGG arrays are involved in telomere capping and satellite DNAs are located in heterochromatin, which is implicated in transcription silencing by gene repositioning (relocalization). It is unknown whether interstitial TTAGGG sequences, which are present in many vertebrates, have a function. Here, evidence will be presented that retroelements and TTAGGG arrays are involved in regulation of gene expression. Retroelements can provide binding sites for transcription factors and protect promoter CpG islands from repressive chromatin modifications, and may be also involved in nuclear compartmentalization of transcriptionally active and inactive domains. Interstitial telomere-like sequences can form dynamically maintained three-dimensional nuclear networks of transcriptionally inactive domains, which may be involved in transcription silencing like classic heterochromatin.

Adaptive evolution and functional redesign of core metabolic proteins in snakes

BACKGROUND

Adaptive evolutionary episodes in core metabolic proteins are uncommon, and are even more rarely linked to major macroevolutionary shifts.

METHODOLOGY/PRINCIPAL FINDINGS

We conducted extensive molecular evolutionary analyses on snake mitochondrial proteins and discovered multiple lines of evidence suggesting that the proteins at the core of aerobic metabolism in snakes have undergone remarkably large episodic bursts of adaptive change. We show that snake mitochondrial proteins experienced unprecedented levels of positive selection, coevolution, convergence, and reversion at functionally critical residues. We examined Cytochrome C oxidase subunit I (COI) in detail, and show that it experienced extensive modification of normally conserved residues involved in proton transport and delivery of electrons and oxygen. Thus, adaptive changes likely altered the flow of protons and other aspects of function in CO, thereby influencing fundamental characteristics of aerobic metabolism. We refer to these processes as "evolutionary redesign" because of the magnitude of the episodic bursts and the degree to which they affected core functional residues.

CONCLUSIONS/SIGNIFICANCE

The evolutionary redesign of snake COI coincided with adaptive bursts in other mitochondrial proteins and substantial changes in mitochondrial genome structure. It also generally coincided with or preceded major shifts in ecological niche and the evolution of extensive physiological adaptations related to lung reduction, large prey consumption, and venom evolution. The parallel timing of these major evolutionary events suggests that evolutionary redesign of metabolic and mitochondrial function may be related to, or underlie, the extreme changes in physiological and metabolic efficiency, flexibility, and innovation observed in snake evolution.

Thalamic nuclei in the opossum Monodelphis domestica

We investigated nuclear divisions of the thalamus in the gray short-tailed opossum (Monodelphis domestica) to gain detailed information for further developmental and comparative studies. Nissl and myelin staining, histochemistry for acetylcholinesterase and immunohistochemistry for calretinin and parvalbumin were performed on parallel series of sections. Many features of the Monodelphis opossum thalamus resemble those in Didelphis and small eutherians showing no particular sensory specializations, particularly in small murid rodents. However, several features of thalamic organization in Monodelphis were distinct from those in rodents. In the opossum the anterior and midline nuclear groups are more clearly separated from adjacent structures than in eutherians. The dorsal lateral geniculate nucleus (LGNd) starts more rostrally and occupies a large part of the lateral wall of the thalamus. As in other marsupials, two cytoarchitectonically different parts, alpha and beta are discernible in the LGNd of the opossum. Each of them may be subdivided into two additional bands in acetylcholinesterase staining, while in murid rodents the LGNd consists of a homogeneous mass of cells. Therefore, differentiation of the LGNd of the Monodelphis opossum is more advanced than in murid rodents. The medial geniculate body consists of three nuclei (medial, dorsal and ventral) that are cytoarchitectonically distinct and stain differentially for parvalbumin. The relatively large size of the MG and LGNd points to specialization of the visual and auditory systems in the Monodelphis opossum. In contrast to rodents, the lateral dorsal and lateral posterior nuclei in the opossum are poorly differentiated cytoarchitectonically.

Loss of genes implicated in gastric function during platypus evolution

BACKGROUND

The duck-billed platypus (Ornithorhynchus anatinus) belongs to the mammalian subclass Prototheria, which diverged from the Theria line early in mammalian evolution. The platypus genome sequence provides a unique opportunity to illuminate some aspects of the biology and evolution of these animals.

RESULTS

We show that several genes implicated in food digestion in the stomach have been deleted or inactivated in platypus. Comparison with other vertebrate genomes revealed that the main genes implicated in the formation and activity of gastric juice have been lost in platypus. These include the aspartyl proteases pepsinogen A and pepsinogens B/C, the hydrochloric acid secretion stimulatory hormone gastrin, and the alpha subunit of the gastric H+/K+-ATPase. Other genes implicated in gastric functions, such as the beta subunit of the H+/K+-ATPase and the aspartyl protease cathepsin E, have been inactivated because of the acquisition of loss-of-function mutations. All of these genes are highly conserved in vertebrates, reflecting a unique pattern of evolution in the platypus genome not previously seen in other mammalian genomes.

CONCLUSION

The observed loss of genes involved in gastric functions might be responsible for the anatomical and physiological differences in gastrointestinal tract between monotremes and other vertebrates, including small size, lack of glands, and high pH of the monotreme stomach. This study contributes to a better understanding of the mechanisms that underlie the evolution of the platypus genome, might extend the less-is-more evolutionary model to monotremes, and provides novel insights into the importance of gene loss events during mammalian evolution.

Transposable elements and the evolution of regulatory networks

The control and coordination of eukaryotic gene expression rely on transcriptional and post-transcriptional regulatory networks. Although progress has been made in mapping the components and deciphering the function of these networks, the mechanisms by which such intricate circuits originate and evolve remain poorly understood. Here I revisit and expand earlier models and propose that genomic repeats, and in particular transposable elements, have been a rich source of material for the assembly and tinkering of eukaryotic gene regulatory systems.

PiRNA-like RNAs in the marsupial Monodelphis domestica identify transcription clusters and likely marsupial transposon targets

PIWI-interacting RNAs (piRNAs) are a recently discovered class of small noncoding RNAs that have been detected in human, mouse, rat, zebrafish, and Drosophila genomes. We have utilized a size-directed small-RNA cloning procedure to clone and map more than 300 candidate piRNA-like small RNAs in the genome of the marsupial species Monodelphis domestica. Our results are consistent with those from other species in that the piRNA-like candidate sequences range in size from 28 to 31 nucleotides, show a pronounced preference for uridine at the 5' end, are transcribed from a few large clusters, appear to target transposons, and display virtually no sequence conservation.

A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach

MicroRNA (miRNA) and other types of small regulatory RNAs play a crucial role in the regulation of gene expression in eukaryotes. Several distinct classes of small regulatory RNAs have been discovered in recent years. To extend the repertoire of small regulatory RNAs characterized in chickens we used a deep sequencing approach developed by Solexa (now Illumina Inc.). We sequenced three small RNA libraries prepared from different developmental stages of the chicken embryo (days five, seven, and nine) to produce over 9.5 million short sequence reads. We developed a bioinformatics pipeline to distinguish authentic mature miRNA sequences from other classes of small RNAs and short RNA fragments represented in the sequencing data. Using this approach we detected almost all of the previously known chicken miRNAs and their respective miRNA* sequences. In addition we discovered 449 new chicken miRNAs including 88 miRNA candidates. Of these, 430 miRNAs appear to be specific to the avian lineage. Another six new miRNAs had evidence of evolutionary conservation in at least one vertebrate species outside of the bird lineage. The remaining 13 putative miRNAs appear to represent chicken orthologs of known vertebrate miRNAs. We discovered 39 additional putative miRNA candidates originating from miRNA generating intronic sequences known as mirtrons.

Genome analysis of the platypus reveals unique signatures of evolution

We present a draft genome sequence of the platypus, Ornithorhynchus anatinus. This monotreme exhibits a fascinating combination of reptilian and mammalian characters. For example, platypuses have a coat of fur adapted to an aquatic lifestyle; platypus females lactate, yet lay eggs; and males are equipped with venom similar to that of reptiles. Analysis of the first monotreme genome aligned these features with genetic innovations. We find that reptile and platypus venom proteins have been co-opted independently from the same gene families; milk protein genes are conserved despite platypuses laying eggs; and immune gene family expansions are directly related to platypus biology. Expansions of protein, non-protein-coding RNA and microRNA families, as well as repeat elements, are identified. Sequencing of this genome now provides a valuable resource for deep mammalian comparative analyses, as well as for monotreme biology and conservation.

Discordant evolution of the adjacent antiretroviral genes TRIM22 and TRIM5 in mammals

TRIM5α provides a cytoplasmic block to retroviral infection, and orthologs encoded by some primates are active against HIV. Here, we present an evolutionary comparison of the TRIM5 gene to its closest human paralogs: TRIM22, TRIM34, and TRIM6. We show that TRIM5 and TRIM22 have a dynamic history of gene expansion and loss during the evolution of mammals. The cow genome contains an expanded cluster of TRIM5 genes and no TRIM22 gene, while the dog genome encodes TRIM22 but has lost TRIM5. In contrast, TRIM6 and TRIM34 have been strictly preserved as single gene orthologs in human, dog, and cow. A more focused analysis of primates reveals that, while TRIM6 and TRIM34 have evolved under purifying selection, TRIM22 has evolved under positive selection as was previously observed for TRIM5. Based on TRIM22 sequences obtained from 27 primate genomes, we find that the positive selection of TRIM22 has occurred episodically for approximately 23 million years, perhaps reflecting the changing pathogenic landscape. However, we find that the evolutionary episodes of positive selection that have acted on TRIM5 and TRIM22 are mutually exclusive, with generally only one of these genes being positively selected in any given primate lineage. We interpret this to mean that the positive selection of one gene has constrained the adaptive flexibility of its neighbor, probably due to genetic linkage. Finally, we find a striking congruence in the positions of amino acid residues found to be under positive selection in both TRIM5α and TRIM22, which in both proteins fall predominantly in the β2-β3 surface loop of the B30.2 domain. Astonishingly, this same loop is under positive selection in the multiple cow TRIM5 genes as well, indicating that this small structural loop may be a viral recognition motif spanning a hundred million years of mammalian evolution.

The intrinsic immunity protein TRIM5α provides a post-entry defense against retroviral infection, which depends on its specific ability to recognize retroviral capsids. TRIM5α has been locked in genetic conflict with retroviruses throughout most of primate evolution, characterized by a higher than expected rate of amino acid change, referred to as positive selection. Here, we find that one of TRIM5's closest human paralogs, TRIM22, has also undergone positive selection in primates. However, we find that its close linkage to TRIM5 has resulted in an anti-correlated pattern of positive selection, with primate lineages generally showing positive selection in either TRIM5 or TRIM22, but not both. Amino acid positions in TRIM22 found to be under positive selection are in remarkable proximity to the “antiviral specificity patch” previously described for TRIM5α. TRIM5 and TRIM22 evolution appears to be equally discordant in other mammals; the cow genome contains an expanded cluster of TRIM5 genes and no TRIM22 gene, while the dog genome encodes TRIM22 but has lost TRIM5. Our analyses highlight TRIM22 as bearing all the evolutionary hallmarks of a candidate intrinsic immunity gene.

Addressing chromosome evolution in the whole-genome sequence era

The evolution of karyotypes has been the subject of intensive study since the middle of the 20th century. This was motivated by the observation that the karyotypes of related species showed remarkable conservation. The recent emergence of whole-genome sequencing projects gives the opportunity to complement the cytogenetic approaches by addressing the conservation of karyotypes using chromosome sequence comparison. In this short review we present a description of recent advances in computational biology methods dedicated to the study of chromosome evolution and more specifically ancestral karyotype reconstruction in an attempt to provide an integrated overview of both cytogenetic and computational approaches.

Comparative genomics search for losses of long-established genes on the human lineage

Taking advantage of the complete genome sequences of several mammals, we developed a novel method to detect losses of well-established genes in the human genome through syntenic mapping of gene structures between the human, mouse, and dog genomes. Unlike most previous genomic methods for pseudogene identification, this analysis is able to differentiate losses of well-established genes from pseudogenes formed shortly after segmental duplication or generated via retrotransposition. Therefore, it enables us to find genes that were inactivated long after their birth, which were likely to have evolved nonredundant biological functions before being inactivated. The method was used to look for gene losses along the human lineage during the approximately 75 million years (My) since the common ancestor of primates and rodents (the euarchontoglire crown group). We identified 26 losses of well-established genes in the human genome that were all lost at least 50 My after their birth. Many of them were previously characterized pseudogenes in the human genome, such as GULO and UOX. Our methodology is highly effective at identifying losses of single-copy genes of ancient origin, allowing us to find a few well-known pseudogenes in the human genome missed by previous high-throughput genome-wide studies. In addition to confirming previously known gene losses, we identified 16 previously uncharacterized human pseudogenes that are definitive losses of long-established genes. Among them is ACYL3, an ancient enzyme present in archaea, bacteria, and eukaryotes, but lost approximately 6 to 8 Mya in the ancestor of humans and chimps. Although losses of well-established genes do not equate to adaptive gene losses, they are a useful proxy to use when searching for such genetic changes. This is especially true for adaptive losses that occurred more than 250,000 years ago, since any genetic evidence of the selective sweep indicative of such an event has been erased.

Non-coding RNAs in imprinted gene clusters

Imprinted ncRNA (non-coding RNA) genes represent a family of untranslated transcripts that are mono-allelically expressed in a parent-of-origin manner (their expression is restricted to either the maternal or the paternal allele). Although the expression of a few long imprinted ncRNAs act as cis-acting silencers in the epigenetic regulation of chromatin at imprinted gene clusters, many of them fall into the growing class of small regulatory RNAs, namely C/D small nucleolar RNAs, microRNAs and also likely piRNAs (Piwi-interacting RNAs), which are known to act as antisense trans-acting regulators of gene expression (for example, site-specific RNA modifications and RNA-mediated gene silencing). Although their biological functions remain elusive, recent studies have pointed to their functional importance in development, in brain plasticity and also perhaps in some pathological situations, such as cancers or Prader-Willi syndrome. In this review, we summarize our current understanding of the molecular and biological roles of these ncRNAs, both in terms of their contribution to genomic imprinting control, as well as in terms of cellular RNA targets they might interact with.

Genome assembly forensics: finding the elusive mis-assembly

A collection of software tools is combined for the first time in an automated pipeline for detecting large-scale genome assembly errors and for validating genome assemblies.

We present the first collection of tools aimed at automated genome assembly validation. This work formalizes several mechanisms for detecting mis-assemblies, and describes their implementation in our automated validation pipeline, called amosvalidate. We demonstrate the application of our pipeline in both bacterial and eukaryotic genome assemblies, and highlight several assembly errors in both draft and finished genomes. The software described is compatible with common assembly formats and is released, open-source, at .

Possible involvement of SINEs in mammalian-specific brain formation

Retroposons, such as short interspersed elements (SINEs) and long interspersed elements (LINEs), are the major constituents of higher vertebrate genomes. Although there are many examples of retroposons' acquiring function, none has been implicated in the morphological innovations specific to a certain taxonomic group. We previously characterized a SINE family, AmnSINE1, members of which constitute a part of conserved noncoding elements (CNEs) in mammalian genomes. We proposed that this family acquired genomic functionality or was exapted after retropositioning in a mammalian ancestor. Here we identified 53 new AmnSINE1 loci and refined 124 total loci, two of which were further analyzed. Using a mouse enhancer assay, we demonstrate that one SINE locus, AS071, 178 kbp from the gene FGF8 (fibroblast growth factor 8), is an enhancer that recapitulates FGF8 expression in two regions of the developing forebrain, namely the diencephalon and the hypothalamus. Our gain-of-function analysis revealed that FGF8 expression in the diencephalon controls patterning of thalamic nuclei, which act as a relay center of the neocortex, suggesting a role for FGF8 in mammalian-specific forebrain patterning. Furthermore, we demonstrated that the locus, AS021, 392 kbp from the gene SATB2, controls gene expression in the lateral telencephalon, which is thought to be a signaling center during development. These results suggest important roles for SINEs in the development of the mammalian neuronal network, a part of which was initiated with the exaptation of AmnSINE1 in a common mammalian ancestor.

Eutherians intrinsically run a higher risk of replication deficiency

Females inherit two X chromosomes from each parent, but males inherit only one from the mother. The POLA gene encodes the catalytic subunit of DNA polymerase alpha (pol-alpha), synthesizing a short DNA primer after RNA priming for DNA chain elongation in eukaryotic DNA replication. The POLA gene is established in the X-added region (XRA) of an X chromosome, translocated from a second autosome about 105 million years ago (mya). Dosage compensation for the mammalian X chromosome is accomplished by the silencing of one X chromosome to equalize gene expression between the male and female. We eutherians must replicate our DNA using pol-alpha expressed from a 'single copy' of genes, whose expression is probably reduced when chromosome translocation occurs. Spontaneous mutation occurring on the 'single' POLA gene, which decreases its replication fidelity, probably accounts for the increase of nucleotide substitution rates more effectively than in other organisms that possess the POLA gene on an autosome pair. Moreover, translocation of the POLA gene from an autosome to the X chromosome might reduce expression of pol-alpha, resulting in promotion of chromosome rearrangement. Therefore, the X-linked POLA gene possibly contributes to greater diversification of eutherian mammals.

Building up the inactive X chromosome

The compensation of the different level of transcripts of X-linked genes in male and female mammals is achieved through X chromosome inactivation, a complex process that differentially regulates the sex chromosomes of female cells. This mechanism has been dissected at evolutionary, genetic and molecular levels: here, we discuss some of the latest examples that illustrate better these intricate connections, focusing particularly on the emerging role of spatial and three-dimensional chromatin arrangements in the building of this special chromosome, the inactive X chromosome.

The Opossum genome reveals further evidence for regulatory evolution in mammalian diversification

The first marsupial genome to be sequenced emphasizes the importance of noncoding sequences in mammalian evolution.

The sequencing of the euchromatic genome of a marsupial, the opossum Monodelphis domestica, identifies shared and unique features of marsupial and placental genomes and reveals a prominent role for the evolution of non-protein-coding elements.

Comparative genomic analysis and evolution of the T cell receptor loci in the opossum Monodelphis domestica

BACKGROUND

All jawed-vertebrates have four T cell receptor (TCR) chains: alpha (TRA), beta (TRB), gamma (TRG) and delta (TRD). Marsupials appear unique by having an additional TCR: mu (TRM). The evolutionary origin of TRM and its relationship to other TCR remain obscure, and is confounded by previous results that support TRM being a hybrid between a TCR and immunoglobulin locus. The availability of the first marsupial genome sequence allows investigation of these evolutionary relationships.

RESULTS

The organization of the conventional TCR loci, encoding the TRA, TRB, TRG and TRD chains, in the opossum Monodelphis domestica are highly conserved with and of similar complexity to that of eutherians (placental mammals). There is a high degree of conserved synteny in the genomic regions encoding the conventional TCR across mammals and birds. In contrast the chromosomal region containing TRM is not well conserved across mammals. None of the conventional TCR loci contain variable region gene segments with homology to those found in TRM; rather TRM variable genes are most similar to that of immunoglobulin heavy chain genes.

CONCLUSION

Complete genomic analyses of the opossum TCR loci continue to support an origin of TRM as a hybrid between a TCR and immunoglobulin locus. None of the conventional TCR loci contain evidence that such a recombination event occurred, rather they demonstrate a high degree of stability across distantly related mammals. TRM, therefore, appears to be derived from receptor genes no longer extant in placental mammals. These analyses provide the first genomic scale structural detail of marsupial TCR genes, a lineage of mammals used as models of early development and human disease.

Molecular and cytogenetic analysis of the spreading of X inactivation in a girl with microcephaly, mild dysmorphic features and t(X;5)(q22.1;q31.1)

X chromosome inactivation involves initiation, propagation, and maintenance of gene inactivation. Studies of replication pattern and timing in X;autosome translocations have suggested that X inactivation may spread to autosomal DNA. To examine this phenomenon at the molecular level, we have tested the transcriptional activity of a number of chromosome 5 loci in a female subject with microcephaly, mild dysmorphic features and 46,X,der(X)t(X;5)(q22.1;q31.1) karyotype. RT-PCR analysis of 20 transcribed sequences spanning 5q31.1-qter revealed that nine of them were not expressed in somatic cell hybrid clones carrying the translocated chromosome. However, eight genes were expressed and therefore escaped inactivation. This direct expression test demonstrates that spreading of inactivation from the X chromosome to the adjoining autosomal DNA was incomplete and 'patchy'. Inactivation was associated in most instances to methylation of the CpG sequences in genes containing CpG islands, but was also present in CpG islandless genes. These results agree with those obtained for other X;autosome translocations and demonstrate that autosomes are partially resistant to Xist-mediated spreading and/or maintenance of inactivation. Repeat distribution analysis does not suggest an association between L1 and LINE repeat density on chromosome 5 and gene inactivation. The expression data may also explain why the proband manifests an attenuated clinical phenotype compared to subjects with partial chromosome 5 trisomy.

Opossum carboxylesterases: sequences, phylogeny and evidence for CES gene duplication events predating the marsupial-eutherian common ancestor

Background

Carboxylesterases (CES) perform diverse metabolic roles in mammalian organisms in the detoxification of a broad range of drugs and xenobiotics and may also serve in specific roles in lipid, cholesterol, pheromone and lung surfactant metabolism. Five CES families have been reported in mammals with human CES1 and CES2 the most extensively studied. Here we describe the genetics, expression and phylogeny of CES isozymes in the opossum and report on the sequences and locations of CES1, CES2 and CES6 'like' genes within two gene clusters on chromosome one. We also discuss the likely sequence of gene duplication events generating multiple CES genes during vertebrate evolution.

Results

We report a cDNA sequence for an opossum CES and present evidence for CES1 and CES2 like genes expressed in opossum liver and intestine and for distinct gene locations of five opossum CES genes,CES1, CES2.1, CES2.2, CES2.3 and CES6, on chromosome 1. Phylogenetic and sequence alignment studies compared the predicted amino acid sequences for opossum CES with those for human, mouse, chicken, frog, salmon and Drosophila CES gene products. Phylogenetic analyses produced congruent phylogenetic trees depicting a rapid early diversification into at least five distinct CES gene family clusters: CES2, CES1, CES7, CES3, and CES6. Molecular divergence estimates based on a Bayesian relaxed clock approach revealed an origin for the five mammalian CES gene families between 328–378 MYA.

Conclusion

The deduced amino acid sequence for an opossum cDNA was consistent with its identity as a mammalian CES2 gene product (designated CES2.1). Distinct gene locations for opossum CES1 (1: 446,222,550–446,274,850), three CES2 genes (1: 677,773,395–677,927,030) and a CES6 gene (1: 677,585,520–677,730,419) were observed on chromosome 1. Opossum CES1 and multiple CES2 genes were expressed in liver and intestine. Amino acid sequences for opossum CES1 and three CES2 gene products revealed conserved residues previously reported for human CES1 involved in catalysis, ligand binding, tertiary structure and organelle localization. Phylogenetic studies indicated the gene duplication events which generated ancestral mammalian CES genes predated the common ancestor for marsupial and eutherian mammals, and appear to coincide with the early diversification of tetrapods.

Ancora: a web resource for exploring highly conserved noncoding elements and their association with developmental regulatory genes

Ancora is a web resource that provides data and tools for exploring genomic organization of highly conserved noncoding elements for multiple genomes.

Metazoan genomes contain arrays of highly conserved noncoding elements (HCNEs) that span developmental regulatory genes and define regulatory domains. We describe Ancora , a web resource that provides data and tools for exploring genomic organization of HCNEs for multiple genomes. Ancora includes a genome browser that shows HCNE locations and features novel HCNE density plots as a powerful tool to discover developmental regulatory genes and distinguish their regulatory elements and domains.

The genome of a Gondwanan mammal

Australia is thought of as the home of marsupials, but South America has 60 or so species of these interesting mammals. The genome of one of these, the South American grey short-tailed opossum, Monodelphis domestica, has just been sequenced and published in June.1 The high quality 6x coverage is the first marsupial genome completed, pipping the 2x coverage of the Australian tammar wallaby at the post by half a year. The opossum genome has an unusual structure with fewer chromosomes than the human genome (9 pairs versus 23 pairs) but a longer total length (3.4 billion versus 3 billion bases). The opossum autosomes, like those of all marsupials, are extremely large but, in contrast, the X chromosome is only 76 Mb long. The opossum genome has turned up several surprises and provided critical new information on the evolution of mammalian genomes.

Analysis of sequence variability in the macronuclear DNA of Paramecium tetraurelia: a somatic view of the germline

Ciliates are the only unicellular eukaryotes known to separate germinal and somatic functions. Diploid but silent micronuclei transmit the genetic information to the next sexual generation. Polyploid macronuclei express the genetic information from a streamlined version of the genome but are replaced at each sexual generation. The macronuclear genome of Paramecium tetraurelia was recently sequenced by a shotgun approach, providing access to the gene repertoire. The 72-Mb assembly represents a consensus sequence for the somatic DNA, which is produced after sexual events by reproducible rearrangements of the zygotic genome involving elimination of repeated sequences, precise excision of unique-copy internal eliminated sequences (IES), and amplification of the cellular genes to high copy number. We report use of the shotgun sequencing data (>10(6) reads representing 13 x coverage of a completely homozygous clone) to evaluate variability in the somatic DNA produced by these developmental genome rearrangements. Although DNA amplification appears uniform, both of the DNA elimination processes produce sequence heterogeneity. The variability that arises from IES excision allowed identification of hundreds of putative new IESs, compared to 42 that were previously known, and revealed cases of erroneous excision of segments of coding sequences. We demonstrate that IESs in coding regions are under selective pressure to introduce premature termination of translation in case of excision failure.

Mammalian non-LTR retrotransposons: for better or worse, in sickness and in health

Transposable elements (TEs) have shared an exceptionally long coexistence with their host organisms and have come to occupy a significant fraction of eukaryotic genomes. The bulk of the expansion occurring within mammalian genomes has arisen from the activity of type I retrotransposons, which amplify in a "copy-and-paste" fashion through an RNA intermediate. For better or worse, the sequences of these retrotransposons are now wedded to the genomes of their mammalian hosts. Although there are several reported instances of the positive contribution of mobile elements to their host genomes, these discoveries have occurred alongside growing evidence of the role of TEs in human disease and genetic instability. Here we examine, with a particular emphasis on human retrotransposon activity, several newly discovered aspects of mammalian retrotransposon biology. We consider their potential impact on host biology as well as their ultimate implications for the nature of the TE-host relationship.

Rats in the genomic era

The rat genome project and the resources that it has generated are transforming the translation of rat biology to human medicine. The rat genome was sequenced to a high quality “draft,” the structure and location of the genes were predicted, and a global assessment was published (Gibbs RA et al., Nature 428: 493–521, 2004). Since that time, researchers have made use of the genome sequence and annotations and related resources. We take this opportunity to review the currently available rat genome resources and to discuss the progress and future plans for the rat genome.

A machine-learning approach to combined evidence validation of genome assemblies

MOTIVATION

While it is common to refer to 'the genome sequence' as if it were a single, complete and contiguous DNA string, it is in fact an assembly of millions of small, partially overlapping DNA fragments. Sophisticated computer algorithms (assemblers and scaffolders) merge these DNA fragments into contigs, and place these contigs into sequence scaffolds using the paired-end sequences derived from large-insert DNA libraries. Each step in this automated process is susceptible to producing errors; hence, the resulting draft assembly represents (in practice) only a likely assembly that requires further validation. Knowing which parts of the draft assembly are likely free of errors is critical if researchers are to draw reliable conclusions from the assembled sequence data.

RESULTS

We develop a machine-learning method to detect assembly errors in sequence assemblies. Several in silico measures for assembly validation have been proposed by various researchers. Using three benchmarking Drosophila draft genomes, we evaluate these techniques along with some new measures that we propose, including the good-minus-bad coverage (GMB), the good-to-bad-ratio (RGB), the average Z-score (AZ) and the average absolute Z-score (ASZ). Our results show that the GMB measure performs better than the others in both its sensitivity and its specificity for assembly error detection. Nevertheless, no single method performs sufficiently well to reliably detect genomic regions requiring attention for further experimental verification. To utilize the advantages of all these measures, we develop a novel machine learning approach that combines these individual measures to achieve a higher prediction accuracy (i.e. greater than 90%). Our combined evidence approach avoids the difficult and often ad hoc selection of many parameters the individual measures require, and significantly improves the overall precisions on the benchmarking data sets.

Non-coding RNAs, epigenetics and complexity

Several aspects of epigenetics are strongly linked to non-coding RNAs, especially small RNAs that can direct the cytosine methylation and histone modifications that are implicated in gene expression regulation in complex organisms. A fundamental characteristic of epigenetics is that the same genome can show alternative phenotypes, which are based in different epigenetic states. Some of the most studied complex epigenetic phenomena including transposon activity and silencing recently exemplified by piRNAs (piwi-interacting RNAs), position effect variegation, X-chromosome inactivation, parental imprinting, and paramutation have direct or indirect participation of an RNA component. Conceivably, most of the non-coding RNAs with no described function yet, are players in epigenetic mechanisms that are still not completely understood. In that regard, RNAs were recently implicated in new mechanisms of genetic information transfer in yeast, plants and mice. In this review article, the hypothesis that non-coding RNAs might be the main component of complex organisms acquired during evolution will be explored. The question of how evolutionary theories have been challenged by these molecules in association with epigenetic mechanisms will also be discussed here.

Sources for comparative studies of placentation. II. Genomic resources

The genomes of over a dozen placental mammal species are now publicly available. These genome sequences have the potential to provide insight into the development and evolution of the placenta. In particular, the variable anatomy of the placenta has likely been affected by natural selection on the genomes of living and extinct mammals. In this note the current availability of mammal genome sequences is reviewed, and strengths and limitations of these data are discussed. Additionally, museums, zoos, and commercial entities are available to provide genomic resources to the placental research community. Recommendations for tissue storage conditions of placentas in genomic research are given.

Maintenance of transposon-free regions throughout vertebrate evolution

Background

We recently reported the existence of large numbers of regions up to 80 kb long that lack transposon insertions in the human, mouse and opossum genomes. These regions are significantly associated with loci involved in developmental and transcriptional regulation.

Results

Here we report that transposon-free regions (TFRs) are prominent genomic features of amphibian and fish lineages, and that many have been maintained throughout vertebrate evolution, although most transposon-derived sequences have entered these lineages after their divergence. The zebrafish genome contains 470 TFRs over 10 kb and a further 3,951 TFRs over 5 kb, which is comparable to the number identified in mammals. Two thirds of zebrafish TFRs over 10 kb are orthologous to TFRs in at least one mammal, and many have orthologous TFRs in all three mammalian genomes as well as in the genome of Xenopus tropicalis. This indicates that the mechanism responsible for the maintenance of TFRs has been active at these loci for over 450 million years. However, the majority of TFR bases cannot be aligned between distantly related species, demonstrating that TFRs are not the by-product of strong primary sequence conservation. Syntenically conserved TFRs are also more enriched for regulatory genes compared to lineage-specific TFRs.

Conclusion

We suggest that TFRs contain extended regulatory sequences that contribute to the precise expression of genes central to early vertebrate development, and can be used as predictors of important regulatory regions.

28-way vertebrate alignment and conservation track in the UCSC Genome Browser

This article describes a set of alignments of 28 vertebrate genome sequences that is provided by the UCSC Genome Browser. The alignments can be viewed on the Human Genome Browser (March 2006 assembly) at http://genome.ucsc.edu, downloaded in bulk by anonymous FTP from http://hgdownload.cse.ucsc.edu/goldenPath/hg18/multiz28way, or analyzed with the Galaxy server at http://g2.bx.psu.edu. This article illustrates the power of this resource for exploring vertebrate and mammalian evolution, using three examples. First, we present several vignettes involving insertions and deletions within protein-coding regions, including a look at some human-specific indels. Then we study the extent to which start codons and stop codons in the human sequence are conserved in other species, showing that start codons are in general more poorly conserved than stop codons. Finally, an investigation of the phylogenetic depth of conservation for several classes of functional elements in the human genome reveals striking differences in the rates and modes of decay in alignability. Each functional class has a distinctive period of stringent constraint, followed by decays that allow (for the case of regulatory regions) or reject (for coding regions and ultraconserved elements) insertions and deletions.

MANTIS: a phylogenetic framework for multi-species genome comparisons

MOTIVATION

Practitioners of comparative genomics face huge analytical challenges as whole genome sequences and functional/expression data accumulate. Furthermore, the field would greatly benefit from a better integration of this wealth of data with evolutionary concepts.

RESULTS

Here, we present MANTIS, a relational database for the analysis of (i) gains and losses of genes on specific branches of the metazoan phylogeny, (ii) reconstructed genome content of ancestral species and (iii) over- or under-representation of functions/processes and tissue specificity of gained, duplicated and lost genes. MANTIS estimates the most likely positions of gene losses on the true phylogeny using a maximum-likelihood function. A user-friendly interface and an extensive query system allow to investigate questions pertaining to gene identity, phylogenetic mapping and function/expression parameters.

AVAILABILITY

MANTIS is freely available at http://www.mantisdb.org and constitutes the missing link between multi-species genome comparisons and functional analyses.

Species-specific endogenous retroviruses shape the transcriptional network of the human tumor suppressor protein p53

The evolutionary forces that establish and hone target gene networks of transcription factors are largely unknown. Transposition of retroelements may play a role, but its global importance, beyond a few well described examples for isolated genes, is not clear. We report that LTR class I endogenous retrovirus (ERV) retroelements impact considerably the transcriptional network of human tumor suppressor protein p53. A total of 1,509 of approximately 319,000 human ERV LTR regions have a near-perfect p53 DNA binding site. The LTR10 and MER61 families are particularly enriched for copies with a p53 site. These ERV families are primate-specific and transposed actively near the time when the New World and Old World monkey lineages split. Other mammalian species lack these p53 response elements. Analysis of published genomewide ChIP data for p53 indicates that more than one-third of identified p53 binding sites are accounted for by ERV copies with a p53 site. ChIP and expression studies for individual genes indicate that human ERV p53 sites are likely part of the p53 transcriptional program and direct regulation of p53 target genes. These results demonstrate how retroelements can significantly shape the regulatory network of a transcription factor in a species-specific manner.