Runs of homozygosity reveal contrasting histories of inbreeding across global lineages of the edible porcini mushroom, Boletus edulis.

Inbreeding, the mating of individuals that are related through common ancestry, is of central importance in evolutionary and conservation biology due to its impacts on individual fitness and population dynamics. However, while advanced genomic approaches have revolutionised the study of inbreeding in animals, genomic studies of inbreeding are rare in plants and lacking in fungi. We investigated global patterns of inbreeding in the prized edible porcini mushroom Boletus edulis using 225 whole genomes from seven lineages distributed across the northern hemisphere. Genomic inbreeding was quantified using runs of homozygosity (ROHs). We found appreciable variation both among and within lineages, with some individuals having over 20% of their genomes in ROHs. Much of this variation could be explained by a combination of elevation and latitude, and to a lesser extent by predicted habitat suitability during the last glacial maximum. In line with this, the majority of ROHs were short, reflecting ancient common ancestry dating back approximately 200-1700 generations ago, while longer ROHs indicative of recent common ancestry (less than approximately 50 generations ago) were infrequent. Our study reveals the inbreeding legacy of major climatic events in a widely distributed forest mutualist, aligning with prevailing theories and empirical studies of the impacts of historical glaciation events on the dominant forest tree species of the northern hemisphere.

Hybrid speciation driven by multilocus introgression of ecological traits.

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

Complexity of Chemical Emissions Increases Concurrently with Sexual Maturity in Heliconius Butterflies.

Pheromone communication is widespread among animals. Since it is often involved in mate choice, pheromone production is often tightly controlled. Although male sex pheromones (MSPs) and anti-aphrodisiacs have been studied in some Heliconius butterfly species, little is known about the factors affecting their production and release in these long-lived butterflies. Here, we investigate the effect of post-eclosion age on chemical blends from pheromone-emitting tissues in Heliconius atthis and Heliconius charithonia, exhibiting respectively free-mating and pupal-mating strategies that are hypothesised to differently affect the timing of their pheromone emissions. We focus on two different tissues: the wing androconia, responsible for MSPs used in courtship, and the genital tip, the production site for anti-aphrodisiac pheromones that affect post-mating behaviour. Gas chromatography-mass spectrometric analysis of tissue extracts from virgin males and females of both species from day 0 to 8 post-eclosion demonstrates the following. Some ubiquitous fatty acid precursors are already detectable at day 0. The complexity of the chemical blends increases with age regardless of tissue or sex. No obvious difference in the time course of blend production was evident between the two species, but female tissues in H. charithonia were more affected by age than in H. atthis. We suggest that compounds unique to male androconia and genitals and whose amount increases with age are potential candidates for future investigation into their roles as pheromones. While this analysis revealed some of the complexity in Heliconius chemical ecology, the effects of other factors, such as the time of day, remain unknown.

Pervasive mimicry in flight behavior among aposematic butterflies.

Flight was a key innovation in the adaptive radiation of insects. However, it is a complex trait influenced by a large number of interacting biotic and abiotic factors, making it difficult to unravel the evolutionary drivers. We investigate flight patterns in neotropical heliconiine butterflies, well known for mimicry of their aposematic wing color patterns. We quantify the flight patterns (wing beat frequency and wing angles) of 351 individuals representing 29 heliconiine and 9 ithomiine species belonging to ten color pattern mimicry groupings. For wing beat frequency and up wing angles, we show that heliconiine species group by color pattern mimicry affiliation. Convergence of down wing angles to mimicry groupings is less pronounced, indicating that distinct components of flight are under different selection pressures and constraints. The flight characteristics of the Tiger mimicry group are particularly divergent due to convergence with distantly related ithomiine species. Predator-driven selection for mimicry also explained variation in flight among subspecies, indicating that this convergence can occur over relatively short evolutionary timescales. Our results suggest that the flight convergence is driven by aposematic signaling rather than shared habitat between comimics. We demonstrate that behavioral mimicry can occur between lineages that have separated over evolutionary timescales ranging from <0.5 to 70 My.

