Quartet inference from SNP data under the coalescent model

Underground speciation: Unraveling the systematics and evolution of the highly diverse tuco-tucos (genus Ctenomys) with genomic data

Subterranean rodents of the genus Ctenomys (tuco-tucos) are endemic to South America and have experienced relatively recent radiation. There are about 67 recognized species that originated in approximately 1-2 MY. They stand out for their species richness, extraordinary chromosomal diversity, and wide range of habitat they occupy in the continent. Phylogenetic relationships among species of tuco-tucos have been challenging to resolve. Groups of closely-related species have been suggested, but their relationships must be resolved. This study estimates the phylogeny of the genus using massive sequencing, generating thousands of independent molecular markers obtained by RADseq, with a taxonomic sampling that includes 66% of the recognized species. The sequences obtained were mapped against the C. sociabilis genome, recovering up to 1,215 widely shared RAD loci with more than 19,000 polymorphic sites. Our new phylogenetic hypothesis corroborated the species groups previously proposed with cytochrome b gene sequences and provided a much greater resolution of the relationships among species groups. The frater group is sister to all other tuco-tucos, whereas some of the earlierliest proposals placed the sociabilis group as sister to all other tuco-tucos. Ctenomys leucodon, previously proposed as an independent lineage, is associated with the frater group with moderate statistical support. The magellanicus and mendocinus are sister groups in a major clade formed by the boliviensis, talarum, tucumanus, torquatus, and opimus groups. Ctenomys viperinus, included in the phylogeny for the first time, belongs to the tucumanus group. This multi-locus phylogenetic hypothesis provides insights into the historical biogeography of understanding this highly diverse genus.

Phylogenomics and biogeography of sawflies and woodwasps (Hymenoptera, Symphyta)

Phylogenomic approaches have recently helped elucidate various insect relationships, but large-scale comprehensive analyses on relationships within sawflies and woodwasps are still lacking. Here, we infer the relationships and long-term biogeographic history of these hymenopteran groups using a large dataset of 354 UCE loci collected from 385 species that represent all major lineages. Early Hymenoptera started diversifying during the Early Triassic ∼249 Ma and spread all over the ancient supercontinent Pangaea. We recovered Xyeloidea as a monophyletic sister group to other Hymenoptera and Pamphilioidea as sister to Unicalcarida. Within the diverse family Tenthredinidae, our taxonomically and geographically expanded taxon sampling highlights the non-monophyly of several traditionally defined subfamilies. In addition, the recent removal of Athalia and related genera from the Tenthredinidae into the separate family Athaliidae is supported. The deep historical biogeography of the group is characterised by independent dispersals and re-colonisations between the northern (Laurasia) and southern (Gondwana) palaeocontinents. The breakup of these landmasses led to ancient vicariance in several Gondwanan lineages, while interchange across the Northern Hemisphere has continued until the Recent. The little-studied African sawfly fauna is likewise a diverse mixture of groups with varying routes of colonization. Our results reveal interesting parallels in the evolution and biogeography of early hymenopterans and other ancient insect groups.

Genomic SNPs resolve the phylogeny of an ancient amphibian island radiation from the Seychelles

Unusually for oceanic islands, the granitic Seychelles host multiple lineages of endemic amphibians. This includes an ancient (likely ca. 60 million years) radiation of eight caecilian species, most of which occur on multiple islands.These caecilians have a complicated taxonomic history and their phylogenetic inter-species relationships have been difficult to resolve. Double-digest RAD sequencing (ddRADseq) has been applied extensively to phylogeography and increasingly to phylogenetics but its utility for resolving ancient divergences is less well established. To address this, we applied ddRADseq to generate a genome-wide SNP panel for phylogenomic analyses of the Seychelles caecilians, whose phylogeny has so far not been satisfactorily resolved with traditional DNA markers. Based on 129,154 SNPs, we resolved deep and shallow splits, with strong support. Our findings demonstrate the capability of genome-wide SNPs for evolutionary inference at multiple taxonomic levels and support the recently proposed synonymy of Grandisonia Taylor, 1968 with Hypogeophis Peters, 1879. We revealed three clades of Hypogeophis (large-, medium- and short-bodied) and identify a single origin of the diminutive, stocky-bodied and pointy-snouted phenotype.

Incomplete lineage sorting and hybridization underlie tree discordance in Petunia and related genera (Petunieae, Solanaceae)

Despite the overarching history of species divergence, phylogenetic studies often reveal distinct topologies across regions of the genome. The sources of these gene tree discordances are variable, but incomplete lineage sorting (ILS) and hybridization are among those with the most biological importance. Petunia serves as a classic system for studying hybridization in the wild. While field studies suggest that hybridization is frequent, the extent of reticulation within Petunia and its closely related genera has never been examined from a phylogenetic perspective. In this study, we used transcriptomic data from 11 Petunia, 16 Calibrachoa, and 10 Fabiana species to illuminate the relationships between these species and investigate whether hybridization played a significant role in the diversification of the clade. We inferred that gene tree discordance within genera is linked to hybridization events along with high levels of ILS due to their rapid diversification. Moreover, network analyses estimated deeper hybridization events between Petunia and Calibrachoa, genera that have different chromosome numbers. Although these genera cannot hybridize at the present time, ancestral hybridization could have played a role in their parallel radiations, as they share the same habitat and life history.

Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China

Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.

Cryptic diversity and phylogeographic patterns of Mediodactylus species in the Eastern Mediterranean region

Cryptic diversity poses a great obstacle in our attempts to assess the current biodiversity crisis and may hamper conservation efforts. The gekkonid genus Mediodactylus, a well-known case of hidden species and genetic diversity, has been taxonomically reclassified several times during the last decade. Focusing on the Mediterranean populations, a recent study within the M. kotschyi species complex using classic mtDNA/nuDNA markers suggested the existence of five distinct species, some being endemic and some possibly threatened, yet their relationships have not been fully resolved. Here, we generated genome-wide SNPs (using ddRADseq) and applied molecular species delimitation approaches and population genomic analyses to further disentangle these relationships. Τhe most extensive nuclear dataset, so far, encompassing 2,360 loci and ∼ 699,000 bp from across the genome of Mediodactylus gecko, enabled us to resolve previously obscure phylogenetic relationships among the five, recently elevated, Mediodactylus species and to support the hypothesis that the taxon includes several new, undescribed species. Population genomic analyses within each of the proposed species showed strong genetic structure and high levels of genetic differentiation among populations.

Evidence for gene flow and trait reversal during radiation of Mexican Goodeid fish

Understanding the phylogeographic history of a group and identifying the factors contributing to speciation is an important challenge in evolutionary biology. The Goodeinae are a group of live-bearing fishes endemic to Mexico. Here, we develop genomic resources for species within the Goodeinae and use phylogenomic approaches to characterise their evolutionary history. We sequenced, assembled and annotated the genomes of four Goodeinae species, including Ataeniobius toweri, the only matrotrophic live-bearing fish without a trophotaenia in the group. We estimated timings of species divergence and examined the extent and timing of introgression between the species to assess if this may have occurred during an early radiation, or in more recent episodes of secondary contact. We used branch-site models to detect genome-wide positive selection across Goodeinae, and we specifically asked whether this differs in A. toweri, where loss of placental viviparity has recently occurred. We found evidence of gene flow between geographically isolated species, suggesting vicariant speciation was supplemented by limited post-speciation gene flow, and gene flow may explain previous uncertainties about Goodeid phylogeny. Genes under positive selection in the group are likely to be associated with the switch to live-bearing. Overall, our studies suggest that both volcanism-driven vicariance and changes in reproductive mode influenced radiation in the Goodeinae.

Dragons in the tropics - Phylogeography and speciation in Diporiphora lizards and common geographic breaks in co-distributed taxa

