Dsuite - Fast D-statistics and related admixture evidence from VCF files

Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini)

The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.

Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China

Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.

Disentangling sources of gene tree discordance for Hordeum species via target-enriched sequencing assays

Hordeum is an economically and evolutionarily important genus within the Triticeae tribe of the family Poaceae, and contains 33 widely distributed and diverse species which cytologically represent four subgenomes (H, Xa, Xu and I). These wild species (except Hordeum spontaneum, which is the primary gene pool of barley) are secondary or tertiary gene-pool germplasms for barley and wheat improvement, and uncovering their complicated evolutionary relationships would benefit for future breeding programs. Here, we developed a complexity-reduced pipeline via capturing genome-wide distributed fragments via two novel target-enriched assays (HorCap v1.0 and BarPlex v1.0) in conjugation with high-throughput sequencing of the enrichments. Both assays were tested for genotyping 40 species from three genera (Hordeum, Triticum, and Aegilops) containing 82 samples 67 accessions. Either of both assays worked efficiently in genotyping, while integration of both assays can significantly improve the robustness and resolution of the Hordeum phylogenetic trees. Interestingly, the incomplete lineage sorting (ILS) was inferred for the first time as the major factor causing phylogenetic discordance among the four subgenomes, whereas in New World species (carrying I genome) post-speciation introgression events were revealed. Through revising the evolutionary relationships of the Hordeum species based on an ancestral state reconstruction for the diploids and parental donor inference for the polyploids, our results raised new queries about the Hordeum phylogeny. Moreover, both newly-developed assays are applicable in genotyping and phylogenetic analysis of Hordeum and other Triticeae wild species.

Museum genomics approach to study the taxonomy and evolution of Woolly-necked storks using historic specimens

The accessibility of genomic tools in evolutionary biology has allowed for a thorough exploration of various evolutionary processes associated with adaptation and speciation. However, genomic studies in natural systems present numerous challenges, reflecting the inherent complexities of studying organisms in their native habitats. The utilization of museum specimens for genomics research has received increased attention in recent times, facilitated by advancements in ancient DNA techniques. In this study, we have utilized a museum genomics approach to analyze historic specimens of Woolly-necked storks (Ciconia spp.) and examine their genetic composition and taxonomic status and explore the evolutionary and adaptive trajectories of populations over the years. The Woolly-necked storks are distributed in Asia and Africa with a taxonomic classification that has been a matter of ambiguity. Asian and African Woollynecks were recently recognized as different species based on their morphological differences; however, their genomic validation was lacking. In this study, we have used ∼70-year-old museum samples for whole-genome population-scale sequencing. Our study has revealed that Asian and African Woollynecks are genetically distinct, consistent with the current taxonomic classification based on morphological features. However, we also found a high genetic divergence between the Asian subspecies Ciconia episcopus neglecta and Ciconia episcopus episcopus, suggesting this classification requires a detailed examination to explore processes of ongoing speciation. Because taxonomic classification directly impacts conservation efforts, and there is evidence of declining populations of Asian Woollynecks in Southeast Asia, our results highlight that population-scale studies are urgent to determine the genetic, ecological, and phylogenetic diversity of these birds.

Genome-wide data reveal bi-direction and asymmetrical hybridization origin of a fern species Microlepia matthewii

Introduction

Natural hybridization is common and plays a crucial role in driving biodiversity in nature. Despite its significance, the understanding of hybridization in ferns remains inadequate. Therefore, it is imperative to study fern hybridization to gain a more comprehensive understanding of fern biodiversity. Our study delves into the role of hybridization in shaping fern species, employing Microlepia matthewii as a case study to investigate its origins of hybridization.

Methods

We performed double digest Genotyping-by-sequencing (dd-GBS) on M. matthewii and its potential parent species, identifying nuclear and chloroplast SNPs. Initially, nuclear SNPs were employed to construct the three cluster analysis: phylogenetic tree, principal component analysis, and population structure analysis. Subsequently, to confirm whether the observed genetic mixture pattern resulted from hybridization, we utilized two methods: ABBA-BABA statistical values in the D-suite program and gene frequency covariance in the Treemix software to detect gene flow. Finally, we employed chloroplast SNPs to construct a phylogenetic tree, tracing the maternal origin.

Results and discussion

The analysis of the nuclear SNP cluster revealed that M. matthewii possesses a genetic composition that is a combination of M. hancei and M. calvescens. Furthermore, the analysis provided strong evidence of significant gene flow signatures from the parental species to the hybrid, as indicated by the two gene flow analyses. The samples of M. matthewii cluster separately with M. hancei or M. calvescens on the chloroplast systematic tree. However, the parentage ratio significantly differs from 1:1, suggesting that M. matthewii is a bidirectional and asymmetrical hybrid offspring of M. hancei and M. calvescens.

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

Ancient and Recent Hybridization in the Oreochromis Cichlid Fishes

Cichlid fishes of the genus Oreochromis (tilapia) are among the most important fish for inland capture fisheries and global aquaculture. Deliberate introductions of non-native species for fisheries improvement and accidental escapees from farms have resulted in admixture with indigenous species. Such hybridization may be detrimental to native biodiversity, potentially leading to genomic homogenization of populations and the loss of important genetic material associated with local adaptation. By contrast, introgression may fuel diversification when combined with ecological opportunity, by supplying novel genetic combinations. To date, the role of introgression in the evolutionary history of tilapia has not been explored. Here we studied both ancient and recent hybridization in tilapia, using whole genome resequencing of 575 individuals from 23 species. We focused on Tanzania, a natural hotspot of tilapia diversity, and a country where hybridization between exotic and native species in the natural environment has been previously reported. We reconstruct the first genome-scale phylogeny of the genus and reveal prevalent ancient gene flow across the Oreochromis phylogeny. This has likely resulted in the hybrid speciation of one species, O. chungruruensis. We identify multiple cases of recent hybridization between native and introduced species in the wild, linked to the use of non-native species in both capture fisheries improvement and aquaculture. This has potential implications for both conservation of wild populations and the development of the global tilapia aquaculture industry.

Inference of Locus-Specific Population Mixtures from Linked Genome-Wide Allele Frequencies

Admixture between populations and species is common in nature. Since the influx of new genetic material might be either facilitated or hindered by selection, variation in mixture proportions along the genome is expected in organisms undergoing recombination. Various graph-based models have been developed to better understand these evolutionary dynamics of population splits and mixtures. However, current models assume a single mixture rate for the entire genome and do not explicitly account for linkage. Here, we introduce TreeSwirl, a novel method for inferring branch lengths and locus-specific mixture proportions by using genome-wide allele frequency data, assuming that the admixture graph is known or has been inferred. TreeSwirl builds upon TreeMix that uses Gaussian processes to estimate the presence of gene flow between diverged populations. However, in contrast to TreeMix, our model infers locus-specific mixture proportions employing a hidden Markov model that accounts for linkage. Through simulated data, we demonstrate that TreeSwirl can accurately estimate locus-specific mixture proportions and handle complex demographic scenarios. It also outperforms related D- and f-statistics in terms of accuracy and sensitivity to detect introgressed loci.

Sirenian genomes illuminate the evolution of fully aquatic species within the mammalian superorder afrotheria

Sirenians of the superorder Afrotheria were the first mammals to transition from land to water and are the only herbivorous marine mammals. Here, we generated a chromosome-level dugong (Dugong dugon) genome. A comparison of our assembly with other afrotherian genomes reveals possible molecular adaptations to aquatic life by sirenians, including a shift in daily activity patterns (circadian clock) and tolerance to a high-iodine plant diet mediated through changes in the iodide transporter NIS (SLC5A5) and its co-transporters. Functional in vitro assays confirm that sirenian amino acid substitutions alter the properties of the circadian clock protein PER2 and NIS. Sirenians show evidence of convergent regression of integumentary system (skin and its appendages) genes with cetaceans. Our analysis also uncovers gene losses that may be maladaptive in a modern environment, including a candidate gene (KCNK18) for sirenian cold stress syndrome likely lost during their evolutionary shift in daily activity patterns. Genomes from nine Australian locations and the functionally extinct Okinawan population confirm and date a genetic break ~10.7 thousand years ago on the Australian east coast and provide evidence of an associated ecotype, and highlight the need for whole-genome resequencing data from dugong populations worldwide for conservation and genetic management.

