A genomic basis of vocal rhythm in birds.

Vocal rhythm plays a fundamental role in sexual selection and species recognition in birds, but little is known of its genetic basis due to the confounding effect of vocal learning in model systems. Uncovering its genetic basis could facilitate identifying genes potentially important in speciation. Here we investigate the genomic underpinnings of rhythm in vocal non-learning Pogoniulus tinkerbirds using 135 individual whole genomes distributed across a southern African hybrid zone. We find rhythm speed is associated with two genes that are also known to affect human speech, Neurexin-1 and Coenzyme Q8A. Models leveraging ancestry reveal these candidate loci also impact rhythmic stability, a trait linked with motor performance which is an indicator of quality. Character displacement in rhythmic stability suggests possible reinforcement against hybridization, supported by evidence of asymmetric assortative mating in the species producing faster, more stable rhythms. Because rhythm is omnipresent in animal communication, candidate genes identified here may shape vocal rhythm across birds and other vertebrates.

Social antagonism facilitates supergene expansion in ants.

Antagonistic selection has long been considered a major driver of the formation and expansion of sex chromosomes. For example, sexually antagonistic variation on an autosome can select for suppressed recombination between that autosome and the sex chromosome, leading to a neo-sex chromosome. Autosomal supergenes, chromosomal regions containing tightly linked variants affecting the same complex trait, share similarities with sex chromosomes, raising the possibility that sex chromosome evolution models can explain the evolution of genome structure and recombination in other contexts. We tested this premise in a Formica ant species, wherein we identified four supergene haplotypes on chromosome 3 underlying colony social organization and sex ratio. We discovered a novel rearranged supergene variant (9r) on chromosome 9 underlying queen miniaturization. The 9r is in strong linkage disequilibrium with one chromosome 3 haplotype (P2) found in multi-queen (polygyne) colonies. We suggest that queen miniaturization is strongly disfavored in the single-queen (monogyne) background and is thus socially antagonistic. As such, divergent selection experienced by ants living in alternative social "environments" (monogyne and polygyne) may have contributed to the emergence of a genetic polymorphism on chromosome 9 and associated queen-size dimorphism. Consequently, an ancestral polygyne-associated haplotype may have expanded to include the polymorphism on chromosome 9, resulting in a larger region of suppressed recombination spanning two chromosomes. This process is analogous to the formation of neo-sex chromosomes and consistent with models of expanding regions of suppressed recombination. We propose that miniaturized queens, 16%-20% smaller than queens without 9r, could be incipient intraspecific social parasites.

Sampling affects population genetic inference: A case study of the Allen's (Selasphorus sasin) and rufous hummingbird (Selasphorus rufus).

Gene flow can affect evolutionary inference when species are undersampled. Here, we evaluate the effects of gene flow and geographic sampling on demographic inference of 2 hummingbirds that hybridize, Allen's hummingbird (Selasphorus sasin) and rufous hummingbird (Selasphorus rufus). Using whole-genome data and extensive geographic sampling, we find widespread connectivity, with introgression far beyond the Allen's × rufous hybrid zone, although the Z chromosome resists introgression beyond the hybrid zone. We test alternative hypotheses of speciation history of Allen's, rufous, and Calliope (S. calliope) hummingbird and find that rufous hummingbird is the sister taxon to Allen's hummingbird, and Calliope hummingbird is the outgroup. A model treating the 2 subspecies of Allen's hummingbird as a single panmictic population fit observed genetic data better than models treating the subspecies as distinct populations, in contrast to morphological and behavioral differences and analyses of spatial population structure. With additional sampling, our study builds upon recent studies that came to conflicting conclusions regarding the evolutionary histories of these 2 species. Our results stress the importance of thorough geographic sampling when assessing demographic history in the presence of gene flow.

Hybridization in late stages of speciation: Strong but incomplete genome-wide reproductive isolation and 'large Z-effect' in a moving hybrid zone.