Co-distributed taxa can respond both similarly or differently to the same climatic and geological events, resulting in a range of phylogeographic patterns across the region. Using a nested approach on a taxonomically diverse yet morphologically conservative group of agamid lizards, we first aimed to evaluate more precisely the extent of phylogeographic structuring within the genus. Then, focusing on four lineages within the more widespread species, we assessed the impact of biogeographic barriers on phylogeographic structuring and demographic history of species, comparing to patterns previously observed in co-distributed taxa. These species occur in the Australian Monsoonal Tropics, a vast tropical savanna system with high richness and endemism associated with environmental heterogeneity and past climate fluctuations. The employment of genomic data helped to determine the relationships between specific taxa that were previously difficult to place. We found a local influence of biogeographic and climatic breaks on population dynamics, analogous to other species. We detected high levels of population structure in the West Kimberley and Arnhem Plateau, which are already known for high endemism. However, we also highlighted unique lineages in areas that have been overlooked until recently, in the South Kimberley and West Top End. Climatic and geographical features in the Arnhem Plateau act as a soft barrier between populations in the east and west regions of the Top End. These observations reflect patterns observed for other vertebrates across this rich biome, indicating how climatic variation, species' ecology, and landscape features interact to shape regional diversity and endemism.

Phylogeny and historical biogeography of the swallow family (Hirundinidae) inferred from comparisons of thousands of UCE loci

Swallows (Hirundinidae) are a globally distributed family of passerine birds that exhibit remarkable similarity in body shape but tremendous variation in plumage, sociality, nesting behavior, and migratory strategies. As a result, swallow species have become models for empirical behavioral ecology and evolutionary studies, and variation across the Hirundinidae presents an excellent opportunity for comparative analyses of trait evolution. Exploiting this potential requires a comprehensive and well-resolved phylogenetic tree of the family. To address this need, we estimated swallow phylogeny using genetic data from thousands of ultraconserved element (UCE) loci sampled from nearly all recognized swallow species. Maximum likelihood, coalescent-based, and Bayesian approaches yielded a well-resolved phylogenetic tree to the generic level, with minor disagreement among inferences at the species level, which likely reflect ongoing population genetic processes. The UCE data were particularly useful in helping to resolve deep nodes, which previously confounded phylogenetic reconstruction efforts. Divergence time estimates from the improved swallow tree support a Miocene origin of the family, roughly 13 million years ago, with subsequent diversification of major groups in the late Miocene and Pliocene. Our estimates of historical biogeography support the hypothesis that swallows originated in the Afrotropics and have subsequently expanded across the globe, with major in situ diversification in Africa and a secondary major radiation following colonization of the Neotropics. Initial examination of nesting and sociality indicates that the origin of mud nesting - a relatively rare nest construction phenotype in birds - was a major innovation coincident with the origin of a clade giving rise to over 40% of extant swallow diversity. In contrast, transitions between social and solitary nesting appear less important for explaining patterns of diversification among swallows.

Gene Flow and Isolation in the Arid Nearctic Revealed by Genomic Analyses of Desert Spiny Lizards

The opposing forces of gene flow and isolation are two major processes shaping genetic diversity. Understanding how these vary across space and time is necessary to identify the environmental features that promote diversification. The detection of considerable geographic structure in taxa from the arid Nearctic has prompted research into the drivers of isolation in the region. Several geographic features have been proposed as barriers to gene flow, including the Colorado River, Western Continental Divide (WCD), and a hypothetical Mid-Peninsular Seaway in Baja California. However, recent studies suggest that the role of barriers in genetic differentiation may have been overestimated when compared to other mechanisms of divergence. In this study, we infer historical and spatial patterns of connectivity and isolation in Desert Spiny Lizards (Sceloporus magister) and Baja Spiny Lizards (Sceloporus zosteromus), which together form a species complex composed of parapatric lineages with wide distributions in arid western North America. Our analyses incorporate mitochondrial sequences, genomic-scale data, and past and present climatic data to evaluate the nature and strength of barriers to gene flow in the region. Our approach relies on estimates of migration under the multispecies coalescent to understand the history of lineage divergence in the face of gene flow. Results show that the S. magister complex is geographically structured, but we also detect instances of gene flow. The WCD is a strong barrier to gene flow, while the Colorado River is more permeable. Analyses yield conflicting results for the catalyst of differentiation of peninsular lineages in S. zosteromus. Our study shows how large-scale genomic data for thoroughly sampled species can shed new light on biogeography. Furthermore, our approach highlights the need for the combined analysis of multiple sources of evidence to adequately characterize the drivers of divergence.

Waves of Colonization and Gene Flow in a Great Speciator

Secondary contact between previously allopatric lineages offers a test of reproductive isolating mechanisms that may have accrued in isolation. Such instances of contact can produce stable hybrid zones-where reproductive isolation can further develop via reinforcement or phenotypic displacement-or result in the lineages merging. Ongoing secondary contact is most visible in continental systems, where steady input from parental taxa can occur readily. In oceanic island systems, however, secondary contact between closely related species of birds is relatively rare. When observed on sufficiently small islands, relative to population size, secondary contact likely represents a recent phenomenon. Here, we examine the dynamics of a group of birds whose apparent widespread hybridization influenced Ernst Mayr's foundational work on allopatric speciation: the whistlers of Fiji (Aves: Pachycephala). We demonstrate two clear instances of secondary contact within the Fijian archipelago, one resulting in a hybrid zone on a larger island, and the other resulting in a wholly admixed population on a smaller, adjacent island. We leveraged low genome-wide divergence in the hybrid zone to pinpoint a single genomic region associated with observed phenotypic differences. We use genomic data to present a new hypothesis that emphasizes rapid plumage evolution and post-divergence gene flow.

Phylogenomic analyses re-examine the evolution of reinforcement and hypothesized hybrid speciation in Phlox wildflowers

The tree of life is riddled with reticulate evolutionary histories, and some clades, such as the eastern standing Phlox, appear to be hotspots of hybridization. In this group, there are two cases of reinforcement and nine hypothesized hybrid species. Given their historical importance in our understanding of plant speciation, the relationships between these taxa and the role of hybridization in their diversification require genomic validation. Using phylogenomic analyses, we resolve the evolutionary relationships of the eastern standing Phlox and evaluate hypotheses about whether and how hybridization and gene flow played a role in their diversification. Our results provide novel resolution of the phylogenetic relationships in this group, including paraphyly across some taxa. We identify gene flow during one case of reinforcement and find genomic support for a hybrid lineage underlying one of the five hypothesized homoploid hybrid speciation events. Additionally, we estimate the ancestries of four allotetraploid hybrid species. Our results are consistent with hybridization contributing to diverse evolutionary outcomes within this group; although, not as extensively as previously hypothesized. This study demonstrates the importance of phylogenomics in evaluating hypothesized evolutionary histories of non-model systems and adds to the growing support of interspecific genetic exchange in the generation of biodiversity.

Differential expression analyses and detection of SNP loci associated with environmental variables: Are salinity and temperature factors involved in population differentiation and speciation in Odontesthes?

