Gene Tree Conflict Research Articles

Repeated shifts out of tropical climates preceded by whole genome duplication.

While flowering plants have diversified in virtually every terrestrial clime, climate constrains the distribution of individual lineages. Overcoming climatic constraints may be associated with diverse evolutionary phenomena including whole genome duplication (WGD), gene-tree conflict, and life-history changes. Climatic shifts may also have facilitated increases in flowering plant diversification rates. We investigate climatic shifts in the flowering plant order Ericales, which consists of c. 14 000 species with diverse climatic tolerances. We estimate phylogenetic trees from transcriptomic data, 64 chloroplast loci, and Angiosperms353 nuclear loci that, respectively, incorporate 147, 4508, and 2870 Ericales species. We use these phylogenetic trees to analyse how climatic shifts are associated with WGD, gene-tree conflict, life-history, and diversification rates. Early branches in the phylogenetic trees are extremely short, and have high levels of gene-tree conflict and at least one WGD. On lineages descended from these early branches, there is a significant association between climatic shifts (primarily out of tropical climates), further WGDs, and life-history. Extremely short early branches, and their associated gene-tree conflict and WGDs, appear to underpin the explosive origin of numerous species rich Ericales clades. The evolution of diverse climatic tolerances in these species rich clades is tightly associated with WGD and life-history.

Phylogenomics, reticulation, and biogeographical history of Elaeagnaceae

The angiosperm family Elaeagnaceae comprises three genera and ca. 100 species distributed mainly in Eurasia and North America. Little family-wide phylogenetic and biogeographic research on Elaeagnaceae has been conducted, limiting the application and preservation of natural genetic resources. Here, we reconstructed a strongly supported phylogenetic framework of Elaeagnaceae to better understand inter- and intrageneric relationships, as well as the origin and biogeographical history of the family. For this purpose, we used both nuclear and plastid sequences from Hyb-Seq and genome skimming approaches to reconstruct a well-supported phylogeny and, along with current distributional data, infer historical biogeographical processes. Our phylogenetic analyses of both nuclear and plastid data strongly support the monophyly of Elaeagnaceae and each of the three genera. Elaeagnus was resolved as sister to the well-supported clade of Hippophae and Shepherdia. The intrageneric relationships of Elaeagnus and Hippophae were also well resolved. High levels of nuclear gene tree conflict and cytonuclear discordance were detected within Elaeagnus, and our analyses suggest putative ancient and recent hybridization. We inferred that Elaeagnaceae originated at ca. 90.48 Ma (95% CI = 89.91–91.05 Ma), and long-distance dispersal likely played a major role in shaping its intercontinentally disjunct distribution. This work presents the most comprehensive phylogenetic framework for Elaeagnaceae to date, offers new insights into previously unresolved relationships in Elaeagnus, and provides a foundation for further studies on classification, evolution, biogeography, and conservation of Elaeagnaceae.

A new nuclear phylogeny of the tea family (Theaceae) unravels rapid radiations in genus Camellia

Molecular analyses of rapidly radiating groups often reveal incongruence between gene trees. This mainly results from incomplete lineage sorting, introgression, and gene tree estimation error, which complicate the estimation of phylogenetic relationships. In this study, we reconstruct the phylogeny of Theaceae using 348 nuclear loci from 68 individuals and two outgroup taxa. Sequence data were obtained by target enrichment using the recently released Angiosperm 353 universal probe set applied to herbarium specimens. The robustness of the topologies to variation in data quality was established under a range of different filtering schemes, using both coalescent and concatenation approaches. Our results confirmed most of the previously hypothesized relationships among tribes and genera, while clarifying additional interspecific relationships within the rapidly radiating genus Camellia. We recovered a remarkably high degree of gene tree heterogeneity indicative of rapid radiation in the group and observed cytonuclear conflicts, especially within Camellia. This was especially pronounced around short branches, which we primarily associate with gene tree estimation error. Our analysis also indicates that incomplete lineage sorting (ILS) contributed to gene-tree conflicts and accounted for approximately 14 % of the explained variation, whereas inferred introgression levels were low. Our study advances the understanding of the evolution of this important plant family and provides guidance on the application of target capture methods and the evaluation of key processes that influence phylogenetic discordances.

