Evolution across the adaptive landscape in a hyperdiverse beetle radiation

The extraordinary diversification of beetles on Earth is a textbook example of adaptive evolution. Yet, the tempo and drivers of this super-radiation remain largely unclear. Here, we address this problem by investigating macroevolutionary dynamics in darkling beetles (Coleoptera: Tenebrionidae), one of the most ecomorphologically diverse beetle families (with over 30,000 species). Using multiple genomic datasets and analytical approaches, we resolve the long-standing inconsistency over deep relationships in the family. In conjunction with a landmark-based dataset of body shape morphology, we show that the evolutionary history of darkling beetles is marked by ancient rapid radiations, frequent ecological transitions, and rapid bursts of morphological diversification. On a global scale, our analyses uncovered a significant pulse of phenotypic diversification proximal to the Cretaceous-Palaeogene (K/Pg) mass extinction and convergence of body shape associated with recurrent ecological specializations. On a regional scale, two major Australasian radiations, the Adeliini and the Heleine clade, exhibited contrasting patterns of ecomorphological diversification, representing phylogenetic niche conservatism versus adaptive radiation. Our findings align with the Simpsonian model of adaptive evolution across the macroevolutionary landscape and highlight a significant role of ecological opportunity in driving the immense ecomorphological diversity in a hyperdiverse beetle group.

GALV-KoRV-related retroviruses in diverse Australian and African rodent species

The enigmatic origins and transmission events of the gibbon ape leukemia virus (GALV) and its close relative the koala retrovirus (KoRV) have been a source of enduring debate. Bats and rodents are each proposed as major reservoirs of interspecies transmission, with ongoing efforts to identify additional animal hosts of GALV-KoRV-related retroviruses. In this study, we identified nine rodent species as novel hosts of GALV-KoRV-related retroviruses. Included among these hosts are two African rodents, revealing the first appearance of this clade beyond the Australian and Southeast Asian region. One of these African rodents, Mastomys natalensis, carries an endogenous GALV-KoRV-related retrovirus that is fully intact and potentially still infectious. Our findings support the hypothesis that rodents are the major carriers of GALV-KoRV-related retroviruses.

Discovering Fragile Clades and Causal Sequences in Phylogenomics by Evolutionary Sparse Learning

Phylogenomic analyses of long sequences, consisting of many genes and genomic segments, reconstruct organismal relationships with high statistical confidence. But, inferred relationships can be sensitive to excluding just a few sequences. Currently, there is no direct way to identify fragile relationships and the associated individual gene sequences in species. Here, we introduce novel metrics for gene-species sequence concordance and clade probability derived from evolutionary sparse learning models. We validated these metrics using fungi, plant, and animal phylogenomic datasets, highlighting the ability of the new metrics to pinpoint fragile clades and the sequences responsible. The new approach does not necessitate the investigation of alternative phylogenetic hypotheses, substitution models, or repeated data subset analyses. Our methodology offers a streamlined approach to evaluating major inferred clades and identifying sequences that may distort reconstructed phylogenies using large datasets.

Phylogenomic analysis and molecular identification of true fruit flies

The family Tephritidae in the order Diptera, known as true fruit flies, are agriculturally important insect pests. However, the phylogenetic relationships of true fruit flies, remain controversial. Moreover, rapid identification of important invasive true fruit flies is essential for plant quarantine but is still challenging. To this end, we sequenced the genome of 16 true fruit fly species at coverage of 47-228×. Together with the previously reported genomes of nine species, we reconstructed phylogenetic trees of the Tephritidae using benchmarking universal single-copy ortholog (BUSCO), ultraconserved element (UCE) and anchored hybrid enrichment (AHE) gene sets, respectively. The resulting trees of 50% taxon-occupancy dataset for each marker type were generally congruent at 88% nodes for both concatenation and coalescent analyses. At the subfamily level, both Dacinae and Trypetinae are monophyletic. At the species level, Bactrocera dorsalis is more closely related to Bactrocera latifrons than Bactrocera tryoni. This is inconsistent with previous conclusions based on mitochondrial genes but consistent with recent studies based on nuclear data. By analyzing these genome data, we screened ten pairs of species-specific primers for molecular identification of ten invasive fruit flies, which PCR validated. In summary, our work provides draft genome data of 16 true fruit fly species, addressing the long-standing taxonomic controversies and providing species-specific primers for molecular identification of invasive fruit flies.

Incomplete lineage sorting and gene flow within Allium (Amayllidaceae)

The phylogeny and systematics of the genus Allium have been studied with a variety of diverse data types, including an increasing amount of molecular data. However, strong phylogenetic discordance and high levels of uncertainty have prevented the identification of a consistent phylogeny. The difficulty in establishing phylogenetic consensus and evidence for genealogical discordance make Allium a compelling test case to assess the relative contribution of incomplete lineage sorting (ILS), gene flow and gene tree estimation error on phylogenetic reconstruction. In this study, we obtained 75 transcriptomes of 38 Allium species across 10 subgenera. Whole plastid genome, single copy genes and consensus CDS were generated to estimate phylogenetic trees both using coalescence and concatenation methods. Multiple approaches including coalescence simulation, quartet sampling, reticulate network inference, sequence simulation, theta of ILS and reticulation index were carried out across the CDS gene trees to investigate the degrees of ILS, gene flow and gene tree estimation error. Afterward, a regression analysis was used to test the relative contributions of each of these forms of uncertainty to the final phylogeny. Despite extensive topological discordance among gene trees, we found a fully supported species tree that agrees with the most of well-accepted relationships and establishes monophyly of the genus Allium. We presented clear evidence for substantial ILS across the phylogeny of Allium. Further, we identified two ancient hybridization events for the formation of the second evolutionary line and subg. Butomissa as well as several introgression events between recently diverged species. Our regression analysis revealed that gene tree inference error and gene flow were the two most dominant factors explaining for the overall gene tree variation, with the difficulty in disentangling the effects of ILS and gene tree estimation error due to a positive correlation between them. Based on our efforts to mitigate the methodological errors in reconstructing trees, we believed ILS and gene flow are two principal reasons for the oft-reported phylogenetic heterogeneity of Allium. This study presents a strongly-supported and well-resolved phylogenetic backbone for the sampled Allium species, and exemplifies how to untangle heterogeneity in phylogenetic signal and reconstruct the true evolutionary history of the target taxa.

