Recent Adaptive Radiation Research Articles

Phylogenomic analysis of the Lake Kronotskoe species flock of Dolly Varden charr reveals genetic and developmental signatures of sympatric radiation.

Recent adaptive radiations provide experimental opportunities to parse the relationship between genomic variation and the origins of distinct phenotypes. Sympatric radiations of the charr complex (genus Salvelinus) present a trove for phylogenetic analyses as charrs have repeatedly diversified into multiple morphs with distinct feeding specializations. However, charr species flocks normally comprise only two to three lineages. Dolly Varden charr inhabiting Lake Kronotskoe represent the most extensive radiation described for the genus, containing at least seven lineages, each with defining morphological and ecological traits. Here, we perform the first genome-wide analysis of this species flock to parse the foundations of adaptive change. Our data support distinct, reproductively isolated lineages within the clade. We find that changes in genes associated with thyroid signaling and craniofacial development provided a foundational shift in evolution to the lake. The thyroid axis is further implicated in subsequent lineage partitioning events. These results delineate a genetic scenario for the diversification of specialized lineages and highlight a common axis of change biasing the generation of specific forms during adaptive radiation.

Diversification of spider silk properties in an adaptive radiation of Hawaiian orb-weaving spiders

IntroductionThe design of biological structures and the materials composing those structures are intimately connected to performance in biological systems. Spider webs present an excellent example of how design and materials interact during their function in capturing prey. Major shifts in how spider webs capture prey have occurred due to evolutionary changes in both web architecture and silk properties. However, these shifts are mostly described for long timescales deep within the spider’s tree of life. Hawaiian Tetragnatha presents an opportunity to ask if such shifts can occur at much shorter timescales because web design diverges significantly among closely related species on the same island while also converging with more distant relatives on other islands. Here, we provide an initial test of whether or not silk properties diversified during the recent adaptive radiation of Hawaiian Tetragnatha.MethodsWe obtained radial and capture spiral silk from orb webs for spiders on two islands and tested their tensile and adhesive properties. We also used solution-state NMR to compare the composition of low molecular weight compounds in the glue because of their influence on capture spiral stickiness.ResultsResults showed differences in the stiffness of radial silk among four populations of Hawaiian Tetragnatha, while extensibility remained unchanged. Although not statistically different, radial strength and toughness varied twofold among species. Stickiness varied threefold among the four populations of orb weavers. No conspicuous qualitative differences in the low molecular weight compound composition of aggregate glue were found, suggesting that differences in capture spiral stickiness were due to the amount or arrangement of glue droplets on threads.DiscussionWhile our sampling is modest, our data provide the first evidence that silk properties can evolve measurably over the relatively short timescales of the adaptive radiation of Hawaiian Tetragnatha spiders.

The Bemisia tabaci cryptic (sibling) species group -imperative for a taxonomic reassessment.

The taxonomy of the Bemisia tabaci cryptic species group remains a challenge due to the lack of morphological differentiation and porous species boundaries among its members. Additionally, it is unclear whether B. tabaci consists of several species in evolutionary stasis with limited morphological change or is the result of a recent adaptive radiation characterized by great ecological diversity but little morphological divergence. Here, a historical overview of the development of the nomenclature used to classify B. tabaci is provided covering changes after synonymizing several species in 1957 until recent insights gained from whole-genome sequencing data. The article discusses the limitations of using a 3.5% mtCOI threshold and argues that a 1% nuclear divergence cutoff better reflects ecological and biogeographic species boundaries. Finally, a plan of action is outlined for naming B. tabaci species using a Latin binomial system in accordance with the International Comission on Zoological Nomenclature (ICZN) regulations.

Testing hypotheses of hybrid taxon formation in the shrubby beardtongues (Penstemon subgenus Dasanthera).

