Reticulate evolution: Detection and utility in the phylogenomics era.

Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.

Mitochondrial phylogeography of grassland caterpillars (Lepidoptera: Lymantriinae: Gynaephora) endemic to the Qinghai-Tibetan plateau.

Grassland caterpillars (Lepidoptera: Lymantriinae: Gynaephora) are the most damaging pests to alpine meadows in the Qinghai-Tibetan Plateau (QTP). Here, we conducted extensive sampling from 39 geographic populations covering almost the entire distribution of the eight QTP Gynaephora (Hübner) species to investigate phylogeographic patterns and speciation based on two mitochondrial genes (COI and ND5). A total of 40 haplotypes were detected in the 39 populations, with >70% of all haplotypes not shared between populations. The monophyletic QTP Gynaephora migrated from non-QTP regions during the Pliocene, corresponding to the uplift of the QTP, suggesting a mode of transport into the QTP. Among the eight QTP Gynaephora species described by morphological characteristics, two species (G. alpherakii and G. menyuanensis) were recovered as monophyletic groups (Clades B and C), while the remaining six formed two monophyletic clades: Clade A (G. qinghaiensis, G. jiuzhiensis, and G. qumalaiensis) and Clade D (G. aureata, G. ruoergensis, and G. minora). These results suggested that the number of the QTP Gynaephora species may be overestimated and further studies based on both morphological and nuclear gene data are needed. Genetic differentiation and speciation of the QTP Gynaephora were likely driven by the QTP uplifts and associated climate fluctuations during the Pleistocene, indicated by divergence time estimation, suggesting that isolation and subsequent divergence was the dominant mode of speciation. The Sanjiangyuan region (i.e., Clade A, characterized by high genetic diversity) may have been a glacial refugium of the QTP Gynaephora, as supported by analyses of gene flow and biogeography. High levels of genetic diversity were found in QTP Gynaephora, without population expansion, which may explain the high-altitude adaptation and outbreaks of grassland caterpillars in alpine meadows of the QTP. This study provides the largest phylogeographic analysis of QTP Gynaephora and improves our understanding of the diversity and speciation of QTP insects.

A new discovery of cryptic species in barred sand burrower Limnichthys fasciatus (Pisces: Creediidae) from the western Pacific, with evolutionary perspectives of anti-equatorial species.

Although the barred sand burrower, Limnichthys fasciatus, is widely distributed throughout the western Pacific, including Japan, Taiwan, and Australia, its morphology and genetics are poorly known. We discovered four cryptic species of Limnichthys from the western Pacific based on mtDNA cytochrome c oxidase subunit I (COI) and 16S ribosomal RNA (rRNA) sequences. Genetic distances showed remarkably large differences (12.7%-24.3% in COI and 7.9%-19.6% in 16S rRNA) between true L. fasciatus (type locality: southeastern Australia) and the others. A relaxed clock model with optimized selected substitution models showed that their deep divergence began in the middle Miocene epoch and subsequently diverged into the current cryptic species in the Plio-Pleistocene. A eurythermal common ancestor may have evolved independently in each region due to geographical events and paleoclimatic fluctuations, which made it possible for L. fasciatus complex to be an anti-equatorial species. Despite their deep genetic divergence, they showed marked phenotypic similarity, suggesting that they have experienced similar selective pressures related to their specific behavior.

Characterization of the complete mitochondrial genome of a newly discovered torrent catfish, Liobagrus geumgangensis, and their phylogenetic relationships.

