Skip to main content
Advertisement
  • Loading metrics

Retrotransposon Silencing by DNA Methylation Can Drive Mammalian Genomic Imprinting

Abstract

Among mammals, only eutherians and marsupials are viviparous and have genomic imprinting that leads to parent-of-origin-specific differential gene expression. We used comparative analysis to investigate the origin of genomic imprinting in mammals. PEG10 (paternally expressed 10) is a retrotransposon-derived imprinted gene that has an essential role for the formation of the placenta of the mouse. Here, we show that an orthologue of PEG10 exists in another therian mammal, the marsupial tammar wallaby (Macropus eugenii), but not in a prototherian mammal, the egg-laying platypus (Ornithorhynchus anatinus), suggesting its close relationship to the origin of placentation in therian mammals. We have discovered a hitherto missing link of the imprinting mechanism between eutherians and marsupials because tammar PEG10 is the first example of a differentially methylated region (DMR) associated with genomic imprinting in marsupials. Surprisingly, the marsupial DMR was strictly limited to the 5′ region of PEG10, unlike the eutherian DMR, which covers the promoter regions of both PEG10 and the adjacent imprinted gene SGCE. These results not only demonstrate a common origin of the DMR-associated imprinting mechanism in therian mammals but provide the first demonstration that DMR-associated genomic imprinting in eutherians can originate from the repression of exogenous DNA sequences and/or retrotransposons by DNA methylation.

Author Summary

Genomic imprinting is a gene regulatory mechanism controlling parent-of-origin-dependent expression of genes. In eutherians, imprinting is essential for fetal and placental development and defects in this mechanism are the cause of several genetic disorders. In eutherian mammals, genomic imprinting is controlled by differential methylation of the DNA. However, no such methylation-dependent mechanism had been previously identified in association with marsupial imprinting. By comparing the genome of all three extant classes of mammals (eutherians, marsupials, and monotremes), we have investigated the evolution of PEG10 (paternally expressed 10), a retrotransposon-derived imprinted gene that is essential for the formation of the placenta in the mouse. PEG10 was present in a marsupial species, the tammar wallaby, but absent from an egg-laying monotreme species, the platypus. Therefore, PEG10 was inserted into the genome at the time when the placenta and viviparity were evolving in therian mammals. This study has shown that PEG10 is not only imprinted in a marsupial, but that its imprint is regulated by differential methylation, suggesting a common origin for methylation in the therian ancestor. These results provide direct evidence that retrotransposon insertion can drive the evolution of genomic imprinting in mammals.

Introduction

Genomic imprinting, or parent-of-origin-specific gene silencing, has been observed in both eutherian and marsupial, but not monotreme mammals. In eutherians, more than 80 imprinted genes have been found and differential DNA methylation plays a crucial role in regulating their imprinted expression patterns [15]. In marsupials, three genes—IGF2, IGF2R, and PEG1/MEST—are imprinted, as they are in eutherians, but no differentially methylated regions (DMRs) have been found in marsupials [68]. Consequently, the regulatory mechanisms of genomic imprinting were thought to have evolved differently between marsupials and eutherians [6,8,9].

It has been hypothesized that genomic imprinting arose as a by-product of a DNA methylation mechanism that silences foreign DNAs [10], such as retrotransposons [11,12]. Similarly, transgenes can also become methylated, depending on parent of origin, further supporting a link between genomic imprinting and silencing of foreign DNAs [1315]. PEG10 is an imprinted gene sharing homology with the sushi-ichi retrotransposon, and in humans and mice it has a clear DMR in its promoter region. Interestingly, PEG10 is conserved in eutherian mammals but not in nonmammalian vertebrates, such as birds and fish [11,1618]. Therefore, investigating the origin of the retrotransposon-derived PEG10 locus would clarify the relationship between retrotransposon (or exogenous DNA sequence) insertion and genomic imprinting.

PEG10 is an essential placental gene in eutherians, since knock-out mice have severe placental defects with loss of spongiotrophoblast and labyrinth layers leading to early embryonic lethality [18]. The origin of PEG10 is therefore of interest in view of its possible contribution to the evolution of mammalian placentation. Sequence identified as PEG10 has recently been reported in the South American marsupial the grey short-tailed opossum (Monodelphis domestica) [16] but its precise location, genetic structure, and imprint status remains unknown. Here, we examined the PEG10 locus in two Australian mammals by isolating bacterial artificial chromosome (BAC) clones from a marsupial, the tammar wallaby, and from a monotreme, the platypus.

