Pleistocene glacial cycles drive isolation, gene flow and speciation in the high-elevation Andes

Incomplete lineage sorting and hybridization underlie tree discordance in Petunia and related genera (Petunieae, Solanaceae)

Despite the overarching history of species divergence, phylogenetic studies often reveal distinct topologies across regions of the genome. The sources of these gene tree discordances are variable, but incomplete lineage sorting (ILS) and hybridization are among those with the most biological importance. Petunia serves as a classic system for studying hybridization in the wild. While field studies suggest that hybridization is frequent, the extent of reticulation within Petunia and its closely related genera has never been examined from a phylogenetic perspective. In this study, we used transcriptomic data from 11 Petunia, 16 Calibrachoa, and 10 Fabiana species to illuminate the relationships between these species and investigate whether hybridization played a significant role in the diversification of the clade. We inferred that gene tree discordance within genera is linked to hybridization events along with high levels of ILS due to their rapid diversification. Moreover, network analyses estimated deeper hybridization events between Petunia and Calibrachoa, genera that have different chromosome numbers. Although these genera cannot hybridize at the present time, ancestral hybridization could have played a role in their parallel radiations, as they share the same habitat and life history.

Botany and geogenomics: Constraining geological hypotheses in the neotropics with large-scale genetic data derived from plants

Decades of empirical research have revealed how the geological history of our planet shaped plant evolution by establishing well-known patterns (e.g., how mountain uplift resulted in high rates of diversification and replicate radiations in montane plant taxa). This follows a traditional approach where botanical data are interpreted in light of geological events. In this synthesis, I instead describe how by integrating natural history, phylogenetics, and population genetics, botanical research can be applied alongside geology and paleontology to inform our understanding of past geological and climatic processes. This conceptual shift aligns with the goals of the emerging field of geogenomics. In the neotropics, plant geogenomics is a powerful tool for the reciprocal exploration of two long standing questions in biology and geology: how the dynamic landscape of the region came to be and how it shaped the evolution of the richest flora. Current challenges that are specific to analytical approaches for plant geogenomics are discussed. I describe the scale at which various geological questions can be addressed from biological data and what makes some groups of plants excellent model systems for geogenomics research. Although plant geogenomics is discussed with reference to the neotropics, the recommendations given here for approaches to plant geogenomics can and should be expanded to exploring long-standing questions on how the earth evolved with the use of plant DNA.

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.

Pleistocene diversification of unifoliolate-leaved Lupinus (Leguminosae: Papilionoideae) in Florida

The importance and prevalence of recent ice-age and post-glacial speciation and species diversification during the Pleistocene across many organismal groups and physiographic settings are well established. However, the extent to which Pleistocene diversification can be attributed to climatic oscillations and their effects on distribution ranges and population structure remains debatable. In this study, we use morphologic, geographic and genetic (RADseq) data to document Pleistocene speciation and intra-specific diversification of the unifoliolate-leaved clade of Florida Lupinus, a small group of species largely restricted to inland and coastal sand ridges across the Florida peninsula and panhandle. Phylogenetic and demographic analyses alongside morphological and geographic evidence suggest that recent speciation and intra-specific divergence within this clade were driven by a combination of non-adaptive allopatric divergence caused by edaphic niche conservatism and opportunities presented by the emergence of new post-glacial sand ridge habitats. These results highlight the central importance of even modest geographic isolation and short periods of allopatric divergence following range expansion in the emergence of new taxa and add to the growing evidence that Pleistocene climatic oscillations may contribute to rapid diversification in a myriad of physiographic settings. Furthermore, our results shed new light on long-standing taxonomic debate surrounding the number of species in the Florida unifoliate Lupinus clade providing support for recognition of five species and a set of intra-specific variants. The important conservation implications for the narrowly restricted, highly endangered species Lupinus aridorum, which we show to be genetically distinct from its sister species Lupinus westianus, are discussed.

Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex

As genomic data proliferates, the prevalence of post-speciation gene flow is making species boundaries and relationships increasingly ambiguous. Although current approaches inferring fully bifurcating phylogenies based on concatenated datasets provide simple and robust answers to many species relationships, they may be inaccurate because the models ignore inter-specific gene flow and incomplete lineage sorting. To examine the potential error resulting from ignoring gene flow, we generated both a RAD-seq and a 500 protein-coding loci highly multiplexed amplicon (HiMAP) dataset for a monophyletic group of 12 species defined as the Bactrocera dorsalis sensu lato clade. With some of the world's worst agricultural pests, the taxonomy of the B. dorsalis s.l. clade is important for trade and quarantines. However, taxonomic confusion confounds resolution due to intra- and interspecific phenotypic variation and convergence, mitochondrial introgression across half of the species, and viable hybrids. We compared the topological convergence of our datasets using concatenated phylogenetic and various multispecies coalescent approaches, some of which account for gene flow. All analyses agreed on species delimitation, but there was incongruence between species relationships. Under concatenation, both datasets suggest identical species relationships with mostly high statistical support. However, multispecies coalescent and multispecies network approaches suggest markedly different hypotheses and detected significant gene flow. We suggest that the network approaches are likely more accurate because gene flow violates the assumptions of the concatenated phylogenetic analyses, but the data-reductive requirements of network approaches resulted in reduced statistical support and could not unambiguously resolve gene flow directions. Our study highlights the importance of testing for gene flow, particularly with phylogenomic datasets, even when concatenated approaches receive high statistical support.

Strategies of diaspore dispersal investment in Compositae: the case of the Andean highlands

BACKGROUND AND AIMS

Understanding diaspore morphology and how much a species invests on dispersal appendages is key for improving our knowledge of dispersal in fragmented habitats. We investigate diaspore morphological traits in high-Andean Compositae and their main abiotic and biotic drivers and test whether they play a role in species distribution patterns across the naturally fragmented high-Andean grasslands.

METHODS

We collected diaspore trait data for 125 Compositae species across 47 tropical high-Andean summits, focusing on achene length and pappus-to-achene length ratio, with the latter as a proxy of dispersal investment. We analysed the role of abiotic (temperature, elevation and latitude) and biotic factors (phylogenetic signal and differences between tribes) on diaspore traits and whether they are related to distribution patterns across the Andes, using phylogenomics, distribution modelling and community ecology analyses.

KEY RESULTS

Seventy-five percent of the studied species show small achenes (length <3.3 mm) and 67% have high dispersal investment (pappus length at least two times the achene length). Dispersal investment increases with elevation, possibly to compensate for lower air density, and achene length increases towards the equator, where non-seasonal climate prevails. Diaspore traits show significant phylogenetic signal, and higher dispersal investment is observed in Gnaphalieae, Astereae and Senecioneae, which together represent 72% of our species. High-Andean-restricted species found across the tropical Andes have, on average, the pappus four times longer than the achene, a significantly higher dispersal investment than species present only in the northern Andes or only in the central Andes.

CONCLUSIONS

Small achenes and high diaspore dispersal investment dominate among high-Andean Compositae, traits typical of mostly three tribes of African origin; but traits are also correlated with the environmental gradients within the high-Andean grasslands. Our results also suggest that diaspore dispersal investment is likely to shape species distribution patterns in naturally fragmented habitats.

Generic concepts and species diversity within the Gynoxyoid clade (Senecioneae, Compositae)

