Museum skins enable identification of introgression associated with cytonuclear discordance.

Increased sampling of genomes and populations across closely related species has revealed that levels of genetic exchange during and after speciation are higher than previously thought. One obvious manifestation of such exchange is strong cytonuclear discordance, where the divergence in mitochondrial DNA (mtDNA) differs from that for nuclear genes more (or less) than expected from differences between mtDNA and nuclear DNA (nDNA) in population size and mutation rate. Given genome-scale datasets and coalescent modelling, we can now confidently identify cases of strong discordance and test specifically for historical or recent introgression as the cause. Using population sampling, combining exon capture data from historical museum specimens and recently collected tissues we showcase how genomic tools can resolve complex evolutionary histories in the brachyotis group of rock-wallabies (Petrogale). In particular, applying population and phylogenomic approaches we can assess the role of demographic processes in driving complex evolutionary patterns and assess a role of ancient introgression and hybridisation. We find that described species are well supported as monophyletic taxa for nDNA genes, but not for mtDNA, with cytonuclear discordance involving at least four operational taxonomic units (OTUs) across four species which diverged 183-278 kya. ABC modelling of nDNA gene trees supports introgression during or after speciation for some taxon pairs with cytonuclear discordance. Given substantial differences in body size between the species involved, this evidence for gene flow is surprising. Heterogenous patterns of introgression were identified but do not appear to be associated with chromosome differences between species. These and previous results suggest that dynamic past climates across the monsoonal tropics could have promoted reticulation among related species.

Functionally mediated cranial allometry evidenced in a genus of rock-wallabies.

In assessments of skeletal variation, allometry (disproportionate change of shape with size) is often corrected to examine size-independent variation for hypotheses relating to function. However, size-related trade-offs in functional demands may themselves be an underestimated driver of mammalian cranial diversity. Here, we use geometric morphometrics alongside dental measurements to assess craniodental allometry in the rock-wallaby genus Petrogale (all 17 species, 370 individuals). We identified functional aspects of evolutionary allometry that can be both extensions of, and correlated negatively with, static or ontogenetic allometric patterns. Regarding constraints, larger species tended to have relatively smaller braincases and more posterior orbits, the former of which might represent a constraint on jaw muscle anatomy. However, they also tended to have more anterior dentition and smaller posterior zygomatic arches, both of which support the hypothesis of relaxed bite force demands and accommodation of different selective pressures that favour facial elongation. By contrast, two dwarf species had stouter crania with divergent dental adaptations that together suggest increased relative bite force capacity. This likely allows them to feed on forage that is mechanically similar to that consumed by larger relatives. Our results highlight a need for nuanced considerations of allometric patterns in future research of mammalian cranial diversity.

Limited Introgression between Rock-Wallabies with Extensive Chromosomal Rearrangements.

Chromosome rearrangements can result in the rapid evolution of hybrid incompatibilities. Robertsonian fusions, particularly those with monobrachial homology, can drive reproductive isolation amongst recently diverged taxa. The recent radiation of rock-wallabies (genus Petrogale) is an important model to explore the role of Robertsonian fusions in speciation. Here, we pursue that goal using an extensive sampling of populations and genomes of Petrogale from north-eastern Australia. In contrast to previous assessments using mitochondrial DNA or nuclear microsatellite loci, genomic data are able to separate the most closely related species and to resolve their divergence histories. Both phylogenetic and population genetic analyses indicate introgression between two species that differ by a single Robertsonian fusion. Based on the available data, there is also evidence for introgression between two species which share complex chromosomal rearrangements. However, the remaining results show no consistent signature of introgression amongst species pairs and where evident, indicate generally low introgression overall. X-linked loci have elevated divergence compared with autosomal loci indicating a potential role for genic evolution to produce reproductive isolation in concert with chromosome change. Our results highlight the value of genome scale data in evaluating the role of Robertsonian fusions and structural variation in divergence, speciation, and patterns of molecular evolution.

Future-proofing the koala: Synergising genomic and environmental data for effective species management.

