Genome-wide introgression among distantly related Heliconius butterfly species

Disentangling a genome-wide mosaic of conflicting phylogenetic signals in Western Rattlesnakes

Species tree inference is often assumed to be more accurate as datasets increase in size, with whole genomes representing the best-case-scenario for estimating a single, most-likely speciation history with high confidence. However, genomes may harbor a complex mixture of evolutionary histories among loci, which amplifies the opportunity for model misspecification and impacts phylogenetic inference. Accordingly, multiple distinct and well-supported phylogenetic trees are often recovered from genome-scale data, and approaches for biologically interpreting these distinct signatures are a major challenge for evolutionary biology in the age of genomics. Here, we analyze 32 whole genomes of nine taxa and two outgroups from the Western Rattlesnake species complex. Using concordance factors, topology weighting, and concatenated and species tree analyses with a chromosome-level reference genome, we characterize the distribution of phylogenetic signal across the genomic landscape. We find that concatenated and species tree analyses of autosomes, the Z (sex) chromosome, and mitochondrial genome yield distinct, yet strongly supported phylogenies. Analyses of site-specific likelihoods show additional patterns consistent with rampant model misspecification, a likely consequence of several evolutionary processes. Together, our results suggest that a combination of historic and recent introgression, along with natural selection, recombination rate variation, and cytonuclear co-evolution of nuclear-encoded mitochondrial genes, underlie genome-wide variation in phylogenetic signal. Our results highlight both the power and complexity of interpreting whole genomes in a phylogenetic context and illustrate how patterns of phylogenetic discordance can reveal the impacts of different evolutionary processes that contribute to genome-wide variation in phylogenetic signal.

Species delimitation in the Populus laurifolia complex (Salicaceae) based on phylogenetic and morphometric evidence

Due to significant morphological differences and extensive interspecific hybridization, there are numerous species complexes with taxonomic challenges in the genus Populus. Integrative taxonomy, which combines evidence of morphology, molecular phylogeny, niche differentiation, and reproductive isolation, provides the most effective approaches for species delimitation. The Populus laurifolia complex, which belongs to Populus subg. Tacamahaca (Salicaceae), is distributed in the Altai Mountains and Tianshan Mountains. This complex exhibits morphological variability, making species delimitation challenging. Due to limited sampling and systematic studies, its taxonomy has remained unresolved. In this study, 337 specimens, along with online digital samples representing nearly all wild populations, were collected. Morphological analyses were performed to evaluate key traits and clarify species boundaries. Phylogenetic relationships were reconstructed using concatenation and coalescent methods based on 566,375 nuclear single-nucleotide polymorphisms (SNPs). Ecological niche differentiation was assessed, and ABBA-BABA analysis was used to examine interspecific hybridization. The results revealed that this complex, based on a series of significant character states, could be morphologically distinguished into three species-P. laurifolia (Populus pilosa considered a synonym of P. laurifolia), Populus talassica, and Populus pamirica-which also correspond to three well-supported clades in the phylogenetic trees. P. pamirica exhibits some degree of ecological niche differentiation from P. talassica and P. laurifolia, whereas the latter two show minimal differentiation. Gene flow within the complex remains limited. This research underscores the importance of integrating multiple lines of evidence in the classification of Populus, providing a framework for future taxonomic studies.

Chromosome-Level Genome Assembly and Whole-Genome Resequencing Revealed Contrasting Population Genetic Differentiation of Black Bream (Megalobrama skolkovii) (Teleostei: Cyprinidae) Allopatric and Sympatric to Its Kin Species

The black bream (Megalobrama skolkovii) is an economically important species widely distributed in China, with its geographic populations potentially having undergone differentiations and local adaptations. In this study, we presented a chromosome-level genome assembly of this species and investigated genetic differentiations of its populations that are allopatric (the northern one) and sympatric (the Poyang Lake) to its kin species, the blunt-snout bream (M. amblycephala), using whole genome resequencing analysis. The results showed that the genome size of black bream was 1.13 Gb, very similar to its kin species but larger than its close relatives, the four Chinese major carps. By resequencing individuals from the northern and Poyang Lake populations, we found that the northern population showed lower genetic diversity, larger genetic differentiation, and two sharp historical declines in population size through demographic analysis, indicating the possible bottlenecks after the allopatric isolation. In contrast, the Poyang Lake population, with its higher genetic diversity, higher Tajima's D value, and lower levels of linkage disequilibrium, reflects the ancestral state of black bream. In addition, we also found that the northern population shared more alleles with its kin species, indicating it may retain more ancestral variations. This was further analyzed to be caused by incomplete lineage sorting and ancient introgression. Some key genes related to reproductive processes, body size development, and muscle metabolism were found under selection in the northern population, possibly responsible for its local adaptation. Our findings that the black bream allopatric population had a loss of genetic diversity but retained more ancestral variations can expand our knowledge on population genetic differentiation and give us hints for future genetic conservation.

Multigenerational hybridisation results in heterosis and facilitates adaptive introgression, with no evidence of outbreeding depression in a pair of marine gastropods

Anthropogenic environmental changes continue to threaten species globally. For example, translocation of species has caused unintentional hybridisation, which has contributed to species declines. On the other hand, hybridisation can be used to increase the evolutionary potential of species vulnerable to rapid environmental change, although the benefits of mixing genetically divergent lineages do not come without risks to individual fitness and the long-term viability of populations. Here, we use a combination of genome-wide Single Nucleotide Polymorphism (SNP) markers, mitochondrial DNA sequencing and measurements of growth rate to determine the genetic consequences of hybridisation between two congeneric marine gastropods across 27 years (~18 generations). Multigeneration hybridisation resulted from the introduction of the intertidal periwinkle Bembicium vittatum (a direct developer) into the native range of its congener Bembicium auratum (a species with planktotrophic larval dispersal). Despite significant genetic divergence between the species, we found no direct evidence of outbreeding depression in the admixed population. Instead, we found evidence for heterosis, which dissipated over time. After an initial lag, the frequency of introduced B. vittatum alleles declined dramatically in the hybrid population. However, a few B. vittatum alleles (3.18%) increased significantly in frequency against the overall trend, providing evidence of adaptive introgression. In the context of hybridisation as a conservation management tool, our results provide some evidence of the potential benefits that can be gained, and suggest that the costs due to outbreeding depression can be small.

The long and short of hyperdivergent regions

The increasing prevalence of genome sequencing and assembly has uncovered evidence of hyperdivergent genomic regions - loci with excess genetic diversity - in species across the tree of life. Hyperdivergent regions are often enriched for genes that mediate environmental responses, such as immunity, parasitism, and sensory perception. Especially in self-fertilizing species where the majority of the genome is homozygous, the existence of hyperdivergent regions might imply the historical action of evolutionary forces such as introgression and/or balancing selection. We anticipate that the application of new sequencing technologies, broader taxonomic sampling, and evolutionary modeling of hyperdivergent regions will provide insights into the mechanisms that generate and maintain genetic diversity within and between species.

Introgression drives adaptation to the plateau environment in a subterranean rodent

BACKGROUND

