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.

Population genetics assessment of two pocilloporid coral species from the northern red sea: Implications for urbanized reef sustainability.

Understanding the genetic makeup of key coral species is vital for effective coral reef management, as heightened genetic diversity directly influences long-term survival and resilience against environmental changes. This study focused on two widespread Indo-Pacific branching corals, Pocillopora damicornis (referred as Pocillopora cf. damicornis (as identified only morphologically) and Seriatopora hystrix, by genotyping 222 and 195 colonies, respectively, from 10 sites in the northern Gulf of Eilat, Red Sea, using six and five microsatellite markers, respectively. Both species exhibited low observed heterozygosity (0.47 for P. cf. damicornis, 0.32 for S. hystrix) and similar expected heterozygosity (0.576 for P. cf. damicornis, 0.578 for S. hystrix). Pocillopora cf. damicornis showed minimal deviations from Hardy-Weinberg equilibrium (HWE) and low but positive F values, indicating high gene flow, while S. hystrix exhibited higher diversion from HWE and positive F values, suggesting isolation by distance and possible non-random mating or genetic drift. As the Gulf of Eilat undergoes rapid urbanization, this study highlights the anthropogenic impacts on the population genetics of key ecosystem engineering species and emphasizes the importance of managing genetics of Marine Protected Areas while implementing active coral reef restoration. The differences in reproductive traits between the two species (S. hystrix being a brooder, while P. cf. damicornis a broadcast spawner), underscore the need for sustainable population genetics management of the coral reefs for the future and resilience of the coral reef ecosystem of the northern Red Sea region.

A theory of measuring natural selection and genetic monitoring.

Two methods have been compared for determining the value of natural selection in the natural populations. The first method, based on the FST-statistics, employs the dependence of genetic diversity of a species on the value of gene flow between subpopulations of the species, derived from the assumption that all the mutations are close to selective neutrality, and subpopulations effect each other equally. Susceptibility to selection is estimated by the degree of deviation from this relationship between genetic diversity and gene flow in certain species. The second method is based on the probability theory and involves comparison between stabilities of the forms, competing in the population, which is computed using the data about fluctuations in their occurrence in several generations. As applied to the problems of genetic monitoring of rare and valuable species, the first method can be employed for express-assessment of susceptibility of a species to rapid intraspecific changes. The second method is suitable for a long-term and in-depth genetic monitoring of the species subjected to extremely intense natural selection of a disruptive or stabilizing form, which were revealed using the first method. There is a lack of long-term observations of intraspecific genetic variation of rare and protected species. The need for funds that finance long-term genetic research is substantiated.

Genetic diversity and population structure of Physella acuta (Gastropoda: Physidae) in Thailand using mitochondrial gene markers: COI and 16S rDNA.

Physella acuta is a freshwater snail native to North America. Understanding the phylogeography and genetic structure of P. acuta will help elucidate its evolution. In this study, we used mitochondrial (COI and 16S rDNA) and nuclear (ITS1) markers to identify the species and examine its genetic diversity, population structure, and demographic history of P. acuta in Thailand. Phylogenetic and network analyses of P. acuta in Thailand pertained to clade A, which exhibits a global distribution. Analysis of the genetic structure of the population revealed that the majority of pairwise comparisons showed no genetic dissimilarity. An isolation-by-distance test indicates no significant correlation between genetic and geographical distances among P. acuta populations, suggesting that gene flow is not restricted by distance. Demographic history and haplotype network analyses suggest a population expansion of P. acuta, as evidenced by the star-like structure detected in the median-joining network. Based on these results, we concluded that P. acuta in Thailand showed gene flow and recent population expansion. Our findings provide fundamental insights into the genetic variation of P. acuta in Thailand.

Uncovering ghost introgression through genomic analysis of a distinct eastern Asian hickory species.

Ghost introgression, or the transfer of genetic material from extinct or unsampled lineages to sampled species, has attracted much attention. However, conclusive evidence for ghost introgression, especially in plant species, remains scarce. Here, we newly assembled chromosome-level genomes for both Carya sinensis and Carya cathayensis, and additionally re-sequenced the whole genomes of 43 C. sinensis individuals as well as 11 individuals representing 11 diploid hickory species. These genomic datasets were used to investigate the reticulation and bifurcation patterns within the genus Carya (Juglandaceae), with a particular focus on the beaked hickory C. sinensis. By combining the D-statistic and BPP methods, we obtained compelling evidence that supports the occurrence of ghost introgression in C. sinensis from an extinct ancestral hickory lineage. This conclusion was reinforced through the phylogenetic network analysis and a genome scan method VolcanoFinder, the latter of which can detect signatures of adaptive introgression from unknown donors. Our results not only dispel certain misconceptions about the phylogenetic history of C. sinensis but also further refine our understanding of Carya's biogeography via divergence estimates. Moreover, the successful integration of the D-statistic and BPP methods demonstrates their efficacy in facilitating a more precise identification of introgression types.

Cryptic diversity within two widespread diadromous freshwater fishes (Teleostei: Galaxiidae).

Identification of taxonomically cryptic species is essential for the effective conservation of biodiversity. Freshwater-limited organisms tend to be genetically isolated by drainage boundaries, and thus may be expected to show substantial cryptic phylogenetic and taxonomic diversity. By comparison, populations of diadromous taxa, that migrate between freshwater and marine environments, are expected to show less genetic differentiation. Here we test for cryptic diversity in Australasian populations (both diadromous and non-diadromous) of two widespread Southern Hemisphere fish species, Galaxias brevipinnis and Galaxias maculatus. Both mtDNA and nuclear markers reveal putative cryptic species within these taxa. The substantial diversity detected within G. brevipinnis may be explained by its strong climbing ability which allows it to form isolated inland populations. In island populations, G. brevipinnis similarly show deeper genetic divergence than those of G. maculatus, which may be explained by the greater abundance of G. maculatus larvae in the sea allowing more ongoing dispersal. Our study highlights that even widespread, 'high-dispersal' species can harbour substantial cryptic diversity and therefore warrant increased taxonomic and conservation attention.

Machine learning provides insights for spatially explicit pest management strategies by integrating information on population connectivity and habitat use in a key agricultural pest.

BACKGROUND: Insect pests have garnered increasing interest because of anthropogenic global change, and their sustainable management requires knowledge of population habitat use and spread patterns. To enhance this knowledge for the prevalent tea pest Empoasca onukii, we utilized a random forest algorithm and a bivariate map to develop and integrate models of its habitat suitability and genetic connectivity across China. RESULTS: Our modeling revealed heterogeneous spatial patterns in suitability and connectivity despite the common key environmental predictor of isothermality. Analyses indicated that tea cultivation in areas surrounding the Tibetan Plateau and the southern tip of China may be at low risk of population outbreaks because of their predicted low suitability and connectivity. However, regions along the middle and lower reaches of the Yangtze River should consider the high abundance and high recolonization potential of E. onukii, and thus the importance of control measures. Our results also emphasized the need to prevent dispersal from outside regions in the areas north of the Yangtze River and highlighted the effectiveness of internal management efforts in southwestern China and along the southeastern coast. Further projections under future conditions suggested the potential for increased abundance and spread in regions north of the Yangtze River and the southern tip of China, and indicated the importance of long-term monitoring efforts in these areas. CONCLUSION: These findings highlighted the significance of combining information on habitat use and spread patterns for spatially explicit pest management planning. In addition, the approaches we used have potential applications in the management of other pest systems and the conservation of endangered biological resources. © 2024 Society of Chemical Industry.

The role of neutral and adaptive evolutionary processes on patterns of genetic diversity across small cave-dwelling populations of Icelandic Arctic charr (Salvelinus alpinus).

Understanding the adaptability of small populations in the face of environmental change is a central problem in evolutionary biology. Solving this problem is challenging because neutral evolutionary processes that operate on historical and contemporary timescales can override the effects of selection in small populations. We assessed the effects of isolation by colonization (IBC), isolation by dispersal limitation (IBDL) as reflected by a pattern of isolation by distance (IBD), and isolation by adaptation (IBA) and the roles of genetic drift and gene flow on patterns of genetic differentiation among 19 cave-dwelling populations of Icelandic Arctic charr (Salvelinus alpinus). We detected evidence of IBC based on the genetic affinity of nearby cave populations and the genetic relationships between the cave populations and the presumed ancestral population in the lake. A pattern of IBD was evident regardless of whether high-level genetic structuring (IBC) was taken into account. Genetic signatures of bottlenecks and lower genetic diversity in smaller populations indicate the effect of drift. Estimates of gene flow and fish movement suggest that gene flow is limited to nearby populations. In contrast, we found little evidence of IBA as patterns of local ecological and phenotypic variation showed little association with genetic differentiation among populations. Thus, patterns of genetic variation in these small populations likely reflect localized gene flow and genetic drift superimposed onto a larger-scale structure that is largely a result of colonization history. Our simultaneous assessment of the effects of neutral and adaptive processes in a tractable and replicated system has yielded novel insights into the evolution of small populations on both historical and contemporary timescales and over a smaller spatial scale than is typically studied.

Phylogenomic analyses re-examine the evolution of reinforcement and hypothesized hybrid speciation in Phlox wildflowers.

The tree of life is riddled with reticulate evolutionary histories, and some clades, such as the eastern standing Phlox, appear to be hotspots of hybridization. In this group, there are two cases of reinforcement and nine hypothesized hybrid species. Given their historical importance in our understanding of plant speciation, the relationships between these taxa and the role of hybridization in their diversification require genomic validation. Using phylogenomic analyses, we resolve the evolutionary relationships of the eastern standing Phlox and evaluate hypotheses about whether and how hybridization and gene flow played a role in their diversification. Our results provide novel resolution of the phylogenetic relationships in this group, including paraphyly across some taxa. We identify gene flow during one case of reinforcement and find genomic support for a hybrid lineage underlying one of the five hypothesized homoploid hybrid speciation events. Additionally, we estimate the ancestries of four allotetraploid hybrid species. Our results are consistent with hybridization contributing to diverse evolutionary outcomes within this group; although, not as extensively as previously hypothesized. This study demonstrates the importance of phylogenomics in evaluating hypothesized evolutionary histories of non-model systems and adds to the growing support of interspecific genetic exchange in the generation of biodiversity.

Locally adaptive inversions in structured populations.

Inversions have been proposed to facilitate local adaptation, by linking together locally coadapted alleles at different loci. Prior work addressing this question theoretically has considered the spread of inversions in "continent-island" scenarios in which there is a unidirectional flow of maladapted migrants into the island population. In this setting, inversions capturing locally adaptive haplotypes are most likely to invade when selection is weak, because stronger local selection i) more effectively purges maladaptive alleles, and ii) generates linkage disequilibrium between adaptive alleles, thus lessening the advantage of inversions. We show this finding only holds under limited conditions by studying the establishment of inversions in a more general two-deme model, which explicitly considers the dynamics of allele frequencies in both populations linked by bidirectional migration. In this model, the level of symmetry between demes can be varied from complete asymmetry (continent-island) to complete symmetry. For symmetric selection and migration, strong selection increases the allele frequency divergence between demes thereby increasing the frequency of maladaptive alleles in migrants, favouring inversions-the opposite of the pattern seen in the asymmetric continent-island scenario. We also account for the likelihood that a new inversion captures an adaptive haplotype in the first instance. When considering the combined process of capture and invasion in "continent island" and symmetric scenarios, relatively strong selection increases inversion establishment probability. Migration must also be low enough that the inversion is likely to capture an adaptive allele combination, but not so low as to eliminate the inversion's advantage. Overall, our analysis suggests that inversions are likely to harbour larger effect alleles that experience relatively strong selection.

