Population genomics and sexual signals support reproductive character displacement in Uperoleia (Anura: Myobatrachidae) in a contact zone.

When closely related species come into contact via range expansion, both may experience reduced fitness as a result of the interaction. Selection is expected to favour traits that minimize costly interspecies reproductive interactions (such as mismating) via a phenomenon called reproductive character displacement (RCD). Research on RCD frequently assumes secondary contact between species, but the geographical history of species interactions is often unknown. Population genomic data permit tests of geographical hypotheses about species origins and secondary contact through range expansion. We used population genomic data from single nucleotide polymorphisms (SNPs), mitochondrial sequence data, advertisement call data and morphological data to investigate a species complex of toadlets (Uperoleia borealis, U. crassa, U. inundata) from northern Australia. Although the three species of frogs were morphologically indistinguishable in our analysis, we determined that U. crassa and U. inundata form a single species (synonymized here) based on an absence of genomic divergence. SNP data identified the phylogeographical origin of U. crassa as the Top End, with subsequent westward invasion into the range of U. borealis in the Kimberley. We identified six F1 hybrids, all of which had the U. borealis mitochondrial haplotype, suggesting unidirectional hybridization. Consistent with the RCD hypothesis, U. borealis and U. crassa sexual signals differ more in sympatry than in allopatry. Hybrid males have intermediate calls, which probably reduces attractiveness to females. Integrating population genomic data, mitochondrial sequencing, morphology and behavioural approaches provides an unusually detailed collection of evidence for reproductive character displacement following range expansion and secondary contact.

On and off the rocks: persistence and ecological diversification in a tropical Australian lizard radiation.

BACKGROUND: Congruent patterns in the distribution of biodiversity between regions or habitats suggest that key factors such as climatic and topographic variation may predictably shape evolutionary processes. In a number of tropical and arid biomes, genetic analyses are revealing deeper and more localised lineage diversity in rocky ranges than surrounding habitats. Two potential drivers of localised endemism in rocky areas are refugial persistence through climatic change, or ecological diversification and specialisation. Here we examine how patterns of lineage and phenotypic diversity differ across two broad habitat types (rocky ranges and open woodlands) in a small radiation of gecko lizards in the genus Gehyra (the australis group) from the Australian Monsoonal Tropics biome. RESULTS: Using a suite of approaches for delineating evolutionarily independent lineages, we find between 26 and 41 putative evolutionary units in the australis group (versus eight species currently recognised). Rocky ranges are home to a greater number of lineages that are also relatively more restricted in distribution, while lineages in open woodland habitats are fewer, more widely distributed, and, in one case, show evidence of range expansion. We infer at least two shifts out of rocky ranges and into surrounding woodlands. Phenotypic divergence between rocky ranges specialist and more generalist taxa is detected, but no convergent evolutionary regimes linked to ecology are inferred. CONCLUSIONS: In climatically unstable biomes such as savannahs, rocky ranges have functioned as zones of persistence, generators of diversity and a source of colonists for surrounding areas. Phenotypic divergence can also be linked to the use of differing habitat types, however, the extent to which ecological specialisation is a primary driver or secondary outcome of localised diversification remains uncertain.

The diverse liver viromes of Australian geckos and skinks are dominated by hepaciviruses and picornaviruses and reflect host taxonomy and habitat.

Lizards have diverse ecologies and evolutionary histories, and represent a promising group to explore how hosts shape virome structure and virus evolution. Yet, little is known about the viromes of these animals. In Australia, squamates (lizards and snakes) comprise the most diverse order of vertebrates, and Australia hosts the highest diversity of lizards globally, with the greatest breadth of habitat use. We used meta-transcriptomic sequencing to determine the virome of nine co-distributed, tropical lizard species from three taxonomic families in Australia and analyzed these data to identify host traits associated with viral abundance and diversity. We show that lizards carry a large diversity of viruses, identifying more than thirty novel, highly divergent vertebrate-associated viruses. These viruses were from nine viral families, including several that contain well known pathogens, such as the Flaviviridae, Picornaviridae, Bornaviridae, Iridoviridae, and Rhabdoviridae. Members of the Flaviviridae were particularly abundant across species sampled here, largely belonging to the genus Hepacivirus: fourteen novel hepaciviruses were identified, broadening the known diversity of this group and better defining its evolution by uncovering new reptilian clades. The evolutionary histories of the viruses studied here frequently aligned with the biogeographic and phylogenetic histories of the hosts, indicating that exogenous viruses may help infer host evolutionary history if sampling is strategic and sampling density high enough. Notably, analysis of alpha and beta diversity revealed that virome composition and richness in the animals sampled here was shaped by host taxonomy and habitat. In sum, we identified a diverse range of reptile viruses that broadly contributes to our understanding of virus-host ecology and evolution.

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.

