The receptor-like kinase ARK controls symbiotic balance across land plants.

The mutualistic arbuscular mycorrhizal (AM) symbiosis arose in land plants more than 450 million years ago and is still widely found in all major land plant lineages. Despite its broad taxonomic distribution, little is known about the molecular components underpinning symbiosis outside of flowering plants. The ARBUSCULAR RECEPTOR-LIKE KINASE (ARK) is required for sustaining AM symbiosis in distantly related angiosperms. Here, we demonstrate that ARK has an equivalent role in symbiosis maintenance in the bryophyte Marchantia paleacea and is part of a broad AM genetic program conserved among land plants. In addition, our comparative transcriptome analysis identified evolutionarily conserved expression patterns for several genes in the core symbiotic program required for presymbiotic signaling, intracellular colonization, and nutrient exchange. This study provides insights into the molecular pathways that consistently associate with AM symbiosis across land plants and identifies an ancestral role for ARK in governing symbiotic balance.

Plastid genome evolution in leafless members of the orchid subfamily Orchidoideae, with a focus on Degranvillea dermaptera.

PREMISE: Leafless, heterotrophic plants are prime examples of organismal modification, the genomic consequences of which have received considerable interest. In particular, plastid genomes (plastomes) are being sequenced at a high rate, allowing continual refinement of conceptual models of reductive evolution in heterotrophs. However, numerous sampling gaps exist, hindering the ability to conduct comprehensive phylogenomic analyses in these plants. METHODS: Using floral tissue from an herbarium specimen, we sequenced and analyzed the plastome of Degranvillea dermaptera, a rarely collected, leafless orchid species from South America about which little is known, including its phylogenetic affinities. RESULTS: The plastome is the most reduced of those sequenced among the orchid subfamily Orchidoideae. In Degranvillea, it has lost the majority of genes found in leafy autotrophic species, is structurally rearranged, and has similar gene content to the most reduced plastomes among the orchids. We found strong evidence for the placement of Degranvillea within the subtribe Spiranthinae using models that explicitly account for heterotachy, or lineage-specific evolutionary rate variation over time. We further found evidence of relaxed selection on several genes and of correlations among substitution rates and several other "traits" of the plastome among leafless members of orchid subfamily Orchidoideae. CONCLUSIONS: Our findings advance knowledge on the phylogenetic relationships and paths of plastid genome evolution among the orchids, which have experienced more independent transitions to heterotrophy than any other plant family. This study demonstrates the importance of herbarium collections in comparative genomics of poorly known species of conservation concern.

Detecting diversifying selection for a trait from within and between-species genotypes and phenotypes.

To quantify selection acting on a trait, methods have been developed using either within or between-species variation. However, methods using within-species variation do not integrate the changes at the macroevolutionary scale. Conversely, current methods using between-species variation usually discard within-species variation, thus not accounting for processes at the micro-evolutionary scale. The main goal of this study is to define a neutrality index for a quantitative trait, by combining withinand between-species variation. This neutrality index integrates nucleotide polymorphism and divergence for normalizing trait variation. As such, it does not require estimation of population size nor of time of speciation for normalization. Our index can be used to seek deviation from the null model of neutral evolution, and test for diversifying selection. Applied to brain mass and body mass at the mammalian scale, we show that brain mass is under diversifying selection. Finally, we show that our test is not sensitive to the assumption that population sizes, mutation rates and generation time are constant across the phylogeny, and automatically adjust for it.

First Insight into the Phylogenetic Diversity of Bovicola caprae Infesting Goats of Different Agro-climatic Locations in India.

Bovicola caprae is an important obligate ectoparasite of goats worldwide including India. The present study aimed at the molecular confirmation, phylogenetics and population structure analyses of B. caprae infesting goats of three different agro-climatic locations in India, by targeting the mitochondrial cytochrome C oxidase subunit 1 (cox1) genetic marker. The phylogenetic tree exhibited the presence of two different lineages of B. caprae. The sequences generated herein clustered in lineage 2 along with the GenBank™ archived sequences from China and Iran. The sequences generated herein also showed the circulation of sub-lineages of B. caprae in India based on the analysis of pairwise genetic distances between sequences and median-joining haplotype network. The population structure analyses revealed low nucleotide (0.00353 ± 0.00291 and 0.02694 ± 0.00363) and high haplotype (0.667 ± 0.314 and 0.618 ± 0.104) diversities for the present study isolates as well as for the complete dataset, respectively, which evinced a recent demographic expansion. High genetic differentiation (FST value = 0.97826) and low gene flow (Nm = 0.00556) were also recorded in the different lineages/populations. In conclusion, the present study addressed the research gap and provided the first insight into the phylogenetics of the goat louse B. caprae and highlighted the circulation of sub-lineages of the ectoparasite in India.