Results/Discussion

We used the tammar wallaby and the platypus as the representative species of marsupials and monotremes, respectively, to compare with representative eutherians, the human and the mouse. Several BAC clones were isolated from tammar and platypus containing SGCE (sarcoglycan epsilon), the neighboring gene of PEG10 in eutherians. DNA sequencing of one tammar BAC clone demonstrated the existence of PEG10 and its conserved location adjacent to SGCE with transcription occurring in a head-to-head manner (Figure 1A). The genetic structure of tammar PEG10 was also the same as that of eutherians, with two open reading frames (ORFs) related to the sushi-ichi retrotransposon GAG and POL proteins, respectively [11]. The CX2CX4HX4C RNA-binding motif of GAG and the DSG sequences of the proteinase activation site of POL were also conserved (Figure S1). Tammar PEG10 was localised close to the telomere of Chromosome 3q, consistent with its autosomal location on proximal mouse Chromosome 6 (Figure 1B). However, in the platypus, there were no PEG10 homologous sequences between SGCE and PPP1R9A (also called NEURABIN1) that flank PEG10 in other mammals (Figure 1A).

thumbnail
Figure 1. Comparison of the Genomic sStructures between SGCE and PPP1R9A, and Chromosomal Location for Tammar PEG10

(A) The green bars on the horizontal lines indicate the exons of each gene. Orthologous exons are connected between species by the broken lines. The grey bars on the upper side of the horizontal lines represent LINEs and SINEs, and the blue bars on the underside represent LTR elements. The arrowheads indicate the direction of transcription. The published opossum genome sequence used was PEG10 downstream, and was combined with our tammar genome sequence at the homologous point in the 3′ UTR of PEG10.

(B) Fluorescence in situ hybridization of PEG10 was localized to Chromosome 3q. There is a nucleus of another cell in the lower left.

https://doi.org/10.1371/journal.pgen.0030055.g001

Using our tammar sequence and the published opossum genome sequence, we compared the entire region between SGCE and PPP1R9A with the equivalent region in several vertebrates from fish (fugu) to mammals. The size of this region in the platypus was similar to that of the chicken and was smaller than in other mammals. There were numerous long interspersed nuclear elements (LINEs) and short interspersed nuclear elements (SINEs) (grey bars in Figure 1A) in all mammalian groups. A previous report suggests that there are less long terminal repeat (LTR)-type retrotransposon-derived sequences in the opossum and wallaby genomes [19], but a large number of these sequences were found in this region as well as in mouse and human (blue bars in Figure 1A). Consistent with the previous report, these were absent in the platypus [19], as was PEG10. Most of the LTR-type retrotransposon-derived sequences observed in these regions are species specific. This suggests that most of these insertions occurred after species diversification. The presence of PEG10 sharing homology with one of the LTR-type retrotransposons, sushi-ichi [11], in both marsupials and eutherians suggests that the original PEG10 sequence insertion in the common therian ancestor was an early event in the therian-specific expansion of LTR-type retrotransposons. These results indicate that the original PEG10 was inserted into the genome of the therian ancestor and evolved to its present function as an essential placental gene after divergence from the monotremes. The acquisition of a new function for an existing character during evolution, a process termed “exaptation” by Gould and colleagues, would be enhanced by the provision of novel genetic materials, such as retrotransposons [20,21]. Thus, the requirement for PEG10 in placental function is a clear example of “exaptation.” The fossil record shows that there was extensive radiation of therian mammals after their split from the Prototheria. New LTR-type retrotranposon-derived sequences might therefore have contributed novel genetic resources to this radiation.

In mice, there is a large imprinting cluster near Peg10 which includes the paternally expressed Sgce gene [22,23] Ppp1r9a, which is maternally biased in extraembryonic tissues [23] and Asb4, which is completely maternally expressed in both embryos and the extraembryonic tissues [23,24]. Tammar PEG10 showed almost complete monoallelic expression in all individuals. Paternal expression was confirmed in two embryos and one yolk sac placenta sample (Figure 2). Unexpectedly, tammar SGCE showed predominantly biallelic expression with only a small paternal bias, despite a short 200-bp distance between the transcription start sites of PEG10 and SGCE (Figure 2). PPP1R9A and ASB4 showed biallelic expression without parental bias (Figure 2). These results clearly demonstrate that imprinting in this region is restricted to the PEG10 gene in the tammar. As described above, the eutherian PEG10 imprinted region includes several neighboring genes, suggesting that the imprinted region expanded in the eutherians while in marsupials imprinting was restricted to PEG10.