Hybridization is increasingly being identified in the genomes of species across the tree of life. This has led to a general recognition that hybridization plays an important role in the generation of species diversity. While hybridization may increase species diversity directly via the formation of new taxa through hybrid speciation, it may also act indirectly via the exchange of phenotypic and genetic variance between species, which may in turn stimulate future speciation events. Using high-throughput sequence data, we resolve phylogenetic relationships and investigate the role of hybridization as a diversification mechanism in the shrubby beardtongues (Penstemon subgenus Dasanthera), a group of North American wildflowers that has undergone a recent and rapid adaptive radiation. Specifically, we test four hypotheses of hybrid taxon formation resulting from hybridization between P. davidsonii and P. fruticosus. Species tree inference supports the monophyly of subgenus Dasanthera and elucidates relationships between taxa distributed in the Cascades and Sierra Nevada Mountains. Results also provide evidence of gene flow between P. davidsonii and P. fruticosus and support at least one hybrid origin hypothesis (P. davidsonii var. menziesii) in a region of contemporary distributional overlap. Hybridization may have also been facilitated by historical overlap in geographic distribution caused by species' responses to climatic changes during the Pleistocene. Our results support a history of hybridization between focal taxa in a rapidly radiating clade of plants and more broadly contribute to our growing understanding of the role of hybridization as a diversification mechanism in plants. This article is protected by copyright. All rights reserved.

Resolution, conflict and rate shifts: insights from a densely sampled plastome phylogeny for Rhododendron (Ericaceae).

Rhododendron is a species-rich and taxonomically challenging genus due to recent adaptive radiation and frequent hybridization. A well-resolved phylogenetic tree would help to understand the diverse history of Rhododendron in the Himalaya-Hengduan Mountains where the genus is most diverse. We reconstructed the phylogeny based on plastid genomes with broad taxon sampling, covering 161 species representing all eight subgenera and all 12 sections, including ~45 % of the Rhododendron species native to the Himalaya-Hengduan Mountains. We compared this phylogeny with nuclear phylogenies to elucidate reticulate evolutionary events and clarify relationships at all levels within the genus. We also estimated the timing and diversification history of Rhododendron, especially the two species-rich subgenera Rhododendron and Hymenanthes that comprise >90 % of Rhododendron species in the Himalaya-Hengduan Mountains. The full plastid dataset produced a well-resolved and supported phylogeny of Rhododendron. We identified 13 clades that were almost always monophyletic across all published phylogenies. The conflicts between nuclear and plastid phylogenies suggested strongly that reticulation events may have occurred in the deep lineage history of the genus. Within Rhododendron, subgenus Therorhodion diverged first at 56 Mya, then a burst of diversification occurred from 23.8 to 17.6 Mya, generating ten lineages among the component 12 clades of core Rhododendron. Diversification in subgenus Rhododendron accelerated c. 16.6 Mya and then became fairly continuous. Conversely, Hymenanthes diversification was slow at first, then accelerated very rapidly around 5 Mya. In the Himalaya-Hengduan Mountains, subgenus Rhododendron contained one major clade adapted to high altitudes and another to low altitudes, whereas most clades in Hymenanthes contained both low- and high-altitude species, indicating greater ecological plasticity during its diversification. The 13 clades proposed here may help to identify specific ancient hybridization events. This study will help to establish a stable and reliable taxonomic framework for Rhododendron, and provides insight into what drove its diversification and ecological adaption. Denser sampling of taxa, examining both organelle and nuclear genomes, is needed to better understand the divergence and diversification history of Rhododendron.

Using functional traits to identify conservation priorities for the world's crocodylians