Populating a Continent: Phylogenomics Reveal the Timing of Australian Frog Diversification

The Australian continent's size and isolation make it an ideal place for studying the accumulation and evolution of biodiversity. Long separated from the ancient supercontinent Gondwana, most of Australia's plants and animals are unique and endemic, including the continent's frogs. Australian frogs comprise a remarkable ecological and morphological diversity categorized into a small number of distantly related radiations. We present a phylogenomic hypothesis based on an exon-capture dataset that spans the main clades of Australian myobatrachoid, pelodryadid hyloid, and microhylid frogs. Our time-calibrated phylogenomic-scale phylogeny identifies great disparity in the relative ages of these groups that vary from Gondwanan relics to recent immigrants from Asia and include arguably the continent's oldest living vertebrate radiation. This age stratification provides insight into the colonization of, and diversification on, the Australian continent through deep time, during periods of dramatic climatic and community changes. Contemporary Australian frog diversity highlights the adaptive capacity of anurans, particularly in response to heat and aridity, and explains why they are one of the continent's most visible faunas. [Anuran; adaptive radiation; Gondwana; phylogenetics].

Speciation across biomes: Rapid diversification with reproductive isolation in the Australian delicate mice

Phylogeographic studies of continental clades, especially when combined with palaeoclimate modelling, provide powerful insight into how environment drives speciation across climatic contexts. Australia, a continent characterized by disparate modern biomes and dynamic climate change, provides diverse opportunity to reconstruct the impact of past and present environments on diversification. Here, we use genomic-scale data (1310 exons and whole mitogenomes from 111 samples) to investigate Pleistocene diversification, cryptic diversity, and secondary contact in the Australian delicate mice (Hydromyini: Pseudomys), a recent radiation spanning almost all Australian environments. Across northern Australia, we find no evidence for introgression between cryptic lineages within Pseudomys delicatulus sensu lato, with palaeoclimate models supporting contraction and expansion of suitable habitat since the last glacial maximum. Despite multiple contact zones, we also find little evidence of introgression at a continental scale, with the exception of a potential hybrid zone in the mesic biome. In the arid zone, combined insights from genetic data and palaeomodels support a recent expansion in the arid specialist P. hermannsburgensis and contraction in the semi-arid P. bolami. In the face of repeated secondary contact, differences in sperm morphology and chromosomal rearrangements are potential mechanisms that maintain species boundaries in these recently diverged species. Additionally, we describe the western delicate mouse as a new species and recommend taxonomic reinstatement of the eastern delicate mouse. Overall, we show that speciation in an evolutionarily young and widespread clade has been driven by environmental change, and potentially maintained by divergence in reproductive morphology and chromosome rearrangements.

OrthoMaM v12: a database of curated single-copy ortholog alignments and trees to study mammalian evolutionary genomics

To date, the databases built to gather information on gene orthology do not provide end-users with descriptors of the molecular evolution information and phylogenetic pattern of these orthologues. In this context, we developed OrthoMaM, a database of ORTHOlogous MAmmalian Markers describing the evolutionary dynamics of coding sequences in mammalian genomes. OrthoMaM version 12 includes 15,868 alignments of orthologous coding sequences (CDS) from the 190 complete mammalian genomes currently available. All annotations and 1-to-1 orthology assignments are based on NCBI. Orthologous CDS can be mined for potential informative markers at the different taxonomic levels of the mammalian tree. To this end, several evolutionary descriptors of DNA sequences are provided for querying purposes (e.g. base composition and relative substitution rate). The graphical web interface allows the user to easily browse and sort the results of combined queries. The corresponding multiple sequence alignments and ML trees, inferred using state-of-the art approaches, are available for download both at the nucleotide and amino acid levels. OrthoMaM v12 can be used by researchers interested either in reconstructing the phylogenetic relationships of mammalian taxa or in understanding the evolutionary dynamics of coding sequences in their genomes. OrthoMaM is available for browsing, querying and complete or filtered download at https://orthomam.mbb.cnrs.fr/.

Phylogenomics reveals extreme gene tree discordance in a lineage of dominant trees: hybridization, introgression, and incomplete lineage sorting blur deep evolutionary relationships despite clear species groupings in Eucalyptus subgenus Eudesmia

Eucalypts are a large and ecologically important group of plants on the Australian continent, and understanding their evolution is important in understanding evolution of the unique Australian flora. Previous phylogenies using plastome DNA, nuclear-ribosomal DNA, or random genome-wide SNPs, have been confounded by limited genetic sampling or by idiosyncratic biological features of the eucalypts, including widespread plastome introgression. Here we present phylogenetic analyses of Eucalyptus subgenus Eudesmia (22 species from western, northern, central and eastern Australia), in the first study to apply a target-capture sequencing approach using custom, eucalypt-specific baits (of 568 genes) to a lineage of Eucalyptus. Multiple accessions of all species were included, and target-capture data were supplemented by separate analyses of plastome genes (average of 63 genes per sample). Analyses revealed a complex evolutionary history likely shaped by incomplete lineage sorting and hybridization. Gene tree discordance generally increased with phylogenetic depth. Species, or groups of species, toward the tips of the tree are mostly supported, and three major clades are identified, but the branching order of these clades cannot be confirmed with confidence. Multiple approaches to filtering the nuclear dataset, by removing genes or samples, could not reduce gene tree conflict or resolve these relationships. Despite inherent complexities in eucalypt evolution, the custom bait kit devised for this research will be a powerful tool for investigating the evolutionary history of eucalypts more broadly.

Genetic mixing in conservation translocations increases diversity of a keystone threatened species, Bettongia lesueur