Environmental factors play a key role in individual adaptation to different local conditions. Because of this, studies about the physiological and genetic responses of individuals exposed to different natural environments offer clues about mechanisms involved in population differentiation, and as a subsequent result, speciation. Marine environments are especially suited to survey this kind of phenomena because they commonly harbor species adapted to different local conditions along a geographic continuum. Silversides belonging to Odontesthes are commonly distributed in tropical and temperate regions of South America and exhibit noticeable phenotypic plasticity, which allows them to adapt to contrasting environments. In this study, the genetic expression of O. argentinensis sampled along the Uruguayan Atlantic coast and estuarine adjacent areas was investigated. In addition, the correlation between individual genotypes and environmental variables was also analysed in O. argentinensis and O. bonariensis. Results obtained suggest a differential expression pattern of low magnitude among individuals from the different areas sampled and a correlation between several SNP loci and environmental variables. The analyses carried out did not show a clear differentiation among individuals sampled along different salinity regimens, but enriched GOTerms seem to be driven by water oxygen content. On the other hand, a total of 46 SNPs analysed in O. argentinensis and O. bonariensis showed a correlation with salinity and temperature. Although none of the correlated SNPs and corresponding genes from our both analyses were directly associated with hypoxia, genes related to the cardiovascular system and muscle cell differentiation were found. All these genes are interesting candidates for future studies since they are closely related to the differentially expressed genes. Although salinity was also mentioned as an important parameter limiting introgression between O. argentinensis and O. bonariensis, it was found that salinity does not drive differential expression in O. argentinensis, but rather oxygen levels. Moreover, salinity does not directly affect the structure and genetic divergence of the populations, they appear to be structured based on their degree of isolation and geographical distance between them. Further studies, like genome-wide analyses, could help to elucidate additional genes adapted to the different environments in these silverside species.

A historical stepping-stone path for an island-colonizing cactus across a submerged "bridge" archipelago

Here we use population genomic data (ddRAD-Seq) and ecological niche modeling to test biogeographic hypotheses for the divergence of the island-endemic cactus species Cereus insularis Hemsl. (Cereeae; Cactaceae) from its sister species C. fernambucensis Lem. The Cereus insularis grows in the Fernando de Noronha Islands (FNI), a Neotropical archipelago located 350 km off the Brazilian Atlantic Forest (BAF) coast. Phylogeographic reconstructions support a northward expansion by the common ancestor of C. insularis and C. fernambucensis along the mainland BAF coast, with C. insularis diverging from the widespread mainland taxon C. fernambucensis after colonizing FNI in the late Pleistocene. The morphologically distinct C. insularis is monophyletic and nested within C. fernambucensis, as expected from a progenitor-derivative speciation model. We tested alternative biogeographic and demographic hypotheses for the colonization of the FNI using Approximate Bayesian Computation. We found the greatest support for a stepping-stone path that emerged during periods of decreased sea level (the "bridge" hypothesis), in congruence with historical ecological niche modeling that shows highly suitable habitats on stepping-stone islands during glacial periods. The outlier analyses reveal signatures of selection in C. insularis, suggesting a putative role of adaptation driving rapid anagenic differentiation of this species in FNI.

Severe Bottleneck Impacted the Genomic Structure of Egg-Eating Cichlids in Lake Victoria

Within 15,000 years, the explosive adaptive radiation of haplochromine cichlids in Lake Victoria, East Africa, generated 500 endemic species. In the 1980s, the upsurge of Nile perch, a carnivorous fish artificially introduced to the lake, drove the extinction of more than 200 endemic cichlids. The Nile perch predation particularly harmed piscivorous cichlids, including paedophages, cichlids eat eggs and fries, which is an example of the unique trophic adaptation seen in African cichlids. Here, aiming to investigate past demographic events possibly triggered by the invasion of Nile perch and the subsequent impacts on the genetic structure of cichlids, we conducted large-scale comparative genomics. We discovered evidence of recent bottleneck events in 4 species, including 2 paedophages, which began during the 1970s to 1980s, and population size rebounded during the 1990s to 2000s. The timing of the bottleneck corresponded to the historical records of endemic haplochromines" disappearance and later resurgence, which is likely associated with the introduction of Nile perch by commercial demand to Lake Victoria in the 1950s. Interestingly, among the 4 species that likely experienced bottleneck, Haplochromis sp. "matumbi hunter," a paedophagous cichlid, showed the most severe bottleneck signatures. The components of shared ancestry inferred by ADMIXTURE suggested a high genetic differentiation between matumbi hunter and other species. In contrast, our phylogenetic analyses highly supported the monophyly of the 5 paedophages, consistent with the results of previous studies. We conclude that high genetic differentiation of matumbi hunter occurred due to the loss of shared genetic components among haplochromines in Lake Victoria caused by the recent severe bottleneck.

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.

Divergence and reticulation in the Mexican white oaks: ecological and phylogenomic evidence on species limits and phylogenetic networks in the Quercus laeta complex (Fagaceae)

BACKGROUND AND AIMS

Introgressive hybridization poses a challenge to taxonomic and phylogenetic understanding of taxa, particularly when there are high numbers of co-occurring, intercrossable species. The genus Quercus exemplifies this situation. Oaks are highly diverse in sympatry and cross freely, creating syngameons of interfertile species. Although a well-resolved, dated phylogeny is available for the American oak clade, evolutionary relationships within many of the more recently derived clades remain to be defined, particularly for the young and exceptionally diverse Mexican white oak clade. Here, we adopted an approach bridging micro- and macroevolutionary scales to resolve evolutionary relationships in a rapidly diversifying clade endemic to Mexico.

METHODS

Ecological data and sequences of 155 low-copy nuclear genes were used to identify distinct lineages within the Quercus laeta complex. Concatenated and coalescent approaches were used to assess the phylogenetic placement of these lineages relative to the Mexican white oak clade. Phylogenetic network methods were applied to evaluate the timing and genomic significance of recent or historical introgression among lineages.

KEY RESULTS

The Q. laeta complex comprises six well-supported lineages, each restricted geographically and with mostly divergent climatic niches. Species trees corroborated that the different lineages are more closely related to other species of Mexican white oaks than to each other, suggesting that this complex is polyphyletic. Phylogenetic networks estimated events of ancient introgression that involved the ancestors of three present-day Q. laeta lineages.

CONCLUSIONS

The Q. laeta complex is a morphologically and ecologically related group of species rather than a clade. Currently, oak phylogenetics is at a turning point, at which it is necessary to integrate phylogenetics and ecology in broad regional samples to figure out species boundaries. Our study illuminates one of the more complicated of the Mexican white oak groups and lays groundwork for further taxonomic study.

The sweet tabaiba or there and back again: phylogeographical history of the Macaronesian Euphorbia balsamifera

BACKGROUND AND AIMS

Biogeographical relationships between the Canary Islands and north-west Africa are often explained by oceanic dispersal and geographical proximity. Sister-group relationships between Canarian and eastern African/Arabian taxa, the 'Rand Flora' pattern, are rare among plants and have been attributed to the extinction of north-western African populations. Euphorbia balsamifera is the only representative species of this pattern that is distributed in the Canary Islands and north-west Africa; it is also one of few species present in all seven islands. Previous studies placed African populations of E. balsamifera as sister to the Canarian populations, but this relationship was based on herbarium samples with highly degraded DNA. Here, we test the extinction hypothesis by sampling new continental populations; we also expand the Canarian sampling to examine the dynamics of island colonization and diversification.

METHODS

Using target enrichment with genome skimming, we reconstructed phylogenetic relationships within E. balsamifera and between this species and its disjunct relatives. A single nucleotide polymorphism dataset obtained from the target sequences was used to infer population genetic diversity patterns. We used convolutional neural networks to discriminate among alternative Canary Islands colonization scenarios.

KEY RESULTS

The results confirmed the Rand Flora sister-group relationship between western E. balsamifera and Euphorbia adenensis in the Eritreo-Arabian region and recovered an eastern-western geographical structure among E. balsamifera Canarian populations. Convolutional neural networks supported a scenario of east-to-west island colonization, followed by population extinctions in Lanzarote and Fuerteventura and recolonization from Tenerife and Gran Canaria; a signal of admixture between the eastern island and north-west African populations was recovered.

CONCLUSIONS