Phylogenomics and the rise of the angiosperms

Angiosperms are the cornerstone of most terrestrial ecosystems and human livelihoods1,2. A robust understanding of angiosperm evolution is required to explain their rise to ecological dominance. So far, the angiosperm tree of life has been determined primarily by means of analyses of the plastid genome3,4. Many studies have drawn on this foundational work, such as classification and first insights into angiosperm diversification since their Mesozoic origins5–7. However, the limited and biased sampling of both taxa and genomes undermines confidence in the tree and its implications. Here, we build the tree of life for almost 8,000 (about 60%) angiosperm genera using a standardized set of 353 nuclear genes8. This 15-fold increase in genus-level sampling relative to comparable nuclear studies9 provides a critical test of earlier results and brings notable change to key groups, especially in rosids, while substantiating many previously predicted relationships. Scaling this tree to time using 200 fossils, we discovered that early angiosperm evolution was characterized by high gene tree conflict and explosive diversification, giving rise to more than 80% of extant angiosperm orders. Steady diversification ensued through the remaining Mesozoic Era until rates resurged in the Cenozoic Era, concurrent with decreasing global temperatures and tightly linked with gene tree conflict. Taken together, our extensive sampling combined with advanced phylogenomic methods shows the deep history and full complexity in the evolution of a megadiverse clade.

Distinct hybridization modes in wide- and narrow-ranged lineages of Causonis (Vitaceae)

BackgroundExplaining contrasting patterns of distribution between related species is crucial for understanding the dynamics of biodiversity. Despite instances where hybridization and whole genome duplication (WGD) can yield detrimental outcomes, a role in facilitating the expansion of distribution range has been proposed. The Vitaceae genus Causonis exhibits great variations in species’ distribution ranges, with most species in the derived lineages having a much wider range than those in the early-diverged lineages. Hybridization and WGD events have been suggested to occur in Causonis based on evidence of phylogenetic discordance. The genus, therefore, provides us with an opportunity to for explore different hybridization and polyploidization modes in lineages with contrasting species’ distribution ranges. However, the evolutionary history of Causonis incorporating potential hybridization and WGD events remains to be explored.ResultsWith plastid and nuclear data from dense sampling, this study resolved the phylogenetic relationships within Causonis and revealed significant cyto-nuclear discordance. Nuclear gene tree conflicts were detected across the genus, especially in the japonica-corniculata clade, which were mainly attributed to gene flow. This study also inferred the allopolyploid origin of the core Causonis species, which promoted the accumulation of stress-related genes. Causonis was estimated to have originated in continental Asia in the early Eocene, and experienced glaciation in the early Oligocene, shortly after the divergence of the early-divergent lineages. The japonica-corniculata clade mainly diversified in the Miocene, followed by temperature declines that may have facilitated secondary contact. Species distribution modeling based on current climate change predicted that the widespread C. japonica tends to be more invasive, while the endemic C. ciliifera may be at risk of extinction.ConclusionsThis study presents Causonis, a genus with complex reticulate evolutionary history, as a model of how hybridization and WGD modes differ in lineages of contrasting species’ geographic ranges. It is important to consider specific evolutionary histories and genetic properties of the focal species within conservation strategies.

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

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

Phylogenomic analyses of Camellia support reticulate evolution among major clades