The Gynoxyoid clade of the Senecioneae (Asteraceae) until now included the five genera Aequatorium, Gynoxys, Nordenstamia, Paracalia and Paragynoxys as diagnosed using selected morphological characters. In their pre-phylogenetic circumscription, the genera Aequatorium and Paragynoxys were considered to inhabit the northern Andes in contrast to Nordenstamia and Paracalia that occur in the central Andes. The most species-rich genus, Gynoxys, was believed to be distributed throughout the Andes. We use a recently established plastid phylogenomic framework that rendered Gynoxys paraphyletic to further evaluate the delimitation of genera in the Gynoxyoid clade. We examine the morphological variation of all members of the Gynoxyoid to identify characters potentially informative at genus level. This results in a matrix of eleven, mostly multistate characters, including those originally used to diagnose these genera. The ancestral character state inference displays a high level of homoplasy, but nevertheless supports the recognition of four genera. Aequatorium is characterised by white radiate capitula. Paracalia and Paragynoxys share white flowers and floral characteristics, such as flower opening and length of disc flowers lobes, as plesiomorphic states, but differ in habit (scandent shrubs vs. trees). Paracalia also retained white flowers, but its two species are characterised by the absence of outer phyllaries. The genera Gynoxys and Nordenstamia comprise species with yellow capitula which appear to be a derived feature in the Gynoxyoids. The genus Nordenstamia, with eight species, is synonymised under Gynoxys since molecular evidence shows its species nested within various parts of the Gynoxys subclade and the morphological variation of Nordenstamia falls well within that of Gynoxys. With the goal to assign all species to four genera (Aequatorium, Gynoxys, Paracalia and Paragynoxys), we assess the states for the eleven characters for all members of the Gynoxyoids and generate new ETS and ITS sequences for 171 specimens belonging to 49 species to further support their generic placement. We provide a taxonomic treatment for the four genera recognised here including amended diagnoses and morphological descriptions. Furthermore, a species-level taxonomic backbone is elaborated for all genera using electronic tools that list 158 currently accepted names and synonyms (209 names in total) with the respective protologue and type information, as well as notes on the current understanding of species limits. Eleven names are newly synonymised, two are lectotypified and eight are newly transferred to other genera.

Are 150 km of open sea enough? Gene flow and population differentiation in a bat-pollinated columnar cactus

Genetic differentiations and phylogeographical patterns are controlled by the interplay between spatial isolation and gene flow. To assess the extent of gene flow across an oceanic barrier, we explored the effect of the separation of the peninsula of Baja California on the evolution of mainland and peninsular populations of the long-lived columnar cactus Stenocereus thurberi. We analyzed twelve populations throughout the OPC distribution range to assess genetic diversity and structure using chloroplast DNA sequences. Genetic diversity was higher (Hd = 0.81), and genetic structure was lower (GST = 0.143) in mainland populations vs peninsular populations (Hd = 0.71, GST = 0.358 respectively). Genetic diversity was negatively associated with elevation but positively with rainfall. Two mainland and one peninsular ancestral haplotypes were reconstructed. Peninsular populations were as isolated among them as with mainland populations. Peninsular haplotypes formed a group with one mainland coastal population, and populations across the gulf shared common haplotypes giving support to regular gene flow across the Gulf. Gene flow is likely mediated by bats, the main pollinators and seed dispersers. Niche modeling suggests that during the Last Glacial Maximum (c. 130 ka), OPC populations shrank to southern locations. Currently, Stenocereus thurberi populations are expanding, and the species is under population divergence despite ongoing gene flow. Ancestral populations are located on the mainland and although vicariant peninsular populations cannot be ruled out, they are likely the result of gene flow across the seemingly formidable barrier of the Gulf of California. Still, unique haplotypes occur in the peninsula and the mainland, and peninsular populations are more structured than those on the mainland.

Allopatric speciation is more prevalent than parapatric ecological divergence in a recent high-Andean diversification (Linochilus: Asteraceae)

Elucidating how species accumulate in diversity hotspots is an ongoing debate in evolutionary biology. The páramo, in the Northern Andes, has remarkably high indices of plant diversity, endemicity, and diversification rates. A hypothesis for explaining such indices is that allopatric speciation is high in the páramo given its island-like distribution. An alternative hypothesis is that the altitudinal gradient of the Andean topography provides a variety of niches that drive vertical parapatric ecological speciation. A formal test for evaluating the relative roles of allopatric and parapatric ecological speciation is lacking. The main aim of our study is to test which kind of speciation is more common in an endemic páramo genus. We developed a framework incorporating phylogenetics, species' distributions, and a morpho-ecological trait (leaf area) to compare sister species and infer whether allopatric or parapatric ecological divergence caused their speciation. We applied our framework to the species-rich genus Linochilus (63 spp.) and found that the majority of recent speciation events in it (12 events, 80%) have been driven by allopatric speciation, while a smaller fraction (one event, 6.7%) is attributed to parapatric ecological speciation; two pairs of sister species produced inconclusive results (13.3%). We conclude that páramo autochthonous (in-situ) diversification has been primarily driven by allopatric speciation.

Phylogenomics and Biogeography of the Mammilloid Clade Revealed an Intricate Evolutionary History Arose in the Mexican Plateau

Mexico harbors ~45% of world’s cacti species richness. Their biogeography and phylogenomics were integrated to elucidate the evolutionary history of the genera Coryphantha, Escobaria, Mammillaria, Mammilloydia, Neolloydia, Ortegocactus, and Pelecyphora (Mammilloid Clade). We analyzed 52 orthologous loci from 142 complete genomes of chloroplast (103 taxa) to generate a cladogram and a chronogram; in the latter, the ancestral distribution was reconstructed with the Dispersal-Extinction-Cladogenesis model. The ancestor of these genera arose ~7 Mya on the Mexican Plateau, from which nine evolutionary lineages evolved. This region was the site of 52% of all the biogeographical processes. The lineages 2, 3 and 6 were responsible for the colonization of the arid southern territories. In the last 4 Mya, the Baja California Peninsula has been a region of prolific evolution, particularly for lineages 8 and 9. Dispersal was the most frequent process and vicariance had relevance in the isolation of cacti distributed in the south of Mexico. The 70 taxa sampled as Mammillaria were distributed in six distinct lineages; one of these presumably corresponded to this genus, which likely had its center of origin in the southern part of the Mexican Plateau. We recommend detailed studies to further determine the taxonomic circumscription of the seven genera.

The importance of the Andes in the evolutionary radiation of Sigmodontinae (Rodentia, Cricetidae), the most diverse group of mammals in the Neotropics

The Andean mountains stand out for their striking species richness and endemicity that characterize many emblematic Neotropical clades distributed in or around these mountains. The radiation of the Sigmodontinae subfamily, the most diversified mammalian group in the Neotropics, has been historically related to Andean orogenesis. We aim to evaluate this interplay between geological processes and biological responses through the diversification dynamics, the biogeographical history, and the range evolution of the subfamily. For these, we built the most comprehensive phylogeny and gathered 14,836 occurrences for the subfamily. We identified one shift in the speciation rate in the genus Akodon, which suffered their Andean radiation after the arrival of non-Andean ancestors. Our biogeographic analyses show multiple dispersal paths throughout the evolution that allowed this subfamily to colonize all Neotropics. The Northern Andes and Central-Southern Andes were the most important sources of diversity. In addition, the Central-Southern Andes were the most relevant sink, receiving the highest number of lineages. The Andean region exhibited higher speciation and turnover rates than non-Andean regions. Thus, our results support the crucial role of the Andean Mountains in the Sigmodontinae radiation, acting as a "macroevolutionary cradle" and "species attractor" for several sigmodontine lineages at different times, and as a "species pump" becoming the biogeographic source of multiple widely distributed neotropical lineages. Then, complex macroevolutionary dynamics would explain these rodents' high extant Andean diversity and their wide distribution in the Neotropics.