thumbnail
Figure 2. Allelic Expression Analysis of Tammar PEG10 and Neighboring Genes

The individual animal numbers are shown above the graphs and the pictures. The pairs of white and grey bars in the graphs represent the two parental alleles in each individual. The vertical axes represent the expression ratios of each allele percentage of total expression, and the raw data are shown under each bar. The characters under the raw data indicate the residues at each polymorphic site, and asterisks are added to the paternal alleles in the informative cases. The gel picture for No. 5 shows the result of RT-PCR for the individual with a PEG10 length polymorphism. Parental origin of both alleles are represented by “P” for paternal and “M” for maternal, respectively. The picture for No. 9 shows the results of direct sequencing for each PCR product. The arrows indicate the polymorphic sites.

https://doi.org/10.1371/journal.pgen.0030055.g002

A CpG island is present in the promoter regions of SGCE and PEG10 in tammar as well as in mouse. To determine why SGCE did not show imprinted expression we examined its methylation status. Surprisingly, we found a DMR with a clear boundary of DNA methylation between the PEG10 side and the SGCE side of the CpG island in both embryos and yolk sac placentas (Figure 3). Furthermore, selective DNA methylation of the maternal allele was confirmed using a DNA polymorphism in this region as predicted by the paternal expression of PEG10. The DNA methylation started about 60 bp downstream from the transcription start site of PEG10, suggesting that maternal transcription is inhibited by methylation of downstream regulatory elements and not by the typical mechanism of promoter methylation (Figure 3). Both LTR [25] and non-LTR retrotransposons [26] are known to have internal transcriptional regulatory elements for their transcription. As PEG10 is a retrotransposon-derived gene, these elements may exist within the DMR not, as is usual, upstream of the transcription start site. A part of the CpG island was possibly derived from LTRs in the original PEG10 sequence and the methylation that originated from a host–defense mechanism may be restricted to the ancient repetitive-element homologous region. Alternatively, a boundary function for DNA methylation spreading and/or transcription regulation may be included in the marsupial CpG island. The retrotransposon LTR sequences are CpG-rich and have such a boundary function in Saccharomyces cerevisiae and Drosophila melanogaster [27,28]. In Drosophila, the boundary function has been attributed to the binding of the SU(HW) protein. A consensus SU(HW) binding site was not found around the methylation boundary in the tammar PEG10. However, the CTCF protein, well known to have a similar insulator/boundary function in mammals, may bind to the possible boundary elements containing CT-rich sequences in this region.

thumbnail
Figure 3. DNA Methylation Status of the CpG Island over the Promoter Regions of SGCE and PEG10

The regions corresponding to SGCE and PEG10 are represented by the boxes. SGCE (yellow box) has biallelic expression whilst PEG10 (blue box) shows paternal expression. Black circles indicate methylated CpGs and white circles indicate unmethylated CpGs. Parental origin of each clone is distinguished by the polymorphism and the mother's genotype, and is represented by “P” for paternal and “M” for maternal, respectively.

https://doi.org/10.1371/journal.pgen.0030055.g003

Even with the presence of a DMR, it is possible that the maternal copy of PEG10 in the tammar is silenced by another mechanism and is only secondarily methylated. We therefore examined whether the imprinted expression of tammar PEG10 was regulated by DNA methylation. A reduced level of DNA methylation was observed in three sites of the CpG island in cells cultured with 5-aza-2′-deoxycytidine, a DNA methylation inhibitor (Figure 4A). Repetitive experiments performed for the most 3′ site using three independent cell lines established from fetal lung and endometrium also showed statistically significant reductions in DNA methylation levels (Figure 4B), and increased PEG10 expression from normally repressed alleles was observed in each case (Figure 4C, black and grey bars), although the expression levels were still much lower than active alleles (Figure 4C, white bars). These results demonstrate the association between imprinted expression of PEG10 and DNA methylation in a marsupial, although it still remains unknown if the differential methylation originates in the germline as does a typical primary DMR in eutherians.

thumbnail
Figure 4. DNA Methylation and Allelic Expression of PEG10 in 5-Aza-2′-deoxycytidine Treated Tammar Cells