Climatic and evolutionary processes are inextricably linked to conservation. Avoiding extinction in rapidly changing environments often depends upon a species' capacity to adapt in the face of extreme selective pressures. Here, we employed exon capture and high-throughput next-generation sequencing to investigate the mechanisms underlying population structure and adaptive genetic variation in the koala (Phascolarctos cinereus), an iconic Australian marsupial that represents a unique conservation challenge because it is not uniformly threatened across its range. An examination of 250 specimens representing 91 wild source locations revealed that five major genetic clusters currently exist on a continental scale. The initial divergence of these clusters appears to have been concordant with the Mid-Brunhes Transition (~430 to 300 kya), a major climatic reorganisation that increased the amplitude of Pleistocene glacial-interglacial cycles. While signatures of polygenic selection and environmental adaptation were detected, strong evidence for repeated, climate-associated range contractions and demographic bottleneck events suggests that geographically isolated refugia may have played a more significant role in the survival of the koala through the Pleistocene glaciation than in situ adaptation. Consequently, the conservation of genome-wide genetic variation must be aligned with the protection of core koala habitat to increase the resilience of vulnerable populations to accelerating anthropogenic threats. Finally, we propose that the five major genetic clusters identified in this study should be accounted for in future koala conservation efforts (e.g., guiding translocations), as existing management divisions in the states of Queensland and New South Wales do not reflect historic or contemporary population structure.

A Framework for Resolving Cryptic Species: A Case Study from the Lizards of the Australian Wet Tropics.

As we collect range-wide genetic data for morphologically-defined species, we increasingly unearth evidence for cryptic diversity. Delimiting this cryptic diversity is challenging, both because the divergences span a continuum and because the lack of overt morphological differentiation suggests divergence has proceeded heterogeneously. Herein, we address these challenges as we diagnose and describe species in three co-occurring species groups of Australian lizards. By integrating genomic and morphological data with data on hybridization and introgression from contact zones, we explore several approaches-and their relative benefits and weaknesses-for testing the validity of cryptic lineages. More generally, we advocate that genetic delimitations of cryptic diversity must consider whether these lineages are likely to be durable and persistent through evolutionary time.

Analysis of Phylogenomic Tree Space Resolves Relationships Among Marsupial Families.

A fundamental challenge in resolving evolutionary relationships across the tree of life is to account for heterogeneity in the evolutionary signal across loci. Studies of marsupial mammals have demonstrated that this heterogeneity can be substantial, leaving considerable uncertainty in the evolutionary timescale and relationships within the group. Using simulations and a new phylogenomic data set comprising nucleotide sequences of 1550 loci from 18 of the 22 extant marsupial families, we demonstrate the power of a method for identifying clusters of loci that support different phylogenetic trees. We find two distinct clusters of loci, each providing an estimate of the species tree that matches previously proposed resolutions of the marsupial phylogeny. We also identify a well-supported placement for the enigmatic marsupial moles (Notoryctes) that contradicts previous molecular estimates but is consistent with morphological evidence. The pattern of gene-tree variation across tree-space is characterized by changes in information content, GC content, substitution-model adequacy, and signatures of purifying selection in the data. In a simulation study, we show that incomplete lineage sorting can explain the division of loci into the two tree-topology clusters, as found in our phylogenomic analysis of marsupials. We also demonstrate the potential benefits of minimizing uncertainty from phylogenetic conflict for molecular dating. Our analyses reveal that Australasian marsupials appeared in the early Paleocene, whereas the diversification of present-day families occurred primarily during the late Eocene and early Oligocene. Our methods provide an intuitive framework for improving the accuracy and precision of phylogenetic inference and molecular dating using genome-scale data.

Marsupial chromosomics: bridging the gap between genomes and chromosomes.

Marsupials have unique features that make them particularly interesting to study, and sequencing of marsupial genomes is helping to understand their evolution. A decade ago, it was a huge feat to sequence the first marsupial genome. Now, the advances in sequencing technology have made the sequencing of many more marsupial genomes possible. However, the DNA sequence is only one component of the structures it is packaged into: chromosomes. Knowing the arrangement of the DNA sequence on each chromosome is essential for a genome assembly to be used to its full potential. The importance of combining sequence information with cytogenetics has previously been demonstrated for rapidly evolving regions of the genome, such as the sex chromosomes, as well as for reconstructing the ancestral marsupial karyotype and understanding the chromosome rearrangements involved in the Tasmanian devil facial tumour disease. Despite the recent advances in sequencing technology assisting in genome assembly, physical anchoring of the sequence to chromosomes is required to achieve a chromosome-level assembly. Once chromosome-level assemblies are achieved for more marsupials, we will be able to investigate changes in the packaging and interactions between chromosomes to gain an understanding of the role genome architecture has played during marsupial evolution.