Hybridization in palms (Arecaceae)

Hybridization has significant evolutionary consequences across the Tree of Life. The process of hybridization has played a major role in plant evolution and has contributed to species richness and trait variation. Since morphological traits are partially a product of their environment, there may be a link between hybridization and ecology. Plant hybrid species richness is noted to be higher in harsh environments, and we explore this hypothesis with a keystone tropical plant lineage, palms (Arecaceae). Leveraging a recent literature review of naturally occurring palm hybrids, we developed a method to calculate hybrid frequency, and then tested if there is phylogenetic signal of hybrids using a phylogeny of all palms. Further, we used phylogenetic comparative methods to examine the interaction between hybrid frequency and presence in dry environments, on islands, and the species richness of genera. Phylogenetic generalized least squares models had stronger support than models of random association, indicating phylogenetic signal for the presence of hybrids in dry and island environments. However, all p-values were >.05 and therefore the correlation was poor between hybridization and the trait frequencies examined. Presence in particular environments are not strongly correlated to hybrid frequency, but phylogenetic signal suggests a role in its distribution in different habitats. Hybridization in palms is not evenly distributed across subfamilies, tribes, subtribes yet plays an important role in palm diversity, nonetheless. Increasing our understanding hybridization in this economically and culturally important plant family is essential, particularly since rates are projected to increase with climate change, reconfiguring the dynamics and distribution of biodiversity.

Broken Ring Speciation in California Mygalomorph Spiders (Nemesiidae, Calisoga)

AbstractIdealized ring species, with approximately continuous gene flow around a geographic barrier but singular reproductive isolation at a ring terminus, are rare in nature. A broken ring species model preserves the geographic setting and fundamental features of an idealized model but accommodates varying degrees of gene flow restriction over complex landscapes through evolutionary time. Here we examine broken ring species dynamics in Calisoga spiders, which, like the classic ring species Ensatina salamanders, are distributed around the Central Valley of California. Using nuclear and mitogenomic data, we test key predictions of common ancestry, ringlike biogeography, biogeographic timing, population connectivity, and terminal overlap. We show that a ring complex of populations shares a single common ancestor, and from an ancestral area in the Sierra Nevada mountains, two distributional and phylogenomic arms encircle the Central Valley. Isolation by distance occurs along these distributional arms, although gene flow restriction is also evident. Where divergent lineages meet in the South Coast Ranges, we find rare lineage sympatry, without evidence for nuclear gene flow and with clear evidence for morphological and ecological divergence. We discuss general insights provided by broken ring species and how such a model could be explored and extended in other systems and future studies.

Population genomics highlights structural variations in local adaptation to saline coastal environments in woolly grape

Structural variations (SVs) are a feature of plant genomes that has been largely unexplored despite their significant impact on plant phenotypic traits and local adaptation to abiotic and biotic stress. In this study, we employed woolly grape (Vitis retordii), a species native to the tropical and subtropical regions of East Asia with both coastal and inland habitats, as a valuable model for examining the impact of SVs on local adaptation. We assembled a haplotype-resolved chromosomal reference genome for woolly grape, and conducted population genetic analyses based on whole-genome sequencing (WGS) data from coastal and inland populations. The demographic analyses revealed recent bottlenecks in all populations and asymmetric gene flow from the inland to the coastal population. In total, 1,035 genes associated with plant adaptive regulation for salt stress, radiation, and environmental adaptation were detected underlying local selection by SVs and SNPs in the coastal population, of which 37.29% and 65.26% were detected by SVs and SNPs, respectively. Candidate genes such as FSD2, RGA1, and AAP8 associated with salt tolerance were found to be highly differentiated and selected during the process of local adaptation to coastal habitats in SV regions. Our study highlights the importance of SVs in local adaptation; candidate genes related to salt stress and climatic adaptation to tropical and subtropical environments are important genomic resources for future breeding programs of grapevine and its rootstocks.

Long-term climatic stability drives accumulation and maintenance of divergent freshwater fish lineages in a temperate biodiversity hotspot

Anthropogenic climate change is forecast to drive regional climate disruption and instability across the globe. These impacts are likely to be exacerbated within biodiversity hotspots, both due to the greater potential for species loss but also to the possibility that endemic lineages might not have experienced significant climatic variation in the past, limiting their evolutionary potential to respond to rapid climate change. We assessed the role of climatic stability on the accumulation and persistence of lineages in an obligate freshwater fish group endemic to the southwest Western Australia (SWWA) biodiversity hotspot. Using 19,426 genomic (ddRAD-seq) markers and species distribution modelling, we explored the phylogeographic history of western (Nannoperca vittata) and little (Nannoperca pygmaea) pygmy perches, assessing population divergence and phylogenetic relationships, delimiting species and estimating changes in species distributions from the Pliocene to 2100. We identified two deep phylogroups comprising three divergent clusters, which showed no historical connectivity since the Pliocene. We conservatively suggest these represent three isolated species with additional intraspecific structure within one widespread species. All lineages showed long-term patterns of isolation and persistence owing to climatic stability but with significant range contractions likely under future climate change. Our results highlighted the role of climatic stability in allowing the persistence of isolated lineages in the SWWA. This biodiversity hotspot is under compounding threat from ongoing climate change and habitat modification, which may further threaten previously undetected cryptic diversity across the region.

Long-term balancing selection for pathogen resistance maintains trans-species polymorphisms in a planktonic crustacean

Balancing selection is an evolutionary process that maintains genetic polymorphisms at selected loci and strongly reduces the likelihood of allele fixation. When allelic polymorphisms that predate speciation events are maintained independently in the resulting lineages, a pattern of trans-species polymorphisms may occur. Trans-species polymorphisms have been identified for loci related to mating systems and the MHC, but they are generally rare. Trans-species polymorphisms in disease loci are believed to be a consequence of long-term host-parasite coevolution by balancing selection, the so-called Red Queen dynamics. Here we scan the genomes of three crustaceans with a divergence of over 15 million years and identify 11 genes containing identical-by-descent trans-species polymorphisms with the same polymorphisms in all three species. Four of these genes display molecular footprints of balancing selection and have a function related to immunity. Three of them are located in or close to loci involved in resistance to a virulent bacterial pathogen, Pasteuria, with which the Daphnia host is known to coevolve. This provides rare evidence of trans-species polymorphisms for loci known to be functionally relevant in interactions with a widespread and highly specific parasite. These findings support the theory that specific antagonistic coevolution is able to maintain genetic diversity over millions of years.

Evidence for gene flow and trait reversal during radiation of Mexican Goodeid fish

Understanding the phylogeographic history of a group and identifying the factors contributing to speciation is an important challenge in evolutionary biology. The Goodeinae are a group of live-bearing fishes endemic to Mexico. Here, we develop genomic resources for species within the Goodeinae and use phylogenomic approaches to characterise their evolutionary history. We sequenced, assembled and annotated the genomes of four Goodeinae species, including Ataeniobius toweri, the only matrotrophic live-bearing fish without a trophotaenia in the group. We estimated timings of species divergence and examined the extent and timing of introgression between the species to assess if this may have occurred during an early radiation, or in more recent episodes of secondary contact. We used branch-site models to detect genome-wide positive selection across Goodeinae, and we specifically asked whether this differs in A. toweri, where loss of placental viviparity has recently occurred. We found evidence of gene flow between geographically isolated species, suggesting vicariant speciation was supplemented by limited post-speciation gene flow, and gene flow may explain previous uncertainties about Goodeid phylogeny. Genes under positive selection in the group are likely to be associated with the switch to live-bearing. Overall, our studies suggest that both volcanism-driven vicariance and changes in reproductive mode influenced radiation in the Goodeinae.

