Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “plant 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.

At least 19 recordsLinked to original sources
PubMed2026

Nucleotide Combination Proportions Across Algae, Monocotyledons and Dicotyledons: Insights into Plant Genome Evolution.

Plant evolution started with unicellular algae, gradually evolving multicellularity and terrestrial colonization. These evolutionary events were accompanied by the interplay of chromosome polyploidization, rearrangement, gene loss, and point mutation. We counted the proportion of nucleotide combinations in the genome sequences of 64 sequenced plants, and analyzed the significant difference in these nucleotide combination proportions among algae, monocotyledons and dicotyledons. The correlation of highly significant different and no significant different nucleotide combinations was analyzed respectively. Nucleotide combinations and their reverse complementary sequence proportions were analyzed in different functional regions of the genome. These results reveal that some nucleotide combinations are subject to strict selection, and these combinations have a higher proportion in the CDS regions and lower proportion in the intergenic regions. Meanwhile, there are some nucleotide combinations that are under less selective pressure, and these combinations have a higher proportion in the intergenic regions and lower proportion in the CDS regions. Cluster analysis based on trinucleotide to octanucleotide combination proportions reveals that plant genome evolution is accompanied by clade-wide differentiation of genome-wide nucleotide composition patterns, in addition to well-documented chromosomal polyploidization, structural rearrangement and gene loss events. We analyzed the changes in the proportion of nucleotide combinations at the genome level in 64 sequenced plants, providing a new idea for studying genome evolution in the plant kingdom.

comparative genomics↗
PubMed2025

Genomic analyses of three Acanthus L. species provide insight into polyploidization-driven speciation and evolution.

Allopolyploidy fundamentally influences plant evolution, yet the genomic dynamics of allotetraploidization remain incompletely understood. We investigated Acanthus tetraploideus (2n = 4x = 96), an ecologically significant allotetraploid true mangrove from Indo-West Pacific intertidal zones. Our prior integrative investigations indicate that A. tetraploideus originated through hybridization of the diploid species A. ilicifolius and A. ebracteatus with subsequent chromosome doubling. Here, we present complete chromosome-scale genome assemblies for all three species, representing the first genomic resources for true mangrove polyploid research. Our analysis reveals that the three species have experienced at least four rounds of polyploidization events, with the most recent, approximately 53 mya, possibly an Acanthus-specific event. The allotetraploid A. tetraploideus, which emerged between 1.5 and 2.2 mya, has A. ebracteatus as its maternal progenitor and A. ilicifolius as its paternal one. Through a comprehensive genomic comparison and analysis of homoeologous gene expression, we propose a gradual evolutionary trajectory for allotetraploidy in A. tetraploideus. Despite the allotetraploidization event dating back to around 2 mya, A. tetraploideus retains a high degree of colinearity with its ancestral genomes, with the majority (76.2%) of duplicated genes retained and no significant sub-genome bias in gene expression. Furthermore, we have identified positive selection in specific genes that may facilitate the adaptation of Acanthus mangrove species to their intertidal habitats. These findings establish A. tetraploideus as a model for studying allopolyploid evolution while providing new insights into mangrove speciation processes.

Genome, Plant↗
PubMed2026

Recent gene duplication and structural remodeling drive rapid lineage-specific gene family evolution in plants.

Gene duplication promotes the generation of novel gene functions and trait diversity across species. Here, we present DupHIST, a computational pipeline that reconstructs the hierarchical timing of gene duplications by integrating maximum likelihood (ML)-based phylogeny with substitution-derived timing via statistical smoothing. Applied to over 4.5 million genes from 114 plant genomes, we successfully inferred duplication histories across nearly 130,000 orthogroups. This large-scale analysis showed that 53.0% of genes arose from recent, lineage-specific duplications, with high concentrations in particular multi-copy families. Among these, NLR, C48, and P450 families exemplified how recently duplicated genes undergo rapid stepwise structural remodeling. This process was primarily driven by small-scale mutations, including insertions, deletions, and frameshifts, that rapidly accumulated shortly after duplication. By resolving the precise duplication order, we reconstructed these architectural changes, thereby enabling both the inference of putative ancestral structures and the exploration of functional diversification arising from structural remodeling. Structure-based clustering further uncovered that recently duplicated, uncharacterized genes retain core domain structures resembling known functional proteins even across phylogenetically distant species lacking sequence homology. Our findings reveal that recent gene duplications and subsequent structural remodeling represent a widespread and lineage-specific force driving rapid diversification of gene families in plants.