AbstractOver half of all crocodylians are threatened with extinction, with many species considered to be functionally extinct throughout their historical ranges. How such losses impact crocodylian functional diversity (FD) and its conservation is poorly understood, but can be investigated using measurable traits. Where such information is unavailable, phylogenetic diversity (PD) has been proposed as an effective strategy to capture FD, but this assumption is not well tested.We constructed a global trait database for 28 extant crocodylians, and quantified their FD, distinctiveness and groupings. Functionally distinctive crocodylians include the most and least evolutionarily distinct species, likely because unusual trait values and combinations can be produced by both evolutionary isolation and evolutionarily recent adaptive radiation. We identified four functional groups of crocodylians, within which species have similar ecological roles. Highly distinctive species, such as the gharial, fit less well into functional groups.We combined functional distinctiveness (FUD) and threat ranking scores to prioritise species for conservation. The metric is termed Ecologically Distinct and Globally Endangered (EcoDGE), and we tested how well it operates to preserve crocodylian FD among extant Crocodylia. Under current extinction scenarios, 32%–38% of crocodylian FD might be lost within 100 years. Losses increase to 48% if all threatened species go extinct, with Asia identified as a hotspot of threatened FD in crocodylians. Highly threatened distinctive species lead this extinction scenario, exacerbated by clusters of species sharing traits linked to high extinction risk. These traits include habitat specialism and lower reproductive investment; in contrast, species that tolerate climate extremes are more resilient.Prioritising at‐risk species based on FUD and threat status appears to be an effective way to conserve present‐day crocodylian diversity. Conservation prioritisation based on the EcoDGE metric performed well to minimise projected losses in FD. Among extant crocodylians, FD losses are projected when high diversity overlaps with high extrinsic threats.We then examined the utility of PD‐based measures as surrogates for FD conservation, to better understand the advantages and drawbacks of such substitutions. A PD‐based measure such as the Evolutionarily Distinct and Globally Endangered metric provided an effective tool to capture FD in our tests.Read the freePlain Language Summaryfor this article on the Journal blog.

Allopatric origin of sympatric whitefish morphs with insights on the genetic basis of their reproductive isolation.

The European whitefish (Coregonus lavaretus) species complex is a classic example of recent adaptive radiation. Here, we examine a whitefish population introduced to northern Finnish Lake Tsahkal in the late 1960s, where three divergent morphs (viz. littoral, pelagic, and profundal feeders) were found 10 generations after. Using demographic modeling based on genomic data, we show that whitefish morphs evolved during a phase of strict isolation, refuting a rapid sympatric divergence scenario. The lake is now an artificial hybrid zone between morphs originated in allopatry. Despite their current syntopy, clear genetic differentiation remains between two of the three morphs. Using admixture mapping, we identify five SNPs associated with gonad weight variation, a proxy for sexual maturity and spawning time. We suggest that ecological adaptations in spawning time evolved in allopatry are currently maintaining partial reproductive isolation in the absence of other barriers to gene flow.

Rapid adaptive radiation of Darwin's finches depends on ancestral genetic modules.

Recent adaptive radiations are models for investigating mechanisms contributing to the evolution of biodiversity. An unresolved question is the relative importance of new mutations, ancestral variants, and introgressive hybridization for phenotypic evolution and speciation. Here, we address this issue using Darwin’s finches and investigate the genomic architecture underlying their phenotypic diversity. Admixture mapping for beak and body size in the small, medium, and large ground finches revealed 28 loci showing strong genetic differentiation. These loci represent ancestral haplotype blocks with origins predating speciation events during the Darwin’s finch radiation. Genes expressed in the developing beak are overrepresented in these genomic regions. Ancestral haplotypes constitute genetic modules for selection and act as key determinants of the unusual phenotypic diversity of Darwin’s finches. Such ancestral haplotype blocks can be critical for how species adapt to environmental variability and change.

Hybridization alters the shape of the genotypic fitness landscape, increasing access to novel fitness peaks during adaptive radiation.