Cladistics of Rhipicephalus microplus and laboratory assessment of acaricidal efficacy of different Cassia fistula extracts against R. microplus larvae.

The present study aimed to analyze the cladistics and population structure analysis of Rhipicephalus microplus ticks infesting buffaloes in Haryana, India, as well as the assessment of the anti-tick efficacy of the ethanolic extracts of Cassia fistula (bark, pod pulp, and flowers) against R. microplus larvae. The molecular characterization and population structure analysis were performed by targeting the amplification of the partial mitochondrial cytochrome C oxidase subunit 1 (cox1) gene, whereas anti-tick efficacy was evaluated using a larval packet test. The sequences generated herein revealed a homology of 98.26-100% to the GenBank-archived R. microplus sequences. In population structure analysis, high haplotype (0.500 ± 0.265) and low nucleotide (0.002 ± 0.001) diversities were recorded for the sequences generated in this study. Significantly negative neutrality indices were recorded for the overall dataset. The extracts were found to significantly affect mortality rates in a dose-dependent manner, and the ethanolic extracts of the bark, pod pulp, and flowers of C. fistula exhibited median lethal concentration (LC50) values of 27.989, 40.457, and 49.43 mg/mL, respectively. The LC50 value recorded for the combination of the ethanolic extracts of the bark, flower, and pod pulp of C. fistula was 19.724 mg/mL, whereas the synthetic acaricide ivermectin had an LC50 value of 351.56 mg/mL. In conclusion, R. microplus populations infesting cattle and buffalo hosts in India exhibited negligible genetic differentiation and high gene flow between them. Additionally, the combination of all C. fistula extracts could serve as a potential substitute for the synthetic acaricide.

Phylogenetic assessment within a species complex of a subterranean rodent (Geomys bursarius) with conservation implications for isolated subspecies.

Range contraction and expansion from glaciation has led to genetic divergence that may be particularly pronounced in fossorial species with low dispersal. The plains pocket gopher (Geomys bursarius) is a fossorial species that ranges widely across North America but has a poorly understood phylogeny. We used mitogenomes (14,996 base pairs) from 56 individuals across seven subspecies, plus two outgroup species, to assess genetic divergence from minimum spanning trees, measure genetic distances, and infer phylogenetic trees using BEAST. We found G. b. wisconsinensis was monophyletic with recent divergence. Further assessment is needed for G. b. major because it was paraphyletic and exhibited inconsistent groupings with other clades. Importantly, we identified G. b. illinoensis as being genetically distinct and monophyletic likely due to a unique colonization event eastward across the Mississippi River. Because G. b. illinoensis faces continued pressures from niche reduction and habitat loss, we recommend that G. b. illinoensis be considered an evolutionary significant unit warranting conservation actions to promote connectivity and restore suitable habitat. Such conservation efforts should benefit other grassland species including those originating from clades west of the Mississippi River that may also be evolutionary significant units.

HIV Subtypes and Drug-resistance-associated Mutations in US Blood Donors, 2015-2020.

Background: Monitoring genotypes of HIV infections in blood donors may provide insights into infection trends in the general population. Methods: HIV RNA was extracted from plasma samples of blood donors confirmed as HIV positive by blood screening nucleic acid and antibody tests. HIV genome target regions were amplified using nested real time-polymerase chain reaction followed by next-generation sequencing. Sequences were compared to those in the Los Alamos National Laboratory (LANL) database. Sequences were also assessed for drug resistance mutations (DRM) using the Stanford HIV DRM Database. Results: From available HIV-positive donations collected between 1 September 2015 and 31 December 2020, 563 of 743 (75.8%) were successfully sequenced; 4 were subtype A, 543 subtype B, 5 subtype C, 1 subtype G, 5 circulating recombinant forms (CRF), and 2 were subtype B and D recombinants. Overall, no significant differences between blood donor and available LANL genotypes were found, and the genotypes of newly acquired versus prevalent HIV infections in donors were similar. The proportion of non-B subtypes and CRF remained a small fraction, with no other subtype or CRF representing more than 1% of the total. DRM were identified in 122 (21.6%) samples with protease inhibitor, nucleoside reverse transcriptase inhibitor and non-nucleoside reverse transcriptase inhibitor DRMs identified in 4.9%, 4.6% and 14.0% of samples, respectively. Conclusions: HIV genetic diversity and DRM in blood donors appear representative of circulating HIV infections in the US general population and may provide more information on infection diversity than sequences reported to LANL, particularly for recently transmitted infections.