Gene duplication history↗
PubMed1976

The evolution of plants and animals under domestication: the contribution of studies at the molecular level.

Protein molecules are essential catalysts in life processes and also form much of the substance of living material. Their three dimensional structures determine their biological function. Their biosynthesis is primarily determined by arrays of nucleic acid macromolecules (DNA and RNA), and the amino acid sequences that constitute their long spatially organized peptide-chain molecules reflect at one remove this DNA coding system, and thus record a step-by-step history of some of the viable genetic events (natural or man-controlled) that have created the organism and the breed. Amino acid sequences can be used to trace the progress of controlled breeding in two ways: by extrapolation back from living breeds, and by analysis of ancient protein material. Of the latter, bone or tendon or skin collagens and hair keratins are the most perfectly preserved as molecular structures through 20,000 years and indeed much longer. Amino acid sequences are expensive to determine (collagen has 1052 amino acid residues), and the potential of this palaeobiological information has been as yet little exploited. The first approach has, however, been more explored, in both plants and animals. Several protein systems must be studied in conjunction to reveal the phylogenetic threads in any one breed. As the three dimensional quaternary structure of protein molecules becomes more appreciated in relation to biological function, and as new techniques and procedures are developed, amino acid sequence data can become more informative in our ultimate understanding of early selective breeding.

Amino Acid Sequence↗
PubMed2025

Evolutionary dynamics and genetic diversity of transposable elements revealed by resequencing data in maize population.

Zea mays (maize) is a globally significant crop with a complex genome enriched with transposable elements (TEs), which are crucial drivers of genomic diversity and plant evolution. In this study, we identified the TE insertion loci (TILs) from resequencing data of 103 maize accessions with the developed pipeline, and 64 293 non-redundant unique TILs were obtained in 82 maize accessions after filtering; approximately 80% (51 361) of loci showed insertion polymorphisms within the population. All TE superfamilies have low frequency in the maize population except for short interspersed nuclear elements, while some TE families have high fixed TE insertions, revealing distinct evolutionary dynamics among TE superfamilies and families. Genetic analysis using the transposon insertion polymorphism information from the maize population showed that the TE polymorphism loci can reflect their geographical origin and evolutionary relationships. Furthermore, TE insertions could also significantly impact gene expression, implying functional consequences for maize phenotypes and adaptation. These findings provide valuable insights into the evolutionary dynamics and genetic diversity of maize genomes, offering a valuable resource for molecular markers and association studies.

Zea mays↗
PubMed2025

Subgenomic divergence and functional innovation following whole-genome duplication in Maleae species of Rosaceae.

Whole-genome duplication (WGD) drives plant evolution by inducing karyotype rearrangements and gene loss through subgenome fractionation. In this study, we investigate post-WGD evolutionary dynamics in Rosaceae, focusing on Maleae species, which uniquely experienced an additional WGD. Using phylogenetic and synteny analyses, we reveal that chromosomal breakpoints act as hotspots for localized fractionation, contributing to blurred homoeologous origins and influencing gene retention patterns. Here, we reconstruct karyotype evolution across Rosaceae subfamilies, highlighting chromosome reductions and lineage-specific rearrangements in Dryadoideae, Rosoideae, and Amygdaloideae. We also identify a bias for retaining transcription factors and hormone-related genes from older WGDs in subsequent polyploidy events. Transcriptome analysis classifies WGD-derived genes in Maleae species, such as apple and loquat, into three expression groups, with hormone-enriched genes playing roles in lignification and fruit-related innovations. These findings demonstrate the interplay between chromosomal breakpoints, biased retention, and functional divergence, revealing their contributions to genomic and phenotypic evolution in Maleae and their adaptive success within Rosaceae.

Genome, Plant↗
PubMed2026

Profile of C. Robin Buell.

C. Robin Buell has been a leading figure in plant genomics since the advent of DNA sequencing technology. She helped lead multiple consortia to sequence some of the first crop genomes at the turn of the millennium. She has since used the genomes to tackle questions in fundamental biology, plant evolution, and agriculture. Recently, she applied single-cell technologies to uncover how complex biosynthetic pathways are compartmentalized across different rare plant cell types. Now at the University of Georgia, Buell explores in her Inaugural Article how tubers arose repeatedly in the plant family tree.

Genome, Plant↗
PubMed2026

Generational variation and stabilization in resynthesized allotetraploid Brassica juncea derived from diploid progenitors B. rapa and B. nigra.

