A maximum pseudo-likelihood approach for phylogenetic networks

Peering through the hedge: Multiple datasets yield insights into the phylogenetic relationships and incongruences in the tribe Lilieae (Liliaceae)

The increasing use of genome-scale data has significantly facilitated phylogenetic analyses, contributing to the dissection of the underlying evolutionary mechanisms that shape phylogenetic incongruences, such as incomplete lineage sorting (ILS) and hybridization. Lilieae, a prominent member of the Liliaceae family, comprises four genera and approximately 260 species, representing 43% of all species within Liliaceae. They possess high ornamental, medicinal and edible values. Yet, no study has explored the validity of various genome-scale data in phylogenetic analyses within this tribe, nor have potential evolutionary mechanisms underlying its phylogenetic incongruences been investigated. Here, transcriptome, Angiosperms353, plastid and mitochondrial data, were collected from 50 to 93 samples of Lilieae, covering all four recognized genera. Multiple datasets were created and used for phylogenetic analyses based on concatenated and coalescent-based methods. Evolutionary rates of different datasets were calculated, and divergence times were estimated. Various approaches, including coalescence simulation, Quartet Sampling (QS), calculation of concordance factors (gCF and sCF), as well as MSCquartets and reticulate network inference, were carried out to infer the phylogenetic discordances and analyze their underlying mechanisms using a reduced 33-taxon dataset. Despite extensive phylogenetic discordances among gene trees, robust phylogenies were inferred from nuclear and plastid data compared to mitochondrial data, with lower synonymous substitution detected in mitochondrial genes than in nuclear and plastid genes. Significant ILS was detected across the phylogeny of Lilieae, with clear evidence of reticulate evolution identified. Divergence time estimation indicated that most of lineages in Lilieae diverged during a narrow time frame (ranging from 5.0 Ma to 10.0 Ma), consistent with the notion of rapid radiation evolution. Our results suggest that integrating transcriptomic and plastid data can serve as cost-effective and efficient tools for phylogenetic inference and evolutionary analysis within Lilieae, and Angiosperms353 data is also a favorable choice. Mitochondrial data are more suitable for phylogenetic analyses at higher taxonomic levels due to their stronger conservation and lower synonymous substitution rates. Significant phylogenetic incongruences detected in Lilieae were caused by both incomplete lineage sorting (ILS) and reticulate evolution, with hybridization and "ghost introgression" likely prevalent in the evolution of Lilieae species. Our findings provide new insights into the phylogeny of Lilieae, enhancing our understanding of the evolution of species in this tribe.

Phylotranscriptomic and ecological analyses reveal the evolution and morphological adaptation of Abies

Coniferous forests are under severe threat of the rapid anthropogenic climate warming. Abies (firs), the fourth-largest conifer genus, is a keystone component of the boreal and temperate dark-coniferous forests and harbors a remarkably large number of relict taxa. However, the uncertainty of the phylogenetic and biogeographic history of Abies significantly impedes our prediction of future dynamics and efficient conservation of firs. In this study, using 1,533 nuclear genes generated from transcriptome sequencing and a complete sampling of all widely recognized species, we have successfully reconstructed a robust phylogeny of global firs, in which four clades are strongly supported and all intersectional relationships are resolved, although phylogenetic discordance caused mainly by incomplete lineage sorting and hybridization was detected. Molecular dating and ancestral area reconstruction suggest a Northern Hemisphere high-latitude origin of Abies during the Late Cretaceous, but all extant firs diversified during the Miocene to the Pleistocene, and multiple continental and intercontinental dispersals took place in response to the late Neogene climate cooling and orogenic movements. Notably, four critically endangered firs endemic to subtropical mountains of China, including A. beshanzuensis, A. ziyuanensis, A. fanjingshanensis and A. yuanbaoshanensis from east to west, have different origins and evolutionary histories. Moreover, three hotspots of species richness, including western North America, central Japan, and the Hengduan Mountains, were identified in Abies. Elevation and precipitation, particularly precipitation of the coldest quarter, are the most significant environmental factors driving the global distribution pattern of fir species diversity. Some morphological traits are evolutionarily constrained, and those linked to elevational variation (e.g., purple cone) and cold resistance (e.g., pubescent branch and resinous bud) may have contributed to the diversification of global firs. Our study sheds new light on the spatiotemporal evolution of global firs, which will be of great help to forest management and species conservation in a warming world.

Identifiability of Level-1 Species Networks from Gene Tree Quartets

When hybridization or other forms of lateral gene transfer have occurred, evolutionary relationships of species are better represented by phylogenetic networks than by trees. While inference of such networks remains challenging, several recently proposed methods are based on quartet concordance factors-the probabilities that a tree relating a gene sampled from the species displays the possible 4-taxon relationships. Building on earlier results, we investigate what level-1 network features are identifiable from concordance factors under the network multispecies coalescent model. We obtain results on both topological features of the network, and numerical parameters, uncovering a number of failures of identifiability related to 3-cycles in the network. Addressing these identifiability issues is essential for designing statistically consistent inference methods.

Phylogenomic data resolved the deep relationships of Gymnogynoideae (Selaginellaceae)

The unresolved phylogenetic framework within the Selaginellaceae subfamily Gymnogynoideae (ca. 130 species) has hindered our comprehension of the diversification and evolution of Selaginellaceae, one of the most important lineages in land plant evolution. Here, based on plastid and nuclear data extracted from genomic sequencing of more than 90% species of all genera except two in Gymnogynoideae, a phylogenomic study focusing on the contentious relationships among the genera in Gymnogynoideae was conducted. Our major results included the following: (1) Only single-copy region (named NR) and only one ribosomal operon was firstly found in Afroselaginella among vascular plants, the plastome structure of Gymnogynoideae is diverse among the six genera, and the direct repeats (DR) type is inferred as the ancestral state in the subfamily; (2) The first strong evidence was found to support Afroselaginella as a sister to Megaloselaginella. Alternative placements of Ericetorum and Gymnogynum were detected, and their relationships were investigated by analyzing the variation of phylogenetic signals; and (3) The most likely genus-level relationships in Gymnogynoideae might be: ((Bryodesma, Lepidoselaginella), (((Megaloselaginella, Afroselaginella), Ericetorum), Gymnogynum)), which was supported by maximum likelihood phylogeny based on plastid datasets, maximum likelihood, and Bayesian inference based on SCG dataset and concatenated nuclear and plastid datasets and the highest proportion of phylogenetic signals of plastid genes.

Phylogenomics of Brachystegia: Insights into the origin of African miombo woodlands

PREMISE

Phylogenetic approaches can provide valuable insights on how and when a biome emerged and developed using its structuring species. In this context, Brachystegia Benth, a dominant genus of trees in miombo woodlands, appears as a key witness of the history of the largest woodland and savanna biome of Africa.

METHODS

We reconstructed the evolutionary history of the genus using targeted-enrichment sequencing on 60 Brachystegia specimens for a nearly complete species sampling. Phylogenomic inferences used supermatrix (RAxML-NG) and summary-method (ASTRAL-III) approaches. Conflicts between species and gene trees were assessed, and the phylogeny was time-calibrated in BEAST. Introgression between species was explored using Phylonet.

RESULTS

The phylogenies were globally congruent regardless of the method used. Most of the species were recovered as monophyletic, unlike previous plastid phylogenetic reconstructions where lineages were shared among geographically close individuals independently of species identity. Still, most of the individual gene trees had low levels of phylogenetic information and, when informative, were mostly in conflict with the reconstructed species trees. These results suggest incomplete lineage sorting and/or reticulate evolution, which was supported by network analyses. The BEAST analysis supported a Pliocene origin for current Brachystegia lineages, with most of the diversification events dated to the Pliocene-Pleistocene.

CONCLUSIONS