BACKGROUND: A new Liobagrus fish was reported from the Korean Peninsula, but research on this taxon is lacking. Moreover, existing research on the mitogenome of the genus Liobagrus in Korea is very limited, and no studies have been conducted on structural characteristics of transfer RNA (tRNA) or gene order comparisons between taxa; instead, research has been restricted to basic phylogeny. OBJECTIVE: The complete mitochondrial genome of Liobagrus geumgangensis was analyzed for the first time. We then aimed to reconstruct the phylogenetic relationships of the genus Liobagrus and estimate the divergence time of speciation events. METHODS: We used a dissected fin clip from an adult of Liobagrus geumgangensis. Genomic DNA was extracted and analyzed with whole genome sequencing (WGS) and assembled by the NOVOPlasty method. The mitogenome sequence was annotated, and a genome map, tRNA structure, and phylogenetic tree were constructed using maximum likelihood analysis. In addition, divergence time was estimated. RESULTS: The mitochondrial genome was 16,522 bp in length and comprised 37 genes. The overall base composition was 30.5% A, 25.5% T, 28.4% C, and 15.7% G. Most tRNAs exhibited the typical clover leaf shape, except trnS1. Phylogenetic analysis revealed that Liobagrus geumgangensis clustered within a clade with four other Liobagrus species exclusive to the southern region of the Korean Peninsula. Its divergence was estimated to have occurred during the late Miocene. CONCLUSION: Characteristics of Liobagrus geumgangensis mitogenome were consistent with those of other torrent catfish species. Time scale estimation revealed distinct groupings, with some distributed across mainland Asia and others in the southern region of the Korean Peninsula. Notably, the Korean Peninsula group was identified as its own lineage, comprising entirely endemic species.

Phylogeography of Parnassius citrinarius Based on Mitochondrial DNA Reveals Large Differences in Genetic Structure between the Eastern and Western Japan.

The Japanese Archipelago hosts a rich butterfly fauna, and elucidating the genetic structures of multiple species is necessary to clarify their formation processes. This study aimed to reveal the genetic structure and distribution formation process of Parnassius citrinarius, which is widely distributed across the Japanese Archipelago from Hokkaido to Shikoku, through phylogeographic analysis based on the mitochondrial cytochrome c oxidase subunit I (COI) gene sequence. Thirty haplotypes were revealed from 311 individuals from 47 sites, indicating significant differences in the genetic structures between the eastern and western parts of the Japanese Archipelago. In Eastern Japan, multiple genetic clusters were found, with some sites harboring two clusters. The divergence times among populations in Eastern Japan were relatively recent, and no genetic differentiation was observed between regions, including between Hokkaido and Honshu, which are separated by a narrow strait. In contrast, in Western Japan, including Shikoku, unique genetic clusters were observed in each region. The phylogenetic relationships among populations were regionally clustered, and the divergence times were relatively ancient. The distribution and genetic structure of P. citrinarius in the Japanese Archipelago have been significantly influenced by temperature fluctuations and the presence of geographical barriers during the Pleistocene glacial-interglacial cycles, including the potential formation of refugia in Western Japan.

A molecular phylogeny of the Petaluridae (Odonata: Anisoptera): A 160-Million-Year-Old story of drift and extinction.

Petaluridae (Odonata: Anisoptera) is a relict dragonfly family, having diverged from its sister family in the Jurassic, of eleven species that are notable among odonates (dragonflies and damselflies) for their exclusive use of fen and bog habitats, their burrowing behavior as nymphs, large body size as adults, and extended lifespans. To date, several nodes within this family remain unresolved, limiting the study of the evolution of this peculiar family. Using an anchored hybrid enrichment dataset of over 900 loci we reconstructed the species tree of Petaluridae. To estimate the temporal origin of the genera within this family, we used a set of well-vetted fossils and a relaxed molecular clock model in a divergence time estimation analysis. We estimate that Petaluridae originated in the early Cretaceous and confirm the existence of monophyletic Gondwanan and Laurasian clades within the family. Our relaxed molecular clock analysis estimated that these clades diverged from their MRCA approximately 160 mya. Extant lineages within this family were identified to have persisted from 6 (Uropetala) to 120 million years (Phenes). Our biogeographical analyses focusing on a set of key regions suggest that divergence within Petaluridae is largely correlated with continental drift, the exposure of land bridges, and the development of mountain ranges. Our results support the hypothesis that species within Petaluridae have persisted for tens of millions of years, with little fossil evidence to suggest widespread extinction in the family, despite optimal conditions for the fossilization of nymphs. Petaluridae appear to be a rare example of habitat specialists that have persisted for tens of millions of years.