Introgression has repeatedly been shown to play an important role in the adaptation of species to extreme environments, yet how introgression enables rodents with specialized subterranean lifestyle to acclimatize to high altitudes is still unclear. Myospalacinae is a group of subterranean rodents, among which the high-altitude plateau zokors (Eospalax baileyi) and the low-altitude Gansu zokors (E. cansus) are sympatrically distributed in the grassland ecosystems of the Qinghai-Tibet Plateau (QTP). Together, they provide a model for the study of the role of introgression in the adaptation of low-altitude subterranean rodents to high altitudes.

RESULTS

Applying low-coverage whole-genome resequencing and population genetics analyses, we identified evidence of adaptive introgression from plateau zokors into Gansu zokors, which likely facilitated the adaptation of the latter to the high-altitude environment of the QTP. We identified positively selected genes with functions related to energy metabolism, cardiovascular system development, calcium ion transport, and response to hypoxia which likely made critical contributions to adaptation to the plateau environment in both plateau zokors and high-altitude populations of Gansu zokors.

CONCLUSIONS

Introgression of genes associated with hypoxia adaptation from plateau zokors may have played a role in the adaptation of Gansu zokors to the plateau environment. Our study provides new insights into the understanding of adaptive evolution of species on the QTP and the importance of introgression in the adaptation of species to high-altitude environments.

Contemporary intergeneric hybridization and backcrossing among birds-of-paradise

Despite large differences in morphology, behavior and lek-mating strategies the birds-of-paradise are known to hybridize occasionally, even across different genera. Many of these bird-of-paradise hybrids were originally described as distinct species based on large morphological differences when compared to recognized species. Nowadays, these specimens are generally recognized as hybrids based on morphological assessments. Having fascinated naturalists for centuries, hybrid specimens of birds-of-paradise have been collected and the specimens kept in Natural History Collections. In the present study, we utilize this remarkable resource in a museomics framework and evaluate the genomic composition of most described intergeneric hybrids and some intrageneric hybrids. We show that the majority of investigated specimens are first-generation hybrids and that the parental species, in most cases, are in line with prior morphological assessments. We also identify two specimens that are the result of introgressive hybridization between different genera. Additionally, two specimens exhibit hybrid morphologies but have no identifiable signals of hybridization, which may indicate that minor levels of introgression can have large morphological effects. Our findings provide direct evidence of contemporary introgressive hybridization taking place between genera of birds-of-paradise in nature, despite markedly different morphologies and lek-mating behaviors.

Genomic evidence for human-mediated introgressive hybridization and selection in the developed breed

BACKGROUND

The pig (Sus Scrofa) is one of the oldest domesticated livestock species that has undergone extensive improvement through modern breeding. European breeds have advantages in lean meat development and highly-productive body type, whereas Asian breeds possess extraordinary fat deposition and reproductive performance. Consequently, Eurasian breeds have been extensively used to develop modern commercial breeds for fast-growing and high prolificacy. However, limited by the sequencing technology, the genome architecture of some nascent developed breeds and the human-mediated impact on their genomes are still unknown.

RESULTS

Through whole-genome analysis of 178 individuals from an Asian locally developed pig breed, Beijing Black pig, and its two ancestors from two different continents, we found the pervasive inconsistent gene trees and species trees across the genome of Beijing Black pig, which suggests its introgressive hybrid origin. Interestingly, we discovered that this developed breed has more genetic relationships with European pigs and an unexpected introgression from Asian pigs to this breed, which indicated that human-mediated introgression could form the porcine genome architecture in a completely different type compared to native introgression. We identified 554 genomic regions occupied 63.30 Mb with signals of introgression from the Asian ancestry to Beijing Black pig, and the genes in these regions enriched in pathways associated with meat quality, fertility, and disease-resistant. Additionally, a proportion of 7.77% of genomic regions were recognized as regions that have been under selection. Moreover, combined with the results of a genome-wide association study for meat quality traits in the 1537 Beijing Black pig population, two important candidate genes related to meat quality traits were identified. DNAJC6 is related to intramuscular fat content and fat deposition, and RUFY4 is related to meat pH and tenderness.

CONCLUSIONS

Our research provides insight for analyzing the origins of nascent developed breeds and genome-wide selection remaining in the developed breeds mediated by humans during modern breeding.

Potential and progress of studying mountain biodiversity by means of butterfly genetics and genomics

Mountains are rich in biodiversity, and butterflies are species-rich and have a good ecological and evolutionary research foundation. This review addresses the potential and progress of studying mountain biodiversity using butterflies as a model. We discuss the uniqueness of mountain ecosystems, factors influencing the distribution of mountain butterflies, representative genetic and evolutionary models in butterfly research, and evolutionary studies of mountain biodiversity involving butterfly genetics and genomics. Finally, we demonstrate the necessity of studying mountain butterflies and propose future perspectives. This review provides insights for studying the biodiversity of mountain butterflies as well as a summary of research methods for reference.

Secondary Contact, Introgressive Hybridization, and Genome Stabilization in Sticklebacks

Advances in genomic studies have revealed that hybridization in nature is pervasive and raised questions about the dynamics of different genetic and evolutionary factors following the initial hybridization event. While recent research has proposed that the genomic outcomes of hybridization might be predictable to some extent, many uncertainties remain. With comprehensive whole-genome sequence data, we investigated the genetic introgression between 2 divergent lineages of 9-spined sticklebacks (Pungitius pungitius) in the Baltic Sea. We found that the intensity and direction of selection on the introgressed variation has varied across different genomic elements: while functionally important regions displayed reduced rates of introgression, promoter regions showed enrichment. Despite the general trend of negative selection, we identified specific genomic regions that were enriched for introgressed variants, and within these regions, we detected footprints of selection, indicating adaptive introgression. Geographically, we found the selection against the functional changes to be strongest in the vicinity of the secondary contact zone and weaken as a function of distance from the initial contact. Altogether, the results suggest that the stabilization of introgressed variation in the genomes is a complex, multistage process involving both negative and positive selection. In spite of the predominance of negative selection against introgressed variants, we also found evidence for adaptive introgression variants likely associated with adaptation to Baltic Sea environmental conditions.

Leveraging shared ancestral variation to detect local introgression

Introgression is a common evolutionary phenomenon that results in shared genetic material across non-sister taxa. Existing statistical methods such as Patterson's D statistic can detect introgression by measuring an excess of shared derived alleles between populations. The D statistic is effective to detect genome-wide patterns of introgression but can give spurious inferences of introgression when applied to local regions. We propose a new statistic, D+, that leverages both shared ancestral and derived alleles to infer local introgressed regions. Incorporating both shared derived and ancestral alleles increases the number of informative sites per region, improving our ability to identify local introgression. We use a coalescent framework to derive the expected value of this statistic as a function of different demographic parameters under an instantaneous admixture model and use coalescent simulations to compute the power and precision of D+. While the power of D and D+ is comparable, D+ has better precision than D. We apply D+ to empirical data from the 1000 Genome Project and Heliconius butterflies to infer local targets of introgression in humans and in butterflies.

Major patterns in the introgression history of Heliconius butterflies

