Phylogenetic Inference Research Articles

Sympatric diversity pattern driven by the secondary contact of two deeply divergent lineages of the soybean pod borer Leguminivora glycinivorella.

The soybean pod borer, Leguminivora glycinivorella (Matsumura), is an important tortricid pest species widely distributed in most parts of China and its adjacent regions. Here, we analyzed the genetic diversity and population differentiation of L. glycinivorella using diverse genetic information including the standard cox1 barcode sequences, mitochondrial genomes (mitogenomes), and single-nucleotide polymorphisms (SNPs) from genotyping-by-sequencing. Based on a comprehensive sampling (including adults or larvae of L. glycinivorella newly collected at 22 of the total 30 localities examined) that covers most of the known distribution range of this pest, analyses of 543 cox1 barcode sequences and 60 mitogenomes revealed that the traditionally recognized and widely distributed L. glycinivorella contains two sympatric and widely distributed genetic lineages (A and B) that were estimated to have diverged ∼1.14 million years ago during the middle Pleistocene. Moreover, low but statistically significant correlations were recognized between genetic differentiation and geographic or environmental distances, indicating the existence of local adaptation to some extent. Based on SNPs, phylogenetic inference, principal component analysis, fixation index, and admixture analysis all confirm the two divergent sympatric lineages. Compared with the stable demographic history of Lineage B, the expansion of Lineage A had possibly made the secondary contact of the two lineages probable, and this process may be driven by the climate fluctuation during the late Pleistocene as revealed by ecological niche modeling.

Phylogenetic tree instability after taxon addition: empirical frequency, predictability, and consequences for online inference.

Online phylogenetic inference methods add sequentially arriving sequences to an inferred phylogeny without the need to recompute the entire tree from scratch. Some online method implementations exist already, but there remains concern that additional sequences may change the topological relationship among the original set of taxa. We call such a change in tree topology a lack of stability for the inferred tree. In this paper, we analyze the stability of single taxon addition in a Maximum Likelihood framework across 1, 000 empirical datasets. We find that instability occurs in almost 90% of our examples, although observed topological differences do not always reach significance under the AU-test. Changes in tree topology after addition of a taxon rarely occur close to its attachment location, and are more frequently observed in more distant tree locations carrying low bootstrap support. To investigate whether instability is predictable, we hypothesize sources of instability and design summary statistics addressing these hypotheses. Using these summary statistics as input features for machine learning under random forests, we are able to predict instability and can identify the most influential features. In summary, it does not appear that a strict insertion-only online inference method will deliver globally optimal trees, although relaxing insertion strictness by allowing for a small number of final tree rearrangements or accepting slightly suboptimal solutions appears feasible.

Novel genetic variants of banana streak MY virus and banana streak IM virus naturally infecting banana in Northeast India.