These results suggest a recent origin of species of the miombo, congruently with their spatial expansion documented from plastid data. Brachystegia species appear to behave potentially as a syngameon, a group of interfertile but still relatively well-delineated species, an aspect that deserves further investigations.

River Drainage Reorganization and Reticulate Evolution in the Two-Lined Salamander (Eurycea bislineata) Species Complex

The origin and eventual loss of biogeographic barriers can create alternating periods of allopatry and secondary contact, facilitating gene flow among distinct metapopulations and generating reticulate evolutionary histories that are not adequately described by a bifurcating evolutionary tree. One such example may exist in the two-lined salamander (Eurycea bislineata) species complex, where discordance among morphological and molecular datasets has created a "vexing taxonomic challenge." Previous phylogeographic analyses of mitochondrial DNA (mtDNA) suggested that the reorganization of Miocene paleodrainages drove vicariance and dispersal, but the inherent limitations of a single-locus dataset precluded the evaluation of subsequent gene flow. Here, we generate triple-enzyme restriction site-associated DNA sequencing (3RAD) data for > 100 individuals representing all major mtDNA lineages and use a suite of complementary methods to demonstrate that discordance among earlier datasets is best explained by a reticulate evolutionary history influenced by river drainage reorganization. Systematics of such groups should acknowledge these complex histories and relationships that are not strictly hierarchical. [Amphibian; hybridization; introgression; Plethodontidae; stream capture.].

Detection of Ghost Introgression Requires Exploiting Topological and Branch Length Information

In recent years, the study of hybridization and introgression has made significant progress, with ghost introgression-the transfer of genetic material from extinct or unsampled lineages to extant species-emerging as a key area for research. Accurately identifying ghost introgression, however, presents a challenge. To address this issue, we focused on simple cases involving 3 species with a known phylogenetic tree. Using mathematical analyses and simulations, we evaluated the performance of popular phylogenetic methods, including HyDe and PhyloNet/MPL, and the full-likelihood method, Bayesian Phylogenetics and Phylogeography (BPP), in detecting ghost introgression. Our findings suggest that heuristic approaches relying on site-pattern counts or gene-tree topologies struggle to differentiate ghost introgression from introgression between sampled non-sister species, frequently leading to incorrect identification of donor and recipient species. The full-likelihood method BPP uses multilocus sequence alignments directly-hence taking into account both gene-tree topologies and branch lengths, by contrast, is capable of detecting ghost introgression in phylogenomic datasets. We analyzed a real-world phylogenomic dataset of 14 species of Jaltomata (Solanaceae) to showcase the potential of full-likelihood methods for accurate inference of introgression.

The rise of baobab trees in Madagascar

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

Phylogenomics of Psammodynastes and Buhoma (Elapoidea: Serpentes), with the description of a new Asian snake family

Asian mock vipers of the genus Psammodynastes and African forest snakes of the genus Buhoma are two genera belonging to the snake superfamily Elapoidea. The phylogenetic placements of Psammodynastes and Buhoma within Elapoidea has been extremely unstable which has resulted in their uncertain and debated taxonomy. We used ultraconserved elements and traditional nuclear and mitochondrial markers to infer the phylogenetic relationships of these two genera with other elapoids. Psammodynastes, for which a reference genome has been sequenced, were found, with strong branch support, to be a relatively early diverging split within Elapoidea that is sister to a clade consisting of Elapidae, Micrelapidae and Lamprophiidae. Hence, we allocate Psammodynastes to its own family, Psammodynastidae new family. However, the phylogenetic position of Buhoma could not be resolved with a high degree of confidence. Attempts to identify the possible sources of conflict in the rapid radiation of elapoid snakes suggest that both hybridisation/introgression during the rapid diversification, including possible ghost introgression, as well as incomplete lineage sorting likely have had a confounding role. The usual practice of combining mitochondrial loci with nuclear genomic data appears to mislead phylogeny reconstructions in rapid radiation scenarios, especially in the absence of genome scale data.

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.

Ultrafast learning of four-node hybridization cycles in phylogenetic networks using algebraic invariants

Motivation

The abundance of gene flow in the Tree of Life challenges the notion that evolution can be represented with a fully bifurcating process which cannot capture important biological realities like hybridization, introgression, or horizontal gene transfer. Coalescent-based network methods are increasingly popular, yet not scalable for big data, because they need to perform a heuristic search in the space of networks as well as numerical optimization that can be NP-hard. Here, we introduce a novel method to reconstruct phylogenetic networks based on algebraic invariants. While there is a long tradition of using algebraic invariants in phylogenetics, our work is the first to define phylogenetic invariants on concordance factors (frequencies of four-taxon splits in the input gene trees) to identify level-1 phylogenetic networks under the multispecies coalescent model.

Results

Our novel hybrid detection methodology is optimization-free as it only requires the evaluation of polynomial equations, and as such, it bypasses the traversal of network space, yielding a computational speed at least 10 times faster than the fastest-to-date network methods. We illustrate our method's performance on simulated and real data from the genus Canis.

Availability and implementation

We present an open-source publicly available Julia package PhyloDiamond.jl available at https://github.com/solislemuslab/PhyloDiamond.jl with broad applicability within the evolutionary community.

IDENTIFIABILITY OF LEVEL-1 SPECIES NETWORKS FROM GENE TREE QUARTETS

When hybridization or other forms of lateral gene transfer have occurred, evolutionary relationships of species are better represented by phylogenetic networks than by trees. While inference of such networks remains challenging, several recently proposed methods are based on quartet concordance factors - the probabilities that a tree relating a gene sampled from the species displays the possible 4-taxon relationships. Building on earlier results, we investigate what level-1 network features are identifiable from concordance factors under the network multispecies coalescent model. We obtain results on both topological features of the network, and numerical parameters, uncovering a number of failures of identifiability related to 3-cycles in the network.

Phylogenomic analyses and chromosome ploidy identification reveal multiple cryptic species in Allium sikkimense complex (Amaryllidaceae)

Polyploidization is a process that typically leads to instantaneous reproductive isolation and has, therefore, been considered as one of the major evolutionary forces in the species-rich Hengduan Mountains (HM), yet this topic remains poorly studied in the region. Allium sikkimense and its relatives (about eight species) compose a natural diploid-polyploid complex with the highest diversity in the HM and adjacent areas. A combination of nuclear ribosomal DNA (nrDNA), plastome, transcriptome, and ploidy identification through chromosome counting and flow cytometry is employed to reconstruct the phylogenetic relationships in this complex and to investigate the frequency and the evolutionary significance of polyploidy in the complex. The plastome failed to resolve the phylogenetic relationships of the different species in the A. sikkimense complex, and the phylogenetic tree based on nrDNA also has limited resolution. However, our study reveals a well-resolved phylogenetic framework for species in the A. sikkimense complex using more than 1,000 orthologous genes from the transcriptome data. Previously recognized morphospecies A. sikkimense are non-monophyletic and comprise at least two independently evolved lineages (i.e., cryptic species), each forming a clade with different diploid species in this complex. The embedded pattern of octoploid A. jichouense and tetraploid A. sp. nov. within different polyploid samples of A. sikkimense supports a possible scenario of budding speciation (via niche divergence). Furthermore, our results reveal that co-occurring species in the A. sikkimense complex usually have different ploidy levels, suggesting that polyploidy is an important process for reproductive isolation of sympatric Allium species. Phylogenetic network analyses suggested that the phylogenetic relationships of the A. sikkimense complex, allowing for reticulation events, always fit the dataset better than a simple bifurcating tree. In addition, the included or exserted filaments, which have long been used to delimit species, are highly unreliable taxonomically due to their extensive parallel and convergent evolution.

Summary Tests of Introgression Are Highly Sensitive to Rate Variation Across Lineages

