Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Mosaic Retroposon Insertion Patterns in Placental Mammals

Gennady Churakov1*, Jan Ole Kriegs1,2*, Robert Baertsch3, Anja Zemann1, Jürgen Brosius1 and Jürgen Schmitz1

1Institute of Experimental Pathology, ZMBE, University of Münster, Von-Esmarch-Str. 56,

48149 Münster, Germany

2current address: LWL-Museum für Naturkunde, Sentruper Str. 285, 48161 Münster,

Germany

3Department of Biomolecular Engineering, University of California, Santa Cruz, CA 95064,

USA

*These authors contributed equally.

Corresponding authors:

Jürgen Schmitz ([email protected]) and

Gennady Churakov ([email protected])

Institute of Experimental Pathology (ZMBE)

University of Münster, Von-Esmarch-Str. 56, 48149 Münster, Germany

Tel +49 251 8352133

Fax +49 251 8352134

Running title: Basal Trifurcation of Placental Mammals

Keywords: retroelements, ancestral hybridization, lineage sorting, mammalian phylogeny, placental origin, soft polytomy

1 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Abstract

One and a half centuries after Charles Darwin and Alfred Russel Wallace outlined our current understanding of evolution, a new scientific era is dawning that enables direct observations of genetic variation. However, pure sequence-based molecular attempts to resolve the basal origin of placental mammals has so far resulted only in apparently conflicting hypotheses. By contrast, in the mammalian genomes where they were highly active, the insertion of retroelements and their comparative insertion patterns constitute a neutral, virtually homoplasy-free archive of evolutionary histories. The “presence” of a retroelement at an orthologous genomic position in two species indicates their common ancestry in contrast to its

“absence” in more distant species. To resolve the placental origin controversy we extracted approximately 2 million potentially phylogenetic informative, retroposon-containing loci from representatives of the major placental mammalian lineages and found highly significant evidence challenging all current single hypotheses of their basal origin. The Exafroplacentalia hypothesis (Afrotheria as the sister group to all remaining placentals) is significantly supported by five retroposon insertions, the Epitheria (Xenarthra as the sister group to all remaining placentals) by nine insertion patterns, and the Atlantogenata hypothesis (a monophyletic clade comprising Xenarthra and Afrotheria as the sistergroup to Boreotheria comprising all remaining placentals) by eight insertion patterns. These findings provide significant support for a ‘soft’ polytomy of the major mammalian clades. Ancestral successive hybridization events and/or incomplete lineage sorting associated with short speciation intervals are viable explanations for the mosaic retroposon insertion patterns of recent placental mammals and for the futile search for a clear root dichotomy.

2 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Introduction

Genomic variation is based on genetic convertibility and is amplified by gene flow and sexual recombination. While only a small fraction of genetic variation is directly exposed to the natural selection driving the evolution of species the majority of changes are neutral (Kimura

1968). In the Molecular Darwinian and post-genomic age when such variations are now directly observable in prodigious numbers, one still needs only compare the volumes of controversy in scientific reports to dispel the notion that all the mysteries of evolutionary history are now an open book. The basal node of the placental mammals is a prime example of such a controversy. After continuous morphological inconsistencies and a lack of resolution due to narrow temporal speciation events at the deeper divergences of the placental mammalian tree, the hope was that large-scale DNA sequencing and an increasing number of sampled species might resolve the higher-level mammalian relationships by reducing systematic biases. But the extraction and analysis of molecular traces of evolutionary history has been no less sensitive to misinterpretation than was the understanding of dusty specimens

150 years ago. The reliability of sequence data in genome-scale phylogenetic approaches is subject to unequal evolutionary rates among lineages, regional- or lineage-specific compositional biases, shifts in site-specific evolutionary rates over time and sequence regions

(see for example Nishihara et al. 2007), variable accuracies of multi-sequence alignments and analytical tools, and last but not least, the general liability of such data to homoplasy.

In a seminal work, Murphy et al. (2001a) investigated approximately 10 kb of sequence information to substantiate the previously proposed (Waddell et al. 1999) major mammalian groups Afrotheria, Xenarthra, Laurasiatheria and Euarchontoglires. Moderate support was found for Afrotheria as the basal split of placentals (Exafroplacentalia;

Figure 1A). A subsequent combination and expansion of previously published data (Madsen et al. 2001; Murphy et al. 2001a) found significant support for a basal split between

Afrotheria and other placentals (Murphy et al. 2001b). Also Nikolaev et al. (2007) found

3 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

support for the Afrotheria as the first placental split drawing on approximately 200 kb of protein-coding sequences from the 1% of the human genome studied in the ENCODE pilot project. However, other studies of large-scale sequences (e.g., Hallström et al. 2007) and previous studies of smaller data sets could not significantly confirm an Afrotherian root

(e.g., Amrine-Madsen et al. 2003; Delsuc et al. 2002; Madsen et al. 2001; Waddell et al.

2001). A sister group relationship of Afrotheria and Xenarthra (Atlantogenata) also received significant support in analyses of sequences from approximately 1,700 conserved genome loci

(Wildman et al. 2007), 2840 protein-coding genes (Hallström et al. 2007), mitochondrial genomes (Kjer and Honeycutt 2007) and was recently confirmed by large-scale analyses of