Divergent banana streak viruses (BSV) were characterized from banana plants exhibiting diverse symptoms in the Northeast region (NER) of India. Using rolling circle amplification (RCA), the complete genome sequences of seven episomal banana streak MY virus (BSMYV) isolates, including two novel variants, and two new banana streak IM virus (BSIMV) isolates were characterized. The novel BSMYV genetic variants were associated with conspicuous necrosis on newly emerged leaves, peduncle distortion, pseudostem internal necrosis, in addition to common streak symptoms. For complete genome nucleotide sequences, BSMYV-IN4 and IN5 shared 77-79% identity with other BSMYVs, while BSMYV-IN7 and IN8 exhibited identities of 77-97%. This study reports for the first time, the complete genomes of two banana streak IM virus (BSIMV-IN1 and -IN2) infecting triploid banana hybrids exhibiting leaf distortion, stunted rosette-like growth, and necrosis, sharing 87% sequence identity with reference BSIMV genome (GenBank accession no. HQ593112). Phylogenetic inference based on complete genomes revealed the distinct and congruent placement of BSMYV-IN4 and IN5 within the BSMYV cluster. Pairwise sequence comparisons of the conserved RT/RNase H nucleotide (nt) sequences revealed that the BSMYV-IN7 and IN4 isolates showed 85% and 97% identity to BSMYV (AY805074), respectively, which shared highest nt identity with BSMYV-IN6, IN9, and IN10, at 100%. The RT/RNase H nt sequences of BSIMV-IN1 and IN2 had 98% identity with the BSIMV (HQ593112), but were characterized as novel variants of BSIMV based on complete genomes. An analysis of relative synonymous codon usage (RSCU) pattern in the ORFIII polyprotein of BSMYV and BSIMV isolates revealed AGA and AGG (arginine) as the most frequently overrepresented codons (>1.5), evolutionary conserved in the genome of both species. A total of 14 recombination events were detected among the 36 BSV genomes, with recombination breakpoints mainly located in the ORFI, III, and IGR genomic regions. A novel phylogenetic cluster, comprised of BSMYV-IN4 and IN5 within the clade I was probably derived from heterologous recombination between parents resembling banana streak VN virus (BSVNV; AY750155) and banana streak GF virus (BSGFV; KJ013507) isolates. The present study conclusively reports the infection of genetically and symptomatically distinct variants of BSMYV and BSIMV infecting banana hybrids in NER India.

Taxonomic and phylogenetic evidence reveal two new Volvariella species (Agaricales, Volvariellaceae) from Denmark

Two new species of Volvariella from Denmark, V. similis and V. sulcata, are described with phylogenetic inference based on multiple loci (ITS1-5.8S-ITS2 and 28S). Volvariella similis is characterized by its brownish pileus, grey stipe without caulocystidia, clavate to broadly subutriform pleurocystidia, and cylindrical, tapering pileipellis elements. Volvariella sulcata is distinguished by its smaller basidiomata with a short stipe, sulcate-striate pileus margin, 35–45 lamellae, lageniform pleurocystidia with a long neck and pedicel, and a preference for man-made habitats. Comprehensive descriptions, illustrations, line drawings, and comparisons with morphologically similar and phylogenetically related species are provided.

Evidence for a Miocene pulse of diversification of the tropical American clade of the Brazil nut family (Lecythidaceae)

ABSTRACT The tropical American species of the Brazil nut family, Lecythidaceae, are abundant in the tree flora, with some economically important species. Yet, the phylogenetic relationships and pace of diversification within this family are understudied. Here, we used shotgun sequencing data for 86 of the 228 currently accepted species in the exclusively tropical American subfamily Lecythidoideae in a phylogenomic context. For each sampled species, we built the full plastid DNA sequence, and also extracted nuclear DNA for 571 regions using the new bioinformatics pipeline REFMAKER, which we used to produce a time-calibrated phylogenetic hypothesis. Our analysis shows that phylogenetic inference from plastid and nuclear DNA alignments resulted in different topologies. The nuclear DNA topology strongly suggests that genus Lecythis should be split into at least four genera. Samples of the genus Eschweilera formed a monophyletic group, with the exception of one sampled species (Eschweilera amazoniciformis S.A.Mori). The Bertholletia clade, which contains the majority of the Lecythidoideae species, and all Lecythis and Eschweilera species, started diversifying around 27.5 Ma, with an accelerated rate of diversification starting in the middle Miocene (c. 12 Ma). The clade sister to Bertholletia (including Corythophora, Eschweilera, and Corrugata, Chartacea and Poiteaui clades) includes at least 124 species and it has diversified less than 10 Ma. This time frame of diversification coincides with major changes in tropical American landscapes and climate associated with the Andean uplift.

Do morphometric data improve phylogenetic reconstruction? A systematic review and assessment

