The Tasmanian tiger or thylacine (Thylacinus cynocephalus) was the largest carnivorous Australian marsupial to survive into the modern era. Despite last sharing a common ancestor with the eutherian canids ~160 million years ago, their phenotypic resemblance is considered the most striking example of convergent evolution in mammals. The last known thylacine died in captivity in 1936 and many aspects of the evolutionary history of this unique marsupial apex predator remain unknown. Here we have sequenced the genome of a preserved thylacine pouch young specimen to clarify the phylogenetic position of the thylacine within the carnivorous marsupials, reconstruct its historical demography and examine the genetic basis of its convergence with canids. Retroposon insertion patterns placed the thylacine as the basal lineage in Dasyuromorphia and suggest incomplete lineage sorting in early dasyuromorphs. Demographic analysis indicated a long-term decline in genetic diversity starting well before the arrival of humans in Australia. In spite of their extraordinary phenotypic convergence, comparative genomic analyses demonstrated that amino acid homoplasies between the thylacine and canids are largely consistent with neutral evolution. Furthermore, the genes and pathways targeted by positive selection differ markedly between these species. Together, these findings support models of adaptive convergence driven primarily bycis-regulatory evolution.

Lateral flight membranes, or patagia, have evolved repeatedly in diverse mammalian lineages. While little is known about patagium development, its recurrent evolution may suggest a shared molecular basis. By combining transcriptomics, developmental experiments, and mouse transgenics, we demonstrate that lateral Wnt5a expression in the marsupial sugar glider ( Petaurus breviceps ) promotes the differentiation of its patagium primordium. We further show that this function of Wnt5a reprises ancestral roles in skin morphogenesis predating mammalian flight and has been convergently used during patagium evolution in eutherian bats. Moreover, we find that many genes involved in limb development have been redeployed during patagium outgrowth in both the sugar glider and bat. Together, our findings reveal that deeply conserved genetic toolkits contribute to the evolutionary transition to flight in mammals.

The extinct marsupial Tasmanian tiger, or thylacine, and the eutherian gray wolf are among the most widely recognized examples of convergent evolution in mammals. Despite being distantly related, these large predators independently evolved extremely similar craniofacial morphologies, and evidence suggests that they filled similar ecological niches. Previous analyses revealed little evidence of adaptive convergence between their protein-coding genes. Thus, the genetic basis of their convergence is still unclear. Here, we identified candidate craniofacial cis-regulatory elements across vertebrates and compared their evolutionary rates in the thylacine and wolf, revealing abundant signatures of convergent positive selection. Craniofacial thylacine–wolf accelerated regions were enriched near genes involved in TGF beta (TGFB) and BMP signaling, both of which are key morphological signaling pathways with critical roles in establishing the identities and boundaries between craniofacial tissues. Similarly, enhancers of genes involved in craniofacial nerve development showed convergent selection and involvement in these pathways. Taken together, these results suggest that adaptation in cis-regulators of TGF beta and BMP signaling may provide a mechanism to explain the coevolution of developmentally and functionally integrated craniofacial structures in these species. We also found that despite major structural differences in marsupial and eutherian brains, accelerated regions in both species were common near genes with roles in brain development. Our findings support the hypothesis that, relative to protein-coding genes, positive selection on cis-regulatory elements is likely to be an essential driver of adaptive convergent evolution and may underpin thylacine–wolf phenotypic similarities.

Phenotypic variation among species is a product of evolutionary changes to developmental programs^1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium—an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

The extinct Tasmanian tiger or thylacine (Thylacinus cynocephalus) was a large marsupial carnivore native to Australia. Once ranging across parts of the mainland, the species remained only on the island of Tasmania by the time of European colonization. It was driven to extinction in the early 20th century and is an emblem of native species loss in Australia. The thylacine was a striking example of convergent evolution with placental canids, with which it shared a similar skull morphology. Consequently, it has been the subject of extensive study. While the original thylacine assemblies published in 2018 enabled the first exploration of the species’ genome biology, further progress is hindered by the lack of high-quality genomic resources. Here, we present a new chromosome-scale hybrid genome assembly for the thylacine, which compares favorably with many recent de novo marsupial genomes. In addition, we provide homology-based gene annotations, characterize the repeat content of the thylacine genome, and show that consistent with demographic decline, the species possessed a low rate of heterozygosity even compared to extant, threatened marsupials.