Our findings support the Surfing Syngameon Hypothesis for the colonization of the Canary Islands by E. balsamifera, but also a recent back-colonization to the continent. Populations of E. balsamifera from northwest Africa are not the remnants of an ancestral continental stock, but originated from migration events from Lanzarote and Fuerteventura. This is further evidence that oceanic archipelagos are not a sink for biodiversity, but may be a source of new genetic variability.

Genomic signatures of strawberry domestication and diversification

Cultivated strawberry (Fragaria × ananassa) has a brief history of less than 300 yr, beginning with the hybridization of octoploids Fragaria chiloensis and Fragaria virginiana. Here we explored the genomic signatures of early domestication and subsequent diversification for different climates using whole-genome sequences of 289 wild, heirloom, and modern varieties from two major breeding programs in the United States. Four nonadmixed wild octoploid populations were identified, with recurrent introgression among the sympatric populations. The proportion of F. virginiana ancestry increased by 20% in modern varieties over initial hybrids, and the proportion of F. chiloensis subsp. pacifica rose from 0% to 3.4%. Effective population size rapidly declined during early breeding. Meanwhile, divergent selection for distinct environments reshaped wild allelic origins in 21 out of 28 chromosomes. Overlapping divergent selective sweeps in natural and domesticated populations revealed 16 convergent genomic signatures that may be important for climatic adaptation. Despite 20 breeding cycles since initial hybridization, more than half of loci underlying yield and fruit size are still not under artificial selection. These insights add clarity to the domestication and breeding history of what is now the most widely cultivated fruit in the world.

Genetic legacies of mega-landslides: Cycles of isolation and contact across flank collapses in an oceanic island

Catastrophic flank collapses are recognized as important drivers of insular biodiversity dynamics, through the disruption of species ranges and subsequent allopatric divergence. However, little empirical data supports this conjecture, with their evolutionary consequences remaining poorly understood. Using genome-wide data within a population genomics and phylogenomics framework, we evaluate how mega-landslides have impacted evolutionary and demographic history within a species complex of weevils (Curculionidae) within the Canary Island of Tenerife. We reveal a complex genomic landscape, within which individuals of single ancestry were sampled in areas characterized by long-term geological stability, relative to the timing of flank collapses. In contrast, individuals of admixed ancestry were almost exclusively sampled within the boundaries of flank collapses. Estimated divergence times among ancestral populations aligned with the timings of mega-landslide events. Our results provide first evidence for a cyclical dynamic of range fragmentation and secondary contact across flank collapse landscapes, with support for a model where this dynamic is mediated by Quaternary climate oscillations. The context within which we reveal climate and topography to interact cyclically through time to shape the geographic structure of genetic variation, together with related recent work, highlights the importance of topoclimatic phenomena as an agent of diversification within insular invertebrates.

Parallel patterns of genetic diversity and structure in circumboreal species of the Sphagnum capillifolium complex

PREMISE

Shared geographical patterns of population genetic variation among related species is a powerful means to identify the historical events that drive diversification. The Sphagnum capillifolium complex is a group of closely related peat mosses within the Sphagnum subgenus Acutifolia and contains several circumboreal species whose ranges encompass both glaciated and unglaciated regions across the northern hemisphere. In this paper, we (1) inferred the phylogeny of subg. Acutifolia and (2) investigated patterns of population structure and genetic diversity among five circumboreal species within the S. capillifolium complex.

METHODS

We generated RAD sequencing data from most species of the subg. Acutifolia and samples from across the distribution ranges of circumboreal species within the S. capillifolium complex.

RESULTS

We resolved at least 14 phylogenetic clusters within the S. capillifolium complex. Five circumboreal species show some common patterns: One population system comprises plants in eastern North America and Europe, and another comprises plants in the Pacific Northwest or around the Beringian and Arctic regions. Alaska appears to be a hotspot for genetic admixture, genetic diversity, and sometimes endemic subclades.

CONCLUSIONS

Our results support the hypothesis that populations of five circumboreal species within the S. capillifolium complex survived in multiple refugia during the last glacial maximum. Long-distance dispersal out of refugia, population bottlenecks, and possible adaptations to conditions unique to each refugium could have contributed to current geographic patterns. These results indicate the important role of historical events in shaping the complex population structure of plants with broad distribution ranges.

Implications of gene tree heterogeneity on downstream phylogenetic analyses: A case study employing the Fair Proportion index

Many questions in evolutionary biology require the specification of a phylogeny for downstream phylogenetic analyses. However, with the increasingly widespread availability of genomic data, phylogenetic studies are often confronted with conflicting signal in the form of genomic heterogeneity and incongruence between gene trees and the species tree. This raises the question of determining what data and phylogeny should be used in downstream analyses, and to what extent the choice of phylogeny (e.g., gene trees versus species trees) impacts the analyses and their outcomes. In this paper, we study this question in the realm of phylogenetic diversity indices, which provide ways to prioritize species for conservation based on their relative evolutionary isolation on a phylogeny, and are thus one example of downstream phylogenetic analyses. We use the Fair Proportion (FP) index, also known as the evolutionary distinctiveness score, and explore the variability in species rankings based on gene trees as compared to the species tree for several empirical data sets. Our results indicate that prioritization rankings among species vary greatly depending on the underlying phylogeny, suggesting that the choice of phylogeny is a major influence in assessing phylogenetic diversity in a conservation setting. While we use phylogenetic diversity conservation as an example, we suspect that other types of downstream phylogenetic analyses such as ancestral state reconstruction are similarly affected by genomic heterogeneity and incongruence. Our aim is thus to raise awareness of this issue and inspire new research on which evolutionary information (species trees, gene trees, or a combination of both) should form the basis for analyses in these settings.

Phylogenomics of American pika (Ochotona princeps) lineage diversification

Quaternary climate oscillations have profoundly influenced current species distributions. For many montane species, these fluctuations were a prominent driver in species range shifts, often resulting in intraspecific diversification, as has been the case for American pikas (Ochotona princeps). Range shifts and population declines in this thermally-sensitive lagomorph have been linked to historical and contemporary environmental changes across its western North American range, with previous research reconstructing five mitochondrial DNA lineages. Here, we paired genome-wide data (25,244 SNPs) with range-wide sampling to re-examine the number and distribution of intra-specific lineages, and investigate patterns of within- and among-lineage divergence and diversity. Our results provide genomic evidence of O. princeps monophyly, reconstructing six distinct lineages that underwent multiple rounds of divergence (0.809-2.81 mya), including a new Central Rocky Mountain lineage. We further found evidence for population differentiation across multiple spatial scales, and reconstructed levels of standing variation comparable to those found in other small mammals. Overall, our findings demonstrate the influence of past glacial cycles on O. princeps lineage diversification, suggest that current subspecific taxonomy may need to be revisited, and provide an important framework for investigations of American pika adaptive potential in the face of anthropogenic climate change.

Whole genomes show contrasting trends of population size changes and genomic diversity for an Amazonian endemic passerine over the late quaternary

The "Amazon tipping point" is a global change scenario resulting in replacement of upland terra-firme forests by large-scale "savannization" of mostly southern and eastern Amazon. Reduced rainfall accompanying the Last Glacial Maximum (LGM) has been proposed to have acted as such a tipping point in the past, with the prediction that terra-firme inhabiting species should have experienced reductions in population size as drier habitats expanded. Here, we use whole-genomes of an Amazonian endemic organism (Scale-backed antbirds - Willisornis spp.) sampled from nine populations across the region to test this historical demography scenario. Populations from southeastern Amazonia and close to the Amazon-Cerrado ecotone exhibited a wide range of demographic patterns, while most of those from northern and western Amazonia experienced uniform expansions between 400 kya and 80-60 kya, with gradual declines toward 20 kya. Southeastern populations of Willisornis were the last to diversify and showed smaller heterozygosity and higher runs of homozygosity values than western and northern populations. These patterns support historical population declines throughout the Amazon that affected more strongly lineages in the southern and eastern areas, where historical "tipping point" conditions existed due to the widespread replacement of humid forest by drier and open vegetation during the LGM.

