Species delimitation using mitochondrial DNA (mtDNA) remains an important and accessible approach for discovering and delimiting species. However, delimiting species with a single locus (e.g. DNA barcoding) is biased towards overestimating species diversity. The highly diverse gecko genus Cyrtodactylus is one such group where delimitation using mtDNA remains the paradigm. In this study, we use genomic data to test putative species boundaries established using mtDNA within three recognized species of Cyrtodactylus on the island of Borneo. We predict that multi-locus genomic data will estimate fewer species than mtDNA, which could have important ramifications for the species diversity within the genus. We aim to (i) investigate the correspondence between species delimitations using mtDNA and genomic data, (ii) infer species trees for each target species, and (iii) quantify gene flow and identify migration patterns to assess population connectivity. We find that species diversity is overestimated and that species boundaries differ between mtDNA and nuclear data. This underscores the value of using genomic data to reassess mtDNA-based species delimitations for taxa lacking clear species boundaries. We expect the number of recognized species within Cyrtodactylus to continue increasing, but, when possible, genomic data should be included to inform more accurate species boundaries.

Biogeographic barriers to gene flow are central to studies of plant phylogeography. There are many physical and geographic barriers in China, but few studies have used molecular ecological evidence to investigate the natural geographic isolation barrier of the Qinling Mountains-Huaihe River Line (QHL). Allium macrostemon is a precious Chinese perennial herb belonging to the Amaryllidaceae family. It is used as a food and medicine, with a variety of health and healing properties. Five SSR markers, three chloroplast DNA (cpDNA) markers (psbA-trnH, rps16 and trnL-F), one nuclear ribosomal DNA (nrDNA) marker (ITS), and simplified genome GBS sequencing were used to analyse the genetic diversity and structure of A. macrostemon. Combining SSR, cpDNA, nrDNA ITS data and GBS analysis results, we divided A. macrostemon populations into northern and southern groups, with the southern group further divided into southwestern and central-southeastern groups. Niche simulation results reveal that the distribution area of A. macrostemon will reach its maximum in the future. These data indicate that the regional separation of A. macrostemon has been maintained by the combined influence of a geographical barrier and Quaternary climate, and that the back-and-forth fluctuations of QHL and Quaternary climate have played an important role in this process. QHL acts as a north-south dividing line in phylogeography and population genetic structure, promoting physical geographic isolation. This study provides a theoretical basis for the conservation, development, and utilization of A. macrostemon resources. It further provides a reference for understanding the systematic geographical pattern of the large-scale spatial distribution of plants in China and enriches our understanding of Quaternary plant evolution in areas with complex terrain.

Coldwater-adapted freshwater fishes, especially their populations along warm-range margins, are most vulnerable to the climate oscillations associated with global warming. Stocking is a major strategy for avoiding the extinction of these species. However, while stocking can reverse the decline of isolated populations, it may also result in a loss of genetic diversity in the native local population due to the introgressive replacement of hatchery genes. To plan an adequate strategy for conserving locally adapted populations, the genetic impacts of stocking on native lineages should be evaluated from small river branches to wide-ranging drainage areas. We investigated the population genetic structure of white-spotted charr (Salvelinus leucomaenis) within its southern range (Lake Biwa basin, Japan). By applying genome-wide SNP analysis to the population\'s genetic structure, we assessed the extent of genetic introgression resulting from stocking. White-spotted charr in the Lake Biwa watershed constitutes a distinctive genetic group, within which apparent genetic differentiation was observed. The hatchery-reared fish line commonly used for supplementation stocking in the catchment was discernable from the native population, enabling us to analyze genetic introgression across the entire drainage area. Admixed individuals resulting from hatchery introgression were observed in most of the stocked sites that showed relatively high heterozygosity and nucleotide diversity. However, their genetic differentiation was much lower than that of native populations. The supplementation history as well as the road availability contributed substantially to the introgression of hatchery genes. Populations with the native genetic structure remained in the upstream regions of the tested rivers. However, their heterozygosity and nucleotide diversity were low when compared with that of the populations with hatchery supplementation. Our results shed light on the genetic impacts of stocking on isolated native populations and suggest that conventional supplementation methods cannot preserve a unique biodiversity in the distribution margin.