Estimating the complex relationship between fitness and genotype or phenotype (i.e. the adaptive landscape) is one of the central goals of evolutionary biology. However, adaptive walks connecting genotypes to organismal fitness, speciation, and novel ecological niches are still poorly understood and processes for surmounting fitness valleys remain controversial. One outstanding system for addressing these connections is a recent adaptive radiation of ecologically and morphologically novel pupfishes (a generalist, molluscivore, and scale-eater) endemic to San Salvador Island, Bahamas. We leveraged whole-genome sequencing of 139 hybrids from two independent field fitness experiments to identify the genomic basis of fitness, estimate genotypic fitness networks, and measure the accessibility of adaptive walks on the fitness landscape. We identified 132 single nucleotide polymorphisms (SNPs) that were significantly associated with fitness in field enclosures. Six out of the 13 regions most strongly associated with fitness contained differentially expressed genes and fixed SNPs between trophic specialists; one gene (mettl21e) was also misexpressed in lab-reared hybrids, suggesting a potential intrinsic genetic incompatibility. We then constructed genotypic fitness networks from adaptive alleles and show that scale-eating specialists are the most isolated of the three species on these networks. Intriguingly, introgressed and de novo variants reduced fitness landscape ruggedness as compared to standing variation, increasing the accessibility of genotypic fitness paths from generalist to specialists. Our results suggest that adaptive introgression and de novo mutations alter the shape of the fitness landscape, providing key connections in adaptive walks circumventing fitness valleys and triggering the evolution of novelty during adaptive radiation.

Extensions and limitations of MacArthur (1958): A review of ecological and evolutionary approaches to competition and diet in the New World wood warblers (Parulidae)

AbstractThe extent to which interspecific competition structures species interactions and coexistence within communities, and the relevant mechanisms, are still debated. We focus on New World wood warblers (Parulidae), beginning with Robert MacArthur’s iconic 1958 paper in which he shows how subtle foraging behaviors, purportedly linked to dietary differences, within spruce trees contribute to the coexistence of 5 spruce-woods warbler species. MacArthur coined the phrase “resource partitioning”, and profoundly impacted the field of Ecology for subsequent decades in diverse ways. To understand what MacArthur got right and what he missed, we reviewed both ecological and evolutionary approaches to questions of the origin and coexistence of competing species in the context of diet. We argue that an important, underappreciated, mechanism of competition among coexisting migratory warbler species, particularly in winter, is diffuse exploitation competition, based in part on our own studies of warbler diets in relation to foraging behavior, substrate use, bird morphology, and other traits. Our review and synthesis of interspecific competition and coexistence in warblers have important consequences, including our questioning of the importance and effectiveness of resource partitioning in birds. We also suggest a novel hypothesis for the success of warblers today in the Caribbean and other habitats, beginning with their relatively recent adaptive radiation and the ecological opportunity on Caribbean islands.

The Build-Up of Population Genetic Divergence along the Speciation Continuum during a Recent Adaptive Radiation of Rhagoletis Flies.

New species form through the evolution of genetic barriers to gene flow between previously interbreeding populations. The understanding of how speciation proceeds is hampered by our inability to follow cases of incipient speciation through time. Comparative approaches examining different diverging taxa may offer limited inferences, unless they fulfill criteria that make the comparisons relevant. Here, we test for those criteria in a recent adaptive radiation of the Rhagoletis pomonella species group (RPSG) hypothesized to have diverged in sympatry via adaptation to different host fruits. We use a large-scale population genetic survey of 1568 flies across 33 populations to: (1) detect on-going hybridization, (2) determine whether the RPSG is derived from the same proximate ancestor, and (3) examine patterns of clustering and differentiation among sympatric populations. We find that divergence of each in-group RPSG taxon is occurring under current gene flow, that the derived members are nested within the large pool of genetic variation present in hawthorn-infesting populations of R. pomonella, and that sympatric population pairs differ markedly in their degree of genotypic clustering and differentiation across loci. We conclude that the RPSG provides a particularly robust opportunity to make direct comparisons to test hypotheses about how ecological speciation proceeds despite on-going gene flow.

Chromosomal-Level Assembly of Antarctic Scaly Rockcod, Trematomus loennbergii Genome Using Long-Read Sequencing and Chromosome Conformation Capture (Hi-C) Technologies