BackgroundIsolating phylogenetic signal from morphological data is crucial for accurately merging fossils into the tree of life and for calibrating molecular dating. However, subjective character definition is a major limitation which can introduce biases that mislead phylogenetic inferences and divergence time estimation. The use of quantitative data, e.g., geometric morphometric (GMM; shape) data can allow for more objective integration of morphological data into phylogenetic inference. This systematic review describes the current state of the field in using continuous morphometric data (e.g., GMM data) for phylogenetic reconstruction and assesses the efficacy of these data compared to discrete characters using the PRISMA-EcoEvo v1.0. reporting guideline, and offers some pathways for approaching this task with GMM data. A comprehensive search string yielded 11,123 phylogenetic studies published in English up to Oct 2023 in the Web of Science database. Title and abstract screening removed 10,975 articles, and full-text screening was performed for 132 articles. Of these, a total of twelve articles met final inclusion criteria and were used for downstream analyses.ResultsPhylogenetic performance was compared between approaches that employed continuous morphometric and discrete morphological data. Overall, the reconstructed phylogenies did not show increased resolution or accuracy (i.e., benchmarked against molecular phylogenies) as continuous data alone or combined with discrete morphological datasets.ConclusionsAn exhaustive search of the literature for existing empirical continuous data resulted in a total of twelve articles for final inclusion following title/abstract, and full-text screening. Our study was performed under a rigorous framework for systematic reviews, which showed that the lack of available comparisons between discrete and continuous data hinders our understanding of the performance of continuous data. Our study demonstrates the problem surrounding the efficacy of continuous data as remaining relatively intractable despite an exhaustive search, due in part to the difficulty in obtaining relevant comparisons from the literature. Thus, we implore researchers to address this issue with studies that collect discrete and continuous data sets with directly comparable properties (i.e., describing shape, or size).

On the importance of assessing topological convergence in Bayesian phylogenetic inference

Abstract Modern phylogenetics research is often performed within a Bayesian framework, using sampling algorithms such as Markov chain Monte Carlo (MCMC) to approximate the posterior distribution. These algorithms require careful evaluation of the quality of the generated samples. Within the field of phylogenetics, one frequently adopted diagnostic approach is to evaluate the effective sample size (ESS) and to investigate trace graphs of the sampled parameters. A major limitation of these approaches is that they are developed for continuous parameters and therefore incompatible with a crucial parameter in these inferences: the tree topology. Several recent advancements have aimed at extending these diagnostics to topological space. In this reflection paper, we present two case studies – one on Ebola virus and one on HIV – illustrating how these topological diagnostics can contain information not found in standard diagnostics, and how decisions regarding which of these diagnostics to compute can impact inferences regarding MCMC convergence and mixing. Our results show the importance of running multiple replicate analyses and of carefully assessing topological convergence using the output of these replicate analyses. To this end, we illustrate different ways of assessing and visualizing the topological convergence of these replicates. Given the major importance of detecting convergence and mixing issues in Bayesian phylogenetic analyses, the lack of a unified approach to this problem warrants further action, especially now that additional tools are becoming available to researchers.

Power and sample size calculations for testing the ratio of reproductive values in phylogenetic samples.

The quality of the inferences we make from pathogen sequence data is determined by the number and composition of pathogen sequences that make up the sample used to drive that inference. However, there remains limited guidance on how to best structure and power studies when the end goal is phylogenetic inference. One question that we can attempt to answer with molecular data is whether some people are more likely to transmit a pathogen than others. Here we present an estimator to quantify differential transmission, as measured by the ratio of reproductive numbers between people with different characteristics, using transmission pairs linked by molecular data, along with a sample size calculation for this estimator. We also provide extensions to our method to correct for imperfect identification of transmission linked pairs, overdispersion in the transmission process, and group imbalance. We validate this method via simulation and provide tools to implement it in an R package, phylosamp.

A step in the deep evolution of Alvinellidae (Annelida: Polychaeta): A phylogenomic comparative approach based on transcriptomes

