Snow carps as an emerging model for polyploid genome evolution.
Explore the source record for details and available documents.
Editorial↗
SEARCH · Search PubMed
Results for “genome evolution”
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.
Cytonuclear conflict and reticulate evolution in the Morelloid clade (Solanum, Solanaceae): Insights from genome skimming and network Phylogenomics.
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.
Phylogeny↗
Correction for Grünig et al., Evolution after whole-genome duplication (WGD) drives phenotypic and transcriptomic divergence more than WGD in an autopolyploid herb.
Explore the source record for details and available documents.
Journal Article↗
Evidence of genome-wide relaxed selection on mildly deleterious mutations in an ancient subterranean catfish.
About one hundred subterranean catfish species have been described, resulting from repeated colonization of cave environments by multiple surface lineages. Most cave-dwelling species are found in the Americas, in particular in South America, but a few species also live in Central and North America. Despite the availability of high-quality genome assemblies for two cave species, the Mexican blind catfish Prietella phreatophila and the Colombian blind catfish Trichomycterus rosablanca, genomic approaches to investigate genetic changes associated with subterranean life or to estimate cave colonization times remain largely unexplored. To fill this gap, we additionally sequenced the genomes of four blind and depigmented subterranean catfishes from Peru (three Trichomycterus and one Astroblepus), as well as the genomes of four close surface relatives. We first extracted a large set of light-related genes, such as phototransduction and crystallin genes, and found contrasting decays of these sequences in different cave species, from 1% of pseudogenes in T. rosablanca to 48% in P. phreatophila. Two independent molecular dating methods gave congruent ages, indicating that these catfishes colonized subterranean habitats at different times, ranging from Early Pliocene to Late Pleistocene, supporting the hypothesis that surface catfishes repeatedly and rapidly adapted to subterranean habitats. The oldest cave species, P. phreatophila, appears to have been thriving in the dark for over 3.5 million years. Moreover, a genome-wide analysis of protein-coding genes suggests weaker purifying selection on mildly deleterious mutations in this cavefish than in other catfish lineages, likely reflecting a long-term small effective population size.
cavefishes↗
Divergent trajectories of genome architecture and chromosome evolution in ferns and angiosperms.
Ferns and angiosperms represent the two largest vascular plant lineages but exhibit striking genomic and ecological contrasts. We investigated whether differences in genome size, chromosome architecture, GC content, and stomatal traits reveal divergent evolutionary trajectories between these lineages. We assembled the most comprehensive dataset to date, integrating genome size, chromosome number and size, GC content, and stomatal traits for over 1100 fern species and compared it with an extensive angiosperm dataset. Ferns exhibited markedly lower variability and c. 16-fold slower rates of chromosome size evolution than angiosperms. A persistent positive relationship between genome size and chromosome number in ferns suggests limited cytological post-polyploid diploidization. While ferns generally possess larger stomata, this difference disappears after accounting for genome size, indicating that nucleotypic constraints, rather than lineage-specific physiology, dictate stomatal dimensions. Both groups share a unimodal GC-genome size relationship peaking at c. 14 Gbp. Larger fern chromosomes imply lower genome-wide recombination rates, potentially limiting genetic reshuffling and adaptive potential. Our results highlight fundamentally divergent evolutionary trajectories, likely shaped by meiotic symmetry in ferns and meiotic asymmetry, possibly centromere drive, and post-polyploid diploidization in angiosperms, defining the functional and genomic landscapes of these lineages across deep evolutionary timescales.
Genome, Plant↗
Comparative genomic and proteomic analysis reveals orthogroup structured evolution of tick protease inhibitors.
Protease inhibitors (PIs) play central roles in regulating endogenous proteolysis and host-parasite interactions in ticks. However, the evolutionary architecture underlying their diversification across tick lineages remains insufficiently resolved. Here, we performed a genome-wide comparative analysis of predicted proteomes from 14 tick species to systematically characterize PI repertoires. In total, 4931 putative PIs were identified and grouped into 20 families using the MEROPS classification system. Further, PI families such as Antistasin, WAP-type, and Pacifastin, which have not previously been systematically reported in tick genomes, were classified. Orthogroup inference demonstrated that PI expansion is structured at the level of evolutionary lineages rather than uniformly across families. By stratifying orthogroups according to duplication burden and taxonomic conservation, we identified a broadly conserved single-copy core under strong purifying selection. Motif level analysis of serpin reactive center loops further revealed conservation of inhibitory specificity within single copy orthogroups and diversification of key functional residues in duplication-associated lineages. Integration of secretion prediction and tissue-resolved proteomics from Hyalomma anatolicum and Rhipicephalus microplus demonstrated that evolutionary stratification is reflected at the protein level. Together, these findings provide an orthogroup-resolved evolutionary framework linking duplication dynamics, molecular evolution, and tissue-level protein deployment. This integrative approach offers a systematic basis for prioritizing conserved and diversified PI lineages for future functional and anti-tick intervention studies.
Animals↗
Comparative analyses of olfactory receptor repertoires in Schizothorax fish based on the chromosome-level genomes: Implications for regulatory roles of dietary differentiation and ploidy variation.
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.
Animals↗
Fructophilic lactic acid bacteria as a window into multi-scale convergent evolution.
Fructophilic lactic acid bacteria (FLAB) are a group of lactic acid bacteria with unique growth characteristics, that is, poor growth on glucose. Their growth is enhanced in the presence of fructose or external electron acceptors. These organisms inhabit fructose-rich environments such as flowers, fruits, and pollinating insects, particularly honey bees. Apilactobacillus spp. and Fructobacillus spp. are representatives of FLAB, although they belong to phylogenetically distant clades. These organisms commonly possess markedly small genomes with a low number of coding DNA sequences. Furthermore, their genomes are characterized by a markedly reduced number of genes involved in carbohydrate transport and metabolism. Genome reduction in FLAB reflects convergent adaptation to fructose-rich environments rather than general genome streamlining. The two distinct FLAB genera, Fructobacillus and Apilactobacillus, independently lost more than 100 genes in statistically similar orders. In contrast, genes involved in carbohydrate and amino acid metabolism exhibited reversed orders of loss between the two genera. Furthermore, FLAB genomes lack an intact bifunctional alcohol/aldehyde dehydrogenase gene (adhE), which causes their poor growth on glucose. A comparative genomic study suggested the evolutionary process underlying adhE gene decay during adaptation to the fructose-rich environments, including pollinating insects. In conclusion, FLAB represent a unique example of habitat-driven convergent reductive evolution that can be investigated across multiple biological scales - from individual genes to whole genomes - in the diverse LAB group with a wide range of habitats, and partially share the fructophilic evolution with eukaryotic yeasts found in fructose-rich habitats.
Fructose↗
Evolution of maize recombination landscape during domestication.
Despite the plethora of knowledge about the benefits of meiotic recombination and numerous theoretical studies examining how recombination rates evolve, there is a general lack of empirical support and consensus across species. To fill this knowledge gap, we characterized the evolution of recombination landscape in maize during its domestication from teosinte and related the observed changes to established theoretical frameworks. Through examining recombination in experimental populations of maize and teosinte and the population genomics approach of identifying historical recombination events using ancestral recombination graph inference to generate saturated maize and teosinte recombination maps, we found that during domestication, maize experienced a 12% increase in its genome-wide recombination rate. Furthermore, maize evolved higher recombination rates on the long arms of chromosomes in regions closer to centromeres, where recombination is generally very low. The repatterning of crossover events came from changes in global crossover positioning rather than alterations in cis-acting chromatin factors. Consequently, we found evidence of selection acting on trans-acting recombination modifiers affecting crossover interference and controlling the interference-dependent class I crossover pathway. We show that CO repatterning was likely beneficial for maize fitness, as significant recombination rate increases were predominantly in gene-rich regions, which harbor domestication-related variation. This work suggests genomic and mechanistic processes leading to the evolution of meiotic recombination landscape in response to directional selection pressure and provides evidence for the evolutionary advantage of recombination.
Zea mays↗
Predictive evolutionary genomics: principles, validation, and practice.
Climate change and habitat loss are driving rapid evolutionary responses in populations world-wide, which creates an urgent need for evolutionary forecasting in conservation and agriculture. Such forecasting can be categorized into three time scales: trait-based models that use multivariate quantitative genetic equations to project correlated phenotypic responses up to c. 20 generations, allele-based analyses that model allele frequency dynamics up to 100 generations, and composite adaptation scores that aggregate many small effects to yield predictions across longer horizons. However, these approaches have remained largely disconnected. Here, we present a Bayesian framework that integrates these three complementary approaches for evolutionary prediction. Our framework combines genomic, phenotypic, and environmental data to yield probabilistic predictions with explicit uncertainty. We show how predictive evolutionary forecasts can be validated with experimental evolution, field experimentation, historical specimens, and reciprocal transplants. These validated forecasts can help advance conservation and agricultural programmes by helping predict which populations are at risk of future extinction, optimizing breeding programmes for future climates, and planning ecosystem management under environmental change. By supporting a shift towards more predictive approaches in evolutionary biology, this framework may help improve our ability to manage biodiversity and food security in a changing world.
Genomics↗
Experimental evolution reveals contrasting adaptive landscapes in lab and field environments.
Experimental evolution is widely used to infer microbial responses to environmental change, yet most laboratory studies impose constant, well-mixed conditions that differ fundamentally from fluctuating, spatially structured field environments. We compared genomic evolution in the leaf litter-associated bacterium Curtobacterium strain MMLR14_002 under control and warming treatments in laboratory culture and in a complementary field experiment. Laboratory-derived isolates accumulated more mutations per genome and exhibited stronger locus-level parallelism, with mutations recurring in a small number of coding loci. Field-derived isolates accumulated fewer mutations per genome, and these mutations rarely occurred in the same coding loci across replicate populations. Instead, field isolates exhibited a higher proportion of intergenic mutations, with mutations recurring in the same intergenic regions across independent field deployments. When coding mutations were detected in the field, they were distributed across functionally diffuse targets and more often involved metabolic pathways than the core cellular processes repeatedly targeted during laboratory evolution. Warming itself did not consistently influence mutation accumulation or the genomic distribution of mutations; instead, laboratory and field contexts primarily shaped the accumulation, targets, and repeatability of genomic change. These results suggest that laboratory thermal evolution identifies adaptive routes favored under sustained selection but may overestimate coding-level parallelism under heterogeneous field conditions. Bridging laboratory and field evolution will likely require experimental designs that incorporate temporal variability and spatial heterogeneity characteristic of natural systems.IMPORTANCEA central goal of experimental evolution is to infer how microbes evolve in nature from laboratory studies. Here, we evaluate this assumption by comparing genomic evolution of a leaf litter-associated Curtobacterium strain in laboratory and field warming experiments to identify broad patterns rather than isolate the contribution of any single environmental factor. We find that the strong parallelism at coding loci observed under laboratory conditions is reduced in the field, while mutations recurring in the same intergenic regions across field deployments suggest that parallel evolution in nature may more often involve regulatory noncoding regions rather than coding targets. These results show that environmental context reshapes adaptive landscapes and may limit the parallelism of coding-level genomic responses inferred from homogeneous laboratory conditions.
experimental evolution↗
Integrating genomic distance analyses in the description of a new family, genus, and species of sponge-associated antipatharians (black corals).
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.
Animals↗
Introgression shapes the genomic conflict landscape of Malus, providing evidence for a reticulate backbone in a woody crop lineage.
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.
Malus↗
Evolutionary architecture and lineage-specific diversification of Forkhead box transcription factors in Perna viridis.
The Forkhead box (Fox) transcription factors are evolutionarily conserved regulators of development, cell cycle, and apoptosis across metazoans. This study provides the first comprehensive genome-wide analysis of the Fox gene family in the Asian green mussel (Perna viridis). We identified 28 Fox genes distributed across 10 chromosomes. Comparative analysis reveals the absence of the FoxI, FoxQ1, FoxR and FoxS subfamily, consistent with other bivalves and indicative of lineage-specific gene loss during molluscan evolution. Notably, gene duplications in the FoxAB, FoxD, FoxH, FoxN1-4, FoxQ2 and FoxQD subfamilies may reflect functional diversification associated with environmental adaptation. Exon-intron structural variability, including intron loss in several paralogues, suggests structural diversification and potential regulatory variation. Phylogenetic reconstruction confirmed the monophyly of core Fox classes while highlighting divergent expansion patterns in lophotrochozoans. Selection analyses showed strong purifying selection across duplicated Fox paralogs, supporting functional conservation after lineage-specific expansion. Gene Ontology enrichment linked Fox genes to stress response, apoptosis, and transcriptional regulation. By integrating phylogenetic, structural, and transcriptomic analyses, this study provides a genomic framework for understanding Fox gene organisation, evolution, and tissue-associated expression patterns in Perna viridis and establishes a comparative resource for future functional studies in bivalves.
Animals↗
A system-level metastable model of cancer evolution: integrating replication stress, cell cycle deregulation and chromosomal instability.
INTRODUCTION: Cancer cell proliferation occurs within the context of persistent genomic instability. In this review, we propose the RS-CCD-CIN axis as a systems-level framework in which replication stress (RS), cell cycle deregulation (CCD) and chromosomal instability (CIN) form an interdependent triad that shapes tumour evolution. This axis represents a constrained metastable state in which genomic instability is tolerated and buffered. The objective of this review is to synthesize the current understanding of how the RS-CCD-CIN axis contributes to tumour heterogeneity, adaptability and therapy response. DISCUSSION: Evidence indicates that RS, CCD and CIN operate as a dynamic, interconnected network rather than as independent processes. Replication stress induces DNA damage and mutagenesis, while partial checkpoint disruption permits cells with unresolved lesions to proliferate. Chromosomal instability generates both structural and numerical alterations, contributing to intratumoural heterogeneity. Together, these processes facilitate adaptation to environmental and therapeutic pressures. Extrachromosomal DNA, micronuclei formation and cytosolic DNA signalling, including the cGAS-STING pathway, connect genomic instability to adaptive responses and immune modulation. Single-cell and spatial profiling reveal temporal and spatial variability in RS, CCD and CIN states, highlighting the limitations of static biomarkers. Therapeutically, targeting individual components often yields limited durability, whereas approaches that simultaneously perturb multiple aspects of the RS-CCD-CIN axis may improve clinical outcomes. CONCLUSIONS: This review highlights the RS-CCD-CIN axis as a fragile and metastable architecture that supports cancer evolution, while also being susceptible to collapse. A deeper understanding of this interconnected framework may inform the development of therapeutic strategies and enhance the management of resistance.
Humans↗
Genome-wide insights into the evolutionary and demographic history of the red alga Mazzaella laminarioides: Evidence for speciation with ancient migration along the southeast Pacific coast.
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.
Rhodophyta↗
Comprehensive identification and evolutionary analysis of the Wnt gene family in bivalves: Insights into the larval development of the noble scallop Chlamys nobilis.
The Wnt gene family regulates fundamental developmental processes in metazoans, but its evolutionary composition and developmental deployment in bivalves remain largely unresolved. Here, we performed a comparative genomic analysis of Wnt genes in 19 bivalve species and examined developmental expression profiles in the noble scallop Chlamys nobilis, with Crassostrea gigas and Chlamys farreri used for cross-species comparison. A total of 235 Wnt genes were identified and assigned to 12 subfamilies. No reliable Wnt3 ortholog was detected in any analyzed bivalve, supporting the view that Wnt3 loss occurred early during lophotrochozoan evolution rather than representing a lineage-specific absence. Most Wnt proteins retained the conserved WNT domain, indicating strong structural conservation, whereas lineage-specific copy-number variation and gene loss were observed among species. C. farreri and C. gigas each retained 12 Wnt genes and lacked Wnt3, whereas C. nobilis lacked Wnt3, Wnt7, and Wnt16. Developmental transcriptome analysis and RT-qPCR revealed clear stage-specific expression patterns. In C. gigas, Wnt2/10/A were highly expressed during earlydevelopment and peaked around the D-shaped larval stage, while Wnt8 and Wnt11 showed distinct stage-specific peaks. By contrast, Wnt1/5/6/9 were more active during later larval development or juvenile formation. These results provide a comparative framework for bivalve Wnt evolution and identify candidate Wnt genes potentially involved in larval development and aquaculture-relevant developmental transitions.
Animals↗
Beyond species trees: pervasive gene flow limits phylogenomic resolution in the diversification of Juniperus from the Qinghai-Tibet Plateau.
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.
Gene Flow↗