The evolutionary implications and frequency of hybridization and introgression are increasingly being recognized across the tree of life. To detect hybridization from multi-locus and genome-wide sequence data, a popular class of methods are based on summary statistics from subsets of 3 or 4 taxa. However, these methods often carry the assumption of a constant substitution rate across lineages and genes, which is commonly violated in many groups. In this work, we quantify the effects of rate variation on the D test (also known as ABBA-BABA test), the D3 test, and HyDe. All 3 tests are used widely across a range of taxonomic groups, in part because they are very fast to compute. We consider rate variation across species lineages, across genes, their lineage-by-gene interaction, and rate variation across gene-tree edges. We simulated species networks according to a birth-death-hybridization process, so as to capture a range of realistic species phylogenies. For all 3 methods tested, we found a marked increase in the false discovery of reticulation (type-1 error rate) when there is rate variation across species lineages. The D3 test was the most sensitive, with around 80% type-1 error, such that D3 appears to more sensitive to a departure from the clock than to the presence of reticulation. For all 3 tests, the power to detect hybridization events decreased as the number of hybridization events increased, indicating that multiple hybridization events can obscure one another if they occur within a small subset of taxa. Our study highlights the need to consider rate variation when using site-based summary statistics, and points to the advantages of methods that do not require assumptions on evolutionary rates across lineages or across genes.

Diploid species phylogeny and evolutionary reticulation indicate early radiation of Ephedra in the Tethys coast

Reconstructing a robust species phylogeny and disentangling the evolutionary and biogeographic history of the gymnosperm genus Ephedra, which has a large genome and rich polyploids, remain a big challenge. Here we reconstructed a transcriptome-based phylogeny of 19 diploid Ephedra species, and explored evolutionary reticulations in this genus represented by 50 diploid and polyploid species, using four low-copy nuclear and nine plastid genes. The diploid species phylogeny indicates that the Mediterranean species diverged first, and the remaining species split into three clades, including the American species (Clade A), E. rhytidosperma, and all other Asian species (Clade B). The single-gene trees placed E. rhytidosperma sister to Clade A, Clade B, or Clades A + B in similar proportions, suggesting that radiation and gene flow likely occurred in the early evolution of Ephedra. In addition, reticulate evolution occurred not only among the deep nodes, but also in the recently evolved South American species, which further caused difficulty in phylogenetic reconstruction. Moreover, we found that allopolyploid speciation was pervasive in Ephedra. Our study also suggests that Ephedra very likely originated in the Tethys coast during the late Cretaceous, and the South American Ephedra species have a single origin by dispersal from Mexico or North America.

The Utilization of Reference-Guided Assembly and In Silico Libraries Improves the Draft Genome of Clarias batrachus and Culter alburnus

Long-read sequencing technologies can generate highly contiguous genome assemblies compared to short-read methods. However, their higher cost often poses a significant barrier. To address this, we explore the utilization of mapping-based genome assembly and reference-guided assembly as cost-effective alternative approaches. We assess the efficacy of these approaches in improving the contiguity of Clarias batrachus and Culter alburnus draft genomes. Our findings demonstrate that employing an iterative mapping strategy leads to a reduction in assembly errors. Specifically, after three iterations, the Mismatches per 100 kbp value for the C. batrachus genome decreased from 2447.20 to 2432.67, reaching a minimum of 2422.67 after two iterations. Additionally, the N50 value for the C. batrachus genome increased from 362,143 to 1,315,126 bp, with a maximum of 1,315,403 bp after two iterations. Furthermore, we achieved Mismatches per 100 kbp values of 3.70 for the reference-guided assembly of C. batrachus and 0.34 for C. alburnus. Correspondingly, the N50 value for the C. batrachus and C. alburnus genomes increased from 362,143 bp and 3,686,385 bp to 2,026,888 bp and 43,735,735 bp, respectively. Finally, we successfully utilized the improved C. batrachus and C. alburnus genomes to compare genome studies using the combined approach of Ragout and Ragtag. Through a comprehensive comparative analysis of mapping-based and reference-guided genome assembly methods, we shed light on the specific contributions of reference-guided assembly in reducing assembly errors and improving assembly continuity and integrity. These advancements establish reference-guided assembly and the utilization of in silico libraries as a promising and suitable approach for comparative genomics studies.

Evolutionary Insights into the Relationship of Frogs, Salamanders, and Caecilians and Their Adaptive Traits, with an Emphasis on Salamander Regeneration and Longevity

The extant amphibians have developed uncanny abilities to adapt to their environment. I compared the genes of amphibians to those of other vertebrates to investigate the genetic changes underlying their unique traits, especially salamanders' regeneration and longevity. Using the well-supported Batrachia tree, I found that salamander genomes have undergone accelerated adaptive evolution, especially for development-related genes. The group-based comparison showed that several genes are under positive selection, rapid evolution, and unexpected parallel evolution with traits shared by distantly related species, such as the tail-regenerative lizard and the longer-lived naked mole rat. The genes, such as EEF1E1, PAFAH1B1, and OGFR, may be involved in salamander regeneration, as they are involved in the apoptotic process, blastema formation, and cell proliferation, respectively. The genes PCNA and SIRT1 may be involved in extending lifespan, as they are involved in DNA repair and histone modification, respectively. Some genes, such as PCNA and OGFR, have dual roles in regeneration and aging, which suggests that these two processes are interconnected. My experiment validated the time course differential expression pattern of SERPINI1 and OGFR, two genes that have evolved in parallel in salamanders and lizards during the regeneration process of salamander limbs. In addition, I found several candidate genes responsible for frogs' frequent vocalization and caecilians' degenerative vision. This study provides much-needed insights into the processes of regeneration and aging, and the discovery of the critical genes paves the way for further functional analysis, which could open up new avenues for exploiting the genetic potential of humans and improving human well-being.

Ancient hybridization in Curcuma (Zingiberaceae)-Accelerator or brake in lineage diversifications?

Hybridization is a widespread phenomenon in the evolution of plants and exploring its role is crucial to understanding diversification processes of many taxonomic groups. Recently, more attention is focused on the role of ancient hybridization that has repeatedly been shown as triggers of evolutionary radiation, although in some cases, it can prevent further diversification. The causes, frequency, and consequences of ancient hybridization remain to be explored. Here, we present an account of several events of ancient hybridization in turmeric, the economically important plant genus Curcuma (Zingiberaceae), which harbors about 130 known species. We analyzed 1094 targeted low-copy genes and plastomes obtained by next-generation sequencing of 37 species of Curcuma, representing the known genetic diversity and spanning the geographical distribution of the genus. Using phylogenetic network analysis, we show that the entire genus Curcuma as well as its most speciose lineage arose via introgression from the genus Pyrgophyllum and one of the extinct lineages, respectively. We also document a single event of ancient hybridization, with C. vamana as a product, that represents an evolutionary dead end. We further discuss distinct circumstances of those hybridization events that deal mainly with (in)congruence in chromosome counts of the parental lineages.

Phylogenomic Analyses Reveal an Allopolyploid Origin of Core Didymocarpinae (Gesneriaceae) Followed by Rapid Radiation