The Potential for Rapid Evolution under Anthropogenic Climate Change.

Understanding how natural populations will respond to rapid anthropogenic climate change is one of the greatest challenges for ecologists and evolutionary biologists. Much research has focussed on whether physiological traits can evolve quickly enough under rapidly increasing temperatures. While the simple Breeder's equation helps to understand how extreme temperatures and genetic variation might drive within-population evolution under climate change, it does not consider two key areas: how different forms of phenotypic plasticity interact and variation among populations. Plasticity can modify the exposure to climatic extremes and the strength of selection from those extremes, while differences among populations provide adaptive diversity not apparent within them. Here, we focus on terrestrial vertebrates and, with a case study on a tropical lizard, demonstrate the complex interplay between spatial, genetic and plastic contributions to variation in climate-relevant physiological traits. We identify several problems that need to be better understood: which traits are under selection in a changing climate; the different forms of plasticity relevant to population persistence and rapid evolution; plastic versus genetic contributions to geographic variation in climate-associated traits and whether plasticity can be harnessed to promote persistence of species. Given ongoing uncertainties around whether natural populations can evolve rapidly enough to persist, we advocate the use of field trials aimed at increasing rates of adaptation, especially in systems known to be strongly impacted by human-driven changes in climate.

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.

Phylodiversity to inform conservation policy: An Australian example.

Phylodiversity measures summarise the phylogenetic diversity patterns of groups of organisms. By using branches of the tree of life, rather than its tips (e.g., species), phylodiversity measures provide important additional information about biodiversity that can improve conservation policy and outcomes. As a biodiverse nation with a strong legislative and policy framework, Australia provides an opportunity to use phylogenetic information to inform conservation decision-making. We explored the application of phylodiversity measures across Australia with a focus on two highly biodiverse regions, the south west of Western Australia (SWWA) and the South East Queensland bioregion (SEQ). We analysed seven diverse groups of organisms spanning five separate phyla on the evolutionary tree of life, the plant genera Acacia and Daviesia, mammals, hylid frogs, myobatrachid frogs, passerine birds, and camaenid land snails. We measured species richness, weighted species endemism (WE) and two phylodiversity measures, phylogenetic diversity (PD) and phylogenetic endemism (PE), as well as their respective complementarity scores (a measure of gains and losses) at 20 km resolution. Higher PD was identified within SEQ for all fauna groups, whereas more PD was found in SWWA for both plant groups. PD and PD complementarity were strongly correlated with species richness and species complementarity for most groups but less so for plants. PD and PE were found to complement traditional species-based measures for all groups studied: PD and PE follow similar spatial patterns to richness and WE, but highlighted different areas that would not be identified by conventional species-based biodiversity analyses alone. The application of phylodiversity measures, particularly the novel weighted complementary measures considered here, in conservation can enhance protection of the evolutionary history that contributes to present day biodiversity values of areas. Phylogenetic measures in conservation can include important elements of biodiversity in conservation planning, such as evolutionary potential and feature diversity that will improve decision-making and lead to better biodiversity conservation outcomes.

High genetic diversity despite the potential for stepping-stone colonizations in an invasive species of gecko on Moorea, French Polynesia.

Invasive species often have reduced genetic diversity, but the opposite can be true if there have been multiple introductions and genetic admixture. Reduced diversity is most likely soon after establishment, in remote locations, when there is lower propagule pressure and with stepping-stone colonizations. The common house gecko (Hemidactylus frenatus) was introduced to Moorea, French Polynesia in the remote eastern Pacific within the last two decades and accordingly is expected to exhibit low diversity. In contrast, we show that H. frenatus on Moorea has exceptionally high genetic diversity, similar to that near the native range in Asia and much higher than reported for other Pacific island reptiles. The high diversity in this recently founded population likely reflects extensive genetic admixture in source population(s) and a life history that promotes retention of diversity. These observations point to the importance of understanding range-wide dynamics of genetic admixture in highly invasive species.

Museum genomics: low-cost and high-accuracy genetic data from historical specimens.

