Interspecific transfer of genetic information through polyploid bridges.

Hybridization blurs species boundaries and leads to intertwined lineages resulting in reticulate evolution. Polyploidy, the outcome of whole genome duplication (WGD), has more recently been implicated in promoting and facilitating hybridization between polyploid species, potentially leading to adaptive introgression. However, because polyploid lineages are usually ephemeral states in the evolutionary history of life it is unclear whether WGD-potentiated hybridization has any appreciable effect on their diploid counterparts. Here, we develop a model of cytotype dynamics within mixed-ploidy populations to demonstrate that polyploidy can in fact serve as a bridge for gene flow between diploid lineages, where introgression is fully or partially hampered by the species barrier. Polyploid bridges emerge in the presence of triploid organisms, which despite critically low levels of fitness, can still allow the transfer of alleles between diploid states of independently evolving mixed-ploidy species. Notably, while marked genetic divergence prevents polyploid-mediated interspecific gene flow, we show that increased recombination rates can offset these evolutionary constraints, allowing a more efficient sorting of alleles at higher-ploidy levels before introgression into diploid gene pools. Additionally, we derive an analytical approximation for the rate of gene flow at the tetraploid level necessary to supersede introgression between diploids with nonzero introgression rates, which is especially relevant for plant species complexes, where interspecific gene flow is ubiquitous. Altogether, our results illustrate the potential impact of polyploid bridges on the (re)distribution of genetic material across ecological communities during evolution, representing a potential force behind reticulation.

Gene flow from Fraxinus cultivars into natural stands of Fraxinus pennsylvanica occurs range-wide, is regionally extensive, and is associated with a loss of allele richness.

In North America, a comparatively small number of Fraxinus (ash) cultivars were planted in large numbers in both urban and rural environments across the entire range of Fraxinus pennsylvanica Marsh (green ash) over the last 80 years. Undetected cultivar gene flow, if extensive, could significantly lower genetic diversity within populations, suppress differentiation between populations, generate interspecific admixture not driven by long-standing natural processes, and affect the impact of abiotic and biotic threats. In this investigation we generated the first range-wide genetic assessment of F. pennsylvanica to detect the extent of cultivar gene flow into natural stands. We used 16 EST-SSR markers to genotype 48 naturally regenerated populations of F. pennsylvanica distributed across the native range (1291 trees), 19 F. pennsylvanica cultivars, and one F. americana L. (white ash) cultivar to detect cultivar propagule dispersal into these populations. We detected first generation cultivar parentage with high confidence in 171 individuals in 34 of the 48 populations and extensive cultivar parentage (23-50%) in eight populations. The incidence of cultivar parentage was negatively associated with allele richness (R2 = 0.151, p = 0.006). The evidence for a locally high frequency of cultivar propagule dispersal and the interspecific admixture in eastern populations will inform Fraxinus gene pool conservation strategies and guide the selection of individuals for breeding programs focused on increasing resistance to the emerald ash borer (Agrilus planipennis Fairmaire), an existential threat to the Fraxinus species of North America.

Barriers and corridors: Assessment of gene flow and movement among red panda populations in eastern Himalayas.

Landscape features can impede dispersal, gene flow, and population demography, resulting in the formation of several meta-populations within a continuous landscape. Understanding a species' ability to overcome these barriers is critical for predicting genetic connectivity and population persistence, and implementing effective conservation strategies. In the present study, we conducted a fine-scale spatial genetic analysis to understand the contemporary gene flow within red panda populations in the Eastern Himalayas. Employing geometric aspects of reserve design, we delineated the critical core habitats for red pandas, which comprise 14.5 % of the landscape (12,189.75 Km2), with only a mere 443 Km2 falling within the protected areas. We identified corridors among the core habitats, which may be vital for the species' long-term genetic viability. Furthermore, we identified substantial landscape barriers, including Sela Pass in the western region, Siang river in the central region, and the Dibang river, Lohit river, along with Dihang, Dipher, and Kumjawng passes in the eastern region, which hinder gene flow. We suggest managing red panda populations through the creation of Community Conservation Reserves in the identified core habitats, following landscape-level management planning based on the core principles of geometric reserve design. This includes a specific emphasis on identified core habitats of red panda (CH-RP 5 and CH-RP 8) to facilitate corridors and implement meta-population dynamics. We propose the development of a comprehensive, long-term conservation and management plan for red pandas in the transboundary landscape, covering China, Nepal, and Bhutan.

The role of historical biogeography in shaping colour morph diversity in the common wall lizard.