Allopolyploid plants have long been regarded as possessing genetic advantages under certain circumstances due to the combined effects of their hybrid origins and duplicated genomes. However, the evolutionary consequences of allopolyploidy in lineage diversification remain to be fully understood. Here, we investigate the evolutionary consequences of allopolyploidy using 138 transcriptomic sequences of Gesneriaceae, including 124 newly sequenced, focusing particularly on the largest subtribe Didymocarpinae. We estimated the phylogeny of Gesneriaceae using concatenated and coalescent-based methods based on five different nuclear matrices and 27 plastid genes, focusing on relationships among major clades. To better understand the evolutionary affinities in this family, we applied a range of approaches to characterize the extent and cause of phylogenetic incongruence. We found that extensive conflicts between nuclear and chloroplast genomes and among nuclear genes were caused by both incomplete lineage sorting (ILS) and reticulation, and we found evidence of widespread ancient hybridization and introgression. Using the most highly supported phylogenomic framework, we revealed multiple bursts of gene duplication throughout the evolutionary history of Gesneriaceae. By incorporating molecular dating and analyses of diversification dynamics, our study shows that an ancient allopolyploidization event occurred around the Oligocene-Miocene boundary, which may have driven the rapid radiation of core Didymocarpinae.

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 32-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.

Hybridization and Transgressive Evolution Generate Diversity in an Adaptive Radiation of Anolis Lizards

Interspecific hybridization may act as a major force contributing to the evolution of biodiversity. Although generally thought to reduce or constrain divergence between 2 species, hybridization can, paradoxically, promote divergence by increasing genetic variation or providing novel combinations of alleles that selection can act upon to move lineages toward new adaptive peaks. Hybridization may, then, play a key role in adaptive radiation by allowing lineages to diversify into new ecological space. Here, we test for signatures of historical hybridization in the Anolis lizards of Puerto Rico and evaluate 2 hypotheses for the role of hybridization in facilitating adaptive radiation-the hybrid swarm origins hypothesis and the syngameon hypothesis. Using whole genome sequences from all 10 species of Puerto Rican anoles, we calculated D and f-statistics (from ABBA-BABA tests) to test for introgression across the radiation and employed multispecies network coalescent methods to reconstruct phylogenetic networks that allow for hybridization. We then analyzed morphological data for these species to test for patterns consistent with transgressive evolution, a phenomenon in which the trait of a hybrid lineage is found outside of the range of its 2 parents. Our analyses uncovered strong evidence for introgression at multiple stages of the radiation, including support for an ancient hybrid origin of a clade comprising half of the extant Puerto Rican anole species. Moreover, we detected significant signals of transgressive evolution for 2 ecologically important traits, head length and toepad width, the latter of which has been described as a key innovation in Anolis. [Adaptive radiation; introgression; multispecies network coalescent; phenotypic evolution; phylogenetic network; reticulation; syngameon; transgressive segregation.].

Phylotranscriptomics supports numerous polyploidization events and phylogenetic relationships in Nicotiana

Introduction

Nicotiana L. (Solanaceae) is of great scientific and economic importance, and polyploidization has been pivotal in shaping this genus. Despite many previous studies on the Nicotiana phylogenetic relationship and hybridization, evidence from whole genome data is still lacking.

Methods

In this study, we obtained 995 low-copy genes and plastid transcript fragments from the transcriptome datasets of 26 Nicotiana species, including all sections. We reconstructed the phylogenetic relationship and phylogenetic network of diploid species.

Results

The incongruence among gene trees showed that the formation of N. sylvestris involved incomplete lineage sorting. The nuclear-plastid discordance and nuclear introgression absence indicated that organelle capture from section Trigonophyllae was involved in forming section Petunioides. Furthermore, we analyzed the evolutionary origin of polyploid species and dated the time of hybridization events based on the analysis of PhyloNet, sequence similarity search, and phylogeny of subgenome approaches. Our results highly evidenced the hybrid origins of five polyploid sections, including sections Nicotiana, Repandae, Rusticae, Polydicliae, and Suaveolentes. Notably, we provide novel insights into the hybridization event of section Polydicliae and Suaveolentes. The section Polydicliae formed from a single hybridization event between maternal progenitor N. attenuata and paternal progenitor N. undulata; the N. sylvestris (paternal progenitor) and the N. glauca (maternal progenitor) were involved in the formation of section Suaveolentes.

Discussion

This study represents the first exploration of Nicotiana polyploidization events and phylogenetic relationships using the high-throughput RNA-seq approach. It will provide guidance for further studies in molecular systematics, population genetics, and ecological adaption studies in Nicotiana and other related species.

Phylogenomics reveals widespread hybridization and polyploidization in Henckelia (Gesneriaceae)

BACKGROUND AND AIMS

Hybridization has long been recognized as an important process for plant evolution and is often accompanied by polyploidization, another prominent force in generating biodiversity. Despite its pivotal importance in evolution, the actual prevalence and distribution of hybridization across the tree of life remain unclear.

METHODS

We used whole-genome shotgun (WGS) sequencing and cytological data to investigate the evolutionary history of Henckelia, a large genus in the family Gesneriaceae with a high frequency of suspected hybridization and polyploidization events. We generated WGS sequencing data at about 10× coverage for 26 Chinese Henckelia species plus one Sri Lankan species. To untangle the hybridization history, we separately extracted whole plastomes and thousands of single-copy nuclear genes from the sequencing data, and reconstructed phylogenies based on both nuclear and plastid data. We also explored sources of both genealogical and cytonuclear conflicts and identified signals of hybridization and introgression within our phylogenomic dataset using several statistical methods. Additionally, to test the polyploidization history, we evaluated chromosome counts for 45 populations of the 27 Henckelia species studied.

KEY RESULTS

We obtained well-supported phylogenetic relationships using both concatenation- and coalescent-based methods. However, the nuclear phylogenies were highly inconsistent with the plastid phylogeny, and we observed intensive discordance among nuclear gene trees. Further analyses suggested that both incomplete lineage sorting and gene flow contributed to the observed cytonuclear and genealogical discordance. Our analyses of introgression and phylogenetic networks revealed a complex history of hybridization within the genus Henckelia. In addition, based on chromosome counts for 27 Henckelia species, we found independent polyploidization events occurred within Henckelia after different hybridization events.

CONCLUSIONS

Our findings demonstrated that hybridization and polyploidization are common in Henckelia. Furthermore, our results revealed that H. oblongifolia is not a member of the redefined Henckelia and they suggested several other taxonomic treatments in this genus.

Comparative genomics reveals the hybrid origin of a macaque group

Although species can arise through hybridization, compelling evidence for hybrid speciation has been reported only rarely in animals. Here, we present phylogenomic analyses on genomes from 12 macaque species and show that the fascicularis group originated from an ancient hybridization between the sinica and silenus groups ~3.45 to 3.56 million years ago. The X chromosomes and low-recombination regions exhibited equal contributions from each parental lineage, suggesting that they were less affected by subsequent backcrossing and hence could have played an important role in maintaining hybrid integrity. We identified many reproduction-associated genes that could have contributed to the development of the mixed sexual phenotypes characteristic of the fascicularis group. The phylogeny within the silenus group was also resolved, and functional experimentation confirmed that all extant Western silenus species are susceptible to HIV-1 infection. Our study provides novel insights into macaque evolution and reveals a hybrid speciation event that has occurred only very rarely in primates.

Integrated phylogenomic analyses unveil reticulate evolution in Parthenocissus (Vitaceae), highlighting speciation dynamics in the Himalayan-Hengduan Mountains

Hybridization caused by frequent environmental changes can lead both to species diversification (speciation) and to speciation reversal (despeciation), but the latter has rarely been demonstrated. Parthenocissus, a genus with its trifoliolate lineage in the Himalayan-Hengduan Mountains (HHM) region showing perplexing phylogenetic relationships, provides an opportunity for investigating speciation dynamics based on integrated evidence. We investigated phylogenetic discordance and reticulate evolution in Parthenocissus based on rigorous analyses of plastome and transcriptome data. We focused on reticulations in the trifoliolate lineage in the HHM region using a population-level genome resequencing dataset, incorporating evidence from morphology, distribution, and elevation. Comprehensive analyses confirmed multiple introgressions within Parthenocissus in a robust temporal-spatial framework. Around the HHM region, at least three hybridization hot spots were identified, one of which showed evidence of ongoing speciation reversal. We present a solid case study using an integrative methodological approach to investigate reticulate evolutionary history and its underlying mechanisms in plants. It demonstrates an example of speciation reversal through frequent hybridizations in the HHM region, which provides new perspectives on speciation dynamics in mountainous areas with strong topographic and environmental heterogeneity.