The Alvinellidae are a family of worms that are endemic to deep-sea hydrothermal vents in the Pacific and Indian Oceans. These annelid worms, a sister group to the Ampharetidae, occupy a wide range of thermal habitats. The family includes the most thermotolerant marine animals described to date such as the Pompeii worm Alvinella pompejana, and other species living at much lower temperatures such as Paralvinella grasslei or Paralvinella pandorae. The phylogeny of this family has not been studied extensively. It is, however, a complex case where molecular phylogenies have given conflicting results, especially concerning the monophyletic or polyphyletic nature of the genus Paralvinella.We carried out a comprehensive study of the phylogeny of this family using the best molecular data currently available from RNAseq datasets, leading to several hundred orthologous transcripts for 11 of the 14 species currently described or in description. The results obtained by the two most popular phylogenetic inference models (using either gene concatenation with maximum likelihood, or a coalescent-based model from gene trees) were compared using a series of ampharetid and terebellid outgroups.Our study shows that the global phylogenetic signal favors the hypothesis of paraphyly for the Paralvinella genus, with P. pandorae being sister species of the other Alvinellidae.However, a high number of gene trees also supports the hypothesis of alternative trees in which the monophyly of the Paralvinella genus, as initially proposed by Desbruyères and Laubier, is valid with the species P. pandorae and Paralvinella unidentata being sister species. According to molecular dating, the radiation of the Alvinellidae was rapid and took place in a short period of time between 80 and 110 million years ago. This is reflected at the genomic scale by gene trees equally separated between different phylogenetic hypothesis, showing high rates of incomplete lineage sorting between the first lineages of the Alvinellidae and probable gene transfers. Although different genomic regions seem to have different phylogenetic stories in the early step of the alvinellid radiation, our study supports the view that the two P. pandorae species can be grouped into a separate genus (possibly Nautalvinella) and that the Miralvinella subgenus, defined by Desbruyères and Laubier, is not valid anymore.

The Impact of Model Misspecification on Phylogenetic Network Inference

The development of statistical methods to infer species phylogenies with reticulations (species networks) has led to many discoveries of gene flow between distinct species. These methods typically assume only incomplete lineage sorting and introgression. Given that phylogenetic networks can be arbitrarily complex, these methods might compensate for model misspecification by increasing the number of dimensions beyond the true value. Herein, we explore the effect of potential model misspecification, including the negligence of gene tree estimation error (GTEE) and assumption of a single substitution rate for all genomic loci, on the accuracy of phylogenetic network inference using both simulated and biological data. In particular, we assess the accuracy of estimated phylogenetic networks as well as test statistics for determining whether a network is the correct evolutionary history, as opposed to the simpler model that is a tree. We found that while GTEE negatively impacts the performance of test statistics to determine the “treeness” of the evolutionary history of a data set, running those tests on triplets of taxa and correcting for multiple-testing significantly ameliorates the problem. We also found that accounting for substitution rate heterogeneity improves the reliability of full Bayesian inference methods of phylogenetic networks, whereas summary statistic methods are robust to GTEE and rate heterogeneity, though currently require manual inspection to determine the network complexity.

Genomes of Alphanucleorhabdovirus Physostegiae Isolates from Two Different Cultivar Groups of Solanum melongena

Plant rhabdoviruses cause considerable economic losses and are a threat to the agriculture of Solanaceae plants. Two novel virus isolates belonging to the family Rhabdoviridae are identified by high-throughput sequencing (HTS) in Russian eggplant cultivars grown in the Volga river delta region for the first time. The phylogenetic inference of L protein (polymerase) shows that these virus isolates belong to Alphanucleorhabdovirus physostegia (Physostegia chlorotic mottle virus—PhCMoV), and their minus-sense RNA genomes have the typical gene order 3′-nucleocapsid (N)—X protein (X)—phosphoprotein (P)—Y protein (Y)—matrix protein (M)—glycoprotein (G)—polymerase (L)-5′ observed in some plant-infecting alphanucleorhabdoviruses. One of the PhCMoV isolates from the eggplant cultivar Almaz is genetically very similar to the Russian PhCMoV isolate from tomato and grouped in a subclade together with four isolates from Belgium, Germany, the Netherlands, and France. However, another eggplant-infecting isolate from the Russian cultivar Boggart is the most divergent compared with the other 45 virus genomes of European PhCMoV isolates. Thus, our comparative analysis reveals that two virus isolates from Russia may either share a close evolutionary relationship with European isolates or significantly diverge from all known virus isolates. The potential to use the protein sequence comparative analysis of accessory polypeptides, along with the early developed strategy of the nucleotide sequence comparison of the RNA genomes, is shown.

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