Globally, the majority of habitat loss is irreversible, and most species will never recover their former ranges. We have learned a great deal about what leads to population decline and extinction, but less about recovery. The recently downlisted giant panda provides a unique opportunity to understand the mechanisms of species recovery. In our study, we estimate giant panda suitable habitats, population density, and gene flow across landscapes to fully investigate the direct and indirect ecological mechanisms underlying bold conservation strategies. We found that the Giant Panda National Survey has modestly but systematically underestimated population size. China\'s effort to mitigate anthropogenic disturbances was associated with increased panda population density through improving habitat quality and reducing habitat fragmentation. Enhanced landscape connectivity reduced inbreeding via gene flow but indirectly increased inbreeding temporarily due to high local panda density. Although the panda\'s recovery has been geographically uneven, we provide evidence for improving connectivity and gene flow resulting from conservation efforts. If these processes can be sustained and improved, the panda\'s path to recovery will be less encumbered by loss of genetic diversity, fostering hope that the present rate of recovery will not be stalled. Findings from this study will not only help guide future giant panda conservation management but also provide a model for how a more mechanistic examination of the genetic processes underlying species recovery can foster the development of more effective strategies for endangered species recovery.

The time that elapsed between the initial introduction and the proliferation of an invasive species is referred to as the lag phase. The identification of the lag phase is critical for generating plans for pest management and for the prevention of biosecurity failure. However, lag phases have been identified mostly through retrospective searches of historical records. The agricultural pest fall armyworm (FAW; Spodoptera frugiperda) is native to the New World. FAW invasion was first reported from West Africa in 2016, then it spread quickly through Africa, Asia, and Oceania. Here, using population genomics approaches, we demonstrate that the FAW invasion involved an undocumented lag phase. Invasive FAW populations have negative signs of genomic Tajima\'s D, and invasive population-specific genetic variations have particularly decreased Tajima\'s D, supporting a substantial amount of time for the generation of new mutations in introduced FAW populations. Model-based diffusion approximations support the existence of a period with a cessation of gene flow between native and invasive FAW populations. Taken together, these results provide strong support for the presence of a lag phase during the FAW invasion. These results show the usefulness of using population genomics analyses to identify lag phases in biological invasions.

To study the interspecific differentiation characteristics of species originating from recent radiation, the genotyping-by-sequencing (GBS) technique was used to explore the kinship, population structure, gene flow, genetic variability, genotype-environment association and selective sweeps of Picea asperata complex with similar phenotypes from a genome-wide perspective. The following results were obtained: 14 populations of P. asperata complex could be divided into 5 clades; P. wilsonii and P. neoveitchii diverged earlier and were more distantly related to the remaining 6 spruce species. Various geological events have promoted the species differentiation of P. asperata complex. There were four instances of gene flow among P. koraiensis, P. meyeri, P. asperata, P. crassifolia and P. mongolica. The population of P. mongolica had the highest level of nucleotide diversity, and P. neoveitchii may have experienced a bottleneck recently. Genotype-environment association found that a total of 20,808 genes were related to the environmental variables, which enhanced the adaptability of spruce in different environments. Genes that were selectively swept in the P. asperata complex were primarily associated with plant stress resistance. Among them were some genes involved in plant growth and development, heat stress, circadian rhythms and flowering. In addition to the commonly selected genes, different spruce species also displayed unique genes subjected to selective sweeps that improved their adaptability to different habitats. Understanding the interspecific gene flow and adaptive evolution of Picea species is beneficial to further understanding the species relationships of spruce and can provide a basis for studying spruce introgression and functional genomics.

Designated silkworms (Antheraea mylitta D.) are found in the Smilipal Biosphere Reserve (SBR) and consist of four ecological groups. This study examined the genetics of eight ethnic groups from four ecological groups (Modal, Nalia, Bogai and Jata) using RAPD and ISSR markers. The RAPD primer produced 127 bands, of which 41 out of 116 polymorphic regions had a PIC value above 0.45. ISSR markers assigned PIC data values to 57 of 127 polymorphic sites. Cluster analysis revealed genetic differences between ecological populations, with markers separating ecological populations into groups. Low testicles and high genetic diversity indicate that these individuals are at different levels. If adequate protection is not provided, the risk of extinction of ecospecies designed in SBR increases. This study shows that the combined use of RAPD + ISSR markers can improve the results of the technology in identifying different diseases in situations where genomic information is limited.