Camellia (Theaceae) is a morphologically highly diverse genus of flowering plants and includes many famous species with high economic value, and the phylogeny of this genus is not fully resolved. We used 95 transcriptomes from 87 Camellia species and identified 1481 low-copy genes to conduct a detailed analysis of the phylogeny of this genus according to various data-screening criteria. The results show that, very different from the two existing classification systems of Camellia, 87 species are grouped into 8 main clades and two independent species, and that all 8 clades except Clade 8 were strongly supported by almost all the coalescent or concatenated trees using different gene subsets. However, the relationships among these clades were weakly supported and different from analyses using different gene subsets; furthermore, they do not agree with the phylogeny from chloroplast genomes of Camellia. Additional analyses support reticulate evolution (probably resulting from introgression or hybridization) among some major Camellia lineages, providing explanation for extensive gene tree conflicts. Furthermore, we inferred that together with the formation of East Asian subtropical evergreen broad-leaved forests, Camellia underwent a radiative divergence of major clades at 23 ∼ 19 Ma in the late Miocene then had a subsequent species burst at 10 ∼ 5 Ma. Principal component and cluster analyses provides new insights into morphological changes underlying the evolution of Camellia and a reference to further clarify subgenus and sections of this genus. The comprehensive study here including a nuclear phylogeny and other analyses reveal the rapid evolutionary history of Camellia.

Evolutionary relationships of mitogenomes in a recently radiated Old World avian family

Environmentally heterogeneous mountains provide opportunities for rapid diversification and speciation. The family Prunellidae (accentors) is a group of birds comprising primarily mountain specialists that have recently radiated across the Palearctic region. This rapid diversification poses challenges to resolving their phylogeny. Herein we sequenced the complete mitogenomes and estimated the phylogeny using all 12 (including 28 individuals) currently recognized species of Prunellidae. We reconstructed the mitochondrial genome phylogeny using 13 protein-coding genes of 12 species and 2 Eurasian Tree Sparrows (Passer montanus). Phylogenetic relationships were estimated using a suite of analyses: maximum likelihood, maximum parsimony and the coalescent-based SVDquartets. Divergence times were estimated by implementing a Bayesian relaxed clock model in BEAST2. Based on the BEAST time-calibrated tree, we implemented an ancestral area reconstruction using RASP v.4.3. Our phylogenies based on the maximum likelihood, maximum parsimony and SVDquartets approaches support a clade of large-sized accentors (subgenus Laiscopus) to be sister to all other accentors with small size (subgenus Prunella). In addition, the trees also support the sister relationship of P. immaculata and P. rubeculoides ​+ ​P.atrogularis with 100% bootstrap support, but the relationships among the remaining eight species in the Prunella clade are poorly resolved. These species cluster in different positions in the three phylogenetic trees and the nodes are often poorly supported. The five nodes separating the seven species diverged simultaneously within less than half million years (i.e., between 2.71 and 3.15 million years ago), suggesting that the recent radiation is likely responsible for rampant incomplete lineage sorting and gene tree conflicts. Ancestral area reconstruction indicates a central Palearctic region origin for Prunellidae. Our study highlights that whole mitochondrial genome phylogeny can resolve major lineages within Prunellidae but is not sufficient to fully resolve the relationship among the species in the Prunella clade that almost simultaneously diversify during a short time period. Our results emphasize the challenge to reconstruct reliable phylogenetic relationship in a group of recently radiated species.

Target capture data resolve recalcitrant relationships in the coffee family (Rubioideae, Rubiaceae).

Subfamily Rubioideae is the largest of the main lineages in the coffee family (Rubiaceae), with over 8,000 species and 29 tribes. Phylogenetic relationships among tribes and other major clades within this group of plants are still only partly resolved despite considerable efforts. While previous studies have mainly utilized data from the organellar genomes and nuclear ribosomal DNA, we here use a large number of low-copy nuclear genes obtained via a target capture approach to infer phylogenetic relationships within Rubioideae. We included 101 Rubioideae species representing all but two (the monogeneric tribes Foonchewieae and Aitchinsonieae) of the currently recognized tribes, and all but one non-monogeneric tribe were represented by more than one genus. Using data from the 353 genes targeted with the universal Angiosperms353 probe set we investigated the impact of data type, analytical approach, and potential paralogs on phylogenetic reconstruction. We inferred a robust phylogenetic hypothesis of Rubioideae with the vast majority (or all) nodes being highly supported across all analyses and datasets and few incongruences between the inferred topologies. The results were similar to those of previous studies but novel relationships were also identified. We found that supercontigs [coding sequence (CDS) + non-coding sequence] clearly outperformed CDS data in levels of support and gene tree congruence. The full datasets (353 genes) outperformed the datasets with potentially paralogous genes removed (186 genes) in levels of support but increased gene tree incongruence slightly. The pattern of gene tree conflict at short internal branches were often consistent with high levels of incomplete lineage sorting (ILS) due to rapid speciation in the group. While concatenation- and coalescence-based trees mainly agreed, the observed phylogenetic discordance between the two approaches may be best explained by their differences in accounting for ILS. The use of target capture data greatly improved our confidence and understanding of the Rubioideae phylogeny, highlighted by the increased support for previously uncertain relationships and the increased possibility to explore sources of underlying phylogenetic discordance.

