Proteo-Transcriptomic Analysis Identifies Potential Novel Toxins Secreted by the Predatory, Prey-Piercing Ribbon Worm Amphiporus lactifloreus.

Nemerteans (ribbon worms) employ toxins to subdue their prey, but research thus far has focused on the small-molecule components of mucus secretions and few protein toxins have been characterized. We carried out a preliminary proteotranscriptomic analysis of putative toxins produced by the hoplonemertean Amphiporus lactifloreus (Hoplonemertea, Amphiporidae). No variants were found of known nemertean-specific toxin proteins (neurotoxins, cytotoxins, parbolysins or nemertides) but several toxin-like transcripts were discovered, expressed strongly in the proboscis, including putative metalloproteinases and sequences resembling sea anemone actitoxins, crown-of-thorn sea star plancitoxins, and multiple classes of inhibitor cystine knot/knottin family proteins. Some of these products were also directly identified in the mucus proteome, supporting their preliminary identification as secreted toxin components. Two new nemertean-typical toxin candidates could be described and were named U-nemertotoxin-1 and U-nemertotoxin-2. Our findings provide insight into the largely overlooked venom system of nemerteans and support a hypothesis in which the nemertean proboscis evolved in several steps from a flesh-melting organ in scavenging nemerteans to a flesh-melting and toxin-secreting venom apparatus in hunting hoplonemerteans.

Composition of marsupial zona pellucida: a molecular and phylogenetic approach.

The zona pellucida (ZP) is an extracellular matrix that surrounds mammalian oocytes. In eutherians it is formed from three or four proteins (ZP1, ZP2, ZP3, ZP4). In the few marsupials that have been studied, however, only three of these have been characterised (ZP2, ZP3, ZP4). Nevertheless, the composition in marsupials may be more complex, since a duplication of the ZP3 gene was recently described in one species. The aim of this work was to elucidate the ZP composition in marsupials and relate it to the evolution of the ZP gene family. For that, an in silico and molecular analysis was undertaken, focusing on two South American species (gray short-tailed opossum and common opossum) and five Australian species (brushtail possum, koala, Bennett's wallaby, Tammar wallaby and Tasmanian devil). This analysis identified the presence of ZP1 mRNA and mRNA from two or three paralogues of ZP3 in marsupials. Furthermore, evidence for ZP1 and ZP4 pseudogenes in the South American subfamily Didelphinae and for ZP3 pseudogenes in two marsupials is provided. In conclusion, two different composition models are proposed for marsupials: a model with four proteins (ZP1, ZP2 and ZP3 (two copies)) for the South American species and a model with six proteins (ZP1, ZP2, ZP3 (three copies) and ZP4) for the Australasian species.

Retroviral envelope gene captures and syncytin exaptation for placentation in marsupials.

Syncytins are genes of retroviral origin captured by eutherian mammals, with a role in placentation. Here we show that some marsupials-which are the closest living relatives to eutherian mammals, although they diverged from the latter ∼190 Mya-also possess a syncytin gene. The gene identified in the South American marsupial opossum and dubbed syncytin-Opo1 has all of the characteristic features of a bona fide syncytin gene: It is fusogenic in an ex vivo cell-cell fusion assay; it is specifically expressed in the short-lived placenta at the level of the syncytial feto-maternal interface; and it is conserved in a functional state in a series of Monodelphis species. We further identify a nonfusogenic retroviral envelope gene that has been conserved for >80 My of evolution among all marsupials (including the opossum and the Australian tammar wallaby), with evidence for purifying selection and conservation of a canonical immunosuppressive domain, but with only limited expression in the placenta. This unusual captured gene, together with a third class of envelope genes from recently endogenized retroviruses-displaying strong expression in the uterine glands where retroviral particles can be detected-plausibly correspond to the different evolutionary statuses of a captured retroviral envelope gene, with only syncytin-Opo1 being the present-day bona fide syncytin active in the opossum and related species. This study would accordingly recapitulate the natural history of syncytin exaptation and evolution in a single species, and definitely extends the presence of such genes to all major placental mammalian clades.

Evolutionary histories of transposable elements in the genome of the largest living marsupial carnivore, the Tasmanian devil.

The largest living carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii), is the sole survivor of a lineage originating about 12 Ma. We set out to investigate the spectrum of transposable elements found in the Tasmanian devil genome, the first high-coverage genome of an Australian marsupial. Marsupial genomes have been shown to have the highest amount of transposable elements among vertebrates. We analyzed the horizontally transmitted DNA transposons OC1 and hAT-1_MEu in the Tasmanian devil genome. OC1 is present in all carnivorous marsupials, while having a very limited distribution among the remaining Australian marsupial orders. In contrast, hAT-1_MEu is present in all Australian marsupial orders, and has so far only been identified in a few placental mammals. We screened 158 introns for phylogenetically informative retrotransposons in the order Dasyuromorphia, and found that the youngest SINE (Short INterspersed Element), WSINE1, is no longer active in the subfamily Dasyuridae. The lack of detectable WSINE1 activity in this group may be due to a retrotransposon inactivation event approximately 30 Ma. We found that the Tasmanian devil genome contains a relatively low number of continuous full-length LINE-1 (Long INterspersed Element 1, L1) retrotransposons compared with the opossum genome. Furthermore, all L1 elements in the Tasmanian devil appeared to be nonfunctional. Hidden Markov Model approaches suggested that other potential sources of functional reverse transcriptase are absent from the genome. We discuss the issues associated with assembling long, highly similar L1 copies from short read Illumina data and describe how assembly artifacts can potentially lead to erroneous conclusions.