Translocation programmes are increasingly being informed by genetic data to monitor and enhance conservation outcomes for both natural and established populations. These data provide a window into contemporary patterns of genetic diversity, structure and relatedness that can guide managers in how to best source animals for their translocation programmes. The inclusion of historical samples, where possible, strengthens monitoring by allowing assessment of changes in genetic diversity over time and by providing a benchmark for future improvements in diversity via management practices. Here, we used reduced representation sequencing (ddRADseq) data to report on the current genetic health of three remnant and seven translocated boodie (Bettongia lesueur) populations, now extinct on the Australian mainland. In addition, we used exon capture data from seven historical mainland specimens and a subset of contemporary samples to compare pre-decline and current diversity. Both data sets showed the significant impact of population founder source (whether multiple or single) on the genetic diversity of translocated populations. Populations founded by animals from multiple sources showed significantly higher genetic diversity than the natural remnant and single-source translocation populations, and we show that by mixing the most divergent populations, exon capture heterozygosity was restored to levels close to that observed in pre-decline mainland samples. Relatedness estimates were surprisingly low across all contemporary populations and there was limited evidence of inbreeding. Our results show that a strategy of genetic mixing has led to successful conservation outcomes for the species in terms of increasing genetic diversity and provides strong rationale for mixing as a management strategy.

Molecular evolution of male reproduction across species with highly divergent sperm morphology in diverse murine rodents

Sperm competition can drive rapid evolution of male reproductive traits, but it remains unclear how variation in sperm competition intensity shapes phenotypic and molecular diversity across clades. Old World mice and rats (subfamily Murinae) comprise a rapid radiation and exhibit incredible diversity in sperm morphology and production. We combined phenotype and sequence data to model the evolution of reproductive traits and genes across 78 murine species. We identified several shifts towards smaller relative testes mass, a trait reflective of reduced sperm competition. Several sperm traits were associated with relative testes mass, suggesting that mating system evolution likely selects for convergent traits related to sperm competitive ability. Molecular evolutionary rates of spermatogenesis proteins also correlated with relative testes mass, but in an unexpected direction. We predicted that sperm competition would result in rapid divergence among species with large relative testes mass, but instead found that many spermatogenesis genes evolve more rapidly in species with smaller relative testes mass due to relaxed purifying selection. While some reproductive genes evolved under positive selection, relaxed selection played a greater role underlying rapid evolution in small testes species. Our work demonstrates that sexual selection can impose strong purifying selection shaping the evolution of male reproduction.

The genomic landscape, causes, and consequences of extensive phylogenomic discordance in Old World mice and rats

A species tree is a central concept in evolutionary biology whereby a single branching phylogeny reflects relationships among species. However, the phylogenies of different genomic regions often differ from the species tree. Although tree discordance is often widespread in phylogenomic studies, we still lack a clear understanding of how variation in phylogenetic patterns is shaped by genome biology or the extent to which discordance may compromise comparative studies. We characterized patterns of phylogenomic discordance across the murine rodents (Old World mice and rats) - a large and ecologically diverse group that gave rise to the mouse and rat model systems. Combining new linked-read genome assemblies for seven murine species with eleven published rodent genomes, we first used ultra-conserved elements (UCEs) to infer a robust species tree. We then used whole genomes to examine finer-scale patterns of discordance and found that phylogenies built from proximate chromosomal regions had similar phylogenies. However, there was no relationship between tree similarity and local recombination rates in house mice, suggesting that genetic linkage influences phylogenetic patterns over deeper timescales. This signal may be independent of contemporary recombination landscapes. We also detected a strong influence of linked selection whereby purifying selection at UCEs led to less discordance, while genes experiencing positive selection showed more discordant and variable phylogenetic signals. Finally, we show that assuming a single species tree can result in high error rates when testing for positive selection under different models. Collectively, our results highlight the complex relationship between phylogenetic inference and genome biology and underscore how failure to account for this complexity can mislead comparative genomic studies.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Genomic data provide a robust phylogeny backbone for Rhodiola L. (Crassulaceae) and reveal extensive reticulate evolution during its rapid radiation

The Tibetan Plateau and adjacent mountain regions (TP; including the Tibetan Plateau, Himalaya, Hengduan Mountains and Mountains of Central Asia) harbor great biodiversity, some lineages on which may have undergone rapid radiations. However, only a few studies have investigated the evolutionary pattern of such diversification in depth using genomic data. In this study, we reconstructed a robust phylogeny backbone of Rhodiola, a lineage that may have undergone rapid radiation in the TP, using Genotyping-by-sequencing data, and conducted a series of gene flow and diversification analyses. The concatenation and coalescent-based methods yield similar tree topologies, and five well-supported clades were revealed. Potential gene flow and introgression events were detected, both between species from different major clades and closely related species, suggesting pervasive hybridization and introgression. An initial rapid and later slowdown of the diversification rate was revealed, indicating niche filling. Molecular dating and correlation analyses showed that the uplift of TP and global cooling in the mid-Miocene might have played an important role in promoting the rapid radiation of Rhodiola. Our work demonstrates that gene flow and introgression might be an important contributor to rapid radiation possibly by quickly reassembling old genetic variation into new combinations.

The Phylogenomic and Biogeographic History of the Gartersnakes, Watersnakes, and Allies (Natricidae: Thamnophiini)

North American Thamnophiini (gartersnakes, watersnakes, brownsnakes, and swampsnakes) are an ecologically and phenotypically diverse temperate clade of snakes representing 61 species across 10 genera. In this study, we estimate phylogenetic trees using ∼3,700 ultraconserved elements (UCEs) for 76 specimens representing 75% of all Thamnophiini species. We infer phylogenies using multispecies coalescent methods and time calibrate them using the fossil record. We also conducted ancestral area estimation to identify how major biogeographic boundaries in North America affect broadscale diversification in the group. While most nodes exhibited strong statistical support, analysis of concordant data across gene trees reveals substantial heterogeneity. Ancestral area estimation demonstrated that the genus Thamnophis was the only taxon in this subfamily to cross the Western Continental Divide, even as other taxa dispersed southward toward the tropics. Additionally, levels of gene tree discordance are overall higher in transition zones between bioregions, including the Rocky Mountains. Therefore, the Western Continental Divide may be a significant transition zone structuring the diversification of Thamnophiini during the Neogene and Pleistocene. Here we show that despite high levels of discordance across gene trees, we were able to infer a highly resolved and well-supported phylogeny for Thamnophiini, which allows us to understand broadscale patterns of diversity and biogeography.

Climate-influenced boreotropical survival and rampant introgressions explain the thriving of New World grapes in the north temperate zone