Phylotranscriptomics of liverworts: revisiting the backbone phylogeny and ancestral gene duplications.

With some 7300 extant species, liverworts (Marchantiophyta) represent one of the major land plant lineages. The backbone relationships, such as the phylogenetic position of Ptilidiales, and the occurrence and timing of whole-genome duplications, are still contentious. Based on analyses of the newly generated transcriptome data for 38 liverworts and complemented with those publicly available, we reconstructed the evolutionary history of liverworts and inferred gene duplication events along the 55 taxon liverwort species tree. Our phylogenomic study provided an ordinal-level liverwort nuclear phylogeny and identified extensive gene tree conflicts and cyto-nuclear incongruences. Gene duplication analyses based on integrated phylogenomics and Ks distributions indicated no evidence of whole-genome duplication events along the backbone phylogeny of liverworts. With a broadened sampling of liverwort transcriptomes, we re-evaluated the backbone phylogeny of liverworts, and provided evidence for ancient hybridizations followed by incomplete lineage sorting that shaped the deep evolutionary history of liverworts. The lack of whole-genome duplication during the deep evolution of liverworts indicates that liverworts might represent one of the few major embryophyte lineages whose evolution was not driven by whole-genome duplications.

Dilemmas in generic delimitation of Senegalia and allies (Caesalpinioideae, mimosoid clade): how to reconcile phylogenomic evidence with morphology and taxonomy?

Senegalia comprises 219 species distributed in tropical and subtropical regions of North and South America, Africa, Asia and Australia. Two sections are currently recognised within Senegalia and these are most readily distinguished by the differences in disposition of their cauline prickles, i.e. sect. Senegalia with prickles at or near leaf nodes and sect. Monacanthea with mostly internodal prickles. Previous phylogenetic studies, based primarily on small numbers of plastid DNA loci, found Senegalia to be monophyletic with two large subclades corresponding to the sections. Here, we present new phylogenomic evidence from 997 single-copy nuclear gene sequences for a small, but representative set of species. These new analyses show that Senegalia is non-monophyletic, but instead, forms a grade that is paraphyletic with respect to the remainder of the ingoid clade (i.e. Ingeae + Acacia s.s. + Acaciella), comprising two well-supported subclades most likely representing the same clades as found in previous phylogenetic studies of the genus and, interspersed between these, a third, moderately supported clade, comprising the genera Mariosousa, Pseudosenegalia and Parasenegalia. In marked contrast to the nuclear phylogeny, the two Senegalia clades are sister groups in the plastid phylogeny, based on analyses of 72 chloroplast genes, rendering the genus monophyletic, based on plastid data alone. We discuss this new evidence that Senegalia is non-monophyletic in relation to the marked cytonuclear discordance, high gene tree conflict and lack of resolution across this senegalioid grade and review the consistency of the key morphological characters distinguishing the two sections of Senegalia. We conclude that it is likely that Senegalia will need to be split into two (or possibly more) genera: a re-circumscribed Senegalia s.s. that corresponds to the existing Senegaliasect.Senegalia plus the S.ataxacantha group (Senegaliasect.Monacanthea s.s.; future studies may show that this group warrants generic status) and a new genus corresponding to the remainder of sect. Monacanthea (here designated as Senegaliasect.Monacanthea p.p.). However, re-delimiting Senegalia now would be premature given that the key morphological characters are not fully congruent with the two sections and pending denser phylogenetic sampling of taxa. A judiciously selected list of critical taxa is presented to facilitate future phylogenomic studies. Finally, we discuss the identity of Albizialeonardii, which is also placed in this senegalioid grade in these new phylogenomic analyses and place it in synonymy with Parasenegaliavogeliana.