The Presence of Two MyoD Genes in a Subset of Acanthopterygii Fish Is Associated with a Polyserine Insert in MyoD1.

The MyoD gene was duplicated during the teleost whole genome duplication and, while a second MyoD gene (MyoD2) was subsequently lost from the genomes of some lineages (including zebrafish), many fish lineages (including Alcolapia species) have retained both MyoD paralogues. Here we reveal the expression patterns of the two MyoD genes in Oreochromis (Alcolapia) alcalica using in situ hybridisation. We report our analysis of MyoD1 and MyoD2 protein sequences from 54 teleost species, and show that O. alcalica, along with some other teleosts, include a polyserine repeat between the amino terminal transactivation domains (TAD) and the cysteine-histidine rich region (H/C) in MyoD1. The evolutionary history of MyoD1 and MyoD2 is compared to the presence of this polyserine region using phylogenetics, and its functional relevance is tested using overexpression in a heterologous system to investigate subcellular localisation, stability, and activity of MyoD proteins that include and do not include the polyserine region.

Two closely related ureotelic fish species of the genus Alcolapia express different levels of ammonium transporters in gills.

Most fish excrete their nitrogenous waste across the gills as ammonia through the activity of the Rhesus glycoprotein ammonium transporters. In contrast, fish of the subgenus Alcolapia (Oreochromis) are the only vertebrates that survive the extreme conditions of the soda lakes of Natron and Magadi in East Africa and have evolved adaptations to the highly alkaline waters including the ability to excrete their nitrogenous waste as urea. Nevertheless, Alcolapia retain the Rhesus glycoprotein genes in their genomes and using two heterologous expression systems, we demonstrate that Alcolapia Rhbg is capable of moving ammonia. Comparing ammonia and urea excretion from two closely related Alcolapia species from the same aquarium, we found that while Alcolapia grahami remains fully ureotelic after many generations in lab conditions, Alcolapia alcalica excretes some of its nitrogenous waste as ammonia. Using in situ hybridisation, we demonstrate robust, localised gene expression of Rhbg, rhcg1 and rhcg2 in the gill tissue in both A. alcalica embryos and adults, similar to that in other ammoniotelic fish. In contrast, the expression of these genes in A. grahami gills is much lower than in A. alcalica, suggesting the rapid evolution of a molecular mechanism underlying the complete ureotelism of A. grahami.

Exploitation of an ancestral pheromone biosynthetic pathway contributes to diversification in Heliconius butterflies.

During courtship, male butterflies of many species produce androconial secretions containing male sex pheromones (MSPs) that communicate species identity and affect female choice. MSPs are thus likely candidates as reproductive barriers, yet their role in speciation remains poorly studied. Although Heliconius butterflies are a model system in speciation, their MSPs have not been investigated from a macroevolutionary perspective. We use GC/MS to characterize male androconial secretions in 33 of the 69 species in the Heliconiini tribe. We found these blends to be species-specific, consistent with a role in reproductive isolation. We detected a burst in blend diversification rate at the most speciose genus, Heliconius; a consequence of Heliconius and Eueides species using a fatty acid (FA) metabolic pathway to unlock more complex blends than basal Heliconiini species, whose secretions are dominated by plant-like metabolites. A comparison of 10 sister species pairs demonstrates a striking positive correlation between blend dissimilarity and range overlap, consistent with character displacement or reinforcement in sympatry. These results demonstrate for the first time that MSP diversification can promote reproductive isolation across this group of butterflies, showcasing how implementation of an ancestral trait, the co-option of the FA metabolic pathway for pheromone production, can facilitate rapid speciation.

Demographic Reconstruction of Antarctic Fur Seals Supports the Krill Surplus Hypothesis.