Complex cycles of divergence and migration shape lineage structure in the common kingsnake species complex

Aim

The Nearctic is a complex patchwork of habitats and geologic features that form barriers to gene flow resulting in phylogeographic structure and speciation in many lineages. Habitats are rarely stable over geologic time, and the Nearctic has undergone major climatic changes in the past few million years. We use the common kingsnake species complex to study how climate, geography, and history influence lineage formation over a large, complex landscape.

Location

Nearctic/North America.

Taxon

Common kingsnake, Lampropeltis getula, species complex.

Methods

We analyzed genome-wide sequence data from 51 snakes spanning the majority of the species complex's range. We used population clustering, generalized dissimilarity modeling, and coalescent methods to identify the number of genetic clusters within the L. getula complex, infer the environmental correlates of genetic differentiation, and estimate models of divergence and gene flow among lineages.

Results

We identified three major lineages within the L. getula complex and further continuous spatial structure within lineages. The most important ecological correlates of genetic distance in the complex are related to aridity and precipitation, consistent with lineage breaks at the Great Plains/Desert ecotone and the Cochise Filter Barrier. Lineages are estimated to have undergone multiple rounds of isolation and secondary contact, with highly asymmetric migration occurring at present.

Main conclusions

Changing climates combined with a large and geologically complex landscape have resulted in a mosaic of discrete and spatially continuous genetic structure. Multiple rounds of isolation and secondary contact as climate fluctuated over the past ~4.4 My have likely driven the evolution of discrete lineages that maintain high levels of gene flow. Continuous structure is strongly shaped by aridity and precipitation, suggesting roles for major precipitation gradients in helping to maintain lineage identity in the face of gene flow when lineages are in geographic contact.

The chromosome-scale genome of Kobresia myosuroides sheds light on karyotype evolution and recent diversification of a dominant herb group on the Qinghai-Tibet Plateau

Kobresia species are common in meadows on the Qinghai-Tibet Plateau. They are important food resources for local livestock, and serve a critical foundation for ecosystem integration. Genetic resources of Kobresia species are scarce. Here, we generated a chromosome-level genome assembly for K. myosuroides (Cyperaceae), using PacBio long-reads, Illumina short-reads, and Hi-C technology. The final assembly had a total size of 399.9 Mb with a contig N50 value of 11.9 Mb. The Hi-C result supported a 29 pseudomolecules model which was in consistent with cytological results. A total of 185.5 Mb (44.89% of the genome) transposable elements were detected, and 26,748 protein-coding genes were predicted. Comparative analysis revealed that Kobresia plants have experienced recent diversification events during the late Miocene to Pliocene. Karyotypes analysis indicated that the fission and fusion of chromosomes have been a major driver of speciation, which complied with the lack of whole-genome duplication (WGD) in K. myosuroides genome. Generally, this high-quality reference genome provides insights into the evolution of alpine sedges, and may be helpful to endemic forage improvement and alpine ecosystem preservation.

Pervasive hybridization during evolutionary radiation of Rhododendron subgenus Hymenanthes in mountains of southwest China

Radiations are especially important for generating species biodiversity in mountainous ecosystems. The contribution of hybridization to such radiations has rarely been examined. Here, we use extensive genomic data to test whether hybridization was involved in evolutionary radiation within Rhododendron subgenus Hymenanthes, whose members show strong geographic isolation in the mountains of southwest China. We sequenced genomes for 143 species of this subgenus and 93 species of four other subgenera, and found that Hymenanthes was monophyletic and radiated during the late Oligocene to middle Miocene. Widespread hybridization events were inferred within and between the identified clades and subclades. This suggests that hybridization occurred both early and late during diversification of subgenus Hymenanthes, although the extent to which hybridization, speciation through mixing-isolation-mixing or hybrid speciation, accelerated the diversification needs further exploration. Cycles of isolation and contact in such and other montane ecosystems may have together promoted species radiation through hybridization between diverging populations and species. Similar radiation processes may apply to other montane floras in this region and elsewhere.

Considerations on the Systematics of Neotropical Pierina, with the Description of Two New Species of Phulia Herrich-Schäffer from the Peruvian Andes (Lepidoptera: Pieridae, Pierinae, Pierini)

A comparative analysis of high-Andean Pierina was carried out, including a total of 25 species. Based on morphological evidence, with an emphasis on venation and genitalia and molecular data, using three genetic markers, we confirm the recent subjective synonymy of the generic names Tatochila Butler, 1870, Piercolias, Staudinger, 1894, Hypsochila Ureta, 1955, Infraphulia Field, 1958, Pierphulia Field, 1958, and Theochila Field, 1958 with Phulia Herrich-Schäffer, 1867. Two new species are described, namely Phulia stoddardi Pyrcz & Cerdeña n. sp., from the Andes of Central Peru, which occurs at an unusually high altitude of close to 5000 m a.s.l. in dry puna habitat, and Phulia phantasma Lamas, Willmott & Boyer n. sp., from dry montane forests in northern Peru and southern Ecuador. An overview of high-elevation butterflies is presented, with some discussion on adaptations to this environment.

Hyb-Seq provides new insights into the phylogeny and evolution of the Chrysanthemum zawadskii species complex in China

A species complex is an assemblage of closely related species with blurred boundaries, and from which species could arise from different speciation processes and/or a speciation continuum. Such a complex can provide an opportunity to investigate evolutionary mechanisms acting on speciation. The Chrysanthemum zawadskii species complex in China, a monophyletic group of Chrysanthemum, consists of seven species with considerable morphological variation, diverse habitats and different distribution patterns. Here, we used Hyb-Seq data to construct a well-resolved phylogeny of the C. zawadskii complex. Then, we performed comparative analyses of variation patterns in morphology, ecology and distribution to investigate the roles of geography and ecology in this complex's diversification. Lastly, we implemented divergence time estimation, species distribution modelling and ancestral area reconstruction to trace the evolutionary history of this complex. We concluded that the C. zawadskii complex originated in the Qinling-Daba mountains during the early Pliocene and then spread west and northward along the mountain ranges to northern China. During this process, geographical and ecological factors imposing different influences resulted in the current diversification and distribution patterns of this species complex, which is composed of both well-diverged species and diverging lineages on the path of speciation.

Linked selection, differential introgression and recombination rate variation promote heterogeneous divergence in a pair of yellow croakers

Understanding the mechanisms underlying heterogeneous genomic divergence is of particular interest in evolutionary biology. Highly differentiated genomic regions, known as genomic islands, often evolve between diverging lineages. These genomic islands may be related to selection promoting adaptation or reproductive isolation. Based on whole genome assembly and genome-wide RAD sequencing in a pair of yellow croakers (genus: Larimichthys), we investigated the evolutionary processes shaping genomic landscapes of divergence. Demographic modelling indicated that the two species diverged following a secondary contact scenario, where differential introgression and linked selection were suggested to be involved in heterogeneous genomic divergence. We identified reduced recombination rate in genomic islands and a relatively good conservation of both genetic diversity and recombination landscapes between species, which highlight the roles of linked selection and recombination rate variation in promoting heterogeneous divergence in the common ancestral lineage of the two species. In addition, we found a positive correlation between differentiation (F_ST ) and absolute sequence divergence (D_xy ), and elevated D_xy in genomic islands, indicating that the genomic landscape of divergence was not shaped by linked selection alone. Restricted gene flow in highly differentiated regions has probably remodelled the landscape of heterogeneous genomic divergence. This study highlights that highly differentiated genomic regions can also arise from a combination of linked selection and differential gene flow in interaction with varying recombination rates.