The north temperate region was characterized by a warm climate and a rich thermophilic flora before the Eocene, but early diversifications of the temperate biome under global climate change and biome shift remain uncertain. Moreover, it is becoming clear that hybridization/introgression is an important driving force of speciation in plant diversity. Here, we applied analyses from biogeography and phylogenetic networks to account for both introgression and incomplete lineage sorting based on genomic data from the New World Vitis, a charismatic component of the temperate North American flora with known and suspected gene flow among species. Biogeographic inference and fossil evidence suggest that the grapes were widely distributed from North America to Europe during the Paleocene to the Eocene, followed by widespread extinction and survival of relicts in the tropical New World. During the climate warming in the early Miocene, a Vitis ancestor migrated northward from the refugia with subsequent diversification in the North American region. We found strong evidence for widespread incongruence and reticulate evolution among nuclear genes within both recent and ancient lineages of the New World Vitis. Furthermore, the organellar genomes showed strong conflicts with the inferred species tree from the nuclear genomes. Our phylogenomic analyses provided an important assessment of the wide occurrence of reticulate introgression in the New World Vitis, which potentially represents one of the most important mechanisms for the diversification of Vitis species in temperate North America and even the entire temperate Northern Hemisphere. The scenario we report here may be a common model of temperate diversification of flowering plants adapted to the global climate cooling and fluctuation in the Neogene.

Population genomics and conservation management of the threatened black-footed tree-rat (Mesembriomys gouldii) in northern Australia

Genomic diversity is a fundamental component of Earth's total biodiversity, and requires explicit consideration in efforts to conserve biodiversity. To conserve genomic diversity, it is necessary to measure its spatial distribution, and quantify the contribution that any intraspecific evolutionary lineages make to overall genomic diversity. Here, we describe the range-wide population genomic structure of a threatened Australian rodent, the black-footed tree-rat (Mesembriomys gouldii), aiming to provide insight into the timing and extent of population declines across a large region with a dearth of long-term monitoring data. By estimating recent trajectories in effective population sizes at four localities, we confirm widespread population decline across the species' range, but find that the population in the peri-urban area of the Darwin region has been more stable. Based on current sampling, the Melville Island population made the greatest contribution to overall allelic richness of the species, and the prioritisation analysis suggested that conservation of the Darwin and Cobourg Peninsula populations would be the most cost-effective scenario to retain more than 90% of all alleles. Our results broadly confirm current sub-specific taxonomy, and provide crucial data on the spatial distribution of genomic diversity to help prioritise limited conservation resources. Along with additional sampling and genomic analysis from the far eastern and western edges of the black-footed tree-rat distribution, we suggest a range of conservation and research priorities that could help improve black-footed tree-rat population trajectories at large and fine spatial scales, including the retention and expansion of structurally complex habitat patches.

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.

Evolutionary models demonstrate rapid and adaptive diversification of Australo-Papuan pythons

Lineages may diversify when they encounter available ecological niches. Adaptive divergence by ecological opportunity often appears to follow the invasion of a new environment with open ecological space. This evolutionary process is hypothesized to explain the explosive diversification of numerous Australian vertebrate groups following the collision of the Eurasian and Australian plates 25 Mya. One of these groups is the pythons, which demonstrate their greatest phenotypic and ecological diversity in Australo-Papua (Australia and New Guinea). Here, using an updated and near complete time-calibrated phylogenomic hypothesis of the group, we show that following invasion of this region, pythons experienced a sudden burst of speciation rates coupled with multiple instances of accelerated phenotypic evolution in head and body shape and body size. These results are consistent with adaptive radiation theory with an initial rapid niche-filling phase and later slow-down approaching niche saturation. We discuss these findings in the context of other Australo-Papuan adaptive radiations and the importance of incorporating adaptive diversification systems that are not extraordinarily species-rich but ecomorphologically diverse to understand how biodiversity is generated.

Evolution: The evolutionary rat race in New Guinea and Australia

A new molecular phylogeny of a remarkable radiation of New Guinean and Australian rodents indicates multiple transitions between biomes and biogeographical regions within the group, and suggests that a key role was played by the geological history of New Guinea.

New Guinea uplift opens ecological opportunity across a continent

Sahul unites the world's largest and highest tropical island and the oldest and most arid continent on the backdrop of dynamic environmental conditions. Massive geological uplift in New Guinea is predicted to have acted as a species pump from the late Miocene onward, but the impact of this process on biogeography and diversification remains untested across Sahul as a whole. To address this, we reconstruct the assembly of a recent and diverse radiation of rodents (Murinae: Hydromyini) spanning New Guinea, Australia, and oceanic islands. Using phylogenomic data from 270 specimens, including many recently extinct and highly elusive species, we find that the orogeny and expansion of New Guinea opened ecological opportunity and triggered diversification across a continent. After a single over-water colonization from Asia ca. 8.5 Ma, ancestral Hydromyini were restricted to the tropical rainforest of proto-New Guinea for 3.5 million years. Following a shift in diversification coincident with the orogeny of New Guinea ca. 5 Ma and subsequent colonization of Australia, transitions between geographic regions (n = 24) and biomes (n = 34) become frequent. Recurrent over-water colonization between mainland and islands demonstrate how islands can play a substantial role in the assembly of continental fauna. Our results are consistent with a model of increased ecological opportunity across Sahul following major geological uplift in New Guinea ca. 5 Ma, with sustained diversification facilitated by over-water colonization from the Pleistocene to present. We show how geological processes, biome transitions, and over-water colonization collectively drove the diversification of an expansive continental radiation.

The contribution of incomplete lineage sorting and introgression to the evolutionary history of the fast-evolving genus Ctenomys (Rodentia, Ctenomyidae)