Gene flow between species, although usually deleterious, is an important evolutionary process that can facilitate adaptation and lead to species diversification. It also makes estimation of species relationships difficult. Here, we use the full-likelihood multispecies coalescent (MSC) approach to estimate species phylogeny and major introgression events in Heliconius butterflies from whole-genome sequence data. We obtain a robust estimate of species branching order among major clades in the genus, including the 'melpomene-silvaniform' group, which shows extensive historical and ongoing gene flow. We obtain chromosome-level estimates of key parameters in the species phylogeny, including species divergence times, present-day and ancestral population sizes, as well as the direction, timing, and intensity of gene flow. Our analysis leads to a phylogeny with introgression events that differ from those obtained in previous studies. We find that Heliconius aoede most likely represents the earliest-branching lineage of the genus and that 'silvaniform' species are paraphyletic within the melpomene-silvaniform group. Our phylogeny provides new, parsimonious histories for the origins of key traits in Heliconius, including pollen feeding and an inversion involved in wing pattern mimicry. Our results demonstrate the power and feasibility of the full-likelihood MSC approach for estimating species phylogeny and key population parameters despite extensive gene flow. The methods used here should be useful for analysis of other difficult species groups with high rates of introgression.

Inferring Historical Introgression with Deep Learning

Resolving phylogenetic relationships among taxa remains a challenge in the era of big data due to the presence of genetic admixture in a wide range of organisms. Rapidly developing sequencing technologies and statistical tests enable evolutionary relationships to be disentangled at a genome-wide level, yet many of these tests are computationally intensive and rely on phased genotypes, large sample sizes, restricted phylogenetic topologies, or hypothesis testing. To overcome these difficulties, we developed a deep learning-based approach, named ERICA, for inferring genome-wide evolutionary relationships and local introgressed regions from sequence data. ERICA accepts sequence alignments of both population genomic data and multiple genome assemblies, and efficiently identifies discordant genealogy patterns and exchanged regions across genomes when compared with other methods. We further tested ERICA using real population genomic data from Heliconius butterflies that have undergone adaptive radiation and frequent hybridization. Finally, we applied ERICA to characterize hybridization and introgression in wild and cultivated rice, revealing the important role of introgression in rice domestication and adaptation. Taken together, our findings demonstrate that ERICA provides an effective method for teasing apart evolutionary relationships using whole genome data, which can ultimately facilitate evolutionary studies on hybridization and introgression.

Introgression between highly divergent fungal sister species

To understand how species evolve and adapt to changing environments, it is important to study gene flow and introgression due to their influence on speciation and radiation events. Here, we apply a novel experimental system for investigating these mechanisms using natural populations. The system is based on two fungal sister species with morphological and ecological similarities occurring in overlapping habitats. We examined introgression between these species by conducting whole genome sequencing of individuals from populations in North America and Europe. We assessed genome-wide nucleotide divergence and performed crossing experiments to study reproductive barriers. We further used ABBA-BABA statistics together with a network analysis to investigate introgression, and conducted demographic modelling to gain insight into divergence times and introgression events. The results revealed that the species are highly divergent and incompatible in vitro. Despite this, small regions of introgression were scattered throughout the genomes and one introgression event likely involves a ghost population (extant or extinct). This study demonstrates that introgression can be found among divergent species and that population histories can be studied without collections of all the populations involved. Moreover, the experimental system is shown to be a useful tool for research on reproductive isolation in natural populations.

Mitochondrial introgression and mito-nuclear discordance obscured the closely related species boundaries in Cletus Stål from China (Heteroptera: Coreidae)

Accurate taxonomy and delimitation are of great importance for pest control strategies and management programs. Here, we focus on Cletus (Insecta: Hemiptera: Coreidae), which includes many crop pests. The species boundaries still conflict and only cytochrome c oxidase subunit I (COI) barcoding has been previously used for molecular studies. We generated new mitochondrial genome and nuclear genome-wide SNPs to explore the species boundaries of 46 Cletus samples from China using multiple species delimitation approaches. All results recovered a monophyly with high support, except for two closely related species in clade I - C. punctiger and C. graminis. Mitochondrial data demonstrated admixture in clade I, while genome-wide SNPs unambiguously identified two separate species, which were confirmed by morphological classification. Inconsistent nuclear and mitochondrial data indicated mito-nuclear discordance. Mitochondrial introgression is the most likely explanation, and more extensive sampling and more comprehensive data are needed to ascertain a pattern. Accurate species delimitation will shed light on species status; thus, an accurate taxonomy is of particular concern, as there is a pressing need to implement precise control of agricultural pests and to perform further research on diversification.

A caste differentiation mutant elucidates the evolution of socially parasitic ants

Most ant species have two distinct female castes-queens and workers-yet the developmental and genetic mechanisms that produce these alternative phenotypes remain poorly understood. Working with a clonal ant, we discovered a variant strain that expresses queen-like traits in individuals that would normally become workers. The variants show changes in morphology, behavior, and fitness that cause them to rely on workers in wild-type (WT) colonies for survival. Overall, they resemble the queens of many obligately parasitic ants that have evolutionarily lost the worker caste and live inside colonies of closely related hosts. The prevailing theory for the evolution of these workerless social parasites is that they evolve from reproductively isolated populations of facultative intermediates that acquire parasitic phenotypes in a stepwise fashion. However, empirical evidence for such facultative ancestors remains weak, and it is unclear how reproductive isolation could gradually arise in sympatry. In contrast, we isolated these variants just a few generations after they arose within their WT parent colony, implying that the complex phenotype reported here was induced in a single genetic step. This suggests that a single genetic module can decouple the coordinated mechanisms of caste development, allowing an obligately parasitic variant to arise directly from a free-living ancestor. Consistent with this hypothesis, the variants have lost one of the two alleles of a putative supergene that is heterozygous in WTs. These findings provide a plausible explanation for the evolution of ant social parasites and implicate new candidate molecular mechanisms for ant caste differentiation.

Prevalent Introgression Underlies Convergent Evolution in the Diversification of Pungitius Sticklebacks

New mutations and standing genetic variations contribute significantly to repeated phenotypic evolution in sticklebacks. However, less is known about the role of introgression in this process. We analyzed taxonomically and geographically comprehensive genomic data from Pungitius sticklebacks to decipher the extent of introgression and its consequences for the diversification of this genus. Our results demonstrate that introgression is more prevalent than suggested by earlier studies. Although gene flow was generally bidirectional, it was often asymmetric and left unequal genomic signatures in hybridizing species, which might, at least partly, be due to biased hybridization and/or population size differences. In several cases, introgression of variants from one species to another was accompanied by transitions of pelvic and/or lateral plate structures-important diagnostic traits in Pungitius systematics-and frequently left signatures of adaptation in the core gene regulatory networks of armor trait development. This finding suggests that introgression has been an important source of genetic variation and enabled phenotypic convergence among Pungitius sticklebacks. The results highlight the importance of introgression of genetic variation as a source of adaptive variation underlying key ecological and taxonomic traits. Taken together, our study indicates that introgression-driven convergence likely explains the long-standing challenges in resolving the taxonomy and systematics of this small but phenotypically highly diverse group of fish.

Divergence and introgression among the virilis group of Drosophila