Genetic markers in Andean Puya species (Bromeliaceae) with implications on plastome evolution and phylogeny

The Andean plant endemic Puya is a striking example of recent and rapid diversification from central Chile to the northern Andes, tracking mountain uplift. This study generated 12 complete plastomes representing nine Puya species and compared them to five published plastomes for their features, genomic evolution, and phylogeny. The total size of the Puya plastomes ranged from 159,542 to 159,839 bp with 37.3%–37.4% GC content. The Puya plastomes were highly conserved in organization and structure with a typical quadripartite genome structure. Each of the 17 consensus plastomes harbored 133 genes, including 87 protein‐coding genes, 38 tRNA (transfer RNA) genes, and eight rRNA (ribosomal RNA) genes; we found 69–78 tandem repeats, 45–60 SSRs (simple sequence repeats), and 8–22 repeat structures among 13 species. Four protein‐coding genes were identified under positive site‐specific selection in Puya. The complete plastomes and hypervariable regions collectively provided pronounced species discrimination in Puya and a practical tool for future phylogenetic studies. The reconstructed phylogeny and estimated divergence time for the lineage suggest that the diversification of Puya is related to Andean orogeny and Pleistocene climatic oscillations. This study provides plastome resources for species delimitation and novel phylogenetic and biogeographic studies.

The role of multiple Pleistocene refugia in promoting diversification in the Pacific Northwest

Pleistocene glacial cycles drastically changed the distributions of taxa endemic to temperate rainforests in the Pacific Northwest, with many experiencing reduced habitat suitability during glacial periods. In this study, we investigate whether glacial cycles promoted intraspecific divergence and whether subsequent range changes led to secondary contact and gene flow. For seven invertebrate species endemic to the PNW, we estimated Species Distribution Models (SDMs) and projected them onto current and historical climate conditions to assess how habitat suitability changed during glacial cycles. Using single nucleotide polymorphism (SNP) data from these species, we assessed population genetic structure and used a machine-learning approach to compare models with and without gene flow between populations upon secondary contact after the Last Glacial Maximum (LGM). Finally, we estimated divergence times and rates of gene flow between populations. SDMs suggest that there was less suitable habitat in the North Cascades and Northern Rocky Mountains during glacial compared to interglacial periods, resulting in reduced habitat suitability and habitat fragmentation during the LGM. Our genomic data identify population structure in all taxa, and support gene flow upon secondary contact in five of the seven taxa. Parameter estimates suggest that population divergences date to the later Pleistocene for most populations. Our results support a role of refugial dynamics in driving intraspecific divergence in the Cascades Range. In these invertebrates, population structure often does not correspond to current biogeographic or environmental barriers. Rather, population structure may reflect refugial lineages that have since expanded their ranges, often leading to secondary contact between once isolated lineages.

Unlocking Andean sigmodontine diversity: five new species of Chilomys (Rodentia: Cricetidae) from the montane forests of Ecuador

The Andean cloud forests of Ecuador are home to several endemic mammals. Members of the Thomasomyini rodents are well represented in the Andes, with Thomasomys being the largest genus (47 species) of the subfamily Sigmodontinae. Within this tribe, however, there are genera that have escaped a taxonomic revision, and Chilomys Thomas, 1897, constitutes a paradigmatic example of these "forgotten" Andean cricetids. Described more than a century ago, current knowledge of this externally unmistakable montane rodent is very limited, and doubts persist as to whether or not it is monotypic. After several years of field efforts in Ecuador, a considerable quantity of specimens of Chilomys were collected from various localities representing both Andean chains. Based on an extensive genetic survey of the obtained material, we can demonstrate that what is currently treated as C. instans in Ecuador is a complex comprising at least five new species which are described in this paper. In addition, based on these noteworthy new evidence, we amend the generic diagnosis in detail, adding several key craniodental traits such as incisor procumbency and microdonty. These results indicate that Chilomys probably has a hidden additional diversity in large parts of the Colombian and Peruvian territories, inviting a necessary revision of the entire genus.

The Andes through time: evolution and distribution of Andean floras

The Andes are the world's most biodiverse mountain chain, encompassing a complex array of ecosystems from tropical rainforests to alpine habitats. We provide a synthesis of Andean vascular plant diversity by estimating a list of all species with publicly available records, which we integrate with a phylogenetic dataset of 14 501 Neotropical plant species in 194 clades. We find that (i) the Andean flora comprises at least 28 691 georeferenced species documented to date, (ii) Northern Andean mid-elevation cloud forests are the most species-rich Andean ecosystems, (iii) the Andes are a key source and sink of Neotropical plant diversity, and (iv) the Andes, Amazonia, and other Neotropical biomes have had a considerable amount of biotic interchange through time.

Population structure in Neotropical plants: integrating pollination biology, topography and climatic niches

Animal pollinators mediate gene flow among plant populations, but in contrast to well‐studied topographic and (Pleistocene) environmental isolating barriers, their impact on population genetic differentiation remains largely unexplored. Comparing how these multifarious factors drive microevolutionary histories is, however, crucial for better resolving macroevolutionary patterns of plant diversification. Here we combined genomic analyses with landscape genetics and niche modelling across six related Neotropical plant species (424 individuals across 33 localities) differing in pollination strategy to test the hypothesis that highly mobile (vertebrate) pollinators more effectively link isolated localities than less mobile (bee) pollinators. We found consistently higher genetic differentiation (F_ST) among localities of bee‐ than vertebrate‐pollinated species with increasing geographical distance, topographic barriers and historical climatic instability. High admixture among montane populations further suggested relative climatic stability of Neotropical montane forests during the Pleistocene. Overall, our results indicate that pollinators may differentially impact the potential for allopatric speciation, thereby critically influencing diversification histories at macroevolutionary scales.

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

What we mean by species and whether they have any biological reality has been debated since the early days of evolutionary biology. Some biologists even suggest that plant species are created by taxonomists as a subjective, artificial division of nature. However, the nature of plant species has been rarely tested critically with data while ignoring taxonomy. We integrate phenomic and genomic data collected across hundreds of individuals at a continental scale to investigate this question in Escallonia (Escalloniaceae), a group of plants which includes 40 taxonomic species (the species proposed by taxonomists). We first show that taxonomic species may be questionable as they match poorly to patterns of phenotypic and genetic variation displayed by individuals collected in nature. We then use explicit statistical methods for species delimitation designed for phenotypic and genomic data, and show that plant species do exist in Escallonia as an objective, discrete property of nature independent of taxonomy. We show that such species correspond poorly to current taxonomic species (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< 20\%$$\end{document}<20%) and that phenomic and genomic data seldom delimit congruent entities (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< 20\%$$\end{document}<20%). These discrepancies suggest that evolutionary forces additional to gene flow can maintain the cohesion of species. We propose that phenomic and genomic data analyzed on an equal footing build a broader perspective on the nature of plant species by helping delineate different ‘types of species’. Our results caution studies which take the accuracy of taxonomic species for granted and challenge the notion of plant species without empirical evidence. Note: A version of the complete manuscript in Spanish is available in the Supplemental Materials.