2.8 Mbp of protein-coding data (Hallström and Janke 2008) and large-scale genomic sequences (Prasad et al. 2008). Thus, despite many different molecular studies, the basal origin of placental mammals remains controversial and unresolved. Unfortunately, pure sequence data cannot resolve these apparent contradictions.

By contrast, the presence/absence patterns of inserted retroelements constitute a virtually homoplasy-free marker system with theoretically infinite character states (Steel and

Penny 2000). Specific genomic insertions of such elements in the ancestor of two species reliably document their common ancestry. The very few known examples of discordance in retroelement presence/absence data can be explained by deletion via illegitimate recombination between perfect direct repeats flanking each insertion (van de Lagemaat et al.

2005), exact parallel insertions (Cantrell et al. 2001), or lineage sorting, a phenomenon related to incomplete allele fixation and frequent intervals of speciation events (Ray et al. 2006;

Shedlock et al. 2004).

Analyses of a smaller number of retroelements as phylogenetic markers in mammals validated the four superordinal mammalian clades (Kriegs et al. 2006; Möller-Krull et al.

2007; Nishihara et al. 2005). However, in support of Shoshani and McKenna (1998), but in contradiction to most other molecular investigations, we found the first evidence, two L1MB5

4 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

retroelements, that placed Xenarthra (Epitheria hypothesis) (Figure 1B) instead of Afrotheria at the base of the placental tree (Kriegs et al. 2006). Meanwhile, Murphy et al. (2007) presented two L1MB5 retroelements and four additional indels that favored a sister group relationship of Afrotheria and Xenarthra (Atlantogenata) (Figure 1C) at the base of the placental tree. Thus, to date even retroposon insertion patterns have not satisfactorily resolved the basal origin of placental mammals.

With the intention of resolving the conflicting positions of placental lineages, we aimed to examine a much larger sample of retroposed elements and in more species. For this study we took a slightly different tack, and aligned whole genomes of three representative placental species, including additional species for informative loci to represent a maximum divergence within each of the groups, and scanned them for diagnostic retroelement insertions. In an analysis of approximately 2 million potential phylogenetically informative loci, we found multiple retroposon insertion patterns significantly supporting all three placental speciation hypotheses, indicating a complex ancestral speciation scenario including successive early divergences in close temporal proximity and hybridization and/or lineage sorting events as viable sources for an effective ‘soft’ polytomy. This very rare example of multiple retroposon incongruence goes a long way towards explaining the decades-long stream of apparently conflicting evidence for one or the other placental evolutionary history based on multiple marker systems including both morphological and molecular sequence evidence.

5 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Results

We focused our computational screening for phylogenetically informative markers on the insertion patterns of L1MB elements that were active during the critical period of early placental divergence ~100 million years ago (Kriegs et al. 2006; Murphy et al. 2007), and independently tested the three viable scenarios of placental origin: the Exafroplacentalia,

Epitheria, and Atlantogenata hypotheses (Fig. 1). We screened the first of the two whole genome 3-way alignments (armadillo-elephant-human) for elements supporting the

Exafroplacentalia (those present in armadillo and human but absent in elephant) and

Atlantogenata (those present in armadillo and elephant but absent in human) hypotheses, and the second alignment (elephant-human-armadillo) for elements supporting the Epitheria hypothesis (those present in elephant and human but absent in armadillo; Figure 2A). Using the criteria presented in Methods, we computationally pre-selected 1227 candidate loci and manually re-inspected them using the UCSC sequence browser in combination with sequences from the trace archives. Manual screening of these loci revealed 22 conserved, potentially informative loci that we then re-aligned to verify the identities and boundaries of the inserted retroelements and to clearly characterize their respective flanking regions. To validate the presence/absence patterns of retroposed elements in these 22 loci in an expanded species sampling and to present as consistent a species representation for each marker as possible, we then retrieved orthologous loci in two distant supraprimates (human and mouse or guinea pig, rabbit), two laurasiatherians (dog and horse or cow), two Afrotheria (elephant and tenrec or hyrax), two Xenarthra (armadillo and sloth), and, as far as alignable, opossum as the outgroup. In some cases sufficient sequence data were not yet available from members of the Afrotheria clade (tenrec) or Xenarthra (sloth) to positively confirm the presence and or absence of all retroposons in a second representative species of the given clades. Some of those gaps were filled by experimental PCR amplification and sequencing of the corresponding loci in tenrec and sloth. Of the 22 phylogenetically informative markers

6 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

analyzed in this study, 20 are highly reliable in that they provide recognizable information of potential direct repeats plus the corresponding unoccupied target site in the absence cases (for examples see Figure 3). The remaining two markers are classified as “good” with some restrictions regarding the recognition of clear direct repeats. Their broad distribution over 15 of the 23 human chromosomal pairs indicates the independent integration of the 22 phylogenetically informative elements (Supplemental Fig. S1).

In strong contrast to most of the mammalian branches that we and others investigated previously with this method (e.g., Huchon et al. 2007; Kriegs et al. 2007; Möller-Krull et al.

2007; Nishihara et al. 2005; Nishihara et al. 2006; Warren et al. 2008; Xing et al. 2007), we found a rare example of virtually equal support for three competing hypotheses (Figure 2B,

Supplemental Table S1, Kriegs J.O. (2007) Dissertation, Univ. Muenster, Germany). We found five independent L1 insertions (2 L1MB4, 2 L1MB5, 1 L1MB7) shared by Boreotheria