The maintenance of polymorphisms often depends on multiple selective forces, but less is known on the role of stochastic or historical processes in maintaining variation. The common wall lizard (Podarcis muralis) is a colour polymorphic species in which local colour morph frequencies are thought to be modulated by natural and sexual selection. Here, we used genome-wide single-nucleotide polymorphism data to investigate the relationships between morph composition and population biogeography at a regional scale, by comparing morph composition with patterns of genetic variation of 54 populations sampled across the Pyrenees. We found that genetic divergence was explained by geographic distance but not by environmental features. Differences in morph composition were associated with genetic and environmental differentiation, as well as differences in sex ratio. Thus, variation in colour morph frequencies could have arisen via historical events and/or differences in the permeability to gene flow, possibly shaped by the complex topography and environment. In agreement with this hypothesis, colour morph diversity was positively correlated with genetic diversity and rates of gene flow and inversely correlated with the likelihood of the occurrence of bottlenecks. Concurrently, we did not find conclusive evidence for selection in the two colour loci. As an illustration of these effects, we observed that populations with higher proportions of the rarer yellow and yellow-orange morphs had higher genetic diversity. Our results suggest that processes involving a decay in overall genetic diversity, such as reduced gene flow and/or bottleneck events have an important role in shaping population-specific morph composition via non-selective processes.

Genetic barriers more than environmental associations explain Serratia marcescens population structure.

Bacterial species often comprise well-separated lineages, likely emerged and maintained by genetic isolation and/or ecological divergence. How these two evolutionary actors interact in the shaping of bacterial population structure is currently not fully understood. In this study, we investigate the genetic and ecological drivers underlying the evolution of Serratia marcescens, an opportunistic pathogen with high genomic flexibility and able to colonise diverse environments. Comparative genomic analyses reveal a population structure composed of five deeply-demarcated genetic clusters with open pan-genome but limited inter-cluster gene flow, partially explained by Restriction-Modification (R-M) systems incompatibility. Furthermore, a large-scale research on hundred-thousands metagenomic datasets reveals only a partial habitat separation of the clusters. Globally, two clusters only show a separate gene composition coherent with ecological adaptations. These results suggest that genetic isolation has preceded ecological adaptations in the shaping of the species diversity, an evolutionary scenario coherent with the Evolutionary Extended Synthesis.

The influence of gene flow on population viability in an isolated urban caracal population.

Wildlife populations are becoming increasingly fragmented by anthropogenic development. Small and isolated populations often face an elevated risk of extinction, in part due to inbreeding depression. Here, we examine the genomic consequences of urbanization in a caracal (Caracal caracal) population that has become isolated in the Cape Peninsula region of the City of Cape Town, South Africa, and is thought to number ~50 individuals. We document low levels of migration into the population over the past ~75 years, with an estimated rate of 1.3 effective migrants per generation. As a consequence of this isolation and small population size, levels of inbreeding are elevated in the contemporary Cape Peninsula population (mean FROH = 0.20). Inbreeding primarily manifests as long runs of homozygosity >10 Mb, consistent with the effects of isolation due to the rapid recent growth of Cape Town. To explore how reduced migration and elevated inbreeding may impact future population dynamics, we parameterized an eco-evolutionary simulation model. We find that if migration rates do not change in the future, the population is expected to decline, though with a low projected risk of extinction. However, if migration rates decline or anthropogenic mortality rates increase, the potential risk of extinction is greatly elevated. To avert a population decline, we suggest that translocating migrants into the Cape Peninsula to initiate a genetic rescue may be warranted in the near future. Our analysis highlights the utility of genomic datasets coupled with computational simulation models for investigating the influence of gene flow on population viability.

The genetic identity of the Vedda: A language isolate of South Asia.

Linguistic data from South Asia identified several language isolates in the subcontinent. The Vedda, an indigenous population of Sri Lanka, are the least studied amongst them. Therefore, to understand the initial peopling of Sri Lanka and the genetic affinity of the Vedda with other populations in Eurasia, we extensively studied the high-resolution autosomal and mitogenomes from the Vedda population of Sri Lanka. Our autosomal analyses suggest a close genetic link of Vedda with the tribal populations of India despite no evidence of close linguistic affinity, thus suggesting a deep genetic link of the Vedda with these populations. The mitogenomic analysis supports this association by pointing to an ancient link with Indian populations. We suggest that the Vedda population is a genetically drifted group with limited gene flow from neighbouring Sinhalese and Sri Lankan Tamil populations. Interestingly, the genetic ancestry sharing of Vedda neglects the isolation-by-distance model. Collectively, the demography of Sri Lanka is unique, where Sinhalese and Sri Lankan Tamil populations excessively admixed, whilst Vedda largely preserved their isolation and deep genetic association with India.

