Divergence and gene flow among Darwin's finches: A genome-wide view of adaptive radiation driven by interspecies allele sharing

Endless forms most frustrating: disentangling species boundaries in the Ramalina decipiens group (Lecanoromycetes, Ascomycota), with the description of six new species and a key to the group

Oceanic islands have been recognized as natural laboratories in which to study a great variety of evolutionary processes. One such process is evolutionary radiations, the diversification of a single ancestor into a number of species that inhabit different environments and differ in the traits that allow them to exploit those environments. The factors that drive evolutionary radiations have been studied for decades in charismatic organisms such as birds or lizards, but are lacking in lichen-forming fungi, despite recent reports of some lineages showing diversification patterns congruent with radiation. Here we propose the Ramalina decipiens group as a model system in which to carry out such studies. This group is currently thought to be comprised of five saxicolous species, all of them endemic to the Macaronesian region (the Azores, Madeira, Selvagens, Canary and Cape Verde islands). Three species are single-island endemics (a rare geographic distribution pattern in lichens), whereas two are widespread and show extreme morphological variation. The latter are suspected to harbor unrecognized species-level lineages. In order to use the Ramalina decipiens group as a model system it is necessary to resolve the group's phylogeny and to clarify its species boundaries. In this study we attempt to do so following an integrative taxonomy approach. We constructed a phylogenetic tree based on six molecular markers, four of which are newly developed and generated competing species hypotheses based on molecular (species discovery strategies based on both single locus and multilocus datasets) and phenotypic data (unsupervised clustering algorithms based on morphology, secondary chemistry and geographic origin). We found that taxonomic diversity in the Ramalina decipiens group has been highly underestimated in previous studies. In consequence, we describe six new species, most of them single-island endemics and provide a key to the group. Phylogenetic relationships among species have been reconstructed with almost full support which, coupled with the endemic character of the group, makes it an excellent system for the study of island radiations in lichen-forming fungi. Citation: Blázquez M, Pérez-Vargas I, Garrido-Benavent I, et al. 2024. Endless forms most frustrating: disentangling species boundaries in the Ramalina decipiens group (Lecanoromycetes, Ascomycota), with the description of six new species and a key to the group. Persoonia 52: 44-93. https://doi.org/10.3767/persoonia.2024.52.03 .

Evolutionary genomics of oceanic island radiations

A recurring feature of oceanic archipelagos is the presence of adaptive radiations that generate endemic, species-rich clades that can offer outstanding insight into the links between ecology and evolution. Recent developments in evolutionary genomics have contributed towards solving long-standing questions at this interface. Using a comprehensive literature search, we identify studies spanning 19 oceanic archipelagos and 110 putative adaptive radiations, but find that most of these radiations have not yet been investigated from an evolutionary genomics perspective. Our review reveals different gaps in knowledge related to the lack of implementation of genomic approaches, as well as undersampled taxonomic and geographic areas. Filling those gaps with the required data will help to deepen our understanding of adaptation, speciation, and other evolutionary processes.

DNA methylation differences during development distinguish sympatric morphs of Arctic charr (Salvelinus alpinus)

Changes in DNA methylation in specific coding or non-coding regions can influence development and potentially divergence in traits within species and groups. While the impact of epigenetic variation on developmental pathways associated with evolutionary divergence is the focus of intense investigation, few studies have looked at recently diverged systems. Phenotypic diversity between closely related populations of Arctic charr (Salvelinus alpinus), which diverged within the last 10 000 years, offers an interesting ecological model to address such effects. Using bisulfite sequencing, we studied general DNA methylation patterns during development in the four sympatric morphs of Arctic charr from Lake Thingvallavatn. The data revealed strong differences between developmental timepoints and between morphs (mainly along the benthic - limnetic axis), both at single CpG sites and in 1,000bp-regions. Genes located close to differentially methylated CpG sites were involved in nucleosome assembly, regulation of osteoclast differentiation, and cell-matrix adhesion. Differentially methylated regions were enriched in tRNA and rRNA sequences, and half of them were located close to transcription start sites. The expression of 14 genes showing methylation differences over time or between morphs was further investigated by qPCR and nine of these were found to be differentially expressed between morphs. Four genes (ARHGEF37-like, H3-like, MPP3 and MEGF9) showed a correlation between methylation and expression. Lastly, histone gene clusters displayed interesting methylation differences between timepoints and morphs, as well as intragenic methylation variation. The results presented here provide a motivation for further studies on the contribution of epigenetic traits, such as DNA methylation, to phenotypic diversity and developmental mechanisms.

Low genome-wide divergence between two lizard populations with high adaptive phenotypic differentiation