Genetic signatures of adaptation revealed from transcriptome sequencing of Arctic and red foxes.

BACKGROUND: The genus Vulpes (true foxes) comprises numerous species that inhabit a wide range of habitats and climatic conditions, including one species, the Arctic fox (Vulpes lagopus) which is adapted to the arctic region. A close relative to the Arctic fox, the red fox (Vulpes vulpes), occurs in subarctic to subtropical habitats. To study the genetic basis of their adaptations to different environments, transcriptome sequences from two Arctic foxes and one red fox individual were generated and analyzed for signatures of positive selection. In addition, the data allowed for a phylogenetic analysis and divergence time estimate between the two fox species. RESULTS: The de novo assembly of reads resulted in more than 160,000 contigs/transcripts per individual. Approximately 17,000 homologous genes were identified using human and the non-redundant databases. Positive selection analyses revealed several genes involved in various metabolic and molecular processes such as energy metabolism, cardiac gene regulation, apoptosis and blood coagulation to be under positive selection in foxes. Branch site tests identified four genes to be under positive selection in the Arctic fox transcriptome, two of which are fat metabolism genes. In the red fox transcriptome eight genes are under positive selection, including molecular process genes, notably genes involved in ATP metabolism. Analysis of the three transcriptomes and five Sanger re-sequenced genes in additional individuals identified a lower genetic variability within Arctic foxes compared to red foxes, which is consistent with distribution range differences and demographic responses to past climatic fluctuations. A phylogenomic analysis estimated that the Arctic and red fox lineages diverged about three million years ago. CONCLUSIONS: Transcriptome data are an economic way to generate genomic resources for evolutionary studies. Despite not representing an entire genome, this transcriptome analysis identified numerous genes that are relevant to arctic adaptation in foxes. Similar to polar bears, fat metabolism seems to play a central role in adaptation of Arctic foxes to the cold climate, as has been identified in the polar bear, another arctic specialist.

Genomic resources and genetic diversity of captive lesser kudu (Tragelaphus imberbis).

The lesser kudu (Tragelaphus imberbis) is a spiral-horned antelope native to northeastern Africa. Individuals kept in zoological gardens are suspected to be highly inbred due to few founder individuals and a small breeding stock. A morphological study suggested two distinct subspecies of the lesser kudu. However, subspecies designation and population structure in zoological gardens has not been analyzed using molecular markers. We analyzed one mitochondrial marker and two nuclear intron loci (total: 2,239 nucleotides) in 52 lesser kudu individuals. Of these, 48 individuals were bred in captivity and sampled from seven different zoos. The four remaining individuals were recently captured in Somalia and are currently held in the Maktoum zoo. Maternally inherited mitochondrial sequences indicate substantial amounts of genetic variation in the zoo populations, while the biparentally inherited intron sequences are, as expected, less variable. The analyzed individuals show 10 mitochondrial haplotypes with a maximal distance of 10 mutational steps. No prominent subspecies structure is detectable in this study. For further studies of the lesser kudu population genetics, we present microsatellite markers from a low-coverage genome survey using 454 sequencing technology.

Near chromosome-level and highly repetitive genome assembly of the snake pipefish Entelurus aequoreus (Syngnathiformes: Syngnathidae).