The evolutionary assembly of forest communities along environmental gradients: recent diversification or sorting of pre-adapted clades?

Recent studies have demonstrated that ecological processes that shape community structure and dynamics change along environmental gradients. However, much less is known about how the emergence of the gradients themselves shape the evolution of species that underlie community assembly. In this study, we address how the creation of novel environments leads to community assembly via two nonmutually exclusive processes: immigration and ecological sorting of pre-adapted clades (ISPC), and recent adaptive diversification (RAD). We study these processes in the context of the elevational gradient created by the uplift of the Central Andes. We develop a novel approach and method based on the decomposition of species turnover into within- and among-clade components, where clades correspond to lineages that originated before mountain uplift. Effects of ISPC and RAD can be inferred from how components of turnover change with elevation. We test our approach using data from over 500 Andean forest plots. We found that species turnover between communities at different elevations is dominated by the replacement of clades that originated before the uplift of the Central Andes. Our results suggest that immigration and sorting of clades pre-adapted to montane habitats is the primary mechanism shaping tree communities across elevations.

Evolutionary history of two rare endemic conifer species from the eastern Qinghai-Tibet Plateau

BACKGROUND AND AIMS

Understanding the population genetics and evolutionary history of endangered species is urgently needed in an era of accelerated biodiversity loss. This knowledge is most important for regions with high endemism that are ecologically vulnerable, such as the Qinghai-Tibet Plateau (QTP).

METHODS

The genetic variation of 84 juniper trees from six populations of Juniperus microsperma and one population of Juniperus erectopatens, two narrow-endemic junipers from the QTP that are sister to each other, was surveyed using RNA-sequencing data. Coalescent-based analyses were used to test speciation, migration and demographic scenarios. Furthermore, positively selected and climate-associated genes were identified, and the genetic load was assessed for both species.

KEY RESULTS

Analyses of 149 052 single nucleotide polymorphisms showed that the two species are well differentiated and monophyletic. They diverged around the late Pliocene, but interspecific gene flow continued until the Last Glacial Maximum. Demographic reconstruction by Stairway Plot detected two severe bottlenecks for J. microsperma but only one for J. erectopatens. The identified positively selected genes and climate-associated genes revealed habitat adaptation of the two species. Furthermore, although J. microsperma had a much wider geographical distribution than J. erectopatens, the former possesses lower genetic diversity and a higher genetic load than the latter.

CONCLUSIONS

This study sheds light on the evolution of two endemic juniper species from the QTP and their responses to Quaternary climate fluctuations. Our findings emphasize the importance of speciation and demographic history reconstructions in understanding the current distribution pattern and genetic diversity of threatened species in mountainous regions.

Population dynamics linked to glacial cycles in Cercis chuniana F. P. Metcalf (Fabaceae) endemic to the montane regions of subtropical China

The mountains of subtropical China are an excellent system for investigating the processes driving the geographical distribution of biodiversity and radiation of plant populations in response to Pleistocene climate fluctuations. How the major mountain ranges in subtropical China have affected the evolution of plant species in the subtropical evergreen broadleaved forest is an issue with long-term concern. Here, we focused on Cercis chuniana, a woody species endemic to the southern mountain ranges in subtropical China, to elucidate its population dynamics. We used genotyping by sequencing (GBS) to investigate the spatial pattern of genetic variation among 11 populations. Geographical isolation was detected between the populations located in adjacent mountain ranges, thought to function as geographical barriers due to their complex physiography. Bayesian time estimation revealed that population divergence occurred in the middle Pleistocene, when populations in the Nanling Mts. separated from those to the east. The orientation and physiography of the mountain ranges of subtropical China appear to have contributed to the geographical pattern of genetic variation between the eastern and western populations of C. chuniana. Complex physiography plus long-term stable ecological conditions across glacial cycles facilitated the demographic expansion in the Nanling Mts., from which contemporary migration began. The Nanling Mts. are thus considered as a suitable area for preserving population diversity and large population sizes of C. chuniana compared with other regions. As inferred by ecological niche modeling and coalescent simulations, secondary contact occurred during the warm Lushan-Tali Interglacial period, with intensified East Asia summer monsoon and continuous habitat available for occupation. Our data support the strong influence of both climatic history and topographic characteristics on the high regional phytodiversity of the subtropical evergreen broadleaved forest in subtropical China.

Spatiotemporal evolution of the global species diversity of Rhododendron

Evolutionary radiation is a widely recognized mode of species diversification, but its underlying mechanisms have not been unambiguously resolved for species-rich cosmopolitan plant genera. In particular, it remains largely unknown how biological and environmental factors have jointly driven its occurrence in specific regions. Here, we use Rhododendron, the largest genus of woody plants in the Northern Hemisphere, to investigate how geographic and climatic factors, as well as functional traits, worked together to trigger plant evolutionary radiations and shape the global patterns of species richness based on a solid species phylogeny. Using 3,437 orthologous nuclear genes, we reconstructed the first highly supported and dated backbone phylogeny of Rhododendron comprising 200 species that represent all subgenera, sections, and nearly all multispecies subsections, and found that most extant species originated by evolutionary radiations when the genus migrated southward from circumboreal areas to tropical/subtropical mountains, showing rapid increases of both net diversification rate and evolutionary rate of environmental factors in the Miocene. We also found that the geographically uneven diversification of Rhododendron led to a much higher diversity in Asia than in other continents, which was mainly driven by two environmental variables, that is, elevation range and annual precipitation, and were further strengthened by the adaptation of leaf functional traits. Our study provides a good example of integrating phylogenomic and ecological analyses in deciphering the mechanisms of plant evolutionary radiations, and sheds new light on how the intensification of the Asian monsoon has driven evolutionary radiations in large plant genera of the Himalaya-Hengduan Mountains.

Geographical isolation versus dispersal: Relictual alpine grasshoppers support a model of interglacial diversification with limited hybridization

Alpine biotas are paradigmatic of the countervailing roles of geographical isolation and dispersal during diversification. In temperate regions, repeated distributional shifts driven by Pleistocene climatic oscillations produced both recurrent pulses of population fragmentation and opportunities for gene flow during range expansions. Here, we test whether a model of divergence in isolation vs. with gene flow is more likely in the diversification of flightless alpine grasshoppers of the genus Podisma from the Iberian Peninsula. The answer to this question can also provide key insights about the pace of evolution. Specifically, if the data fit a divergence in isolation model, this suggests rapid evolution of reproductive isolation. Genomic data confirm a Pleistocene origin of the species complex, and multiple analytical approaches revealed limited asymmetric historical hybridization between two taxa. Genomic-based demographic reconstructions, spatial patterns of genetic structure and range shifts inferred from palaeodistribution modelling suggest severe range contraction accompanied by declines in effective population sizes during interglacials (i.e., contemporary populations confined to sky islands are relicts) and expansions during the coldest stages of the Pleistocene in each taxon. Although limited hybridization during secondary contact leads to phylogenetic uncertainty if gene flow is not accommodated when estimating evolutionary relationships, all species exhibit strong genetic cohesiveness. Our study lends support to the notion that the accumulation of incipient differences during periods of isolation were sufficient to lead to lineage persistence, but also that the demographic changes, dispersal constraints and spatial distribution of the sky islands themselves mediated species diversification in temperate alpine biotas.