Phylogenetic reconciliation: making the most of genomes to understand microbial ecology and evolution.

In recent years, phylogenetic reconciliation has emerged as a promising approach for studying microbial ecology and evolution. The core idea is to model how gene trees evolve along a species tree, and to explain differences between them via evolutionary events including gene duplications, transfers, and losses. Here, we describe how phylogenetic reconciliation provides a natural framework for studying genome evolution, and highlight recent applications including ancestral gene content inference, the rooting of species trees, and the insights into metabolic evolution and ecological transitions they yield. Reconciliation analyses have elucidated the evolution of diverse microbial lineages, from Chlamydiae to Asgard archaea, shedding light on ecological adaptation, host-microbe interactions, and symbiotic relationships. However, there are many opportunities for broader application of the approach in microbiology. Continuing improvements to make reconciliation models more realistic and scalable, and integration of ecological metadata such as habitat, pH, temperature and oxygen use, offer enormous potential for understanding the rich tapestry of microbial life.

Characterization of the family-level Borreliaceae pan-genome and development of an episomal typing protocol.

Background: The Borreliaceae family includes many obligate parasitic bacterial species which are etiologically associated with a myriad of zoonotic borrelioses including Lyme disease and vector-borne relapsing fevers. Infections by the Borreliaceae are difficult to detect by both direct and indirect methods, often leading to delayed and missed diagnoses. Efforts to improve diagnoses center around the development of molecular diagnostics (MDx), but due to deep tissue sequestration of the causative spirochaetes and the lack of persistent bacteremias, even MDx assays suffer from a lack of sensitivity. Additionally, the highly extensive genomic heterogeneity among isolates, even within the same species, contributes to the lack of assay sensitivity as single target assays cannot provide universal coverage. This within-species heterogeneity is partly due to differences in replicon repertoires and genomic structures that have likely arisen to support the complex Borreliaceae lifecycle in which these parasites have to survive in multiple hosts each with unique immune responses. Results: We constructed a Borreliaceae family-level pangenome and characterized the phylogenetic relationships among the constituent taxa which supports the recent taxonomy of splitting the family into at least two genera. Gene content pro les were created for the majority of the Borreliaceae replicons, providing for the first time their unambiguous molecular typing. Conclusion: Our characterization of the Borreliaceae pan-genome supports the splitting of the former Borrelia genus into two genera and provides for the phylogenetic placement of several non-species designated isolates. Mining this family-level pangenome will enable precision diagnostics corresponding to gene content-driven clinical outcomes while also providing targets for interventions.

Maximum Likelihood Estimation for Unrooted 3-Leaf Trees: An Analytic Solution for the CFN Model.

Maximum likelihood estimation is among the most widely-used methods for inferring phylogenetic trees from sequence data. This paper solves the problem of computing solutions to the maximum likelihood problem for 3-leaf trees under the 2-state symmetric mutation model (CFN model). Our main result is a closed-form solution to the maximum likelihood problem for unrooted 3-leaf trees, given generic data; this result characterizes all of the ways that a maximum likelihood estimate can fail to exist for generic data and provides theoretical validation for predictions made in Parks and Goldman (Syst Biol 63(5):798-811, 2014). Our proof makes use of both classical tools for studying group-based phylogenetic models such as Hadamard conjugation and reparameterization in terms of Fourier coordinates, as well as more recent results concerning the semi-algebraic constraints of the CFN model. To be able to put these into practice, we also give a complete characterization to test genericity.

An updated LSU database and pipeline for environmental DNA identification of arbuscular mycorrhizal fungi.