The snake pipefish, Entelurus aequoreus (Linnaeus, 1758), is a northern Atlantic fish inhabiting open seagrass environments that recently expanded its distribution range. Here, we present a highly contiguous, near chromosome-scale genome of E. aequoreus. The final assembly spans 1.6 Gbp in 7,391 scaffolds, with a scaffold N50 of 62.3 Mbp and L50 of 12. The 28 largest scaffolds (>21 Mbp) span 89.7% of the assembly length. A BUSCO completeness score of 94.1% and a mapping rate above 98% suggest a high assembly completeness. Repetitive elements cover 74.93% of the genome, one of the highest proportions identified in vertebrates. Our demographic modeling identified a peak in population size during the last interglacial period, suggesting the species might benefit from warmer water conditions. Our updated snake pipefish assembly is essential for future analyses of the morphological and molecular changes unique to the Syngnathidae.

Tracking marsupial evolution using archaic genomic retroposon insertions.

The Australasian and South American marsupial mammals, such as kangaroos and opossums, are the closest living relatives to placental mammals, having shared a common ancestor around 130 million years ago. The evolutionary relationships among the seven marsupial orders have, however, so far eluded resolution. In particular, the relationships between the four Australasian and three South American marsupial orders have been intensively debated since the South American order Microbiotheria was taxonomically moved into the group Australidelphia. Australidelphia is significantly supported by both molecular and morphological data and comprises the four Australasian marsupial orders and the South American order Microbiotheria, indicating a complex, ancient, biogeographic history of marsupials. However, the exact phylogenetic position of Microbiotheria within Australidelphia has yet to be resolved using either sequence or morphological data analysis. Here, we provide evidence from newly established and virtually homoplasy-free retroposon insertion markers for the basal relationships among marsupial orders. Fifty-three phylogenetically informative markers were retrieved after in silico and experimental screening of approximately 217,000 retroposon-containing loci from opossum and kangaroo. The four Australasian orders share a single origin with Microbiotheria as their closest sister group, supporting a clear divergence between South American and Australasian marsupials. In addition, the new data place the South American opossums (Didelphimorphia) as the first branch of the marsupial tree. The exhaustive computational and experimental evidence provides important insight into the evolution of retroposable elements in the marsupial genome. Placing the retroposon insertion pattern in a paleobiogeographic context indicates a single marsupial migration from South America to Australia. The now firmly established phylogeny can be used to determine the direction of genomic changes and morphological transitions within marsupials.

Expansion of CORE-SINEs in the genome of the Tasmanian devil.

BACKGROUND: The genome of the carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii, Order: Dasyuromorphia), was sequenced in the hopes of finding a cure for or gaining a better understanding of the contagious devil facial tumor disease that is threatening the species' survival. To better understand the Tasmanian devil genome, we screened it for transposable elements and investigated the dynamics of short interspersed element (SINE) retroposons. RESULTS: The temporal history of Tasmanian devil SINEs, elucidated using a transposition in transposition analysis, indicates that WSINE1, a CORE-SINE present in around 200,000 copies, is the most recently active element. Moreover, we discovered a new subtype of WSINE1 (WSINE1b) that comprises at least 90% of all Tasmanian devil WSINE1s. The frequencies of WSINE1 subtypes differ in the genomes of two of the other Australian marsupial orders. A co-segregation analysis indicated that at least 66 subfamilies of WSINE1 evolved during the evolution of Dasyuromorphia. Using a substitution rate derived from WSINE1 insertions, the ages of the subfamilies were estimated and correlated with a newly established phylogeny of Dasyuromorphia. Phylogenetic analyses and divergence time estimates of mitochondrial genome data indicate a rapid radiation of the Tasmanian devil and the closest relative the quolls (Dasyurus) around 14 million years ago. CONCLUSIONS: The radiation and abundance of CORE-SINEs in marsupial genomes indicates that they may be a major player in the evolution of marsupials. It is evident that the early phases of evolution of the carnivorous marsupial order Dasyuromorphia was characterized by a burst of SINE activity. A correlation between a speciation event and a major burst of retroposon activity is for the first time shown in a marsupial genome.

Speciation and population divergence in a mutualistic seed dispersing bird.