Gene flow assessment helps to distinguish strong genomic structure from speciation in an Iberian ant-eating spider

Although genomic data is boosting our understanding of evolution, we still lack a solid framework to perform reliable genome-based species delineation. This problem is especially critical in the case of phylogeographically structured organisms, with allopatric populations showing similar divergence patterns as species. Here, we assess the species limits and phylogeography of Zodarion alacre, an ant-eating spider widely distributed across the Iberian Peninsula. We first performed species delimitation based on genome-wide data and then validated these results using additional evidence. A commonly employed species delimitation strategy detected four distinct lineages with almost no admixture, which present allopatric distributions. These lineages showed ecological differentiation but no clear morphological differentiation, and evidence of introgression in a mitochondrial barcode. Phylogenomic networks found evidence of substantial gene flow between lineages. Finally, phylogeographic methods highlighted remarkable isolation by distance and detected evidence of range expansion from south-central Portugal to central-north Spain. We conclude that despite their deep genomic differentiation, the lineages of Z. alacre do not show evidence of complete speciation. Our results likely shed light on why Zodarion is among the most diversified spider genera despite its limited distribution and support the use of gene flow evidence to inform species boundaries.

Hybridization and introgression are prevalent in Southern European Erysimum (Brassicaceae) species

BACKGROUND AND AIMS

Hybridization is a common and important force in plant evolution. One of its outcomes is introgression - the transfer of small genomic regions from one taxon to another by hybridization and repeated backcrossing. This process is believed to be common in glacial refugia, where range expansions and contractions can lead to cycles of sympatry and isolation, creating conditions for extensive hybridization and introgression. Polyploidization is another genome-wide process with a major influence on plant evolution. Both hybridization and polyploidization can have complex effects on plant evolution. However, these effects are often difficult to understand in recently evolved species complexes.

METHODS

We combined flow cytometry, analyses of transcriptomic sequences and pollen tube growth assays to investigate the consequences of polyploidization, hybridization and introgression on the recent evolution of several Erysimum (Brassicaceae) species from the South of the Iberian Peninsula, a well-known glacial refugium. This species complex differentiated in the last 2 million years, and its evolution has been hypothesized to be determined mainly by polyploidization, interspecific hybridization and introgression.

KEY RESULTS

Our results support a scenario of widespread hybridization involving both extant and 'ghost' taxa. Several taxa studied here, most notably those with purple corollas, are polyploids, probably of allopolyploid origin. Moreover, hybridization in this group might be an ongoing phenomenon, as pre-zygotic barriers appeared weak in many cases.

CONCLUSIONS

The evolution of Erysimum spp. has been determined by hybridization to a large extent. Species with purple (polyploids) and yellow flowers (mostly diploid) exhibit a strong signature of introgression in their genomes, indicating that hybridization occurred regardless of colour and across ploidy levels. Although the adaptive value of such genomic exchanges remains unclear, our results demonstrate the significance of hybridization for plant diversification, which should be taken into account when studying plant evolution.

Bayesian inference of admixture graphs on Native American and Arctic populations

Admixture graphs are mathematical structures that describe the ancestry of populations in terms of divergence and merging (admixing) of ancestral populations as a graph. An admixture graph consists of a graph topology, branch lengths, and admixture proportions. The branch lengths and admixture proportions can be estimated using numerous numerical optimization methods, but inferring the topology involves a combinatorial search for which no polynomial algorithm is known. In this paper, we present a reversible jump MCMC algorithm for sampling high-probability admixture graphs and show that this approach works well both as a heuristic search for a single best-fitting graph and for summarizing shared features extracted from posterior samples of graphs. We apply the method to 11 Native American and Siberian populations and exploit the shared structure of high-probability graphs to characterize the relationship between Saqqaq, Inuit, Koryaks, and Athabascans. Our analyses show that the Saqqaq is not a good proxy for the previously identified gene flow from Arctic people into the Na-Dene speaking Athabascans.

Phylogeny of Saxifraga section Saxifraga subsection Arachnoideae (Saxifragaceae) and the origin of low elevation shade-dwelling species

Saxifraga section Saxifraga subsection Arachnoideae is a lineage of 12 species distributed mainly in the European Alps. It is unusual in terms of ecological diversification by containing both high elevation species from exposed alpine habitats and low elevation species from shady habitats such as overhanging rocks and cave entrances. Our aims are to explore which of these habitat types is ancestral, and to identify the possible drivers of this remarkable ecological diversification. Using a Hybseq DNA-sequencing approach and a complete species sample we reconstructed and dated the phylogeny of subsection Arachnoideae. Using Landolt indicator values, this phylogenetic tree was used for the reconstruction of the evolution of temperature, light and soil pH requirements in this lineage. Diversification of subsection Arachnoideae started in the late Pliocene and continued through the Pleistocene. Both diversification among and within clades was largely allopatric, and species from shady habitats with low light requirements are distributed in well-known refugia. We hypothesize that low light requirements evolved when species persisting in cold-stage refugia were forced into marginal habitats by more competitive warm-stage vegetation. While we do not claim that such competition resulted in speciation, it very likely resulted in adaptive evolution.

The tree of blobs of a species network: identifiability under the coalescent

Inference of species networks from genomic data under the Network Multispecies Coalescent Model is currently severely limited by heavy computational demands. It also remains unclear how complicated networks can be for consistent inference to be possible. As a step toward inferring a general species network, this work considers its tree of blobs, in which non-cut edges are contracted to nodes, so only tree-like relationships between the taxa are shown. An identifiability theorem, that most features of the unrooted tree of blobs can be determined from the distribution of gene quartet topologies, is established. This depends upon an analysis of gene quartet concordance factors under the model, together with a new combinatorial inference rule. The arguments for this theoretical result suggest a practical algorithm for tree of blobs inference, to be fully developed in a subsequent work.

Phylotranscriptomics reveals the reticulate evolutionary history of a widespread diatom species complex

In contrast to surveys based on a few genes that often provide limited taxonomic resolution, transcriptomes provide a wealth of genomic loci that can resolve relationships among taxonomically challenging lineages. Diatoms are a diverse group of aquatic microalgae that includes important bioindicator species and many such lineages. One example is Nitzschia palea, a widespread species complex with several morphologically defined taxonomic varieties, some of which are critical pollution indicators. Morphological differences among the varieties are subtle and phylogenetic studies based on a few genes fail to resolve their evolutionary relationships. We conducted morphometric and transcriptome analyses of 10 Nitzschia palea strains to resolve the relationships among strains and taxonomic varieties. Nitzschia palea was resolved into three clades, one of which corresponds to a group of strains with narrow linear-lanceolate valves. The other morphological group recovered in the shape outline analysis was not monophyletic and consisted of two clades. Gene-tree concordance analyses and phylogenetic network estimations revealed patterns of incomplete lineage sorting and gene flow between intraspecific lineages. We detected reticulated evolutionary patterns among lineages with different morphologies, resulting in a putative recent hybrid. Our study shows that phylogenomic analyses of unlinked nuclear loci, complemented with morphometrics, can resolve complex evolutionary histories of recently diverged species complexes.

Genomic resolution of cryptic species diversity in chipmunks

