Initiation of speciation across multiple dimensions in a rock-restricted, tropical lizard.

Population isolation and concomitant genetic divergence, resulting in strong phylogeographical structure, is a core aspect of speciation initiation. If and how speciation then proceeds and ultimately completes depends on multiple factors that mediate reproductive isolation, including divergence in genomes, ecology and mating traits. Here we explored these multiple dimensions in two young (Plio-Pleistocene) species complexes of gekkonid lizards (Heteronotia) from the Kimberley-Victoria River regions of tropical Australia. Using mitochondrial DNA screening and exon capture phylogenomics, we show that the rock-restricted Heteronotia planiceps exhibits exceptional fine-scale phylogeographical structure compared to the codistributed habitat generalist Heteronotia binoei. This indicates pervasive population isolation and persistence in the rock-specialist, and thus a high rate of speciation initiation across this geographically complex region, with levels of genomic divergence spanning the "grey zone" of speciation. Proximal lineages of H. planiceps were often separated by different rock substrates, suggesting a potential role for ecological isolation; however, phylogenetic incongruence and historical introgression were inferred between one such pair. Ecomorphological divergence among lineages within both H. planiceps and H. binoei was limited, except that limestone-restricted lineages of H. planiceps tended to be larger than rock-generalists. By contrast, among-lineage divergence in the chemical composition of epidermal pore secretions (putative mating trait) exceeded ecomorphology in both complexes, but with less trait overlap among lineages in H. planiceps. This system-particularly the rock-specialist H. planiceps-highlights the role of multidimensional divergence during incipient speciation, with divergence in genomes, ecomorphology and chemical signals all at play at very fine spatial scales.

Composition of a chemical signalling trait varies with phylogeny and precipitation across an Australian lizard radiation.

The environment presents challenges to the transmission and detection of animal signalling systems, resulting in selective pressures that can drive signal divergence amongst populations in disparate environments. For chemical signals, climate is a potentially important selective force because factors such as temperature and moisture influence the persistence and detection of chemicals. We investigated an Australian lizard radiation (Heteronotia) to explore relationships between a sexually dimorphic chemical signalling trait (epidermal pore secretions) and two key climate variables: temperature and precipitation. We reconstructed the phylogeny of Heteronotia with exon capture phylogenomics, estimated phylogenetic signal in amongst-lineage chemical variation and assessed how chemical composition relates to temperature and precipitation using multivariate phylogenetic regressions. High estimates of phylogenetic signal indicate that the composition of epidermal pore secretions varies amongst lineages in a manner consistent with Brownian motion, although there are deviations to this, with stark divergences coinciding with two phylogenetic splits. Accounting for phylogenetic non-independence, we found that amongst-lineage chemical variation is associated with geographic variation in precipitation but not temperature. This contrasts somewhat with previous lizard studies, which have generally found an association between temperature and chemical composition. Our results suggest that geographic variation in precipitation can affect the evolution of chemical signalling traits, possibly influencing patterns of divergence amongst lineages and species.

Phylogenomics reveals accelerated late Cretaceous diversification of bee flies (Diptera: Bombyliidae).

Bombyliidae is a very species-rich and widespread family of parasitoid flies with more than 250 genera classified into 17 extant subfamilies. However, little is known about their evolutionary history or how their present-day diversity was shaped. Transcriptomes of 15 species and anchored hybrid enrichment (AHE) sequence captures of 86 species, representing 94 bee fly species and 14 subfamilies, were used to reconstruct the phylogeny of Bombyliidae. We integrated data from transcriptomes across each of the main lineages in our AHE tree to build a data set with more genes (550 loci versus 216 loci) and higher support levels. Our overall results show strong congruence with the current classification of the family, with 11 out of 14 included subfamilies recovered as monophyletic. Heterotropinae and Mythicomyiinae are successive sister groups to the remainder of the family. We examined the evolution of key morphological characters through our phylogenetic hypotheses and show that neither the "sand chamber subfamilies" nor the "Tomophthalmae" are monophyletic in our phylogenomic analyses. Based on our results, we reinstate two tribes at the subfamily level (Phthiriinae stat. rev. and Ecliminae stat. rev.) and we include the genus Sericosoma Macquart (previously incertae sedis) in the subfamily Oniromyiinae, bringing the total number of bee fly subfamilies to 19. Our dating analyses indicate a Jurassic origin of the family (165-194 Ma), with the sand chamber evolving early in bee fly evolution, in the late Jurassic or mid-Cretaceous (100-165 Ma). We hypothesize that the angiosperm radiation and the hothouse climate established during the late Cretaceous accelerated the diversification of bee flies, by providing an expanded range of resources for the parasitoid larvae and nectarivorous adults.