BACKGROUND: Polyploidy is a major driver of plant evolution and crop improvement, generating novel variation in morphology, physiology, and agronomic traits. Brassica juncea (AABB, 2n = 36), a natural allotetraploid derived from B. rapa (AA) and B. nigra (BB), is an important oilseed and vegetable crop; however, its narrow genetic base limits further breeding gains. Resynthesized B. juncea (RBJ), developed from known progenitors, provides a tractable system to investigate polyploid stabilization, trait diversification, and generational variation. This study evaluated RBJ across nine generations (F1-S8) to elucidate generational variation in morphological, molecular, cytological, and oil content traits during progressive stabilization. RESULTS: Substantial variation was observed for key yield-related traits, including siliqua length, seeds per siliqua, and thousand-seed weight. High estimates of heritability, genotypic variance, and genetic advance indicated their potential utility in selection based improvement. Comparative analyses revealed a clear generational progression, characterized by relatively enhanced performance in early generations, increased recombination-driven variability in intermediate generations, and the partial stabilization of several traits in later generations. Generation mean analysis suggested the involvement of additive, dominance, and epistatic gene effects in trait inheritance. Molecular analysis using SSR markers confirmed the amphidiploid origin and genomic integrity of RBJ generations. Cytological assessments, pollen viability assays, and flow cytometric analysis collectively demonstrated stable chromosome numbers, improved fertility, and maintenance of ploidy stability across successive generations. CONCLUSIONS: The study provides valuable insights into the generational variation and stabilization of morphological, molecular, and oil content traits in resynthesized B. juncea. The findings suggest that variability arising from polyploidization and interspecific hybridization undergoes gradual reorganization across successive generations, leading to increased trait stabilization and more consistent expression of selected agronomic characteristics. Collectively, these results contribute to the understanding of early stabilization processes in RBJ, highlighting resynthesized polyploids as useful systems for studying variation and stabilization in allopolyploid crops.

Mustard Plant↗
PubMed2025

MEANtools integrates multi-omics data to identify metabolites and predict biosynthetic pathways.

During evolution, plants have developed the ability to produce a vast array of specialized metabolites, which play crucial roles in helping plants adapt to different environmental niches. However, their biosynthetic pathways remain largely elusive. In the past decades, increasing numbers of plant biosynthetic pathways have been elucidated based on approaches utilizing genomics, transcriptomics, and metabolomics. These efforts, however, are limited by the fact that they typically adopt a target-based approach, requiring prior knowledge. Here, we present MEANtools, a systematic and unsupervised computational integrative omics workflow to predict candidate metabolic pathways de novo by leveraging knowledge of general reaction rules and metabolic structures stored in public databases. In our approach, possible connections between metabolites and transcripts that show correlated abundance across samples are identified using reaction rules linked to the transcript-encoded enzyme families. MEANtools thus assesses whether these reactions can connect transcript-correlated mass features within a candidate metabolic pathway. We validate MEANtools using a paired transcriptomic-metabolomic dataset recently generated to reconstruct the falcarindiol biosynthetic pathway in tomato. MEANtools correctly anticipated five out of seven steps of the characterized pathway and also identified other candidate pathways involved in specialized metabolism, which demonstrates its potential for hypothesis generation. Altogether, MEANtools represents a significant advancement to integrate multi-omics data for the elucidation of biochemical pathways in plants and beyond.

Metabolomics↗
PubMed2025

Comprehensive identification and analysis of clusters of tandemly duplicated genes reveal their contributions to adaptive evolution of green plants.

Tandem gene duplication occurred more frequently compared with the episodic whole-genome duplication (WGD), providing a continuous supply of genetic material for evolutionary innovation and adaptation to changing environments. The rising roles of clusters of tandemly duplicated genes (CTDGs) in the evolution of phenotypic diversity have been unraveled in mammals. However, the content and biological roles of CTDGs remain largely unknown in plants. Here, we comprehensively identified CTDGs in 220 published plant genomes representing major lineages of green plants. The number of CTDGs showed great variation across taxa, ranging from 0 to 6028. The size of CTDGs varied from 2 to 47 genes, with small clusters containing two members predominating. Interestingly, significant expansion of CTDGs was found in early-diverging land plants and is closely associated with the evolution of key traits (e.g., ABA response, plant cuticle, UV-B resistance) required for plants to conquer terrestrial environments. Functional enrichment analysis revealed conserved and specialized functional profiles among different sizes of CTDGs in both Arabidopsis thaliana and the bryophyte Physcomitrium patens. Small CTDGs were enriched in fundamental stress responses, including protein modification, signal transduction, and responses to diverse stress stimuli, while large CTDGs were enriched in more sophisticated processes such as plant hormone biosynthesis and signaling, plant-microbe interactions, and reproductive processes. Expression pattern analyses of CTDGs under different stress conditions in A. thaliana and P. patens revealed that the highest number of CTDGs showed differential expression under drought stress, suggesting important roles of CTDGs in the evolution of desiccation tolerance in early land plants. The results of this study provide new additions to our knowledge about the abundance of CTDGs across green plants and reveal their important contributions to enable plants to overcome stressful environments on land.