Incomplete lineage sorting and introgression have been increasingly recognized as important processes involved in biological differentiation. Both incomplete lineage sorting and introgression result in incongruences between gene trees and species trees, consequently causing difficulties in phylogenetic reconstruction. This is particularly the case for rapid radiations, as short internodal distances and incomplete reproductive isolation increase the likelihood of both ILS and introgression. Estimation of the relative frequency of these processes requires assessments across many genomic regions. We use transcriptomics to test for introgression and estimate the frequency of incomplete lineage sorting in a set of three closely related and geographically adjacent South American tuco-tucos species (Ctenomys), a genus comprising 64 species resulting from recent, rapid radiation. After cleaning and filtering, 5764 orthologous genes strongly support paraphyly of C. pearsoni relative to C. brasiliensis (putatively represented by the population of Villa Serrana). In line with earlier phylogenetic work, the C. pearsoni - C. brasiliensis pair is closely related to C. torquatus, whereas C. rionegrensis is more distantly related to these three nominal species. Classical Patterson's D-statistic shows significant signals of introgression from C. torquatus into C. brasiliensis. However, a 5-taxon test shows no significant results. Incomplete lineage sorting was estimated to have involved about 9% of the loci, suggesting it represents an important process in the incipient diversification of tuco-tucos.

Sequence Capture From Historical Museum Specimens: Maximizing Value for Population and Phylogenomic Studies

The application of high-throughput, short-read sequencing to degraded DNA has greatly increased the feasibility of generating genomic data from historical museum specimens. While many published studies report successful sequencing results from historical specimens; in reality, success and quality of sequence data can be highly variable. To examine predictors of sequencing quality, and methodological approaches to improving data accuracy, we generated and analyzed genomic sequence data from 115 historically collected museum specimens up to 180 years old. Data span both population genomic and phylogenomic scales, including historically collected specimens from 34 specimens of four species of Australian rock-wallabies (genus Petrogale) and 92 samples from 79 specimens of Australo-Papuan murine rodents (subfamily Murinae). For historical rodent specimens, where the focus was sampling for phylogenomics, we found that regardless of specimen age, DNA sequence libraries prepared from toe pad or bone subsamples performed significantly better than those taken from the skin (in terms of proportion of reads on target, number of loci captured, and data accuracy). In total, 93% of DNA libraries from toe pad or bone subsamples resulted in reliable data for phylogenetic inference, compared to 63% of skin subsamples. For skin subsamples, proportion of reads on target weakly correlated with collection year. Then using population genomic data from rock-wallaby skins as a test case, we found substantial improvement in final data quality by mapping to a high-quality “closest sister” de novo assembly from fresh tissues, compared to mapping to a sample-specific historical de novo assembly. Choice of mapping approach also affected final estimates of the number of segregating sites and Watterson's θ, both important parameters for population genomic inference. The incorporation of accurate and reliable sequence data from historical specimens has important outcomes for evolutionary studies at both population and phylogenomic scales. By assessing the outcomes of different approaches to specimen subsampling, library preparation and bioinformatic processing, our results provide a framework for increasing sequencing success for irreplaceable historical specimens.

Synthesis of Nuclear and Chloroplast Data Combined With Network Analyses Supports the Polyploid Origin of the Apple Tribe and the Hybrid Origin of the Maleae-Gillenieae Clade

Plant biologists have debated the evolutionary origin of the apple tribe (Maleae; Rosaceae) for over a century. The "wide-hybridization hypothesis" posits that the pome-bearing members of Maleae (base chromosome number x = 17) resulted from a hybridization and/or allopolyploid event between progenitors of other tribes in the subfamily Amygdaloideae with x = 8 and x = 9, respectively. An alternative "spiraeoid hypothesis" proposed that the x = 17 of Maleae arose via the genome doubling of x = 9 ancestors to x = 18, and subsequent aneuploidy resulting in x = 17. We use publicly available genomic data-448 nuclear genes and complete plastomes-from 27 species representing all major tribes within the Amygdaloideae to investigate evolutionary relationships within the subfamily containing the apple tribe. Specifically, we use network analyses and multi-labeled trees to test the competing wide-hybridization and spiraeoid hypotheses. Hybridization occurred between an ancestor of the tribe Spiraeeae (x = 9) and an ancestor of the clade Sorbarieae (x = 9) + Exochordeae (x = 8) + Kerrieae (x = 9), giving rise to the clade Gillenieae (x = 9) + Maleae (x = 17). The ancestor of the Maleae + Gillenieae arose via hybridization between distantly related tribes in the Amygdaloideae (i.e., supporting the wide hybridization hypothesis). However, some evidence supports an aspect of the spiraeoid hypothesis-the ancestors involved in the hybridization event were likely both x = 9, so genome doubling was followed by aneuploidy to result in x = 17 observed in Maleae. By synthesizing existing genomic data with novel analyses, we resolve the nearly century-old mystery regarding the origin of the apple tribe. Our results also indicate that nuclear gene tree-species tree conflict and/or cytonuclear conflict are pervasive at several other nodes in subfamily Amygdaloideae of Rosaceae.

The dynamics of introgression across an avian radiation

Hybridization and resulting introgression can play both a destructive and a creative role in the evolution of diversity. Thus, characterizing when and where introgression is most likely to occur can help us understand the causes of diversification dynamics. Here, we examine the prevalence of and variation in introgression using phylogenomic data from a large (1300+ species), geographically widespread avian group, the suboscine birds. We first examine patterns of gene tree discordance across the geographic distribution of the entire clade. We then evaluate the signal of introgression in a subset of 206 species triads using Patterson's D‐statistic and test for associations between introgression signal and evolutionary, geographic, and environmental variables. We find that gene tree discordance varies across lineages and geographic regions. The signal of introgression is highest in cases where species occur in close geographic proximity and in regions with more dynamic climates since the Pleistocene. Our results highlight the potential of phylogenomic datasets for examining broad patterns of hybridization and suggest that the degree of introgression between diverging lineages might be predictable based on the setting in which they occur.

Whole-Genome Phylogenetic Reconstruction as a Powerful Tool to Reveal Homoplasy and Ancient Rapid Radiation in Waterflea Evolution

