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.

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.

The Tasmanian tiger or thylacine (Thylacinus cynocephalus) was an iconic Australian marsupial predator that was hunted to extinction in the early 1900s. Despite sharing striking similarities with canids, they failed to evolve many of the specialized anatomical features that characterize carnivorous placental mammals. These evolutionary limitations are thought to arise from functional constraints associated with the marsupial mode of reproduction, in which otherwise highly altricial young use their well-developed forelimbs to climb to the pouch and mouth to suckle. Here we present the first three-dimensional digital developmental series of the thylacine throughout its pouch life using X-ray computed tomography on all known ethanol-preserved specimens. Based on detailed skeletal measurements, we refine the species growth curve to improve age estimates for the individuals. Comparison of allometric growth trends in the appendicular skeleton (fore- and hindlimbs) with that of other placental and marsupial mammals revealed that despite their unique adult morphologies, thylacines retained a generalized early marsupial ontogeny. Our approach also revealed mislabelled specimens that possessed large epipubic bones (vestigial in thylacine) and differing vertebral numbers. All of our generated CT models are publicly available, preserving their developmental morphology and providing a novel digital resource for future studies of this unique marsupial.

Runt-related transcription factor 2 (RUNX2) is critical for the development of the vertebrate bony skeleton. Unlike other RUNX family members, RUNX2 possesses a variable poly-glutamine, poly-alanine (QA) repeat domain. Natural variation within this repeat is able to alter the transactivation potential of RUNX2, acting as an evolutionary ‘tuning knob’ suggested to influence mammalian skull shape. However, the broader role of the RUNX2 QA repeat throughout vertebrate evolution is unknown. In this perspective, we examine the role of the RUNX2 QA repeat during skeletal development and discuss how its emergence and expansion may have facilitated the evolution of morphological novelty in vertebrates.

The lateral plate mesoderm (LPM) is a transient tissue that produces a diverse range of differentiated structures, including the limbs. However, the molecular mechanisms that drive early LPM specification and development are poorly understood. In this study, we use single-cell transcriptomics to define the cell-fate decisions directing LPM specification, subdivision and early initiation of the forelimb mesenchyme in chicken embryos. We establish a transcriptional atlas and global cell-cell signalling interactions in progenitor, transitional and mature cell types throughout the developing forelimb field. During LPM subdivision, somatic and splanchnic LPM fate is achieved through activation of lineage-specific gene modules. During the earliest stages of limb initiation, we identify activation of TWIST1 in the somatic LPM as a putative driver of limb bud epithelial-to-mesenchymal transition. Furthermore, we define a new role for BMP signalling during early limb development, revealing that it is necessary for inducing a somatic LPM fate and initiation of limb outgrowth, potentially through activation of TBX5. Together, these findings provide new insights into the mechanisms underlying LPM development, somatic LPM fate choice and early initiation of the vertebrate limb.

Runt-related transcription factor 2 (RUNX2) is a transcription factor essential for skeletal development. Variation within the RUNX2 polyglutamine / polyalanine (QA) repeat is correlated with facial length within orders of placental mammals and is suggested to be a major driver of craniofacial diversity. However, it is not known if this correlation exists outside of the placental mammals. Here we examined the correlation between the RUNX2 QA repeat ratio and facial length in the naturally evolving sister group to the placental mammals, the marsupials. Marsupials have a diverse range of facial lengths similar to that seen in placental mammals. Despite their diversity there was almost no variation seen in the RUNX2 QA repeat across individuals spanning the entire marsupial infraclass. The extreme conservation of the marsupial RUNX2 QA repeat indicates it is under strong purifying selection. Despite this, we observed an unexpectedly high level of repeat purity. Unlike within orders of placental mammals, RUNX2 repeat variation cannot drive craniofacial diversity in marsupials. We propose conservation of the marsupial RUNX2 QA repeat is driven by the constraint of accelerated ossification of the anterior skeleton to facilitate life in the pouch. Thus, marsupials must utilize alternate pathways to placental mammals to drive craniofacial evolution.