Ad fontes: divergence-time estimation and the age of angiosperms.

Accurate divergence times are essential for interpreting and understanding the context in which lineages have evolved. Over the past several decades, debates have surrounded the discrepancies between the inferred molecular ages of crown angiosperms, often estimated from the Late Jurassic into the Permian, and the fossil record, placing angiosperms in the Early Cretaceous. That crown angiosperms could have emerged as early as the Permian or even the Triassic would have major implications for the paleoecological context of the origin of one of the most consequential clades in the tree of life. Here, we argue, and demonstrate through simulations, that the older ages inferred from molecular data and relaxed-clock models are misled by lineage-specific rate heterogeneity resulting from life history changes that occurred several times throughout the evolution of vascular plants. To overcome persistent discrepancies in age estimates, more biologically informed and realistic models should be developed, and our results should be considered in the context of their biological implications before we accept inferences that are a major departure from our strongest evidence.

An updated phylogeny and adaptive evolution within Amaranthaceae s.l. inferred from multiple phylogenomic datasets.

Amaranthaceae s.l. is a widely distributed family consisting of over 170 genera and 2000 species. Previous molecular phylogenetic studies have shown that Amaranthaceae s.s. and traditional Chenopodiaceae form a monophyletic group (Amaranthaceae s.l.), however, the relationships within this evolutionary branch have yet to be fully resolved. In this study, we assembled the complete plastomes and full-length ITS of 21 Amaranthaceae s.l. individuals and compared them with 38 species of Amaranthaceae s.l. Through plastome structure and sequence alignment analysis, we identified a reverse complementary region approximately 5200 bp long in the genera Atriplex and Chenopodium. Adaptive evolution analysis revealed significant positive selection in eight genes, which likely played a driving role in the evolution of Amaranthaceae s.l., as demonstrated by partitioned evolutionary analysis. Furthermore, we found that about two-thirds of the examined species lack the ycf15 gene, potentially associated with natural selection pressures from their adapted habitats. The phylogenetic tree indicated that some genera (Chenopodium, Halogeton, and Subtr. Salsolinae) are paraphyletic lineages. Our results strongly support the clustering of Amaranthaceae s.l. with monophyletic traditional Chenopodiaceae (Clades I and II) and Amaranthaceae s.s. After a comprehensive analysis, we determined that cytonuclear conflict, gene selection by adapted habitats, and incomplete lineage sorting (ILS) events were the primary reasons for the inconsistent phylogeny of Amaranthaceae s.l. During the last glacial period, certain species within Amaranthaceae s.l. underwent adaptations to different environments and began to differentiate rapidly. Since then, these species may have experienced morphological and genetic changes distinct from those of other genera due to intense selection pressure.

Characteristics of the chloroplast genome and genetic divergence of Tamarix hispida Willd. 1816 (Tamaricaceae).

Tamarix hispida Willd. 1816, a crucial native plant species in the arid desert region of northwestern China, plays a significant role in maintaining ecological stability. It is instrumental in addressing soil salinity-alkalinity and heavy metal pollution. This research aims to analyze the phylogenetic divergence pattern and evolutionary history of T. hispida by comparing chloroplast genome structures across different populations. Despite the minimal differences in chloroplast genome structure due to conserved genes and junction regions, sequencing was conducted using the Illumina NovaSeq platform to verify the historical evolutionary processes between different populations, followed by assembly and annotation. The results revealed that the T. hispida chloroplast genome is approximately 156,164-156,186 bp in length, with a quadripartite structure and 131 annotated genes. Phylogenetic analysis indicated two lineages within T. hispida, with a divergence time of 3.15 Ma. These findings emphasize the low genetic diversity in T. hispida and offer valuable insights into its evolutionary past. To effectively protect and manage this species, increased scientific research and monitoring of its genetic diversity are necessary. This study underscores the importance of comprehending the genetic mechanisms behind species divergence to develop informed conservation strategies.