Speciation with gene flow is now widely regarded as common. However, the frequency of introgression between recently diverged species and the evolutionary consequences of gene flow are still poorly understood. The virilis group of Drosophila contains 12 species that are geographically widespread and show varying levels of prezygotic and postzygotic isolation. Here, we use de novo genome assemblies and whole-genome sequencing data to resolve phylogenetic relationships and describe patterns of introgression and divergence across the group. We suggest that the virilis group consists of three, rather than the traditional two, subgroups. Some genes undergoing rapid sequence divergence across the group were involved in chemical communication and desiccation tolerance, and may be related to the evolution of sexual isolation and adaptation. We found evidence of pervasive phylogenetic discordance caused by ancient introgression events between distant lineages within the group, and more recent gene flow between closely related species. When assessing patterns of genome-wide divergence in species pairs across the group, we found no consistent genomic evidence of a disproportionate role for the X chromosome as has been found in other systems. Our results show how ancient and recent introgressions confuse phylogenetic reconstruction, but may play an important role during early radiation of a group.

Genetic introgression between different groups reveals the differential process of Asian cultivated rice

The Asian cultivated rice consists of two major subspecies: indica and japonica. There are already many reports about the existence of genetic introgression between the two subspecies. They propose some possible introgression-related genes from the comparison of population parameters. This study uses the genome-wide variation data of Asia cultivated rice to investigate their genetic introgression on the whole genome level. We detect a total of 13 significantly high introgression loci between the tropical japonica and indica populations. Two different methods are used to identify the genetic introgression regions. For most of the detected introgression regions, they generally get consistent results. Some previously known introgression genes are detected in the identified introgression loci, such as heat resistance gene TT1 and GLW7. The biological functions for these genetic introgression regions are annotated by the published QTL mapping results. We find that genetic introgression plays a vital role in the determination of both the phenotype and the domestication process of different groups. Our study also provides useful information and resources for the study of rice gene function and the domestication process.

Early stages of sympatric homoploid hybrid speciation in crater lake cichlid fishes

Homoploid hybrid speciation (i.e., hybrid speciation without a change in ploidy) has traditionally been considered to be rare in animals. Only few accepted empirical examples of homoploid hybrid speciation in nature exist, and in only one previous case (insects) was it convincingly shown that this process occurred in complete sympatry. Here, we report an instance of sympatric homoploid hybrid speciation in Midas cichlid fishes in Crater Lake Xiloá, Nicaragua. The hybrid lineage, albeit at an early stage of speciation, has genomically and phenotypically diverged from both of its two parental species. Together with a distinct stable isotope signature this suggests that this hybrid lineages occupies a different trophic niche compared to the other sympatric Midas cichlid species in Crater Lake Xiloá.

Here, the authors present a rare example of sympatric homoploid hybrid speciation, without ploidy changes, in the Midas cichlid fishes from Nicaragua. Midas cichlid hybrids occupy a different ecological niche, likely facilitated by body shape adaptations.

The evolution and diversification of oakleaf butterflies

Oakleaf butterflies in the genus Kallima have a polymorphic wing phenotype, enabling these insects to masquerade as dead leaves. This iconic example of protective resemblance provides an interesting evolutionary paradigm that can be employed to study biodiversity. We integrated multi-omic data analyses and functional validation to infer the evolutionary history of Kallima species and investigate the genetic basis of their variable leaf wing patterns. We find that Kallima butterflies diversified in the eastern Himalayas and dispersed to East and Southeast Asia. Moreover, we find that leaf wing polymorphism is controlled by the wing patterning gene cortex, which has been maintained in Kallima by long-term balancing selection. Our results provide macroevolutionary and microevolutionary insights into a model species originating from a mountain ecosystem.

Interpreting phylogenetic conflict: Hybridization in the most speciose genus of lichen-forming fungi

While advances in sequencing technologies have been invaluable for understanding evolutionary relationships, increasingly large genomic data sets may result in conflicting evolutionary signals that are often caused by biological processes, including hybridization. Hybridization has been detected in a variety of organisms, influencing evolutionary processes such as generating reproductive barriers and mixing standing genetic variation. Here, we investigate the potential role of hybridization in the diversification of the most speciose genus of lichen-forming fungi, Xanthoparmelia. As Xanthoparmelia is projected to have gone through recent, rapid diversification, this genus is particularly suitable for investigating and interpreting the origins of phylogenomic conflict. Focusing on a clade of Xanthoparmelia largely restricted to the Holarctic region, we used a genome skimming approach to generate 962 single-copy gene regions representing over 2 Mbp of the mycobiont genome. From this genome-scale dataset, we inferred evolutionary relationships using both concatenation and coalescent-based species tree approaches. We also used three independent tests for hybridization. Although different species tree reconstruction methods recovered largely consistent and well-supported trees, there was widespread incongruence among individual gene trees. Despite challenges in differentiating hybridization from ILS in situations of recent rapid radiations, our genome-wide analyses detected multiple potential hybridization events in the Holarctic clade, suggesting one possible source of trait variability in this hyperdiverse genus. This study highlights the value in using a pluralistic approach for characterizing genome-scale conflict, even in groups with well-resolved phylogenies, while highlighting current challenges in detecting the specific impacts of hybridization.

Genome-wide analyses of introgression between two sympatric Asian oak species

Introgression can be an important source of new alleles for adaption under rapidly changing environments, perhaps even more important than standing variation. Though introgression has been extensively studied in many plants and animals, key questions on the underlying mechanisms of introgression still remain unanswered. In particular, we are yet to determine the genomic distribution of introgressed regions along the genome; whether the extent and patterns of introgression are influenced by ecological factors; and when and how introgression contributes to adaptation. Here, we generated high-quality genomic resources for two sympatric widespread Asian oak species, Quercus acutissima and Q. variabilis, sampled in multiple forests to study introgression between them. We show that introgressed regions are broadly distributed across the genome. Introgression was affected by genetic divergence between pairs of populations and by the similarity of the environments in which they live-populations occupying similar ecological sites tended to share the same introgressed regions. Introgressed genomic footprints of adaptation were preferentially located in regions with suppressed recombination rate. Introgression probably confers adaptation in these oak populations by introducing allelic variation in cis-regulatory elements, in particular through transposable element insertions, thereby altering the regulation of genes related to stress. Our results provide new avenues of research for uncovering mechanisms of adaptation due to hybridization in sympatric species.

Phylogenomics and Biogeography of Populus Based on Comprehensive Sampling Reveal Deep-Level Relationships and Multiple Intercontinental Dispersals

Populus not only has significant economic and ecological values, but also serves as a model tree that is widely used in the basic research of tree growth, physiology, and genetics. However, high levels of morphological variation and extensive interspecific hybridization of Populus pose an obstacle for taxonomy, and also to the understanding of phylogenetic interspecific relationships and biogeographical history. In this study, a total of 103 accessions representing almost all the wild species of Populus were collected and whole-genome re-sequenced to examine the phylogenetic relationships and biogeography history. On the basis of 12,916,788 nuclear single nucleotide polymorphisms (SNPs), we reconstructed backbone phylogenies using concatenate and coalescent methods, we highly disentangled the species relationships of Populus, and several problematic taxa were treated as species complexes. Furthermore, the phylogeny of the chloroplast genome showed extensive discordance with the trees from the nuclear genome data, and due to extensive chloroplast capture and hybridization of Populus species, plastomes could not accurately evaluate interspecies relationships. Ancient gene flow between clades and some hybridization events were also identified by ABBA-BABA analysis. The reconstruction of chronogram and ancestral distributions suggested that North America was the original region of this genus, and subsequent long dispersal and migration across land bridges were contributed to the modern range of Populus. The diversification of Populus mainly occurred in East Asia in recent 15 Ma, possibly promoted by the uplift of the Tibetan Plateau. This study provided comprehensive evidence on the phylogeny of Populus and proposed a four-subgeneric classification and a new status, subgenus Abaso. Meanwhile, ancestral distribution reconstruction with nuclear data advanced the understanding of the biogeographic history of Populus.