Recent work established a backbone reference tree and phylogenetic placement pipeline for identification of arbuscular mycorrhizal fungal (AMF) large subunit (LSU) rDNA environmental sequences. Our previously published pipeline allowed any environmental sequence to be identified as putative AMF or within one of the major families. Despite this contribution, difficulties in implementation of the pipeline remain. Here, we present an updated database and pipeline with (1) an expanded backbone tree to include four newly described genera and (2) several changes to improve ease and consistency of implementation. In particular, packages required for the pipeline are now installed as a single folder (conda environment) and the pipeline has been tested across three university computing clusters. This updated backbone tree and pipeline will enable broadened adoption by the community, advancing our understanding of these ubiquitous and ecologically important fungi.

Integrons and multidrug resistance across phylogenetic groups of clinical isolates of Escherichia coli.

Objective: This study was aimed to investigate the multidrug resistance patterns in clinical isolates of Escherichia coli and their correlation with integrons and phylogenetic groupings. Methods: A total of 37 clinical E. coli isolates were evaluated for drug resistance patterns by disk diffusion method. Phylogenetic groupings and the presence of integrons among E. coli were determined by multiplex PCR assays. Results: Multidrug resistance was identified in 84% of the clinical isolates of E. coli with higher resistance found against cephalosporins (94.6%) and fluoroquinolones (83.8%), while lower resistance was observed against polymyxins (24.3%) and carbapenems (29.7%). Metallo-ß-lactamases were found in all carbapenem resistant isolates. The phylogenetic group B2 was the most dominant (40.5%), followed by groups A (35.1%), D (13.5%) and B1 (10.8%). Integrons were detected in 25 (67.6%) isolates and intI1, intI2, and intI3 genes were found in 62.2%, 18.9% and 10.8% of isolates respectively. Conclusion: Our results show that phylogenetic classification of E. coli is not relevant with antimicrobial resistance. However, there was strong association between the integron classes and resistance against ß-lactam and fluoroquinolones antimicrobials. Additionally, this study highlighted that the presence of integrons plays a crucial role in the development of multidrug resistance in clinical isolates of E. coli. Most significantly, this is the first report of detection of three classes of integron among clinical isolates of E. coli in Pakistan.

Sources of Human Immunodeficiency Virus Infections Among Men Who Have Sex With Men With a Migration Background: A Viral Phylogenetic Case Study in Amsterdam, The Netherlands.

BACKGROUND: Men and women with a migration background comprise an increasing proportion of incident human immunodeficiency virus (HIV) cases across Western Europe. METHODS: To characterize sources of transmission in local transmission chains, we used partial HIV consensus sequences with linked demographic and clinical data from the opt-out AIDS Therapy Evaluation in the Netherlands (ATHENA) cohort of people with HIV in the Netherlands and identified phylogenetically and epidemiologically possible HIV transmission pairs in Amsterdam. We interpreted these in the context of estimated infection dates, and quantified population-level sources of transmission to foreign-born and Dutch-born Amsterdam men who have sex with men (MSM) within Amsterdam transmission chains. RESULTS: We estimate that Dutch-born MSM were the predominant sources of infections among all Amsterdam MSM who acquired their infection locally in 2010-2021, and among almost all foreign-born Amsterdam MSM subpopulations. Stratifying by 2-year intervals indicated time trends in transmission dynamics, with a majority of infections originating from foreign-born MSM since 2016, although uncertainty ranges remained wide. CONCLUSIONS: Native-born MSM have predominantly driven HIV transmissions in Amsterdam in 2010-2021. However, in the context of rapidly declining incidence in Amsterdam, the contribution from foreign-born MSM living in Amsterdam is increasing, with some evidence that most local transmissions have been from foreign-born Amsterdam MSM since 2016.

Embedding of Markov matrices for d ⩽ 4.

The embedding problem of Markov matrices in Markov semigroups is a classic problem that regained a lot of impetus and activities through recent needs in phylogeny and population genetics. Here, we give an account for dimensions d ⩽ 4 , including a complete and simplified treatment of the case d = 3 , and derive the results in a systematic fashion, with an eye on the potential applications. Further, we reconsider the setup of the corresponding problem for time-inhomogeneous Markov chains, which is needed for real-world applications because transition rates need not be constant over time. Additional cases of this more general embedding occur for any d ⩾ 3 . We review the known case of d = 3 and describe the setting for future work on d = 4 .

Prorocentrum canariense sp. nov., a case of pseudo-cryptic speciation in the cosmopolitan dinoflagellate P. compressum (Prorocentrales, Dinophyceae).