Phylogenetic insights into the Salicaceae: The evolution of willows and beyond.

The Salicaceae includes approximately 54 genera and over 1,400 species with a cosmopolitan distribution. Members of the family are well-known for their diverse secondary plant metabolites, and they play crucial roles in tropical and temperate forest ecosystems. Phylogenetic reconstruction of the Salicaceae has been historically challenging due to the limitations of molecular markers and the extensive history of hybridization and polyploidy within the family. Our study employs whole-genome sequencing of 74 species to generate an extensive phylogeny of the Salicaceae. We generated two RAD-Seq enriched whole-genome sequence datasets and extracted two additional gene sets corresponding to the universal Angiosperms353 and Salicaceae-specific targeted-capture arrays. We reconstructed maximum likelihood-based molecular phylogenies using supermatrix and coalescent-based supertree approaches. Our fossil-calibrated phylogeny estimates that the Salicaceae originated around 128 million years ago and unravels the complex taxonomic relationships within the family. Our findings confirm the non-monophyly of the subgenus Salix s.l. and further support the merging of subgenera Chamaetia and Vetrix, both of which exhibit intricate patterns within and among different sections. Overall, our study not only enhances our understanding of the evolution of the Salicaceae, but also provides valuable insights into the complex relationships within the family.

Structure, Evolution, and Mitochondrial Genome Analysis of Mussel Species (Bivalvia, Mytilidae).

Based on the nucleotide sequences of the mitochondrial genome (mitogenome) of specimens taken from two mussel species (Arcuatula senhousia and Mytilus coruscus), an investigation was performed by means of the complex approaches of the genomics, molecular phylogenetics, and evolutionary genetics. The mitogenome structure of studied mussels, like in many other invertebrates, appears to be much more variable than in vertebrates and includes changing gene order, duplications, and deletions, which were most frequent for tRNA genes; the mussel species' mitogenomes also have variable sizes. The results demonstrate some of the very important properties of protein polypeptides, such as hydrophobicity and its determination by the purine and pyrimidine nucleotide ratio. This fact might indirectly indicate the necessity of purifying natural selection for the support of polypeptide functionality. However, in accordance with the widely accepted and logical concept of natural cutoff selection for organisms living in nature, which explains its action against deleterious nucleotide substitutions in the nonsynonymous codons (mutations) and its holding of the active (effective) macromolecules of the polypeptides in a population, we were unable to get unambiguous evidence in favor of this concept in the current paper. Here, the phylogeny and systematics of mussel species from one of the largest taxons of bivalve mollusks are studied, the family known as Mytilidae. The phylogeny for Mytilidae (order Mytilida), which currently has no consensus in terms of systematics, is reconstructed using a data matrix of 26-27 mitogenomes. Initially, a set of 100 sequences from GenBank were downloaded and checked for their gender: whether they were female (F) or male (M) in origin. Our analysis of the new data confirms the known drastic differences between the F/M mitogenome lines in mussels. Phylogenetic reconstructions of the F-lines were performed using the combined set of genetic markers, reconstructing only protein-coding genes (PCGs), only rRNA + tRNA genes, and all genes. Additionally, the analysis includes the usage of nucleotide sequences composed of other data matrices, such as 20-68 mitogenome sequences. The time of divergence from MRCA, estimated via BEAST2, for Mytilidae is close to 293 Mya, suggesting that they originate in the Silurian Period. From all these data, a consensus for the phylogeny of the subfamily of Mytilinae and its systematics is suggested. In particular, the long-debated argument on mussel systematics was resolved as to whether Mytilidae, and the subfamily of Mytilinae, are monophyletic. The topology signal, which was strongly resolved in this paper and in the literature, has refuted the theory regarding the monophyly of Mytilinae.