Reticulation, divergence, and the phylogeography-phylogenetics continuum.

Phylogeography, and its extensions into comparative phylogeography, have their roots in the layering of gene trees across geography, a paradigm that was greatly facilitated by the nonrecombining, fast evolution provided by animal mtDNA. As phylogeography moves into the era of next-generation sequencing, the specter of reticulation at several levels-within loci and genomes in the form of recombination and across populations and species in the form of introgression-has raised its head with a prominence even greater than glimpsed during the nuclear gene PCR era. Here we explore the theme of reticulation in comparative phylogeography, speciation analysis, and phylogenomics, and ask how the centrality of gene trees has fared in the next-generation era. To frame these issues, we first provide a snapshot of multilocus phylogeographic studies across the Carpentarian Barrier, a prominent biogeographic barrier dividing faunas spanning the monsoon tropics in northern Australia. We find that divergence across this barrier is evident in most species, but is heterogeneous in time and demographic history, often reflecting the taxonomic distinctness of lineages spanning it. We then discuss a variety of forces generating reticulate patterns in phylogeography, including introgression, contact zones, and the potential selection-driven outliers on next-generation molecular markers. We emphasize the continued need for demographic models incorporating reticulation at the level of genomes and populations, and conclude that gene trees, whether explicit or implicit, should continue to play a role in the future of phylogeography.

Residents' Reactions to Earthquake Early Warnings in Japan.

This article empirically examines the effectiveness of earthquake early warning (EEW) in Japan based on experiences of residents who received warnings before earthquake shaking occurred. In Study 1, a survey (N = 299) was conducted to investigate residents' experiences of, and reactions to, an EEW issued in Gunma and neighboring regions on June 17, 2018. The main results were as follows. (1) People's primary reactions to the EEW were mental, not physical, and thus motionless. Most residents stayed still, not for safety reasons, but because they were focusing on mentally bracing themselves. (2) Residents perceived the EEW to be effective because it enabled them to mentally prepare, rather than take physical protective actions, before strong shaking arrived. (3) In future, residents anticipate that on receipt of an EEW they would undertake mental preparation as opposed to physical protective actions. In Study 2, a survey (N = 450) was conducted on another EEW issued for an earthquake offshore of Chiba Prefecture on July 7, 2018. Results were in line with those of Study 1, suggesting that the findings described above are robust. Finally, given people's lack of impetus to undertake protective action on receipt of an EEW, this article discusses ways to enhance such actions.

Phylogenomics of a rapid radiation: the Australian rainbow skinks.

BACKGROUND: The application of target capture with next-generation sequencing now enables phylogenomic analyses of rapidly radiating clades of species. But such analyses are complicated by extensive incomplete lineage sorting, demanding the use of methods that consider this process explicitly, such as the multispecies coalescent (MSC) model. However, the MSC makes strong assumptions about divergence history and population structure, and when using the full Bayesian implementation, current computational limits mean that relatively few loci and samples can be analysed for even modest sized radiations. We explore these issues through analyses of an extensive (> 1000 loci) dataset for the Australian rainbow skinks. This group consists of 3 genera and 41 described species, which likely diversified rapidly in Australia during the mid-late Miocene to occupy rainforest, woodland, and rocky habitats with corresponding diversity of morphology and breeding colouration. Previous phylogenetic analyses of this group have revealed short inter-nodes and high discordance among loci, limiting the resolution of inferred trees. A further complication is that many species have deep phylogeographic structure - this poses the question of how to sample individuals within species for analyses using the MSC. RESULTS: Phylogenies obtained using concatenation and summary coalescent species tree approaches to the full dataset are well resolved with generally consistent topology, including for previously intractable relationships near the base of the clade. As expected, branch lengths at the tips are substantially overestimated using concatenation. Comparisons of different strategies for sampling haplotypes for full Bayesian MSC analyses (for one clade and using smaller sets of loci) revealed, unexpectedly, that combining haplotypes across divergent phylogeographic lineages yielded consistent species trees. CONCLUSIONS: This study of more than 1000 loci provides a strongly-supported estimate of the phylogeny of the Australian rainbow skinks, which will inform future research on the evolution and taxonomy of this group. Our analyses suggest that species tree estimation with the MSC can be quite robust to violation of the assumption that the individuals representing a taxon are sampled from a panmictic population.