In organisms reproducing sexually, speciation occurs when increasing divergence results in pre- or post-zygotic reproductive isolation between lineages. Studies focusing on reproductive isolation origin in early stages of speciation are common and many rely on genomic scans to infer introgression providing limited information on the genomic architecture of reproductive isolation long-term maintenance. This study analyses a natural hybrid zone between two species in a late stage of speciation. We used ddRADseq genotyping in the contact between Podarcis bocagei and P. carbonelli to examine admixture extent, analyse hybrid zone stability and assess genome-wide variation in selection against introgression. We confirmed strong but incomplete reproductive isolation in a bimodal hybrid zone. New findings revealed population genetic structure within P. carbonelli in the contact zone; geographical and genomic clines analysis suggested strong selection against gene flow, but a relatively small proportion of the loci can introgress, mostly within the narrow contact zone. However, geographical clines revealed that a few introgressed loci show signs of potential positive selection, particularly into P. bocagei. Geographical clines also detected a signal of hybrid zone movement towards P. bocagei distribution. Genomic cline analysis revealed heterogeneous patterns of introgression among loci within the syntopy zone, but the majority maintain a strong association with the genomic background of origin. However, incongruences between both cline approaches were found, potentially driven by confounding effects on genomic clines. Last, an important role of the Z chromosome in reproductive isolation is suggested. Importantly, overall patterns of restricted introgression seem to result from numerous strong intrinsic barriers across the genome.

Inferring genetic structure when there is little: population genetics versus genomics of the threatened bat Miniopterus schreibersii across Europe.

Despite their paramount importance in molecular ecology and conservation, genetic diversity and structure remain challenging to quantify with traditional genotyping methods. Next-generation sequencing holds great promises, but this has not been properly tested in highly mobile species. In this article, we compared microsatellite and RAD-sequencing (RAD-seq) analyses to investigate population structure in the declining bent-winged bat (Miniopterus schreibersii) across Europe. Both markers retrieved general patterns of weak range-wide differentiation, little sex-biased dispersal, and strong isolation by distance that associated with significant genetic structure between the three Mediterranean Peninsulas, which could have acted as glacial refugia. Microsatellites proved uninformative in individual-based analyses, but the resolution offered by genomic SNPs illuminated on regional substructures within several countries, with colonies sharing migrators of distinct ancestry without admixture. This finding is consistent with a marked philopatry and spatial partitioning between mating and rearing grounds in the species, which was suspected from marked-recaptured data. Our study advocates that genomic data are necessary to properly unveil the genetic footprints left by biogeographic processes and social organization in long-distant flyers, which are otherwise rapidly blurred by their high levels of gene flow.

Species recognition limits mating between hybridizing ant species.

Identifying mechanisms limiting hybridization is a central goal of speciation research. Here, we studied premating and postmating barriers to hybridization between two ant species, Formica selysi and Formica cinerea. These species hybridize in the Rhône valley in Switzerland, where they form a mosaic hybrid zone, with limited introgression from F. selysi into F. cinerea. There was no sign of temporal isolation between the two species in the production of queens and males. With choice experiments, we showed that queens and males strongly prefer to mate with conspecifics. Yet, we did not detect postmating barriers caused by genetic incompatibilities. Specifically, hybrids of all sexes and castes were found in the field and F1 hybrid workers did not show reduced viability compared to nonhybrid workers. To gain insights into the cues involved in species recognition, we analyzed the cuticular hydrocarbons (CHCs) of queens, males, and workers and staged dyadic encounters between workers. CHC profiles differed markedly between species, but were similar in F. cinerea and hybrids. Accordingly, workers also discriminated species, but they did not discriminate F. cinerea and hybrids. We discuss how the CHC-based recognition system of ants may facilitate the establishment of premating barriers to hybridization, independent of hybridization costs.

A socially polymorphic Formica ant species exhibits a novel distribution of social supergene genotypes.