Phylogeography of the endemic red-tailed cicadas of New Zealand (Hemiptera: Cicadidae: Rhodopsalta), and molecular, morphological and bioacoustical confirmation of the existence of Hudson’s Rhodopsalta microdora

Why do some genera radiate, whereas others do not? The genetic structure of present-day populations can provide clues for developing hypotheses. In New Zealand, three Cicadidae genera are depauperate [Amphipsalta (three species), Notopsalta (one species) and Rhodopsalta (three species)], whereas two have speciated extensively [Kikihia (~30 species/subspecies) and Maoricicada (~20 species/subspecies). Here, we examine the evolution of Rhodopsalta, the last New Zealand genus to be studied phylogenetically and phylogeographically. We use Bayesian and maximum-likelihood analyses of mitochondrial cox1 and nuclear EF1α gene sequences. Concatenated and single-gene phylogenies for 70 specimens (58 localities) support its monophyly and three described species: Rhodopsalta cruentata, Rhodopsalta leptomera and Rhodopsalta microdora, the last taxon previously regarded as uncertain. We provide distribution maps, biological notes and the first descriptions of diagnostic songs. We show that both R. cruentata and R. microdora exhibit northern and southern genetic subclades. Subclades of the dry-adapted R. microdora clade show geographical structure, whereas those of the mesic R. cruentata and sand-dune specialist R. leptomera have few discernible patterns. Genetic, bioacoustical and detailed distributional evidence for R. microdora add to the known biodiversity of New Zealand. We designate a lectotype for Tettigonia cruentata  Fabricius, 1775, the type species of Rhodopsalta.

Dynamic Capabilities and Their Impact on Intellectual Capital and Innovation Performance

There is a high tendency for conversion from a statistical economy based on measuring tangible assets into investigating non-tangible capital drive in the present economic status worldwide. The implications of intellectual capital on innovation performance have widely attracted attention among researchers in the global arena. The present study investigated the impacts of intellectual capital on innovation performance in the banking sector as influencing non-tangible assets. Besides, the role of dynamic capabilities in moderating the relationship between intellectual capital and innovation performance was examined. A purposive sampling technique was applied to 364 participants from Iraqi commercial banks as the research context. Thereafter, structural equation modelling (SEM) was utilised to analyse the collected data from the survey questionnaire using SPSS.v25 and AMOS.v24. The study found that the employees’ levels of intellectual capital significantly increased toward innovativeness through the moderating role of dynamic capabilities between intellectual capital and innovation performance in the commercial banking sector for better competitive advantages. Consequently, the study provides valuable insights and guidance for academicians and practitioners on the impacts of developing intellectual capital on enhancing competitive performance, especially in the context of Iraqi commercial banks.

Floral Scent Evolution in the Genus Jaborosa (Solanaceae): Influence of Ecological and Environmental Factors

Floral scent is a key communication channel between plants and pollinators. However, the contributions of environment and phylogeny to floral scent composition remain poorly understood. In this study, we characterized interspecific variation of floral scent composition in the genus Jaborosa Juss. (Solanaceae) and, using an ecological niche modelling approach (ENM), we assessed the environmental variables that exerted the strongest influence on floral scent variation, taking into account pollination mode and phylogenetic relationships. Our results indicate that two major evolutionary themes have emerged: (i) a 'warm Lowland Subtropical nectar-rewarding clade' with large white hawkmoth pollinated flowers that emit fragrances dominated by oxygenated aromatic or sesquiterpenoid volatiles, and (ii) a 'cool-temperate brood-deceptive clade' of largely fly-pollinated species found at high altitudes (Andes) or latitudes (Patagonian Steppe) that emit foul odors including cresol, indole and sulfuric volatiles. The joint consideration of floral scent profiles, pollination mode, and geoclimatic context helped us to disentangle the factors that shaped floral scent evolution across "pollinator climates" (geographic differences in pollinator abundance or preference). Our findings suggest that the ability of plants in the genus Jaborosa to colonize newly formed habitats during Andean orogeny was associated with striking transitions in flower scent composition that trigger specific odor-driven behaviors in nocturnal hawkmoths and saprophilous fly pollinators.

The evolutionary history of vines in a neotropical biodiversity hotspot: Phylogenomics and biogeography of a large passion flower clade (Passiflora section Decaloba)

Because of their extraordinary flower and leaf morphology, passion flowers (Passifloraceae) have fascinated naturalists since their discovery. Within the large, diverse (600 species) genus Passiflora is an especially enigmatic and species-rich (120 spp.) subclade, Section Decaloba, which occurs in the Neotropics and has its center of diversity in Andean montane forests. A recent phylogenetic study of Passifloraceae showed that Section Decaloba was monophyletic, but was unable to resolve relationships within the clade, thus preventing inferences of evolutionary history and biogeography. The goal of this study was to elucidate the phylogeny and biogeography of Section Decaloba. We sampled 206 accessions representing 91 of the ~ 120 known species in section Decaloba and four outgroups, with samples derived predominantly from herbarium specimens. We generated DNA sequences using a high-throughput DNA sequencing technique called 2b-RAD, reconstructed the phylogeny, and conducted ancestral area reconstructions to infer the biogeographic history of the group. We recovered predominantly well-supported trees in which species were grouped into two main clades: 1) the Central American clade, within which the majority of nodes well supported and species were monophyletic and 2) the South American clade, a large clade that showed overall lower resolution and included several polyphyletic species and species complexes that need additional research. RASP analysis showed that section Decaloba originated in Central America around 10.4 Ma, and then dispersed to South America, the Greater Antilles, and the Bahamas. The South American clade diversified in the Northern Andes and then dispersed to the rest of South America, and Lesser Antilles. Results suggest that both long-distance dispersal and colonization of newly available habitats (i.e., in the Andes) likely promoted diversification of this clade. This study also illustrates how using herbarium specimens and a RAD-seq approach can produce phylogenies for broadly distributed, highly diverse, and poorly accessible groups of plants where field collections would be unfeasible.

Interaction of climate change, potentially toxic elements (PTEs), and topography on plant diversity and ecosystem functions in a high-altitude mountainous region of the Tibetan Plateau

Potentially toxic elements (PTEs) generated from mining activities have affected ecological diversity and ecosystem functions around the world. Accurately assessing the long-term effects of PTEs is critical to classifying recoverable areas and proposing management strategies. Mining activities that shape geographical patterns of biodiversity in individual regions are increasingly understood, but the complex interactions on broad scales and in changing environments are still unclear. In this study, we developed a series of empirical models that simulate the changes in biodiversity and ecosystem functions in mine-affected regions along elevation gradients (1500-3600 m a.s.l) in the metal-rich Qilian Mountains (∼800 km) on the northeastern Tibetan Plateau (China). Our results confirmed the crucial role of PTEs dispersal, topography, and climatic heterogeneity in the diversification of plant community composition. On average, 54% of the changes in ecosystem functions were explained by the interactions among topography, climate, and PTEs. However, merely 30% of the changes were correlated with a single driver. The changes in species composition (explained variables = 94.8%) in the PTE-polluted habitats located in the warm and humid low-elevation deserts and grasslands were greater than those occurring in the dry alpine deserts and grasslands. The ecosystem functions (soil characteristics, nutrient migration, and plant biomass) experienced greater changes in the humid low-elevation grasslands and alpine deserts. Our results suggest that the processes driven by climate or other factors can result in high-altitude PTE-affected habitat facing greater threats.

Draft genome of Puya raimondii (Bromeliaceae), the Queen of the Andes