Usually, adaptive phenotypic differentiation is paralleled by genetic divergence between locally adapted populations. However, adaptation can also happen in a scenario of nonsignificant genetic divergence due to intense gene flow and/or recent differentiation. While this phenomenon is rarely published, findings on incipient ecologically driven divergence or isolation by adaptation are relatively common, which could confound our understanding about the frequency at which they actually occur in nature. Here, we explore genome-wide traces of divergence between two populations of the lacertid lizard Psammodromus algirus separated by a 600 m elevational gradient. These populations seem to be differentially adapted to their environments despite showing low levels of genetic differentiation (according to previously studies of mtDNA and microsatellite data). We performed a search for outliers (i.e., loci subject to selection) trying to identify specific loci with FST statistics significantly higher than those expected on the basis of overall, genome-wide estimates of genetic divergence. We find that local phenotypic adaptation (in terms of a wide diversity of characters) was not accompanied by genome-wide differentiation, even when we maximized the chances of unveiling such differentiation at particular loci with FST-based outlier detection tests. Instead, our analyses confirmed the lack of genome-wide differentiation on the basis of more than 70,000 SNPs, which is concordant with a scenario of local adaptation without isolation by environment. Our results add evidence to previous studies in which local adaptation does not lead to any kind of isolation (or early stages of ecological speciation), but maintains phenotypic divergence despite the lack of a differentiated genomic background.

Hybridization fluctuates with rainfall in Darwin’s tree finches

Hybridization in natural populations may be an adaptive response to shifting climatic regimes, but understanding this can be limited by the timing of sampling effort and confident identification of hybrids. On the Galapagos Islands, Darwin’s finches regularly hybridize; the islands also show extreme annual variation in rainfall, but the effect of annual rainfall on the frequency of finch hybridization is little known. Across a 20-year period on Floreana Island, we compare patterns of hybridization in sympatric Darwin’s tree finches (N = 425; Camaryhnchus spp.) and test for an effect of annual rainfall on (1) the frequency of hybrids (C. pauper × C. parvulus) and (2) the percentage of male hybrid birds produced per year (hybrid recruitment). Annual rainfall correlated with recruitment positively for hybrids, negatively for C. parvulus and not at all for C. pauper. Furthermore, the percentage of hybrids (range: 12–56%) and C. parvulus did not change with sampling year, but the critically endangered C. pauper declined. Our findings indicate that hybrid recruitment is recurring and variable according to annual rainfall in Camarhynchus Darwin’s finches.

Female-biased dispersal and non-random gene flow of MC1R variants do not result in a migration load in barn owls

Non-random gene flow is a widely neglected force in evolution and ecology. This genotype-dependent dispersal is difficult to assess, yet can impact the genetic variation of natural populations and their fitness. In this work, we demonstrate a high immigration rate of barn owls (Tyto alba) inside a Swiss population surveyed during 15 years. Using ten microsatellite loci as an indirect method to characterize dispersal, two-third of the genetic tests failed to detect a female-biased dispersal, and Monte Carlo simulations confirmed a low statistical power to detect sex-biased dispersal in case of high dispersal rate of both sexes. The capture-recapture data revealed a female-biased dispersal associated with an excess of heterozygote for the melanocortin-1 receptor gene (MC1R), which is responsible for their ventral rufous coloration. Thus, female homozygotes for the MC1R_WHITE allele might be negatively selected during dispersal. Despite the higher immigration of females that are heterozygote at MC1R, non-random gene flow should not lead to a migration load regarding this gene because we did not detect an effect of MC1R on survival and reproductive success in our local population. The present study highlights the usefulness of using multiple methods to correctly decrypt dispersal and gene flow. Moreover, despite theoretical expectations, we show that non-random dispersal of particular genotypes does not necessarily lead to migration load in recipient populations.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Biological action in Read-Write genome evolution

Many of the most important evolutionary variations that generated phenotypic adaptations and originated novel taxa resulted from complex cellular activities affecting genome content and expression. These activities included (i) the symbiogenetic cell merger that produced the mitochondrion-bearing ancestor of all extant eukaryotes, (ii) symbiogenetic cell mergers that produced chloroplast-bearing ancestors of photosynthetic eukaryotes, and (iii) interspecific hybridizations and genome doublings that generated new species and adaptive radiations of higher plants and animals. Adaptive variations also involved horizontal DNA transfers and natural genetic engineering by mobile DNA elements to rewire regulatory networks, such as those essential to viviparous reproduction in mammals. In the most highly evolved multicellular organisms, biological complexity scales with 'non-coding' DNA content rather than with protein-coding capacity in the genome. Coincidentally, 'non-coding' RNAs rich in repetitive mobile DNA sequences function as key regulators of complex adaptive phenotypes, such as stem cell pluripotency. The intersections of cell fusion activities, horizontal DNA transfers and natural genetic engineering of Read-Write genomes provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Adaptive introgression from distant Caribbean islands contributed to the diversification of a microendemic adaptive radiation of trophic specialist pupfishes