Genomics of hybrid parallel origin in Aquilegia ecalcarata

BACKGROUND

The parallel evolution of similar traits or species provides strong evidence for the role of natural selection in evolution. Traits or species that evolved repeatedly can be driven by separate de novo mutations or interspecific gene flow. Although parallel evolution has been reported in many studies, documented cases of parallel evolution caused by gene flow are scarce by comparison. Aquilegia ecalcarata and A. kansuensis belong to the genus of Aquilegia, and are the closest related sister species. Mutiple origins of A. ecalcarata have been reported in previous studies, but whether they have been driven by separate de novo mutations or gene flow remains unclear.

RESULTS

In this study, We conducted genomic analysis from 158 individuals of two repeatedly evolving pairs of A. ecalcarata and A. kansuensis. All samples were divided into two distinct clades with obvious geographical distribution based on phylogeny and population structure. Demographic modeling revealed that the origin of the A. ecalcarata in the Eastern of China was caused by gene flow, and the Eastern A. ecalcarata occurred following introgression from Western A. ecalcarata population. Analysis of Treemix and D-statistic also revealed that a strong signal of gene flow was detected from Western A. ecalcarata to Eastern A. ecalcarata. Genetic divergence and selective sweep analyses inferred parallel regions of genomic divergence and identified many candidate genes associated with ecologically adaptive divergence between species pair. Comparative analysis of parallel diverged regions and gene introgression confirms that gene flow contributed to the parallel evolution of A. ecalcarata.

CONCLUSIONS

Our results further confirmed the multiple origins of A. ecalcarata and highlighted the roles of gene flow. These findings provide new evidence for parallel origin after hybridization as well as insights into the ecological adaptation mechanisms underlying the parallel origins of species.

Evolutionary genomics of climatic adaptation and resilience to climate change in alfalfa

Given the escalating impact of climate change on agriculture and food security, gaining insights into the evolutionary dynamics of climatic adaptation and uncovering climate-adapted variation can empower the breeding of climate-resilient crops to face future climate change. Alfalfa (Medicago sativa subsp. sativa), the queen of forages, shows remarkable adaptability across diverse global environments, making it an excellent model for investigating species responses to climate change. In this study, we performed population genomic analyses using genome resequencing data from 702 accessions of 24 Medicago species to unravel alfalfa's climatic adaptation and genetic susceptibility to future climate change. We found that interspecific genetic exchange has contributed to the gene pool of alfalfa, particularly enriching defense and stress-response genes. Intersubspecific introgression between M. sativa subsp. falcata (subsp. falcata) and alfalfa not only aids alfalfa's climatic adaptation but also introduces genetic burden. A total of 1671 genes were associated with climatic adaptation, and 5.7% of them were introgressions from subsp. falcata. By integrating climate-associated variants and climate data, we identified populations that are vulnerable to future climate change, particularly in higher latitudes of the Northern Hemisphere. These findings serve as a clarion call for targeted conservation initiatives and breeding efforts. We also identified pre-adaptive populations that demonstrate heightened resilience to climate fluctuations, illuminating a pathway for future breeding strategies. Collectively, this study enhances our understanding about the local adaptation mechanisms of alfalfa and facilitates the breeding of climate-resilient alfalfa cultivars, contributing to effective agricultural strategies for facing future climate change.

Integrative taxonomic analyses reveal that rapid genetic divergence drives Acropora speciation

Reef-building corals provide the structural basis for one of Earth's most spectacular and diverse but increasingly threatened ecosystems. The reef-building coral genus Acropora may have undergone substantial speciation during the Pleistocene climate and sea-level changes. Here, we aimed to evaluate the speciation history of four morphologically similar tabular Acropora species (Acropora aff. hyacinthus, A. cf. bifurcata, A. cf. cytherea, and A. cf. subulata) using an integrative approach with morphology, genetic, and reproduction methodology. Extensive morphological analyses showed that these four species are distinct and exhibited high gamete incompatibility, preventing hybridization. Furthermore, population structure and principal component analyses with SNPs (>60,000) indicated that these species were genetically distinct, and the ABBA-BABA test did not support introgression among these species. Many of their coding and noncoding RNA sequences showed high genetic variance at loci with high Fst values along the genome. Comparison of these orthologs with those of other Acropora species suggested that many of these genes are under positive selection, which could be associated with spawning time, gamete, and morphological divergence. Our findings show that the speciation of tabular Acropora occurred without hybridization, and the divergence accompanying the rapid evolution of genes in species-rich Acropora could be associated with speciation.

Comparative population genomics reveals convergent and divergent selection in the apricot-peach-plum-mei complex

The economically significant genus Prunus includes fruit and nut crops that have been domesticated for shared and specific agronomic traits; however, the genomic signals of convergent and divergent selection have not been elucidated. In this study, we aimed to detect genomic signatures of convergent and divergent selection by conducting comparative population genomic analyses of the apricot-peach-plum-mei (APPM) complex, utilizing a haplotype-resolved telomere-to-telomere (T2T) genome assembly and population resequencing data. The haplotype-resolved T2T reference genome for the plum cultivar was assembled through HiFi and Hi-C reads, resulting in two haplotypes 251.25 and 251.29 Mb in size, respectively. Comparative genomics reveals a chromosomal translocation of ~1.17 Mb in the apricot genomes compared with peach, plum, and mei. Notably, the translocation involves the D locus, significantly impacting titratable acidity (TA), pH, and sugar content. Population genetic analysis detected substantial gene flow between plum and apricot, with introgression regions enriched in post-embryonic development and pollen germination processes. Comparative population genetic analyses revealed convergent selection for stress tolerance, flower development, and fruit ripening, along with divergent selection shaping specific crop, such as somatic embryogenesis in plum, pollen germination in mei, and hormone regulation in peach. Notably, selective sweeps on chromosome 7 coincide with a chromosomal collinearity from the comparative genomics, impacting key fruit-softening genes such as PG, regulated by ERF and RMA1H1. Overall, this study provides insights into the genetic diversity, evolutionary history, and domestication of the APPM complex, offering valuable implications for genetic studies and breeding programs of Prunus crops.

Phylogenomics of Afrotherian mammals and improved resolution of extant Paenungulata

Molecular investigations have gathered a diverse set of mammals-predominantly African natives like elephants, hyraxes, and aardvarks-into a clade known as Afrotheria. Nevertheless, the precise phylogenetic relationships among these species remain contentious. Here, we sourced orthologous markers and ultraconserved elements to discern the interordinal connections among Afrotherian mammals. Our phylogenetic analyses bolster the common origin of Afroinsectiphilia and Paenungulata, and propose Afrosoricida as the closer relative to Macroscelidea rather than Tubulidentata, while also challenging the notion of Sirenia and Hyracoidea as sister taxa. The approximately unbiased test and the gene concordance factor uniformly recognized the alliance of Proboscidea with Hyracoidea as the dominant topology within Paenungulata. Investigation into sites with extremly high phylogenetic signal unveiled their potential to intensify conflicts in the Paenungulata topology. Subsequent exploration suggested that incomplete lineage sorting was predominantly responsible for the observed contentious relationships, whereas introgression exerted a subsidiary influence. The divergence times estimated in our study hint at the Cretaceous-Paleogene (K-Pg) extinction event as a catalyst for Afrotherian diversification. Overall, our findings deliver a tentative but insightful overview of Afrotheria phylogeny and divergence, elucidating these relationships through the lens of phylogenomics.

A chromosome-level genome of Chenghua pig provides new insights into the domestication and local adaptation of pigs