Much debate surrounds the importance of top-down and bottom-up effects in the Southern Ocean, where the harvesting of over two million whales in the mid twentieth century is thought to have produced a massive surplus of Antarctic krill. This excess of krill may have allowed populations of other predators, such as seals and penguins, to increase, a top-down hypothesis known as the 'krill surplus hypothesis'. However, a lack of pre-whaling population baselines has made it challenging to investigate historical changes in the abundance of the major krill predators in relation to whaling. Therefore, we used reduced representation sequencing and a coalescent-based maximum composite likelihood approach to reconstruct the recent demographic history of the Antarctic fur seal, a pinniped that was hunted to the brink of extinction by 18th and 19th century sealers. In line with the known history of this species, we found support for a demographic model that included a substantial reduction in population size around the time period of sealing. Furthermore, maximum likelihood estimates from this model suggest that the recovered, post-sealing population at South Georgia may have been around two times larger than the pre-sealing population. Our findings lend support to the krill surplus hypothesis and illustrate the potential of genomic approaches to shed light on long-standing questions in population biology.

Complex basis of hybrid female sterility and Haldane's rule in Heliconius butterflies: Z-linkage and epistasis.

Hybrids between species are often sterile or inviable. Hybrid unfitness usually evolves first in the heterogametic sex-a pattern known as Haldane's rule. The genetics of Haldane's rule have been extensively studied in species where the male is the heterogametic (XX/XY) sex, but its basis in taxa where the female is heterogametic (ZW/ZZ), such as Lepidoptera and birds, is largely unknown. Here, we analyse a new case of female hybrid sterility between geographic subspecies of Heliconius pardalinus. The two subspecies mate freely in captivity, but female F1 hybrids in both directions of cross are sterile. Sterility is due to arrested development of oocytes after they become differentiated from nurse cells, but before yolk deposition. We backcrossed fertile male F1 hybrids to parental females and mapped quantitative trait loci (QTLs) for female sterility. We also identified genes differentially expressed in the ovary as a function of oocyte development. The Z chromosome has a major effect, similar to the 'large X effect' in Drosophila, with strong epistatic interactions between loci at either end of the Z chromosome, and between the Z chromosome and autosomal loci on chromosomes 8 and 20. By intersecting the list of genes within these QTLs with those differentially expressed in sterile and fertile hybrids, we identified three candidate genes with relevant phenotypes. This study is the first to characterize hybrid sterility using genome mapping in the Lepidoptera and shows that it is produced by multiple complex epistatic interactions often involving the sex chromosome, as predicted by the dominance theory of Haldane's rule.

Adaptation of the carbamoyl-phosphate synthetase enzyme in an extremophile fish.

Tetrapods and fish have adapted distinct carbamoyl-phosphate synthase (CPS) enzymes to initiate the ornithine urea cycle during the detoxification of nitrogenous wastes. We report evidence that in the ureotelic subgenus of extremophile fish Oreochromis Alcolapia, CPS III has undergone convergent evolution and adapted its substrate affinity to ammonia, which is typical of terrestrial vertebrate CPS I. Unusually, unlike in other vertebrates, the expression of CPS III in Alcolapia is localized to the skeletal muscle and is activated in the myogenic lineage during early embryonic development with expression remaining in mature fish. We propose that adaptation in Alcolapia included both convergent evolution of CPS function to that of terrestrial vertebrates, as well as changes in development mechanisms redirecting CPS III gene expression to the skeletal muscle.

Past, current, and potential future distributions of unique genetic diversity in a cold-adapted mountain butterfly.

AIM: Climatic changes throughout the Pleistocene have strongly modified species distributions. We examine how these range shifts have affected the genetic diversity of a montane butterfly species and whether the genetic diversity in the extant populations is threatened by future climate change. LOCATION: Europe. TAXON: Erebia epiphron Lepidoptera: Nymphalidae. METHODS: We analyzed mtDNA to map current genetic diversity and differentiation of E. epiphron across Europe to identify population refugia and postglacial range shifts. We used species distribution modeling (SDM) to hindcast distributions over the last 21,000 years to identify source locations of extant populations and to project distributions into the future (2070) to predict potential losses in genetic diversity. RESULTS: We found substantial genetic diversity unique to specific regions within Europe (total number of haplotypes = 31, number of unique haplotypes = 27, H d = 0.9). Genetic data and SDM hindcasting suggest long-term separation and survival of discrete populations. Particularly, high rates of unique diversity in postglacially colonized sites in England (H d = 0.64) suggest this population was colonized from a now extinct cryptic refugium. Under future climate change, SDMs predict loss of climate suitability for E. epiphron, particularly at lower elevations (<1,000 meters above sea level) equating to 1 to 12 unique haplotypes being at risk under climate scenarios projecting 1°C and 2-3°C increases respectfully in global temperature by 2070. MAIN CONCLUSIONS: Our results suggest that historical range expansion and retraction processes by a cold-adapted mountain species caused diversification between populations, resulting in unique genetic diversity which may be at risk if distributions of cold-adapted species shrink in future. Assisted colonizations of individuals from at-risk populations into climatically suitable unoccupied habitat might help conserve unique genetic diversity, and translocations into remaining populations might increase their genetic diversity and hence their ability to adapt to future climate change.