Phylogenomic approaches to detecting and characterizing introgression

Phylogenomics has revealed the remarkable frequency with which introgression occurs across the tree of life. These discoveries have been enabled by the rapid growth of methods designed to detect and characterize introgression from whole-genome sequencing data. A large class of phylogenomic methods makes use of data across species to infer and characterize introgression based on expectations from the multispecies coalescent. These methods range from simple tests, such as the D-statistic, to model-based approaches for inferring phylogenetic networks. Here, we provide a detailed overview of the various signals that different modes of introgression are expected leave in the genome, and how current methods are designed to detect them. We discuss the strengths and pitfalls of these approaches and identify areas for future development, highlighting the different signals of introgression, and the power of each method to detect them. We conclude with a discussion of current challenges in inferring introgression and how they could potentially be addressed.

The effects of introgression across thousands of quantitative traits revealed by gene expression in wild tomatoes

It is now understood that introgression can serve as powerful evolutionary force, providing genetic variation that can shape the course of trait evolution. Introgression also induces a shared evolutionary history that is not captured by the species phylogeny, potentially complicating evolutionary analyses that use a species tree. Such analyses are often carried out on gene expression data across species, where the measurement of thousands of trait values allows for powerful inferences while controlling for shared phylogeny. Here, we present a Brownian motion model for quantitative trait evolution under the multispecies network coalescent framework, demonstrating that introgression can generate apparently convergent patterns of evolution when averaged across thousands of quantitative traits. We test our theoretical predictions using whole-transcriptome expression data from ovules in the wild tomato genus Solanum. Examining two sub-clades that both have evidence for post-speciation introgression, but that differ substantially in its magnitude, we find patterns of evolution that are consistent with histories of introgression in both the sign and magnitude of ovule gene expression. Additionally, in the sub-clade with a higher rate of introgression, we observe a correlation between local gene tree topology and expression similarity, implicating a role for introgressed cis-regulatory variation in generating these broad-scale patterns. Our results reveal a general role for introgression in shaping patterns of variation across many thousands of quantitative traits, and provide a framework for testing for these effects using simple model-informed predictions.

It is now known from studying large genetic datasets that species often hybridize and cross with each other over many generations – a phenomenon known as introgression. Introgression introduces new genetic variation into a population, and this variation can cause traits to be shared among the introgressing species. When researchers study the evolution of trait variation among species, this source of trait sharing is rarely accounted for. Here, we present a statistical model of the effects of introgression on trait variation. This model predicts that, when averaged across many thousands of traits, introgressing species are consistently more similar than expected from standard approaches. Researchers studying gene expression often consider the expression of many thousands of genes, making this a case where the expected effects of introgression are likely to manifest. We tested our model prediction using ovule gene expression data from the wild tomato genus Solanum, in two groups of species with evidence of historical introgression. We found that patterns of expression similarity in both groups are consistent with their histories of introgression and the predictions from our model. Our results highlight the importance of accounting for introgression as a source of trait variation among species.

Spatial and host related genomic variation in partially sympatric cactophagous moth species

Surveys of patterns of genetic variation in natural sympatric and allopatric populations of recently diverged species are necessary to understand the processes driving intra- and interspecific diversification. The South American moths Cactoblastis cactorum, Cactoblastis doddi and Cactoblastis bucyrus are specialized in the use of cacti as host plants. These species have partially different geographic ranges and differ in patterns of host plant use. However, there are areas that overlap, particularly, in northwestern Argentina, where they are sympatric. Using a combination of genome-wide SNPs and mitochondrial data we assessed intra and interspecific genetic variation and investigated the relative roles of geography and host plants on genetic divergence. We also searched for genetic footprints of hybridization between species. We identified three well delimited species and detected signs of hybridization in the area of sympatry. Our results supported a hypothetical scenario of allopatric speciation in the generalist C. cactorum and genetic interchange during secondary geographic contact with the pair of specialists C. bucyrus and C. doddi that probably speciated sympatrically. In both cases, adaptation to new host plants probably played an important role in speciation. The results also suggested the interplay of geography and host plant use as drivers of divergence and limiting gene flow at intra and interspecific levels.

Variational inference using approximate likelihood under the coalescent with recombination

Coalescent methods are proven and powerful tools for population genetics, phylogenetics, epidemiology, and other fields. A promising avenue for the analysis of large genomic alignments, which are increasingly common, is coalescent hidden Markov model (coalHMM) methods, but these methods have lacked general usability and flexibility. We introduce a novel method for automatically learning a coalHMM and inferring the posterior distributions of evolutionary parameters using black-box variational inference, with the transition rates between local genealogies derived empirically by simulation. This derivation enables our method to work directly with three or four taxa and through a divide-and-conquer approach with more taxa. Using a simulated data set resembling a human–chimp–gorilla scenario, we show that our method has comparable or better accuracy to previous coalHMM methods. Both species divergence times and population sizes were accurately inferred. The method also infers local genealogies, and we report on their accuracy. Furthermore, we discuss a potential direction for scaling the method to larger data sets through a divide-and-conquer approach. This accuracy means our method is useful now, and by deriving transition rates by simulation, it is flexible enough to enable future implementations of various population models.

A widely diverged locus involved in locomotor adaptation in Heliconius butterflies

Heliconius butterflies have undergone adaptive radiation and therefore serve as an excellent system for exploring the continuum of speciation and adaptive evolution. However, there is a long-lasting paradox between their convergent mimetic wing patterns and rapid divergence in speciation. Here, we characterize a locus that consistently displays high divergence among Heliconius butterflies and acts as an introgression hotspot. We further show that this locus contains multiple genes related to locomotion and conserved in Lepidoptera. In light of these findings, we consider that locomotion traits may be under selection, and if these are heritable traits that are selected for, then they might act as species barriers.

A widely diverged locus involved in locomotor adaptation in Heliconius butterflies

Heliconius butterflies have undergone adaptive radiation and therefore serve as an excellent system for exploring the continuum of speciation and adaptive evolution. However, there is a long-lasting paradox between their convergent mimetic wing patterns and rapid divergence in speciation. Here, we characterize a locus that consistently displays high divergence among Heliconius butterflies and acts as an introgression hotspot. We further show that this locus contains multiple genes related to locomotion and conserved in Lepidoptera. In light of these findings, we consider that locomotion traits may be under selection, and if these are heritable traits that are selected for, then they might act as species barriers.

Genomic introgression from a distant congener in the Levant fritillary butterfly, Melitaea acentria