As a country with abundant genetic resources of pigs, the domestication history of pigs in China and the adaptive evolution of Chinese pig breeds at different latitudes have rarely been elucidated at the genome-wide level. To fill this gap, we first assembled a high-quality chromosome-level genome of the Chenghua pig and used it as a benchmark to analyse the genomes of 272 samples from three genera of three continents. The divergence of the three species belonging to three genera, Phacochoerus africanus, Potamochoerus porcus, and Sus scrofa, was assessed. The introgression of pig breeds redefined that the migration routes were basically from southern China to central and southwestern China, then spread to eastern China, arrived in northern China, and finally reached Europe. The domestication of pigs in China occurred ∼12,000 years ago, earlier than the available Chinese archaeological domestication evidence. In addition, FBN1 and NR6A1 were identified in our study as candidate genes related to extreme skin thickness differences in Eurasian pig breeds and adaptive evolution at different latitudes in Chinese pig breeds, respectively. Our study provides a new resource for the pig genomic pool and refines our understanding of pig genetic diversity, domestication, migration, and adaptive evolution at different latitudes.

Genomic analyses provide insights into sex differentiation of tetraploid strawberry (Fragaria moupinensis)

The strawberry genus, Fragaria, exhibits a wide range of sexual systems and natural ploidy variation. Nearly, all polyploid strawberry species exhibit separate sexes (dioecy). Research has identified the sex-determining sequences as roughly conserved but with repeatedly changed genomic locations across octoploid strawberries. However, it remains unclear whether tetraploid wild strawberries evolved dioecy independently or shared a common origin with octoploid strawberries. In this study, we investigated the sex determinants of F. moupinensis, a dioecious plant with heterogametic females (ZW). Utilizing a combination of haplotype-resolved genome sequencing of the female F. moupinensis, k-mer-based and coverage-based genome-wide association studies (GWAS), and transcriptomic analysis, we discovered a non-recombining, approximately 33.6 kb W-specific region on chromosome 2a. Within this region, only one candidate sex-determining gene (FmoAFT) was identified. Furthermore, an extensive resequencing of the entire Fragaria genus indicated that the W-specific region displays conservative female specificity across all tetraploid species. This observation suggests that dioecy evolved independently in tetraploid and octoploid strawberries. Moreover, employing virus-induced gene silencing (VIGS), we knocked down the expression of the FmoAFT homologue transcript in cultivated strawberries, revealing its potential role in promoting female functions during early carpel development. We also applied DNA affinity purification sequencing (DAP-seq) and yeast one-hybrid assays to identify potential direct targets of FmoAFT. These insights shed new light on the genetic basis and evolutionary history of sex determination in strawberries, thereby facilitating the formulation of strategies to manipulate sex determination in breeding programs.

Adaptive divergence, historical population dynamics, and simulation of suitable distributions for Picea Meyeri and P. Mongolica at the whole-genome level

The taxonomic classification of Picea meyeri and P. mongolica has long been controversial. To investigate the genetic relatedness, evolutionary history, and population history dynamics of these species, genotyping-by-sequencing (GBS) technology was utilized to acquire whole-genome single nucleotide polymorphism (SNP) markers, which were subsequently used to assess population structure, population dynamics, and adaptive differentiation. Phylogenetic and population structural analyses at the genomic level indicated that although the ancestor of P. mongolica was a hybrid of P. meyeri and P. koraiensis, P. mongolica is an independent Picea species. Additionally, P. mongolica is more closely related to P. meyeri than to P. koraiensis, which is consistent with its geographic distribution. There were up to eight instances of interspecific and intraspecific gene flow between P. meyeri and P. mongolica. The P. meyeri and P. mongolica effective population sizes generally decreased, and Maxent modeling revealed that from the Last Glacial Maximum (LGM) to the present, their habitat areas decreased initially and then increased. However, under future climate scenarios, the habitat areas of both species were projected to decrease, especially under high-emission scenarios, which would place P. mongolica at risk of extinction and in urgent need of protection. Local adaptation has promoted differentiation between P. meyeri and P. mongolica. Genotype‒environment association analysis revealed 96,543 SNPs associated with environmental factors, mainly related to plant adaptations to moisture and temperature. Selective sweeps revealed that the selected genes among P. meyeri, P. mongolica and P. koraiensis are primarily associated in vascular plants with flowering, fruit development, and stress resistance. This research enhances our understanding of Picea species classification and provides a basis for future genetic improvement and species conservation efforts.

Introgression Underlies Phylogenetic Uncertainty But Not Parallel Plumage Evolution in a Recent Songbird Radiation

Instances of parallel phenotypic evolution offer great opportunities to understand the evolutionary processes underlying phenotypic changes. However, confirming parallel phenotypic evolution and studying its causes requires a robust phylogenetic framework. One such example is the "black-and-white wagtails," a group of 5 species in the songbird genus Motacilla: 1 species, Motacilla alba, shows wide intra-specific plumage variation, while the 4r others form 2 pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the 2 species in each of these pairs were not recovered as sisters in previous phylogenetic inferences. Their relationships varied depending on the markers used, suggesting that gene tree heterogeneity might have hampered accurate phylogenetic inference. Here, we use whole genome resequencing data to explore the phylogenetic relationships within this group, with a special emphasis on characterizing the extent of gene tree heterogeneity and its underlying causes. We first used multispecies coalescent methods to generate a "complete evidence" phylogenetic hypothesis based on genome-wide variants, while accounting for incomplete lineage sorting (ILS) and introgression. We then investigated the variation in phylogenetic signal across the genome to quantify the extent of discordance across genomic regions and test its underlying causes. We found that wagtail genomes are mosaics of regions supporting variable genealogies, because of ILS and inter-specific introgression. The most common topology across the genome, supporting M. alba and M. aguimp as sister species, appears to be influenced by ancient introgression. Additionally, we inferred another ancient introgression event, between M. alba and M. grandis. By combining results from multiple analyses, we propose a phylogenetic network for the black-and-white wagtails that confirms that similar phenotypes evolved in non-sister lineages, supporting parallel plumage evolution. Furthermore, the inferred reticulations do not connect species with similar plumage coloration, suggesting that introgression does not underlie parallel plumage evolution in this group. Our results demonstrate the importance of investing genome-wide patterns of gene tree heterogeneity to help understand the mechanisms underlying phenotypic evolution. [Gene tree heterogeneity; incomplete lineage sorting; introgression; parallel evolution; phylogenomics; plumage evolution; wagtails.].

Recent selection and introgression facilitated high-altitude adaptation in cattle

During the past 3000 years, cattle on the Qinghai-Xizang Plateau have developed adaptive phenotypes under the selective pressure of hypoxia, ultraviolet (UV) radiation, and extreme cold. The genetic mechanism underlying this rapid adaptation is not yet well understood. Here, we present whole-genome resequencing data for 258 cattle from 32 cattle breeds/populations, including 89 Tibetan cattle representing eight populations distributed at altitudes ranging from 3400 m to 4300 m. Our genomic analysis revealed that Tibetan cattle exhibited a continuous phylogeographic cline from the East Asian taurine to the South Asian indicine ancestries. We found that recently selected genes in Tibetan cattle were related to body size (HMGA2 and NCAPG) and energy expenditure (DUOXA2). We identified signals of sympatric introgression from yak into Tibetan cattle at different altitudes, covering 0.64%-3.26% of their genomes, which included introgressed genes responsible for hypoxia response (EGLN1), cold adaptation (LRP11), DNA damage repair (LATS1), and UV radiation resistance (GNPAT). We observed that introgressed yak alleles were associated with noncoding variants, including those in present EGLN1. In Tibetan cattle, three yak introgressed SNPs in the EGLN1 promoter region reduced the expression of EGLN1, suggesting that these genomic variants enhance hypoxia tolerance. Taken together, our results indicated complex adaptation processes in Tibetan cattle, where recently selected genes and introgressed yak alleles jointly facilitated rapid adaptation to high-altitude environments.

Extreme elevational migration spurred cryptic speciation in giant hummingbirds