Most supergenes discovered so far are young, occurring in one species or a few closely related species. An ancient supergene in the ant genus Formica presents an unusual opportunity to compare supergene-associated phenotypes and the factors that influence the persistence of polymorphism in different species. We investigate the genetic architecture of social organization in Formica francoeuri, an ant species native to low- and mid-elevation semiarid regions of southern California, and describe an efficient technique for estimating mode of social organization using population genomic data. Using this technique, we show that F. francoeuri exhibits polymorphism in colony social organization and that the phenotypic polymorphism is strongly associated with genotypes within the Formica social supergene region. The distribution of supergene haplotypes in F. francoeuri differs from that of related species Formica selysi in that colonies with multiple queens contain almost exclusively workers that are heterozygous for alternative supergene haplotypes. Moreover, heterozygous workers exhibit allele-specific expression of the polygyne-associated haplotype at the candidate gene Knockout, which is thought to influence social organization. We also report geographic population structure and variation in worker size across a large fraction of the species range. Our results suggest that, although the Formica supergene is conserved within the genus, the mechanisms that maintain the supergene and its associated polymorphisms differ among species.

Linked supergenes underlie split sex ratio and social organization in an ant.

Sexually reproducing organisms usually invest equally in male and female offspring. Deviations from this pattern have led researchers to new discoveries in the study of parent-offspring conflict, genomic conflict, and cooperative breeding. Some social insect species exhibit the unusual population-level pattern of split sex ratio, wherein some colonies specialize in the production of future queens and others specialize in the production of males. Theoretical work predicted that worker control of sex ratio and variation in relatedness asymmetry among colonies would cause each colony to specialize in the production of one sex. While some empirical tests supported theoretical predictions, others deviated from them, leaving many questions about how split sex ratio emerges. One factor yet to be investigated is whether colony sex ratio may be influenced by the genotypes of queens or workers. Here, we sequence the genomes of 138 Formica glacialis workers from 34 male-producing and 34 gyne-producing colonies to determine whether split sex ratio is under genetic control. We identify a supergene spanning 5.5 Mbp that is closely associated with sex allocation in this system. Strikingly, this supergene is adjacent to another supergene spanning 5 Mbp that is associated with variation in colony queen number. We identify a similar pattern in a second related species, Formica podzolica. The discovery that split sex ratio is determined, at least in part, by a supergene in two species opens future research on the evolutionary drivers of split sex ratio.

The maintenance of polymorphism in an ancient social supergene.

Supergenes, regions of the genome with suppressed recombination between sets of functional mutations, contribute to the evolution of complex phenotypes in diverse systems. Excluding sex chromosomes, most supergenes discovered so far appear to be young, being found in one species or a few closely related species. Here, we investigate how a chromosome harbouring an ancient supergene has evolved over about 30 million years (Ma). The Formica supergene underlies variation in colony queen number in at least five species. We expand previous analyses of sequence divergence on this chromosome to encompass about 90 species spanning the Formica phylogeny. Within the nonrecombining region, the gene knockout contains 22 single nucleotide polymorphisms (SNPs) that are consistently differentiated between two alternative supergene haplotypes in divergent European Formica species, and we show that these same SNPs are present in most Formica clades. In these clades, including an early diverging Nearctic Formica clade, individuals with alternative genotypes at knockout also have higher differentiation in other portions of this chromosome. We identify hotspots of SNPs along this chromosome that are present in multiple Formica clades to detect genes that may have contributed to the emergence and maintenance of the genetic polymorphism. Finally, we infer three gene duplications on one haplotype, based on apparent heterozygosity within these genes in the genomes of haploid males. This study strengthens the evidence that this supergene originated early in the evolution of Formica and that just a few loci in this large region of suppressed recombination retain strongly differentiated alleles across contemporary Formica lineages.

Mass of genes rather than master genes underlie the genomic architecture of amphibian speciation.