Multiple transitions between realms shape relict lineages of Proteus cave salamanders

In comparison to biodiversity on Earth's surface, subterranean biodiversity has largely remained concealed. The olm (Proteus anguinus) is one of the most enigmatic extant cave inhabitants, and until now little was known regarding its genetic structure and evolutionary history. Olms inhabit subterranean waters throughout the Dinaric Karst of the western Balkans, with a seemingly uniform phenotypic appearance of cave-specialized traits: an elongate body, snout and limbs, degenerated eyes and loss of pigmentation ("white olm"). Only a single small region in southeastern Slovenia harbours olms with a phenotype typical of surface animals: pigmented skin, eyes, a blunt snout and short limbs ("black olm"). We used a combination of mitochondrial DNA and genome-wide single nucleotide polymorphism data to investigate the molecular diversity, evolutionary history and biogeography of olms along the Dinaric Karst. We found nine deeply divergent species-level lineages that separated between 17 and 4 million years ago, while molecular diversity within lineages was low. We detected no signal of recent admixture between lineages and only limited historical gene flow. Biogeographically, the contemporaneous distribution of lineages mostly mirrors hydrologically separated subterranean environments, while the historical separation of olm lineages follows microtectonic and climatic changes in the area. The reconstructed phylogeny suggests at least four independent transitions to the cave phenotype. Two of the species-level lineages have miniscule ranges and may represent Europe's rarest amphibians. Their rarity and the decline in other lineages call for protection of their subterranean habitats.

Hybridization, polyploidization, and morphological convergence make dozens of taxa into one chaotic genetic pool: a phylogenomic case of the Ficus erecta species complex (Moraceae)

The Ficus erecta complex, characterized by its morphological diversity and frequent interspecific overlap, shares pollinating fig wasps among several species. This attribute, coupled with its intricate phylogenetic relationships, establishes it as an exemplary model for studying speciation and evolutionary patterns. Extensive researches involving RADseq (Restriction-site associated DNA sequencing), complete chloroplast genome data, and flow cytometry methods were conducted, focusing on phylogenomic analysis, genetic structure, and ploidy detection within the complex. Significantly, the findings exposed a pronounced nuclear-cytoplasmic conflict. This evidence, together with genetic structure analysis, confirmed that hybridization within the complex is a frequent occurrence. The ploidy detection revealed widespread polyploidy, with certain species exhibiting multiple ploidy levels, including 2×, 3×, and 4×. Of particular note, only five species (F. abelii, F. erecta, F. formosana, F. tannoensis and F. vaccinioides) in the complex were proved to be monophyletic. Species such as F. gasparriniana, F. pandurata, and F. stenophylla were found to encompass multiple phylogenetically distinct lineages. This discovery, along with morphological comparisons, suggests a significant underestimation of species diversity within the complex. This study also identified F. tannoensis as an allopolyploid species originating from F. vaccinioide and F. erecta. Considering the integration of morphological, molecular systematics, and cytological evidences, it is proposed that the scope of the F. erecta complex should be expanded to the entire subsect. Frutescentiae. This would redefine the complex as a continuously evolving group comprising at least 33 taxa, characterized by blurred species boundaries, frequent hybridization and polyploidization, and ambiguous genetic differentiation.

Median quartet tree search algorithms using optimal subtree prune and regraft

Gene trees can be different from the species tree due to biological processes and inference errors. One way to obtain a species tree is to find one that maximizes some measure of similarity to a set of gene trees. The number of shared quartets between a potential species tree and gene trees provides a statistically justifiable score; if maximized properly, it could result in a statistically consistent estimator of the species tree under several statistical models of discordance. However, finding the median quartet score tree, one that maximizes this score, is NP-Hard, motivating several existing heuristic algorithms. These heuristics do not follow the hill-climbing paradigm used extensively in phylogenetics. In this paper, we make theoretical contributions that enable an efficient hill-climbing approach. Specifically, we show that a subtree of size m can be placed optimally on a tree of size n in quasi-linear time with respect to n and (almost) independently of m. This result enables us to perform subtree prune and regraft (SPR) rearrangements as part of a hill-climbing search. We show that this approach can slightly improve upon the results of widely-used methods such as ASTRAL in terms of the optimization score but not necessarily accuracy.

Japanese wolves are most closely related to dogs and share DNA with East Eurasian dogs

Although the domestic dog's origin is still unclear, this lineage is believed to have been domesticated from an extinct population of gray wolves, which is expected to be more closely related to dogs than to other populations of gray wolves. Here, we sequence the whole genomes of nine Japanese wolves (7.5-100x: Edo to Meiji periods) and 11 modern Japanese dogs and analyze them together with those from other populations of dogs and wolves. A phylogenomic tree shows that, among the gray wolves, Japanese wolves are closest to the dog, suggesting that the ancestor of dogs is closely related to the ancestor of the Japanese wolf. Based on phylogenetic and geographic relationships, the dog lineage has most likely originated in East Asia, where it diverged from a common ancestor with the Japanese wolf. Since East Eurasian dogs possess Japanese wolf ancestry, we estimate an introgression event from the ancestor of the Japanese wolf to the ancestor of the East Eurasian dog that occurred before the dog's arrival in the Japanese archipelago.

Evolutionary relationships, hybridization and diversification under domestication of the locoto chile (Capsicum pubescens) and its wild relatives

Patterns of genetic variation in crops are the result of multiple processes that have occurred during their domestication and improvement, and are influenced by their wild progenitors that often remain understudied. The locoto chile, Capsicum pubescens, is a crop grown mainly in mid-highlands of South-Central America. This species is not known from the wild and exists only as a cultigen. The evolutionary affinities and exact origin of C. pubescens have still not been elucidated, with hypotheses suggesting its genetic relatedness and origin to two wild putative ancestral Capsicum species from the Central Andes, C. eximium and C. cardenasii. In the current study, RAD-sequencing was applied to obtain genome-wide data for 48 individuals of C. pubescens and its wild allies representing different geographical areas. Bayesian, Maximum Likelihood and coalescent-based analytical approaches were used to reconstruct population genetic patterns and phylogenetic relationships of the studied species. The results revealed that C. pubescens forms a well-defined monotypic lineage closely related to wild C. cardenasii and C. eximium, and also to C. eshbaughii. The primary lineages associated with the diversification under domestication of C. pubescens were also identified. Although direct ancestor-descendant relationship could not be inferred within this group of taxa, hybridization events were detected between C. pubescens and both C. cardenasii and C. eximium. Therefore, although hybrid origin of C. pubescens could not be inferred, gene flow involving its wild siblings was shown to be an important factor contributing to its contemporary genetic diversity. The data allowed for the inference of the center of origin of C. pubescens in central-western Bolivia highlands and for better understanding of the dynamics of its gene pool. The results of this study are essential for germplasm conservation and breeding purposes, and provide excellent basis for further research of the locoto chile and its wild relatives.

The natural and human-mediated expansion of a human-commensal lizard into the fringes of Southeast Asia

BACKGROUND

Human-commensal species often display deep ancestral genetic structure within their native range and founder-effects and/or evidence of multiple introductions and admixture in newly established areas. We investigated the phylogeography of Eutropis multifasciata, an abundant human-commensal scincid lizard that occurs across Southeast Asia, to determine the extent of its native range and to assess the sources and signatures of human introduction outside of the native range. We sequenced over 350 samples of E. multifasciata for the mitochondrial ND2 gene and reanalyzed a previous RADseq population genetic dataset in a phylogenetic framework.