The ecoevolutionary drivers of species niche expansion or contraction are critical for biodiversity but challenging to infer. Niche expansion may be promoted by local adaptation or constrained by physiological performance trade-offs. For birds, evolutionary shifts in migratory behavior permit the broadening of the climatic niche by expansion into varied, seasonal environments. Broader niches can be short-lived if diversifying selection and geography promote speciation and niche subdivision across climatic gradients. To illuminate niche breadth dynamics, we can ask how "outlier" species defy constraints. Of the 363 hummingbird species, the giant hummingbird (Patagona gigas) has the broadest climatic niche by a large margin. To test the roles of migratory behavior, performance trade-offs, and genetic structure in maintaining its exceptional niche breadth, we studied its movements, respiratory traits, and population genomics. Satellite and light-level geolocator tracks revealed an >8,300-km loop migration over the Central Andean Plateau. This migration included a 3-wk, ~4,100-m ascent punctuated by upward bursts and pauses, resembling the acclimatization routines of human mountain climbers, and accompanied by surging blood-hemoglobin concentrations. Extreme migration was accompanied by deep genomic divergence from high-elevation resident populations, with decisive postzygotic barriers to gene flow. The two forms occur side-by-side but differ almost imperceptibly in size, plumage, and respiratory traits. The high-elevation resident taxon is the world's largest hummingbird, a previously undiscovered species that we describe and name here. The giant hummingbirds demonstrate evolutionary limits on niche breadth: when the ancestral niche expanded due to evolution (or loss) of an extreme migratory behavior, speciation followed.

The rise of baobab trees in Madagascar

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

Genomic signatures of strawberry domestication and diversification

Cultivated strawberry (Fragaria × ananassa) has a brief history of less than 300 yr, beginning with the hybridization of octoploids Fragaria chiloensis and Fragaria virginiana. Here we explored the genomic signatures of early domestication and subsequent diversification for different climates using whole-genome sequences of 289 wild, heirloom, and modern varieties from two major breeding programs in the United States. Four nonadmixed wild octoploid populations were identified, with recurrent introgression among the sympatric populations. The proportion of F. virginiana ancestry increased by 20% in modern varieties over initial hybrids, and the proportion of F. chiloensis subsp. pacifica rose from 0% to 3.4%. Effective population size rapidly declined during early breeding. Meanwhile, divergent selection for distinct environments reshaped wild allelic origins in 21 out of 28 chromosomes. Overlapping divergent selective sweeps in natural and domesticated populations revealed 16 convergent genomic signatures that may be important for climatic adaptation. Despite 20 breeding cycles since initial hybridization, more than half of loci underlying yield and fruit size are still not under artificial selection. These insights add clarity to the domestication and breeding history of what is now the most widely cultivated fruit in the world.

The genetics of niche-specific behavioral tendencies in an adaptive radiation of cichlid fishes

Behavior is critical for animal survival and reproduction, and possibly for diversification and evolutionary radiation. However, the genetics behind adaptive variation in behavior are poorly understood. In this work, we examined a fundamental and widespread behavioral trait, exploratory behavior, in one of the largest adaptive radiations on Earth, the cichlid fishes of Lake Tanganyika. By integrating quantitative behavioral data from 57 cichlid species (702 wild-caught individuals) with high-resolution ecomorphological and genomic information, we show that exploratory behavior is linked to macrohabitat niche adaptations in Tanganyikan cichlids. Furthermore, we uncovered a correlation between the genotypes at a single-nucleotide polymorphism upstream of the AMPA glutamate-receptor regulatory gene cacng5b and variation in exploratory tendency. We validated this association using behavioral predictions with a neural network approach and CRISPR-Cas9 genome editing.

Balancing selection and the functional effects of shared polymorphism in cryptic Daphnia species

The patterns of genetic variation within and between related taxa represent the genetic history of a species. Shared polymorphisms, loci with identical alleles across species, are of unique interest as they may represent cases of ancient selection maintaining functional variation post-speciation. In this study, we investigate the abundance of shared polymorphism in the Daphnia pulex species complex. We test whether shared mutations are consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting, or convergent evolution. We analyzed over 2,000 genomes from North American and European D. pulex and several outgroup species to examine the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that while North American and European D. pulex diverged over ten million years ago, they retained tens of thousands of shared alleles. We found that the number of shared polymorphisms between North American and European D. pulex cannot be explained by hybridization or incomplete lineage sorting alone. Instead, we show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting young and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild. Ultimately, our findings provide insights into the genetic basis of adaptation and the maintenance of genetic diversity between species.

From sub-Saharan Africa to China: Evolutionary history and adaptation of Drosophila melanogaster revealed by population genomics

Drosophila melanogaster is a widely used model organism for studying environmental adaptation. However, the genetic diversity of populations in Asia is poorly understood, leaving a notable gap in our knowledge of the global evolution and adaptation of this species. We sequenced genomes of 292 D. melanogaster strains from various ecological settings in China and analyzed them along with previously published genome sequences. We have identified six global genetic ancestry groups, despite the presence of widespread genetic admixture. The strains from China represent a unique ancestry group, although detectable differentiation exists among populations within China. We deciphered the global migration and demography of D. melanogaster, and identified widespread signals of adaptation, including genetic changes in response to insecticides. We validated the effects of insecticide resistance variants using population cage trials and deep sequencing. This work highlights the importance of population genomics in understanding the genetic underpinnings of adaptation, an effort that is particularly relevant given the deterioration of ecosystems.

Introgression and disruption of migration routes have shaped the genetic integrity of wildebeest populations

The blue wildebeest (Connochaetes taurinus) is a keystone species in savanna ecosystems from southern to eastern Africa, and is well known for its spectacular migrations and locally extreme abundance. In contrast, the black wildebeest (C. gnou) is endemic to southern Africa, barely escaped extinction in the 1900s and is feared to be in danger of genetic swamping from the blue wildebeest. Despite the ecological importance of the wildebeest, there is a lack of understanding of how its unique migratory ecology has affected its gene flow, genetic structure and phylogeography. Here, we analyze whole genomes from 121 blue and 22 black wildebeest across the genus' range. We find discrete genetic structure consistent with the morphologically defined subspecies. Unexpectedly, our analyses reveal no signs of recent interspecific admixture, but rather a late Pleistocene introgression of black wildebeest into the southern blue wildebeest populations. Finally, we find that migratory blue wildebeest populations exhibit a combination of long-range panmixia, higher genetic diversity and lower inbreeding levels compared to neighboring populations whose migration has recently been disrupted. These findings provide crucial insights into the evolutionary history of the wildebeest, and tangible genetic evidence for the negative effects of anthropogenic activities on highly migratory ungulates.

Population genomics of the Isoetes appalachiana (Isoetaceae) complex supports a 'diploids-first' approach to conservation

BACKGROUND AND AIMS

Allopolyploidy is an important driver of diversification and a key contributor to genetic novelty across the tree of life. However, many studies have questioned the importance of extant polyploid lineages, suggesting that the vast majority may constitute evolutionary 'dead ends'. This has important implications for conservation efforts where polyploids and diploid progenitors often compete for wildlife management resources. Isoetes appalachiana is an allotetraploid that is broadly distributed throughout the eastern USA alongside its diploid progenitors, I. valida and I. engelmannii. As such, this species complex provides an excellent opportunity to investigate the processes that underpin the formation and survival of allopolyploid lineages.

METHODS

Here we utilized RADseq and whole-chloroplast sequencing to unravel the demographic and evolutionary history of hybridization in this widespread species complex. We developed a modified protocol for phasing RADseq loci from an allopolyploid in order to examine each progenitor's genetic contribution independently in a phylogenetic context. Additionally, we conducted population-level analyses to examine genetic diversity and evidence of gene flow within species.

KEY RESULTS