Although phylogeny estimation is notoriously difficult in radiations that occurred several hundred million years ago, phylogenomic approaches offer new ways to examine relationships among ancient lineages and evaluate hypotheses that are key to evolutionary biology. Here, we reconstruct the deep-rooted relationships of one of the oldest living arthropod clades, the branchiopod crustaceans, using a kaleidoscopic approach. We use concatenation and coalescent tree-building methods to analyze a large multigene data set at the nucleotide and amino acid level and examine gene tree versus species tree discordance. We unequivocally resolve long-debated relationships among extant orders of the Cladocera, the waterfleas, an ecologically relevant zooplankton group in global aquatic and marine ecosystems that is famous for its model systems in ecology and evolution. To build the data set, we assembled eight de novo genomes of key taxa including representatives of all extant cladoceran orders and suborders. Our phylogenetic analysis focused on a BUSCO-based set of 823 conserved single-copy orthologs shared among 23 representative taxa spanning all living branchiopod orders, including 11 cladoceran families. Our analysis supports the monophyly of the Cladocera and reveals remarkable homoplasy in their body plans. We found large phylogenetic distances between lineages with similar ecological specializations, indicating independent evolution in major body plans, such as in the pelagic predatory orders Haplopoda and Onychopoda (the “Gymnomera”). In addition, we assessed rapid cladogenesis by estimating relative timings of divergence in major lineages using reliable fossil-calibrated priors on eight nodes in the branchiopod tree, suggesting a Paleozoic origin around 325 Ma for the cladoceran ancestor and an ancient rapid radiation around 252 Ma at the Perm/Triassic boundary. These findings raise new questions about the roles of homoplasy and rapid radiation in the diversification of the cladocerans and help examine trait evolution from a genomic perspective in a functionally well understood, ancient arthropod group. [Cladocera; Daphnia; evolution; homoplasy; molecular clock; phylogenomics; systematics; waterfleas.]

Phylogenomic and Macroevolutionary Evidence for an Explosive Radiation of a Plant Genus in the Miocene

Mountain systems harbor a substantial fraction of global biodiversity and, thus, provide excellent opportunities to study rapid diversification and to understand the historical processes underlying the assembly of biodiversity hotspots. The rich biodiversity in mountains is widely regarded as having arisen under the influence of geological and climatic processes as well as the complex interactions among them. However, the relative contribution of geology and climate in driving species radiation is seldom explored. Here, we studied the evolutionary radiation of Oreocharis (Gesneriaceae), which has diversified extensively throughout East Asia, especially within the Hengduan Mountains (HDM), using transcriptomic data and a time calibrated phylogeny for 88% (111/126) of all species of the genus. In particular, we applied phylogenetic reconstructions to evaluate the extent of incomplete lineage sorting accompanying the early and rapid radiation in the genus. We then fit macroevolutionary models to explore its spatial and diversification dynamics in Oreocharis and applied explicit birth–death models to investigate the effects of past environmental changes on its diversification. Evidence from 574 orthologous loci suggest that Oreocharis underwent an impressive early burst of speciation starting ca. 12 Ma in the Miocene, followed by a drastic decline in speciation toward the present. Although we found no evidence for a shift in diversification rate across the phylogeny of Oreocharis, we showed a difference in diversification dynamics between the HDM and non-HDM lineages, with higher diversification rates in the HDM. The diversification dynamic of Oreocharis is most likely positively associated with temperature-dependent speciation and dependency on the Asian monsoons. We suggest that the warm and humid climate of the mid-Miocene was probably the primary driver of the rapid diversification in Oreocharis, while mountain building of the HDM might have indirectly affected species diversification of the HDM lineage. This study highlights the importance of past climatic changes, combined with mountain building, in creating strong environmental heterogeneity and driving diversification of mountain plants, and suggests that the biodiversity in the HDM cannot directly be attributed to mountain uplift, contrary to many recent speculations.[East Asian monsoons; environmental heterogeneity; Hengduan Mountains; incomplete lineage sorting; Oreocharis; past climate change; rapid diversification; transcriptome.]

Museum genomics reveals the rapid decline and extinction of Australian rodents since European settlement

Australia has the highest historically recorded rate of mammalian extinction in the world, with 34 terrestrial species declared extinct since European colonization in 1788. Among Australian mammals, rodents have been the most severely affected by these recent extinctions; however, given a sparse historical record, the scale and timing of their decline remain unresolved. Using museum specimens up to 184 y old, we generate genomic-scale data from across the entire assemblage of Australian hydromyine rodents (i.e., eight extinct species and their 42 living relatives). We reconstruct a phylogenomic tree for these species spanning ∼5.2 million years, revealing a cumulative total of 10 million years (>10%) of unique evolutionary history lost to extinction within the past ∼150 y. We find no evidence for reduced genetic diversity in extinct species just prior to or during decline, indicating that their extinction was extremely rapid. This suggests that populations of extinct Australian rodents were large prior to European colonization, and that genetic diversity does not necessarily protect species from catastrophic extinction. In addition, comparative analyses suggest that body size and biome interact to predict extinction and decline, with larger species more likely to go extinct. Finally, we taxonomically resurrect a species from extinction, Gould's mouse (Pseudomys gouldii Waterhouse, 1839), which survives as an island population in Shark Bay, Western Australia (currently classified as Pseudomys fieldi Waite, 1896). With unprecedented sampling across a radiation of extinct and living species, we unlock a previously inaccessible historical perspective on extinction in Australia. Our results highlight the capacity of collections-based research to inform conservation and management of persisting species.

A genomic perspective on an old question: Salmo trouts or Salmo trutta (Teleostei: Salmonidae)?