The Net of Life, a short story: Intricate patterns of gene flows across hundreds of extant genomes, all the way to LUCA.

Over the past quarter-century, the field of evolutionary biology has been transformed by the emergence of complete genome sequences and the conceptual framework known as the 'Net of Life.' This paradigm shift challenges traditional notions of evolution as a tree-like process, emphasizing the complex, interconnected network of gene flow that may blur the boundaries between distinct lineages. In this context, gene loss, rather than horizontal gene transfer, is the primary driver of gene content, with vertical inheritance playing a principal role. The 'Net of Life' not only impacts our understanding of genome evolution but also has profound implications for classification systems, the rapid appearance of new traits, and the spread of diseases. Here, we explore the core tenets of the 'Net of Life' and its implications for genome-scale phylogenetic divergence, providing a comprehensive framework for further investigations in evolutionary biology.

Temporal analysis of genetic diversity and gene flow in the threatened catfish Pseudoplatystoma magdaleniatum from a dammed neotropical river.

The striped catfish Pseudoplatystoma magdaleniatum is a large-sized migratory species from the north Andes region, endemic to Magdalena basin and one of the major fishery resources. Despite the estimated reduction of over 80% of the fisheries production of this species throughout the basin in recent decades, its population in the lower Magdalena-Cauca basin showed healthy genetics after molecular analyses. However, the current conservation status of this species and several habitat disturbances demand the re-evaluation of its population genetics to infer evolutionary risks and assess potential changes. This work analyzed a total of 164 samples from the Cauca River collected downstream the Ituango Dam between 2019-2021 using species-specific microsatellite markers to compare the genetic diversity and structure in samples collected between 2010-2014 from the lower Magdalena-Cauca basin, previously analyzed. Our results showed a relatively stable panmictic population over time (4 to 10 years), with high genetic diversity and evidence of recent bottleneck. Promoting habitat connectivity to conserve gene flow, characterizing diversity and genetic structure over the entire basin, and integrating the results with future monitoring are important aspects for the management planning for P. magdaleniatum in the Magdalena-Cauca basin.

DNA barcoding, phylogeography and evolutionary dynamics of Chrysichthys auratus.

This study embarked on an exploration into the genetic structure and evolutionary history of the Chrysichthys auratus species, leveraging PCR amplification, phylogenetic trees, and haplotype networks. Specific DNA segments were successfully amplified and visualized through electrophoresis. Newly obtained sequences were Bank into GenBank and given accession numbers (OR730807-OR730808-OR730809). The Neighbor-Joining method provided insights into the evolutionary relationships among taxa, further augmented by bootstrap values and the Tamura 3-parameter method. A comprehensive geographical haplotype network showcased pronounced genetic differentiation, especially between remote populations. Nonetheless, shared haplotypes between proximate regions indicated either ancestral genetic connections or ongoing gene flow. Employing the COI-DNA barcodes, an in-depth understanding of intra- and inter-populational genetic diversity was achieved. The study's findings unravel the intricate genetic landscape and evolutionary dynamics of C. auratus, offering novel perspectives into its demographic history across its vast native habitat.

Global assessment of effective population sizes: Consistent taxonomic differences in meeting the 50/500 rule.

Effective population size (Ne) is a particularly useful metric for conservation as it affects genetic drift, inbreeding and adaptive potential within populations. Current guidelines recommend a minimum Ne of 50 and 500 to avoid short-term inbreeding and to preserve long-term adaptive potential respectively. However, the extent to which wild populations reach these thresholds globally has not been investigated, nor has the relationship between Ne and human activities. Through a quantitative review, we generated a dataset with 4610 georeferenced Ne estimates from 3829 populations, extracted from 723 articles. These data show that certain taxonomic groups are less likely to meet 50/500 thresholds and are disproportionately impacted by human activities; plant, mammal and amphibian populations had a <54% probability of reaching N ̂ e = 50 and a <9% probability of reaching N ̂ e = 500. Populations listed as being of conservation concern according to the IUCN Red List had a smaller median N ̂ e than unlisted populations, and this was consistent across all taxonomic groups. N ̂ e was reduced in areas with a greater Global Human Footprint, especially for amphibians, birds and mammals, however relationships varied between taxa. We also highlight several considerations for future works, including the role that gene flow and subpopulation structure plays in the estimation of N ̂ e in wild populations, and the need for finer-scale taxonomic analyses. Our findings provide guidance for more specific thresholds based on Ne and help prioritise assessment of populations from taxa most at risk of failing to meet conservation thresholds.