Enhanced and asymmetric signatures of hybridization at climatic margins: Evidence from closely related dioecious fig species.

Hybridization plays a significant role in biological evolution. However, it is not clear whether ecological contingency differentially influences likelihood of hybridization, particularly at ecological margins where parental species may exhibit reduced fitnesses. Moreover, it is unknown whether future ecosystem change will increase the prevalence of hybridization. Ficus heterostyla and F. squamosa are closely related species co-distributed from southern Thailand to southwest China where hybridization, yielding viable seeds, has been documented. As a robust test of ecological factors driving hybridization, we investigated spatial hybridization signatures based on nuclear microsatellites from extensive population sampling across a widespread contact range. Both species showed high population differentiation and strong patterns of isolation by distance. Admixture estimates exposed asymmetric interspecific gene flow. Signatures of hybridization increase significantly towards higher latitude zones, peaking at the northern climatic margins. Geographic variation in reproductive phenology combined with ecologically challenging marginal habitats may promote this phenomenon. Our work is a first systematic evaluation of such patterns in a comprehensive, latitudinally-based clinal context, and indicates that tendency to hybridize appears strongly influenced by environmental conditions. Moreover, that future climate change scenarios will likely alter and possibly augment cases of hybridization at ecosystem scales.

Divergence of the freshwater sleeper, Neodontobutishainanensis (Chen, 1985) (Teleostei, Odontobutidae), in the Pearl River basin and on Hainan Island of southern China.