There are particular challenges in defining the taxonomic status of recently radiated groups due to the low level of phylogenetic signal. Members of the Salmo trutta species-complex, which mostly evolved during and following the Pleistocene, show high morphological and ecological diversity that, along with their very wide geographic distribution, have led to morphological description of 47 extant nominal species. However, many of these species have not been supported by previous phylogenetic studies, which could be partly due to lack of significant genetic differences among them, the limited resolution offered by molecular methods previously used, as well as the often local scale of these studies. The development of next-generation sequencing (NGS) and related analytical tools have enhanced our ability to address such challenging questions. In this study, Genotyping-by-Sequencing (GBS) of 15,169 filtered SNPs and mitochondrial DNA (mtDNA) D-loop sequences were combined to assess the phylogenetic relationships among 166 brown trouts representing 21 described species and three undescribed groups collected from 84 localities throughout their natural distribution in Europe, west Asia, and North Africa. The data were analysed using different clustering algorithms (admixture analysis and discriminant analysis of principal components-DAPC), a Bayes Factor Delimitation (BFD) test, species tree reconstruction, gene flow tests (three- and four-population tests), and Rogue taxa identification tests. Genomic contributions of the Atlantic lineage brown trout were found in all major sea basins excluding the North African and Aral Sea basins, suggesting introgressive hybridization of native brown trouts driven by stocking using strains of the Atlantic lineage. After removing the phylogenetic noise caused by the Atlantic brown trout, admixture clusters and DAPC clustering based on GBS data, respectively, resolved 11 and 13 clusters among the previously described brown trout species, which were also supported by BFD test results. Our results suggest that natural hybridization between different brown trout lineages has probably played an important role in the origin of several of the putative species, including S. marmoratus, S. carpio, S. farioides, S. pellegrini, S. caspius (in the Kura River drainage) and Salmo sp. in the Danube River basin. Overall, our results support a multi-species taxonomy for brown trouts. They also resolve some species in the Adriatic-Mediterranean and Black Sea drainages as members of very closely related genomic clusters that may need taxonomic revision. However, any final conclusions pertaining to the taxonomy of the brown trout complex should be based on an integrative approach combining genomic, morphological, and ecological data. To avoid challenges in taxonomy and conservation of species complexes like brown trouts, it is suggested to describe species based on genomic clusters of populations instead of describing species based only on morphologically differentiated single type populations.

Molecular Evolution of Ecological Specialisation: Genomic Insights from the Diversification of Murine Rodents

Adaptive radiations are characterized by the diversification and ecological differentiation of species, and replicated cases of this process provide natural experiments for understanding the repeatability and pace of molecular evolution. During adaptive radiation, genes related to ecological specialization may be subject to recurrent positive directional selection. However, it is not clear to what extent patterns of lineage-specific ecological specialization (including phenotypic convergence) are correlated with shared signatures of molecular evolution. To test this, we sequenced whole exomes from a phylogenetically dispersed sample of 38 murine rodent species, a group characterized by multiple, nested adaptive radiations comprising extensive ecological and phenotypic diversity. We found that genes associated with immunity, reproduction, diet, digestion, and taste have been subject to pervasive positive selection during the diversification of murine rodents. We also found a significant correlation between genome-wide positive selection and dietary specialization, with a higher proportion of positively selected codon sites in derived dietary forms (i.e., carnivores and herbivores) than in ancestral forms (i.e., omnivores). Despite striking convergent evolution of skull morphology and dentition in two distantly related worm-eating specialists, we did not detect more genes with shared signatures of positive or relaxed selection than in a nonconvergent species comparison. Although a small number of the genes we detected can be incidentally linked to craniofacial morphology or diet, protein-coding regions are unlikely to be the primary genetic basis of this complex convergent phenotype. Our results suggest a link between positive selection and derived ecological phenotypes, and highlight specific genes and general functional categories that may have played an integral role in the extensive and rapid diversification of murine rodents.

Population genomics and conservation management of a declining tropical rodent

Conservation management is improved by incorporating information about the spatial distribution of population genetic diversity into planning strategies. Northern Australia is the location of some of the world's most severe ongoing declines of endemic mammal species, yet we have little genetic information from this regional mammal assemblage to inform a genetic perspective on conservation assessment and planning. We used next-generation sequencing data from remnant populations of the threatened brush-tailed rabbit-rat (Conilurus penicillatus) to compare patterns of genomic diversity and differentiation across the landscape and investigate standardised hierarchical genomic diversity metrics to better understand brush-tailed rabbit-rat population genomic structure. We found strong population structuring, with high levels of differentiation between populations (F_ST = 0.21-0.78). Two distinct genomic lineages between the Tiwi Islands and mainland are also present. Prioritisation analysis showed that one population in both lineages would need to be conserved to retain at least ~80% of alleles for the species. Analysis of standardised genomic diversity metrics showed that approximately half of the total diversity occurs among lineages (δ = 0.091 from grand total γ = 0.184). We suggest that a focus on conserving remnant island populations may not be appropriate for the preservation of species-level genomic diversity and adaptive potential, as these populations represent a small component of the total diversity and a narrow subset of the environmental conditions in which the species occurs. We also highlight the importance of considering both genomic and ecological differentiation between source and receiving populations when considering translocations for conservation purposes.

Testing the utility of alternative metrics of branch support to address the ancient evolutionary radiation of tunas, stromateoids, and allies (Teleostei: Pelagiaria)

The use of high-throughput sequencing technologies to produce genome-scale datasets was expected to settle some long-standing controversies across the Tree of Life, particularly in areas where short branches occur at deep timescales. Instead, these datasets have often yielded many well-supported but conflicting topologies, and highly variable gene-tree distributions. A variety of branch-support metrics beyond the nonparametric bootstrap are now available to assess how robust a phylogenetic hypothesis may be, as well as new methods to quantify gene-tree discordance. We applied multiple branch support metrics to an ancient group of marine fishes (Teleostei: Pelagiaria) whose interfamilial relationships have proven difficult to resolve due to a rapid accumulation of lineages very early in its history. We analyzed hundreds of loci including published UCE data and newly generated exonic data along with their flanking regions to represent all 16 extant families for more than 150 out of 284 valid species in the group. Branch support was lower for interfamilial relationships (except the SH-like aLRT and aBayes methods) regardless of the type of marker used. Several nodes that were highly supported with bootstrap had very low site and gene-tree concordance, revealing underlying conflict. Despite this conflict, we were able to identify four consistent interfamilial clades, each comprised of two or three families. Combining exons with their flanking regions also produced increased branch lengths in the deep branches of the pelagiarian tree. Our results demonstrate the limitations of employing current metrics of branch support and species-tree estimation when assessing the confidence of ancient evolutionary radiations and emphasize the necessity to embrace alternative measurements to explore phylogenetic uncertainty and discordance in phylogenomic datasets.

Diversification of the mygalomorph spider genus Aname (Araneae: Anamidae) across the Australian arid zone: Tracing the evolution and biogeography of a continent-wide radiation