Discovery of cryptic species is essential to understand the process of speciation and assessing the impacts of anthropogenic stressors. Here, we used genomic data to test for cryptic species diversity within an ecologically well-known radiation of North American rodents, western chipmunks (Tamias). We assembled a de novo reference genome for a single species (Tamias minimus) combined with new and published targeted sequence-capture data for 21,551 autosomal and 493 X-linked loci sampled from 121 individuals spanning 22 species. We identified at least two cryptic lineages corresponding with an isolated subspecies of least chipmunk (T. minimus grisescens) and with a restricted subspecies of the yellow-pine chipmunk (Tamias amoenus cratericus) known only from around the extensive Craters of the Moon lava flow. Additional population-level sequence data revealed that the so-called Crater chipmunk is a distinct species that is abundant throughout the coniferous forests of southern Idaho. This cryptic lineage does not appear to be most closely related to the ecologically and phenotypically similar yellow-pine chipmunk but does show evidence for recurrent hybridization with this and other species.

Comparative transcriptome analysis reveals molecular adaptations underlying distinct immunity and inverted resting posture in bats

Understanding how natural selection shapes unique traits in mammals is a central topic in evolutionary biology. The mammalian order Chiroptera (bats) is attractive for biologists as well as the general public due to their specific traits of extraordinary immunity and inverted resting posture. However, genomic resources for bats that occupy key phylogenetic positions are not sufficient, which hinders comprehensive investigation of the molecular mechanisms underpinning the origin of specific traits in bats. Here, we sequenced the transcriptomes of five bats that are phylogenetically divergent and occupy key positions in the phylogenetic tree of bats. In combination with the available genomes of 19 bats and 21 other mammals, we built a database consisting of 10,918 one-to-one ortholog genes and reconstructed phylogenetic relationships of these mammals. We found that genes related to immunity, bone remodeling and cardiovascular system are targets of natural selection along the ancestral branch of bats. Further analyses revealed that the T cell receptor signaling pathway involved in immune adaptation is specifically enriched in bats. Moreover, molecular adaptations of bone remodeling, cardiovascular system, and balance sensing may help to explain the reverted resting posture in bats. Our study provides valuable transcriptome resources, enabling us to tentatively identify genetic changes associated with bat-specific traits. This work is among the first to advance our understanding of molecular underpinnings of inverted resting posture in bats, which could provide insight into healthcare applications such as hypertension in humans. This article is protected by copyright. All rights reserved.

Phylogeny and evolution of Cupressaceae: Updates on intergeneric relationships and new insights on ancient intergeneric hybridization

After the merger of the former Taxodiaceae and Cupressaceae s.s., currently the conifer family Cupressaceae (sensu lato) comprises seven subfamilies and 32 genera, most of which are important components of temperate and mountainous forests. With the exception of a recently published genus-level phylogeny of gymnosperms inferred from sequence analysis of 790 orthologs, previous phylogenetic studies of Cupressaceae were based mainly on morphological characters or a few molecular markers, and did not completely resolve the intergeneric relationships. In this study, we reconstructed a robust and well-resolved phylogeny of Cupressaceae represented by all 32 genera, using 1944 genes (Orthogroups) generated from transcriptome sequencing. Reticulate evolution analyses detected a possible ancient hybridization that occurred between ancestors of two subclades of Cupressoideae, including Microbiota-Platycladus-Tetraclinis (MPT) and Juniperus-Cupressus-Hesperocyparis-Callitropsis-Xanthocyparis (JCHCX), although both concatenation and coalescent trees are highly supported. Moreover, divergence time estimation and ancestral area reconstruction indicate that Cupressaceae very likely originated in Asia in the Triassic, and geographic isolation caused by continental separation drove the vicariant evolution of the two subfamilies Cupressoideae and Callitroideae in the northern and southern hemispheres, respectively. Evolutionary analyses of some morphological characters suggest that helically arranged linear-acicular leaves and imbricate bract-scale complexes represent ancestral states, and the shift from linear-acicular leaves to scale-like leaves was associated with the shift from helical to decussate arrangement. Our study sheds new light on phylogeny and evolutionary history of Cupressaceae, and strongly suggests that both dichotomous phylogenetic and reticulate evolution analyses be conducted in phylogenomic studies.

Impact of Ghost Introgression on Coalescent-based Species Tree Inference and Estimation of Divergence Time

The species studied in any evolutionary investigation generally constitute a small proportion of all the species currently existing or that have gone extinct. It is therefore likely that introgression, which is widespread across the tree of life, involves "ghosts," i.e., unsampled, unknown, or extinct lineages. However, the impact of ghost introgression on estimations of species trees has rarely been studied and is poorly understood. Here, we use mathematical analysis and simulations to examine the robustness of species tree methods based on the multispecies coalescent model to introgression from a ghost or extant lineage. We found that many results originally obtained for introgression between extant species can easily be extended to ghost introgression, such as the strongly interactive effects of incomplete lineage sorting (ILS) and introgression on the occurrence of anomalous gene trees (AGTs). The relative performance of the summary species tree method (ASTRAL) and the full-likelihood method (*BEAST) varies under different introgression scenarios, with the former being more robust to gene flow between non-sister species whereas the latter performing better under certain conditions of ghost introgression. When an outgroup ghost (defined as a lineage that diverged before the most basal species under investigation) acts as the donor of the introgressed genes, the time of root divergence among the investigated species generally was overestimated, whereas ingroup introgression, as commonly perceived, can only lead to underestimation. In many cases of ingroup introgression that may or may not involve ghost lineages, the stronger the ILS, the higher the accuracy achieved in estimating the time of root divergence, although the topology of the species tree is more prone to be biased by the effect of introgression.

Interpreting phylogenetic conflict: Hybridization in the most speciose genus of lichen-forming fungi

While advances in sequencing technologies have been invaluable for understanding evolutionary relationships, increasingly large genomic data sets may result in conflicting evolutionary signals that are often caused by biological processes, including hybridization. Hybridization has been detected in a variety of organisms, influencing evolutionary processes such as generating reproductive barriers and mixing standing genetic variation. Here, we investigate the potential role of hybridization in the diversification of the most speciose genus of lichen-forming fungi, Xanthoparmelia. As Xanthoparmelia is projected to have gone through recent, rapid diversification, this genus is particularly suitable for investigating and interpreting the origins of phylogenomic conflict. Focusing on a clade of Xanthoparmelia largely restricted to the Holarctic region, we used a genome skimming approach to generate 962 single-copy gene regions representing over 2 Mbp of the mycobiont genome. From this genome-scale dataset, we inferred evolutionary relationships using both concatenation and coalescent-based species tree approaches. We also used three independent tests for hybridization. Although different species tree reconstruction methods recovered largely consistent and well-supported trees, there was widespread incongruence among individual gene trees. Despite challenges in differentiating hybridization from ILS in situations of recent rapid radiations, our genome-wide analyses detected multiple potential hybridization events in the Holarctic clade, suggesting one possible source of trait variability in this hyperdiverse genus. This study highlights the value in using a pluralistic approach for characterizing genome-scale conflict, even in groups with well-resolved phylogenies, while highlighting current challenges in detecting the specific impacts of hybridization.

A phylotranscriptome study using silica gel-dried leaf tissues produces an updated robust phylogeny of Ranunculaceae