RESULTS

Nuclear and mitochondrial trees are concordant and show that E. multifasciata has retained high levels of genetic structure across Southeast Asia despite being frequently moved by humans. Lineage boundaries in the native range roughly correspond to several major biogeographic barriers, including Wallace's Line and the Isthmus of Kra. Islands at the outer fringe of the range show evidence of founder-effects and multiple introductions.

CONCLUSIONS

Most of enormous range of E. multifasciata across Southeast Asia is native and it only displays signs of human-introduction or recent expansion along the eastern and northern fringe of its range. There were at least three events of human-introductions to Taiwan and offshore islands, and several oceanic islands in eastern Indonesia show a similar pattern. In Myanmar and Hainan, there is a founder-effect consistent with post-warming expansion after the last glacial maxima or human introduction.

Anomalous networks under the multispecies coalescent: theory and prevalence

Reticulations in a phylogenetic network represent processes such as gene flow, admixture, recombination and hybrid speciation. Extending definitions from the tree setting, an anomalous network is one in which some unrooted tree topology displayed in the network appears in gene trees with a lower frequency than a tree not displayed in the network. We investigate anomalous networks under the Network Multispecies Coalescent Model with possible correlated inheritance at reticulations. Focusing on subsets of 4 taxa, we describe a new algorithm to calculate quartet concordance factors on networks of any level, faster than previous algorithms because of its focus on 4 taxa. We then study topological properties required for a 4-taxon network to be anomalous, uncovering the key role of [Formula: see text]-cycles: cycles of 3 edges parent to a sister group of 2 taxa. Under the model of common inheritance, that is, when each gene tree coalesces within a species tree displayed in the network, we prove that 4-taxon networks are never anomalous. Under independent and various levels of correlated inheritance, we use simulations under realistic parameters to quantify the prevalence of anomalous 4-taxon networks, finding that truly anomalous networks are rare. At the same time, however, we find a significant fraction of networks close enough to the anomaly zone to appear anomalous, when considering the quartet concordance factors observed from a few hundred genes. These apparent anomalies may challenge network inference methods.

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

Glacial vicariance and secondary contact shape demographic histories in a freshwater mussel species complex

Characterizing the mechanisms influencing the distribution of genetic variation in aquatic species can be difficult due to the dynamic nature of hydrological landscapes. In North America's Central Highlands, a complex history of glacial dynamics, long-term isolation, and secondary contact have shaped genetic variation in aquatic species. Although the effects of glacial history have been demonstrated in many taxa, responses are often lineage- or species-specific and driven by organismal ecology. In this study, we reconstruct the evolutionary history of a freshwater mussel species complex using a suite of mitochondrial and nuclear loci to resolve taxonomic and demographic uncertainties. Our findings do not support Pleurobema rubrum as a valid species, which is proposed for listing as threatened under the U.S. Endangered Species Act. We synonymize P. rubrum under Pleurobema sintoxia-a common and widespread species found throughout the Mississippi River Basin. Further investigation of patterns of genetic variation in P. sintoxia identified a complex demographic history, including ancestral vicariance and secondary contact, within the Eastern Highlands. We hypothesize these patterns were shaped by ancestral vicariance driven by the formation of Lake Green and subsequent secondary contact after the last glacial maximum. Our inference aligns with demographic histories observed in other aquatic taxa in the region and mirrors patterns of genetic variation of a freshwater fish species (Erimystax dissimilis) confirmed to serve as a parasitic larval host for P. sintoxia. Our findings directly link species ecology to observed patterns of genetic variation and may have significant implications for future conservation and recovery actions of freshwater mussels.

Rapid diversification of grey mangroves (Avicennia marina) driven by geographic isolation and extreme environmental conditions in the Arabian Peninsula

Biological systems occurring in ecologically heterogeneous and spatially discontinuous habitats provide an ideal opportunity to investigate the relative roles of neutral and selective factors in driving lineage diversification. The grey mangroves (Avicennia marina) of Arabia occur at the northern edge of the species' range and are subject to variable, often extreme, environmental conditions, as well as historic large fluctuations in habitat availability and connectivity resulting from Quaternary glacial cycles. Here, we analyse fully sequenced genomes sampled from 19 locations across the Red Sea, the Arabian Sea and the Persian/Arabian Gulf (PAG) to reconstruct the evolutionary history of the species in the region and to identify adaptive mechanisms of lineage diversification. Population structure and phylogenetic analyses revealed marked genetic structure correlating with geographic distance and highly supported clades among and within the seas surrounding the Arabian Peninsula. Demographic modelling showed times of divergence consistent with recent periods of geographic isolation and low marine connectivity during glaciations, suggesting the presence of (cryptic) glacial refugia in the Red Sea and the PAG. Significant migration was detected within the Red Sea and the PAG, and across the Strait of Hormuz to the Arabian Sea, suggesting gene flow upon secondary contact among populations. Genetic-environment association analyses revealed high levels of adaptive divergence and detected signs of multi-loci local adaptation driven by temperature extremes and hypersalinity. These results support a process of rapid diversification resulting from the combined effects of historical factors and ecological selection and reveal mangrove peripheral environments as relevant drivers of lineage diversity.

Phylogenetic, morphological and niche differentiation unveil new species limits for the big brown bat (Eptesicus fuscus)

Phylogeographic accounts of mammals across fragmented landscapes show high levels of genetic, morphological and ecological variation. The big brown bat (Eptesicus fuscus) widely spans mainland landmasses from Canada to Ecuador and Colombia, and the insular Caribbean through The Bahamas and Greater and Lesser Antilles. Given the distribution of E. fuscus, we hypothesized that insular lineages could represent a different species aided by isolation. We assessed species limits by capitalizing on available mitochondrial and genomic data. Novel morphological and spatial datasets were produced to examine limits phenotypically and whether ecological niches could be associated with differences between groups. Phylogenetics strongly supported the Caribbean as unique compared to the mainland. Genomic data indicated high levels of genetic structure within the Caribbean and no detectable admixture of the Caribbean with continental lineages. Similarly, the Caribbean group shows high phenotypic disparity, and niche models revealed differences in habitat suitability between groups, concordant with the phylogenetic results. This study uncovered signals of divergence supporting the Caribbean clade of E. fuscus as unique through an integrative framework. We endorse re-evaluating the taxonomic status of Caribbean big brown bats as Eptesicus dutertreus. This recognition can help promote local conservation plans for insular lineages of big brown bats.

Phylogenomics and species delimitation of an abundant and little-studied Amazonian forest spiny rat

Species delimitation studies based on integrating different datasets such as genomic, morphometric, and cytogenetics data are rare in studies focused on Neotropical rodents. As a consequence, the evolutionary history of most of these genera remains poorly understood. Proechimys is a highly diverse and widely distributed genus of Neotropical spiny rats with unique traits like multiple sympatry, micro-habitat segregation, and fuzzy species limits. Here, we applied RAD-Seq to infer the phylogenetic relationships, estimate the species boundaries, and estimate the divergence times for Proechimys, one of the most common and least studied small mammals in the Amazon. We tested whether inferred lineages in the phylogenetic trees could be considered distinct species based on the genomic dataset and morphometric data. Analyses revealed the genus is not monophyletic, with Proechimys hoplomyoides sister to a group of Hoplomys gymnurus + all other Proechimys species, contesting the generic status of Hoplomys. There are five main clades in Proechimys stricto sensu (excluding H. gymnurus and P. hoplomyoides). Species delimitation analyses supported 25 species within the genus Proechimys. The five main clades in Proechimys stricto sensu also showed similar ages for their origins, and two rapid diversification events were identified in the Early Pliocene and in the Early Pleistocene. Most cases of sympatry in Proechimys occur among species from the different main clades, and although Proechimys is an inhabitant of the Amazon, three species occupied the Cerrado biome during the Pleistocene. We could associate available nominal taxon, cytogenetics information, and DNA sequences in Genbank to most of the 25 species we hypothesized from our delimitation analyses. Based on our analyses, we estimate that eight forms represent putative new species that need a taxonomic revision.

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

