An integrative genomic and phenomic analysis to investigate the nature of plant species in Escallonia (Escalloniaceae)

When the sand blossoms: Phylogeny, trait evolution, and geography of speciation in Linanthus

PREMISE

Understanding how plants successfully diversified in novel environments is a central question in evolutionary biology. Linanthus occurs in arid areas of western North America and exhibits extensive floral trait variation, multiple color polymorphisms, differences in blooming time, and variation in life history strategies. We reconstructed the evolutionary history of this genus.

METHODS

We generated restriction-site associated (ddRAD) sequences for 180 individuals and target capture (TC) sequences for 63 individuals, with complete species sampling. Using maximum likelihood and pseudo-coalescent approaches, we inferred phylogenies of Linanthus and used them to model the evolution of phenotypic traits and investigate the genus's geographic speciation history.

RESULTS

Relationships are consistent and well supported with both ddRAD and TC data. Most species are monophyletic despite extensive local sympatry and range overlap, suggesting strong isolating barriers. The non-monophyly of the night-blooming and perennial species may be due to rapid speciation or cryptic diversity. Perenniality likely evolved from annuality, a rare shift in angiosperms. Night-blooming evolved three times independently. Flower color polymorphism is an evolutionarily labile trait that is likely ancestral. No single geographic mode of speciation characterizes this diversification, but most species overlap in range, which suggests that they evolved in parapatry.

CONCLUSIONS

Our results illustrate the complexity of phylogenetic inference for recent radiations, even with multiple sources of genomic data and extensive sampling. This analysis provides a foundation for understanding aridity adaptations, such as evolution of flower color polymorphisms, night-blooming, and perenniality, as well as speciation mechanisms.

Tracing evolutionary trajectories in the presence of gene flow in South American temperate lizards (Squamata: Liolaemus kingii group)

Evolutionary processes behind lineage divergence often involve multidimensional differentiation. However, in the context of recent divergences, the signals exhibited by each dimension may not converge. In such scenarios, incomplete lineage sorting, gene flow, and scarce phenotypic differentiation are pervasive. Here, we integrated genomic (RAD loci of 90 individuals), phenotypic (linear and geometric traits of 823 and 411 individuals, respectively), spatial, and climatic data to reconstruct the evolutionary history of a speciation continuum of liolaemid lizards (Liolaemus kingii group). Specifically, we (a) inferred the population structure of the group and contrasted it with the phenotypic variability; (b) assessed the role of postdivergence gene flow in shaping phylogeographic and phenotypic patterns; and (c) explored ecogeographic drivers of diversification across time and space. We inferred eight genomic clusters exhibiting leaky genetic borders coincident with geographic transitions. We also found evidence of postdivergence gene flow resulting in transgressive phenotypic evolution in one species. Predicted ancestral niches unveiled suitable areas in southern and eastern Patagonia during glacial and interglacial periods. Our study underscores integrating different data and model-based approaches to determine the underlying causes of diversification, a challenge faced in the study of recently diverged groups. We also highlight Liolaemus as a model system for phylogeographic and broader evolutionary studies.

Draft genome assemblies for two species of Escallonia (Escalloniales)

OBJECTIVES

Escallonia (Escalloniaceae) belongs to the Escalloniales, a diverse clade of flowering plants with unclear placement in the tree of life. Escallonia species show impressive morphological and ecological diversity and are widely distributed across three hotspots of biodiversity in the Neotropics. To shed light on the genomic substrate of this radiation and the phylogenetic placement of Escalloniales as well as to generate useful data for comparative evolutionary genomics across flowering plants, we produced and annotated draft genomes for two species of Escallonia.

DATA DESCRIPTION

Genomic DNA from E. rubra and E. herrerae was sequenced with Oxford Nanopore sequencing chemistry, generating 3.4 and 12 million sequence reads with an average read length of 9.4 and 9.1 Kb (approximately 31 and 111 Gb of sequence data), respectively. In addition, we generated Illumina 100-bp paired-end short read data for E. rubra (approximately 75 Gb of sequence data). The Escallonia rubra genome was 566 Mb, with 3,233 contigs and an N50 of 285 Kb. The assembled genome for E. herrerae was 994 Mp, with 5,760 contigs and an N50 of 317 Kb. The genome sequences were annotated with 31,038 (E. rubra) and 47,905 (E. herrerea) protein-coding gene models supported by transcriptome/protein evidence and/or Pfam domain content. BUSCO assessments indicated completeness levels of approximately 98% for the genome assemblies and 88% for the genome annotations.