Pleistocene climatic changes drive diversification across a tropical savanna.

Spatial responses of species to past climate change depend on both intrinsic traits (climatic niche breadth, dispersal rates) and the scale of climatic fluctuations across the landscape. New capabilities in generating and analysing population genomic data, along with spatial modelling, have unleashed our capacity to infer how past climate changes have shaped populations, and by extension, complex communities. Combining these approaches, we uncover lineage diversity across four codistributed lizards from the Australian Monsoonal Tropics and explore how varying climatic tolerances interact with regional climate history to generate common vs. disparate responses to late Pleistocene change. We find more divergent spatial structuring and temporal demographic responses in the drier Kimberley region compared to the more mesic and consistently suitable Top End. We hypothesize that, in general, the effects of species' traits on sensitivity to climate fluctuation will be more evident in climatically marginal regions. If true, this points to the need in climatically marginal areas to craft more species-(or trait)-specific strategies for persistence under future climate change.

Phylogenetic analysis of the tree-kangaroos (Dendrolagus) reveals multiple divergent lineages within New Guinea.

Amongst the Australasian kangaroos and wallabies (Macropodidae) one anomalous genus, the tree-kangaroos, Dendrolagus, has secondarily returned to arboreality. Modern tree-kangaroos are confined to the wet tropical forests of north Queensland, Australia (2 species) and New Guinea (8 species). Due to their behavior, distribution and habitat most species are poorly known and our understanding of the evolutionary history and systematics of the genus is limited and controversial. We obtained tissue samples from 36 individual Dendrolagus including representatives from 14 of the 17 currently recognised or proposed subspecies and generated DNA sequence data from three mitochondrial (3116 bp) and five nuclear (4097 bp) loci. Phylogenetic analysis of these multi-locus data resolved long-standing questions regarding inter-relationships within Dendrolagus. The presence of a paraphyletic ancestral long-footed and derived monophyletic short-footed group was confirmed. Six major lineages were identified: one in Australia (D. lumholtzi, D. bennettianus) and five in New Guinea (D. inustus, D. ursinus, a Goodfellow's group, D. mbaiso and a Doria's group). Two major episodes of diversification within Dendrolagus were identified: the first during the late Miocene/early Pliocene associated with orogenic processes in New Guinea and the second mostly during the early Pleistocene associated with the intensification of climatic cycling. All sampled subspecies showed high levels of genetic divergence and currently recognized species within both the Doria's and Goodfellow's groups were paraphyletic indicating that adjustments to current taxonomy are warranted.

Accounting for Uncertainty in Gene Tree Estimation: Summary-Coalescent Species Tree Inference in a Challenging Radiation of Australian Lizards.

Accurate gene tree inference is an important aspect of species tree estimation in a summary-coalescent framework. Yet, in empirical studies, inferred gene trees differ in accuracy due to stochastic variation in phylogenetic signal between targeted loci. Empiricists should, therefore, examine the consistency of species tree inference, while accounting for the observed heterogeneity in gene tree resolution of phylogenomic data sets. Here, we assess the impact of gene tree estimation error on summary-coalescent species tree inference by screening ${\sim}2000$ exonic loci based on gene tree resolution prior to phylogenetic inference. We focus on a phylogenetically challenging radiation of Australian lizards (genus Cryptoblepharus, Scincidae) and explore effects on topology and support. We identify a well-supported topology based on all loci and find that a relatively small number of high-resolution gene trees can be sufficient to converge on the same topology. Adding gene trees with decreasing resolution produced a generally consistent topology, and increased confidence for specific bipartitions that were poorly supported when using a small number of informative loci. This corroborates coalescent-based simulation studies that have highlighted the need for a large number of loci to confidently resolve challenging relationships and refutes the notion that low-resolution gene trees introduce phylogenetic noise. Further, our study also highlights the value of quantifying changes in nodal support across locus subsets of increasing size (but decreasing gene tree resolution). Such detailed analyses can reveal anomalous fluctuations in support at some nodes, suggesting the possibility of model violation. By characterizing the heterogeneity in phylogenetic signal among loci, we can account for uncertainty in gene tree estimation and assess its effect on the consistency of the species tree estimate. We suggest that the evaluation of gene tree resolution should be incorporated in the analysis of empirical phylogenomic data sets. This will ultimately increase our confidence in species tree estimation using summary-coalescent methods and enable us to exploit genomic data for phylogenetic inference. [Coalescence; concatenation; Cryptoblepharus; exon capture; gene tree; phylogenomics; species tree.].

