Corrigendum to "Reticulate evolution, introgression, and recent diversification in Epimedium sect. Macroceras" [Molecul. Phylog. Evol. 223 (2026) 108646].
Explore the source record for details and available documents.
SEARCH · Search PubMed
Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
The Morelloid clade (black nightshades) is one of the most strongly supported clades within the megadiverse Solanum genus. It comprises 76 globally distributed, non-spiny herbaceous and suffrutescent species. While often erroneously considered poisonous weeds, several species are economically important as orphan crops. The clade is closely related to tomato and potato but, due to a lack of focused breeding efforts, remains a putative reservoir of genetic diversity for crop improvement. Despite this potential, we lack fundamental knowledge on the evolution of the Morelloid clade. The group includes polyploid species with unknown parental origins-likely reflecting reticulate processes such as hybridization, introgression, and associated backcrossing events. Prior analyses have been unable to disentangle these processes, leaving the mechanisms underlying reticulate evolution in the Morelloid clade poorly understood. Here, we use genome skimming to produce a well-supported maximum likelihood plastid phylogeny from complete circularized plastomes and a coalescent-based species tree from combined Angiosperms353 and conserved ortholog set nuclear markers. Our dataset, composed of previously published data and deep genome skimming from herbarium samples, spans 26 Morelloid species. To investigate phylogenetic discordance, we used a nuclear phylogenetic network, multispecies coalescent simulations, a fused rooted nuclear chloroplast tree, and quantification of nuclear gene tree concordance. We show that incongruence between nuclear and plastid trees is pervasive and cannot be explained by incomplete lineage sorting alone. Instead, our results demonstrate that events consistent with repeated chloroplast capture have shaped the reticulate evolutionary history of the clade, especially among African polyploid and Pan-American diploid lineages.
Understanding how lineages diversify despite persistent ancestral polymorphism and recurrent gene flow remains a central challenge in evolutionary biology. Juniperus distributed across the Qinghai-Tibet Plateau provide an ideal system for addressing this question because repeated geological uplift and climatic oscillations have likely promoted cycles of lineage divergence, range shifts, and secondary contact. Here, we combined approximately 1.08 million genome-wide SNPs from 164 individuals representing thirteen Juniperus lineages with phylogenomic datasets comprising 3,381 nuclear single-copy genes and nearly complete plastomes. We detected extensive phylogenomic discordance and cytonuclear incongruence across genomic datasets. Topology weighting, coalescent simulations, quartet-based tests, and analyses of gene flow and reticulation collectively support the interpretation that these patterns were shaped by the combined effects of prolonged incomplete lineage sorting and gene flow during lineage diversification. Ecological niche analyses further provide a spatial and climatic context in which environmentally similar lineages may have had greater opportunities for secondary contact during historical range shifts. Collectively, our results reveal that the evolutionary history of Qinghai-Tibet Plateau Juniperus is characterized by reticulate diversification rather than strictly bifurcating evolution, and demonstrate how genome-wide discordance can provide biological insights into the evolutionary processes underlying lineage diversification.
Diervilleae (Caprifoliaceae) is a horticulturally important lineage with striking floral diversity and a long history of interspecific crossing, suggesting reticulate evolution. We integrated nuclear SNPs and whole plastome data to reconstruct a phylogenomic backbone for the tribe and to identify hybrids, cultivated accessions, and introgression among lineages. Nuclear and plastid phylogenies consistently recover Weigela and Diervilla as reciprocally monophyletic and resolve four major lineages within Weigela, providing a reproducible framework for revising sectional limits and species boundaries. Cultivated accessions form a well supported clade sister to W. florida and show predominantly W. florida ancestry while retaining contributions from multiple wild lineages, consistent with recurrent crossing, backcrossing, and selection. Analyses of wild populations reveal recurrent hybrids and enable plausible parental combinations to be inferred. Tests across the genome further indicate strong evidence for historical introgression across Diervilleae, with the strongest signals involving W. middendorffiana, W. maximowiczii, and Diervilla. Fossil evidence, divergence time estimation, and paleodistribution modelling together suggest range expansion during the Miocene and Pliocene followed by climate driven contraction, providing a spatiotemporal context for episodic contact, introgression, and the East Asia-North America disjunction.
Paramecium bursaria maintains a stable endosymbiosis with green algae, yet the evolutionary consequences of this association remain unclear. Here, we screened the host genome for algal-derived horizontally transferred genes (HTGs) using a lineage-aware workflow designed to detect horizontal gene transfer (HGT) between two defined lineages. We identified 16 candidate HTGs, including four putative newly transferred genes and 12 homologous transferred genes, most of which were functionally associated with redox homeostasis and metabolism. Five HTGs showed symbiosis-dependent expression. RNAi knockdown of GH32s and SATs reduced host proliferation, total cell area, and motility, while GH32s knockdown also reduced endosymbiont load. Duplication patterns suggest that most transfers may have occurred after the P. bursaria lineage diverged from the sampled Paramecium species but before its lineage-specific whole-genome duplication (WGD). The HTGs also showed host-associated shifts in GC content and gene length, while representative HTGs retained conserved domains and functional motifs. Together, our results support algae-to-host HGT in P. bursaria and suggest that some transferred genes may contribute to metabolic integration during endosymbiosis.
The angiosperm order Rosales still represents a major challenge for phylogenetic reconstruction. Although its circumscription is now well-defined, phylogenetic relationships among families are still uncertain. Here, we used nuclear, plastid, and mitochondrial genomic data from 33 species representing all nine families to further clarify interfamilial relationships and the group's evolutionary history. We detected significant phylogenetic conflict among the three datasets. Further analyses at the nuclear level identified incomplete lineage sorting (ILS) as the main cause of unstable phylogenetic positions among families. The discordant placements of Rhamnaceae and Elaeagnaceae based on plastid and mitochondrial data are caused by ancient hybridization events, potentially involving differences in organellar inheritance. Our molecular dating confirms earlier suggestions that the ancient rapid diversification of the three Rosaceae subfamilies could be the main reason for the difficulties in resolving their phylogenetic relationships. Our findings provide new insights into the interfamilial relationships of Rosales and demonstrate that the evolutionary history of this order was shaped by ancient and rapid radiation as well as extensive ILS and reticulate evolution. They also suggest that previous attempts to clarify interfamilial relationships in this order were hampered by combining nuclear and organellar sequence data, leading to inconsistent topologies observed across earlier studies.
The olfactory receptor (OR) genes constitute the molecular basis of fish olfaction, mediating survival behaviors and environmental adaptation while coevolving with habitat-driven evolution. Schizothorax, a cyprinid genus endemic to the Qinghai-Tibetan Plateau, exhibits remarkable dietary divergence and ploidy variation in response to plateau environmental changes, which presumably facilitates the adaptive evolution of OR genes. However, the evolutionary patterns of OR genes associated with trophic divergence and ploidy variation in this genus remain unclear. In this study, three species were selected: the herbivorous diploid S. macropogon, the carnivorous diploid S. lantsangensis, and the herbivorous tetraploid S. curvilabiatus. S. macropogon possessed 142 OR genes (92.25% functional), primarily located on chromosomes 14 and 24, with the fewest sequence clusters. Such compact gene repertoire and highly overlapping chromosomal clusters indicated specialization for a herbivorous olfactory niche. S. lantsangensis contained 127 OR genes (93.70% functional), concentrated on chromosomes 4 and 5, with fewer sequence clusters and a scattered distribution, reflecting evolution of OR genes under carnivorous feeding habits. The herbivorous tetraploid S. curvilabiatus exhibited striking features: 316 OR genes (94.30% functional), the most subfamilies, unique ε and κ OR subfamilies, and species-specific motifs. These characteristics revealed that ploidy, rather than herbivory, dominated OR gene evolution. In conclusion, dietary differentiation and ploidy variation together drove olfactory adaptive evolution in Schizothorax, providing new insights into vertebrate OR gene ecological adaptation.
Antipatharians (black corals) are among the least studied coral groups, with much of their diversity still undescribed. Here, we present an integrative morphological, phylogenomic and genomic distance study of deep-sea antipatharians sampled in high seas areas of the North Pacific Ocean and from New Zealand's Exclusive Economic Zone. These corals grow on hexactinellid sponges - a unique characteristic in the order Antipatharia. Using a dataset of ultra-conserved elements and exons, combined with morphological analyses, we reconstruct phylogenomic relationships and formally describe a new family (Eidikopathidae fam. nov.), a new genus (Eidikopathesgen. nov.), and two new species (E. korallispongiasp. nov., E. zealandkoralliasp. nov.). Morphologically, the new family is distinguished by a corallum consisting of a network of loose branches that fuse with the sponge skeletal framework. Phylogenomic analyses recovered consistent topologies with strong nodal support, corroborating the distinct evolutionary placement of this sponge-associated lineage. Pairwise genomic distances estimated using the Tamura-Nei model were concordant with patristic genomic distances, identifying Pteridopathidae as the genetically closest family to Eidikopathidae fam. nov., followed by Myriopathidae and Stylopathidae, which were recovered as sister families in the phylogeny. This pattern shows that genomic distance complements, rather than simply mirrors, tree topology by quantifying accumulated sequence divergence among lineages. Together, these results provide the first genomic distance framework for Antipatharia, offering a baseline for future systematic, evolutionary, and biodiversity studies on this fundamental shallow, mesophotic and deep-sea coral group.
The mechanisms driving lineage divergence in red algae remain unexplored, despite the group's remarkable diversity and ancient evolutionary history. The red alga Mazzaella laminarioides, a Chilean intertidal species complex composed of three parapatric cryptic lineages (North, Center, South), offers a valuable system to evaluate these processes, as its life history combines severe dispersal limitation with a haploid-diploid cycle that may influence the emergence of reproductive barriers. We reconstructed its evolutionary history using whole-genome sequencing and nuclear genome assembly of representative individuals from each lineage. Phylogenomic analyses based on 1,507 single-copy orthologs recovered three deeply divergent lineages with limited nuclear discordance consistent with incomplete lineage sorting. For both splits, demographic modelling was most consistent with an Ancient Migration scenario, although support over strict isolation was moderate, suggesting that divergence may have begun with low asymmetric ancestral gene flow followed by subsequent loss of connectivity, demographic bottlenecks, and later population expansion. Coding sequence analyses revealed lineage-specific dN/dS heterogeneity; only one South-lineage locus passed FDR correction (metaxin-1, mitochondrial protein import), with two further South-lineage candidates in chlorophyll and heme biosynthesis falling below the FDR threshold. Together, these signals suggest that divergent selective pressures on energy acquisition may have contributed to divergence at the southern end of the distribution. These results add to the small but growing body of whole-genome data for red algae and, alongside recent macroalgal studies, suggest that ancestral connectivity could be a recurrent feature of lineage divergence even in marine organisms with extremely restricted dispersal.
Fine-resolution information on species relationships and biological diversity is critically needed to guide conservation efforts amidst rapid environmental changes. Systematics, which forms the foundation of this knowledge, has been revolutionized by phylogenomics, utilizing genome-scale datasets. However, the use of diverse marker types, non-comparable taxon sampling, and outgroup selection can lead to conflicting phylogenetic hypotheses. These inconsistencies complicate study comparisons and hinder our ability to assess marker-specific impacts on phylogenetic resolution. The phylogenetic reconstruction of the bat genus Myotis, encompassing over 140 species and characterized by a rapid radiation in the last 20 million years, has been particularly influenced by these challenges. Achieving phylogenetic resolution in Myotis is particularly complex due to subtle interspecific differences in both morphological and molecular traits. Mitochondrial and nuclear markers often produce discordant trees, influenced by hybridization, introgression, and methodological variations. In this study, we employed a consistent taxonomic sample set of 44 Myotis taxa to evaluate the impact of five different genetic marker types on phylogenetic reconstruction. We observed significant discordance between topologies derived from conserved nuclear and mitochondrial markers and found that transposable elements were inadequate for resolving relationships across the entire genus. Our results also clarify the placement of previously problematic taxa within the genus. These findings emphasize the importance of aligning genetic marker choice with specific phylogenetic questions and highlight the influence of taxonomic and methodological variation on phylogenomic outcomes. This work provides a framework for improving phylogenetic inference in rapidly radiating groups and enhances our understanding of evolutionary history in Myotis.
Phylogenomic discordance is widespread across plants, but its evolutionary significance is often obscured when conflict is treated primarily as analytical noise rather than as evidence of underlying processes. In woody lineages in particular, incomplete lineage sorting, introgression, and genome duplication can interact over long timescales to produce complex genomic histories that are not adequately summarized by a strictly bifurcating tree. Here, we use Malus as a model woody genus to investigate how these processes structure conflict across a genus-scale, accession-based phylogenomic framework. Using broad taxon sampling, hundreds of nuclear loci, plastid genomes, and genome-wide SNP summaries, we reconstruct a robust nuclear backbone for sampled Malus lineages and evaluate where discordance is concentrated and which processes best explain it. Nuclear analyses resolve eight major clades, whereas conflict is non-random and localized to recurrent hotspots rather than evenly distributed across the tree. Cytonuclear discordance is similarly concentrated, especially around Clade H, represented by sampled accessions of M. tschonoskii, where localized plastid-nuclear disagreement is consistent with candidate plastid capture or organellar introgression. Multiple complementary analyses further indicate that the strongest conflict is not explained by ILS alone, but instead reflects lineage-structured introgression, while polyploid complexes represent additional localized sources of evolutionary complexity. Together, these results provide evidence for a reticulate genomic backbone in Malus and show how integrating nuclear, plastid, and genome-wide conflict analyses can help distinguish background discordance from process-specific signals in woody plant radiations. Several lineage-level reticulation hypotheses identified here should now be tested with broader population-level sampling and curated reference accessions.
The generic-level classification of the bee tribe Megachilini (Megachilidae) has remained controversial due to poor phylogenetic resolution at the base of the group, particularly among the brood parasitic genera and the numerous dauber ("Chalicodoma s. l.") lineages. We present a phylogenomic analysis of Megachilini based on ultraconserved elements (UCEs), sampling 52 ingroup taxa with emphasis on the dauber lineages. We also present a combined UCE + six-gene analysis to improve taxon coverage, resulting in a dataset with 127 ingroup taxa. Maximum likelihood, coalescent, and Bayesian analyses of multiple UCE matrices recover largely congruent topologies with substantially improved support relative to previous studies. Our results strongly support the monophyly of Megachilini, the early divergence of Noteriades and Gronoceras, and a single origin of brood parasitism. All remaining non-parasitic Megachilini form a moderately supported clade sister to the brood parasitic lineage. The leafcutter bees are monophyletic and nested within dauber lineages. Several major dauber clades are consistently recovered, including an exclusively Australian clade corresponding to the Hackeriapis group of subgenera, while several recognized subgenera are paraphyletic. The lineage known as Morphella, previously placed in synonymy with the subgenus Callomegachile, was not closely related to that subgenus and is here treated as a valid subgenus. Divergence-time analyses place the crown age of Megachilini in the late Eocene to early Oligocene, with major extant lineages diversifying during the Miocene. Limited morphological diagnosability of several clades indicates that splitting non-parasitic lineages into numerous genera would result in an impractical classification that would widen the gap between taxonomists and non-specialists and exacerbate the taxonomic impediment in bees. We therefore advocate retaining a single genus Megachile for non-parasitic Megachilini (excluding Noteriades and Gronoceras), as the classification best supported by phylogenomic evidence and most robust to future taxon sampling.
Target capture is a common method of generating high throughput DNA sequencing data for phylogenetic reconstruction of species relationships, for which single copy genes are usually most informative. However, a pervasive problem with target capture is that putatively single copy genes may in fact be paralogs resulting from gene duplication, which are problematic for phylogenetic inference because their evolutionary history may differ from the divergence history of species. Here, we use as a case study a target enrichment dataset of 88 species of Detarioideae (Leguminosae) with a focus on the Sindora clade to examine approaches for handling paralogs, including the built-in paralog handling functions in HybPiper and CAPTUS, plus subsequent steps using Putative Paralog Detection and the tree-based Yang & Smith orthology inference approach. We compare the paralogs flagged using these methods and verify their performance with BLAST mapping against a reference genome sequence of Sindora glabra, and then subsequently compare the species tree topologies produced across these methods. Our comparisons of paralogs flagged across the Sindora clade show that the Putative Paralog Detection pipeline was the most accurate in identifying paralogs in terms of its similarity to the BLAST mapping, followed by the built-in paralog identification function of CAPTUS. However, the results we recovered for the Detarioideae subfamily suggest that the largest differences in species tree topology resulted from the use of paralog-filtered alignments (such as with the Putative Paralog Detection pipeline and the Yang & Smith orthology inference approaches) rather than just by removing the sequences of identified paralogous genes. This was the true for HybPiper-assembled datasets but was not seen in CAPTUS-assembled datasets. In all comparisons, the topological differences caused by different paralog handling methods tended to be confined to clades where processes such as hybridisation and introgression are prevalent. Our study provides a roadmap to establish the best approach to identify, eliminate or separate paralogs in the absence of a chromosomally contiguous reference genome for a study group, and highlights the importance of careful data inspection and processing in addition to understanding the extent of paralogy and paralog characteristics (e.g. sequence divergence between copies) for their study group.
G protein-coupled receptors (GPCRs) are one of the largest families of membrane proteins in insects, regulating vision, neural signal transduction, and various physiological behaviors. Megalopta genalis exhibits a unique facultatively eusocial lifestyle and possesses adaptations for nocturnal activity; however, its GPCR family has not yet been systematically characterized. In this study, we performed genome-wide identification, phylogenetic analysis, and expression profiling of GPCRs in M. genalis by integrating genomic annotation and transcriptomic analysis. The results showed that a total of 99 GPCRs were identified in the genome of M. genalis, which were classified into four major families. Here, we show that M. genalis has undergone lineage-specific GPCR repertoire remodeling, marked by the expansion of novel orphan receptors and the systematic loss of multiple receptor subtypes, such as the neuropeptide receptors MIP-R and NPFR. Moreover, opsins have formed a diverse array of combinations and non-GPCR odorant receptors have undergone significant expansion via tandem duplication. Together, these features may represent part of the molecular repertoire associated with the adaptation of M. genalis to a nocturnal lifestyle. Furthermore, transcriptomic analysis revealed distinct spatiotemporal expression divergence within each of the Mth/Mthl and Fz GPCR families, suggesting functional specialization across development and adult tissues. This study provides the first systematic identification and initial functional characterization of GPCRs in M. genalis, revealing an evolutionary pattern characterized by the coexistence of contraction and expansion within the GPCR family. These findings lay a foundation for further studies aimed at elucidating the roles of these GPCRs in regulating M. genalis physiology and behavior.
The somatotropic axis encompassing growth hormone (GH), insulin-like growth factor (IGF), myostatin (MSTN), and prolactin (PRL) signalling cascades is the master regulator of somatic growth, metabolism, and development in vertebrates. African catfish (Clarias gariepinus), a commercially pivotal aquaculture species, now possesses a chromosome-level reference genome (CGAR_prim_01v2); however, a systematic, genome-wide characterization spanning all five interconnected growth-related gene families has not previously been undertaken in this species. Here, we identified and characterized 15 growth-related genes spanning gh1, ghra, ghrb, Igf1, Igf2a, Igf2b, igf1ra, Igf1rb, Igf2r, Mstna, Mstnb, prl, prlra, prlrb, and smtlb distributed across 13 chromosomes. Complete one-to-one orthology with zebrafish confirmed strong dosage-balance conservation across >120 million years of teleost divergence. Physicochemical analysis resolved a clear biochemical dichotomy between compact, basic secreted ligands (19.88-45.81 kDa; pI up to 10.02) and large, acidic, heavily glycosylated membrane receptors (56.82-270.80 kDa; pI 4.85-5.97). Phylogenetic analysis confirmed 3R whole-genome duplication origins for all paralog pairs, while synteny analysis revealed a disruption of the ancestral gh1-prl chromosomal block in C. gariepinus, a finding that warrants further comparative and functional investigation. This genomic atlas provides the sequence and structural information including exon-intron boundaries, domain architecture, and chromosomal coordinates needed as a prerequisite for future marker-assisted selection and CRISPR-based myostatin-editing efforts in African catfish aquaculture, though translation into applied breeding outcomes will require subsequent functional and expression studies.
Heat shock proteins (Hsps) are essential molecular chaperones for protein homeostasis and stress responses. However, the Hsp repertoires and functions in Blattodea remain underexplored. Our genome-scale survey of nine Blattodea species revealed 37-46 conserved Hsp90, Hsp70, and DNAJ (Hsp40) genes, with DNAJ the most abundant and Hsp90 the least. Phylogenetic analysis confirmed the evolutionary conservation of three Hsp90, seven Hsp70, and 29 DNAJ subclades in Blattodea. Selection pressure analysis revealed predominant purifying selection (dN/dS ≪ 1) across lineages, strongest in DNAJ and highest in Hsp90 conservation. In Blattella germanica, expression of six representative BgHsp genes progressively increased during development, peaking in fifth-instar nymphs. Tissue expression profiling revealed that BgHspA1-2/3/4 were predominantly expressed in legs, BgDNAJB5 and BgHsp90AB1-2 were enriched in the fat body, and BgHsp90AB1 was highly expressed in the head. dsRNA injection targeting conserved Hsp gene regions achieved 61.9-94.1% knockdown of all six target genes. RNAi knockdown of six BgHsp genes disrupted wing morphogenesis, causing distinct phenotypes: wing whitening (56.7%, dsBgHspA1-4), unequal length (66.7%, dsBgHspA1-3; 76.7%, dsBgDNAJB5), and wing wrinkling (70%, dsBgHspA1-2; 63.3%, dsBgHsp90AB1; 76.7%, dsBgHsp90AB1-2). During ovarian formation, the developmental delay was most severe in the dsBgHsp90AB1 group, moderate in the dsBgHsp90AB1-2 and dsBgHspA1-2/3/4 groups, and weakest in the dsBgDNAJB5 group. Besides, knockdown significantly downregulated key developmental genes (apterous-a, nubbin, scalloped, ultrabithorax, wingless, and vitellogenin). These findings provide a reference for understanding the evolutionary patterns of Hsps in Blattodea, and offer mechanistic insights into the developmental regulation mediated by Hsps in this important public-health pest.
Acetyl-CoA carboxylase (ACC) and fatty acid synthase (FASN) are rate-limiting enzymes in the fatty acid biosynthetic pathway, yet their evolutionary relationships, sequence features, and expression profiles remain poorly understood in crustaceans, particularly in the economically important Chinese mitten crab (Eriocheir sinensis). Here, we identified and systematically analyzed ACC and FASN genes in E. sinensis using comparative genomic analyses across 43 species. ACC was highly conserved as a single-copy gene in invertebrates, in contrast to the multiple paralogs observed in vertebrates. Similarly, FASN was generally maintained as a single-copy gene across most taxa but exhibited lineage-specific expansion in certain insect groups. Phylogenetic and structural analysis revealed strong conservation of both genes within crustaceans, supported by multiple conserved motifs and canonical functional domains. Expression profiling showed predominant expression in the hepatopancreas and midgut, suggesting their potential involvement in crustacean lipid metabolism. During the molting cycle, ACC and FASN exhibited higher expression levels during stages C and D, suggesting an increased capacity for fatty acid biosynthesis before molting. In addition, dietary lipid levels experiment revealed that ACC and FASN expression responded dynamically to dietary lipid availability, with increased expression at moderate lipid levels but reduced expression under excessive lipid supplementation, indicating a possible adaptive transcriptional response to lipid status. Collectively, this study provides insights into the evolutionary conservation and expression dynamics of ACC and FASN and improves our understanding of lipid metabolic adaptation in crustaceans.
Species delimitation in the South American genus Polylepis is notoriously challenging due to high morphological similarity and phenotypic plasticity, likely driven by hybridization and gene flow. Previous phylogenetic studies suggested that genetic structure aligns more strongly with geography than with taxonomy, questioning existing species concepts and hampering conservation efforts. We used double-digest RAD sequencing (ddRADseq) to generate genome-wide SNP data for 11 Polylepis species sampled across multiple localities in Bolivia and Ecuador. Population genetic analyses, phylogenetic inference, and network approaches were combined to assess whether genetic structure aligns more closely with taxonomy or geography. Morphologically defined species formed largely cohesive genetic lineages across regions, with species identity explaining substantially more genetic variation than locality. While localized admixture and reticulation were detected among closely related taxa, widespread species showed strong genetic cohesion and clear separation from congeners. Our results indicate that the sampled Polylepis species from Bolivia and Ecuador maintain distinct genetic identities despite localized signals consistent with gene flow. This genome-wide support for current taxonomy highlights Polylepis as a valuable model for studying speciation under gene flow and indicates that multiple geographic sampling will be essential in reconstructing a robust phylogeny of the genus, with important implications for conservation planning in Andean montane forests.