Study of divergence of freshwater fish populations between island and adjacent mainland areas can shed light on the phylogeographical relationships of these regions. Neodontobutishainanensis is a freshwater fish species restricted to Hainan Island and in Guangdong and Guangxi provinces in the southern mainland China. We examine the phylogenetic relationship and population structure of N.hainanensis based on 3,176 nuclear loci using a gene-capture method. STRUCTURE analysis and principal coordinate analyses (PCA) indicate that populations from Guangdong, Guangxi and Hainan are each distinct, except that some individuals of the Guangdong population share minor genetic components with individuals of the Guangxi population. In the concatenated gene tree, the Hainan population is grouped with the Guangdong population, but the coalescent tree groups the Hainan population as the sister to the Guangxi population. Finally, coalescent simulations confirmed the divergence pattern supported by the coalescent tree and revealed a one-way introgression from the Guangxi population to the Guangdong population, which can explain the discordant results supported by the concatenated and coalescent phylogenetic analyses. Due to recent decline of N.hainanensis populations and the genetic patterns in this species, as revealed in this study, the populations in the three areas should be treated as separate conservation units.

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.

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana.

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

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.

Genetic consequences of improved river connectivity in brown trout (Salmo trutta L.).

Fragmentation of watercourses poses a significant threat to biodiversity, particularly for migratory fish species. Mitigation measures such as fishways, have been increasingly implemented to restore river connectivity and support fish migration. The effects of such restoration efforts are typically tested using telemetry and fisheries methods, which do not fully capture the broader population movements that may have important consequences for population viability. We performed a before-and-after control-impact (BACI) study using genetic tools (SNPs) to investigate the effect of a newly implemented fishway, aiming to enhance upstream spawning migration of brown trout (Salmo trutta Linnaeus) in a reservoir with two headwater tributaries fragmented by man-made weirs. Another reservoir with two barrier-free tributaries was also analysed as a control. Our results showed that the isolated brown trout population was spawning in the reservoir before the installation of the fishway, and we found genetic structuring and differentiation between fragmented headwater tributaries before the fishway construction, but not in the control reservoir. Unexpectedly, after the fishway construction we observed signals consistent with increased genetic differentiation between populations of newly recruited juvenile fish in the reservoir tributary and fish in the reservoir. We propose this was caused by newly enabled philopatric behaviour of brown trout to their natal spawning tributary. In contrast, we did not find any genetic changes in the tributary without a fishway or in the barrier-free reservoir system. Given the scarcity of similar studies, we advocate for an increased use of genetic analyses in BACI studies to monitor and evaluate the effect of efforts to restore habitat connectivity and inform future management strategies.