Trematomus species (suborder Notothenioidei; family Nototheniidae) are widely distributed in the southern oceans near Antarctica. There are 11 recognized species in the genus Trematomus, and notothenioids are known to have high chromosomal diversity (2n = 24–58) because of relatively recent and rapid adaptive radiation. Herein, we report the chromosomal-level genome assembly of T. loennbergii, the first characterized genome representative of the genus Trematomus. The final genome assembly of T. loennbergii was obtained using a Pacific Biosciences long-read sequencing platform and high-throughput chromosome conformation capture technology. Twenty-three chromosomal-level scaffolds were assembled to 940 Mb in total size, with a longest contig size of 48.5 Mb and contig N50 length of 24.7 Mb. The genome contained 42.03% repeat sequences, and a total of 24,525 protein-coding genes were annotated. We produced a high-quality genome assembly of T. loennbergii. Our results provide a first reference genome for the genus Trematomus and will serve as a basis for studying the molecular taxonomy and evolution of Antarctic fish.

Divergent diapause life history timing drives both allochronic speciation and reticulate hybridization in an adaptive radiation of Rhagoletis flies

Divergent adaptation to new ecological opportunities can be an important factor initiating speciation. However, as niches are filled during adaptive radiations, trait divergence driving reproductive isolation between sister taxa may also result in trait convergence with more distantly related taxa, increasing the potential for reticulated gene flow across the radiation. Here, we demonstrate such a scenario in a recent adaptive radiation of Rhagoletis fruit flies, specialized on different host plants. Throughout this radiation, shifts to novel hosts are associated with changes in diapause life history timing, which act as "magic traits" generating allochronic reproductive isolation and facilitating speciation-with-gene-flow. Evidence from laboratory rearing experiments measuring adult emergence timing and genome-wide DNA-sequencing surveys supported allochronic speciation between summer-fruiting Vaccinium spp.-infesting Rhagoletis mendax and its hypothesized and undescribed sister taxon infesting autumn-fruiting sparkleberries. The sparkleberry fly and R.mendax were shown to be genetically discrete sister taxa, exhibiting no detectable gene flow and allochronically isolated by a 2-month average difference in emergence time corresponding to host availability. At sympatric sites across the southern USA, the later fruiting phenology of sparkleberries overlaps with that of flowering dogwood, the host of another more distantly related and undescribed Rhagoletis taxon. Laboratory emergence data confirmed broadly overlapping life history timing and genomic evidence supported on-going gene flow between sparkleberry and flowering dogwood flies. Thus, divergent phenological adaptation can drive the initiation of reproductive isolation, while also enhancing genetic exchange across broader adaptive radiations, potentially serving as a source of novel genotypic variation and accentuating further diversification.

Extraordinary diversification of the “bristle flies” (Diptera: Tachinidae) and its underlying causes

AbstractThe family Tachinidae (“bristle flies”) is the most diverse and ecologically important group of insect parasitoids outside the parasitic wasps. It is among the most species rich families of flies (Diptera) and has experienced a recent adaptive radiation across the globe. We make use of a molecular phylogeny of the family to examine its rapid radiation and explore the traits of tachinid lineages that may have contributed to variation in their diversification. We apply a range of diversification analyses to assess the consistency and robustness of effects. We find that the Tachinidae are among the most rapidly diversifying families of animals. Six to eight clades of bristle flies, distributed across the phylogeny, exhibit strong evidence of accelerated diversification. Our results suggest that the use of holometabolous insect larvae, and specifically caterpillars (Lepidoptera), as hosts, is associated with increased diversification rates. However, these effects were inconsistent across analyses. We detected little influence of oviposition strategy (egg type) or host feeding habit, and we recovered evidence that unmeasured “hidden” traits may explain greater variance in diversification. We evaluated the strengths and weaknesses of different Maximum Likelihood and Bayesian approaches for analysing diversification and the potential for extrinsic factors, such as geography, to influence patterns of richness and diversification. In general, we conclude that although certain traits may provide opportunities for diversification, whether this is capitalized on may depend on additional traits and/or historical contingency.