Isoetes appalachiana is the product of multiple phylogenetic origins, suggesting that formation and establishment of allopolyploids are common in this group. Hybridization appears to be unidirectional, with I. engelmannii consistently being the maternal progenitor. Additionally, we find that polyploid lineages are genetically isolated, rarely if ever experiencing gene flow between geographically distinct populations.

CONCLUSIONS

Allopolyploid lineages of I. appalachiana appear to form frequently and experience a high degree of genetic isolation following formation. Thus, our results appear to corroborate the hypothesis that the vast majority of recently formed polyploids may represent evolutionary dead ends. However, this does not necessarily lessen the evolutionary importance or ecological impact of polyploidy per se. Accordingly, we propose a conservation strategy that prioritizes diploid taxa, thus preserving downstream processes that recurrently generate allopolyploid diversity.

Giraffe lineages are shaped by major ancient admixture events

Strong genetic structure has prompted discussion regarding giraffe taxonomy,1,2,3 including a suggestion to split the giraffe into four species: Northern (Giraffa c. camelopardalis), Reticulated (G. c. reticulata), Masai (G. c. tippelskirchi), and Southern giraffes (G. c. giraffa).4,5,6 However, their evolutionary history is not yet fully resolved, as previous studies used a simple bifurcating model and did not explore the presence or extent of gene flow between lineages. We therefore inferred a model that incorporates various evolutionary processes to assess the drivers of contemporary giraffe diversity. We analyzed whole-genome sequencing data from 90 wild giraffes from 29 localities across their current distribution. The most basal divergence was dated to 280 kya. Genetic differentiation, FST, among major lineages ranged between 0.28 and 0.62, and we found significant levels of ancient gene flow between them. In particular, several analyses suggested that the Reticulated lineage evolved through admixture, with almost equal contribution from the Northern lineage and an ancestral lineage related to Masai and Southern giraffes. These new results highlight a scenario of strong differentiation despite gene flow, providing further context for the interpretation of giraffe diversity and the process of speciation in general. They also illustrate that conservation measures need to target various lineages and sublineages and that separate management strategies are needed to conserve giraffe diversity effectively. Given local extinctions and recent dramatic declines in many giraffe populations, this improved understanding of giraffe evolutionary history is relevant for conservation interventions, including reintroductions and reinforcements of existing populations.

Adaptive introgression reveals the genetic basis of a sexually selected syndrome in wall lizards

The joint expression of particular colors, morphologies, and behaviors is a common feature of adaptation, but the genetic basis for such "phenotypic syndromes" remains poorly understood. Here, we identified a complex genetic architecture associated with a sexually selected syndrome in common wall lizards, by capitalizing on the adaptive introgression of coloration and morphology into a distantly related lineage. Consistent with the hypothesis that the evolution of phenotypic syndromes in vertebrates is facilitated by developmental linkage through neural crest cells, most of the genes associated with the syndrome are involved in neural crest cell regulation. A major locus was a ~400-kb region, characterized by standing structural genetic variation and previously implied in the evolutionary innovation of coloration and beak size in birds. We conclude that features of the developmental and genetic architecture contribute to maintaining trait integration, facilitating the extensive and rapid introgressive spread of suites of sexually selected characters.

Genetic introgression from commercial European pigs to the indigenous Chinese Lijiang breed and associated changes in phenotypes

BACKGROUND

Gene flow is crucial for enhancing economic traits of livestock. In China, breeders have used hybridization strategies for decades to improve livestock performance. Here, we performed whole-genome sequencing of a native Chinese Lijiang pig (LJP) breed. By integrating previously published data, we explored the genetic structure and introgression of genetic components from commercial European pigs (EP) into the LJP, and examined the impact of this introgression on phenotypic traits.

RESULTS

Our analysis revealed significant introgression of EP breeds into the LJP and other domestic pig breeds in China. Using a haplotype-based approach, we quantified introgression levels and compared EP to LJP and other Chinese domestic pigs. The results show that EP introgression is widely prevalent in Chinese domestic pigs, although there are significant differences between breeds. We propose that LJP could potentially act as a mediator for the transmission of EP haplotypes. We also examined the correlation between EP introgression and the number of thoracic vertebrae in LJP and identified VRTN and STUM as candidate genes for this trait.

CONCLUSIONS

Our study provides evidence of introgressed European haplotypes in the LJP breed and describes the potential role of EP introgression on phenotypic changes of this indigenous breed.

Species delimitation in Amana (Liliaceae): transcriptomes battle with evolutionary complexity

Species delimitation has long been a subject of controversy, and there are many alternative concepts and approaches used to define species in plants. The genus Amana (Liliaceae), known as "East Asian tulips" has a number of cryptic species and a huge genome size (1C = 21.48-57.35 pg). It also is intriguing how such a spring ephemeral genus thrives in subtropical areas. However, phylogenetic relationships and species delimitation within Amana are challenging. Here we included all species and 84 populations of Amana, which are collected throughout its distribution range. A variety of methods were used to clarify its species relationships based on a combination of morphological, ecological, genetic, evolutionary and phylogenetic species concepts. This evidence supports the recognition of at least 12 species in Amana. Moreover, we explored the complex evolutionary history within the genus and detected several historical hybridization and introgression events based on phylogenetic trees (transcriptomic and plastid), phylonetworks, admixture and ABBA-BABA analyses. Morphological traits have undergone parallel evolution in the genus. This spring ephemeral genus might have originated from a temperate region, yet finally thrives in subtropical areas, and three hypotheses about its adaptive evolution are proposed for future testing. In addition, we propose a new species, Amana polymorpha, from eastern Zhejiang Province, China. This research also demonstrates that molecular evidence at the genome level (such as transcriptomes) has greatly improved the accuracy and reasonability of species delimitation and taxon classification.

Genomic evidence for evolutionary history and local adaptation of two endemic apricots: Prunus hongpingensis and P. zhengheensis

Apricot, belonging to the Armeniaca section of Rosaceae, is one of the economically important crop fruits that has been extensively cultivated. The natural wild apricots offer valuable genetic resources for crop improvement. However, some of them are endemic, with small populations, and are even at risk of extinction. In this study we unveil chromosome-level genome assemblies for two southern China endemic apricots, Prunus hongpingensis (PHP) and P. zhengheensis (PZH). We also characterize their evolutionary history and the genomic basis of their local adaptation using whole-genome resequencing data. Our findings reveal that PHP and PZH are closely related to Prunus armeniaca and form a distinct lineage. Both species experienced a decline in effective population size following the Last Glacial Maximum (LGM), which likely contributed to their current small population sizes. Despite the observed decrease in genetic diversity and heterozygosity, we do not observe an increased accumulation of deleterious mutations in these two endemic apricots. This is likely due to the combined effects of a low inbreeding coefficient and strong purifying selection. Furthermore, we identify a set of genes that have undergone positive selection and are associated with local environmental adaptation in PHP and PZH, respectively. These candidate genes can serve as valuable genetic resources for targeted breeding and improvement of cultivated apricots. Overall, our study not only enriches our comprehension of the evolutionary history of apricot species but also offers crucial insights for the conservation and future breeding of other endemic species amidst rapid climate changes.

The role of introgressive hybridization in shaping the geographically isolated gene pools of wax palm populations (genus Ceroxylon)