Tropical specialist vs. climate generalist: Diversification and demographic history of sister species of Carlia skinks from northwestern Australia.

Species endemic to the tropical regions are expected to be vulnerable to future climate change due in part to their relatively narrow climatic niches. In addition, these species are more likely to have responded strongly to past climatic change, and this can be explored through phylogeographic analyses. To test the hypothesis that tropical specialists are more sensitive to climate change than climate generalists, we generated and analyse sequence data from mtDNA and ~2500 exons to compare scales of historical persistence and population fluctuation in two sister species of Australian rainbow skinks: the tropical specialist Carlia johnstonei and the climate generalist C. triacantha. We expect the tropical specialist species to have deeper and finer-scale phylogeographic structure and stronger demographic fluctuations relative to the closely related climate generalist species, which should have had more stable populations through periods of harsh climate in the late Quaternary. Within C. johnstonei, we find that some populations from the northern Kimberley islands are highly divergent from mainland populations. In C. triacantha, one major clade occurs across the deserts and into the mesic Top End, and another occurs primarily in the Kimberley with scattered records eastwards. Where their ranges overlap in the Kimberley, both mitochondrial DNA and nuclear DNA suggest stronger phylogeographic structure and range expansion within the tropical specialist, whereas the climate generalist has minimal structuring and no evidence of recent past range expansion. These results are consistent with the hypothesis that tropical specialists are more sensitive to past climatic change.

Phylogenomics at the tips: inferring lineages and their demographic history in a tropical lizard, Carlia amax.

High-throughput sequencing approaches offer opportunities to better understand the evolutionary processes driving diversification, particularly in nonmodel organisms. In particular, the 100-1000's of loci that can now be sequenced are providing unprecedented power in population, speciation and phylogenetic studies. Here, we apply an exon capture approach to generate >99% complete sequence and SNP data across >2000 loci from a tropical skink, Carlia amax, and exploit these data to identify divergent lineages and infer their relationships and demographic histories. This is especially relevant to low-dispersal tropical taxa that often have cryptic diversity and spatially dynamic histories. For C. amax, clustering of nuclear SNPs and coalescent-based species delimitation analyses identify four divergent lineages, one fewer than predicted based on geographically coherent mtDNA clades (>9.4% sequence divergence). Three of these lineages are widespread and parapatric on the mainland, whereas the most divergent is restricted to islands off the northeast Northern Territory. Tests for population expansion reject an equilibrium isolation-by-distance model for two of the three widespread lineages and infer refugial expansion sources in the relatively mesic northeast Top End and northwest Kimberley. The latter is already recognized as a hotspot of endemism, but our results also suggest that a stronger focus on the northeast Top End, and adjacent islands is warranted. More generally, our results show how genome-reduction methods such as exon capture can yield insights into the pattern and dynamics of biodiversity across complex landscapes with as yet poorly understood biogeographic history and how exon data can link between population and phylogenetic questions.

Gene flow despite complex Robertsonian fusions among rock-wallaby (Petrogale) species.

Complex Robertsonian rearrangements, with shared arms in different fusions, are expected to prevent gene flow between hybrids through missegregation during meiosis. Here, we estimate gene flow between recently diverged and chromosomally diverse rock-wallabies (Petrogale) to test for this form of chromosomal speciation. Contrary to expectations, we observe relatively high admixture among species with complex fusions. Our results reinforce the need to consider alternative roles of chromosome change, together with genic divergence, in driving speciation.

The importance of an evolutionary perspective in conservation policy planning.

Prioritization of taxa for conservation must rest on a foundation of correctly identified species boundaries, enhanced by an understanding of evolutionary history and phylogenetic relationships. Therefore, we can incorporate both evolutionary and ecological processes into efforts to sustain biodiversity. In this issue of Molecular Ecology, Malaney & Cook (2013) highlight the critical value of an evolutionary biogeographical approach, combining multilocus phylogeography with climatic niche modelling to infer phylogenetically weighted conservation priorities for evolutionary lineages of jumping mice across North America. Remarkably, they find that the Preble's meadow jumping mouse (Zapus hudsonius preblei), long debated as a threatened taxon, in fact represents the southern terminus of a relatively uniform lineage that expanded well into Alaska during the Holocene. By contrast, some other relictual and phylogenetically divergent taxa of jumping mice likely warrant greater conservation priority. This study highlights the value of integrative approaches that place current taxonomy in a broader evolutionary context to identify taxa for conservation assessment, but also highlights the challenges in maintaining potential for adaptive responses to environmental change.