Introgressive hybridization is more common in nature than previously thought, and its role and creative power in evolution is hotly discussed but not completely understood. Introgression occurs more frequently in sympatry between recently diverged taxa, or when the speciation process has not yet been completed. However, there are relatively few documented cases of hybridization that erodes reproductive barriers between distantly related species. Here, we use whole genome and mitochondrial data to examine how introgression from a distant congener affects pattern of genetic differentiation in the Levant fritillary butterfly Melitaea acentria. We show that this local taxon has evolved as a peripatric geographic isolate of the widespread Melitaea persea, and that there has been significant unidirectional gene flow from the sympatric, nonclosely related Melitaea didyma to M. acentria. We found direct evidence of ongoing sporadic hybridization between M. didyma and M. acentria, which are separated by at least 5 million years of independent evolution. Elevated differentiation and lower level of introgression on the sex Z chromosome compared to autosomes suggest that the Z chromosome has accumulated loci acting as intrinsic postzygotic barriers. Our results show that introgression from M. didyma has been an additional source of nucleotide diversity in the M. acentria population, providing material for drift and selection.

Rampant Genome-Wide Admixture across the Heliconius Radiation

How frequent is gene flow between species? The pattern of evolution is typically portrayed as a phylogenetic tree, yet gene flow between good species may be an important mechanism in diversification, spreading adaptive traits and leading to a complex pattern of phylogenetic incongruence. This process has thus far been studied mainly among a few closely related species, or in geographically restricted areas such as islands, but not on the scale of a continental radiation. Using a genomic representation of 40 out of 47 species in the genus, we demonstrate that admixture has played a role throughout the evolution of the charismatic Neotropical butterflies Heliconius. Modeling of phylogenetic networks based on the exome uncovers up to 13 instances of interspecific gene flow. Admixture is detected among the relatives of Heliconius erato, as well as between the ancient lineages leading to modern clades. Interspecific gene flow played a role throughout the evolution of the genus, although the process has been most frequent in the clade of Heliconius melpomene and relatives. We identify Heliconius hecalesia and relatives as putative hybrids, including new evidence for introgression at the loci controlling the mimetic wing patterns. Models accounting for interspecific gene flow yield a more complete picture of the radiation as a network, which will improve our ability to study trait evolution in a realistic comparative framework.

CARACTERIZACIÓN GENÉTICA DE LA POBLACIÓN DE Heliconius sara (Nymphalidae) EN LA ISLA GORGONA, COLOMBIA

La estructura genética de poblaciones de mariposas con distribución en islas y sus pares continentales ha sido poco documentada para el neotrópico. Este estudio presenta la caracterización de una población de Heliconius sara con distribución en la Isla Gorgona, ubicada en la región del Pacífico Oriental Colombiano. Para esto se examinaron secuencias parciales de un marcador mitocondrial incluyendo información obtenida del GenBank. Se comparó la diversidad y estructura genética con sus conespecíficos continentales y también con congéneres, con los que comparte un ancestro común cercano en el clado Sapho-Sara. Para el análisis de diversidad y estructura genética se realizó un análisis molecular de varianza. Este análisis muestra que la distancia entre la población de la isla y sus pares en el continente es consistente con la variación intraespecífica observada en otras especies del género Heliconius. Para la reconstrucción de la genealogía y datación reciente en el Pleistoceno superior del grupo monofilético de secuencias de H. sara, se realizó un análisis de inferencia bayesiana, así como una de máxima verosimilitud. Del análisis demográfico se seleccionó un modelo histórico de flujo asimétrico desde la isla hacia el continente que sugiere baja resistencia de la discontinuidad geográfica a la dispersión de esta mariposa diurna desde la isla. Este es el primer estudio en examinar un posible evento de aislamiento de una población insular de mariposas en Colombia.

Phylogenomic approach reveals strong signatures of introgression in the rapid diversification of neotropical true fruit flies (Anastrepha: Tephritidae)

New sequencing techniques have allowed us to explore the variation on thousands of genes and elucidate evolutionary relationships of lineages even in complex scenarios, such as when there is rapid diversification. That seems to be the case of species in the genus Anastrepha, which shows great species diversity that has been divided into 21 species groups, several of which show wide geographical distribution. The fraterculus group has several economically important species and it is also an outstanding model for speciation studies, since it includes several lineages that have diverged recently possibly in the presence of interspecific gene flow. Our main goal is to test whether we can infer phylogenetic relationships of recently diverged taxa with gene flow, such as what is expected for the fraterculus group and determine whether certain genes remain informative even in this complex scenario. An analysis of thousands of orthologous genes derived from transcriptome datasets of 10 different lineages across the genus, including some of the economically most important pests, revealed signals of incomplete lineage sorting, vestiges of ancestral introgression between more distant lineages and ongoing gene flow between closely related lineages. Though these patterns affect the phylogenetic signal, the phylogenomic inferences consistently show that the morphologically identified species here investigated are in different evolutionary lineages, with the sole exception involving Brazilian lineages of A. fraterculus, which has been suggested to be a complex assembly of cryptic species. A tree space analysis suggested that genes with greater phylogenetic resolution have evolved under similar selection pressures and are more resilient to intraspecific gene flow, which would make it more likely that these genomic regions may be useful for identifying fraterculus group lineages. Our findings help establish relationships among the most important Anastrepha species groups, as well as bring further data to indicate that the diversification of fraterculus group lineages, and even other lineages in the genus Anastrepha, has been strongly influenced by interspecific gene flow.

Doublesex Mediates the Development of Sex-Specific Pheromone Organs in Bicyclus Butterflies via Multiple Mechanisms

The Bicyclus lineage of satyrid butterflies exhibits male-specific traits, the scent organs, used for chemical communication during courtship. These organs consist of tightly packed brush-like scales (hair-pencils) that rub against scent patches to disperse pheromones, but the evolution and molecular basis of these organ’s male-limited development remains unknown. Here, we examine the evolution of the number and location of the scent patches and hair-pencils within 53 species of Bicyclus butterflies, and the involvement of the sex determinant gene doublesex (dsx) in scent organ development in Bicyclus anynana using CRISPR/Cas9. We show that scent patches and hair-pencils arose via multiple, independent gains, in a correlated manner. Further, an initially nonsex-specific Dsx protein expression pattern in developing wing discs becomes male-specific and spatially refined to areas that develop the scent patches. Functional perturbations of dsx show that this gene activates patch development in males whereas hair-pencils develop in both sexes without Dsx input. Dsx in females is, instead, required to repress hair-pencils whereas Dsx in males regulates minor aspects of its development. These findings suggest that the patches and hair-pencils evolve as correlated composite organs presumably due to their functional integration. Divergence in the function of dsx isoforms occurred across the sexes, where the male isoform promotes patch development in males and the female isoform represses hair-pencil development in females, both leading to the development of male-limited traits. Furthermore, evolution in number of patches in males is due to the evolution of spatial regulation of dsx.

The natural hybridization between species Ligularia nelumbifolia and Cremanthodium stenoglossum (Senecioneae, Asteraceae) suggests underdeveloped reproductive isolation and ambiguous intergeneric boundary