Phylogenomic analysis of 997 nuclear genes reveals the need for extensive generic re-delimitation in Caesalpinioideae (Leguminosae).

Subfamily Caesalpinioideae with ca. 4,600 species in 152 genera is the second-largest subfamily of legumes (Leguminosae) and forms an ecologically and economically important group of trees, shrubs and lianas with a pantropical distribution. Despite major advances in the last few decades towards aligning genera with clades across Caesalpinioideae, generic delimitation remains in a state of considerable flux, especially across the mimosoid clade. We test the monophyly of genera across Caesalpinioideae via phylogenomic analysis of 997 nuclear genes sequenced via targeted enrichment (Hybseq) for 420 species and 147 of the 152 genera currently recognised in the subfamily. We show that 22 genera are non-monophyletic or nested in other genera and that non-monophyly is concentrated in the mimosoid clade where ca. 25% of the 90 genera are found to be non-monophyletic. We suggest two main reasons for this pervasive generic non-monophyly: (i) extensive morphological homoplasy that we document here for a handful of important traits and, particularly, the repeated evolution of distinctive fruit types that were historically emphasised in delimiting genera and (ii) this is an artefact of the lack of pantropical taxonomic syntheses and sampling in previous phylogenies and the consequent failure to identify clades that span the Old World and New World or conversely amphi-Atlantic genera that are non-monophyletic, both of which are critical for delimiting genera across this large pantropical clade. Finally, we discuss taxon delimitation in the phylogenomic era and especially how assessing patterns of gene tree conflict can provide additional insights into generic delimitation. This new phylogenomic framework provides the foundations for a series of papers reclassifying genera that are presented here in Advances in Legume Systematics (ALS) 14 Part 1, for establishing a new higher-level phylogenetic tribal and clade-based classification of Caesalpinioideae that is the focus of ALS14 Part 2 and for downstream analyses of evolutionary diversification and biogeography of this important group of legumes which are presented elsewhere.

Categorical edge-based analyses of phylogenomic data reveal conflicting signals for difficult relationships in the avian tree

Phylogenetic analyses fail to yield a satisfactory resolution of some relationships in the tree of life even with genome-scale datasets, so the failure is unlikely to reflect limitations in the amount of data. Gene tree conflicts are particularly notable in studies focused on these contentious nodes, and taxon sampling, different analytical methods, and/or data type effects can further confound analyses. Although many efforts have been made to incorporate biological conflicts, few studies have curated individual genes for their efficiency in phylogenomic studies. Here, we conduct an edge-based analysis of Neoavian evolution, examining the phylogenetic efficacy of two recent phylogenomic bird datasets and three datatypes (ultraconserved elements [UCEs], introns, and coding regions). We assess the potential causes for biases in signal-resolution for three difficult nodes: the earliest divergence of Neoaves, the position of the enigmatic Hoatzin (Opisthocomus hoazin), and the position of owls (Strigiformes). We observed extensive conflict among genes for all data types and datasets even after meticulous curation. Edge-based analyses (EBA) increased congruence and provided information about the impact of data type, GC content variation (GCCV), and outlier genes on each of nodes we examined. First, outlier gene signals appeared to drive different patterns of support for the relationships among the earliest diverging Neoaves. Second, the placement of Hoatzin was highly variable, although our EBA did reveal a previously unappreciated data type effect with an impact on its position. It also revealed that the resolution with the most support here was Hoatzin + shorebirds. Finally, GCCV, rather than data type (i.e., coding vs non-coding) per se, was correlated with a signal that supports monophyly of owls + Accipitriformes (hawks, eagles, and vultures). Eliminating high GCCV loci increased the signal for owls + mousebirds. Categorical EBA was able to reveal the nature of each edge and provide a way to highlight especially problematic branches that warrant a further examination. The current study increases our understanding about the contentious parts of the avian tree, which show even greater conflicts than appreciated previously.