Biased Learning of Sexual Signals by Female Bengalese Finches

Peak shift in mate preference learning can be a driver of rapid repeated speciation. Therefore, clades that have undergone recent adaptive radiations are predicted to show biased learning of signals from the opposite sex. The estrildid finches are one such clade. In species including the Zebra Finch Taeniopygia guttata and the Bengalese Finch Lonchura striata var. domestica, females choose mates in part based on their songs. Consistent with theory, female Zebra Finches show peak shift in their learned response to male song characteristics. We used operant conditioning to train female Bengalese Finches to respond to songs with trills of one length and to ignore songs with trills of another length. Then, we exposed those females to songs with a range of trill lengths, and we observed their responses. We found that at least some Bengalese Finches also show behaviour consistent with peak shift in their response to male songs. Moreover, females evaluated songs relative to other songs they had recently heard. Our results suggest that females respond to male sexual signals with bias in multiple species in the rapidly speciating estrildid clade.

Ancestral Hybridization Facilitated Species Diversification in the Lake Malawi Cichlid Fish Adaptive Radiation.

The adaptive radiation of cichlid fishes in East African Lake Malawi encompasses over 500 species that are believed to have evolved within the last 800,000 years from a common founder population. It has been proposed that hybridization between ancestral lineages can provide the genetic raw material to fuel such exceptionally high diversification rates, and evidence for this has recently been presented for the Lake Victoria region cichlid superflock. Here, we report that Lake Malawi cichlid genomes also show evidence of hybridization between two lineages that split 3–4 Ma, today represented by Lake Victoria cichlids and the riverine Astatotilapia sp. “ruaha blue.” The two ancestries in Malawi cichlid genomes are present in large blocks of several kilobases, but there is little variation in this pattern between Malawi cichlid species, suggesting that the large-scale mosaic structure of the genomes was largely established prior to the radiation. Nevertheless, tens of thousands of polymorphic variants apparently derived from the hybridization are interspersed in the genomes. These loci show a striking excess of differentiation across ecological subgroups in the Lake Malawi cichlid assemblage, and parental alleles sort differentially into benthic and pelagic Malawi cichlid lineages, consistent with strong differential selection on these loci during species divergence. Furthermore, these loci are enriched for genes involved in immune response and vision, including opsin genes previously identified as important for speciation. Our results reinforce the role of ancestral hybridization in explosive diversification by demonstrating its significance in one of the largest recent vertebrate adaptive radiations.

Phylogenomic Analyses Clarify True Species within the Butterfly Genus Speyeria despite Evidence of a Recent Adaptive Radiation.

When confronted with an adaptive radiation, considerable evidence is needed to resolve the evolutionary relationships of these closely related lineages. The North American genus Speyeria is one especially challenging radiation of butterflies due to potential signs of incomplete lineage sorting, ongoing hybridization, and similar morphological characters between species. Previous studies have found species to be paraphyletic and have been unable to disentangle taxa, often due to a lack of data and/or incomplete sampling. As a result, Speyeria remains unresolved. To achieve phylogenetic resolution of the genus, we conducted phylogenomic and population genomic analyses of all currently recognized North American Speyeria species, as well as several subspecies, using restriction-site-associated DNA sequencing (RADseq). Together, these analyses confirm the 16 canonical species, and clarify many internal relationships. However, a few relationships within Speyeria were poorly supported depending on the evolutionary model applied. This lack of resolution among certain taxa corroborates Speyeria is experiencing an ongoing adaptive radiation, with incomplete lineage sorting and lack of postzygotic reproductive barriers contributing to hybridization and further ambiguity. Given that many Speyeria taxa are under duress from anthropogenic factors, their legal protection must be viewed cautiously and on a case by case basis in order to properly conserve the diversity being generated.