Puya raimondii, the Queen of the Andes, is an endangered high Andean species in the Bromeliaceae family. Here, we report its first genome to promote its conservation and evolutionary study. Comparative genomics showed P. raimondii diverged from Ananas comosus about 14.8 million years ago, and the long terminal repeats were likely to contribute to the genus diversification in last 3.5 million years. The gene families related to plant reproductive development and stress responses significantly expanded in the genome. At the same time, gene families involved in disease defense, photosynthesis and carbohydrate metabolism significantly contracted, which may be an evolutionary strategy to adapt to the harsh conditions in high Andes. The demographic history analysis revealed the P. raimondii population size sharply declined in the Pleistocene and then increased in the Holocene. We also designed and tested 46 pairs of universal primers for amplifying orthologous single-copy nuclear genes in Puya species.

The importance of environmental conditions in maintaining lineage identity in Epithelantha (Cactaceae)

The use of environmental variables to explain the evolution of lineages has gained relevance in recent studies. Additionally, it has allowed the recognition of species by adding more characters to morphological and molecular information. This study focuses on identifying environmental and landscape variables that have acted as barriers that could have influenced the evolution of Epithelantha species and its close genera.Our results show that soil pH, isothermality, temperature seasonality, and annual precipitation have a significant phylogenetic signal for Epithelantha. Soil type and landforms are also relevant as ecological barriers that maintain the identity of Epithelantha species.The variables associated with the soil (pH) have influenced the evolution of Epithelantha and probably in other genera of Cactaceae. Additionally, Epithelantha is frequent in the piedmont and haplic kastanozems. Bioclimatic variables reinforce the recognition of E. micromeris, and E. cryptica as independent species. Therefore, ecology can be considered as a factor to explain the high level of endemism in Cactaceae.

Effect of the Andean Geography and Climate on the Specialized Metabolism of Its Vegetation: The Subtribe Espeletiinae (Asteraceae) as a Case Example

The Andean mountains are 'center stage' to some of the most spectacular examples of plant diversifications, where geographic isolation and past climatic fluctuations have played a major role. However, the influence of Andean geography and climate as drivers of metabolic variation in Andean plants is poorly elucidated. Here, we studied the influence of those factors on the metabolome of the subtribe Espeletiinae (Asteraceae) using liquid chromatography coupled to high-resolution mass spectrometry data of over two hundred samples from multiple locations. Our results demonstrate that metabolic profiles can discriminate Espeletiinae taxa at different geographic scales, revealing inter- and intraspecific metabolic variations: at the country level, segregation between Colombian and Venezuelan taxa was observed; regionally, between páramo massifs; and locally, between páramo complexes. Metabolic differences in Espeletiinae were mainly explained by geographic isolation, although differences in taxonomic genera, temperature, and elevation, were also important. Furthermore, we found that different species inhabiting the same páramo complex showed stronger chemical similarities than the same species at different locations, corroborating that geographic isolation represents the main driver of metabolic change in Espeletiinae. The current study serves as a starting point to fill in the gaps in how Andean geography and climate have shaped the metabolism of its vegetation and reveal the potential of untargeted metabolomics to study the environmental physiology of plants.

Phylogenetic signatures of ecological divergence and leapfrog adaptive radiation in Espeletia

PREMISE

Events of accelerated species diversification represent one of Earth's most celebrated evolutionary outcomes. Northern Andean high-elevation ecosystems, or páramos, host some plant lineages that have experienced the fastest diversification rates, likely triggered by ecological opportunities created by mountain uplifts, local climate shifts, and key trait innovations. However, the mechanisms behind rapid speciation into the new adaptive zone provided by these opportunities have long remained unclear.

METHODS

We address this issue by studying the Venezuelan clade of Espeletia, a species-rich group of páramo-endemics showing a dazzling ecological and morphological diversity. We performed several comparative analyses to study both lineage and trait diversification, using an updated molecular phylogeny of this plant group.

RESULTS

We showed that sets of either vegetative or reproductive traits have conjointly diversified in Espeletia along different vegetation belts, leading to adaptive syndromes. Diversification in vegetative traits occurred earlier than in reproductive ones. The rate of species and morphological diversification showed a tendency to slow down over time, probably due to diversity dependence. We also found that closely related species exhibit significantly more overlap in their geographic distributions than distantly related taxa, suggesting that most events of ecological divergence occurred at close geographic proximity within páramos.

CONCLUSIONS

These results provide compelling support for a scenario of small-scale ecological divergence along multiple ecological niche dimensions, possibly driven by competitive interactions between species, and acting sequentially over time in a leapfrog pattern.

How to Tackle Phylogenetic Discordance in Recent and Rapidly Radiating Groups? Developing a Workflow Using Loricaria (Asteraceae) as an Example

A major challenge in phylogenetics and -genomics is to resolve young rapidly radiating groups. The fast succession of species increases the probability of incomplete lineage sorting (ILS), and different topologies of the gene trees are expected, leading to gene tree discordance, i.e., not all gene trees represent the species tree. Phylogenetic discordance is common in phylogenomic datasets, and apart from ILS, additional sources include hybridization, whole-genome duplication, and methodological artifacts. Despite a high degree of gene tree discordance, species trees are often well supported and the sources of discordance are not further addressed in phylogenomic studies, which can eventually lead to incorrect phylogenetic hypotheses, especially in rapidly radiating groups. We chose the high-Andean Asteraceae genus Loricaria to shed light on the potential sources of phylogenetic discordance and generated a phylogenetic hypothesis. By accounting for paralogy during gene tree inference, we generated a species tree based on hundreds of nuclear loci, using Hyb-Seq, and a plastome phylogeny obtained from off-target reads during target enrichment. We observed a high degree of gene tree discordance, which we found implausible at first sight, because the genus did not show evidence of hybridization in previous studies. We used various phylogenomic analyses (trees and networks) as well as the D-statistics to test for ILS and hybridization, which we developed into a workflow on how to tackle phylogenetic discordance in recent radiations. We found strong evidence for ILS and hybridization within the genus Loricaria. Low genetic differentiation was evident between species located in different Andean cordilleras, which could be indicative of substantial introgression between populations, promoted during Pleistocene glaciations, when alpine habitats shifted creating opportunities for secondary contact and hybridization.

Characterizing the genetic diversity of the Andean blueberry (Vaccinium floribundum Kunth.) across the Ecuadorian Highlands

The Ecuadorian páramo, a high altitude tundra-like ecosystem, is a unique source of various ecosystem services and distinct biodiversity. Anthropogenic activities are associated with its fragmentation, which alters ecological factors and directly threatens resident species. Vaccinium floribundum Kunth., commonly known as Andean blueberry or mortiño, is a wild shrub endemic to the Andean region and highly valued in Ecuador for its berries, which are widely used in food preparations and hold an important cultural value. Since it is a wild species, mortiño could be vulnerable to environmental changes, resulting in a reduction of the size and distribution of its populations. To evaluate the extent of these effects on the mortiño populations, we assessed the genetic diversity and population structure of the species along the Ecuadorian highlands. We designed and developed a set of 30 species-specific SSR (simple sequence repeats) markers and used 16 of these to characterize 100 mortiño individuals from 27 collection sites. Our results revealed a high degree of genetic diversity (HE = 0.73) for the Ecuadorian mortiño, and a population structure analyses suggested the existence of distinct genetic clusters present in the northern, central and southern highlands. A fourth, clearly differentiated cluster was also found and included individuals from locations at higher elevations. We suggest that the population structure of the species could be explained by an isolation-by-distance model and can be associated with the geological history of the Andean region. Our results suggest that elevation could also be a key factor in the differentiation of mortiño populations. This study provides an extensive overview of the species across its distribution range in Ecuador, contributing to a better understanding of its conservation status. These results can assist in the development of conservation programs for this valuable biological and cultural resource and for the páramo ecosystem as a whole.