Mitochondrial Genome Characteristics Reveal Evolution of Acanthopsetta nadeshnyi (Jordan and Starks, 1904) and Phylogenetic Relationships.

In the present study, the mitochondrial genomic characteristics of Acanthopsetta nadeshnyi have been reported and have depicted the phylogenetic relationship among Pleuronectidae. Combined with a comparative analysis of 13 PCGs, the TN93 model was used to review the neutral evolution and habitat evolution catalysis of the mitogenome to verify the distancing and purification selectivity of the mitogenome in Pleuronectidae. At the same time, a species differentiation and classification model based on mitogenome analysis data was established. This study is expected to provide a new perspective on the phylogenetic relationship and taxonomic status of A. nadeshnyi and lay a foundation for further exploration of environmental and biological evolutionary mechanisms.

Morphological observation, molecular identification and evolutionary analysis of Hydatigera kamiyai found in Neodon fuscus from the Qinghai-Tibetan plateau.

Hydatigera kamiyai (H. kamiyai) is a new species within Hydatigera that has recently been resurrected. Voles and cats are hosts of H. kamiyai and have a certain impact on its health and economy. Moreover, the Qinghai-Tibetan plateau (QTP) is a research hotspot representing Earth's biodiversity, as its unique geographical environment and climatic conditions support the growth of a variety of mammals and provide favorable conditions for various parasites to complete their life history. The aim of this study was to reveal the phylogenetic relationships and divergence times of H. kamiyai strains isolated from Neodon fuscus on the QTP using morphological and molecular methods. In this study, we morphologically observed H. kamiyai and sequenced the whole mitochondrial genome. Then, we constructed phylogenetic trees with the maximum likelihood (ML) and Bayesian inference (BI) methods. The GTR alternative model was selected for divergence time analysis. These data demonstrated that the results were consistent with the general morphological characteristics of Hydatigera. The whole genome of H. kamiyai was 13,822 bp in size, and the A + T content (73%) was greater than the G + C content (27%). The Ka/Ks values were all <1, indicating that all 13 protein-coding genes (13 PCGs) underwent purifying selection during the process of evolution. The phylogenetic tree generated based on the 13 PCGs, cytochrom oxidase subunit I (COI), 18S rRNA and 28S rRNA revealed close phylogenetic relationships between H. kamiyai and Hydatigera, with high node support for the relationship. The divergence time based on 13 PCGs indicated that H. kamiyai diverged approximately 11.3 million years ago (Mya) in the Miocene. Interestingly, it diverged later than the period of rapid uplift in the QTP. We also speculated that H. kamiyai differentiation was caused by host differentiation due to the favorable living conditions brought about by the uplift of the QTP. As there have been relatively few investigations on the mitochondrial genome of H. kamiyai, our study could provide factual support for further studies of H. kamiyai on the QTP. We also emphasized the importance of further studies of its hosts, Neodon fuscus and cats, which will be important for further understanding the life cycle of H. kamiyai.

Characterization of the organellar genomes of Mazzaella laminarioides and Mazzaella membranacea (Gigartinaceae, Rhodophyta).