Climate and habitat type interact to influence contemporary dispersal potential in Prairie Smoke (Geum triflorum).

Understanding dispersal potential, or the probability a species will move a given distance, under different environmental conditions is essential to predicting species' ability to move across the landscape and track shifting ecological niches. Two important drivers of dispersal ability are climatic differences and variations in local habitat type. Despite the likelihood these global drivers act simultaneously on plant populations, and thus dispersal potential is likely to change as a result, their combined effects on dispersal are rarely examined. To understand the effect of climate and varying habitat types on dispersal potential, we studied Geum triflorum-a perennial grassland species that spans a wide range of environments, including both prairie and alvar habitats. We explored how the climate of the growing season and habitat type (prairie vs. alvar) interact to alter dispersal potential. We found a consistent interactive effect of climate and habitat type on dispersal potential. Across prairie populations, an increased number of growing degree days favored traits that increase dispersal potential or the probability of dispersing farther distances. However, for alvar populations, dispersal potential tended to decrease as the number of growing degree days increased. Our findings suggest that under continued warming, populations in prairie habitats will benefit from increased gene flow, while alvar populations will become increasingly segregated, with reduced potential to track shifting fitness optima.

Complex Polyploids: Origins, Genomic Composition, and Role of Introgressed Alleles.

Introgression allows polyploid species to acquire new genomic content from diploid progenitors or from other unrelated diploid or polyploid lineages, contributing to genetic diversity and facilitating adaptive allele discovery. In some cases, high levels of introgression elicit the replacement of large numbers of alleles inherited from the polyploid's ancestral species, profoundly reshaping the polyploid's genomic composition. In such complex polyploids it is often difficult to determine which taxa were the progenitor species and which taxa provided additional introgressive blocks through subsequent hybridization. Here, we use population-level genomic data to reconstruct the phylogenetic history of Betula pubescens (downy birch), a tetraploid species often assumed to be of allopolyploid origin and which is known to hybridize with at least four other birch species. This was achieved by modeling of polyploidization and introgression events under the multispecies coalescent and then using an approximate Bayesian computation (ABC) rejection algorithm to evaluate and compare competing polyploidization models. We provide evidence that B. pubescens is the outcome of an autoploid genome doubling event in the common ancestor of B. pendula and its extant sister species, B. platyphylla, that took place approximately 178,000-188,000 generations ago. Extensive hybridization with B. pendula, B. nana, and B. humilis followed in the aftermath of autopolyploidization, with the relative contribution of each of these species to the B. pubescens genome varying markedly across the species' range. Functional analysis of B. pubescens loci containing alleles introgressed from B. nana identified multiple genes involved in climate adaptation, while loci containing alleles derived from B. humilis revealed several genes involved in the regulation of meiotic stability and pollen viability in plant species.