Recurrent gene flow events occurred during the diversification of clownfishes of the skunk complex.

Clownfish (subfamily Amphiprioninae) are an iconic group of coral reef fish that evolved a mutualistic interaction with sea anemones, which triggered the adaptive radiation of the clade. Within clownfishes, the "skunk complex" is particularly interesting. Besides ecological speciation, interspecific gene flow and hybrid speciation are thought to have shaped the evolution of the group. We investigated the mechanisms characterizing the diversification of this complex. By taking advantage of their disjunct geographical distribution, we obtained whole-genome data of sympatric and allopatric populations of the three main species of the complex (Amphiprion akallopisos, A. perideraion and A. sandaracinos). We examined population structure, genomic divergence and introgression signals and performed demographic modelling to identify the most realistic diversification scenario. We excluded scenarios of strict isolation or hybrid origin of A. sandaracinos. We discovered moderate gene flow from A. perideraion to the ancestor of A. akallopisos + A. sandaracinos and weak gene flow between the species in the Indo-Australian Archipelago throughout the diversification of the group. We identified introgressed regions in A. sandaracinos and detected in A. perideraion two large regions of high divergence from the two other species. While we found that gene flow has occurred throughout the species' diversification, we also observed that recent admixture was less pervasive than initially thought, suggesting a role of host repartition or behavioural barriers in maintaining the genetic identity of the species in sympatry.

Colonial-driven extinction of the blue antelope despite genomic adaptation to low population size.

Low genomic diversity is generally indicative of small population size and is considered detrimental by decreasing long-term adaptability.1,2,3,4,5,6 Moreover, small population size may promote gene flow with congeners and outbreeding depression.7,8,9,10,11,12,13 Here, we examine the connection between habitat availability, effective population size (Ne), and extinction by generating a 40× nuclear genome from the extinct blue antelope (Hippotragus leucophaeus). Historically endemic to the relatively small Cape Floristic Region in southernmost Africa,14,15 populations were thought to have expanded and contracted across glacial-interglacial cycles, tracking suitable habitat.16,17,18 However, we found long-term low Ne, unaffected by glacial cycles, suggesting persistence with low genomic diversity for many millennia prior to extinction in ∼AD 1800. A lack of inbreeding, alongside high levels of genetic purging, suggests adaptation to this long-term low Ne and that human impacts during the colonial era (e.g., hunting and landscape transformation), rather than longer-term ecological processes, were central to its extinction. Phylogenomic analyses uncovered gene flow between roan (H. equinus) and blue antelope, as well as between roan and sable antelope (H. niger), approximately at the time of divergence of blue and sable antelope (∼1.9 Ma). Finally, we identified the LYST and ASIP genes as candidates for the eponymous bluish pelt color of the blue antelope. Our results revise numerous aspects of our understanding of the interplay between genomic diversity and evolutionary history and provide the resources for uncovering the genetic basis of this extinct species' unique traits.

The Andes as a semi-permeable geographical barrier: Genetic connectivity between structured populations in a widespread spider.

Geographical barriers like mountain ranges impede genetic exchange among populations, promoting diversification. The effectiveness of these barriers in limiting gene flow varies between lineages due to each species' dispersal modes and capacities. Our understanding of how the Andes orogeny contributes to species diversification comes from well-studied vertebrates and a few arthropods and plants, neglecting organisms unable to fly or walk long distances. Some arachnids, such as Gasteracantha cancriformis, have been hypothesized to disperse long distances via ballooning (i.e. using their silk to interact with the wind). Yet, we do not know how the environment and geography shape its genetic diversity. Therefore, we tested whether the Andes contributed to the diversification of G. cancriformis acting as an absolute or semi-permeable barrier to genetic connectivity between populations of this spider at opposite sides of the mountain range. We sampled thousands of loci across the distribution of the species and implemented population genetics, phylogenetic, and landscape genetic analyses. We identified two genetically distinct groups structured by the Central Andes, and a third less structured group in the Northern Andes that shares ancestry with the previous two. This structure is largely explained by the altitude along the Andes, which decreases in some regions, possibly facilitating cross-Andean dispersal and gene flow. Our findings support that altitude in the Andes plays a major role in structuring populations in South America, but the strength of this barrier can be overcome by organisms with long-distance dispersal modes together with altitudinal depressions.