Exploring the Expression of Cardiac Regulators in a Vertebrate Extremophile: The Cichlid Fish Oreochromis (Alcolapia) alcalica.

Although it is widely accepted that the cellular and molecular mechanisms of vertebrate cardiac development are evolutionarily conserved, this is on the basis of data from only a few model organisms suited to laboratory studies. Here, we investigate gene expression during cardiac development in the extremophile, non-model fish species, Oreochromis (Alcolapia) alcalica. We first characterise the early development of O. alcalica and observe extensive vascularisation across the yolk prior to hatching. We further investigate heart development by identifying and cloning O. alcalica orthologues of conserved cardiac transcription factors gata4, tbx5, and mef2c for analysis by in situ hybridisation. Expression of these three key cardiac developmental regulators also reveals other aspects of O. alcalica development, as these genes are expressed in developing blood, limb, eyes, and muscle, as well as the heart. Our data support the notion that O. alcalica is a direct-developing vertebrate that shares the highly conserved molecular regulation of the vertebrate body plan. However, the expression of gata4 in O. alcalica reveals interesting differences in the development of the circulatory system distinct from that of the well-studied zebrafish. Understanding the development of O. alcalica embryos is an important step towards providing a model for future research into the adaptation to extreme conditions; this is particularly relevant given that anthropogenic-driven climate change will likely result in more freshwater organisms being exposed to less favourable conditions.

Deep Convergence, Shared Ancestry, and Evolutionary Novelty in the Genetic Architecture of Heliconius Mimicry.

Convergent evolution can occur through different genetic mechanisms in different species. It is now clear that convergence at the genetic level is also widespread, and can be caused by either (i) parallel genetic evolution, where independently evolved convergent mutations arise in different populations or species, or (ii) collateral evolution in which shared ancestry results from either ancestral polymorphism or introgression among taxa. The adaptive radiation of Heliconius butterflies shows color pattern variation within species, as well as mimetic convergence between species. Using comparisons from across multiple hybrid zones, we use signals of shared ancestry to identify and refine multiple putative regulatory elements in Heliconius melpomene and its comimics, Heliconius elevatus and Heliconius besckei, around three known major color patterning genes: optix, WntA, and cortex While we find that convergence between H. melpomene and H. elevatus is caused by a complex history of collateral evolution via introgression in the Amazon, convergence between these species in the Guianas appears to have evolved independently. Thus, we find adaptive convergent genetic evolution to be a key driver of regulatory changes that lead to rapid phenotypic changes. Furthermore, we uncover evidence of parallel genetic evolution at some loci around optix and WntA in H. melpomene and its distant comimic Heliconius erato Ultimately, we show that all three of convergence, conservation, and novelty underlie the modular architecture of Heliconius color pattern mimicry.

Genomic architecture and introgression shape a butterfly radiation.

We used 20 de novo genome assemblies to probe the speciation history and architecture of gene flow in rapidly radiating Heliconius butterflies. Our tests to distinguish incomplete lineage sorting from introgression indicate that gene flow has obscured several ancient phylogenetic relationships in this group over large swathes of the genome. Introgressed loci are underrepresented in low-recombination and gene-rich regions, consistent with the purging of foreign alleles more tightly linked to incompatibility loci. Here, we identify a hitherto unknown inversion that traps a color pattern switch locus. We infer that this inversion was transferred between lineages by introgression and is convergent with a similar rearrangement in another part of the genus. These multiple de novo genome sequences enable improved understanding of the importance of introgression and selective processes in adaptive radiation.