The genetic architecture of speciation, i.e., how intrinsic genomic incompatibilities promote reproductive isolation (RI) between diverging lineages, is one of the best-kept secrets of evolution. To directly assess whether incompatibilities arise in a limited set of large-effect speciation genes, or in a multitude of loci, we examined the geographic and genomic landscapes of introgression across the hybrid zones of 41 pairs of frog and toad lineages in the Western Palearctic region. As the divergence between lineages increases, phylogeographic transitions progressively become narrower, and larger parts of the genome resist introgression. This suggests that anuran speciation proceeds through a gradual accumulation of multiple barrier loci scattered across the genome, which ultimately deplete hybrid fitness by intrinsic postzygotic isolation, with behavioral isolation being achieved only at later stages. Moreover, these loci were disproportionately sex linked in one group (Hyla) but not in others (Rana and Bufotes), implying that large X-effects are not necessarily a rule of speciation with undifferentiated sex chromosomes. The highly polygenic nature of RI and the lack of hemizygous X/Z chromosomes could explain why the speciation clock ticks slower in amphibians compared to other vertebrates. The clock-like dynamics of speciation combined with the analytical focus on hybrid zones offer perspectives for more standardized practices of species delimitation.

Nuclear phylogenies and genomics of a contact zone establish the species rank of Podarcis lusitanicus (Squamata, Lacertidae).

Unravelling when divergent lineages constitute distinct species can be challenging, particularly in complex scenarios combining cryptic diversity and phylogenetic discordances between different types of molecular markers. Combining a phylogenetic approach with the study of contact zones can help to overcome such difficulties. The Podarcis hispanicus species complex has proven to be prosperous in independent evolutionary units, sometimes associated with cryptic diversity. Previous studies have revealed that one of the species of this complex, P. guadarramae, comprises two deeply divergent yet morphologically indistinguishable evolutionary units, currently regarded as subspecies (P. g. guadarramae and P. g. lusitanicus). In this study we used molecular data to address the systematics of the two lineages of Podarcis guadarramae and the closely related P. bocagei. Firstly, we reconstructed the species tree of these three and two additional taxa based on 30 nuclear loci using the multispecies coalescent with and without gene flow. Secondly, we used SNPs obtained from RADseq data to analyze the population structure across the distribution limits P. g. lusitanicus and P. g. guadarramae, and for comparison, a contact zone between P. bocagei and P. g. lusitanicus. Nuclear phylogenetic relationships between these three taxa are clearly difficult to determine due to the influence of gene flow, but our results give little support to the monophyly of P. guadarramae, potentially due to a nearly simultaneous divergence between them. Genetic structure and geographic cline analysis revealed that the two lineages of P. guadarramae replace each other abruptly across the sampled region and that gene flow is geographically restricted, implying the existence of strong reproductive isolation. Podarcis bocagei and P. g. lusitanicus show a similar degree of genetic differentiation and reproductive isolation, with very low levels of admixture in syntopy. These results support that all three forms are equally differentiated and reproductively isolated. In consequence, we conclude that the two former subspecies of Podarcis guadarramae constitute valid, yet cryptic species, that should be referred to as P. lusitanicus and P. guadarramae.

Comparing divergence landscapes from reduced-representation and whole genome resequencing in the yellow-rumped warbler (Setophaga coronata) species complex.

Researchers seeking to generate genomic data for non-model organisms are faced with a number of trade-offs when deciding which method to use. The selection of reduced representation approaches versus whole genome resequencing will ultimately affect the marker density, sequencing depth, and the number of individuals that can multiplexed. These factors can affect researchers' ability to accurately characterize certain genomic features, such as landscapes of divergence-how FST varies across the genomes. To provide insight into the effect of sequencing method on the estimation of divergence landscapes, we applied an identical bioinformatic pipeline to three generations of sequencing data (GBS, ddRAD, and WGS) produced for the same system, the yellow-rumped warbler species complex. We compare divergence landscapes generated using each method for the myrtle warbler (Setophaga coronata coronata) and the Audubon's warbler (S. c. auduboni), and for Audubon's warblers with deeply divergent mtDNA resulting from mitochondrial introgression. We found that most high-FST peaks were not detected in the ddRAD data set, and that while both GBS and WGS were able to identify the presence of large peaks, WGS was superior at a finer scale. Comparing Audubon's warblers with divergent mitochondrial haplotypes, only WGS allowed us to identify small (10-20 kb) regions of elevated differentiation, one of which contained the nuclear-encoded mitochondrial gene NDUFAF3. We calculated the cost per base pair for each method and found it was comparable between GBS and WGS, but significantly higher for ddRAD. These comparisons highlight the advantages of WGS over reduced representation methods when characterizing landscapes of divergence.