The speciation continuum is the process by which genetic groups diverge until they reach reproductive isolation. It has become common in the literature to show that this process is gradual and flickering, with possibly many instances of secondary contact and introgression after divergence has started. The level of divergence might vary among genomic regions due to, among others, the different forces and roles of selection played by the shared regions. Through hybrid capture, we sequenced ca. 4,000 nuclear regions in populations of six species of wax palms, five of which form a monophyletic group (genus Ceroxylon, Arecaceae: Ceroxyloideae). We show that in this group, the different populations show varying degrees of introgressive hybridization, and two of them are backcrosses of the other three 'pure' species. This is particularly interesting because these three species are dioecious, have a shared main pollinator, and have slightly overlapping reproductive seasons but highly divergent morphologies. Our work supports shows wax palms diverge under positive and background selection in allopatry, and hybridize due to secondary contact and inefficient reproductive barriers, which sustain genetic diversity. Introgressed regions are generally not under positive selection. Peripheral populations are backcrosses of other species; thus, introgressive hybridization is likely modulated by demographic effects rather than selective pressures. In general, these species might function as an 'evolutionary syngameon' where expanding, peripheral, small, and isolated populations maintain diversity by crossing with available individuals of other wax palms. In the Andean context, species can benefit from gained variation from a second taxon or the enhancement of population sizes by recreating a common genetic pool.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for over a century. However, laboratory mice capture only a subset of the genetic variation found in wild mouse populations, ultimately limiting the potential of classical inbred strains to uncover phenotype-associated variants and pathways. Wild mouse populations are reservoirs of genetic diversity that could facilitate the discovery of new functional and disease-associated alleles, but the scarcity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently developed, sequenced, and phenotyped a set of 11 inbred strains derived from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from one of five environmentally distinct locations across North and South America. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the GRCm39 mouse reference, with 42.5% of variants in the Nachman strain genomes absent from current classical inbred mouse strain panels. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels. These novel wild-derived inbred mouse strain resources are set to empower new discoveries in both basic and preclinical research.

Significant genomic introgression from grey junglefowl (Gallus sonneratii) to domestic chickens (Gallus gallus domesticus)

BACKGROUND

Chicken is one of the most numerous and widely distributed species around the world, and many studies support the multiple ancestral origins of domestic chickens. The research regarding the yellow skin phenotype in domestic chickens (regulated by BCO2) likely originating from the grey junglefowl serves as crucial evidence for demonstrating the multiple origins of chickens. However, beyond the BCO2 gene region, much remains unknown about the introgression from the grey junglefowl into domestic chickens. Therefore, in this study, based on whole-genome data of 149 samples including 4 species of wild junglefowls and 13 local domestic chicken breeds, we explored the introgression events from the grey junglefowl to domestic chickens.

RESULTS

We successfully detected introgression regions besides BCO2, including two associated with growth trait (IGFBP2 and TKT), one associated with angiogenesis (TIMP3) and two members of the heat shock protein family (HSPB2 and CRYAB). Our findings suggest that the introgression from the grey junglefowl may impact the growth performance of chickens. Furthermore, we revealed introgression events from grey junglefowl at the BCO2 region in multiple domestic chicken breeds, indicating a phenomenon where the yellow skin phenotype likely underwent strong selection and was retained. Additionally, our haplotype analysis shed light on BCO2 introgression event from different sources of grey junglefowl into domestic chickens, possibly suggesting multiple genetic flows between the grey junglefowl and domestic chickens.

CONCLUSIONS

In summary, our findings provide evidences of the grey junglefowl contributing to the genetic diversity of domestic chickens, laying the foundation for a deeper understanding of the genetic composition within domestic chickens, and offering new perspectives on the impact of introgression on domestic chickens.

High genetic load without purging in caribou, a diverse species at risk

High intra-specific genetic diversity is associated with adaptive potential, which is key for resilience to global change. However, high variation may also support deleterious alleles through genetic load, thereby increasing the risk of inbreeding depression if population sizes decrease. Purging of deleterious variation has been demonstrated in some threatened species. However, less is known about the costs of declines and inbreeding in species with large population sizes and high genetic diversity even though this encompasses many species globally that are expected to undergo population declines. Caribou is a species of ecological and cultural significance in North America with a wide distribution supporting extensive phenotypic variation but with some populations undergoing significant declines resulting in their at-risk status in Canada. We assessed intra-specific genetic variation, adaptive divergence, inbreeding, and genetic load across populations with different demographic histories using an annotated chromosome-scale reference genome and 66 whole-genome sequences. We found high genetic diversity and nine phylogenomic lineages across the continent with adaptive diversification of genes, but also high genetic load among lineages. We found highly divergent levels of inbreeding across individuals, including the loss of alleles by drift but not increased purging in inbred individuals, which had more homozygous deleterious alleles. We also found comparable frequencies of homozygous deleterious alleles between lineages regardless of nucleotide diversity. Thus, further inbreeding may need to be mitigated through conservation efforts. Our results highlight the "double-edged sword" of genetic diversity that may be representative of other species atrisk affected by anthropogenic activities.

Recent genetic, phenetic and ecological divergence across the Mesoamerican highlands: a study case with Diglossa baritula (Aves: Thraupidae)

The topographical, geological, climatic and biodiversity complexity of Mesoamerica has made it a primary research focus. The Mesoamerican highlands is a region with particularly high species richness and within-species variation. The Cinnamon-bellied Flowerpiercer, Diglossa baritula (Wagler, 1832), is a species endemic to the Mesoamerican highlands, with three allopatric subspecies currently recognized. To characterize divergence within this species, we integrated genomics, morphology, coloration and ecological niche modeling approaches, obtained from sampling individuals across the entire geographic distribution of the species. Our results revealed a clear genomic divergence between the populations to the east versus the west of the Isthmus of Tehuantepec. In contrast to the genomic results, morphology and coloration analyses showed intermediate levels of differentiation, indicating that population groups within D. baritula have probably been under similar selective pressures. Our morphology results indicated that the only sexually dimorphic morphological variable is the wing chord, with males having a longer wing chord than females. Finally, ecological data indicated that there are differences in ecological niche within D. baritula. Our data suggest that D. baritula could contain two or more incipient species at the intermediate phase of the speciation continuum. These results highlight the importance of the geographical barrier of the Isthmus of Tehuantepec and Pleistocene climatic events in driving isolation and population divergence in D. baritula. The present investigation illustrates the speciation potential of the D. baritula complex and the capacity of Mesoamerican highlands to create cryptic biodiversity and endemism.

Evolutionary radiation of the Eurasian Pinus species under pervasive gene flow

Evolutionary radiation, a pivotal aspect of macroevolution, offers valuable insights into evolutionary processes. The genus Pinus is the largest genus in conifers withc . $$ c. $$ 90% of the extant species emerged in the Miocene, which signifies a case of rapid diversification. Despite this remarkable history, our understanding of the mechanisms driving radiation within this expansive genus has remained limited. Using exome capture sequencing and a fossil-calibrated phylogeny, we investigated the divergence history, niche diversification, and introgression among 13 closely related Eurasian species spanning climate zones from the tropics to the boreal Arctic. We detected complex introgression among lineages in subsection Pinus at all stages of the phylogeny. Despite this widespread gene exchange, each species maintained its genetic identity and showed clear niche differentiation. Demographic analysis unveiled distinct population histories among these species, which further influenced the nucleotide diversity and efficacy of purifying and positive selection in each species. Our findings suggest that radiation in the Eurasian pines was likely fueled by interspecific recombination and further reinforced by their adaptation to distinct environments. Our study highlights the constraints and opportunities for evolutionary change, and the expectations of future adaptation in response to environmental changes in different lineages.

Invasive fall armyworms are corn strain

The fall armyworm (Spodoptera frugiperda) is one of the major pest insects in diverse crop plants, including maize, rice, and cotton. While the fall armyworm is native to North and South America, its invasion was first reported in West Africa in 2016. Since then, this species has rapidly spread across Sub-Saharan Africa, Asia, and Oceania, as well as Egypt and Cyprus. The fall armyworm is composed of two sympatric strains, the corn and rice strains, designated to their preferred host plants, in native areas. It remains surprisingly unclear whether invasive fall armyworms belong to the corn strain, rice strain, or hybrids of the two, despite a large number of population genetics studies. In this study, we performed population genomics analyses using globally collected 116 samples to identify the strains of invasive fall armyworms. We observed that invasive fall armyworms are genomically most similar to the corn strain. The reconstructed phylogenetic tree supports the hypothesis that invasive fall armyworms originated from the corn strain. All genomic loci of invasive populations exhibit higher genetic similarity to the corn strains compared to the rice strains. Furthermore, we found no evidence of gene flow from rice strains to invasive populations at any genomic locus. These results demonstrate that invasive fall armyworms belong to the corn strain. These results suggest that invasive fall armyworms likely have very limited potential to infest rice. Therefore, the management plan should primarily focus on crops preferred by the corn strain.