Gene Duplication↗
PubMed2025

NGATHA carpel development genes evolved in the common ancestor of seed plants.

The evolution of the carpel, the defining feature of angiosperms, remains a fundamental question in plant biology. Understanding how this organ originated is crucial because it underpins the reproductive success and diversity of flowering plants. Here, we investigated the functional conservation between gymnosperms and angiosperms of key transcription factors involved in carpel development. We found that Ginkgo biloba homologs can functionally substitute for their angiosperm counterparts in stigma development. We discovered that GbRAV5 is related to angiosperm NGA genes, challenging previous notions that these are exclusive to angiosperms, and we found a parallel loss of the AP2 domain in gymnosperms providing a rare snapshot of how protein families evolve. Conserved protein interactions and overlapping expression patterns of GbRAV5 and GbHEC in Ginkgo ovules suggest that the molecular toolkit for carpel development was largely present in the last common ancestor of seed plants, offering new insights into the evolution of reproductive structures.

Plant Proteins↗
PubMed2024

A gene with a thousand alleles: The hyper-variable effectors of plant-parasitic nematodes.

Pathogens are engaged in a fierce evolutionary arms race with their host. The genes at the forefront of the engagement between kingdoms are often part of diverse and highly mutable gene families. Even in this context, we discovered unprecedented variation in the hyper-variable (HYP) effectors of plant-parasitic nematodes. HYP effectors are single-gene loci that potentially harbor thousands of alleles. Alleles vary in the organization, as well as the number, of motifs within a central hyper-variable domain (HVD). We dramatically expand the HYP repertoire of two plant-parasitic nematodes and define distinct species-specific "rules" underlying the apparently flawless genetic rearrangements. Finally, by analyzing the HYPs in 68 individual nematodes, we unexpectedly found that despite the huge number of alleles, most individuals are germline homozygous. These data support a mechanism of programmed genetic variation, termed HVD editing, where alterations are locus specific, strictly governed by rules, and theoretically produce thousands of variants without errors.

Animals↗
PubMed2026

The evolution of Ca2+-ATPases across plants with profiles in Rhododendron and the function of key members in alleviating high calcium stress.

Ca2+-ATPase (CAP) is a key Ca2+ efflux protein in plants. Our previous research suggests that CAPs may play a crucial role in the adaptation of rhododendrons to high calcium environments. However, the evolution, variation, characteristic expression, and subfunctionalization of this gene family in Rhododendron remain unknown. Through the analysis of pan-genomes and pan-transcriptomes, we elucidated the systematic evolution of CAPs in plants, as well as their characteristic expression patterns in Rhododendron. During the evolutionary process from lower to higher plants, CAPs can be divided into six clades and exhibit structural conservation. CAPs have emerged and differentiated in lower plants such as algae, and they have undergone significant amplification in Eudicots plants like rhododendrons. Three Rhododendron species (Rhododendron bailiense, R. delavayi, and R. irroratum) located in the karst province of Guizhou in Southwest China exhibit the highest copies of CAPs, suggesting a strong association between CAP copy number variation and habitat, particularly in high calcium environments. Through multiple transcriptome analyses, we revealed that CAPs are induced under various environmental/developmental conditions (e.g. karst environments, high altitude, early flower development, hormones, etc.). Co-expression network analysis highlighted key members of calcineurin B-like protein (CBL) and CBL-interacting protein kinases (CIPK) that are associated with the high expression of CAPs. Experimental validation demonstrated that CAPb1 and CAPd1 significantly alleviate high calcium stress, and the CAPb1-CIPK1-CBL1 and CAPd1-CIPK2-CBL1 modules can further enhance the alleviation. These findings provide new insights into the evolution, characteristic expression, and function of CAPs, as well as new perspectives on the high calcium adaptability of rhododendrons.

Journal Article↗
PubMed2026

Comparison of paralog identification methods and their impact on species tree topologies in target capture phylogenomics within the Sindora clade (Detarioideae: Leguminosae).

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.