Cryptic lineage diversity, body size divergence, and sympatry in a species complex of Australian lizards (Gehyra).

Understanding the joint evolutionary and ecological underpinnings of sympatry among close relatives remains a key challenge in biology. This problem can be addressed through joint phylogenomic and phenotypic analysis of complexes of closely related lineages within, and across, species and hence representing the speciation continuum. For a complex of tropical geckos from northern Australia-Gehyra nana and close relatives-we combine mtDNA phylogeography, exon-capture sequencing, and morphological data to resolve independently evolving lineages and infer their divergence history and patterns of morphological evolution. Gehyra nana is found to include nine divergent lineages and is paraphyletic with four other species from the Kimberley region of north-west Australia. Across these 13 taxa, 12 of which are restricted to rocky habitats, several lineages overlap geographically, including on the diverse Kimberley islands. Morphological evolution is dominated by body size shifts, and both body size and shape have evolved gradually across the group. However, larger body size shifts are observed among overlapping taxa than among closely related parapatric lineages of G. nana, and sympatric lineages are more divergent than expected at random. Whether elevated body size differences among sympatric lineages are due to ecological sorting or character displacement remains to be determined.

Multiscale Evaluation of Thermal Dependence in the Glucocorticoid Response of Vertebrates.

Environmental temperature has profound effects on animal physiology, ecology, and evolution. Glucocorticoid (GC) hormones, through effects on phenotypic performance and life history, provide fundamental vertebrate physiological adaptations to environmental variation, yet we lack a comprehensive understanding of how temperature influences GC regulation in vertebrates. Using field studies and meta- and comparative phylogenetic analyses, we investigated how acute change and broadscale variation in temperature correlated with baseline and stress-induced GC levels. Glucocorticoid levels were found to be temperature and taxon dependent, but generally, vertebrates exhibited strong positive correlations with acute changes in temperature. Furthermore, reptile baseline, bird baseline, and capture stress-induced GC levels to some extent covaried with broadscale environmental temperature. Thus, vertebrate GC function appears clearly thermally influenced. However, we caution that lack of detailed knowledge of thermal plasticity, heritability, and the basis for strong phylogenetic signal in GC responses limits our current understanding of the role of GC hormones in species' responses to current and future climate variation.

Identification and qualification of 500 nuclear, single-copy, orthologous genes for the Eupulmonata (Gastropoda) using transcriptome sequencing and exon capture.

The qualification of orthology is a significant challenge when developing large, multiloci phylogenetic data sets from assembled transcripts. Transcriptome assemblies have various attributes, such as fragmentation, frameshifts and mis-indexing, which pose problems to automated methods of orthology assessment. Here, we identify a set of orthologous single-copy genes from transcriptome assemblies for the land snails and slugs (Eupulmonata) using a thorough approach to orthology determination involving manual alignment curation, gene tree assessment and sequencing from genomic DNA. We qualified the orthology of 500 nuclear, protein-coding genes from the transcriptome assemblies of 21 eupulmonate species to produce the most complete phylogenetic data matrix for a major molluscan lineage to date, both in terms of taxon and character completeness. Exon capture targeting 490 of the 500 genes (those with at least one exon >120 bp) from 22 species of Australian Camaenidae successfully captured sequences of 2825 exons (representing all targeted genes), with only a 3.7% reduction in the data matrix due to the presence of putative paralogs or pseudogenes. The automated pipeline Agalma retrieved the majority of the manually qualified 500 single-copy gene set and identified a further 375 putative single-copy genes, although it failed to account for fragmented transcripts resulting in lower data matrix completeness when considering the original 500 genes. This could potentially explain the minor inconsistencies we observed in the supported topologies for the 21 eupulmonate species between the manually curated and 'Agalma-equivalent' data set (sharing 458 genes). Overall, our study confirms the utility of the 500 gene set to resolve phylogenetic relationships at a range of evolutionary depths and highlights the importance of addressing fragmentation at the homolog alignment stage for probe design.