Natural hybridization is frequent in plants and is considered an important factor facilitating speciation. The natural intergeneric hybridization between Ligularia and Cremanthodium was previously confirmed using a couple of DNA markers. However, the mechanism of this intergeneric hybridization and the role of reproductive isolation in the process of hybridization remain unclear. Here we used double-digest restriction site-associated DNA sequencing (ddRAD-seq) to further quantify the occurrence of hybridization, the genetic structure of the hybrid population and the role of reproductive isolation between Ligularia nelumbifolia and Cremanthodium stenoglossum. The results based on the ddRAD-seq SNP data sets indicated that hybridization between L. nelumbifolia and C. stenoglossum was restricted to F₁s, and no gene introgression was identified between these two species. STRUCTURE analysis and maximum likelihood (ML) tree results showed a slightly larger genetic contribution of L. nelumbifolia to putative hybrid F₁s. We deduced that the reproductive isolation between these two parent species is not well-developed but still strong enough to maintain the genetic integrity of the species, and that their F₁s are sterile or with low fertility. Given the poorly resolved phylogenetic relationship between Ligularia and Cremanthodium, the occurrence of natural hybridization between L. nelumbifolia and C. stenoglossum may provide new insights into the re-circumscription and re-delimitation of these two genera.

Disentangling Population History and Character Evolution among Hybridizing Lineages

Understanding the origin and maintenance of adaptive phenotypic novelty is a central goal of evolutionary biology. However, both hybridization and incomplete lineage sorting can lead to genealogical discordance between the regions of the genome underlying adaptive traits and the remainder of the genome, decoupling inferences about character evolution from population history. Here, to disentangle these effects, we investigated the evolutionary origins and maintenance of Batesian mimicry between North American admiral butterflies (Limenitis arthemis) and their chemically defended model (Battus philenor) using a combination of de novo genome sequencing, whole-genome resequencing, and statistical introgression mapping. Our results suggest that balancing selection, arising from geographic variation in the presence or absence of the unpalatable model, has maintained two deeply divergent color patterning haplotypes that have been repeatedly sieved among distinct mimetic and nonmimetic lineages of Limenitis via introgressive hybridization.

Natural hybridization between two butterfly bushes in Tibet: dominance of F1 hybrids promotes strong reproductive isolation

BACKGROUND

It has been recognized that a certain amount of habitat disturbance is a facilitating factor for the occurrence of natural hybridization, yet to date we are unaware of any studies exploring hybridization and reproductive barriers in those plants preferentially occupying disturbed habitats. Buddleja plants (also called butterfly bush) generally do grow in disturbed habitats, and several species with hybrid origin have been proposed, based solely on morphological evidence.

RESULTS

In the present study, we test the hypothesis that B. × wardii is of natural hybridization origin in two sympatric populations of three taxa including B. × wardii and its parents (B. alternifolia and B. crispa) plus 4 referenced parental populations, using four nuclear genes and three chloroplast intergenic spacers, as well as with 10 morphological characters. Our results suggest that at both sites B. × wardii is likely to be a hybrid between B. alternifolia and B. crispa, and moreover, we confirm that most of the hybrids examined are F1s. That these plants are F1s is further supported by morphology, as no transgressive characters were detected. B. crispa was found to be the maternal parent in the Bahe (BH) population, from cpDNA evidence. However, in the Taji (TJ) population, the direction of hybridization was difficult to establish due to the shared cpDNA haplotypes between B. alternifolia and B. crispa, however we still predicted a similar unidirectional hybridization pattern due to results from cross-specific pollination treatments which supported the "SI × SC rule".

CONCLUSIONS

The presence of mainly F1 hybrids can successfully impede gene flow and thus maintain species boundaries in parental species in a typical distribution of Buddleja, i.e. in disturbed habitats.

Rapid colour shift by reproductive character displacement in Cupido butterflies

Reproductive character displacement occurs when competition for successful breeding imposes a divergent selection on the interacting species, causing a divergence of reproductive traits. Here, we show that a disputed butterfly taxon is actually a case of male wing colour shift, apparently produced by reproductive character displacement. Using double digest restriction-site associated DNA sequencing and mitochondrial DNA sequencing we studied four butterfly taxa of the subgenus Cupido (Lepidoptera: Lycaenidae): Cupido minimus and the taxon carswelli, both characterized by brown males and females, plus C. lorquinii and C. osiris, both with blue males and brown females. Unexpectedly, taxa carswelli and C. lorquinii were close to indistinguishable based on our genomic and mitochondrial data, despite displaying strikingly different male coloration. In addition, we report and analysed a brown male within the C. lorquinii range, which demonstrates that the brown morph occurs at very low frequency in C. lorquinii. Such evidence strongly suggests that carswelli is conspecific with C. lorquinii and represents populations with a fixed male brown colour morph. Considering that these brown populations occur in sympatry with or very close to the blue C. osiris, and that the blue C. lorquinii populations never do, we propose that the taxon carswelli could have lost the blue colour due to reproductive character displacement with C. osiris. Since male colour is important for conspecific recognition during courtship, we hypothesize that the observed colour shift may eventually trigger incipient speciation between blue and brown populations. Male colour seems to be an evolutionarily labile character in the Polyommatinae, and the mechanism described here might be at work in the wide diversification of this subfamily of butterflies.

A phylogenomic study of Steganinae fruit flies (Diptera: Drosophilidae): strong gene tree heterogeneity and evidence for monophyly

BACKGROUND

The Drosophilidae family is traditionally divided into two subfamilies: Drosophilinae and Steganinae. This division is based on morphological characters, and the two subfamilies have been treated as monophyletic in most of the literature, but some molecular phylogenies have suggested Steganinae to be paraphyletic. To test the paraphyletic-Steganinae hypothesis, here, we used genomic sequences of eight Drosophilidae (three Steganinae and five Drosophilinae) and two Ephydridae (outgroup) species and inferred the phylogeny for the group based on a dataset of 1,028 orthologous genes present in all species (> 1,000,000 bp). This dataset includes three genera that broke the monophyly of the subfamilies in previous works. To investigate possible biases introduced by small sample sizes and automatic gene annotation, we used the same methods to infer species trees from a set of 10 manually annotated genes that are commonly used in phylogenetics.

RESULTS

Most of the 1,028 gene trees depicted Steganinae as paraphyletic with distinct topologies, but the most common topology depicted it as monophyletic (43.7% of the gene trees). Despite the high levels of gene tree heterogeneity observed, species tree inference in ASTRAL, in PhyloNet, and with the concatenation approach strongly supported the monophyly of both subfamilies for the 1,028-gene dataset. However, when using the concatenation approach to infer a species tree from the smaller set of 10 genes, we recovered Steganinae as a paraphyletic group. The pattern of gene tree heterogeneity was asymmetrical and thus could not be explained solely by incomplete lineage sorting (ILS).

CONCLUSIONS

Steganinae was clearly a monophyletic group in the dataset that we analyzed. In addition to ILS, gene tree discordance was possibly the result of introgression, suggesting complex branching processes during the early evolution of Drosophilidae with short speciation intervals and gene flow. Our study highlights the importance of genomic data in elucidating contentious phylogenetic relationships and suggests that phylogenetic inference for drosophilids based on small molecular datasets should be performed cautiously. Finally, we suggest an approach for the correction and cleaning of BUSCO-derived genomic datasets that will be useful to other researchers planning to use this tool for phylogenomic studies.

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.

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 roles of hybridization and habitat fragmentation in the evolution of Brazil's enigmatic longwing butterflies, Heliconius nattereri and H. hermathena

BACKGROUND