Geographic contrasts between pre- and postzygotic barriers are consistent with reinforcement in Heliconius butterflies.

Identifying the traits causing reproductive isolation and the order in which they evolve is fundamental to understanding speciation. Here, we quantify prezygotic and intrinsic postzygotic isolation among allopatric, parapatric, and sympatric populations of the butterflies Heliconius elevatus and Heliconius pardalinus. Sympatric populations from the Amazon (H. elevatus and H. p. butleri) exhibit strong prezygotic isolation and rarely mate in captivity; however, hybrids are fertile. Allopatric populations from the Amazon (H. p. butleri) and Andes (H. p. sergestus) mate freely when brought together in captivity, but the female F1 hybrids are sterile. Parapatric populations (H. elevatus and H. p. sergestus) exhibit both assortative mating and sterility of female F1s. Assortative mating in sympatric populations is consistent with reinforcement in the face of gene flow, where the driving force, selection against hybrids, is due to disruption of mimicry and other ecological traits rather than hybrid sterility. In contrast, the lack of assortative mating and hybrid sterility observed in allopatric populations suggests that geographic isolation enables the evolution of intrinsic postzygotic reproductive isolation. Our results show how the types of reproductive barriers that evolve between species may depend on geography.

The genetic architecture of adaptation: convergence and pleiotropy in Heliconius wing pattern evolution.

Unravelling the genetic basis of adaptive traits is a major challenge in evolutionary biology. Doing so informs our understanding of evolution towards an adaptive optimum, the distribution of locus effect sizes, and the influence of genetic architecture on the evolvability of a trait. In the Müllerian co-mimics Heliconius melpomene and Heliconius erato some Mendelian loci affecting mimicry shifts are well known. However, several phenotypes in H. melpomene remain to be mapped, and the quantitative genetics of colour pattern variation has rarely been analysed. Here we use quantitative trait loci (QTL) analyses of crosses between H. melpomene races from Peru and Suriname to map, for the first time, the control of the broken band phenotype to WntA and identify a ~100 kb region controlling this variation. Additionally, we map variation in basal forewing red-orange pigmentation to a locus centred around the gene ventral veins lacking (vvl). The locus also appears to affect medial band shape variation as it was previously known to do in H. erato. This adds to the list of homologous regions controlling convergent phenotypes between these two species. Finally we show that Heliconius wing-patterning genes are strikingly pleiotropic among wing pattern traits. Our results demonstrate how genetic architecture can shape, aid and constrain adaptive evolution.

RAD Sequencing and a Hybrid Antarctic Fur Seal Genome Assembly Reveal Rapidly Decaying Linkage Disequilibrium, Global Population Structure and Evidence for Inbreeding.

Recent advances in high throughput sequencing have transformed the study of wild organisms by facilitating the generation of high quality genome assemblies and dense genetic marker datasets. These resources have the potential to significantly advance our understanding of diverse phenomena at the level of species, populations and individuals, ranging from patterns of synteny through rates of linkage disequilibrium (LD) decay and population structure to individual inbreeding. Consequently, we used PacBio sequencing to refine an existing Antarctic fur seal (Arctocephalus gazella) genome assembly and genotyped 83 individuals from six populations using restriction site associated DNA (RAD) sequencing. The resulting hybrid genome comprised 6,169 scaffolds with an N50 of 6.21 Mb and provided clear evidence for the conservation of large chromosomal segments between the fur seal and dog (Canis lupus familiaris). Focusing on the most extensively sampled population of South Georgia, we found that LD decayed rapidly, reaching the background level by around 400 kb, consistent with other vertebrates but at odds with the notion that fur seals experienced a strong historical bottleneck. We also found evidence for population structuring, with four main Antarctic island groups being resolved. Finally, appreciable variance in individual inbreeding could be detected, reflecting the strong polygyny and site fidelity of the species. Overall, our study contributes important resources for future genomic studies of fur seals and other pinnipeds while also providing a clear example of how high throughput sequencing can generate diverse biological insights at multiple levels of organization.