Phylogeny↗
PubMed2026

Haplotype-resolved genome of Forsythia suspensa reveals the reticulate evolution in Oleaceae and a novel gene cluster regulating stamen development.

The olive family (Oleaceae) comprises numerous species of economic, horticultural, and medicinal importance. Despite its significance, the evolutionary history of this complex family remains enigmatic. Here, we generated a high-quality haplotype-resolved genome of Forsythia suspensa, a distylous species that occupies a key phylogenetic position in Oleaceae. The 2 haplotypes exhibit significant allelic divergence with potential allele-specific regulation. We reconstructed the polyploidization history of Oleaceae by confirming and precisely dating a shared whole-genome triplication and an independent whole-genome duplication event. We revealed a complex reticulate evolution that gave rise to the tribe Oleeae: an initial hybridization between Forsythieae (♂) and Jasmineae (♀), a subsequent backcrossing event, and a final whole-genome duplication. We identified a novel tandemly duplicated pectin methylesterase inhibitor gene cluster that regulates filament length and pollen size via restricting cell elongation in the long-styled morph. Dosage augmentation via stepwise cluster formation (0.99 to 3.83 Mya) may contribute to maintaining stamen traits of the long-styled morph. These FsPMEIs are co-expressed with many cell wall-related genes, suggesting a functional link in cell wall modification. Our study reveals the reticulate evolution in Oleaceae and a novel gene cluster controlling stamen development in F. suspensa and provides valuable haplotype-resolved genomic resources for heterostylous species, offering novel framework and molecular pathways to understand plant adaptive evolution.

Forsythia↗
PubMed2026

The genetic control of rapid genome content divergence in Arabidopsis thaliana.

Genome evolution in eukaryotes is predominantly driven by the dynamics of repetitive sequences, which vary widely in both copy number and sequence composition. Rates of repeat evolution differ between and within species and are likely modulated by both genetics and environment. To uncover factors shaping the rate of genome content evolution, we analyzed 1043 resequenced Arabidopsis thaliana genomes using a novel K-mer-based approach to characterize genome content variation and identify hypervariable regions underlying differences in repeat abundance. We next treated repeat abundance as a quantitative trait and performed genome-wide association analyses across more than 400 repeat families to identify the genetic basis of copy number variation. Integrating these results through a meta-GWAS approach revealed both cis-acting variants and more than 50 candidate trans-acting loci associated with repeat abundance genome-wide. Cis-acting variation was predominantly localized to pericentromeric and centromeric regions, whereas trans-acting loci were enriched for candidate genes involved in DNA replication, DNA repair, and DNA methylation regulation. The results are consistent with purifying selection acting against mutations that accelerate genome content divergence, favoring alleles that constrain repeat expansion. Together, these findings provide new insights into the genetic architecture and evolutionary forces shaping genome evolution in A. thaliana and establish a framework for investigating these processes in other plant species.

Arabidopsis↗
PubMed2026

Evolution of virulence of a plant RNA virus in developmental stage-structured host populations.

Natural host populations are age-structured, and developmental stages differ in susceptibility and within-host pathogen dynamics, potentially imposing distinct selective pressures on viruses. However, the evolutionary consequences of host age structure remain poorly understood. We experimentally evolved turnip mosaic potyvirus for 5 passages in Arabidopsis thaliana populations spanning 7 demographic regimes, from juvenile- to mature-dominated cohorts. We quantified disease progression, symptom severity, and viral load, cross-inoculated evolved lineages across host stages to construct infection matrices, and performed whole-population sequencing at passages 1 and 5. Disease traits changed markedly with passage, demography, and their interaction. Disease progression evolved faster in older populations, whereas symptom severity was independent of median age, indicating demographic reweighting of virulence components. Viral load increased across passages and positively correlated with severity, linking within-host fitness to symptoms. Cross-inoculation assays revealed a modular infection network: juvenile-evolved lineages specialized on juvenile hosts, whereas lineages from intermediate and older populations were more generalist. Genomically, we detected both parallel and demography-specific adaptations, including recurrent changes in the viral protein VPg (involved in translation, replication, and host interactions) as well as synonymous variants showing consistent or opposing selection across host population stage structures. Overall, host age structure emerges as a major ecological driver of virulence evolution, shaping tradeoffs between disease progression and severity and determining specialization versus generalism. These results integrate phenotypic and genomic responses and suggest that manipulating crop age structure could steer virus evolution toward less damaging outcomes.

Virulence↗
PubMed2026

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↗