Detection of trait donors and QTL boundaries for early blight resistance using local ancestry inference in a library of genomic sequences for tomato

By conducting hierarchical clustering along a sliding window, we generated haplotypes across hundreds of re-sequenced genomes in a few hours. We leveraged our method to define cryptic introgressions underlying disease resistance in tomato (Solanum lycopersicum L.) and to discover resistant germplasm in the tomato seed bank. The genomes of 9 accessions with early blight (Alternaria linariae) disease resistance were newly sequenced and analyzed together with published sequences for 770 tomato and wild species accessions, most of which are available in germplasm collections. Identification of common ancestral haplotypes among resistant germplasm enabled rapid fine mapping of recently discovered quantitative trait loci (QTL) conferring resistance and the identification of possible causal variants. The source of the early blight QTL EB-9 was traced to a vintage tomato named 'Devon Surprise'. Another QTL, EB-5, as well as resistance to bacterial spot disease (Xanthomonas spp.), was traced to Hawaii 7998. A genomic survey of all accessions forecasted EB-9-derived resistance in several heirloom tomatoes, accessions of S. lycopersicum var. cerasiforme, and S. pimpinellifolium PI 37009. Our haplotype-based predictions were validated by screening the accessions against the causal pathogen. There was little evidence of EB-5 prevalence in surveyed contemporary germplasm, presenting an opportunity to bolster tomato disease resistance by adding this rare locus. Our work demonstrates practical insights that can be derived from the efficient processing of large genome-scale datasets, including rapid functional prediction of disease resistance QTL in diverse genetic backgrounds. Finally, our work finds more efficient ways to leverage public genetic resources for crop improvement.

Reticulate evolution of the tertiary relict Osmanthus

When interspecific gene flow is common, species relationships are more accurately represented by a phylogenetic network than by a bifurcating tree. This study aimed to uncover the role of introgression in the evolution of Osmanthus, the only genus of the subtribe Oleinae (Oleaceae) with its distribution center in East Asia. We built species trees, detected introgression, and constructed networks using multiple kinds of sequencing data (whole genome resequencing, transcriptome sequencing, and Sanger sequencing of nrDNA) combined with concatenation and coalescence approaches. Then, based on well-understood species relationships, historical biogeographic analyses and diversification rate estimates were employed to reveal the history of Osmanthus. Osmanthus originated in mid-Miocene Europe and dispersed to the eastern Tibetan Plateau in the late Miocene. Thereafter, it continued to spread eastwards. Phylogenetic conflict is common within the 'Core Osmanthus' clade and is seen at both early and late stages of diversification, leading to hypotheses of net-like species relationships. Incomplete lineage sorting proved ineffective in explaining phylogenetic conflicts and thus supported introgression as the main cause of conflicts. This study elucidates the diversification history of a relict genus in the subtropical regions of eastern Asia and reveals that introgression had profound effects on its evolutionary history.

Phylogenomics of the narrowly endemic Eurycheilichthys (Siluriformes: Loricariidae): Sympatric species with non-sister relationships suggest mainly allopatric speciation

Armored catfishes of the genus Eurycheilichthys are endemic to Southern Brazil and Misiones (Argentina) comprising nine species of small size, with a high degree of sympatry and species diversity distributed in two river basins. Here we use new genome-wide data to infer a species phylogeny and test species boundaries for this poorly known group. We estimate 1) the phylogenetic relationships of the species of Eurycheilichthys based on 29,350 loci in 65 individuals of nine species plus outgroups, and 2) the population structure and differentiation based on 43,712 loci and 62 individuals to estimate how geography may have acted on speciation and formation of the sympatric species groups. Analyses support the monophyly of the genus and suggest two species-inclusive clades (East and West) with high support and very recently diverged species. Western clade contains E. limulus (from upper Jacuí River basin) that is sister to Western species of the Taquari-Antas basin plus E. paucidens. The Eastern clade contains E. pantherinus (from Uruguay River basin) sister to the Eastern species of the Taquari-Antas basin E. coryphaenus, plus the central-distributed species E. planus and E. vacariensis, and the more widely-distributed species E. luisae. Eurycheilichthys luisae is not monophyletic and may contain one or more cryptic species or hybrid individuals. A stronger diversity on structure of lineages on the Taquari-Antas, when compared to upper Uruguay and Jacuí River basins, and the fact that most of the sympatrically distributed taxa have non-sister relationships suggest a scenario of mainly allopatric speciation and may indicate a more dynamic landscape with headwater capture events among these tributaries.

Phylogenomic analyses reveal reticulate evolution between Neomicrocalamus and Temochloa (Poaceae: Bambusoideae)

Neomicrocalamus and Temochloa are closely related to bamboo genera. However, when considered with newly discovered and morphologically similar material from China and Vietnam, the phylogenetic relationship among these three groups was ambiguous in the analyses based on DNA regions. Here, as a means of investigating the relationships among the three bamboo groups and exploring potential sources of genomic conflicts, we present a phylogenomic examination based on the whole plastome, single-nucleotide polymorphism (SNP), and single-copy nuclear (SCN) gene datasets. Three different phylogenetic hypotheses were found. The inconsistency is attributed to the combination of incomplete lineage sorting and introgression. The origin of newly discovered bamboos is from introgressive hybridization between Temochloa liliana (which contributed 80.7% of the genome) and Neomicrocalamus prainii (19.3%), indicating that the newly discovered bamboos are closer to T. liliana in genetics. The more similar morphology and closer distribution elevation also imply a closer relationship between Temochloa and newly discovered bamboos.

Phylogenomic investigation of lampreys (Petromyzontiformes)

The history of lamprey evolution has been contentious due to limited morphological differentiation and limited genetic data. Available data has produced inconsistent results, including in the relationship among northern and southern species and the monophyly of putative clades. Here we use whole genome sequence data sourced from a public database to identify orthologs for 11 lamprey species from across the globe and build phylogenies. The phylogeny showed a clear separation between northern and southern lamprey species, which contrasts with some prior work. We also find that the phylogenetic relationships of our samples of two genera, Lethenteron and Eudontomyzon, deviate from the taxonomic classification of these species, suggesting that they require reclassification.

Unraveling phylogenetic relationships and species boundaries in the arid adapted Gerbillus rodents (Muridae: Gerbillinae) by RAD-seq data

Gerbillus is one of the most speciose genera among rodents, with ca. 51 recognized species. Previous attempts to reconstruct the evolutionary history of Gerbillus mainly relied on the mitochondrial cyt-b marker as a source of phylogenetic information. In this study, we utilize RAD-seq genomic data from 37 specimens representing 11 species to reconstruct the phylogenetic tree for Gerbillus, applying concatenation and coalescence methods. We identified four highly supported clades corresponding to the traditionally recognized subgenera: Dipodillus, Gerbillus, Hendecapleura and Monodia. Only two uncertain branches were detected in the resulting trees, with one leading to diversification of the main lineages in the genus, recognized by quartet sampling analysis as uncertain due to possible introgression. We also examined species boundaries for four pairs of sister taxa, including potentially new species from Morocco, using SNAPP. The results strongly supported a speciation model in which all taxa are treated as separate species. The dating analyses confirmed the Plio-Pleistocene diversification of the genus, with the uncertain branch coinciding with the beginning of aridification of the Sahara at the the Plio-Pleistocene boundary. This study aligns well with the earlier analyses based on the cyt-b marker, reaffirming its suitability as an adequate marker for estimating genetic diversity in Gerbillus.