AbstractPremiseHybridization capture kits are now commonly used for reduced representation approaches in genomic sequencing, with both universal and clade‐specific kits available. Here, we present a probe kit targeting 799 low‐copy genes for the plant family Annonaceae.MethodsThis new version of the kit combines the original 469 genes from the previous Annonaceae kit with 334 genes from the universal Angiosperms353 kit. We also compare the results obtained using the original Angiosperms353 kit with our custom approach using a subset of specimens. Parsimony‐informative sites and the results of maximum likelihood phylogenetic inference were assessed for combined matrices using the genera Asimina and Deeringothamnus.ResultsThe Annonaceae799 genes derived from the Angiosperms353 kit have extremely high recovery rates. Off‐target reads were also detected. When evaluating size, the proportion of on‐ and off‐target regions, and the number of parsimony‐informative sites, the genes incorporated from the Angiosperms353 panel generally outperformed the genes from the original Annonaceae probe kit.DiscussionWe demonstrated that the new sequences from the Angiosperms353 probe set are variable and relevant for future studies on species‐level phylogenomics and within‐species studies in the Annonaceae. The integration of kits also establishes a connection between projects and makes new genes available for phylogenetic and population studies.

Phylogenomics resolves a 100-year-old debate regarding the evolutionary history of caddisflies (Insecta: Trichoptera)

Trichoptera (caddisfly) phylogeny provides an interesting example of aquatic insect evolution, with rich ecological diversification, especially for underwater architecture. Trichoptera provide numerous critical ecosystem services and are also one of the most important groups of aquatic insects for assessing water quality. The phylogenetic relationships of Trichoptera have been debated for nearly a century. In particular, the phylogenetic position of the “cocoon-makers” within Trichoptera has long been contested. Here, we designed a universal single-copy orthologue and sets of ultraconserved element markers specific for Trichoptera for the first time. Simultaneously, we reconstructed the phylogenetic relationship of Trichoptera based on genome data from 111 species, representing 29 families and 71 genera. Our phylogenetic inference clarifies the probable phylogenetic relationships of “cocoon-makers” within Integripalpia. Hydroptilidae is considered as the basal lineage within Integripalpia, and the families Glossosomatidae, Hydrobiosidae, and Rhyacophilidae formed a monophyletic clade, sister to the integripalpian subterorder Phryganides. The resulting divergence time and ancestral state reconstruction suggest that the most recent common ancestor of Trichoptera appeared in the early Permian and that diversification was strongly correlated with habitat change.

An automated convergence diagnostic for phylogenetic MCMC analyses.