Mazzaella, a genus with no genomic resources available, has extensive distribution in the cold waters of the Pacific, where they represent ecologically and economically important species. In this study, we aimed to sequence, assemble, and annotate the complete mitochondrial and chloroplast genomes from two Mazzaella spp. and characterize the intraspecific variation among them. We report for the first time seven whole organellar genomes (mitochondria: OR915856, OR947465, OR947466, OR947467, OR947468, OR947469, OR947470; chloroplast: OR881974, OR909680, OR909681, OR909682, OR909683, OR909684, OR909685) obtained through high-throughput sequencing for six M. laminarioides sampled from three Chilean regions and one M. membranacea. Sequenced Mazzaella mitogenomes have identical gene number, gene order, and genome structure. The same results were observed for assembled plastomes. A total of 52 genes were identified in mitogenomes, and a total of 235 genes were identified in plastomes. Although the M. membranacea plastome included a full-length pbsA gene, in all M. laminarioides samples, the pbsA gene was split in three open reading frames (ORFs). Within M. laminarioides, we observed important plastome lineage-specific variations, such as the pseudogenization of the two hypothetical protein-coding genes, ycf23 and ycf45. Nonsense mutations in the ycf23 and ycf45 genes were only detected in the northern lineage. These results are consistent with phylogenetic reconstructions and divergence time estimation using concatenated coding sequences that not only support the monophyly of M. laminarioides but also underscore that the three M. laminarioides lineages are in an advanced stage of divergence. These new results open the question of the existence of still undisclosed species in M. laminarioides.

Plastid phylogenomics reveals evolutionary relationships in the mycoheterotrophic orchid genus Dipodium and provides insights into plastid gene degeneration.

The orchid genus Dipodium R.Br. (Epidendroideae) comprises leafy autotrophic and leafless mycoheterotrophic species, with the latter confined to sect. Dipodium. This study examined plastome degeneration in Dipodium in a phylogenomic and temporal context. Whole plastomes were reconstructed and annotated for 24 Dipodium samples representing 14 species and two putatively new species, encompassing over 80% of species diversity in sect. Dipodium. Phylogenomic analysis based on 68 plastid loci including a broad outgroup sampling across Orchidaceae found that sect. Leopardanthus is the sister lineage to sect. Dipodium. Dipodium ensifolium, the only leafy autotrophic species in sect. Dipodium, was found to be a sister to all leafless, mycoheterotrophic species, supporting a single evolutionary origin of mycoheterotrophy in the genus. Divergence-time estimations found that Dipodium arose ca. 33.3 Ma near the lower boundary of the Oligocene and that crown diversification commenced in the late Miocene, ca. 11.3 Ma. Mycoheterotrophy in the genus was estimated to have evolved in the late Miocene, ca. 7.3 Ma, in sect. Dipodium. The comparative assessment of plastome structure and gene degradation in Dipodium revealed that plastid ndh genes were pseudogenised or physically lost in all Dipodium species, including in leafy autotrophic species of both Dipodium sections. Levels of plastid ndh gene degradation were found to vary among species as well as within species, providing evidence of relaxed selection for retention of the NADH dehydrogenase complex within the genus. Dipodium exhibits an early stage of plastid genome degradation, as all species were found to have retained a full set of functional photosynthesis-related genes and housekeeping genes. This study provides important insights into plastid genome degradation along the transition from autotrophy to mycoheterotrophy in a phylogenomic and temporal context.

Phylogenetics and biogeography of the olive family (Oleaceae).

BACKGROUND AND AIMS: Progress in the systematic studies of the olive family (Oleaceae) during the last two decades provides the opportunity to update its backbone phylogeny and to investigate its historical biogeography. We additionally aimed to understand the factors underlying the disjunct distribution pattern between East Asia and both West Asia and Europe that is found more commonly in this family than in any other woody plant families. METHODS: Using a sampling of 298 species out of ca. 750, the largest in a phylogenetic study of Oleaceae thus far, and a set of 36 plastid and nuclear markers, we reconstructed and dated a new phylogenetic tree based on maximum likelihood and Bayesian methods and checked for any reticulation events. We also assessed the relative support of four competing hypotheses [Qinghai-Tibet Plateau uplift (QTP-only hypothesis), climatic fluctuations (Climate-only hypothesis), combined effects of QTP uplift and climate (QTP-Climate hypothesis), and no effects (Null hypothesis)] in explaining these disjunct distributions. KEY RESULTS: We recovered all tribes and subtribes within Oleaceae as monophyletic, but uncertainty in the position of tribe Forsythieae remains. Based on this dataset, no reticulation event was detected. Our biogeographic analyses support the QTP-Climate hypothesis as the likely main explanation for the East-West Eurasian disjunctions in Oleaceae. Our results also show an earlier origin of Oleaceae at ca. 86 Mya and the role of Tropical Asia as a main source of species dispersals. CONCLUSION: Our new family-wide and extensive phylogenetic tree highlights both the stable relationships within Oleaceae, including the polyphyly of the genus Chionanthus, and the need for further systematic studies within the family's largest and most under-sampled genera (Chionanthus and Jasminum). Increased sampling will also help to fine-tune biogeographic analyses across spatial scales and geological times.