Multiple biogeographical barriers identified across the monsoon tropics of northern Australia: phylogeographic analysis of the brachyotis group of rock-wallabies.

The monsoon tropics of northern Australia are a globally significant biodiversity hotspot, but its phylogeography is poorly known. A major challenge for this region is to understand the biogeographical processes that have shaped the distribution and diversity of taxa, without detailed knowledge of past climatic and environmental fluctuations. Although molecular data have great potential to address these questions, only a few species have been examined phylogeographically. Here, we use the widely distributed and abundant short-eared rock-wallaby (Petrogale brachyotis; n = 101), together with the sympatric monjon (P. burbidgei; n = 11) and nabarlek (P. concinna; n = 1), to assess historical evolutionary and biogeographical processes in northern Australia. We sequenced ∼1000 bp of mitochondrial DNA (control region, ND2) and ∼3000 bp of nDNA (BRCA1, ω-globin and two anonymous loci) to investigate phylogeographic structuring and delineate the time-scale of diversification within the region. Our results indicate multiple barriers between the Top End (Northern Territory) and Kimberley (Western Australia), which have caused divergence throughout the Plio-Pleistocene. Eight geographically discrete and genetically distinct lineages within the brachyotis group were identified, five of which are separated by major river valleys (Ord, Victoria, Daly), arid lowlands and discontinuous sandstone ranges. It is likely that these barriers have similarly influenced genetic structure in other monsoonal biota.

Phylogenetic relationships of rock-wallabies, Petrogale (Marsupialia: Macropodidae) and their biogeographic history within Australia.

The rock-wallaby genus Petrogale comprises a group of habitat-specialist macropodids endemic to Australia. Their restriction to rocky outcrops, with infrequent interpopulation dispersal, has been suggested as the cause of their recent and rapid diversification. Molecular phylogenetic relationships within and among species of Petrogale were analysed using mitochondrial (cytochrome oxidase c subunit 1, cytochrome b, NADH dehydrogenase subunit 2) and nuclear (omega-globin intron, breast and ovarian cancer susceptibility gene) sequence data with representatives that encompassed the morphological and chromosomal variation within the genus, including for the first time both Petrogale concinna and Petrogale purpureicollis. Four distinct lineages were identified, (1) the brachyotis group, (2) Petrogale persephone, (3) Petrogalexanthopus and (4) the lateralis-penicillata group. Three of these lineages include taxa with the ancestral karyotype (2n=22). Paraphyletic relationships within the brachyotis group indicate the need for a focused phylogeographic study. There was support for P. purpureicollis being reinstated as a full species and P. concinna being placed within Petrogale rather than in the monotypic genus Peradorcas. Bayesian analyses of divergence times suggest that episodes of diversification commenced in the late Miocene-Pliocene and continued throughout the Pleistocene. Ancestral state reconstructions suggest that Petrogale originated in a mesic environment and dispersed into more arid environments, events that correlate with the timing of radiations in other arid zone vertebrate taxa across Australia.

Temperature predicts the rate of molecular evolution in Australian Eugongylinae skinks.

Temperature differences over time and space have been hypothesized to cause variation in the rate of molecular evolution of species, but empirical evidence is mixed. To further test this hypothesis, we utilized a large exon-capture sequence data of Australian Eugongylinae skinks, exemplifying a radiation of temperature-sensitive ectotherms spanning a large latitudinal gradient. The association between temperature (and other species traits) and long-term substitution rate was assessed based on 1268 sequenced exons of 44 species pairs from the Eugongylinae subfamily using regression analyses. Temperature is the strongest, positively correlated predictor of variation in substitution rate across the Australian Eugongylinae. It explains 45% of variation in synonymous substitution rate, and 11% after controlling for all the other factors. Synonymous substitution rate is also negatively associated with body size, with a 6% variation explained by body size after controlling for the effects of temperature. Other factors are not associated with synonymous substitution rate after controlling for temperature. Overall, this study points to temperature as a strong predictor of the molecular evolution rate in the Eugongylinae subfamily, and demonstrates the power of large-scale exonic data to identify correlates of the rate of molecular evolution.