Origin and diversity of Capsella bursa-pastoris from the genomic point of view

BACKGROUND

Capsella bursa-pastoris, a cosmopolitan weed of hybrid origin, is an emerging model object for the study of early consequences of polyploidy, being a fast growing annual and a close relative of Arabidopsis thaliana. The development of this model is hampered by the absence of a reference genome sequence.

RESULTS

We present here a subgenome-resolved chromosome-scale assembly and a genetic map of the genome of Capsella bursa-pastoris. It shows that the subgenomes are mostly colinear, with no massive deletions, insertions, or rearrangements in any of them. A subgenome-aware annotation reveals the lack of genome dominance-both subgenomes carry similar number of genes. While most chromosomes can be unambiguously recognized as derived from either paternal or maternal parent, we also found homeologous exchange between two chromosomes. It led to an emergence of two hybrid chromosomes; this event is shared between distant populations of C. bursa-pastoris. The whole-genome analysis of 119 samples belonging to C. bursa-pastoris and its parental species C. grandiflora/rubella and C. orientalis reveals introgression from C. orientalis but not from C. grandiflora/rubella.

CONCLUSIONS

C. bursa-pastoris does not show genome dominance. In the earliest stages of evolution of this species, a homeologous exchange occurred; its presence in all present-day populations of C. bursa-pastoris indicates on a single origin of this species. The evidence coming from whole-genome analysis challenges the current view that C. grandiflora/rubella was a direct progenitor of C. bursa-pastoris; we hypothesize that it was an extinct (or undiscovered) species sister to C. grandiflora/rubella.

Cryptic diversity and rampant hybridization in annual gentians on the Qinghai-Tibet Plateau revealed by population genomic analysis

Understanding the evolutionary and ecological processes involved in population differentiation and speciation provides critical insights into biodiversity formation. In this study, we employed 29,865 single nucleotide polymorphisms (SNPs) and complete plastomes to examine genomic divergence and hybridization in Gentiana aristata, which is endemic to the Qinghai-Tibet Plateau (QTP) region. Genetic clustering revealed that G. aristata is characterized by geographic genetic structures with five clusters (West, East, Central, South and North). The West cluster has a specific morphological character (i.e., blue corolla) and higher values of FST compared to the remaining clusters, likely the result of the geological barrier formed by the Yangtze River. The West cluster diverged from the other clusters in the Early Pliocene; these remaining clusters diverged from one another in the Early Quaternary. Phylogenetic reconstructions based on SNPs and plastid data revealed substantial cyto-nuclear conflicts. Genetic clustering and D-statistics demonstrated rampant hybridization between the Central and North clusters, along the Bayankala Mountains, which form the geological barrier between the Central and North clusters. Species distribution modeling demonstrated the range of G. aristata expanded since the Last Interglacial period. Our findings provide genetic and morphological evidence of cryptic diversity in G. aristata, and identified rampant hybridization between genetic clusters along a geological barrier. These findings suggest that geological barriers and climatic fluctuations have an important role in triggering diversification as well as hybridization, indicating that cryptic diversity and hybridization are essential factors in biodiversity formation within the QTP region.

Gene flow and an anomaly zone complicate phylogenomic inference in a rapidly radiated avian family (Prunellidae)

BACKGROUND

Resolving the phylogeny of rapidly radiating lineages presents a challenge when building the Tree of Life. An Old World avian family Prunellidae (Accentors) comprises twelve species that rapidly diversified at the Pliocene-Pleistocene boundary.

RESULTS

Here we investigate the phylogenetic relationships of all species of Prunellidae using a chromosome-level de novo assembly of Prunella strophiata and 36 high-coverage resequenced genomes. We use homologous alignments of thousands of exonic and intronic loci to build the coalescent and concatenated phylogenies and recover four different species trees. Topology tests show a large degree of gene tree-species tree discordance but only 40-54% of intronic gene trees and 36-75% of exonic genic trees can be explained by incomplete lineage sorting and gene tree estimation errors. Estimated branch lengths for three successive internal branches in the inferred species trees suggest the existence of an empirical anomaly zone. The most common topology recovered for species in this anomaly zone was not similar to any coalescent or concatenated inference phylogenies, suggesting presence of anomalous gene trees. However, this interpretation is complicated by the presence of gene flow because extensive introgression was detected among these species. When exploring tree topology distributions, introgression, and regional variation in recombination rate, we find that many autosomal regions contain signatures of introgression and thus may mislead phylogenetic inference. Conversely, the phylogenetic signal is concentrated to regions with low-recombination rate, such as the Z chromosome, which are also more resistant to interspecific introgression.

CONCLUSIONS

Collectively, our results suggest that phylogenomic inference should consider the underlying genomic architecture to maximize the consistency of phylogenomic signal.

Genomic evidence for climate-linked diversity loss and increased vulnerability of wild barley spanning 28 years of climate warming

The information on how plant populations respond genetically to climate warming is scarce. Here, landscape genomic and machine learning approaches were integrated to assess genetic response of 10 wild barley (Hordeum vulgare ssp. spontaneum; WB) populations in the past and future, using whole genomic sequencing (WGS) data. The WB populations were sampled in 1980 and again in 2008. Phylogeny of accessions was roughly in conformity with sampling sites, which accompanied by admixture/introgressions. The 28-y climate warming resulted in decreased genetic diversity, increased selection pressure, and an increase in deleterious single nucleotide polymorphism (dSNP) numbers, heterozygous deleterious and total deleterious burdens for WB. Genome-environment associations identified some candidate genes belonging to peroxidase family (HORVU2Hr1G057450, HORVU4Hr1G052060 and HORVU4Hr1G057210) and heat shock protein 70 family (HORVU2Hr1G112630). The gene HORVU2Hr1G120170 identified by selective sweep analysis was under strong selection during the climate warming of the 28-y, and its derived haplotypes were fixed by WB when faced with the 28-y increasingly severe environment. Temperature variables were found to be more important than precipitation variables in influencing genomic variation, with an eco-physiological index gdd5 (growing degree-days at the baseline threshold temperature of 5 °C) being the most important determinant. Gradient forest modelling revealed higher predicted genomic vulnerability in Sede Boqer under future climate scenarios at 2041-2070 and 2071-2100. Additionally, estimates of effective population size (Ne) tracing back to 250 years indicated a forward decline in all populations over time. Our assessment about past genetic response and future vulnerability of WB under climate warming is crucial for informing conservation efforts for wild cereals and rational use strategies.

The phylogeographic journey of a plant species from lowland to highlands during the Pleistocene

Phylogeographic history refers to how species evolve and diversify in response to historical, ecological, and demographic factors. The climate fluctuation during the Pleistocene period marked a crucial time in shaping many species' distribution and genetic structure, particularly those from southern South American grasslands. This work investigated the phylogeographic history of a highland grassland, Petunia altiplana T. Ando & Hashim. (Solanaceae), its diversity, and geographic distribution using a population genomic approach based on RAD-seq data. Our results indicated that, during the Pleistocene, when the grasslands expanded to highlands, the lowland populations of P. altiplana reached the higher open fields, enlarging their geographic distribution. We found that the P. altiplana genetic diversity followed the geographic division into eastern (E) and western (WE) population groups, with a subtle division in the E group regarding the Pelotas River headwater. The results also showed that isolation by distance was the main divergence pattern, with elevation playing a pivotal role in shaping WE and E groups. Our findings indicated that lowland-adapted populations quickly colonized highlands during the late Pleistocene.