The Arctic is warming four times faster than the rest of the world, threatening the persistence of many Arctic species. It is uncertain if Arctic wildlife will have sufficient time to adapt to such rapidly warming environments. We used genetic forecasting to measure the risk of maladaptation to warming temperatures and sea ice loss in polar bears (Ursus maritimus) sampled across the Canadian Arctic. We found evidence for local adaptation to sea ice conditions and temperature. Forecasting of genome-environment mismatches for predicted climate scenarios suggested that polar bears in the Canadian high Arctic had the greatest risk of becoming maladapted to climate warming. While Canadian high Arctic bears may be the most likely to become maladapted, all polar bears face potentially negative outcomes to climate change. Given the importance of the sea ice habitat to polar bears, we expect that maladaptation to future warming is already widespread across Canada.

Many cooperatively breeding species live in groups with complex structure-large group sizes, low and variable kin structure, and multiple breeding pairs. Since these mixed-kin groups typically form because of immigration of unrelated individuals of both sexes in addition to limited offspring dispersal, differences in patterns of dispersal can generate variation in group structure, even within the same species or population. Here, we examine how environmentally mediated dispersal patterns influence variation in group structure in the plural breeding superb starling (Lamprotornis superbus), an avian cooperative breeder that inhabits a spatiotemporally variable savanna environment and forms mixed-kin groups with variable group sizes and more than one breeding pair per group. Using 4068 genome-wide polymorphic loci and fine-scale, remotely sensed ecological data from 22 groups sampled across a nearly 200 km2 environmental gradient in central Kenya, we find evidence of not only frequent and long-distance dispersal in both sexes (low isolation-by-distance and weak genetic structure), but also directional dispersal from small groups in lower quality habitat with low normalised difference vegetation index (NDVI) to large groups in higher quality habitat with high NDVI. Additionally, we find stronger genetic structure among groups in lower quality habitat, and higher genetic diversity and lower relatedness of groups in higher quality habitat. Previous work using long-term data from groups in the same population has shown that groups with lower relatedness are larger and have more breeding pairs. Long-distance, directional dispersal to maximise individual fitness can thus lead to smaller and simpler kin-based social groups in lower quality habitat, but larger and more complex mixed-kin groups in higher quality habitat. Such intraspecific, within-population variation in group structure, including variation in kin structure of social groups, could have profound implications for the relative importance of the evolutionary mechanisms (i.e. direct vs. indirect fitness benefits) underlying the formation of cooperative societies.

Selaginellaceae, originated in the Carboniferous and survived the Permian-Triassic mass extinction, is the largest family of lycophyte, which is sister to other tracheophytes. It stands out from tracheophytes by exhibiting extraordinary habitat diversity and lacking polyploidization. The organelle genome-based phylogenies confirmed the monophyly of Selaginella, with six or seven subgenera grouped into two superclades, but the phylogenetic positions of the enigmatic Selaginella sanguinolenta clade remained problematic. Here, we conducted a phylogenomic study on Selaginellaceae utilizing large-scale nuclear gene data from RNA-seq to elucidate the phylogeny and explore the causes of the phylogenetic incongruence of the S. sanguinolenta clade. Our phylogenetic analyses resolved three different positions of the S. sanguinolenta clade, which were supported by the sorted three nuclear gene sets, respectively. The results from the gene flow test, species network inference, and plastome-based phylogeny congruently suggested a probable hybrid origin of the S. sanguinolenta clade involving each common ancestor of the two superclades in Selaginellaceae. The hybrid hypothesis is corroborated by the evidence from rhizophore morphology and spore micromorphology. The chromosome observation and Ks distributions further suggested hybridization accompanied by polyploidization. Divergence time estimation based on independent datasets from nuclear gene sets and plastid genome data congruently inferred that allopolyploidization occurred in the Early Triassic. To our best knowledge, the allopolyploidization in the Mesozoic reported here represents the earliest record of tracheophytes. Our study revealed a unique triad of phylogenetic positions for a hybrid-originated group with comprehensive evidence and proposed a hypothesis for retaining both parental alleles through gene conversion.