Natural history collections are unparalleled repositories of geographical and temporal variation in faunal conditions. Molecular studies offer an opportunity to uncover much of this variation; however, genetic studies of historical museum specimens typically rely on extracting highly degraded and chemically modified DNA samples from skins, skulls or other dried samples. Despite this limitation, obtaining short fragments of DNA sequences using traditional PCR amplification of DNA has been the primary method for genetic study of historical specimens. Few laboratories have succeeded in obtaining genome-scale sequences from historical specimens and then only with considerable effort and cost. Here, we describe a low-cost approach using high-throughput next-generation sequencing to obtain reliable genome-scale sequence data from a traditionally preserved mammal skin and skull using a simple extraction protocol. We show that single-nucleotide polymorphisms (SNPs) from the genome sequences obtained independently from the skin and from the skull are highly repeatable compared to a reference genome.

MITOCHONDRIAL-DNA ANALYSES AND THE ORIGIN AND RELATIVE AGE OF PARTHENOGENETIC LIZARDS (GENUS CNEMIDOPHORUS). III. C. VELOX AND C. EXSANGUIS.

Mitochondrial DNAs (mtDNAs) of two unisexual, parthenogenetically reproducing species of whiptail lizards (Cnemidophorus velox and C. exsanguis) and their bisexual relatives were compared by restriction-enzyme analysis to assess levels of mtDNA variation and to establish the maternal ancestry of the unisexuals. No cleavage-site differences were found to be diagnostic between C. velox and C. exsanguis mtDNAs, suggesting an ancestry rooted in the same maternal lineage. The mtDNA of the unisexuals is relatively homogeneous, indicating that these lineages are of recent origin. Phylogenetic analysis revealed that the maternal ancestor of both C. velox and C. exsanguis was most probably C. burti stictogrammus, C. costatus barrancorum, or an unidentified taxon closely related to them. In addition, the mtDNA analyses demonstrate conclusively that the triploid species C. velox could not have been formed by the fertilization of an unreduced (diploid) C. inornatus egg, further strengthening the hypothesis that parthenogenesis in Cnemidophorus results from hybridization.

MITOCHONDRIAL-DNA DIVERSITY AND THE ORIGIN OF THE MENIDIA CLARKHUBBSI COMPLEX OF UNISEXUAL FISHES (ATHERINIDAE).

Restriction-endonuclease analysis of mitochondrial DNA (mtDNA) of the unisexual M. clarkhubbsi complex and close sexual relatives indicated that the unisexual complex arose through multiple, nonreciprocal hybridizations involving females of M. peninsulae. High-resolution analyses using restriction endonucleases that cleave at 4-bp sites revealed mtDNA sequence diversity that was low among unisexuals but high among individuals of M. peninsulae. The identification of M. peninsulae as a parent of the unisexuals conflicts with some details of previous allozyme comparisons. One possibility is that the unisexuals were derived from hybridization involving M. beryllina and a recently extinct form of M. peninsulae. In contrast to the unisexuals, contemporary hybrids of M. peninsulae and M. beryllina are formed by reciprocal matings. The origins of extant unisexual lineages from nonreciprocal hybridizations, together with their low mtDNA diversity relative to the maternal ancestor, implies strong constraints on origins of unisexuality via hybridization. Data on reproduction by contemporary F1 hybrids reveal one form of genetic/developmental constraint: M. peninsulae and M. beryllina may now have diverged beyond the point where the hybrid origin of new unisexual lineages is possible.

MITOCHONDRIAL-DNA ANALYSES AND THE ORIGIN AND RELATIVE AGE OF PARTHENOGENETIC LIZARDS (GENUS CNEMIDOPHORUS). IV. NINE SEXLINEATUS-GROUP UNISEXUALS.

Mitochondrial DNAs (mtDNAs) from nine morphologically distinct unisexual species and five bisexual species of lizards, all from the sexlineatus species-group of Cnemidophorus, were compared using restriction endonucleases. The unisexual lizards have mtDNAs that are identical at all or nearly all of the 128 sites cleaved. Although differing little in sequence, some mtDNAs differed in size due to the presence of tandem sequence duplications. Phylogenetic analysis of cleavage maps indicates that the mtDNAs of the unisexuals are most similar to that of the bisexual species C. inornatus. Considerable mtDNA diversity exists among C. inornatus populations, and one geographically restricted subspecies, C. i. arizonae, was identified as the most probable maternal ancestor of all nine unisexuals. All but one of these are triploid, and all have at least one C. inornatus gene complement. This, together with the homogeneity of their mtDNAs, suggests that all stem from one or a small number of allodiploid females (presumably parthenogenetic) that originated in a restricted geographic area in the recent past. These data, when combined with those from allozyme studies, preclude the possibility that most of the triploid unisexuals could have arisen via fertilization of an unreduced diploid ovum from one species by a haploid sperm from a different species.