Assessing convergence of Markov chain Monte Carlo (MCMC) based analyses is crucial but challenging, especially so in high dimensional and complex spaces such as the space of phylogenetic trees (treespace). In practice, it is assumed that the target distribution is the unique stationary distribution of the MCMC and convergence is achieved when samples appear to be stationary. Here we leverage recent advances in computational geometry of the treespace and introduce a method that combines classical statistical techniques and algorithms with geometric properties of the treespace to automatically evaluate and assess practical convergence of phylogenetic MCMC analyses. Our method monitors convergence across multiple MCMC chains and achieves high accuracy in detecting both practical convergence and convergence issues within treespace. Furthermore, our approach is developed to allow for real-time evaluation during the MCMC algorithm run, eliminating any of the chain post-processing steps that are currently required. Our tool therefore improves reliability and efficiency of MCMC based phylogenetic inference methods and makes analyses easier to reproduce and compare. We demonstrate the efficacy of our diagnostic via a well-calibrated simulation study and provide examples of its performance on real data sets. Although our method performs well in practice, a significant part of the underlying treespace probability theory is still missing, which creates an excellent opportunity for future mathematical research in this area. The open source package for the phylogenetic inference framework BEAST2, called ASM, that implements these methods, making them accessible through a user-friendly GUI, is available from https://github.com/rbouckaert/asm/. The open source Python package, called tetres, that provides an interface for these methods enabling their applications beyond BEAST2 can be accessed at https://github.com/bioDS/tetres/.

Broadening the horizon of Indian mycobiota: two new species and one new generic record in Boletoid mushrooms (Boletaceae)

A combination of morphological studies and multigene molecular phylogenetic analyses unveiled two significant novel species of boletes namely Anthracoporus indohimalayanus sp. nov. and Boletus dhakuricus sp. nov. from the states of Sikkim and Uttarakhand, respectively. The first novel species, i.e. A. indohimalayanus is also a new generic record of the genus Anthracoporus from India, whereas the second novel species i.e. B. dhakuricus is a member of wild edible porcini mushroom. Both the taxa are described with their morphological details and molecular phylogenetic inferences.

Identifying Impacts of Contact Tracing on Epidemiological Inference from Phylogenetic Data.

Robust sampling methods are foundational to inferences using phylogenies. Yet the impact of using contact tracing, a type of non-uniform sampling used in public health applications such as infectious disease outbreak investigations, has not been investigated in the molecular epidemiology field. To understand how contact tracing influences a recovered phylogeny, we developed a new simulation tool called SEEPS (Sequence Evolution and Epidemiological Process Simulator) that allows for the simulation of contact tracing and the resulting transmission tree, pathogen phylogeny, and corresponding virus genetic sequences. Importantly, SEEPS takes within-host evolution into account when generating pathogen phylogenies and sequences from transmission histories. Using SEEPS, we demonstrate that contact tracing can significantly impact the structure of the resulting tree, as described by popular tree statistics. Contact tracing generates phylogenies that are less balanced than the underlying transmission process, less representative of the larger epidemiological process, and affects the internal/external branch length ratios that characterize specific epidemiological scenarios. We also examined real data from a 2007-2008 Swedish HIV-1 outbreak and the broader 1998-2010 European HIV-1 epidemic to highlight the differences in contact tracing and expected phylogenies. Aided by SEEPS, we show that the data collection of the Swedish outbreak was strongly influenced by contact tracing even after downsampling, while the broader European Union epidemic showed little evidence of universal contact tracing, agreeing with the known epidemiological information about sampling and spread. Overall, our results highlight the importance of including possible non-uniform sampling schemes when examining phylogenetic trees. For that, SEEPS serves as a useful tool to evaluate such impacts, thereby facilitating better phylogenetic inferences of the characteristics of a disease outbreak. SEEPS is available at github.com/MolEvolEpid/SEEPS.

Comparative plastome analyses and evolutionary relationships of 25 East Asian species within the medicinal plant genus Scrophularia (Scrophulariaceae)