The planktonic dinoflagellate Prorocentrum compressum is widespread in warm and temperate seas. A strain identified as P. cf. compressum BEA 0681B isolated from the island of Gran Canaria, NE Atlantic Ocean, showed a divergence in rDNA/ITS phylogenies with respect to P. compressum. The Canarian strain was oval, with an average length-to-width ratio of 1.35, smooth thecal surface with less than 150 thecal pores, including oblique pores, sometimes with a bifurcated opening. In contrast, P. compressum was rounder, with a length-to-width ratio < 1.2, with reticulate-foveate ornamentation and 200-300 pores per valve. We propose Prorocentrum canariense sp. nov. These species clustered as the most early-branching lineage in the clade Prorocentrum sensu stricto. Although this clade mainly contains planktonic species, the closer relatives were the benthic species P. tsawwassenense and P. elegans. Interestingly, P. compressum and P. canariense sp. nov. are widely distributed in temperate and warm seas without an apparent morphological adaptation to planktonic life. The formation of two concentric hyaline mucilaginous walls could contribute to this success. We discuss the use of Prorocentrum bidens to solve the nomenclature issue of P. compressum that was described citing a diatom as basionym.

Mitonuclear compatibility is maintained despite relaxed selection on male mitochondrial DNA in bivalves with doubly uniparental inheritance.

Mitonuclear coevolution is common in eukaryotes, but bivalve lineages that have doubly uniparental inheritance (DUI) of mitochondria may be an interesting example. In this system, females transmit mtDNA (F mtDNA) to all offspring, while males transmit a different mtDNA (M mtDNA) solely to their sons. Molecular evolution and functional data suggest oxidative phosphorylation (OXPHOS) genes encoded in M mtDNA evolve under relaxed selection due to their function being limited to sperm only (vs. all other tissues for F mtDNA). This has led to the hypothesis that mitonuclear coevolution is less important for M mtDNA. Here, we use comparative phylogenetics, transcriptomics, and proteomics to understand mitonuclear interactions in DUI bivalves. We found nuclear OXPHOS proteins coevolve and maintain compatibility similarly with both F and M mtDNA OXPHOS proteins. Mitochondrial recombination did not influence mitonuclear compatibility and nuclear-encoded OXPHOS genes were not upregulated in tissues with M mtDNA to offset dysfunction. Our results support that selection maintains mitonuclear compatibility with F and M mtDNA despite relaxed selection on M mtDNA. Strict sperm transmission, lower effective population size, and higher mutation rates may explain the evolution of M mtDNA. Our study highlights that mitonuclear coevolution and compatibility may be broad features of eukaryotes.

Epidemiological, virological and clinical characterization of a Dengue/Zika outbreak in the Caribbean region of Costa Rica 2017-2018.