Phylogenomics of the extinct Heath Hen provides support for sex-biased introgression among extant prairie grouse

Rapid divergence and subsequent reoccurring patterns of gene flow can complicate our ability to discern phylogenetic relationships among closely related species. To what degree such patterns may differ across the genome can provide an opportunity to extrapolate better how life history constraints may influence species boundaries. By exploring differences between autosomal and Z (or X) chromosomal-derived phylogenetic patterns, we can better identify factors that may limit introgression despite patterns of incomplete lineage sorting among closely related taxa. Here, using a whole-genome resequencing approach coupled with an exhaustive sampling of subspecies within the recently divergent prairie grouse complex (genus: Tympanuchus), including the extinct Heath Hen (T. cupido cupido), we show that their phylogenomic history differs depending on autosomal or Z-chromosome partitioned SNPs. Because the Heath Hen was allopatric relative to the other prairie grouse taxa, its phylogenetic signature should not be influenced by gene flow. In contrast, all the other extant prairie grouse taxa, except Attwater's Prairie-chicken (T. c. attwateri), possess overlapping contemporary geographic distributions and have been known to hybridize. After excluding samples that were likely translocated prairie grouse from the Midwest to the eastern coastal states or their resulting hybrids with mainland Heath Hens, species tree analyses based on autosomal SNPs consistently identified a paraphyletic relationship with regard to the Heath Hen with Lesser Prairie-chicken (T. pallidicinctus) sister to Greater Prairie-chicken (T. c. pinnatus) regardless of genic or intergenic partitions. In contrast, species trees based on the Z-chromosome were consistent with Heath Hen sister to a clade that included its conspecifics, Greater and Attwater's Prairie-chickens (T. c. attwateri). These results were further explained by historic gene flow, as shown with an excess of autosomal SNPs shared between Lesser and Greater Prairie-chickens but not with the Z-chromosome. Phylogenetic placement of Sharp-tailed Grouse (T. phasianellus), however, did not differ among analyses and was sister to a clade that included all other prairie grouse despite low levels of autosomal gene flow with Greater Prairie-chicken. These results, along with strong sexual selection (i.e., male hybrid behavioral isolation) and a lek breeding system (i.e., high variance in male mating success), are consistent with a pattern of female-biased introgression between prairie grouse taxa with overlapping geographic distributions. Additional study is warranted to explore how genomic components associated with the Z-chromosome influence the phenotype and thereby impact species limits among prairie grouse taxa despite ongoing contemporary gene flow.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Phylogeography, hybrid zones and contemporary species boundaries in the south-eastern Australian smooth frogs (Anura: Myobatrachidae: Geocrinia)

Paleo-climatic fluctuations have driven episodic changes in species distributions, providing opportunities for populations to diverge in isolation and hybridise following secondary contact. Studies of phylogeographic diversity and patterns of gene flow across hybrid zones can provide insight into contemporary species boundaries and help to inform taxonomic and conservation inferences. Here we explore geographic diversity within the acoustically divergent yet morphologically conserved south-eastern Australian smooth frog complex and assess gene flow across a narrow hybrid zone using mitochondrial nucleotide sequences and nuclear genome-wide single nucleotide polymorphisms. Our analyses reveal the presence of an evolutionarily distinct taxon restricted to the Otway Plains and Ranges, Victoria, which forms a narrow (9-30 km wide), spatiotemporally stable (>50 years) hybrid zone with Geocrinia laevis, which we describe herein as a new species.

Genomic signatures of ecological divergence between savanna and forest populations of a Neotropical tree

BACKGROUND AND AIMS

In eastern Neotropical South America, the Cerrado, a large savanna vegetation, and the Atlantic Forest harbour high biodiversity levels, and their habitats are rather different from each other. The biomes have intrinsic evolutionary relationships, with high lineage exchange that can be attributed, in part, to a large contact zone between them. The genomic study of ecotypes, i.e. populations adapted to divergent habitats, can be a model to study the genomic signatures of ecological divergence. Here, we investigated two ecotypes of the tree Plathymenia reticulata, one from the Cerrado and the other from the Atlantic Forest, which have a hybrid zone in the ecotonal zone of Atlantic Forest-Cerrado.

METHODS

The ecotypes were sampled in the two biomes and their ecotone. The evolutionary history of the divergence of the species was analysed with double-digest restriction site-associated DNA sequencing. The genetic structure and the genotypic composition of the hybrid zone were determined. Genotype-association analyses were performed, and the loci under putative selection and their functions were investigated.

KEY RESULTS

High divergence between the two ecotypes was found, and only early-generation hybrids were found in the hybrid zone, suggesting a partial reproductive barrier. Ancient introgression between the Cerrado and Atlantic Forest was not detected. The soil and climate were associated with genetic divergence in Plathymenia ecotypes and outlier loci were found to be associated with the stress response, with stomatal and root development and with reproduction.

CONCLUSIONS

The high genomic, ecological and morphophysiological divergence between ecotypes, coupled with partial reproductive isolation, indicate that the ecotypes represent two species and should be managed as different evolutionary lineages. We advise that the forest species should be re-evaluated and restated as vulnerable. Our results provide insights into the genomic mechanisms underlying the diversification of species across savanna and forest habitats and the evolutionary forces acting in the species diversification in the Neotropics.

Phylotranscriptomic analyses reveal deep gene tree discordance in Camellia (Theaceae)

Gene tree discordance is a significant legacy of biological evolution. Multiple factors can result in incongruence among genes, such as introgression, incomplete lineage sorting (ILS), gene duplication or loss. Resolving the background of gene tree discordance is a critical way to uncover the process of species diversification. Camellia, the largest genus in Theaceae, has controversial taxonomy and systematics due in part to a complex evolutionary history. We used 60 transcriptomes of 55 species, which represented 15 sections of Camellia to investigate its phylogeny and the possible causes of gene tree discordance. We conducted gene tree discordance analysis based on 1,617 orthologous low-copy nuclear genes, primarily using coalescent species trees and polytomy tests to distinguish hard and soft conflict. A selective pressure analysis was also performed to assess the impact of selection on phylogenetic topology reconstruction. Our results detected different levels of gene tree discordance in the backbone of Camellia, and recovered rapid diversification as one of the possible causes of gene tree discordance. Furthermore, we confirmed that none of the currently proposed sections of Camellia was monophyletic. Comparisons among datasets partitioned under different selective pressure regimes showed that integrating all orthologous genes provided the best phylogenetic resolution of the species tree of Camellia. The findings of this study reveal rapid diversification as a major source of gene tree discordance in Camellia and will facilitate future investigation of reticulate relationships at the species level in this important plant genus.

QR-STAR: A Polynomial-Time Statistically Consistent Method for Rooting Species Trees Under the Coalescent

We address the problem of rooting an unrooted species tree given a set of unrooted gene trees, under the assumption that gene trees evolve within the model species tree under the multispecies coalescent (MSC) model. Quintet Rooting (QR) is a polynomial time algorithm that was recently proposed for this problem, which is based on the theory developed by Allman, Degnan, and Rhodes that proves the identifiability of rooted 5-taxon trees from unrooted gene trees under the MSC. However, although QR had good accuracy in simulations, its statistical consistency was left as an open problem. We present QR-STAR, a variant of QR with an additional step and a different cost function, and prove that it is statistically consistent under the MSC. Moreover, we derive sample complexity bounds for QR-STAR and show that a particular variant of it based on "short quintets" has polynomial sample complexity. Finally, our simulation study under a variety of model conditions shows that QR-STAR matches or improves on the accuracy of QR. QR-STAR is available in open-source form on github.