Rapid diversification often involves complex histories of gene flow that leave variable and conflicting signatures of evolutionary relatedness across the genome. Identifying the extent and source of variation in these evolutionary relationships can provide insight into the evolutionary mechanisms involved in rapid radiations. Here we compare the discordant evolutionary relationships associated with species phenotypes across 42 whole genomes from a sympatric adaptive radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas and several outgroup pupfish species in order to understand the rarity of these trophic specialists within the larger radiation of Cyprinodon. 82% of the genome depicts close evolutionary relationships among the San Salvador Island species reflecting their geographic proximity, but the vast majority of variants fixed between specialist species lie in regions with discordant topologies. Top candidate adaptive introgression regions include signatures of selective sweeps and adaptive introgression of genetic variation from a single population in the northwestern Bahamas into each of the specialist species. Hard selective sweeps of genetic variation on San Salvador Island contributed 5 times more to speciation of trophic specialists than adaptive introgression of Caribbean genetic variation; however, four of the 11 introgressed regions came from a single distant island and were associated with the primary axis of oral jaw divergence within the radiation. For example, standing variation in a proto-oncogene (ski) known to have effects on jaw size introgressed into one San Salvador Island specialist from an island 300 km away approximately 10 kya. The complex emerging picture of the origins of adaptive radiation on San Salvador Island indicates that multiple sources of genetic variation contributed to the adaptive phenotypes of novel trophic specialists on the island. Our findings suggest that a suite of factors, including rare adaptive introgression, may be necessary for adaptive radiation in addition to ecological opportunity.

The effects of Pleistocene climate change on biotic differentiation in a montane songbird clade from Wallacea

The role of Pleistocene Ice Age in tropical diversification is poorly understood, especially in archipelagos, in which glaciation-induced sea level fluctuations may lead to complicated changes in land distribution. To assess how Pleistocene land bridges may have facilitated gene flow in tropical archipelagos, we investigated patterns of diversification in the rarely-collected rusty-bellied fantail Rhipidura teysmanni (Passeriformes: Rhipiduridae) complex from Wallacea using a combination of bioacoustic traits and whole-genome sequencing methods (dd-RADSeq). We report a biogeographic leapfrog pattern in the vocalizations of these birds, and uncover deep genomic divergence among island populations despite the presence of intermittent land connections between some. We demonstrate how rare instances of genetic introgression have affected the evolution of this species complex, and document the presence of double introgressive mitochondrial sweeps, highlighting the dangers of using only mitochondrial DNA in evolutionary research. By applying different tree inference approaches, we demonstrate how concatenation methods can give inaccurate results when investigating divergence in closely-related taxa. Our study highlights high levels of cryptic avian diversity in poorly-explored Wallacea, elucidates complex patterns of Pleistocene climate-mediated diversification in an elusive montane songbird, and suggests that Pleistocene land bridges may have accounted for limited connectivity among montane Wallacean biota.

Genetic divergence with ongoing gene flow is maintained by the use of different hosts in phytophagous ladybird beetles genus Henosepilachna

Adaptation to different environments can promote population divergence via natural selection even in the presence of gene flow - a phenomenon that typically occurs during ecological speciation. To elucidate how natural selection promotes and maintains population divergence during speciation, we investigated the population genetic structure, degree of gene flow and heterogeneous genomic divergence in three closely related Japanese phytophagous ladybird beetles: Henosepilachna pustulosa, H. niponica and H. yasutomii. These species act as a generalist, a wild thistle (Cirsium spp.) specialist and a blue cohosh (Caulophyllum robustum) specialist, respectively, and their ranges differ accordingly. The two specialist species widely co-occur but are reproductively isolated solely due to their high specialization to a particular host plant. Genomewide amplified fragment-length polymorphism (AFLP) markers and mitochondrial cytochrome c oxidase subunit I (COI) gene sequences demonstrated obvious genomewide divergence associated with both geographic distance and ecological divergence. However, a hybridization assessment for both AFLP loci and the mitochondrial sequences revealed a certain degree of unidirectional gene flow between the two sympatric specialist species. Principal coordinates analysis (PCoA) based on all of the variable AFLP loci demonstrated that there are genetic similarities between populations from adjacent localities irrespective of the species (i.e. host range). However, a further comparative genome scan identified a few fractions of loci representing approximately 1% of all loci as different host-associated outliers. These results suggest that these three species had a complex origin, which could be obscured by current gene flow, and that ecological divergence can be maintained with only a small fraction of the genome is related to different host use even when there is a certain degree of gene flow between sympatric species pairs.