Genetic Analysis Reveals Dispersal Patterns of Japanese Serow in Two Different Habitats of a Mountainous Region.

Dispersal increases the costs of feeding and predation risk in the new environment and is reported to be biased toward habitats similar to the natal region in some mammals. The benefits and costs of dispersal often differ between sexes, and most mammals show male-biased dispersal in relation to a polygamous mating system. Japanese serow is generally a solitary and monogamous species. However, recent studies have shown that the sociality of serows on Mt. Asama differs between habitat types. In the mountain forests with low forage availability, solitary habits and social monogamy were observed, while, in alpine grasslands, female grouping and social polygyny were observed, which is probably due to abundant forage availability. We investigated the effects of habitat characteristics and sociality on the dispersal of serows using fecal and tissue samples from two different habitats on Mt. Asama. The Fst value between the two areas was significantly positive, and the mean relatedness within areas was significantly higher than that between areas, which suggests limited gene flow and natal habitat-biased dispersal. Bayesian clustering analysis showed unidirectional gene flow from forest to grassland, which was probably due to the high forage availability of the grassland. Analyses of the assignment index and mean relatedness did not show male-biased dispersal, even in the grassland, where serows were polygynous. Thus, polygyny in the grassland is not linked to male-biased dispersal. In summary, our study suggests that dispersal patterns in Japanese serows are affected by habitat rather than social differences.

Speciation within the Anopheles gambiae complex: high-throughput whole genome sequencing reveals evidence of a putative new cryptic taxon in 'far-west' Africa.

The two main Afrotropical malaria vectors - Anopheles coluzzii and An. gambiae - are genetically distinct and reproductively isolated across West Africa. However, populations at the western extreme of their range are assigned as "intermediate" between the two species by whole genome sequence (WGS) data, and as hybrid forms by conventional molecular diagnostics. By exploiting WGS data from 1,190 specimens collected across west Africa via the Anopheles gambiae 1000 Genomes network, we identify a novel putative taxon in the far-west (provisionally named Bissau molecular form), which did not arise by admixture but rather originated at the same time as the split between An. coluzzii and An. gambiae. Intriguingly, these populations lack insecticide resistance mechanisms commonly observed in the two main species. These findings lead to a change of perspective on malaria vector species in the far-west region with potential for epidemiological implications, and a new challenge for genetic-based mosquito control approaches.