Challenge accepted: Evolutionary lineages versus taxonomic classification of North American shrub willows (Salix).

PREMISE: The huge diversity of Salix subgenus Chamaetia/Vetrix clade in North America and the lack of phylogenetic resolution within this clade has presented a difficult but fascinating challenge for taxonomists to resolve. Here we tested the existing taxonomic classification with molecular tools. METHODS: In this study, 132 samples representing 46 species from 22 described sections of shrub willows from the United States and Canada were analyzed and combined with 67 samples from Eurasia. The ploidy levels of the samples were determined using flow cytometry and nQuire. Sequences were produced using a RAD sequencing approach and subsequently analyzed with ipyrad, then used for phylogenetic reconstructions (RAxML, SplitsTree), dating analyses (BEAST, SNAPPER), and character evolution analyses of 14 selected morphological traits (Mesquite). RESULTS: The RAD sequencing approach allowed the production of a well-resolved phylogeny of shrub willows. The resulting tree showed an exclusively North American (NA) clade in sister position to a Eurasian clade, which included some North American endemics. The NA clade began to diversify in the Miocene. Polyploid species appeared in each observed clade. Character evolution analyses revealed that adaptive traits such as habit and adaxial nectaries evolved multiple times independently. CONCLUSIONS: The diversity in shrub willows was shaped by an evolutionary radiation in North America. Most species were monophyletic, but the existing sectional classification could not be supported by molecular data. Nevertheless, monophyletic lineages share several morphological characters, which might be useful in the revision of the taxonomic classification of shrub willows.

Phylogenetic analysis and divergence time estimation of Lycium species in China based on the chloroplast genomes.

BACKGROUND: Lycium is an economically and ecologically important genus of shrubs, consisting of approximately 70 species distributed worldwide, 15 of which are located in China. Despite the economic and ecological importance of Lycium, its phylogeny, interspecific relationships, and evolutionary history remain relatively unknown. In this study, we constructed a phylogeny and estimated divergence time based on the chloroplast genomes (CPGs) of 15 species, including subspecies, of the genus Lycium from China. RESULTS: We sequenced and annotated 15 CPGs in this study. Comparative analysis of these genomes from these Lycium species revealed a typical quadripartite structure, with a total sequence length ranging from 154,890 to 155,677 base pairs (bp). The CPGs was highly conserved and moderately differentiated. Through annotation, we identified a total of 128-132 genes. Analysis of the boundaries of inverted repeat (IR) regions showed consistent positioning: the junctions of the IRb/LSC region were located in rps19 in all Lycium species, IRb/SSC between the ycf1 and ndhF genes, and SSC/IRa within the ycf1 gene. Sequence variation in the SSC region exceeded that in the IR region. We did not detect major expansions or contractions in the IR region or rearrangements or insertions in the CPGs of the 15 Lycium species. Comparative analyses revealed five hotspot regions in the CPG: trnR(UCU), atpF-atpH, ycf3-trnS(GGA), trnS(GGA), and trnL-UAG, which could potentially serve as molecular markers. In addition, phylogenetic tree construction based on the CPG indicated that the 15 Lycium species formed a monophyletic group and were divided into two typical subbranches and three minor branches. Molecular dating suggested that Lycium diverged from its sister genus approximately 17.7 million years ago (Mya) and species diversification within the Lycium species of China primarily occurred during the recent Pliocene epoch. CONCLUSION: The divergence time estimation presented in this study will facilitate future research on Lycium, aid in species differentiation, and facilitate diverse investigations into this economically and ecologically important genus.