"Bringing Taxonomy to the Service of Genetics": Edgar Anderson and Introgressive Hybridization

In introgressive hybridization (the repeated backcrossing of hybrids with parental populations), Edgar Anderson found a source for variation upon which natural selection could work. In his 1953 review article "Introgressive Hybridization," he asserted that he was "bringing taxonomy to the service of genetics" whereas distinguished colleagues such as Theodosius Dobzhansky and Ernst Mayr did the precise opposite. His work as a geneticist particularly focused on linkage and recombination and was enriched by collaborations with Missouri Botanical Garden colleagues interested in taxonomy as well as with cytologists C.D. Darlington and Karl Sax. As the culmination of a biosystemtatic research program, Anderson's views challenged the mainstream of the Evolutionary Synthesis.

Patterns of gene flow and selection across multiple species of Acrocephalus warblers: footprints of parallel selection on the Z chromosome

Background

Understanding the mechanisms and selective forces leading to adaptive radiations and origin of biodiversity is a major goal of evolutionary biology. Acrocephalus warblers are small passerines that underwent an adaptive radiation in the last approximately 10 million years that gave rise to 37 extant species, many of which still hybridize in nature. Acrocephalus warblers have served as model organisms for a wide variety of ecological and behavioral studies, yet our knowledge of mechanisms and selective forces driving their radiation is limited. Here we studied patterns of interspecific gene flow and selection across three European Acrocephalus warblers to get a first insight into mechanisms of radiation of this avian group.

Results

We analyzed nucleotide variation at eight nuclear loci in three hybridizing Acrocephalus species with overlapping breeding ranges in Europe. Using an isolation-with-migration model for multiple populations, we found evidence for unidirectional gene flow from A. scirpaceus to A. palustris and from A. palustris to A. dumetorum. Gene flow was higher between genetically more closely related A. scirpaceus and A. palustris than between ecologically more similar A. palustris and A. dumetorum, suggesting that gradual accumulation of intrinsic barriers rather than divergent ecological selection are more efficient in restricting interspecific gene flow in Acrocephalus warblers. Although levels of genetic differentiation between different species pairs were in general not correlated, we found signatures of apparently independent instances of positive selection at the same two Z-linked loci in multiple species.

Conclusions

Our study brings the first evidence that gene flow occurred during Acrocephalus radiation and not only between sister species. Interspecific gene flow could thus be an important source of genetic variation in individual Acrocephalus species and could have accelerated adaptive evolution and speciation rate in this avian group by creating novel genetic combinations and new phenotypes. Independent instances of positive selection at the same loci in multiple species indicate an interesting possibility that the same loci might have contributed to reproductive isolation in several speciation events.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0692-2) contains supplementary material, which is available to authorized users.

Evolvability Is an Evolved Ability: The Coding Concept as the Arch-Unit of Natural Selection

Physical processes that characterize living matter are qualitatively distinct in that they involve encoding and transfer of specific types of information. Such information plays an active part in the control of events that are ultimately linked to the capacity of the system to persist and multiply. This algorithmicity of life is a key prerequisite for its Darwinian evolution, driven by natural selection acting upon stochastically arising variations of the encoded information. The concept of evolvability attempts to define the total capacity of a system to evolve new encoded traits under appropriate conditions, i.e., the accessible section of total morphological space. Since this is dependent on previously evolved regulatory networks that govern information flow in the system, evolvability itself may be regarded as an evolved ability. The way information is physically written, read and modified in living cells (the "coding concept") has not changed substantially during the whole history of the Earth's biosphere. This biosphere, be it alone or one of many, is, accordingly, itself a product of natural selection, since the overall evolvability conferred by its coding concept (nucleic acids as information carriers with the "rulebook of meanings" provided by codons, as well as all the subsystems that regulate various conditional information-reading modes) certainly played a key role in enabling this biosphere to survive up to the present, through alterations of planetary conditions, including at least five catastrophic events linked to major mass extinctions. We submit that, whatever the actual prebiotic physical and chemical processes may have been on our home planet, or may, in principle, occur at some time and place in the Universe, a particular coding concept, with its respective potential to give rise to a biosphere, or class of biospheres, of a certain evolvability, may itself be regarded as a unit (indeed the arch-unit) of natural selection.