Phylogenomic discordance suggests polytomies along the backbone of the large genus Solanum.

PremiseEvolutionary studies require solid phylogenetic frameworks, but increased volumes of phylogenomic data have revealed incongruent topologies among gene trees in many organisms both between and within genomes. Some of these incongruences indicate polytomies that may remain impossible to resolve. Here we investigate the degree of gene‐tree discordance in Solanum, one of the largest flowering plant genera that includes the cultivated potato, tomato, and eggplant, as well as 24 minor crop plants.MethodsA densely sampled species‐level phylogeny of Solanum is built using unpublished and publicly available Sanger sequences comprising 60% of all accepted species (742 spp.) and nine regions (ITS, waxy, and seven plastid markers). The robustness of this topology is tested by examining a full plastome dataset with 140 species and a nuclear target‐capture dataset with 39 species of Solanum (Angiosperms353 probe set).ResultsWhile the taxonomic framework of Solanum remained stable, gene tree conflicts and discordance between phylogenetic trees generated from the target‐capture and plastome datasets were observed. The latter correspond to regions with short internodal branches, and network analysis and polytomy tests suggest the backbone is composed of three polytomies found at different evolutionary depths. The strongest area of discordance, near the crown node of Solanum, could potentially represent a hard polytomy.ConclusionsWe argue that incomplete lineage sorting due to rapid diversification is the most likely cause for these polytomies, and that embracing the uncertainty that underlies them is crucial to understand the evolution of large and rapidly radiating lineages.

Increased resolution in the face of conflict: phylogenomics of the Neotropical bellflowers (Campanulaceae: Lobelioideae), a rapid plant radiation.

The centropogonid clade (Lobelioideae: Campanulaceae) is an Andean-centred rapid radiation characterized by repeated convergent evolution of morphological traits, including fruit type and pollination syndromes. While previous studies have resolved relationships of lineages with fleshy fruits into subclades, relationships among capsular species remain unresolved. This lack of resolution has impeded reclassification of non-monophyletic genera, whose current taxonomy relies heavily on traits that have undergone convergent evolution. Targeted sequence capture using a probe-set recently developed for the centropogonid clade was used to obtain phylogenomic data from DNA extracted from both silica-dried and herbarium leaf tissue. These data were used to infer relationships among species using concatenated and partitioned species tree methods, and to quantify gene tree discordance. While silica-dried leaf tissue resulted in longer assembled sequence data, the inclusion of herbarium samples improved taxonomic representation. Relationships among baccate lineages are similar to those inferred in previous studies, although they differ for lineages within and among capsular clades. We improve the phylogenetic resolution of Siphocampylus, which forms ten groups of closely related species which we informally name. Two subclades of Siphocampylus and two individual species are rogue taxa whose placement differs widely across analyses. Gene tree discordance (including cytonuclear discordance) is rampant. This first phylogenomic study of the centropogonid clade considerably improves our understanding of relationships in this rapid radiation. Differences across analyses and the possibility of additional lineage discoveries still hamper a solid and stable reclassification. Rapid morphological innovation corresponds with a high degree of phylogenomic complexity, including cytonuclear discordance, nuclear gene tree conflict and well-supported differences between analyses based on different nuclear loci. Together, these results point to a potential role of hemiplasy underlying repeated convergent evolution. This hallmark of rapid radiations is probably present in many other species-rich Andean plant radiations.

Dense sampling of taxa and characters improves phylogenetic resolution among deltocephaline leafhoppers (Hemiptera: Cicadellidae: Deltocephalinae)