An expanded phylogeny of the genus Pseudamnicola (Gastropoda; Truncatelloidea; Hydrobiidae) across the Mediterranean Basin.

The genus Pseudamnicola Paulucci, 1878, is commonly found throughout the Mediterranean region. The genus displays considerable levels of endemism, accompanied by notable systematic and taxonomic ambiguity. However, the application of molecular data has proven highly effective in clarifying taxonomy and unveiling the diversity of cryptic species within the genus. Therefore, we employed all cytochrome c oxidase subunit I sequence data available and generated new ones from Greece to infer the phylogeny of the genus throughout its Mediterranean range and estimate the divergence times as well as the ancestral area of diversification. Our phylogenetic and time-estimate analyses demonstrate that with 36 to 38 extant Pseudamnicola species and genetic divergences across species ranging from 0.5% to 11.9% on average, the genus underwent relatively recent diversification during late Miocene (6.53 Ma), and the primary speciation events occurred during Plio-Pleistocene. The Italian Peninsula and Islands and the Ionian Drainages as defined by the Freshwater Ecoregions of the World are the ancestral regions of the genus following two different dispersal routes. Our study contributes to deepening our understanding of Pseudamnicola phylogeny by using data from throughout its range for the first time. This phylogeny provides evidence and confirms previous studies that relatively recent habitat isolation, followed by founder and dispersal events, has been one of the primary reasons for the evolution of the genus Pseudamnicola in the Mediterranean basin.

Sequential Bayesian Phylogenetic Inference.

The ideal approach to Bayesian phylogenetic inference is to estimate all parameters of interest jointly in a single hierarchical model. However, this is often not feasible in practice due to the high computational cost. Instead, phylogenetic pipelines generally consist of sequential analyses, whereby a single point estimate from a given analysis is used as input for the next analysis (e.g., a single multiple sequence alignment is used to estimate a gene tree). In this framework, uncertainty is not propagated from step to step, which can lead to inaccurate or spuriously confident results. Here, we formally develop and test a sequential inference approach for Bayesian phylogenetic inference, which uses importance sampling to generate observations for the next step of an analysis pipeline from the posterior distribution produced in the previous step. Our sequential inference approach presented here not only accounts for uncertainty between analysis steps, but also allows for greater flexibility in software choice (and hence model availability) and can be computationally more efficient than the traditional joint inference approach when multiple models are being tested. We show that our sequential inference approach is identical in practice to the joint inference approach only if sufficient information in the data is present (a narrow posterior distribution) and/or sufficiently many importance samples are used. Conversely, we show that the common practice of using a single point estimate can be biased, e.g., a single phylogeny estimate to transform an unrooted phylogeny into a time-calibrated phylogeny. We demonstrate the theory of sequential Bayesian inference using both a toy example and an empirical case study of divergence-time estimation in insects using a relaxed clock model from transcriptome data. In the empirical example, we estimate three posterior distributions of branch lengths from the same data (DNA character matrix with a GTR+Γ+I substitution model, an amino acid data matrix with empirical substitution models, and an amino acid data matrix with the PhyloBayes CAT-GTR model). Finally, we apply three different node-calibration strategies and show that divergence-time estimates are affected by both the data source and underlying substitution process to estimate branch lengths as well as the node-calibration strategies. Thus, our new sequential Bayesian phylogenetic inference provides the opportunity to efficiently test different approaches for divergence time estimation, including branch-length estimation from other software.