The increase in incidence and geographical expansion of viruses transmitted by the Aedes mosquitoes, such as dengue (DENV) and zika (ZIKV) in the Americas, represents a burden for healthcare systems in tropical and subtropical regions. These and other under-detected arboviruses co-circulate in Costa Rica, adding additional complexity to their management due to their shared epidemiological behavior and similarity of symptoms in early stages. Since diagnostics of febrile illness is mostly based on clinical symptoms alone, we gathered acute-phase serum and urine from 399 samples of acute dengue-like cases from two healthcare facilities of Costa Rica, during an outbreak of arboviruses from July 2017 to May 2018, and tested them using molecular and serological methods. The analyses showed that of the clinically presumptive arbovirus cases that were reported, only 39.4% (n=153) of the samples were confirmed positive by RT-PCR to be DENV (DENV (10.3%), CHIKV (0.2%), ZIKV (27.3%), or mixed infections (1.5%). RT-PCR for other alphaviruses and flaviviruses, and PCR for Leptospira sp were negative. Furthermore, to assess flavivirus positivity in post-acute patients, the negative sera were tested against Dengue-IgM. 20% of sera were found positive, confounding even more the definitive number of cases, and emphasizing the need of several distinct diagnostic tools for accurate diagnostics. Molecular characterization of the prM and E genes from isolated viruses revealed that the American/Asian genotype of DENV-2 and the Asian lineage of ZIKV were circulating during this outbreak. Two different clades of DENV-2 American/Asian genotype were identified to co-circulate in the same region and a difference in the platelet and leukocyte count was noted between people infected with each clade, suggesting a putative distinct virulence. Our study sheds light on the necessity for healthcare strategies in managing arbovirus outbreaks, emphasizing the importance of comprehensive molecular and serological diagnostic approaches, as well as molecular characterization. This approach aids in enhancing our understanding of the clinical and epidemiological aspects of arboviral diseases during outbreaks. Our research highlights the need to strengthen training programs for health professionals and the need to increase research-based on laboratory evidence for diagnostic accuracy, guidance, development and implementation of public health interventions and epidemiological surveillance.

Genetic diversity of porcine circoviruses 2 and 3 circulating among wild boars in the Moscow Region of Russia.

Porcine circoviruses (PCVs) are widely distributed in swine herds. PCV2, the significant swine pathogen, causes infections characterized by growth and development disorders, skin lesions, and respiratory distress. PCV3 has been circulating worldwide and can be associated with various clinical signs and disease developments. Wild boars are the main reservoir of these pathogens in wildlife and can create an alarming threat to pig farming. In Russia, three PCV2 genotypes (PCV2a, PCV2b, and PCV2d) were identified in pig farms. Additionally, PCV3 was observed in pig herds during the monitoring studies in the country. However, data considering the circulation of PCVs in herds of wild boars in Russia is scant. For this purpose, we performed PCR assays of the samples from 30 wild boars hunted in the Moscow Region of Russia in 2021-2023. The ratios of wild boars positive for PCV2, PCV3, or coinfected were 50, 10, and 13.3%, respectively. Additionally, we sequenced 15 PCV2 and four PCV3 complete genomes and conducted phylogenetic analysis, which divided PCV2 isolates into two groups: PCV2d and PCV2b. The study showed a high infection rate of PCV2 among wild boars, with PCV2d dominance. Simultaneously, PCV3 also circulates among wild boars. The obtained results can provide a basis for the development of preventive measures to support infection transmission risks between farm and wild animals.

Utility of various genome lengths in diversity and evolution analyses of Hepatitis E virus.

The aim of this study was to investigate to what extent fragments of the HEV genome could be used for accurate diagnostics and inference of viral population-scale processes. For this, we selected all the published whole genome sequences from the NCBI GenBank and trimmed them to various fragment lengths (ORF1,2,3, ORF1, ORF2, ORF3, 493 nt in ORF2 and 148 nt in ORF2). Each of the fragment lengths was used to infer the richness and diversity of the viral sequence types, typing accuracy, and potential use in phylodynamics. The results obtained from the different fragments were compared. We observed that, generally, the longer the nucleic acid fragment used in typing, the better the accuracy in predicting the viral subtype. However, the dominant HEV subtypes circulating in Europe were relatively well classified even by the 493 nt fragment, with false negative rates as low as 8 in 1000 typed sequences. Most fragments also give comparable results in analyses of population size, albeit with shorter fragments showing a broader 95 % highest posterior density interval and less obvious increase of the viral effective population size. The reconstructed phylogenies of a heterochronous subset indicated a good concordance between all the fragments, with the major clades following similar branching patterns. Furthermore, we have used the HEV sequence data from the Netherlands available in the HEVnet database as a case study for reconstruction of population size changes in the past decades. This data showed that molecular and epidemiological results are concordant and point to an increase in the viral effective population size underlying the observed increase in incidence of acute HEV infection cases. In the absence of whole genome sequencing data, the 493 bp fragment can be used for analyzing HEV strains currently circulating in Europe, as it is informative for describing short term population-scale processes.

The polyphyletic Caucasus-centred Campanula subg. Scapiflorae (Campanulaceae) revisited with a newly circumscribed C. sect. Tridentatae for its core clade.

Campanula L. is among the genera with the highest number of endemics in the Caucasus ecoregion. A group of attractive alpine and subalpine perennial rosette plants with short single-flowered stems centred in the Caucasus has been treated as Campanulasubg.Scapiflorae or at other ranks, with considerably varying circumscription and classification. Molecular phylogenetic analysis of three plastid DNA regions (trnK/matK, petD, rpl16) of a strongly extended sampling, comprising 23 of the 27 commonly accepted taxa (85%) with 330 accessions built on and guided by the results of our previous study of the group, confirmed the polyphyly of C.subg.Scapiflorae in any of its circumscriptions. The core clade of the group comprises exclusively endemics and near-endemics of the Caucasus and is treated here as C.sect.Tridentatae in a revised circumscription. The phylogenetic relationships of the disparate other elements of the Scapiflorae group are outlined.