The assembly of the Australian arid zone biota has long fascinated biogeographers. Covering over two-thirds of the continent, Australia's vast arid zone biome is home to a distinctive fauna and flora, including numerous lineages which have diversified since the Eocene. Tracing the origins and speciation history of these arid zone taxa has been an ongoing endeavour since the advent of molecular phylogenetics, and an increasing number of studies on invertebrate animals are beginning to complement a rich history of research on vertebrate and plant taxa. In this study, we apply continent-wide genetic sampling and one of the largest phylogenetic data matrices yet assembled for a genus of Australian spiders, to reconstruct the phylogeny and biogeographic history of the open-holed trapdoor spider genus Aname L. Koch, 1873. This highly diverse lineage of Australian mygalomorph spiders has a distribution covering the majority of Australia west of the Great Dividing Range, but apparently excluding the high rainfall zones of eastern Australia and Tasmania. Original and legacy sequences were obtained for three mtDNA and four nuDNA markers from 174 taxa in seven genera, including 150 Aname specimen terminals belonging to 102 species-level operational taxonomic units, sampled from 32 bioregions across Australia. Reconstruction of the phylogeny and biogeographic history of Aname revealed three radiations (Tropical, Temperate-Eastern and Continental), which could be further broken into eight major inclusive clades. Ancestral area reconstruction revealed the Pilbara, Monsoon Tropics and Mid-West to be important ancestral areas for the genus Aname and its closest relatives, with the origin of Aname itself inferred in the Pilbara bioregion. From these origins in the arid north-west of Australia, our study found evidence for a series of subsequent biome transitions in separate lineages, with at least eight tertiary incursions back into the arid zone from more mesic tropical, temperate or eastern biomes, and only two major clades which experienced widespread (primary) in situ diversification within the arid zone. Based on our phylogenetic results, and results from independent legacy divergence dating studies, we further reveal the importance of climate-driven biotic change in the Miocene and Pliocene in shaping the distribution and composition of the Australian arid zone biota, and the value of continent-wide studies in revealing potentially complex patterns of arid zone diversification in dispersal-limited invertebrate taxa.

Dissecting incongruence between concatenation- and quartet-based approaches in phylogenomic data

Topological conflict or incongruence is widespread in phylogenomic data. Concatenation- and coalescent-based approaches often result in incongruent topologies, but the causes of this conflict can be difficult to characterize. We examined incongruence stemming from conflict between likelihood-based signal (quantified by the difference in gene-wise log likelihood score or ΔGLS) and quartet-based topological signal (quantified by the difference in gene-wise quartet score or ΔGQS) for every gene in three phylogenomic studies in animals, fungi, and plants, which were chosen because their concatenation-based IQ-TREE (T1) and quartet-based ASTRAL (T2) phylogenies are known to produce eight conflicting internal branches (bipartitions). By comparing the types of phylogenetic signal for all genes in these three data matrices, we found that 30% - 36% of genes in each data matrix are inconsistent, that is, each of these genes has higher log likelihood score for T1 versus T2 (i.e., ΔGLS >0) whereas its T1 topology has lower quartet score than its T2 topology (i.e., ΔGQS <0) or vice versa. Comparison of inconsistent and consistent genes using a variety of metrics (e.g., evolutionary rate, gene tree topology, distribution of branch lengths, hidden paralogy, and gene tree discordance) showed that inconsistent genes are more likely to recover neither T1 nor T2 and have higher levels of gene tree discordance than consistent genes. Simulation analyses demonstrate that removal of inconsistent genes from datasets with low levels of incomplete lineage sorting (ILS) and low and medium levels of gene tree estimation error (GTEE) reduced incongruence and increased accuracy. In contrast, removal of inconsistent genes from datasets with medium and high ILS levels and high GTEE levels eliminated or extensively reduced incongruence, but the resulting congruent species phylogenies were not always topologically identical to the true species trees.

Target-capture phylogenomics provide insights on gene and species tree discordances in Old World treefrogs (Anura: Rhacophoridae)

Genome-scale data have greatly facilitated the resolution of recalcitrant nodes that Sanger-based datasets have been unable to resolve. However, phylogenomic studies continue to use traditional methods such as bootstrapping to estimate branch support; and high bootstrap values are still interpreted as providing strong support for the correct topology. Furthermore, relatively little attention has been given to assessing discordances between gene and species trees, and the underlying processes that produce phylogenetic conflict. We generated novel genomic datasets to characterize and determine the causes of discordance in Old World treefrogs (Family: Rhacophoridae)-a group that is fraught with conflicting and poorly supported topologies among major clades. Additionally, a suite of data filtering strategies and analytical methods were applied to assess their impact on phylogenetic inference. We showed that incomplete lineage sorting was detected at all nodes that exhibited high levels of discordance. Those nodes were also associated with extremely short internal branches. We also clearly demonstrate that bootstrap values do not reflect uncertainty or confidence for the correct topology and, hence, should not be used as a measure of branch support in phylogenomic datasets. Overall, we showed that phylogenetic discordances in Old World treefrogs resulted from incomplete lineage sorting and that species tree inference can be improved using a multi-faceted, total-evidence approach, which uses the most amount of data and considers results from different analytical methods and datasets.

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations. [Amaranthaceae; gene tree discordance; hybridization; incomplete lineage sorting; phylogenomics; species network; species tree; transcriptomics.]

Environment predicts repeated body size shifts in a recent radiation of Australian mammals

Closely related species that occur across steep environmental gradients often display clear body size differences, and examining this pattern is crucial to understanding how environmental variation shapes diversity. Australian endemic rodents in the Pseudomys Division (Muridae: Murinae) have repeatedly colonized the arid, monsoon, and mesic biomes over the last 5 million years. Using occurrence records, body mass data, and Bayesian phylogenetic models, we test whether body mass of 31 species in the Pseudomys Division can be predicted by their biome association. We also model the effect of eight environmental variables on body mass. Despite high phylogenetic signal in body mass evolution across the phylogeny, we find that mass predictably increases in the mesic biome and decreases in arid and monsoon biomes. As per Bergmann's rule, temperature is strongly correlated with body mass, as well as several other variables. Our results highlight two important findings. First, body size in Australian rodents has tracked with climate through the Pleistocene, likely due to several environmental variables rather than a single factor. Second, support for both Brownian motion and predictable change at different taxonomic levels in the Pseudomys Division phylogeny demonstrates how the level at which we test hypotheses can alter interpretation of evolutionary processes.