AbstractWe analysed a phylogenomic dataset comprising 730 terminal taxa and >160,000 nucleotide positions obtained using anchored hybrid enrichment of genomic DNA for a sample of deltocephaline leafhoppers and outgroups. Maximum likelihood analyses of concatenated nucleotide and amino acid sequences as well as coalescent gene tree analysis, yielded well‐resolved phylogenetic estimates that were highly congruent with most branches receiving maximum support. Some topological incongruence occurs among the trees resulting from different analyses, mainly distributed among very short branches at intermediate levels in the phylogeny, pertaining to relationships among some tribes and multi‐tribe lineages restricted to particular continents. Coalescent gene tree analysis revealed extensive gene tree conflict at these nodes, suggesting that certain relationships may remain difficult to resolve consistently even with genome‐scale datasets and dense taxon sampling. Ancestral character state reconstruction of feeding preference indicates that grass specialization has been acquired three or fewer times in some highly diverse but relatively derived lineages. Molecular divergence time estimation suggests that the earliest divergences in the subfamily occurred during the Cretaceous but that most modern tribes did not appear until after the Cretaceous–Palaeogene boundary and that grass specialist lineages generally diversified and spread at the same time as grasslands were becoming widespread globally. These analyses also show a high level of global biogeographic structure, with several large lineages of primarily arboreal deltocephalines restricted to particular regions or continents and long‐distance dispersal among continents occurring primarily among grass‐specialist lineages. The results represent a major improvement over previous analyses of this diverse subfamily, which were based on 152 or fewer taxa and data from morphology and partial sequences of two genes.

Calophyllaceae plastomes, their structure and insights in relationships within the clusioids

A complete chloroplast genome is not yet available for numerous species of plants. Among the groups that lack plastome information is the clusioid clade (Malpighiales), which includes five families: Bonnetiaceae, Calophyllaceae, Clusiaceae, Hypericaceae, and Podostemaceae. With around 2200 species, it has few published plastomes and most of them are from Podostemaceae. Here we assembled and compared six plastomes from members of the clusioids: five from Calophyllaceae (newly sequenced) and one from Clusiaceae. Putative regions for evolutionary studies were identified and the newly assembled chloroplasts were analyzed with other available chloroplasts for the group, focusing on Calophyllaceae. Our results mostly agree with recent studies which found a general conserved structure, except for the two Podostemaceae species that have a large inversion (trnK-UUU–rbcL) and lack one intron from ycf3. Within Calophyllaceae we observed a longer LSC and reduced IRs in Mahurea exstipulata, resulting in some genic rearrangement, and a short inversion (psbJ–psbE) in Kielmeyera coriacea. Phylogenetic analyses recovered the clusioids and the five families as monophyletic and revealed that conflicts in relationships reported in the literature for the group agree with nodes concentrating uninformative or conflicting gene trees. Our study brings new insights about clusioid plastome architecture and its evolution.

Gene duplications and phylogenomic conflict underlie major pulses of phenotypic evolution in gymnosperms.

Inferring the intrinsic and extrinsic drivers of species diversification and phenotypic disparity across the tree of life is a major challenge in evolutionary biology. In green plants, polyploidy (or whole-genome duplication, WGD) is known to play a major role in microevolution and speciation, but the extent to which WGD has shaped macroevolutionary patterns of diversification and phenotypic innovation across plant phylogeny remains an open question. Here, we examine the relationship of various facets of genomic evolution-including gene and genome duplication, genome size, and chromosome number-with macroevolutionary patterns of phenotypic innovation, species diversification, and climatic occupancy in gymnosperms. We show that genomic changes, such as WGD and genome-size shifts, underlie the origins of most major extant gymnosperm clades, and notably, our results support an ancestral WGD in the gymnosperm lineage. Spikes of gene duplication typically coincide with major spikes of phenotypic innovation, while increased rates of phenotypic evolution are typically found at nodes with high gene-tree conflict, representing historic population-level dynamics during speciation. Most shifts in gymnosperm diversification since the rise of angiosperms are decoupled from putative WGDs and instead are associated with increased rates of climatic occupancy evolution, particularly in cooler and/or more arid climatic conditions, suggesting that ecological opportunity, especially in the later Cenozoic, and environmental heterogeneity have driven a resurgence of gymnosperm diversification. Our study provides critical insight on the processes underlying diversification and phenotypic evolution in gymnosperms, with important broader implications for the major drivers of both micro- and macroevolution in plants.