Bird-mediated seed dispersal is crucial for the regeneration and viability of ecosystems, often resulting in complex mutualistic species networks. Yet, how this mutualism drives the evolution of seed dispersing birds is still poorly understood. In the present study we combine whole genome re-sequencing analyses and morphometric data to assess the evolutionary processes that shaped the diversification of the Eurasian nutcracker (Nucifraga), a seed disperser known for its mutualism with pines (Pinus). Our results show that the divergence and phylogeographic patterns of nutcrackers resemble those of other non-mutualistic passerine birds and suggest that their early diversification was shaped by similar biogeographic and climatic processes. The limited variation in foraging traits indicates that local adaptation to pines likely played a minor role. Our study shows that close mutualistic relationships between bird and plant species might not necessarily act as a primary driver of evolution and diversification in resource-specialized birds.

Genome Size Dynamics in Marine Ribbon Worms (Nemertea, Spiralia).

Nemertea is a phylum consisting of 1300 mostly marine species. Nemertea is distinguished by an eversible muscular proboscis, and most of the species are venomous. Genomic resources for this phylum are scarce despite their value in understanding biodiversity. Here, we present genome size estimates of Nemertea based on flow cytometry and their relationship to different morphological and developmental traits. Ancestral genome size estimations were done across the nemertean phylogeny. The results increase the available genome size estimates for Nemertea three-fold. Our analyses show that Nemertea has a narrow genome size range (0.43-3.89 pg) compared to other phyla in Lophotrochozoa. A relationship between genome size and evolutionary rate, developmental modes, and habitat was found. Trait analyses show that the highest evolutionary rate of genome size is found in upper intertidal, viviparous species with direct development. Despite previous findings, body size in nemerteans was not correlated with genome size. A relatively small genome (1.18 pg) is assumed for the most recent common ancestor of all extant nemerteans. The results provide an important basis for future studies in nemertean genomics, which will be instrumental to understanding the evolution of this enigmatic and often neglected phylum.

Population analysis of retrotransposons in giraffe genomes supports RTE decline and widespread LINE1 activity in Giraffidae.

BACKGROUND: The majority of structural variation in genomes is caused by insertions of transposable elements (TEs). In mammalian genomes, the main TE fraction is made up of autonomous and non-autonomous non-LTR retrotransposons commonly known as LINEs and SINEs (Long and Short Interspersed Nuclear Elements). Here we present one of the first population-level analysis of TE insertions in a non-model organism, the giraffe. Giraffes are ruminant artiodactyls, one of the few mammalian groups with genomes that are colonized by putatively active LINEs of two different clades of non-LTR retrotransposons, namely the LINE1 and RTE/BovB LINEs as well as their associated SINEs. We analyzed TE insertions of both types, and their associated SINEs in three giraffe genome assemblies, as well as across a population level sampling of 48 individuals covering all extant giraffe species. RESULTS: The comparative genome screen identified 139,525 recent LINE1 and RTE insertions in the sampled giraffe population. The analysis revealed a drastically reduced RTE activity in giraffes, whereas LINE1 is still actively propagating in the genomes of extant (sub)-species. In concert with the extremely low activity of the giraffe RTE, we also found that RTE-dependent SINEs, namely Bov-tA and Bov-A2, have been virtually immobile in the last 2 million years. Despite the high current activity of the giraffe LINE1, we did not find evidence for the presence of currently active LINE1-dependent SINEs. TE insertion heterozygosity rates differ among the different (sub)-species, likely due to divergent population histories. CONCLUSIONS: The horizontally transferred RTE/BovB and its derived SINEs appear to be close to inactivation and subsequent extinction in the genomes of extant giraffe species. This is the first time that the decline of a TE family has been meticulously analyzed from a population genetics perspective. Our study shows how detailed information about past and present TE activity can be obtained by analyzing large-scale population-level genomic data sets.

Transposable Elements in the Genome of the Lichen-Forming Fungus Umbilicaria pustulata and Their Distribution in Different Climate Zones along Elevation.

Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Drift and natural selection are important forces shaping TE distribution and accumulation. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host's ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold temperate. We trace the occurrence of the newly identified TEs in populations along three elevation gradients using a Pool-Seq approach to identify TE insertions of potential adaptive significance. We found that TEs cover 21.26% of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. We identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function. This pioneering study of the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution.

Education in the genomics era: Generating high-quality genome assemblies in university courses.

Recent advances in genome sequencing technologies have simplified the generation of genome data and reduced the costs for genome assemblies, even for complex genomes like those of vertebrates. More practically oriented genomic courses can prepare university students for the increasing importance of genomic data used in biological and medical research. Low-cost third-generation sequencing technology, along with publicly available data, can be used to teach students how to process genomic data, assemble full chromosome-level genomes, and publish the results in peer-reviewed journals, or preprint servers. Here we outline experiences gained from 2 master's-level courses and discuss practical considerations for teaching hands-on genome assembly courses.