BackgroudScrophularia L., a genus of the Scrophulariaceae, is a group of important medicinal plants used for eliminating heat and detoxifying. East Asia has an abundance of potentially medicinal Scrophularia species, and it serves as a secondary diversity center of the genus. However, the genomic resources available for germplasm identification and pharmaceutical exploration of East Asian Scrophularia are insufficient, hindering its commercial and industrial development. Additionally, the interspecific relationships of most East Asian Scrophularia species remain unclear.MethodsIn this study, we sequenced the leaves of 25 East Asian species of the genus Scrophularia, assembled and annotated the complete chloroplast genomes, and subsequently performed comparative and phylogenetic analyses on these genomes.Results and discussionThe conserved plastome length of these 25 species ranged from 151,582 bp to 153,239 bp, containing a total of 132 coding genes, including 18 duplicated genes and 114 unique genes. Through genome alignment of these 25 species, 38-53 repeated sequences and 7 shared SSRs were identified, along with regions with high nucleotide polymorphism (Pi), which could potentially serve as molecular markers for species identification. The genome structure, gene content, and arrangement showed conservation, while variations were observed in the IR boundary regions and IGS. Phylogenetic inferences based on whole plastomes or on coding sequences (CDS) only yielded congruent results. We categorized the 25 East Asian Scrophularia species into six distinct clades and further explored their interspecies relationships using morphological characteristics, such as flower color, the relative position of stamens and corolla, and plant height. This could lay a genetic basis for future resource development of Scrophularia in East Asia.

Epidemiology of Staphylococcus aureus food isolates: Comparison of conventional methods with whole genome sequencing typing methods

A variety of methods exists for typing bacteria. However, guidelines for the application and interpretation of typing tools in epidemiologic investigations of Staphylococcus aureus are lacking. This study aimed to identify appropriate typing methods for S. aureus population studies and outbreak investigation. We compared pulsed-field gel electrophoresis (PFGE), seven loci multi-locus sequence typing (MLST), core genome MLST (cgMLST), core single nucleotide polymorphism (cSNP), and enterotoxin (se/SE) profiles on 351 S. aureus isolates. The discriminatory power, concordance, and congruence of typing results were assessed. cgMLST, cSNP, and PFGE yielded the highest discrimination value, followed by se/SE typing and MLST. The best concordance of results was found between cgMLST and cSNP, while the best congruence was observed for cgMLST and cSNP with all methods, followed by PFGE with MLST. The strengths and weaknesses of each method are highlighted. For population structure, cgMLST and cSNP performed better than PFGE and MLST in terms of resolution of clusters and in phylogenetic inference. Enterotoxin profiles matched with MLST groups, suggesting the use of se/SE typing to predict MLST results. For the retrospective analysis of 31 outbreaks, all methods performed almost equally to discriminate epidemiologically related strains and can be used to unambiguously distinguish outbreak strains.

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.

An integrative approach to the delineation of locally endemic species: A case study of Vitis baihuashanensis (Vitaceae) and its implication for conservation management

Species definition and classification have long been the subject of controversy, and name changes can often be a source of frustration for both conservationists and the public. Vitis baihuashanensis is a critically endangered grapevine endemic to Beijing with only two remaining wild individuals. Although formerly treated as a synonym of V. amurensis, it was later listed in the List of National Key Protected Wild Plants (LNKPWP) in 2021. Here, we evaluated the taxonomic identity of V. baihuashanensis based on phylogenomic inference and population genetic analysis, and tested the species boundaries under various major species concepts (biological, morphological, ecological, phylogenetic, and integrative species concepts). Phylogenetic inference revealed that the amurensis-baihuashanensis clade was monophyletic, samples of V. baihuashanensis nested within the V. amurensis node. Population genetic analysis demonstrated that V. baihuashanensis samples clustered with V. amurensis. Our findings from field observations, morphological comparisons, pollination experiments, and genetic analyses indicate that V. baihuashanensis and V. amurensis are the same species under the integrative species concept, and support the synonymy of V. baihuashanensis under V. amurensis. Given that V. baihuashanensis was listed in the LNKPWP, we suggest that it be removed in the next update of the list because of limited financial resources. Nevertheless, the two wild individuals deserve conservation attention because of their distinct leaf characteristics and potential of cold tolerance, which could have important applications in the horticulture and grape industries. Our study provides multiple lines of evidence for clarifying the identity of a taxonomically questionable species, as well as valuable suggestion with conservation and management implications.