Las barreras geográficas como las cordilleras montañosas impiden el intercambio genético entre poblaciones, promoviendo la diversificación. La efectividad de estas barreras para limitar el flujo genético varía entre linajes debido a los modos y capacidades de dispersión de cada especie. Nuestra comprensión de cómo la orogenia de los Andes contribuye a la diversificación de especies proviene de vertebrados y algunos artrópodos y plantas bien estudiados, descuidando a los organismos incapaces de volar o caminar grandes distancias. Se ha hipotetizado que algunas arañas, como Gasteracantha cancriformis, se dispersan a grandes distancias mediante la técnica de "ballooning" (es decir, utilizando su seda para interactuar con el viento). Sin embargo, no sabemos cómo el entorno y la geografía han dado forma a su diversidad genética. Por lo tanto, probamos si los Andes contribuyeron a la diversificación de G. cancriformis actuando como una barrera absoluta o permeable para la conectividad genética entre poblaciones de esta araña en lados opuestos de la cordillera. Muestreamos miles de loci a través de la distribución de la especie e implementamos análisis de genética de poblaciones, filogenéticos y de genética del paisaje. Identificamos dos grupos genéticamente distintos estructurados por los Andes Centrales, y un tercer grupo menos estructurado en los Andes del Norte que comparte ascendencia con los dos anteriores. Esta estructura se explica en gran medida por la altitud a lo largo de los Andes, que disminuye en algunas regiones, posiblemente facilitando la dispersión y el flujo genético a través de los Andes. Nuestros hallazgos apoyan que la altitud en los Andes juega un papel importante en la estructuración de las poblaciones en América del Sur, pero la fuerza de esta barrera puede ser superada por organismos con modos de dispersión a larga distancia junto con depresiones altitudinales.

Effective dispersal patterns in prairie plant species across human-modified landscapes.

Effective dispersal among plant populations is dependent on vector behaviour, landscape features and availability of adequate habitats. To capture landscape feature effects on dispersal, studies must be conducted at scales reflecting single-generation dispersal events (mesoscale). Many studies are conducted at large scales where genetic differentiation is due to dispersal occurring over multiple generations, making it difficult to interpret the effects of specific landscape features on vector behaviour. Genetic structure at the mesoscale may be determined by ecological and evolutionary processes, such as the consequences of vector behaviour on patterns of gene flow. We used chloroplast haplotypes and nuclear genome SNP surveys to identify landscape features influencing seed and pollen dispersal at a mesoscale within the Rogue River Valley in southern Oregon. We evaluated biotic and abiotic vector behaviour by contrasting two annual species with differing dispersal mechanisms; Achyrachaena mollis (Asteraceae) is a self-pollinating and anemochoric species, and Plectritis congesta (Caprifoliaceae) is biotically pollinated with barochoric seeds. Using landscape genetics methods, we identified features of the study region that conduct or restrict dispersal. We found chloroplast haplotypes were indicative of historic patterns of gene flow prior to human modification of landscapes. Seed dispersal of A. mollis was best supported by models of isolation by distance, while seed-driven gene flow of P. congesta was determined by the distribution of preserved natural spaces and quality habitat. Nuclear genetic structure was driven by both pollen and seed dispersal, and both species responded to contemporary landscape changes, such as urban and agricultural conversion, and habitat availability.

Combining population genomics and biophysical modelling to assess connectivity patterns in an Antarctic fish.

Connectivity is a fundamental process of population dynamics in marine ecosystems. In the last decade, with the emergence of new methods, combining different approaches to understand the patterns of connectivity among populations and their regulation has become increasingly feasible. The Western Antarctic Peninsula (WAP) is characterized by complex oceanographic dynamics, where local conditions could act as barriers to population connectivity. Here, the notothenioid fish Harpagifer antarcticus, a demersal species with a complex life cycle (adults with poor swim capabilities and pelagic larvae), was used to assess connectivity along the WAP by combining biophysical modelling and population genomics methods. Both approaches showed congruent patterns. Areas of larvae retention and low potential connectivity, observed in the biophysical model output, coincide with four genetic groups within the WAP: (1) South Shetland Islands, (2) Bransfield Strait, (3) the central and (4) the southern area of WAP (Marguerite Bay). These genetic groups exhibited limited gene flow between them, consistent with local oceanographic conditions, which would represent barriers to larval dispersal. The joint effect of geographic distance and larval dispersal by ocean currents had a greater influence on the observed population structure than each variable evaluated separately. The combined effect of geographic distance and a complex oceanographic dynamic would be generating limited levels of population connectivity in the fish H. antarcticus along the WAP. Based on this, population connectivity estimations and priority areas for conservation were discussed, considering the marine protected area proposed for this threatened region of the Southern Ocean.