The utility of transcriptome data in plant phylogenetics has gained popularity in recent years. However, because RNA degrades much more easily than DNA, the logistics of obtaining fresh tissues has become a major limiting factor for widely applying this method. Here, we used Ranunculaceae to test whether silica-dried plant tissues could be used for RNA extraction and subsequent phylogenomic studies. We sequenced 27 transcriptomes, 21 from silica gel-dried (SD-samples) and six from liquid nitrogen-preserved (LN-samples) leaf tissues, and downloaded 27 additional transcriptomes from GenBank. Our results showed that although the LN-samples produced slightly better reads than the SD-samples, there were no significant differences in RNA quality and quantity, assembled contig lengths and numbers, and BUSCO comparisons between two treatments. Using these data, we conducted phylogenomic analyses, including concatenated- and coalescent-based phylogenetic reconstruction, molecular dating, coalescent simulation, phylogenetic network estimation, and whole genome duplication (WGD) inference. The resulting phylogeny was consistent with previous studies with higher resolution and statistical support. The 11 core Ranunculaceae tribes grouped into two chromosome type clades (T- and R-types), with high support. Discordance among gene trees is likely due to hybridization and introgression, ancient genetic polymorphism and incomplete lineage sorting. Our results strongly support one ancient hybridization event within the R-type clade and three WGD events in Ranunculales. Evolution of the three Ranunculaceae chromosome types is likely not directly related to WGD events. By clearly resolving the Ranunculaceae phylogeny, we demonstrated that SD-samples can be used for RNA-seq and phylotranscriptomic studies of angiosperms.

On the hybrid origin of the C2 Salsola divaricata agg. (Amaranthaceae) from C3 and C4 parental lineages

C₂ photosynthesis is characterised using recapturing photorespiratory CO₂ by RuBisCo in Kranz-like cells and is therefore physiologically intermediate between C₃ and C₄ photosynthesis. C₂ can be interpreted as an evolutionary precursor of C₄ and/or as the result of hybridisation between a C₃ and C₄ lineage. We compared the expression of photosynthetic traits among populations of the Salsola divaricata agg. (C₂ ) from humid subtropical to arid habitats on the coasts of the Canary Islands and Morocco and subjected them to salt and drought treatments. We screened for enhanced C₄ -like expression of traits related to habitat or treatment. We estimated species trees with a transcriptome dataset of Salsoleae and explored patterns of gene tree discordance. With phylogenetic networks and hybridisation analyses we tested for the hybrid origin of the Salsola divaricata agg. We observed distinct independent variation of photosynthetic traits within and among populations and no clear evidence for selection towards C₄ -like trait expression in more stressful habitats or treatments. We found reticulation and gene tree incongruence in Salsoleae supporting a putative hybrid origin of the Salsola divaricata agg. C₂ photosynthesis in the Salsola divaricata agg. combines traits inherited from its C₃ and C₄ parental lineages and seems evolutionarily stable, possibly well adapted to a wide climatic amplitude.

Classes of explicit phylogenetic networks and their biological and mathematical significance

The evolutionary relationships among organisms have traditionally been represented using rooted phylogenetic trees. However, due to reticulate processes such as hybridization or lateral gene transfer, evolution cannot always be adequately represented by a phylogenetic tree, and rooted phylogenetic networks that describe such complex processes have been introduced as a generalization of rooted phylogenetic trees. In fact, estimating rooted phylogenetic networks from genomic sequence data and analyzing their structural properties is one of the most important tasks in contemporary phylogenetics. Over the last two decades, several subclasses of rooted phylogenetic networks (characterized by certain structural constraints) have been introduced in the literature, either to model specific biological phenomena or to enable tractable mathematical and computational analyses. In the present manuscript, we provide a thorough review of these network classes, as well as provide a biological interpretation of the structural constraints underlying these networks where possible. In addition, we discuss how imposing structural constraints on the network topology can be used to address the scalability and identifiability challenges faced in the estimation of phylogenetic networks from empirical data.

RF-Net 2: fast inference of virus reassortment and hybridization networks

MOTIVATION

A phylogenetic network is a powerful model to represent entangled evolutionary histories with both divergent (speciation) and convergent (e.g. hybridization, reassortment, recombination) evolution. The standard approach to inference of hybridization networks is to (i) reconstruct rooted gene trees and (ii) leverage gene tree discordance for network inference. Recently, we introduced a method called RF-Net for accurate inference of virus reassortment and hybridization networks from input gene trees in the presence of errors commonly found in phylogenetic trees. While RF-Net demonstrated the ability to accurately infer networks with up to four reticulations from erroneous input gene trees, its application was limited by the number of reticulations it could handle in a reasonable amount of time. This limitation is particularly restrictive in the inference of the evolutionary history of segmented RNA viruses such as influenza A virus (IAV), where reassortment is one of the major mechanisms shaping the evolution of these pathogens.

RESULTS

Here, we expand the functionality of RF-Net that makes it significantly more applicable in practice. Crucially, we introduce a fast extension to RF-Net, called Fast-RF-Net, that can handle large numbers of reticulations without sacrificing accuracy. In addition, we develop automatic stopping criteria to select the appropriate number of reticulations heuristically and implement a feature for RF-Net to output error-corrected input gene trees. We then conduct a comprehensive study of the original method and its novel extensions and confirm their efficacy in practice using extensive simulation and empirical IAV evolutionary analyses.

AVAILABILITY AND IMPLEMENTATION

RF-Net 2 is available at https://github.com/flu-crew/rf-net-2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identifiability of species network topologies from genomic sequences using the logDet distance

Inference of network-like evolutionary relationships between species from genomic data must address the interwoven signals from both gene flow and incomplete lineage sorting. The heavy computational demands of standard approaches to this problem severely limit the size of datasets that may be analyzed, in both the number of species and the number of genetic loci. Here we provide a theoretical pointer to more efficient methods, by showing that logDet distances computed from genomic-scale sequences retain sufficient information to recover network relationships in the level-1 ultrametric case. This result is obtained under the Network Multispecies Coalescent model combined with a mixture of General Time-Reversible sequence evolution models across individual gene trees. It applies to both unlinked site data, such as for SNPs, and to sequence data in which many contiguous sites may have evolved on a common tree, such as concatenated gene sequences. Thus under standard stochastic models statistically justifiable inference of network relationships from sequences can be accomplished without consideration of individual genes or gene trees.

Phylotranscriptomics of Theaceae: generic-level relationships, reticulation and whole-genome duplication

BACKGROUND AND AIMS

Theaceae, with three tribes, nine genera and more than 200 species, are of great economic and ecological importance. Recent phylogenetic analyses based on plastomic data resolved the relationships among the three tribes and the intergeneric relationships within two of those tribes. However, generic-level relationships within the largest tribe, Theeae, were not fully resolved. The role of putative whole-genome duplication (WGD) events in the family and possible hybridization events among genera within Theeae also remain to be tested further.

METHODS

Transcriptomes or low-depth whole-genome sequencing of 57 species of Theaceae, as well as additional plastome sequence data, were generated. Using a dataset of low-copy nuclear genes, we reconstructed phylogenetic relationships using concatenated, species tree and phylogenetic network approaches. We further conducted molecular dating analyses and inferred possible WGD events by examining the distribution of the number of synonymous substitutions per synonymous site (Ks) for paralogues in each species. For plastid protein-coding sequences , phylogenies were reconstructed for comparison with the results obtained from analysis of the nuclear dataset.

RESULTS

Based on the 610 low-copy nuclear genes (858 606 bp in length) investigated, Stewartieae was resolved as sister to the other two tribes. Within Theeae, the Apterosperma-Laplacea clade grouped with Pyrenaria, leaving Camellia and Polyspora as sister. The estimated ages within Theaceae were largely consistent with previous studies based mainly on plastome data. Two reticulation events within Camellia and one between the common ancestor of Gordonia and Schima were found. All members of the tea family shared two WGD events, an older At-γ and a recent Ad-β; both events were also shared with the outgroups (Diapensiaceae, Pentaphylacaceae, Styracaceae and Symplocaceae).

CONCLUSIONS

Our analyses using low-copy nuclear genes improved understanding of phylogenetic relationships at the tribal and generic levels previously proposed based on plastome data, but the phylogenetic position of the Apterosperma-Laplacea clade needs more attention. There is no evidence for extensive intergeneric hybridization within Theeae or for a Theaceae-specific WGD event. Land bridges (e.g. the Bering land bridge) during the Late Oligocene may have permitted the intercontinental plant movements that facilitated the putative ancient introgression between the common ancestor of Gordonia and Schima.