Chromosome-level genome assembly of a benthic associated Syngnathiformes species: the common dragonet, Callionymus lyra.

Background: The common dragonet, Callionymus lyra, is one of three Callionymus species inhabiting the North Sea. All three species show strong sexual dimorphism. The males show strong morphological differentiation, e.g., species-specific colouration and size relations, while the females of different species have few distinguishing characters. Callionymus belongs to the 'benthic associated clade' of the order Syngnathiformes. The 'benthic associated clade' so far is not represented by genome data and serves as an important outgroup to understand the morphological transformation in 'long-snouted' syngnatiformes such as seahorses and pipefishes. Findings: Here, we present the chromosome-level genome assembly of C. lyra. We applied Oxford Nanopore Technologies' long-read sequencing, short-read DNBseq, and proximity-ligation-based scaffolding to generate a high-quality genome assembly. The resulting assembly has a contig N50 of 2.2 Mbp and a scaffold N50 of 26.7 Mbp. The total assembly length is 568.7 Mbp, of which over 538 Mbp were scaffolded into 19 chromosome-length scaffolds. The identification of 94.5% complete BUSCO genes indicates high assembly completeness. Additionally, we sequenced and assembled a multi-tissue transcriptome with a total length of 255.5 Mbp that was used to aid the annotation of the genome assembly. The annotation resulted in 19,849 annotated transcripts and identified a repeat content of 27.7%. Conclusions: The chromosome-level assembly of C. lyra provides a high-quality reference genome for future population genomic, phylogenomic, and phylogeographic analyses.

Retrophylogenomics in rorquals indicate large ancestral population sizes and a rapid radiation.

BACKGROUND: Baleen whales (Mysticeti) are the largest animals on earth and their evolutionary history has been studied in detail, but some relationships still remain contentious. In particular, reconstructing the phylogenetic position of the gray whales (Eschrichtiidae) has been complicated by evolutionary processes such as gene flow and incomplete lineage sorting (ILS). Here, whole-genome sequencing data of the extant baleen whale radiation allowed us to identify transposable element (TE) insertions in order to perform phylogenomic analyses and measure germline insertion rates of TEs. Baleen whales exhibit the slowest nucleotide substitution rate among mammals, hence we additionally examined the evolutionary insertion rates of TE insertions across the genomes. RESULTS: In eleven whole-genome sequences representing the extant radiation of baleen whales, we identified 91,859 CHR-SINE insertions that were used to reconstruct the phylogeny with different approaches as well as perform evolutionary network analyses and a quantification of conflicting phylogenetic signals. Our results indicate that the radiation of rorquals and gray whales might not be bifurcating. The morphologically derived gray whales are placed inside the rorqual group, as the sister-species to humpback and fin whales. Detailed investigation of TE insertion rates confirm that a mutational slow down in the whale lineage is present but less pronounced for TEs than for nucleotide substitutions. CONCLUSIONS: Whole genome sequencing based detection of TE insertions showed that the speciation processes in baleen whales represent a rapid radiation. Large genome-scale TE data sets in addition allow to understand retrotransposition rates in non-model organisms and show the potential for TE calling methods to study the evolutionary history of species.

Speciation Generates Mosaic Genomes in Kangaroos.

The iconic Australasian kangaroos and wallabies represent a successful marsupial radiation. However, the evolutionary relationship within the two genera, Macropus and Wallabia, is controversial: mitochondrial and nuclear genes, and morphological data have produced conflicting scenarios regarding the phylogenetic relationships, which in turn impact the classification and taxonomy. We sequenced and analyzed the genomes of 11 kangaroos to investigate the evolutionary cause of the observed phylogenetic conflict. A multilocus coalescent analysis using ∼14,900 genome fragments, each 10 kb long, significantly resolved the species relationships between and among the sister-genera Macropus and Wallabia. The phylogenomic approach reconstructed the swamp wallaby (Wallabia) as nested inside Macropus, making this genus paraphyletic. However, the phylogenomic analyses indicate multiple conflicting phylogenetic signals in the swamp wallaby genome. This is interpreted as at least one introgression event between the ancestor of the genus Wallabia and a now extinct ghost lineage outside the genus Macropus. Additional phylogenetic signals must therefore be caused by incomplete lineage sorting and/or introgression, but available statistical methods cannot convincingly disentangle the two processes. In addition, the relationships inside the Macropus subgenus M. (Notamacropus) represent a hard polytomy. Thus, the relationships between tammar, red-necked, agile, and parma wallabies remain unresolvable even with whole-genome data. Even if most methods resolve bifurcating trees from genomic data, hard polytomies, incomplete lineage sorting, and introgression complicate the interpretation of the phylogeny and thus taxonomy.