Heliconius butterflies are widely distributed across the Neotropics and have evolved a stunning array of wing color patterns that mediate Müllerian mimicry and mating behavior. Their rapid radiation has been strongly influenced by hybridization, which has created new species and allowed sharing of color patterning alleles between mimetic species pairs. While these processes have frequently been observed in widespread species with contiguous distributions, many Heliconius species inhabit patchy or rare habitats that may strongly influence the origin and spread of species and color patterns. Here, we assess the effects of historical population fragmentation and unique biology on the origins, genetic health, and color pattern evolution of two rare and sparsely distributed Brazilian butterflies, Heliconius hermathena and Heliconius nattereri.

RESULTS

We assembled genomes and re-sequenced whole genomes of eight H. nattereri and 71 H. hermathena individuals. These species harbor little genetic diversity, skewed site frequency spectra, and high deleterious mutation loads consistent with recent population bottlenecks. Heliconius hermathena consists of discrete, strongly isolated populations that likely arose from a single population that dispersed after the last glacial maximum. Despite having a unique color pattern combination that suggested a hybrid origin, we found no genome-wide evidence that H. hermathena is a hybrid species. However, H. hermathena mimicry evolved via introgression, from co-mimetic Heliconius erato, of a small genomic region upstream of the color patterning gene cortex.

CONCLUSIONS

Heliconius hermathena and H. nattereri population fragmentation, potentially driven by historical climate change and recent deforestation, has significantly reduced the genetic health of these rare species. Our results contribute to a growing body of evidence that introgression of color patterning alleles between co-mimetic species appears to be a general feature of Heliconius evolution.

High-density genetic variation maps reveal the correlation between asymmetric interspecific introgressions and improvement of agronomic traits in Upland and Pima cotton varieties developed in Xinjiang, China

The two new world tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, are cultivated worldwide and are characterised by a high yield and superior fibre quality, respectively. Historical genetic introgression has been reported between them; however, the existence of introgression and its genetic effects on agronomic traits remain unclear with regard to independent breeding of G. hirsutum (Upland cotton) and G. barbadense (Pima cotton) elite cultivars. We collected 159 G. hirsutum and 70 G. barbadense cultivars developed in Xinjiang, China, along with 30 semi-wild accessions of G. hirsutum, to perform interspecific introgression tests, intraspecific selection analyses and genome-wide association studies (GWAS) with fibre quality and yield component traits in multiple environments. In total, we identified seven interspecific introgression events and 52 selective sweep loci in G. hirsutum, as well as 17 interspecific introgression events and 19 selective sweep loci in G. barbadense. Correlation tests between agronomic traits and introgressions showed that introgression loci were mutually beneficial for the improvement of fibre quality and yield traits in both species. In addition, the phenotypic effects of four interspecific introgression events could be detected by intraspecific GWAS, with Gb_INT13 significantly improving fibre yield in G. barbadense. The present study describes the landscape of genetic introgression and selection between the two species, and highlights the genetic effects of introgression among populations, which can be used for future improvement of fibre yield and quality in G. barbadense and G. hirsutum, respectively.

Genome-wide data reveal discordant mitonuclear introgression in the intermediate horseshoe bat (Rhinolophus affinis)

Closely related taxa often exhibit mitonuclear discordance attributed to introgression of mitochondrial DNA (mtDNA), yet few studies have considered the underlying causes of mtDNA introgression. Here we test for demographic versus adaptive processes as explanations for mtDNA introgression in three subspecies of the intermediate horseshoe bat (Rhinolophus affinis). We generated sequences of 1692 nuclear genes and 13 mitochondrial protein-coding genes for 48 individuals. Phylogenetic reconstructions based on 320 exon sequences and 2217 single nucleotide polymorphisms (SNPs) both revealed conflicts between the species tree and mtDNA tree. These results, together with geographic patterns of mitonuclear discordance, and shared identical or near-identical mtDNA sequences, suggest extensive introgression of mtDNA between the two parapatric mainland subspecies. Under demographic hypotheses, we would also expect to uncover traces of ncDNA introgression, however, population structure and gene flow analyses revealed little nuclear admixture. Furthermore, we found inconsistent estimates of the timing of population expansion and that of the most recent common ancestor for the clade containing introgressed haplotypes. Without a clear demographic explanation, we also examined whether introgression likely arises from adaptation. We found that five mtDNA genes contained fixed amino acid differences between introgressed and non-introgressed individuals, including putative positive selection found in one codon, although this did not show introgression. While our evidence for rejecting demographic hypotheses is arguably stronger than that for rejecting adaptation, we find no definitive support for either explanation. Future efforts will focus on larger-scale resequencing to decipher the underlying causes of discordant mitonuclear introgression in this system.

Genome-wide SNP Data Reveal an Overestimation of Species Diversity in a Group of Hawkmoths

The interface between populations and evolving young species continues to generate much contemporary debate in systematics depending on the species concept(s) applied but which ultimately reduces to the fundamental question of “when do nondiscrete entities become distinct, mutually exclusive evolutionary units”? Species are perceived as critical biological entities, and the discovery and naming of new species is perceived by many authors as a major research aim for assessing current biodiversity before much of it becomes extinct. However, less attention is given to determining whether these names represent valid biological entities because this is perceived as both a laborious chore and an undesirable research outcome. The charismatic spurge hawkmoths (Hyles euphorbiae complex, HEC) offer an opportunity to study this less fashionable aspect of systematics. To elucidate this intriguing systematic challenge, we analyzed over 10,000 ddRAD single nucleotide polymorphisms from 62 individuals using coalescent-based and population genomic methodology. These genome-wide data reveal a clear overestimation of (sub)species-level diversity and demonstrate that the HEC taxonomy has been seriously oversplit. We conclude that only one valid species name should be retained for the entire HEC, namely Hyles euphorbiae, and we do not recognize any formal subspecies or other taxonomic subdivisions within it. Although the adoption of genetic tools has frequently revealed morphologically cryptic diversity, the converse, taxonomic oversplitting of species, is generally (and wrongly in our opinion) accepted as rare. Furthermore, taxonomic oversplitting is most likely to have taken place in intensively studied popular and charismatic organisms such as the HEC.

Gene Flow in the Müllerian Mimicry Ring of a Poisonous Papuan Songbird Clade (Pitohui; Aves)

Müllerian mimicry rings are remarkable symbiotic species assemblages in which multiple members share a similar phenotype. However, their evolutionary origin remains poorly understood. Although gene flow among species has been shown to generate mimetic patterns in some Heliconius butterflies, mimicry is believed to be due to true convergence without gene flow in many other cases. We investigated the evolutionary history of multiple members of a passerine mimicry ring in the poisonous Papuan pitohuis. Previous phylogenetic evidence indicates that the aposematic coloration shared by many, but not all, members of this genus is ancestral and has only been retained by members of the mimicry ring. Using a newly assembled genome and thousands of genomic DNA markers, we demonstrate gene flow from the hooded pitohui (Pitohui dichrous) into the southern variable pitohui (Pitohui uropygialis), consistent with shared patterns of aposematic coloration. The vicinity of putatively introgressed loci is significantly enriched for genes that are important in melanin pigment expression and toxin resistance, suggesting that gene flow may have been instrumental in the sharing of plumage patterns and toxicity. These results indicate that interspecies gene flow may be a more general mechanism in generating mimicry rings than hitherto appreciated.