Introgression from local cultivars is a driver of agricultural adaptation in Argentinian weedy rice.

Weedy rice, a pervasive and troublesome weed found across the globe, has often evolved through fertilization of rice cultivars with little importance of crop-weed gene flow. In Argentina, weedy rice has been reported as an important constraint since the early 1970s, and, in the last few years, strains with herbicide-resistance are suspected to evolve. Despite their importance, the origin and genetic composition of Argentinian weedy rice as well its adaptation to agricultural environments has not been explored so far. To study this, we conducted genotyping-by-sequencing on samples of Argentinian weedy and cultivated rice and compared them with published data from weedy, cultivated and wild rice accessions distributed worldwide. In addition, we conducted a phenotypic characterization for weedy-related traits, a herbicide resistance screening and genotyped accessions for known mutations in the acetolactate synthase (ALS) gene, which confers herbicide resistance. Our results revealed large phenotypic variability in Argentinian weedy rice. Most strains were resistant to ALS-inhibiting herbicides with a high frequency of the ALS mutation (A122T) present in Argentinian rice cultivars. Argentinian cultivars belonged to the three major genetic groups of rice: japonica, indica and aus while weeds were mostly aus or aus-indica admixed, resembling weedy rice strains from the Southern Cone region. Phylogenetic analysis supports a single origin for aus-like South American weeds, likely as seed contaminants from the United States, and then admixture with local indica cultivars. Our findings demonstrate that crop to weed introgression can facilitate rapid adaptation to agriculture environments.

Synthesizing environmental, epidemiological and vector and parasite genetic data to assist decision making for disease elimination.

We present a framework for identifying when conditions are favourable for transmission of vector-borne diseases between communities by incorporating predicted disease prevalence mapping with landscape analysis of sociological, environmental and host/parasite genetic data. We explored the relationship between environmental features and gene flow of a filarial parasite of humans, Onchocerca volvulus, and its vector, blackflies in the genus Simulium. We generated a baseline microfilarial prevalence map from point estimates from 47 locations in the ecological transition separating the savannah and forest in Ghana, where transmission of O. volvulus persists despite onchocerciasis control efforts. We generated movement suitability maps based on environmental correlates with mitochondrial population structure of 164 parasites from 15 communities and 93 vectors from only four sampling sites, and compared these to the baseline prevalence map. Parasite genetic distance between sampling locations was significantly associated with elevation (r = .793, p = .005) and soil moisture (r = .507, p = .002), while vector genetic distance was associated with soil moisture (r = .788, p = .0417) and precipitation (r = .835, p = .0417). The correlation between baseline prevalence and parasite resistance surface maps was stronger than that between prevalence and vector resistance surface maps. The centre of the study area had high prevalence and suitability for parasite and vector gene flow, potentially contributing to persistent transmission and suggesting the importance of re-evaluating transmission zone boundaries. With suitably dense sampling, this framework can help delineate transmission zones for onchocerciasis and would be translatable to other vector-borne diseases.

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.

Study on the origin of the Baise horse based on whole-genome resequencing.

The Baise horse, an indigenous horse breed mainly distributed in the Baise region of Guangxi province in southwest China, has a long history as draft animal. However, there is a lack of research regarding the origin and ancestral composition of the Baise horse. In this study, whole-genome resequencing data from 236 horses of seven Chinese indigenous horse breeds, five foreign horse breeds, and four Przewalski's horses were used to investigate the relationships between the Baise horse and other horse breeds. The results showed that foreign horse breeds had no significant impact on the formation of the Baise horse. The two southwestern horse populations, the Debao pony and the Jinjiang horse, exhibit the closest genetic affinity with the Baise horse. This is consistent with their adjacent geographical distribution. Analysis of the migration route revealed a gene flow from the Chakouyi horse into the Baise horse. In summary, our results confirm that the formation of the Baise horse did not involve participation from foreign breeds. Geographical distance emerges as a crucial factor in determining the genetic relationships with the Baise horse. Gene flows of indigenous horse breeds along ancient routes of trade activities had played a role in the formation of the Baise horse.