(A) The bars on the horizontal line represent each CpG site. Three CpG sites were examined for DNA methylation status using combined bisulphite restriction analysis indicated by the arrowheads. The intensity of the bands marked with the white and black circles show the amount of unmethylated and methylated CpG, respectively.

(B) Repetitive experiments and quantification of the cut and uncut bands of combined bisulphite restriction analysis were performed for the most 3′ site indicated in Figure 4A. One black and two grey bars represent the results of positive control cells and of two independent 5-aza-2′-deoxycytidine treated cells, respectively. The decrease in methylation was statistically significant (**) (p < 0.01). Quantification of each samples was performed three times using independent PCR products.

(C) Relative expression was calculated by quantifying the results of restriction fragment length polymorphism analysis. White and black bars represent the expression from active and inactive alleles of positive control cells, respectively. It should be noted that expression from the inactive alleles was negligible. Two grey bars represent induced expression from the inactive alleles of two independent 5-aza-2′-deoxycytidine treated cells. Statistically significant increase in expression after treatment is shown by * (p < 0.03) or ** (p < 0.01).

https://doi.org/10.1371/journal.pgen.0030055.g004

The DNA methylation status of retrotransposons can differ between male and female germ cells. For example, IAP and LINE1 are more highly methylated in sperm than oocytes, while Alu is less methylated [12]. Mice and humans with paternal disomy that express PEG10 and SGCE biallelically have normal phenotypes, so monoallelic expression of these genes is not essential for development. Therefore, although PEG10 is essential for placental development in the mouse, the imprinting of this locus may be a functionally unimportant inheritance derived from the nature of the original retrotransposition of PEG10.

There were CpG islands in the putative promoter region of SGCE of the chicken, platypus, tammar, mouse, and human (Figure 5A). We hypothesize that insertion of PEG10 after the divergence of therian from prototherian mammals expanded the CpG islands (Figure 5). In the tammar, DNA methylation is restricted to PEG10, but in the mouse and human, the entire region is differentially methylated. These differences may be explained by the presence or absence of a boundary function of the CpG island in these groups as discussed above. However, in both cases, insertion of PEG10, which must have occurred in the therian ancestor, is clearly sufficient to establish imprinting of this region (Figure 5B). Our study confirms that silencing of exogenous DNA after retrotransposon insertion can drive the evolution of genomic imprinting in mammals.

thumbnail
Figure 5. Comparisons of the Genomic Region Surrounding SGCE in Various Species and Phylogeny of PEG10 Imprinting

(A) CpG content of a 1-kb genomic sequence starting in the conserved end of the exon1 in SGCE is shown in each species. The regions corresponding to SGCE and PEG10 are represented by boxes below. Yellow boxes indicate biallelic expression and blue boxes represent paternal expression. Precise start sites of SGCE in chicken and platypus remain unknown.

(B) Illustration of the insertion of the original PEG10 sequence in therian not prototherian lineage at least 125 million years ago [31,32]. The short and long boxes represent SGCE and PEG10, respectively. The black lollipops indicate DNA methylation of CpG sites and the grey boxes represent the silent state of transcription.

https://doi.org/10.1371/journal.pgen.0030055.g005

Materials and Methods

Animals and tissue collection.

Tammar wallabies of Kangaroo Island origin were maintained in our breeding colony in grassy, outdoor enclosures. Lucerne cubes, grass and water were provided ad libitum and supplemented with fresh vegetables. Fetuses and yolk sac placenta tissue were collected between days 22 and 25 of the 26.5-d gestation [29]. Experimental procedures conformed to Australian National Health and Medical Research Council (1990) guidelines and were approved by the Animal Experimentation Ethics Committees of the University of Melbourne.

Isolation of BAC clones and determination of genomic DNA sequences.

Each BAC clone in the tammar and platypus BAC libraries was stored separately and was spotted onto nylon membranes correspondently. These membranes were hybridized with the partial SGCE probes of each species, and the positive clones were identified according to the locus information of the signals on the membranes. The tammar PEG10 sequence was determined by the primer walking method from a partial fragment amplified using cross-species degenerate primers. The platypus sequence between SGCE and PPP1R9A was determined by the shotgun sequencing method, and it was completed using a published database of whole genome shotgun sequences and direct sequencing of PCR products.

Detection of repetitive sequences.