The population genetics of nonmigratory Allen's Hummingbird (Selasphorus sasin sedentarius) following a recent mainland colonization.

Allen's Hummingbird comprises two subspecies, one migratory (Selasphorus sasin sasin) and one nonmigratory (S. s. sedentarius). The nonmigratory subspecies, previously endemic to the California Channel Islands, apparently colonized the California mainland on the Palos Verdes Peninsula some time before 1970 and now breeds throughout coastal southern California. We sequenced and compared populations of mainland nonmigratory Allen's Hummingbird to Channel Island populations from Santa Catalina, San Clemente, and Santa Cruz Island. We found no evidence of founder effects on the mainland population. Values of nucleotide diversity on the mainland were higher than on the Channel Islands. There were low levels of divergence between the Channel Islands and the mainland, and Santa Cruz Island was the most genetically distinct. Ecological niche models showed that rainfall and temperature variables on the Channel Islands are similar in the Los Angeles basin and predicted continued expansion of nonmigratory Allen's Hummingbird north along the coast and inland. We also reviewed previous genetic studies of vertebrate species found on the Channel Islands and mainland and showed that broad conclusions regarding island-mainland patterns remain elusive. Challenges include the idiosyncratic nature of colonization itself as well as the lack of a comprehensive approach that incorporates similar markers and sampling strategies across taxa, which, within the context of a comparative study of island-mainland relationships, may lead to inconsistent results.

Variable levels of introgression between the endangered Podarcis carbonelli and highly divergent congeneric species.

Recent empirical studies have demonstrated that speciation with gene flow is more common than previously thought. From a conservation perspective, the potential negative effects of hybridization raise concerns on the genetic integrity of endangered species. However, introgressive hybridization has also been growingly recognized as a source of diversity and new advantageous alleles. Carbonell's wall lizard (Podarcis carbonelli) is an endangered species whose distribution overlaps with four other congeneric species. Our goal here was to determine whether P. carbonelli is completely reproductively isolated from its congeners and to evaluate the relevance of hybridization and interspecific gene flow for developing a conservation plan. We used restriction site associated DNA (RAD) sequencing to discover SNPs in samples from four contact zones between P. carbonelli and four other species. Principal component analysis, multilocus genotype assignment and interspecific heterozygosity suggest incomplete reproductive isolation and ongoing gene flow between species. However, hybridization dynamics vary across all pairs, suggesting complex interactions between multiple intrinsic and extrinsic barriers. Despite seemingly ubiquitous interspecific gene flow, we found evidence of strong reproductive isolation across most contact zones. Instead, indirect effects of hybridization like waste of reproductive effort in small isolated populations may be more problematic. Our results highlight the need to further evaluate the consequences of introgression for P. carbonelli, both on a geographic and genomic level and included in a comprehensive and urgently needed conservation plan. Besides, those findings will add important insights on the potential effects of hybridization and introgression for endangered species.

When Sex Chromosomes Recombine Only in the Heterogametic Sex: Heterochiasmy and Heterogamety in Hyla Tree Frogs.