Supergene Evolution Triggered by the Introgression of a Chromosomal Inversion.

Supergenes are groups of tightly linked loci whose variation is inherited as a single Mendelian locus and are a common genetic architecture for complex traits under balancing selection [1-8]. Supergene alleles are long-range haplotypes with numerous mutations underlying distinct adaptive strategies, often maintained in linkage disequilibrium through the suppression of recombination by chromosomal rearrangements [1, 5, 7-9]. However, the mechanism governing the formation of supergenes is not well understood and poses the paradox of establishing divergent functional haplotypes in the face of recombination. Here, we show that the formation of the supergene alleles encoding mimicry polymorphism in the butterfly Heliconius numata is associated with the introgression of a divergent, inverted chromosomal segment. Haplotype divergence and linkage disequilibrium indicate that supergene alleles, each allowing precise wing-pattern resemblance to distinct butterfly models, originate from over a million years of independent chromosomal evolution in separate lineages. These "superalleles" have evolved from a chromosomal inversion captured by introgression and maintained in balanced polymorphism, triggering supergene inheritance. This mode of evolution involving the introgression of a chromosomal rearrangement is likely to be a common feature of complex structural polymorphisms associated with the coexistence of distinct adaptive syndromes. This shows that the reticulation of genealogies may have a powerful influence on the evolution of genetic architectures in nature.

Contrasting geographic structure in evolutionarily divergent Lake Tanganyika catfishes.

Geographic isolation is suggested to be among the most important processes in the generation of cichlid fish diversity in East Africa's Great Lakes, both through isolation by distance and fluctuating connectivity caused by changing lake levels. However, even broad scale phylogeographic patterns are currently unknown in many non-cichlid littoral taxa from these systems. To begin to address this, we generated restriction-site-associated DNA sequence (RADseq) data to investigate phylogeographic structure throughout Lake Tanganyika (LT) in two broadly sympatric rocky shore catfish species from independent evolutionary radiations with differing behaviors: the mouthbrooding claroteine, Lophiobagrus cyclurus, and the brood-parasite mochokid, Synodontis multipunctatus. Our results indicated contrasting patterns between these species, with strong lake-wide phylogeographic signal in L. cyclurus including a deep divergence between the northern and southern lake basins. Further structuring of these populations was observed across a heterogeneous habitat over much smaller distances. Strong population growth was observed in L. cyclurus sampled from shallow shorelines, suggesting population growth associated with the colonization of new habitats following lake-level rises. Conversely, S. multipunctatus, which occupies a broader depth range, showed little phylogeographic structure and lower rates of population growth. Our findings suggest that isolation by distance and/or habitat barriers may play a role in the divergence of non-cichlid fishes in LT, but this effect varies by species.

The Scent Chemistry of Heliconius Wing Androconia.

Neotropical Heliconius butterflies are members of various mimicry rings characterized by diverse colour patterns. In the present study we investigated whether a similar diversity is observed in the chemistry of volatile compounds present in male wing androconia. Recent research has shown that these androconia are used during courting of females. Three to five wild-caught male Heliconius individuals of 17 species and subspecies were analyzed by GC/MS. Most of the identified compounds originate from common fatty acids precursors, including aldehydes, alcohols, acetates or esters preferentially with a C18 and C20 chain, together with some alkanes. The compounds occurred in species-specific mixtures or signatures. For example, octadecanal is characteristic for H. melpomene, but variation in composition between the individuals was observed. Cluster analysis of compound occurrence in individual bouquets and analyses based on biosynthetic motifs such as functional group, chain length, or basic carbon-backbone modification were used to reveal structural patterns. Mimetic pairs contain different scent bouquets, but also some compounds in common, whereas sympatric species, both mimetic and non-mimetic, have more distinct compound compositions. The compounds identified here may play a role in mate choice thus helping maintain species integrity in a butterfly genus characterized by pervasive interspecific gene flow.