Abstract: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.

Resolving the Mortierellaceae phylogeny through synthesis of multi-gene phylogenetics and phylogenomics

Back together: Over 1000 single-copy nuclear loci and reproductive features support the holoendoparasitic Apodanthaceae and Rafflesiaceae as sister lineages in the order Malpighiales

Phylogenomics of Mulberries (morus, Moraceae) Inferred from Plastomes and Single Copy Nuclear Genes

A multigene phylogeny toward a new phylogenetic classification of Leotiomycetes

Incomplete Lineage Sorting and Introgression among Genera and Species of Liliaceae Tribe Tulipeae: Insights from Phylogenomics

Multiple historical processes obscure phylogenetic relationships in a taxonomically difficult group (Lobariaceae, Ascomycota)

New insights into the phylogenetic relationships within the Lauraceae from mitogenomes

Optimal data partitioning, multispecies coalescent and Bayesian concordance analyses resolve early divergences of the grape family (Vitaceae)

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

Modeling compositional heterogeneity resolves deep phylogeny of flowering plants

Large scale genome skimming from herbarium material for accurate plant identification and phylogenomics

The polyphyletic nature of Pleosporales: an example from Massariosphaeria based on rDNA and RBP2 gene phylogenies.

Organellar phylogenomics of Ophioglossaceae fern genera

Using target enrichment sequencing to study the higher-level phylogeny of the largest lichen-forming fungi family: Parmeliaceae (Ascomycota)

Phylogenetic reconstruction in the order Nymphaeales: ITS2 secondary structure analysis and in silico testing of maturase k (matK) as a potential marker for DNA bar coding

Genome Report: Pseudomolecule-scale genome assemblies of Drepanocaryum sewerzowii and Marmoritis complanata

Resolving conflict in eutherian mammal phylogeny using phylogenomics and the multispecies coalescent model.

Truly the best of both worlds: Merging lineage‐specific and universal probe kits to maximize phylogenomic inference

Untargeted Metabolomics for Integrative Taxonomy: Metabolomics, DNA Marker-Based Sequencing, and Phenotype Bioimaging

Genome-scale approaches to resolving incongruence in molecular phylogenies

A multi-gene phylogeny of Clavicipitaceae (Ascomycota, Fungi): identification of localized incongruence using a combinational bootstrap approach