Species Tree Inference Methods Intended to Deal with Incomplete Lineage Sorting Are Robust to the Presence of Paralogs

Many recent phylogenetic methods have focused on accurately inferring species trees when there is gene tree discordance due to incomplete lineage sorting (ILS). For almost all of these methods, and for phylogenetic methods in general, the data for each locus are assumed to consist of orthologous, single-copy sequences. Loci that are present in more than a single copy in any of the studied genomes are excluded from the data. These steps greatly reduce the number of loci available for analysis. The question we seek to answer in this study is: what happens if one runs such species tree inference methods on data where paralogy is present, in addition to or without ILS being present? Through simulation studies and analyses of two large biological data sets, we show that running such methods on data with paralogs can still provide accurate results. We use multiple different methods, some of which are based directly on the multispecies coalescent model, and some of which have been proven to be statistically consistent under it. We also treat the paralogous loci in multiple ways: from explicitly denoting them as paralogs, to randomly selecting one copy per species. In all cases, the inferred species trees are as accurate as equivalent analyses using single-copy orthologs. Our results have significant implications for the use of ILS-aware phylogenomic analyses, demonstrating that they do not have to be restricted to single-copy loci. This will greatly increase the amount of data that can be used for phylogenetic inference.[Gene duplication and loss; incomplete lineage sorting; multispecies coalescent; orthology; paralogy.].

Phylogenomic approaches to detecting and characterizing introgression

Phylogenomics has revealed the remarkable frequency with which introgression occurs across the tree of life. These discoveries have been enabled by the rapid growth of methods designed to detect and characterize introgression from whole-genome sequencing data. A large class of phylogenomic methods makes use of data across species to infer and characterize introgression based on expectations from the multispecies coalescent. These methods range from simple tests, such as the D-statistic, to model-based approaches for inferring phylogenetic networks. Here, we provide a detailed overview of the various signals that different modes of introgression are expected leave in the genome, and how current methods are designed to detect them. We discuss the strengths and pitfalls of these approaches and identify areas for future development, highlighting the different signals of introgression, and the power of each method to detect them. We conclude with a discussion of current challenges in inferring introgression and how they could potentially be addressed.

Sage Insights Into the Phylogeny of Salvia: Dealing With Sources of Discordance Within and Across Genomes

Next-generation sequencing technologies have facilitated new phylogenomic approaches to help clarify previously intractable relationships while simultaneously highlighting the pervasive nature of incongruence within and among genomes that can complicate definitive taxonomic conclusions. Salvia L., with ∼1,000 species, makes up nearly 15% of the species diversity in the mint family and has attracted great interest from biologists across subdisciplines. Despite the great progress that has been achieved in discerning the placement of Salvia within Lamiaceae and in clarifying its infrageneric relationships through plastid, nuclear ribosomal, and nuclear single-copy genes, the incomplete resolution has left open major questions regarding the phylogenetic relationships among and within the subgenera, as well as to what extent the infrageneric relationships differ across genomes. We expanded a previously published anchored hybrid enrichment dataset of 35 exemplars of Salvia to 179 terminals. We also reconstructed nearly complete plastomes for these samples from off-target reads. We used these data to examine the concordance and discordance among the nuclear loci and between the nuclear and plastid genomes in detail, elucidating both broad-scale and species-level relationships within Salvia. We found that despite the widespread gene tree discordance, nuclear phylogenies reconstructed using concatenated, coalescent, and network-based approaches recover a common backbone topology. Moreover, all subgenera, except for Audibertia, are strongly supported as monophyletic in all analyses. The plastome genealogy is largely resolved and is congruent with the nuclear backbone. However, multiple analyses suggest that incomplete lineage sorting does not fully explain the gene tree discordance. Instead, horizontal gene flow has been important in both the deep and more recent history of Salvia. Our results provide a robust species tree of Salvia across phylogenetic scales and genomes. Future comparative analyses in the genus will need to account for the impacts of hybridization/introgression and incomplete lineage sorting in topology and divergence time estimation.

Multispecies Coalescent: Theory and Applications in Phylogenetics

Species tree estimation is a basic part of many biological research projects, ranging from answering basic evolutionary questions (e.g., how did a group of species adapt to their environments?) to addressing questions in functional biology. Yet, species tree estimation is very challenging, due to processes such as incomplete lineage sorting, gene duplication and loss, horizontal gene transfer, and hybridization, which can make gene trees differ from each other and from the overall evolutionary history of the species. Over the last 10–20 years, there has been tremendous growth in methods and mathematical theory for estimating species trees and phylogenetic networks, and some of these methods are now in wide use. In this survey, we provide an overview of the current state of the art, identify the limitations of existing methods and theory, and propose additional research problems and directions.

Natural Hybrid Origin of the Controversial "Species" Clematis × pinnata (Ranunculaceae) Based on Multidisciplinary Evidence

Interspecific hybridization is common and has often been viewed as a driving force of plant diversity. However, it raises taxonomic problems and thus impacts biodiversity estimation and biological conservation. Although previous molecular phylogenetic studies suggested that interspecific hybridization may be rather common in Clematis, and artificial hybridization has been widely applied to produce new Clematis cultivars for nearly two centuries, the issue of natural hybridization of Clematis has never been addressed in detail. In this study, we tested the hybrid origin of a mesophytic and cold-adapted vine species, Clematis pinnata, which is a rare and taxonomically controversial taxon endemic to northern China. Using field investigations, flow cytometry (FCM), phylogenomic analysis, morphological statistics, and niche modeling, we tested hybrid origin and species status of C. pinnata. The FCM results showed that all the tested species were homoploid (2n = 16). Phylonet and HyDe analyses based on transcriptome data showed the hybrid origins of C. × pinnata from either C. brevicaudata × C. heracleifolia or C. brevicaudata × C. tubulosa. The plastome phylogeny depicted that C. × pinnata in different sampling sites originated by different hybridization events. Morphological analysis showed intermediacy of C. × pinnata between its putative parental species in many qualitative and quantitative characters. Niche modeling results suggested that C. × pinnata had not been adapted to a novel ecological niche independent of its putative parents. These findings demonstrated that plants of C. × pinnata did not formed a self-evolved clade and should not be treated as a species. The present study also suggests that interspecific hybridization is a common mechanism in Clematis to generate diversity and variation, and it may play an important role in the evolution and diversification of this genus. Our study implies that morphological diversity caused by natural hybridization may overstate the real species diversity in Clematis.

Maximum Parsimony Inference of Phylogenetic Networks in the Presence of Polyploid Complexes

Phylogenetic networks provide a powerful framework for modeling and analyzing reticulate evolutionary histories. While polyploidy has been shown to be prevalent not only in plants but also in other groups of eukaryotic species, most work done thus far on phylogenetic network inference assumes diploid hybridization. These inference methods have been applied, with varying degrees of success, to data sets with polyploid species, even though polyploidy violates the mathematical assumptions underlying these methods. Statistical methods were developed recently for handling specific types of polyploids and so were parsimony methods that could handle polyploidy more generally yet while excluding processes such as incomplete lineage sorting. In this article, we introduce a new method for inferring most parsimonious phylogenetic networks on data that include polyploid species. Taking gene tree topologies as input, the method seeks a phylogenetic network that minimizes deep coalescences while accounting for polyploidy. We demonstrate the performance of the method on both simulated and biological data. The inference method as well as a method for evaluating evolutionary hypotheses in the form of phylogenetic networks are implemented and publicly available in the PhyloNet software package. [Incomplete lineage sorting; minimizing deep coalescences; multilabeled trees; multispecies network coalescent; phylogenetic networks; polyploidy.]