Ecological divergence and hybridization of Neotropical Leishmania parasites

Parasites are interesting models for studying speciation processes because they have a high potential for specialization, thanks to the intimate ecological association with their hosts and vectors. Yet little is known about the circumstances under which new parasite lineages emerge. Here we studied the genome diversity of parasites of the Leishmania braziliensis species complex that inhabit both Amazonian and Andean biotas in Peru. We identify three major parasite lineages that occupy particular ecological niches and show that these emerged during forestation changes over the past 150,000 y. We furthermore discovered that meiotic recombination between Amazonian and Andean lineages resulted in full-genome hybrids presenting mixed mitochondrial genomes, providing insights into the genetic consequences of hybridization in parasitic protozoa.

The tropical Andes are an important natural laboratory to understand speciation in many taxa. Here we examined the evolutionary history of parasites of the Leishmania braziliensis species complex based on whole-genome sequencing of 67 isolates from 47 localities in Peru. We first show the origin of Andean Leishmania as a clade of near-clonal lineages that diverged from admixed Amazonian ancestors, accompanied by a significant reduction in genome diversity and large structural variations implicated in host–parasite interactions. Within the Andean species, patterns of population structure were strongly associated with biogeographical origin. Molecular clock and ecological niche modeling suggested that the history of diversification of the Andean lineages is limited to the Late Pleistocene and intimately associated with habitat contractions driven by climate change. These results suggest that changes in forestation over the past 150,000 y have influenced speciation and diversity of these Neotropical parasites. Second, genome-scale analyses provided evidence of meiotic-like recombination between Andean and Amazonian Leishmania species, resulting in full-genome hybrids. The mitochondrial genome of these hybrids consisted of homogeneous uniparental maxicircles, but minicircles originated from both parental species. We further show that mitochondrial minicircles—but not maxicircles—show a similar evolutionary pattern to the nuclear genome, suggesting that compatibility between nuclear-encoded mitochondrial genes and minicircle-encoded guide RNA genes is essential to maintain efficient respiration. By comparing full nuclear and mitochondrial genome ancestries, our data expand our appreciation on the genetic consequences of diversification and hybridization in parasitic protozoa.

Introgression across evolutionary scales suggests reticulation contributes to Amazonian tree diversity

Hybridization has the potential to generate or homogenize biodiversity and is a particularly common phenomenon in plants, with an estimated 25% of plant species undergoing interspecific gene flow. However, hybridization in Amazonia's megadiverse tree flora was assumed to be extremely rare despite extensive sympatry between closely related species, and its role in diversification remains enigmatic because it has not yet been examined empirically. Using members of a dominant Amazonian tree family (Brownea, Fabaceae) as a model to address this knowledge gap, our study recovered extensive evidence of hybridization among multiple lineages across phylogenetic scales. More specifically, using targeted sequence capture our results uncovered several historical introgression events between Brownea lineages and indicated that gene tree incongruence in Brownea is best explained by reticulation, rather than solely by incomplete lineage sorting. Furthermore, investigation of recent hybridization using ~19,000 ddRAD loci recovered a high degree of shared variation between two Brownea species that co-occur in the Ecuadorian Amazon. Our analyses also showed that these sympatric lineages exhibit homogeneous rates of introgression among loci relative to the genome-wide average, implying a lack of selection against hybrid genotypes and persistent hybridization. Our results demonstrate that gene flow between multiple Amazonian tree species has occurred across temporal scales, and contrasts with the prevailing view of hybridization's rarity in Amazonia. Overall, our results provide novel evidence that reticulate evolution influenced diversification in part of the Amazonian tree flora, which is the most diverse on Earth.

Fast diversification through a mosaic of evolutionary histories characterizes the endemic flora of ancient Neotropical mountains

Mountains are among the most biodiverse areas on the globe. In young mountain ranges, exceptional plant species richness is often associated with recent and rapid radiations linked to the mountain uplift itself. In ancient mountains, however, orogeny vastly precedes the evolution of vascular plants, so species richness has been explained by species accumulation during long periods of low extinction rates. Here we evaluate these assumptions by analysing plant diversification dynamics in the campo rupestre, an ecosystem associated with pre-Cambrian mountaintops and highlands of eastern South America, areas where plant species richness and endemism are among the highest in the world. Analyses of 15 angiosperm clades show that radiations of endemics exhibit fastest rates of diversification during the last 5 Myr, a climatically unstable period. However, results from ancestral range estimations using different models disagree on the age of the earliest in situ speciation events and point to a complex floristic assembly. There is a general trend for higher diversification rates associated with these areas, but endemism may also increase or reduce extinction rates, depending on the group. Montane habitats, regardless of their geological age, may lead to boosts in speciation rates by accelerating population isolation in archipelago-like systems, circumstances that can also result in higher extinction rates and fast species turnover, misleading the age estimates of endemic lineages.

Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S. aethnensis and S. chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome-wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly-differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs' allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16-0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species - pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.

Building mountain biodiversity: Geological and evolutionary processes

Mountain regions are unusually biodiverse, with rich aggregations of small-ranged species that form centers of endemism. Mountains play an array of roles for Earth’s biodiversity and affect neighboring lowlands through biotic interchange, changes in regional climate, and nutrient runoff. The high biodiversity of certain mountains reflects the interplay of multiple evolutionary mechanisms: enhanced speciation rates with distinct opportunities for coexistence and persistence of lineages, shaped by long-term climatic changes interacting with topographically dynamic landscapes. High diversity in most tropical mountains is tightly linked to bedrock geology—notably, areas comprising mafic and ultramafic lithologies, rock types rich in magnesium and poor in phosphate that present special requirements for plant physiology. Mountain biodiversity bears the signature of deep-time evolutionary and ecological processes, a history well worth preserving.

The rise of the Andes promoted rapid diversification in Neotropical Phlegmariurus (Lycopodiaceae)

Tropical mountains are disproportionately biodiverse relative to their surface area, but the processes underlying their exceptional diversity require further study. Here, we use comparative phylogenetic methods to examine the impact of the Andean orogeny on the diversification of Neotropical Phlegmariurus, a species-rich lycophyte clade. We generated a time-calibrated phylogeny of 105 species of Neotropical Phlegmariurus and estimated lineage diversification rates. We tested for correlations between lineage diversification rates and species range size, niche breadth, elevational range amplitude, and mean elevation of occurrence. A recently developed macroevolutionary model was used to incorporate geological data and test for an association between diversification rates and the Andean uplift. Diversification rates of Neotropical Phlegmariurus are negatively correlated with species range size and positively correlated with mean elevation of species occurrence. The rise of the Andes is strongly associated with increased rates of diversification in Neotropical Phlegmariurus during the last 10 Myr. Our study demonstrates the importance of mountain-building events and geographical isolation of alpine populations as drivers of rapid diversification, even in spore-dispersed plants. This work also highlights the usefulness of combined phylogenetic, geological and ecological datasets, and the promise of comparative environment-dependent diversification models in better understanding the evolutionary origins of biodiversity.