Whole-genome sequencing of the blue whale and other rorquals finds signatures for introgressive gene flow.

Reconstructing the evolution of baleen whales (Mysticeti) has been problematic because morphological and genetic analyses have produced different scenarios. This might be caused by genomic admixture that may have taken place among some rorquals. We present the genomes of six whales, including the blue whale (Balaenoptera musculus), to reconstruct a species tree of baleen whales and to identify phylogenetic conflicts. Evolutionary multilocus analyses of 34,192 genome fragments reveal a fast radiation of rorquals at 10.5 to 7.5 million years ago coinciding with oceanic circulation shifts. The evolutionarily enigmatic gray whale (Eschrichtius robustus) is placed among rorquals, and the blue whale genome shows a high degree of heterozygosity. The nearly equal frequency of conflicting gene trees suggests that speciation of rorqual evolution occurred under gene flow, which is best depicted by evolutionary networks. Especially in marine environments, sympatric speciation might be common; our results raise questions about how genetic divergence can be established.

Phylogenetic Conflict in Bears Identified by Automated Discovery of Transposable Element Insertions in Low-Coverage Genomes.

Phylogenetic reconstruction from transposable elements (TEs) offers an additional perspective to study evolutionary processes. However, detecting phylogenetically informative TE insertions requires tedious experimental work, limiting the power of phylogenetic inference. Here, we analyzed the genomes of seven bear species using high-throughput sequencing data to detect thousands of TE insertions. The newly developed pipeline for TE detection called TeddyPi (TE detection and discovery for Phylogenetic Inference) identified 150,513 high-quality TE insertions in the genomes of ursine and tremarctine bears. By integrating different TE insertion callers and using a stringent filtering approach, the TeddyPi pipeline produced highly reliable TE insertion calls, which were confirmed by extensive in vitro validation experiments. Analysis of single nucleotide substitutions in the flanking regions of the TEs shows that these substitutions correlate with the phylogenetic signal from the TE insertions. Our phylogenomic analyses show that TEs are a major driver of genomic variation in bears and enabled phylogenetic reconstruction of a well-resolved species tree, despite strong signals for incomplete lineage sorting and introgression. The analyses show that the Asiatic black, sun, and sloth bear form a monophyletic clade, in which phylogenetic incongruence originates from incomplete lineage sorting. TeddyPi is open source and can be adapted to various TE and structural variation callers. The pipeline makes it possible to confidently extract thousands of TE insertions even from low-coverage genomes (∼10×) of nonmodel organisms. This opens new possibilities for biologists to study phylogenies and evolutionary processes as well as rates and patterns of (retro-)transposition and structural variation.

Resolving kangaroo phylogeny and overcoming retrotransposon ascertainment bias.

Reconstructing phylogeny from retrotransposon insertions is often limited by access to only a single reference genome, whereby support for clades that do not include the reference taxon cannot be directly observed. Here we have developed a new statistical framework that accounts for this ascertainment bias, allowing us to employ phylogenetically powerful retrotransposon markers to explore the radiation of the largest living marsupials, the kangaroos and wallabies of the genera Macropus and Wallabia. An exhaustive in silico screening of the tammar wallaby (Macropus eugenii) reference genome followed by experimental screening revealed 29 phylogenetically informative retrotransposon markers belonging to a family of endogenous retroviruses. We identified robust support for the enigmatic swamp wallaby (Wallabia bicolor) falling within a paraphyletic genus, Macropus. Our statistical approach provides a means to test for incomplete lineage sorting and introgression/hybridization in the presence of the ascertainment bias. Using retrotransposons as "molecular fossils", we reveal one of the most complex patterns of hemiplasy yet identified, during the rapid diversification of kangaroos and wallabies. Ancestral state reconstruction incorporating the new retrotransposon phylogenetic information reveals multiple independent ecological shifts among kangaroos into more open habitats, coinciding with the Pliocene onset of increased aridification in Australia from ~3.6 million years ago.