Sex chromosomes are classically predicted to stop recombining in the heterogametic sex, thereby enforcing linkage between sex-determining (SD) and sex-antagonistic (SA) genes. With the same rationale, a pre-existing sex asymmetry in recombination is expected to affect the evolution of heterogamety, for example, a low rate of male recombination might favor transitions to XY systems, by generating immediate linkage between SD and SA genes. Furthermore, the accumulation of deleterious mutations on nonrecombining Y chromosomes should favor XY-to-XY transitions (which discard the decayed Y), but disfavor XY-to-ZW transitions (which fix the decayed Y as an autosome). Like many anuran amphibians, Hyla tree frogs have been shown to display drastic heterochiasmy (males only recombine at chromosome tips) and are typically XY, which seems to fit the above expectations. Instead, here we demonstrate that two species, H. sarda and H. savignyi, share a common ZW system since at least 11 Ma. Surprisingly, the typical pattern of restricted male recombination has been maintained since then, despite female heterogamety. Hence, sex chromosomes recombine freely in ZW females, not in ZZ males. This suggests that heterochiasmy does not constrain heterogamety (and vice versa), and that the role of SA genes in the evolution of sex chromosomes might have been overemphasized.

Pigmentation Genes Show Evidence of Repeated Divergence and Multiple Bouts of Introgression in Setophaga Warblers.

Species radiations have long served as model systems in evolutionary biology.1,2 However, it has only recently become possible to study the genetic bases of the traits responsible for diversification and only in a small number of model systems.3 Here, we use genomes of 36 species of North, Central, and South American warblers to highlight the role of pigmentation genes-involved in melanin and carotenoid processing-in the diversification of this group. We show that agouti signaling protein (ASIP) and beta-carotene oxygenase 2 (BCO2) are predictably divergent between species that differ in the distribution of melanin and carotenoid in their plumages, respectively. Among species, sequence variation at ASIP broadly mirrors the species' phylogenetic history, consistent with repeated, independent mutations generating melanin-based variation. In contrast, BCO2 variation is highly discordant from the species tree, with evidence of cross-lineage introgression among species like the yellow warbler (Setophaga petechia) and magnolia warbler (S. magnolia) with extensive carotenoid-based coloration. We also detect introgression of a small part of the BCO2 coding region (<3 kb) in S. discolor and S. vitellina, including an amino acid substitution that is unique to warblers but otherwise highly conserved across birds. Lateral transfer of carotenoid-processing genes has been documented in arthropods, but introgression of BCO2 as demonstrated here-presumably adaptive-represents the first example of carotenoid gene transfer among vertebrates. These contrasting genomic patterns show that both independent evolution in a common set of genes and past hybridization have fueled plumage diversification in this colorful avian radiation.

CYP2J19 mediates carotenoid colour introgression across a natural avian hybrid zone.

It has long been of interest to identify the phenotypic traits that mediate reproductive isolation between related species, and more recently, the genes that underpin them. Much work has focused on identifying genes associated with animal colour, with the candidate gene CYP2J19 identified in laboratory studies as the ketolase converting yellow dietary carotenoids to red ketocarotenoids in birds with red pigments. However, evidence that CYP2J19 explains variation between red and yellow feather coloration in wild populations of birds is lacking. Hybrid zones provide the opportunity to identify genes associated with specific traits. Here we investigate genomic regions associated with colour in red-fronted and yellow-fronted tinkerbirds across a hybrid zone in southern Africa. We sampled 85 individuals, measuring spectral reflectance of forecrown feathers and scoring colours from photographs, while testing for carotenoid presence with Raman spectroscopy. We performed a genome-wide association study to identify associations with carotenoid-based coloration, using double-digest RAD sequencing aligned to a short-read whole genome of a Pogoniulus tinkerbird. Admixture mapping using 104,933 single nucleotide polymorphisms (SNPs) identified a region of chromosome 8 that includes CYP2J19 as the only locus with more than two SNPs significantly associated with both crown hue and crown score, while Raman spectra provided evidence of ketocarotenoids in red feathers. Asymmetric backcrossing in the hybrid zone suggests that yellow-fronted females mate more often with red-fronted males than vice versa. Female red-fronted tinkerbirds mating assortatively with red-crowned males is consistent with the hypothesis that converted carotenoids are an honest signal of quality.

Genomic differentiation across the speciation continuum in three hummingbird species pairs.