RepeatMasker (http://www.repeatmasker.org) was used for the detection of LINEs, SINEs, and LTR elements in the genomic region between SGCE and PPP1R9A.

Fluorescence in situ hybridization.

BAC DNA was labeled by nick translation with digoxygenin-11-dUTP. Hybridization was operated with the labeled BAC DNA and C0t-1 DNA at 37 °C overnight. Anti-digoxygenin-Cy3 and DAPI were used for the detection of the signals and for the counterstain, respectively.

Single nucleotide primer extension assay.

The details have been described in our previously published paper [8]. The 3′ UTRs of PEG10 and ASB4 including the polymorphisms were amplified by 30–35 cycles of RT-PCR using the following pairs of primers: PEG10-F1, 5′-CAAATGCCATTGCCGTCT-3′ and PEG10-R1, 5′-GTTAGACGGTCAGCTCCACG-3′; PEG10-F2, 5′-CAACCAGGGGAGCTAGGATT-3′ and PEG10-R2, 5′-GAACATCCATGCACCGTAGA-3′; PEG10-F3, 5′-CTCTCTGGAGCGGTATCCAG-3′ and PEG10-R3, 5′-TGTGAGATTTGGCAATCATACA-3′; and ASB4-F, 5′-AACACCCCGAGGTCTCTCAT-3′ and ASB4-R, 5′-GAGGACCATGGCATTTATTCA-3′. The following correspondent SNuPE primers were used for the single nucleotide extension: PEG10-SN1, 5′-ATTCATTTCCCTTCCCAACAT-3′; PEG10-SN2, 5′-CCCTGGCTGCGAGACCA-3′; PEG10-SN3, 5′-CCCGGGGAGCTCCGAGC-3′; and ASB4-SN, 5′-CAAGAAGCAAGTAGTTCTTCAAAGG-3′

Hot stop PCR.

The labeling of the primers for the last hot cycle was operated using [γ-32P]ATP and T4 DNA kinase. The final PCR products of SGCE and PPP1R9A were digested by Alul and Mspl, respectively. The 3′ UTRs of SGCE and PPP1R9A including the polymorphisms on the recognition sequences of these restriction enzymes were amplified by 30–35 cycles of PCR using the following pairs of primers. Asterisks indicate the labeled primers: SGCE-F, 5′-CAGTGATGGCGTTCTGTACG-3′; *SGCE-R, 5′-GTTGATGACCAGGTTGTGCC-3′; *PPP1R9A-F, 5′-CCAGGAGAAGATGGAGAAGC-3′; and PPP1R9A-R, 5′-GTTGGGGATGAAGGAGTGTG-3′. Ten percent polyacrylamide gels were used for the gel electrophoresis of the digested samples.

Bisulphite sequencing.

After the bisulphite treatment [30] for the genomic DNA of tammar, the region corresponding to the CpG islands over the promoter regions of SGCE and PEG10 was amplified by 35 cycles of PCR using the following pair of primers: CGI-F, 5′-GGAGTGATTGTGGAAATGGAGGTG-3′ and CGI-R, 5′-ATACAAAATCCCCCCCCTAAACCTC-3′. The PCR products were cloned and the clones were analyzed by sequencing.

Cell culture.

Primary culture of tammar fetal lung cells from day 25 of gestation and adult endometrium cells were used in this study. Control cells were cultured in 50% AmnioMAX (Invitrogen, http://www.invitrogen.com) and 50% DMEM supplemented with 10% fetal calf serum and penicillin/streptomycin at 37 °C/5% CO2. Cells for 5-aza-2′-deoxycytidine treatment were cultured in the same media but containing 10 μM of 5-aza-2′-deoxycytidine (Sigma). Fresh media with 5-aza-2′-deoxycytidine were added every 24 h for 6 d.

Combined bisulphite restriction analysis.

Three regions in the PEG10 DMR were amplified by 35 (for the middle and right side in Figure 4A) and 40 (for the left side in Figure 4A) cycles of PCR from the bisulphite treated genomic DNA of cells using the following pair of primers: LEFT-F, 5′-GTATTAGTTTTTTTGTAGTT-3′ and LEFT-R, 5′-CCTAAAAAACTACCCTACTCC-3′; MID-F, 5′- GAGATGGGGAGATTGATATTT-3′ and MID-R, 5′- CCCTATAACTAAACTACAATCTCTCC-3′; and RIGHT-F, 5′- CCTCCCATTAACTTTAAAATCACC-3′ and RIGHT-R, 5′- ATTGTAGTAATGGGGTAGGTTATG-3′.

PCR products were digested by RsaI (for the left side) or Aci I (for the middle and right side) for analyses.

Restriction fragment length polymorphism analysis.

The 3′ UTR of PEG10, including the polymorphism on the TaqI or BceAI recognition sequences, was amplified by 30 cycles of RT-PCR using the PEG10-F1R1 or PEG10-F3R3 primer pairs. PCR products were digested by TaqI or BceAI for analyses.

Supporting Information

Figure S1. Amino Acid Sequence Alignment between Human and Tammar PEG10

Amino acid sequence identity (asterisks), homology (dots), and divergence (no marks) between human and tammar PEG10 (ORF1 and 2 are combined) are shown. The red boxes represent the conserved CCHC and DSG motifs.

https://doi.org/10.1371/journal.pgen.0030055.sg001

(15 KB PDF)

Accession Numbers

The National Center for Biotechnology Information GenBank (http://www.ncbi.nlm.nih.gov/Genbank) sequence accession numbers for tammar and platypus BACs are AB260975 and AB260976, respectively.

Acknowledgments

We thank Arizona Genomic Institute and Clemson University Genomic Institute for tammar and platypus BAC clones. Animals were held in our research colony under permits 1000450 and 10001910 from the Victorian Department of Sustainability and Environment. We thank Sue Osborn for establishing the fetal lung cell line and Kerry Martin, Scott Brownlees, and members of the wallaby research group for assistance in animal husbandry and handling.

Author Contributions

SS, RO, TK, FI, and TKI conceived and designed the experiments. SS performed the experiments. SS, TN, and CW analyzed the data. TN, AJP, GS, AEA, JAMG, YK, and MBR contributed reagents/materials/analysis tools. SS, AJP, GS, FI, and MBR wrote the paper.

References

  1. 1. Bartolomei MS, Zemel S, Tilghman SM (1991) Parental imprinting of the mouse H19 gene. Nature 351: 153–155.
  2. 2. Ferguson-Smith AC, Sasaki H, Cattanach BM, Surani MA (1993) Parental-origin-specific epigenetic modification of the mouse H19 gene. Nature 362: 751–755.
  3. 3. Li E, Beard C, Jaenisch R (1993) Role for DNA methylation in genomic imprinting. Nature 366: 362–365.
  4. 4. Kaneko-Ishino T, Kohda T, Ono R, Ishino F (2006) Complementation hypothesis: The necessity of a monoallelic gene expression mechanism in mammalian development. Cytogenet Genome Res 113: 24–30.
  5. 5. Lewis A, Reik W (2006) How imprinting centres work. Cytogenet Genome Res 113: 81–89.
  6. 6. Killian JK, Byrd JC, Jirtle JV, Munday BL, Stoskopf MK, et al. (2000) M6P/IGF2R imprinting evolution in mammals. Mol Cell 5: 707–716.
  7. 7. O'Neill MJ, Ingram RS, Vrana PB, Tilghman SM (2000) Allelic expression of IGF2 in marsupials and birds. Dev Genes Evol 210: 18–20.
  8. 8. Suzuki S, Renfree MB, Pask AJ, Shaw G, Kobayashi S, et al. (2005) Genomic imprinting of IGF2, p57(KIP2) and PEG1/MEST in a marsupial, the tammar wallaby. Mech Dev 122: 213–222.
  9. 9. Lewis A, Mitsuya K, Umlauf D, Smith P, Dean W, et al. (2004) Imprinting on distal Chromosome 7 in the placenta involves repressive histone methylation independent of DNA methylation. Nat Genet 36: 1291–1295.
  10. 10. Barlow DP (1993) Methylation and imprinting: from host defense to gene regulation? Science 260: 309–310.
  11. 11. Ono R, Kobayashi S, Wagatsuma H, Aisaka K, Kohda T, et al. (2001) A retrotransposon-derived gene, PEG10, is a novel imprinted gene located on human Chromosome 7q21. Genomics 73: 232–237.
  12. 12. Yoder JA, Walsh CP, Bestor TH (1997) Cytosine methylation and the ecology of intragenomic parasites. Trends Genet 13: 335–340.
  13. 13. Sapienza C, Peterson AC, Rossant J, Balling R (1987) Degree of methylation of transgenes is dependent on gamete of origin. Nature 328: 251–254.
  14. 14. Surani MA, Reik W, Allen ND (1988) Transgenes as molecular probes for genomic imprinting. Trends Genet 4: 59–62.
  15. 15. Chaillet JR, Vogt TF, Beier DR, Leder P (1991) Parental-specific methylation of an imprinted transgene is established during gametogenesis and progressively changes during embryogenesis. Cell 66: 77–83.
  16. 16. Brandt J, Schrauth S, Veith AM, Froschauer A, Haneke T, et al. (2005) Transposable elements as a source of genetic innovation: expression and evolution of a family of retrotransposon-derived neogenes in mammals. Gene 345: 101–111.
  17. 17. Youngson NA, Kocialkowski S, Peel N, Ferguson-Smith AC (2005) A small family of sushi-class retrotransposon-derived genes in mammals and their relation to genomic imprinting. J Mol Evol 61: 481–490.
  18. 18. Ono R, Nakamura K, Inoue K, Naruse M, Usami T, et al. (2006) Deletion of Peg10, an imprinted gene acquired from a retrotransposon, causes early embryonic lethality. Nat Genet 38: 101–106.
  19. 19. Margulies EH, Maduro VV, Thomas PJ, Tomkins JP, Amemiya CT, et al. (2005) Comparative sequencing provides insights about the structure and conservation of marsupial and monotreme genomes. Proc Natl Acad Sci U S A 102: 3354–3359.
  20. 20. Gould SJ, Vrba ES (1982) Exaptation–a missing term in the science of form. Palaeobiology 8: 4–15.
  21. 21. Brosius J, Gould SJ (1992) On “genomenclature”: A comprehensive (and respectful) taxonomy for pseudogenes and other “junk DNA.”. Proc Natl Acad Sci U S A 89: 10706–10710.
  22. 22. Piras G, El Kharroubi A, Kozlov S, Escalante-Alcalde D, Hernandez L, et al. (2000) Zac1 (Lot1), a potential tumor suppressor gene, and the gene for epsilon-sarcoglycan are maternally imprinted genes: Identification by a subtractive screen of novel uniparental fibroblast lines. Mol Cell Biol 20: 3308–3315.
  23. 23. Ono R, Shiura H, Aburatani H, Kohda T, Kaneko-Ishino T, et al. (2003) Identification of a large novel imprinted gene cluster on mouse proximal Chromosome 6. Genome Res 13: 1696–1705.
  24. 24. Mizuno Y, Sotomaru Y, Katsuzawa Y, Kono T, Meguro M, et al. (2002) Asb4, Ata3, and Dcn are novel imprinted genes identified by high-throughput screening using RIKEN cDNA microarray. Biochem Biophys Res Commun 290: 1499–1505.
  25. 25. Christy RJ, Huang RC (1988) Functional analysis of the long terminal repeats of intracisternal A-particle genes: Sequences within the U3 region determine both the efficiency and direction of promoter activity. Mol Cell Biol 8: 1093–1102.
  26. 26. Minakami R, Kurose K, Etoh K, Furuhata Y, Hattori M, et al. (1992) Identification of an internal cis-element essential for the human L1 transcription and a nuclear factor(s) binding to the element. Nucleic Acids Res 20: 3139–3145.
  27. 27. Roseman RR, Pirrotta V, Geyer PK (1993) The SU(HW) protein insulates expression of the Drosophila melanogaster white gene from chromosomal position-effects. Embo J 12: 435–442.
  28. 28. Donze D, Adams CR, Rine J, Kamakaka RT (1999) The boundaries of the silenced HMR domain in Saccharomyces cerevisiae. Genes Dev 13: 698–708.
  29. 29. Renfree MB (1973) The composition of fetal fluids of the marsupial, Macropus eugenii. Dev Biol 33: 62–79.
  30. 30. Frommer M, McDonald LE, Millar DS, Collis CM, Watt F, et al. (1992) A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. Proc Natl Acad Sci U S A 89: 1827–1831.
  31. 31. Luo ZX, Ji Q, Wible JR, Yuan CX (2003) An Early Cretaceous tribosphenic mammal and metatherian evolution. Science 302: 1934–1940.
  32. 32. Asher RJ, Horovitz I, Sanchez-Villagra MR (2004) First combined cladistic analysis of marsupial mammal interrelationships. Mol Phylogenet Evol 33: 240–250.