Hybrid gene misregulation in multiple developing tissues within a recent adaptive radiation of Cyprinodon pupfishes.

Genetic incompatibilities constitute the final stages of reproductive isolation and speciation, but little is known about incompatibilities that occur within recent adaptive radiations among closely related diverging populations. Crossing divergent species to form hybrids can break up coadapted variation, resulting in genetic incompatibilities within developmental networks shaping divergent adaptive traits. We crossed two closely related sympatric Cyprinodon pupfish species–a dietary generalist and a specialized molluscivore–and measured expression levels in their F1 hybrids to identify regulatory variation underlying the novel craniofacial morphology found in this recent microendemic adaptive radiation. We extracted mRNA from eight day old whole-larvae tissue and from craniofacial tissues dissected from 17–20 day old larvae to compare gene expression between a total of seven F1 hybrids and 24 individuals from parental species populations. We found 3.9% of genes differentially expressed between generalists and molluscivores in whole-larvae tissues and 0.6% of genes differentially expressed in craniofacial tissue. We found that 2.1% of genes were misregulated in whole-larvae hybrids whereas 19.1% of genes were misregulated in hybrid craniofacial tissues, after correcting for sequencing biases. We also measured allele specific expression across 15,429 heterozygous sites to identify putative compensatory regulatory mechanisms underlying differential expression between generalists and molluscivores. Together, our results highlight the importance of considering misregulation as an early indicator of genetic incompatibilities in the context of rapidly diverging adaptive radiations and suggests that compensatory regulatory divergence drives hybrid gene misregulation in developing tissues that give rise to novel craniofacial traits.

The behavioral origins of novelty: did increased aggression lead to scale-eating in pupfishes?

Behavioral changes in a new environment are often assumed to precede the origins of evolutionary novelties. Here, we examined whether an increase in aggression is associated with a novel scale-eating trophic niche within a recent radiation of Cyprinodon pupfishes endemic to San Salvador Island, Bahamas. We measured aggression using multiple behavioral assays and used transcriptomic analyses to identify differentially expressed genes in aggression and other behavioral pathways across 3 sympatric species in the San Salvador radiation (generalist, snail-eating specialist, and scale-eating specialist) and 2 generalist outgroups. Surprisingly, we found increased behavioral aggression and differential expression of aggression-related pathways in both the scale-eating and snail-eating specialists, despite their independent evolutionary origins. Increased behavioral aggression varied across both sex and stimulus context in both species. Our results indicate that aggression is not unique to scale-eating specialists. Instead, selection may increase aggression in other contexts such as niche specialization in general or mate competition. Alternatively, increased aggression may result from indirect selection on craniofacial traits, pigmentation, or metabolism-all traits which are highly divergent, exhibit signs of selective sweeps, and are affected by aggression-related genetic pathways which are differentially expressed in this system. In conclusion, the evolution of a novel predatory trophic niche within a recent adaptive radiation does not have clear-cut behavioral origins as previously assumed, highlighting the multivariate nature of adaptation and the complex integration of behavior with other phenotypic traits.

Whole-genome sequences of Malawi cichlids reveal multiple radiations interconnected by gene flow

The hundreds of cichlid fish species in Lake Malawi constitute the most extensive recent vertebrate adaptive radiation. Here we characterize its genomic diversity by sequencing 134 individuals covering 73 species across all major lineages. The average sequence divergence between species pairs is only 0.1–0.25%. These divergence values overlap diversity within species, with 82% of heterozygosity shared between species. Phylogenetic analyses suggest that diversification initially proceeded by serial branching from a generalist Astatotilapia-like ancestor. However, no single species tree adequately represents all species relationships, with evidence for substantial gene flow at multiple times. Common signatures of selection on visual and oxygen transport genes shared by distantly related deep-water species point to both adaptive introgression and independent selection. These findings enhance our understanding of genomic processes underlying rapid species diversification, and provide a platform for future genetic analysis of the Malawi radiation.