BACKGROUND: The study of speciation has expanded with the increasing availability and affordability of high-resolution genomic data. How the genome evolves throughout the process of divergence and which regions of the genome are responsible for causing and maintaining that divergence have been central questions in recent work. Here, we use three pairs of species from the recently diverged bee hummingbird clade to investigate differences in the genome at different stages of speciation, using divergence times as a proxy for the speciation continuum. RESULTS: Population measures of relative differentiation between hybridizing species reveal that different chromosome types diverge at different stages of speciation. Using FST as our relative measure of differentiation we found that the sex chromosome shows signs of divergence early in speciation. Next, small autosomes (microchromosomes) accumulate highly diverged genomic regions, while the large autosomes (macrochromosomes) accumulate genomic regions of divergence at a later stage of speciation. CONCLUSIONS: Our finding that genomic windows of elevated FST accumulate on small autosomes earlier in speciation than on larger autosomes is counter to the prediction that FST increases with size of chromosome (i.e. with decreased recombination rate), and is not represented when weighted average FST per chromosome is compared with chromosome size. The results of this study suggest that multiple chromosome characteristics such as recombination rate and gene density combine to influence the genomic locations of signatures of divergence.

Subspecies differentiation in an enigmatic chaparral shrub species.

PREMISE: Delimiting biodiversity units is difficult in organisms in which differentiation is obscured by hybridization, plasticity, and other factors that blur phenotypic boundaries. Such work is more complicated when the focal units are subspecies, the definition of which has not been broadly explored in the era of modern genetic methods. Eastwood manzanita (Arctostaphylos glandulosa Eastw.) is a widely distributed and morphologically complex chaparral shrub species with much subspecific variation, which has proven challenging to categorize. Currently 10 subspecies are recognized, however, many of them are not geographically segregated, and morphological intermediates are common. Subspecies delimitation is of particular importance in this species because two of the subspecies are rare. The goal of this study was to apply an evolutionary definition of "subspecies" to characterize structure within Eastwood manzanita. METHODS: We used publicly available geospatial environmental data and reduced-representation genome sequencing to characterize environmental and genetic differentiation among subspecies. In addition, we tested whether subspecies could be differentiated by environmentally associated genetic variation. RESULTS: Our analyses do not show genetic differentiation among subspecies of Eastwood manzanita, with the exception of one of the two rare subspecies. In addition, our environmental analyses did not show ecological differentiation, though limitations of the analysis prevent strong conclusions. CONCLUSIONS: Genetic structure within Eastwood manzanita does not correspond to current subspecies circumscriptions, but rather reflects geographic distribution. Our study suggests that subspecies concepts need to be reconsidered in long-lived plant species, especially in the age of next-generation sequencing.

An Ancient and Eroded Social Supergene Is Widespread across Formica Ants.

Supergenes, clusters of tightly linked genes, play a key role in the evolution of complex adaptive variation [1, 2]. Although supergenes have been identified in many species, we lack an understanding of their origin, evolution, and persistence [3]. Here, we uncover 20-40 Ma of evolutionary history of a supergene associated with polymorphic social organization in Formica ants [4]. We show that five Formica species exhibit homologous divergent haplotypes spanning 11 Mbp on chromosome 3. Despite the supergene's size, only 142 single nucleotide polymorphisms (SNPs) consistently distinguish alternative supergene haplotypes across all five species. These conserved trans-species SNPs are localized in a small number of disjunct clusters distributed across the supergene. This unexpected pattern of divergence indicates that the Formica supergene does not follow standard models of sex chromosome evolution, in which distinct evolutionary strata reflect an expanding region of suppressed recombination [5]. We propose an alternative "eroded strata model" in which clusters of conserved trans-species SNPs represent functionally important areas maintained by selection in the face of rare recombination between ancestral haplotypes. The comparison of whole-genome sequences across 10 additional Formica species reveals that the most conserved region of the supergene contains a transcription factor essential for motor neuron development in Drosophila [6]. The discovery that a very small portion of this large and ancient supergene harbors conserved trans-species SNPs linked to colony social organization suggests that the ancestral haplotypes have been eroded by recombination, with selection preserving differentiation at one or a few genes generating alternative social organization.