A phylogenomic perspective on gene tree conflict and character evolution in Caprifoliaceae using target enrichment data, with Zabelioideae recognized as a new subfamily

AbstractThe use of diverse data sets in phylogenetic studies aiming for understanding evolutionary histories of species can yield conflicting inference. Phylogenetic conflicts observed in animal and plant systems have often been explained by hybridization, incomplete lineage sorting (ILS), or horizontal gene transfer. Here, we used target enrichment data, species tree, and species network approaches to infer the backbone phylogeny of the family Caprifoliaceae, while distinguishing among sources of incongruence. We used 713 nuclear loci and 46 complete plastome sequence data from 43 samples representing 38 species from all major clades to reconstruct the phylogeny of the family using concatenation and coalescence approaches. We found significant nuclear gene tree conflict as well as cytonuclear discordance. Additionally, coalescent simulations and phylogenetic species network analyses suggested putative ancient hybridization among subfamilies of Caprifoliaceae, which seems to be the main source of phylogenetic discordance. Ancestral state reconstruction of six morphological characters revealed some homoplasy for each character examined. By dating the branching events, we inferred the origin of Caprifoliaceae at approximately 66.65 Ma in the late Cretaceous. By integrating evidence from molecular phylogeny, divergence times, and morphology, we here recognize Zabelioideae as a new subfamily in Caprifoliaceae. This work shows the necessity of using a combination of multiple approaches to identify the sources of gene tree discordance. Our study also highlights the importance of using data from both nuclear and plastid genomes to reconstruct deep and shallow phylogenies of plants.

Concordance-Based Approaches for the Inference of Relationships and Molecular Rates with Phylogenomic Data Sets.

Gene tree conflict is common and finding methods to analyze and alleviate the negative effects that conflict has on species tree analysis is a crucial part of phylogenomics. This study aims to expand the discussion of inferring species trees and molecular branch lengths when conflict is present. Conflict is typically examined in two ways: inferring its prevalence and inferring the influence of the individual genes (how strongly one gene supports any given topology compared to an alternative topology). Here, we examine a procedure for incorporating both conflict and the influence of genes in order to infer evolutionary relationships. All supported relationships in the gene trees are analyzed and the likelihood of the genes constrained to these relationships is summed to provide a likelihood for the relationship. Consensus tree assembly is conducted based on the sum of likelihoods for a given relationship and choosing relationships based on the most likely relationship assuming it does not conflict with a relationship that has a higher likelihood score. If it is not possible for all most likely relationships to be combined into a single bifurcating tree then multiple trees are produced and a consensus tree with a polytomy is created. This procedure allows for more influential genes to have a greater influence on an inferred relationship, does not assume conflict has arisen from any one source and does not force the data set to produce a single bifurcating tree. Using this approach, on three empirical data sets, we examine and discuss the relationship between influence and prevalence of gene tree conflict. We find that in one of the data sets, assembling a bifurcating consensus tree solely composed of the most likely relationships is impossible. To account for conflict in molecular rate analysis we also introduce a concordance-based approach to the summary and estimation of branch lengths suitable for downstream comparative analyses. We demonstrate through simulation that even under high levels of stochastic conflict, the mean and median of the concordant rates recapitulate the true molecular rate better than using a supermatrix approach. Using a large phylogenomic data set, we examine rate heterogeneity across concordant genes with a focus on the branch subtending crown angiosperms. Notably, we find highly variable rates of evolution along the branch subtending crown angiosperms. The approaches outlined here have several limitations, but they also represent some alternative methods for harnessing the complexity of phylogenomic data sets and enrich our inferences of both species relationships and evolutionary processes.[Branch length estimation; consensus tree; gene tree conflict; gene tree filtering; phylogenetics; phylogenomics.].