(supraprimates and laurasiatherians) and Xenarthra, but absent in Afrotheria, that would support the Exafroplacentalia hypothesis with Afrotheria as the sister group to all remaining placentals. Our search returned nine (including one previously reported by Kriegs et al.

(2006)), independent insertions of L1 elements (7 L1MB5, 1 L1MB8, 1 L1MB4) shared by

Boreotheria and Afrotheria but absent in Xenarthra. Together with an additional independent marker reported by Kriegs et al. (2006) there are now ten retroposon markers that would support the Epitheria hypothesis. We also found eight (two of which were previously reported by Murphy et al. (2007) L1 insertions (7 L1MB5 and 1 L1MB8) that would support the

Atlantogenata hypothesis, as they were present in Afrotheria and Xenarthra species and clearly absent in those of Boreotheria. Importantly, for each branching scenario we found 5, 9, and 8 clear markers that would individually translate to probabilities of P<0.0001 that the individual hypotheses (by screening for markers in just one direction) are incorrect ([5 0 0], [9

0 0], [8 0 0] (Waddell et al. 2001).

7 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Of all the phylogenetically informative mammalian markers investigated to date

(Kriegs et al. 2006; Nishihara et al. 2006; Murphy et al. 2007), just the one branch leading to the major placental clades has exhibited such extreme contradictions. The only hypothesis to best explain these apparently conflicting results is that there is no clear dichotomy at the base of placentals, but rather a trifurcation.

8 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Discussion

Comparative analyses of large-scale genome data afford a new era of phylogenetic inference combining the fields of evolution and genomics into one of phylogenomics. A completely objective history of the mammalian orders is one of the main visions of this new field of phylogenomics. Most of the interordinal mammalian branches have already been confirmed with convincing support from genome data, but problematic areas of the placental tree have also been recognized (Hallström and Janke 2008). As an example, narrow ancestral splitting events such as the early divergence of placentals have been contradictorily resolved by genome-wide analyses of protein-coding data (Hallström et al. 2007; Nikolaev et al. 2007).

But, because nucleotide sequences are highly exposed to homoplasy, using them as phylogenetic markers cannot resolve these contradictions.

Our current data, with significant, nearly equivalent support for all three hypotheses in one study, could be explained by three possible scenarios: deletion via illegitimate recombination between perfect direct repeats flanking each insertion, exact parallel insertions, or lineage sorting, a phenomenon related to incomplete allele fixation and frequent intervals of speciation events (Ray et al. 2006; Shedlock et al. 2004). Precise deletion was observed in

0.5% of LINE1-derived SINE insertions in primates (van de Lagemaat et al. 2005).

Statistically speaking, this could affect at most six of our 1227 candidate loci, and the chance of all six being in the 22 markers used for the final analysis is infinitesimally small.

Moreover, the process of precise deletion is restricted to a very small period of time where both element flanking direct repeats must be identical to allow recombination. Lagemaat et al.

(2005) also did not consider the effects of lineage sorting or mosaic evolution in his calculation, hence the expected number of real precise deletions should be much lower than

0.5%. From a previous investigation (Ray et al. 2006) it is known that precise parallel insertion occurs in 0.05% of observed primate LINE1-derived SINE insertions. Considering our total of 1227 candidate loci, this is less than one. Furthermore, all previous investigations

9 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

did not consider random truncations that are an essential, additional factor for characterizing

LINE element orthology.

Therefore, as the homoplasy of retroposon insertion patterns is extremely rare in mammalians and was only described recently for one apparent conflicting marker challenging the newly described Pegasoferae clade (Nishihara et al. 2006), this leaves us with the last explanation, lineage sorting. We feel that the most parsimonious interpretation of the current data is that the ancestral placental populations were characterized by severe ancestral subdivisions and re-joinings, leading to a complex mosaic of phylogenetic relationships in recent species (Figure 4). Effects of alternating divergence, hybridization, introgression, and lineage sorting might complicate our search for a clear dichotomy at the base of this tree and leave us with an indistinct, effective ‘soft’ polytomy, leading sometimes to one or the other solution depending on the size of the data set and the particular markers examined.

Two different scenarios or a mixture of both might be responsible for the mosaic pattern that we found in our screening for informative retroposon insertions and that other groups found in phylogenetic investigations of subsets of genomic sequence data for these species. The earliest eutherian mammals (extant placentals plus closely related extinct mammals) originated about 125 million years ago in the Early Cretaceous, and based on appearance can probably be represented by the recently discovered and reconstructed skeletal remains of Eomaia scansoria (Ji et al. 2002). In the first scenario, the Eomaia-like ancestral population probably diverged into three distinct lineages of preAfrotheria, preXenarthra, and preBoreotheria. Limited gene flow and hybridization that occurred among these populations left behind mosaic signs of intermittent relationships that we identify 100 million years later as a patchwork of Exafroplacentalia (I), Epitheria (II) and Atlantogenata (III) roots (Figure

4A). This scenario requires temporary overlap of the three lineages or parts of their populations before final speciation occurred (Figure 4B). A second scenario that could explain the results involves incomplete lineage sorting among the three populations (Figure

10 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

4C). Because retroposon fixation requires possibly millions of years before being consistently represented in a population (Schmitz and Zischler 2002), alleles with and without certain elements surely coexisted temporarily, and, in the case of rapid speciation, were distributed randomly into different lineages. Such a process could contribute as well to the patchwork presentation of seemingly competing insertion patterns now present in recent species.

A mosaic of evolutionary signals, blueprinted in the human genome, was also recently described by Ebersberger et al. (2007). They reported that 23% of the genomic sequences

(equal representations of genes and intergenic regions; 23,210 alignments, each ≥300 nt) that they analyzed represented an incongruent genealogy in which chimpanzee was not the closest relative to human.

The majority of sequence-based investigations aimed at resolving the placental root, by examining numerous genes, sequence partition, and various analytic tools, also found almost equal support for one or the other of the various hypotheses (e.g., Delsuc et al. 2002;

Madsen et al. 2001; Waddell and Shelley 2003); once more emphasizing the notion that early placental divergence constituted a patchwork evolution of the initial lineages.

That such conflicting results are so very rare requires some discussion of the actual markers that we recovered. The 22 phylogenetically informative markers providing significant support for a trifurcation at the basal node of placental mammals that were retrieved in the current screen included both of the previously described Atlantogenata markers (Murphy et al. 2007), but only one of the two previously described Epitheria markers

(Kriegs et al. 2006). Although the undetected one was not present in the 3-way MAF alignments, we do consider it to be a reliable marker (Supplemental Figure S2). A distinct advantage of the current expanded screening compared to previous searches is a more unbiased evaluation of equal amounts of data for all three hypothetical scenarios of placental origin based on whole-genome alignments. We based our previous strategies (Kriegs et al.

2006; Strategy II) on a random search of human intron data and screens of the limited 11 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

amounts of trace sequences of elephant and armadillo available at the time, while the specifically designed, whole-genome, 3-way alignments currently employed offer a much more comprehensive data source resulting in detection of a high number of potentially informative markers. However, the fact that we recovered only three of the four previously described retroposon markers for early placental splits indicates that the current limitations of an exhaustive and maximized unbiased search is the naturally imperfect and variable quality of available specific whole-genome alignments.

One critical point in using retroposon insertion data, especially in deep phylogenetic splits, is the continuous divergence of the elements and their flanking regions. Nevertheless, we recently demonstrated that retroposon data could provide reliable markers in deep mammalian phylogeny by supporting the more than 130 million year old therian split between marsupials and placentals (Warren et al. 2008). And as stated previously, 20 of the retroposons examined here had clearly recognizable information concerning their direct repeats and the other two were, though recognizable, were not quite as clear.

Direct repeats generated by target site duplications in the process of retroposition are valuable landmarks of an orthologous insertion but are unlikely to remain unchanged throughout deep phylogenetic splits. Hence, there is an additional and more reliable orthology criterion available, particularly for the LINE-related retroposons such as L1MB elements.

These elements insert almost exclusively as randomly truncated forms. Thus, as independent insertions of identically truncated LINE elements of the same subtype at the same genomic locus of two species is highly unlikely, determining the exact point of truncation of the observed orthologous elements will clearly trace them back to a common origin. A genome- wide screening for the 20,700 LINE1 (L1MB) truncations demonstrated a significant random distribution of truncations with no bias for specific truncation points (P<0.002, Supplemental

Figure S3), and all investigated diagnostic LINE1 elements differ in their lengths

12 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

(Supplemental Figure S3). Thus, retroposon truncation points serve as an additional, powerful criterion for orthologous insertions.

The randomness of the insertions of LINE1-related elements is well demonstrated by

Ray et al. (2006) with only a very few cases of observed conflicts. A genome-wide investigation of target site preferences indicates a slight preference for TT/AAAAA sites

(Supplemental Figure S4). This is well known as a kinkable DNA site (Jurka et al. 1998).

However, the full target site duplication is much longer (8-30 nts) and rather individual.

Nevertheless, the multiplicity of criteria, including (1) identical genomic insertion points, (2) element orientation, (3) identical LINE1 subtypes, (4) matching length and sequence of direct repeats, (5) concurrent truncation points, and (6) the complete consistency in diverse members of representative orders, make it improbable to have screened for artifacts.

Investigations of very young retroposon insertion events present another problem, lineage specific exact or nearly exact deletions. For young retroposon insertions flanked by perfect identical direct repeats, precise removal might be possible via recombination involving the direct repeats (van de Lagemaat et al. 2005). Because direct repeats diversify quickly after insertion, this possibility should not be critical for deep phylogenetic splits (after altering of perfect direct repeats), but it is not completely impossible for such random perfect or nearly perfect deletions of retroelements to have accumulated, especially in deep divergences where exact comparisons of diverged direct repeats are difficult. However, the tendency of LINE1-related elements to integrate within A/T-rich sequence regions could promote such deletions. Comparing potential retroposon absence regions with a clear absence in an outgroup species is one way to diminish such potential misinterpretations.

Unfortunately, for placentals the next outgroup are marsupials that diverged substantially.

Due to the evolutionary distance between marsupials and placentals, even if possible, alignments between these species are, particularly for intergenic regions, only approximate and mostly not a critical proof of true absence. On the other hand, the outgroup-comparison

13 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

for informative L1MB elements is not essential, because L1MB elements are placental- specific and therefore necessarily absent in any outgroup species. Thus, the suggestion of

Murphy et al. (2007) that one of our informative L1MB5 markers for Epitheria (Kriegs et al.

2006) is traceable in opossum may be based on a misalignment and misannotation of the corresponding BLAT browser location (Supplemental Figure S2).

Conclusions

Comparing the genome-wide insertion patterns of retroelements is a powerful strategy to gain a virtually homoplasy-free picture of the evolutionary histories of species. Virtually homoplasy-free implies that this marker system is, as any other marker-system, not insensitive to the effects of ancestral hybridization and incomplete lineage sorting combined with short splitting times. Because of the clear character polarity (presence of an element as the derived state) and the low probability of orthologous exact parallel insertions or deletions, the retroposon presence/absence marker system is exceptionally effective for detecting the ancestral exchange of genetic material among lineages and the lineage sorting effects that are difficult to pinpoint by the statistical evaluation of heterogeneous sequence data (Hallström and Janke 2008). Such effects are more difficult to recognize in deep phylogenetic splits and require maximized, unbiased, multidirectional screening of whole-genome phylogenetic signals. They must also be conducted with the notion in mind, that in such rare cases, the minimum of three clear markers proposed by Waddell et al. (2001) is not absolutely sufficient as a significant clade support. The present example of apparently incongruent markers inherent in the early branching of placentals offers a seminal example of conflicting retroposon presence/absence patterns that, at the same time, shed new light on the decades old, controversially discussed scenarios of sequence-based phylogenies. As demonstrated, retroposon screening offers a powerful test of critical signals in narrow splits to detect nodes with more than two immediate descending branches. The basal rodent splits, the Euarchonta

14 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

relationships, and the affiliations among laurasiatherians are only a few additional contentious examples that may possibly find elucidation by comparing the insertion histories of genome- wide, multidirectionally extracted, phylogenetically informative retroelements.

15 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Materials and Methods

Genome alignments

To perform an exhaustive, genome-wide screening for phylogenetically informative retroposon insertions in placentals and to test the three viable phylogenetic hypotheses of the basal origin of placentals, we derived two different 3-way (species) whole-genome sequence alignments in MAF format (UCSC – University of California Santa Cruz): (1) armadillo - elephant - human (3,050 giga bases) and (2) elephant – human – armadillo (2,730 giga bases), the decisive difference between the two alignments being the full retroposon representation of the leading species that was used as the profile for the remaining species.

The 3-way alignments where prepared using pairwise (armadillo-elephant, armadillo- human for the first 3-way alignment and elephant-human and armadillo-elephant for the second 3-way alignment) BLASTZ alignments (Chiaromonte et al. 2002; Kent et al. 2003;

Schwartz et al. 2003) and the standard matrix used in the UCSC browser. The alignments were fed into axtChain with the parameters -chainMinScore=3000 -chainLinearGap=medium, which organizes all alignments by chromosome (or scaffold) and creates a kd-tree from the gapless subsections (blocks) of the alignments. A dynamic program was then run over the kd- trees to find the maximally scoring chains of these blocks. Finally MULTIZ (Blanchette et al.

2004) was used to generate the two 3-way alignments referencing elephant and armadillo.

We used the first alignment to search for elements supporting the Atlantogenata

[(armadillo + elephant) - human] and Exafroplacentalia [(armadillo + human) - elephant] hypotheses, and the second one to test the Epitheria [(elephant + human) - armadillo] hypothesis.

Presence/absence loci

The MAF alignments were organized in blocks of conservation including chromosomal coordinates (Figure 2A). Using our own C-language script we extracted a total of 653,181 16 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

triple blocks, each composed of three continuous conserved regions the middle one of which included sequence regions larger than 50 nts in two of the species that was missing in the third species (corresponding to the potential absence of the corresponding sequence in this species). All triple blocks were converted into FASTA format using another C-language script. In addition, using more relaxed conditions we screened all blocks (1,636,542) for orthologous retroposon regions in two species that corresponded to at least 80% of the gaps in the third species so as to also consider potentially informative but slightly misaligned regions.

Computational screening for potential informative retroposons

With the local RepeatMasker (www.repeatmasker.org) we screened the triple block sequences for mammalian-specific L1MB LINEs. We selected triple blocks in which at least 50% of block 2 was composed of an L1MB element (human + elephant 284 loci; human + armadillo

256 loci; elephant + armadillo 677 loci). In addition, to detect slightly misaligned but informative loci, all conserved blocks were independently screened for retroelements and the corresponding regions compared among all three species for potential absences (10 cases).

Visual screening for informative retroposons

Links to the Genome Browser (http://genome.ucsc.edu/cgi-bin/hgBlat) were generated to filter out loci with less than 70% sequence similarity between multiple species and retroposons that overlapped with blocks 1 and/or 3, leaving us with a total of 22 phylogenetically informative loci.

Experimental amplifiction to detect presence/absence of retroposons in sloth and tenrec.

As sufficient sequence data were not yet available for all the phylogenetically informative loci in sloth and tenrec, we PCR-amplified some of those missing loci to provide a more complete presence/absence pattern. PCR reactions were performed using Phusion DNA Polymerase

17 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

(New England BioLabs, Beverly, MA) for 30 s at 98°C followed by 35 cycles of 10 s at 98°C,

30 s at primer specific temperatures (Supplementary Table S3), and 30 s at 72°C. Amplified

PCR fragments were sequenced directly or purified on agarose gels, ligated into the pDrive

Cloning Vector (Qiagen, Hilden) and electroporated into TOP10 cells (Invitrogen,

Groningen). Sequencing was performed using the Ampli Taq FS Big Dye Terminator Kit (PE

Biosystems, Foster City). A list of the PCR primers used is shown in Supplementary Table

S3.

Comparative analyses

For each locus identified in the MAF alignments, we compiled available sequence information of other selected mammalian species from trace or genomic data sequences from mouse, guinea pig, rabbit, dog, horse, cow, tenrec, hyrax, sloth, and as far as available, from opossum. All inserted L1MB elements were aligned against the consensus repbase sequences

(www.girinst.org/repbase) (Kapitonov and Jurka 2008). Each insertion site was carefully checked to validate the orthology of the insertions based on the following criteria: (1) identical genomic insertion points, (2) element orientation, (3) identical LINE1 subtypes, (4) matching length and sequence of direct repeats, and (5) identical truncation points. We tried as much as was possible to compile a representative and consistent species sampling for each marker (Supplemental Table S1; Supplemental Material 1). This was not always possible because the genomic information from some species remains fragmented. For consistent comparisons for most markers we selected two Supraprimates from human, mouse, guinea pig, and rabbit; two Laurasiatheria from dog, horse, and cow; two Afrotheria from elephant, tenrec, and hyrax; two Xenarthra (armadillo and sloth); and if alignable, the opossum outgroup sequence.

18 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Acknowledgments

We thank Marsha Bundman for editorial assistance and Ronald Blakey for providing the two terrestrial globes in Figure 4. Many thanks go to Webb Miller and Maria Nilsson for providing valuable comments on the manuscript. We thank Agencourt Biosciences, the

Baylor Genome Sequencing Center, and the Washington University Genome Center for providing raw genomic sequences (trace archives). This work was supported by the Deutsche

Forschungsgemeinschaft (SCHM1469/3-1).

19 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Figure legends

Figure 1. Different hypotheses of the placental origin. (A) The Exafroplacentalia hypothesis proposes a mammalian clade merging Boreotheria (Supraprimates plus

Laurasiathera) and Xenarthra, and Afrotheria at the base. (B) The Epitheria merges

Boreotheria and Afrotheria, with Xenarthra at the base. (C) The Atlantogenata hypothesis proposes Xenarthra and Afrotheria in one clade as the base.

Figure 2. Computational strategy to extract phylogenetically informative retroposon presence/absence loci. (A) UCSC 3-way Multiple Alignment Format (MAF) data are organized in conserved alignment blocks with chromosomal coordinates. Three adjacent blocks including an internal gap region in block 2 were transformed in FASTA format.

Sequences above gap regions were analyzed for mammalian specific repeats using the

RepeatMasker. Loci including such repeats were inspected in the UCSC Genome Browser.

All available sequences of phylogenetically informative loci were transformed to FASTA sequences and manually re-aligned. (B) 5, 9, and 8 loci were identified that support the three competing hypotheses of the origin of placental mammals, Exafroplacentalia, Epitheria and

Atlantogenata, respectively. (+) denotes the presence and (-) the absence of a diagnostic retroposon.

Figure 3. Sequence alignments of representative informative loci. Two representative, multispecies alignments of phylogenetically informative L1MB insertions supporting the

Exafroplacentalia (markers 1b and 1d), Epitheria (markers 2a and 2c), and Atlantogenata

(markers 3b and 3c) hypotheses. For full species sampling see Supplemental Material 1.

Potential direct repeats are boxed. The sequences above each set of alignments represent the

RepeatMasker consensus elements; w corresponds to a or t. Subscript numbers indicate serial

20 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

identical nucleotides. Note that the 5’ ends of the L1MB elements are truncated. The 5’ and 3’ ends of the retroposed insertions are partially shown in lower case letters.

Figure 4. Hypothetical scenarios of interspersed ancestry of placentals. (A)

Hybridization/introgression scenario: Initially the three placental lineages diverged and separated from an Eomaia-like (Eomaia scansoria) common ancestor. Alternate gene flow between the lineages led to a mosaic pattern of relationships among the current species. The sequence of genetic flow leading to significant support for Exafroplacentalia, Epitheria and

Atlantogenata was chosen randomly and not meant to indicate a specific time course. The lineages leading to Boreotheria, Afrotheria, and Xenarthra are shown in red, green, and blue, respectively. Gray circles represent diagnostic retroposon insertions for Exafroplacentalia, orange circles for Epitheria, and black circles for Atlantogenata. (B) Geographic transitions at the time of separation of the three lineages during the Cretaceous. Left terrestrial globe:

Possible overlapping geographical distribution of the three placental lineages before continental drift, the expected period of genetic flow between the lineages. Right terrestrial globe: separation of the continents and divergence of placental lineages possibly prohibiting further gene flow (Murphy et al. 2001b; Wildman et al. 2007). Roman numerals indicate exchanges of genetic material between the species shown in (A). (C) Randomly sorted ancestral alleles: Fixation of retroposon alleles was not completed in the ancestral population before the separation of lineages. After the proliferation of inserted retroposon alleles, the distributions to different lineages occurred randomly (dotted lines). After fixation of the alleles in the terminal lineages, a mosaic pattern of relatedness remains. (I) Five orthologous markers represent Exafroplacentalia, (II) 9 markers (plus an additional marker previously described in Kriegs et al. (2006), that was undetected here) Epitheria, and (III) 8

Atlantogenata. Unfilled circles represent alleles without retroelements.

21 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

References

Amrine-Madsen, H., Koepfli, K.P., Wayne, R.K., and Springer, M.S. 2003. A new phylogenetic

marker, apolipoprotein B, provides compelling evidence for eutherian relationships. Mol.

Phylogenet. Evol. 28: 225-240.

Blanchette, M., Kent, W.J., Riemer, C., Elnitski, L., Smit, A.F.A., Roskin, K.M., Baertsch, R.,

Rosenbloom, K., Clawson, H., Green, E.D., Haussler, D., and Miller, W. 2004. Aligning

muliple genomic sequences with the threaded blockset aligner. Genome Res. 14:708-715.

Cantrell, M.A., Filanoski, B.J., Ingermann, A.R., Olsson, K., DiLuglio, N., Lister, Z., and

Wichman, H.A. 2001. An ancient retrovirus-like element contains hot spots for SINE

insertion. Genetics 158: 769-777.

Chiaromonte, F., Yap, V.B., and Miller, W. 2002. Scoring pairwise genomic sequence

alignments. Pac Symp Biocomput 2002: 115-126.

Delsuc, F., Scally, M., Madsen, O., Stanhope, M.J., de Jong, W.W., Catzeflis, F.M., Springer,

M.S., and Douzery, E.J. 2002. Molecular phylogeny of living xenarthrans and the impact

of character and taxon sampling on the placental tree rooting. Mol. Biol. Evol. 19: 1656-

1671.

Ebersberger, I., Galgoczy, P., Taudien, S., Taenzer, S., Platzer, M., and von Haeseler, A. 2007.

Mapping human genetic ancestry. Mol. Biol. Evol. 24: 2266-2276.

Hallström, B.M. Kullberg M., Nilsson M.A., and Janke A. 2007. Phylogenomic data analyses

provide evidence that Xenarthra and Afrotheria are sister groups. Mol. Biol. Evol. 24:

2059-2068.

Hallström, B.M. and Janke, A. 2008. Resolution among major placental mammal interordinal

relationships with genome data imply that speciation influenced their earliest radiations.

BMC Evol. Biol. 8: 162.

22 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Huchon, D., Chevret, P., Jordan, U., Kilpatrick, C.W., Ranwez, V., Jenkins, P.D., Brosius, J., and

Schmitz, J. 2007. Multiple molecular evidences for a living mammalian fossil. Proc. Natl.

Acad. Sci. 104: 7495-7499.

Ji, Q., Luo, Z.X., Yuan, C.X., Wible, J.R., Zhang, J.P., and Georgi, J.A. 2002. The earliest known

eutherian mammal. Nature 416: 816-822.

Kapitonov, V.V. and Jurka, J. 2008. A universal classification of eukaryotic transposable

elements implemented in Repbase. Nat. Rev. Genet. 9: 411-412.

Kent, W.J., Baertsch, R., Hinrichs, A., Miller, W., and Haussler, D. 2003. Evolution's cauldron:

duplication, deletion, and rearrangement in the mouse and human genomes. Proc. Natl.

Acad. Sci. 100: 11484-11489.

Kimura, M. 1968. Evolutionary rate at the molecular level. Nature 217: 624-626.

Kjer, K.M. and Honeycutt, R.L. 2007. Site specific rates of mitochondrial genomes and the

phylogeny of eutheria. BMC Evol. Biol. 7: 8.

Kriegs, J.O., Churakov, G., Jurka, J., Brosius, J., and Schmitz, J. 2007. Evolutionary history of

7SL RNA-derived SINEs in Supraprimates. Trends Genet. 23: 158-161.

Kriegs, J.O., Churakov, G., Kiefmann, M., Jordan, U., Brosius, J., and Schmitz, J. 2006.

Retroposed elements as archives for the evolutionary history of placental mammals. PLoS

Biol. 4: e91.

Madsen, O., Scally, M., Douady, C.J., Kao, D.J., DeBry, R.W., Adkins, R., Amrine, H.M.,

Stanhope, M.J., de Jong, W.W., and Springer, M.S. 2001. Parallel adaptive radiations in

two major clades of placental mammals. Nature 409: 610-614.

Möller-Krull, M., Delsuc, F., Churakov, G., Marker, C., Superina, M., Brosius, J., Douzery, E.J.,

and Schmitz, J. 2007. Retroposed elements and their flanking regions resolve the

23 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

evolutionary history of xenarthran mammals (armadillos, anteaters, and sloths). Mol. Biol.

Evol. 24: 2573-2582.

Murphy, W.J., Eizirik, E., Johnson, W.E., Zhang, Y.P., Ryder, O.A., and O'Brien, S.J. 2001a.

Molecular phylogenetics and the origins of placental mammals. Nature 409: 614-618.

Murphy, W.J., Eizirik, E., O'Brien, S.J., Madsen, O., Scally, M., Douady, C.J., Teeling, E.,

Ryder, O.A., Stanhope, M.J., de Jong, W.W. et al. 2001b. Resolution of the early

placental mammal radiation using Bayesian phylogenetics. Science 294: 2348-2351.

Murphy, W.J., Pringle, T.H., Crider, T.A., Springer, M.S., and Miller, W. 2007. Using genomic

data to unravel the root of the placental mammal phylogeny. Genome Res. 17: 413-421.

Nikolaev, S., Montoya-Burgos, J.I., Margulies, E.H., Rougemont, J., Nyffeler, B., and

Antonarakis, S.E. 2007. Early history of mammals is elucidated with the ENCODE

multiple species sequencing data. PLoS Genet. 3: e2.

Nishihara, H., Satta, Y., Nikaido, M., Thewissen, J.G.M., Stanhope, M.J., and Okada, N. 2005. A

retroposon analysis of afrotherian phylogeny. Mol. Biol. Evol. 22: 1823-1833.

Nishihara, H., Hasegawa, M., and Okada, N. 2006. Pegasoferae, an unexpected mammalian clade

revealed by tracking ancient retroposon insertions. Proc. Natl. Acad. Sci. 103: 9929-9934.

Nishihara, H., Okada, N., and Hasegawa, M. 2007. Rooting the eutherian tree: the power and

pitfalls of phylogenomics. Genome Biol. 8: R199.

Prasad, A.B., Allard, M.W., and Green, E.D. 2008. Confirming the phylogeny of mammals by

use of large comparative sequence data sets. Mol. Biol. Evol. 25: 1795-1808.

Ray, D.A., Xing, J., Salem, A.H., and Batzer, M.A. 2006. SINEs of a nearly perfect character.

Syst. Biol. 55: 928-935.

24 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Schmitz, J. and Zischler, H. 2002. Molecular cladistic markers and the infraordinal phylogenetic

relationships of primates. In Anthropoid Origins: New Visions (eds. R.F. Kay and C.

Ross). Kluwer Academic Press, NY.

Schwartz, S., Kent, W.J., Smit, A., Zhang, Z., Baertsch, R., Hardison, R.C., Haussler, D., and

Miller, W. 2003. Human-Mouse Alignments with BLASTZ. Genome Res. 13: 103-107.

Shedlock, A.M., Takahashi, K., and Okada, N. 2004. SINEs of speciation: tracking lineages with

retroposons. Trends Ecol. Evol. 19: 545-553.

Shoshani, J. and McKenna, M.C. 1998. Higher taxonomic relationships among extant mammals

based on morphology, with selected comparisons of results from molecular data. Mol.

Phylogenet Evol. 9: 572-584.

Steel, M. and Penny, D. 2000. Parsimony, likelihood, and the role of models in molecular

phylogenetics. Mol. Biol. Evol. 17: 839-850. van de Lagemaat, L.N., Gagnier, L., Medstrand, P., and Mager, D.L. 2005. Genomic deletions

and precise removal of transposable elements mediated by short identical DNA segments

in primates. Genome Res. 15: 1243-1249.

Waddell, P.J., Cao, Y., Hasegawa, M., and Mindell, D.P. 1999. Assessing the Cretaceous

superordinal divergence times within birds and placental mammals by using whole

mitochondrial protein sequences and an extended statistical framework. Syst. Biol. 48:

119-137.

Waddell, P.J., Kishino, H., and Ota, R. 2001. A phylogenetic foundation for comparative

mammalian genomics. Genome Inform. 12: 141-154.

Waddell, P.J. and Shelley, S. 2003. Evaluating placental inter-ordinal phylogenies with novel

sequences including RAG1, gamma-fibrinogen, ND6, and mt-tRNA, plus MCMC-driven

nucleotide, amino acid, and codon models. Mol. Phylogenet Evol. 28: 197-224.

25 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Warren, W.C. Hillier, L.W. Marshall Graves, J.A. Birney, E. Ponting, C.P. Grutzner, F. Belov, K.

Miller, W. Clarke, L. Chinwalla, A.T. et al. 2008. Genome analysis of the platypus

reveals unique signatures of evolution. Nature 453: 175-183.

Wildman, D.E., Uddin, M., Opazo, J.C., Liu, G., Lefort, V., Guindon, S., Gascuel, O., Grossman,

L.I., Romero, R., and Goodman, M. 2007. Genomics, biogeography, and the

diversification of placental mammals. Proc. Natl. Acad. Sci. 104: 14395-14400.

Xing, J., Witherspoon, D.J., Ray, D.A., Batzer, M.A., and Jorde, L.B. 2007. Mobile DNA

elements in primate and human evolution. Am. J. Phys. Anthropol Suppl 45: 2-19.

26 Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press Downloaded from genome.cshlp.org on September 30, 2021 - Published by Cold Spring Harbor Laboratory Press

Mosaic Retroposon Insertion Patterns in Placental Mammals

Gennady Churakov, Jan Ole Kriegs, Robert Baertsch, et al.

Genome Res. published online March 4, 2009 Access the most recent version at doi:10.1101/gr.090647.108

Supplemental http://genome.cshlp.org/content/suppl/2009/05/01/gr.090647.108.DC1 Material

P

Accepted Peer-reviewed and accepted for publication but not copyedited or typeset; accepted Manuscript manuscript is likely to differ from the final, published version.

License

Email Alerting Receive free email alerts when new articles cite this article - sign up in the box at the Service top right corner of the article or click here.

Advance online articles have been peer reviewed and accepted for publication but have not yet appeared in the paper journal (edited, typeset versions may be posted when available prior to final publication). Advance online articles are citable and establish publication priority; they are indexed by PubMed from initial publication. Citations to Advance online articles must include the digital object identifier (DOIs) and date of initial publication.

To subscribe to Genome Research go to: https://genome.cshlp.org/subscriptions

Copyright © 2009, Cold Spring Harbor Laboratory Press