Species divergence and environmental adaptation of Picea asperata complex at the whole genome level.

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.

Signatures of selective sweeps in continuous-space populations.

Selective sweeps describe the process by which an adaptive mutation arises and rapidly fixes in the population, thereby removing genetic variation in its genomic vicinity. The expected signatures of selective sweeps are relatively well understood in panmictic population models, yet natural populations often extend across larger geographic ranges where individuals are more likely to mate with those born nearby. To investigate how such spatial population structure can affect sweep dynamics and signatures, we simulated selective sweeps in populations inhabiting a two-dimensional continuous landscape. The maximum dispersal distance of offspring from their parents can be varied in our simulations from an essentially panmictic population to scenarios with increasingly limited dispersal. We find that in low-dispersal populations, adaptive mutations spread more slowly than in panmictic ones, while recombination becomes less effective at breaking up genetic linkage around the sweep locus. Together, these factors result in a trough of reduced genetic diversity around the sweep locus that looks very similar across dispersal rates. We also find that the site frequency spectrum around hard sweeps in low-dispersal populations becomes enriched for intermediate-frequency variants, making these sweeps appear softer than they are. Furthermore, haplotype heterozygosity at the sweep locus tends to be elevated in low-dispersal scenarios as compared to panmixia, contrary to what we observe in neutral scenarios without sweeps. The haplotype patterns generated by these hard sweeps in low-dispersal populations can resemble soft sweeps from standing genetic variation that arose from substantially older alleles. Our results highlight the need for better accounting for spatial population structure when making inferences about selective sweeps.

Adaptive alien genes are maintained amid a vanishing introgression footprint in a sea squirt.

Human transport of species across oceans disrupts natural dispersal barriers and facilitates hybridization between previously allopatric species. The recent introduction of the North Pacific sea squirt, Ciona robusta, into the native range of the North Atlantic sea squirt, Ciona intestinalis, is a good example of this outcome. Recent studies have revealed an adaptive introgression in a single chromosomal region from the introduced into the native species. Here, we monitored this adaptive introgression over time, examining both the frequency of adaptive alleles at the core and the hitchhiking footprint in the shoulders of the introgression island by studying a thousand Ciona spp. individuals collected in 22 ports of the contact zone, 14 of which were sampled 20 generations apart. For that purpose, we developed a KASP multiplex genotyping approach, which proved effective in identifying native, nonindigenous and hybrid individuals and in detecting introgressed haplotypes. We found no early generation hybrids in the entire sample, and field observations suggest a decline in the introduced species. At the core region of the introgression sweep, where the frequency of C. robusta alleles is the highest and local adaptation genes must be, we observed stable frequencies of adaptive alien alleles in both space and time. In contrast, we observed erosion of C. robusta ancestry tracts in flanking chromosomal shoulders on the edges of the core, consistent with the second phase of a local sweep and a purge of hitchhiked incompatible mutations. We hypothesize that adaptive introgression may have modified the competition relationships between the native and invasive species in human-altered environments.

Rapid Adaptation and Interspecific Introgression in the North American Crop Pest Helicoverpa zea.

Insect crop pests threaten global food security. This threat is amplified through the spread of nonnative species and through adaptation of native pests to control measures. Adaptations such as pesticide resistance can result from selection on variation within a population, or through gene flow from another population. We investigate these processes in an economically important noctuid crop pest, Helicoverpa zea, which has evolved resistance to a wide range of pesticides. Its sister species Helicoverpa armigera, first detected as an invasive species in Brazil in 2013, introduced the pyrethroid-resistance gene CYP337B3 to South American H. zea via adaptive introgression. To understand whether this could contribute to pesticide resistance in North America, we sequenced 237 H. zea genomes across 10 sample sites. We report H. armigera introgression into the North American H. zea population. Two individuals sampled in Texas in 2019 carry H. armigera haplotypes in a 4 Mbp region containing CYP337B3. Next, we identify signatures of selection in the panmictic population of nonadmixed H. zea, identifying a selective sweep at a second cytochrome P450 gene: CYP333B3. We estimate that its derived allele conferred a ∼5% fitness advantage and show that this estimate explains independently observed rare nonsynonymous CYP333B3 mutations approaching fixation over a ∼20-year period. We also detect putative signatures of selection at a kinesin gene associated with Bt resistance. Overall, we document two mechanisms of rapid adaptation: the introduction of fitness-enhancing alleles through interspecific introgression, and selection on intraspecific variation.

Phenotypic Identification, Genetic Characterization, and Selective Signal Detection of Huitang Duck.

The Huitang duck (HT), a long-domesticated elite local breed from Hunan Province, China, with excellent meat quality, has not had its population genetic structure and genomic selective sweeps extensively studied to date. This study measured the phenotypic characteristics of HT and conducted comparative analysis between HT and 16 different duck breeds, including wild, indigenous, and meat breeds, to characterize its population structure and genetic potential. The results revealed that HT is a dual-purpose indigenous breed with a genetic background closely related to the Youxian sheldrake and Linwu ducks. In the selective sweep analysis between HT and Linwu ducks, genes such as PLCG2, FN1, and IGF2BP2, which are associated with muscle growth and development, were identified near the 27 selection signals. The comparison between HT and Jinding ducks revealed 68 selective signals that contained important genes associated with ovarian development (GRIK4, MAP3K8, and TGIF1) and egg-laying behaviors (ERBB4). Selective sweep analysis between HT and Youxian sheldrake ducks found 93 selective regions covering genes related to both meat (IGF1R and IGFBP5) and egg-production (FOXO3 and ITPR1) traits. Our study may provide novel knowledge for exploring the population structure and genetic potential of HT, offering a theoretical basis for its breeding strategies in the future.

Whole-genome resequencing reveals genetic diversity and adaptive evolution in Chinese honeybee (Apis cerana cerana) in Guizhou, China.

Introduction: Guizhou Province, characterized by complex and diverse geographic and climatic environments, has rich genetic resources for the Chinese honeybee (Apis cerana cerana) and is one of the main bee-producing areas in China. However, research on the genetic diversity of Chinese honeybee in the Guizhou region is very limited, despite implications for conservation of biodiversity. Methods: In this study, we analyzed the genetic diversity, differentiation, and selection signals based on 116 Chinese honeybees from 12 regions in Guizhou Province using whole-genome sequencing. Results: We identified 1,400,430 high-quality SNPs across all samples. A population structure analysis revealed two independent genetic subgroups of Chinese honeybees in Guizhou, a Yunnan-Guizhou Plateau population in western Guizhou and a hilly-mountainous population in eastern Guizhou. The average nucleotide diversity (Pi) ranged from 0.00138 to 0.00161 and average expected heterozygosity (He) ranged from 0.2592 to 0.2604. The average genetic differentiation index (F ST) for Chinese honeybees in pairwise comparisons of 12 regions ranged from 0.0094 to 0.0293. There was clear genetic differentiation between the western plateau and the eastern hilly mountainous areas of Guizhou; however, F ST values between the eastern and western populations ranged from 0.0170 to 0.0293, indicating a low degree of differentiation. A genome-wide scan revealed a number of genes under selection in the Yunnan-Guizhou Plateau environment. These genes were related to growth and development, reproduction, and cold resistance, and several candidate genes involved in environmental adaptation were identified, including CTR, MAPK, MAST, HSF, and MKKK. Discussion: The results of the present study provide important theoretical bases for the conservation, evaluation, development, and utilization of genetic resources for Chinese honeybees in the Guizhou region and for further investigations of environmental adaptation and underlying mechanisms in the species.

A comparative analysis reveals the genomic diversity among 8 Muscovy duck populations.

The Muscovy duck (Cairina moschata) is a waterfowl indigenous to the neotropical regions of Central and South America. It has low demand for concentrated feed and strong adaptability to different rearing conditions. After introduced to China through Eurasian commercial trade, Muscovy ducks have a domestication history of around 300 years in the Fujian Province of China. In the 1990s, the commodity Muscovy duck breed "Crimo," cultivated in Europe, entered the Chinese market for consumption and breeding purposes. Due to the different selective breeding processes, Muscovy ducks have various populational traits and lack transparency of their genetic background. To remove this burden in the Muscovy duck breeding process, we analyzed genomic data from 8 populations totaling 83 individuals. We identify 11.24 million single nucleotide polymorphisms (SNPs) and categorized these individuals into the Fujian-bred and the Crimo populations according to phylogenetic analyses. We then delved deeper into their evolutionary relationships through assessing population structure, calculating fixation index (FST) values, and measuring genetic distances. Our exploration of runs of homozygosity (ROHs) and homozygous-by-descent (HBD) uncovered genomic regions enriched for genes implicated in fatty acid metabolism, development, and immunity pathways. Selective sweep analyses further indicated strong selective pressures exerted on genes including TECR, STAT2, and TRAF5. These findings provide insights into genetic variations of Muscovy ducks, thus offering valuable information regarding genetic diversity, population conservation, and genome associated with the breeding of Muscovy ducks.

Natural selection on feralization genes contributed to the invasive spread of wild pigs throughout the United States.

Despite a long presence in the contiguous United States (US), the distribution of invasive wild pigs (Sus scrofa × domesticus) has expanded rapidly since the 1980s, suggesting a more recent evolutionary shift towards greater invasiveness. Contemporary populations of wild pigs represent exoferal hybrid descendants of domestic pigs and European wild boar, with such hybridization expected to enrich genetic diversity and increase the adaptive potential of populations. Our objective was to characterize how genetic enrichment through hybridization increases the invasiveness of populations by identifying signals of selection and the ancestral origins of selected loci. Our study focused on invasive wild pigs within Great Smoky Mountains National Park, which represents a hybrid population descendent from the admixture of established populations of feral pigs and an introduction of European wild boar to North America. Accordingly, we genotyped 881 wild pigs with multiple high-density single-nucleotide polymorphism (SNP) arrays. We found 233 markers under putative selection spread over 79 regions across 16 out of 18 autosomes, which contained genes involved in traits affecting feralization. Among these, genes were found to be related to skull formation and neurogenesis, with two genes, TYRP1 and TYR, also encoding for crucial melanogenesis enzymes. The most common haplotypes associated with regions under selection for the Great Smoky Mountains population were also common among other populations throughout the region, indicating a key role of putatively selective variants in the fitness of invasive populations. Interestingly, many of these haplotypes were absent among European wild boar reference genotypes, indicating feralization through genetic adaptation.

Whole-genome resequencing deciphers patterns of genetic diversity, phylogeny, and evolutionary dynamics in Kashmir cattle.

Kashmir cattle, which were kept by local pastoralists for centuries, are exceptionally resilient and adaptive to harsh environments. Despite its significance, the genomic characteristics of this cattle breed remain elusive. This study utilized whole genome sequences of Kashmir cattle (n = 20; newly sequenced) alongside published whole genomes of 32 distinct breeds and seven core cattle populations (n = 135). The analysis identified ~25.87 million biallelic single nucleotide polymorphisms in Kashmir cattle, predominantly in intergenic and intron regions. Population structure analyses revealed distinct clustering patterns of Kashmir cattle with proximity to the South Asian, African and Chinese indicine cattle populations. Genetic diversity analysis of Kashmir cattle demonstrated lower inbreeding and greater nucleotide diversity than analyzed global breeds. Homozygosity runs indicated less consanguineous mating in Kashmir cattle compared with European taurine breeds. Furthermore, six selection sweep detection methods were used within Kashmir cattle and other cattle populations to identify genes associated with vital traits, including immunity (BOLA-DQA5, BOLA-DQB, TNFAIP8L, FCRL4, AOAH, HIF1AN, FBXL3, MPEG1, CDC40, etc.), reproduction (GOLGA4, BRWD1, OSBP2, LEO1 ADCY5, etc.), growth (ADPRHL1, NRG2, TCF12, TMOD4, GBP4, IGF2, RSPO3, SCD, etc.), milk composition (MRPS30 and CSF1) and high-altitude adaptation (EDNRA, ITPR2, AGBL4 and SCG3). These findings provide essential genetic insights into the characteristics and establish the foundation for the scientific conservation and utilization of Kashmir cattle breed.

Whole-genome resequencing uncovers diversity and selective sweep in Kazakh cattle.

The Kazakh cattle in the Xinjiang Uygur Autonomous Region of China are highly adaptable and have multiple uses, including milk and meat production, and use as draft animals. They are an excellent original breed that could be enhanced by breeding and hybrid improvement. However, the genomic diversity and signature of selection underlying the germplasm characteristics require further elucidation. Herein, we evaluated 26 Kazakh cattle genomes in comparison with 103 genomes of seven other cattle breeds from regions around the world to assess the Kazakh cattle genetic variability. We revealed that the relatively low linkage disequilibrium at large SNP distances was strongly correlated with the largest effective population size among Kazakh cattle. Using population structural analysis, we next demonstrated a taurine lineage with restricted Bos indicus introgression among Kazakh cattle. Notably, we identified putative selected genes associated with resistance to disease and body size within Kazakh cattle. Together, our findings shed light on the evolutionary history and breeding profile of Kazakh cattle, as well as offering indispensable resources for germplasm resource conservation and crossbreeding program implementation.

Artificial selection footprints in indigenous and commercial chicken genomes.

BACKGROUND: Although many studies have been done to reveal artificial selection signatures in commercial and indigenous chickens, a limited number of genes have been linked to specific traits. To identify more trait-related artificial selection signatures and genes, we re-sequenced a total of 85 individuals of five indigenous chicken breeds with distinct traits from Yunnan Province, China. RESULTS: We found 30 million non-redundant single nucleotide variants and small indels (< 50 bp) in the indigenous chickens, of which 10 million were not seen in 60 broilers, 56 layers and 35 red jungle fowls (RJFs) that we compared with. The variants in each breed are enriched in non-coding regions, while those in coding regions are largely tolerant, suggesting that most variants might affect cis-regulatory sequences. Based on 27 million bi-allelic single nucleotide polymorphisms identified in the chickens, we found numerous selective sweeps and affected genes in each indigenous chicken breed and substantially larger numbers of selective sweeps and affected genes in the broilers and layers than previously reported using a rigorous statistical model. Consistent with the locations of the variants, the vast majority (~ 98.3%) of the identified selective sweeps overlap known quantitative trait loci (QTLs). Meanwhile, 74.2% known QTLs overlap our identified selective sweeps. We confirmed most of previously identified trait-related genes and identified many novel ones, some of which might be related to body size and high egg production traits. Using RT-qPCR, we validated differential expression of eight genes (GHR, GHRHR, IGF2BP1, OVALX, ELF2, MGARP, NOCT, SLC25A15) that might be related to body size and high egg production traits in relevant tissues of relevant breeds. CONCLUSION: We identify 30 million single nucleotide variants and small indels in the five indigenous chicken breeds, 10 million of which are novel. We predict substantially more selective sweeps and affected genes than previously reported in both indigenous and commercial breeds. These variants and affected genes are good candidates for further experimental investigations of genotype-phenotype relationships and practical applications in chicken breeding programs.

Genomic structure analysis and construction of DNA fingerprint for four sheep populations.

The Yongdeng Qishan sheep (QS) is a sheep population found locally in China. To gain in-depth knowledge of its population characteristics, three control groups were chosen, comprising the Lanzhou fat-tailed sheep (LFT), TAN sheep (TAN), and Minxian black fur sheep (MBF), inhabiting the nearby environments. This study genotyped a total of 120 individuals from four sheep populations: QS, LFT, TAN, and MBF. Using Specific-Locus Amplified Fragment Sequencing, we conducted genetic diversity, population structure, and selective sweep analysis, and constructed the fingerprint of each population. In total, there were 782 535 single nucleotide polymorphism (SNP) variations identified, with most being situated within regions that are intergenic or intronic. The genetic diversity analysis revealed that the QS population exhibited lower genetic diversity compared to the other three populations. Consistent results were obtained from the principal component, phylogenetic tree, and population structure analysis, indicating significant genetic differences between QS and the other three populations. However, a certain degree of differentiation was observed within the QS population. The linkage disequilibrium (LD) patterns among the four populations showed clear distinctions, with the QS group demonstrating the most rapid LD decline. Kinship analysis supported the findings of population structure, dividing the 90 QS individuals into two subgroups consisting of 23 and 67 individuals. Selective sweep analysis identified a range of genes associated with reproduction, immunity, and adaptation to high-altitude hypoxia. These genes hold potential as candidate genes for marker-assisted selection breeding. Additionally, a total of 86 523 runs of homozygosity (ROHs) were detected, showing non-uniform distribution across chromosomes, with chromosome 1 having the highest coverage percentage and chromosome 26 the lowest. In the high-frequency ROH islands, 79 candidate genes were associated with biological processes such as reproduction and fat digestion and absorption. Furthermore, a DNA fingerprint was constructed for the four populations using 349 highly polymorphic SNPs. In summary, our research delves into the genetic diversity and population structure of QS population. The construction of DNA fingerprint profiles for each population can provide valuable references for the identification of sheep breeds both domestically and internationally.

Genomic Profiling of Insecticide Resistance in Malaria Vectors: Insights into Molecular Mechanisms.

Malaria control faces challenges from widespread insecticide resistance in major Anopheles species. This study, employing a cross-species approach, integrates RNA-Sequencing, whole-genome sequencing, and microarray data to elucidate drivers of insecticide resistance in Anopheles gambiae complex and An. funestus. Findings show an inverse relationship between genetic diversity and gene expression, with highly expressed genes experiencing stronger purifying selection. These genes cluster physically in the genome, revealing potential coordinated regulation. We identified known and novel candidate insecticide resistance genes, enriched in metabolic, cuticular, and behavioural functions. We also present AnoExpress, a Python package, and an online interface for user-friendly exploration of resistance candidate expression. Despite millions of years of speciation, convergent gene expression responses to insecticidal selection pressures are observed across Anopheles species, providing crucial insights for malaria vector control. This study culminates in a rich dataset that allows us to understand molecular mechanisms, better enabling us to combat insecticide resistance effectively.

A Novel Approach Utilizing Domain Adversarial Neural Networks for the Detection and Classification of Selective Sweeps.

The identification and classification of selective sweeps are of great significance for improving the understanding of biological evolution and exploring opportunities for precision medicine and genetic improvement. Here, a domain adaptation sweep detection and classification (DASDC) method is presented to balance the alignment of two domains and the classification performance through a domain-adversarial neural network and its adversarial learning modules. DASDC effectively addresses the issue of mismatch between training data and real genomic data in deep learning models, leading to a significant improvement in its generalization capability, prediction robustness, and accuracy. The DASDC method demonstrates improved identification performance compared to existing methods and excels in classification performance, particularly in scenarios where there is a mismatch between application data and training data. The successful implementation of DASDC in real data of three distinct species highlights its potential as a useful tool for identifying crucial functional genes and investigating adaptive evolutionary mechanisms, particularly with the increasing availability of genomic data.

Integrated analysis strategy of genome-wide functional gene mining reveals DKK2 gene underlying meat quality in Shaziling synthesized pigs.

BACKGROUND: Shaziling pig is a well-known indigenous breed in China who has superior meat quality traits. However, the genetic mechanism and genomic evidence underlying meat quality characteristics of Shaziling pigs are still unclear. To explore and investigate the germplasm characteristics of Shaziling pigs, we totally analyzed 67 individual's whole genome sequencing data for the first time (20 Shaziling pigs [S], 20 Dabasha pigs [DBS], 11 Yorkshire pigs [Y], 10 Berkshire pigs [BKX], 5 Basha pigs [BS] and 1 Warthog). RESULTS: A total of 2,538,577 SNPs with high quality were detected and 9 candidate genes which was specifically selected in S and shared in S to DBS were precisely mined and screened using an integrated analysis strategy of identity-by-descent (IBD) and selective sweep. Of them, dickkopf WNT signaling pathway inhibitor 2 (DKK2), the antagonist of Wnt signaling pathway, was the most promising candidate gene which was not only identified an association of palmitic acid and palmitoleic acid quantitative trait locus in PigQTLdb, but also specifically selected in S compared to other 48 Chinese local pigs of 12 populations and 39 foreign pigs of 4 populations. Subsequently, a mutation at 12,726-bp of DKK2 intron 1 (g.114874954 A > C) was identified associated with intramuscular fat content using method of PCR-RFLP in 21 different pig populations. We observed DKK2 specifically expressed in adipose tissues. Overexpression of DKK2 decreased the content of triglyceride, fatty acid synthase and expression of relevant genes of adipogenic and Wnt signaling pathway, while interference of DKK2 got contrary effect during adipogenesis differentiation of porcine preadipocytes and 3T3-L1 cells. CONCLUSIONS: Our findings provide an analysis strategy for mining functional genes of important economic traits and provide fundamental data and molecular evidence for improving pig meat quality traits and molecular breeding.

Whole-genome resequencing of major populations revealed domestication-related genes in yaks.

BACKGROUND: The yak is a symbol of the Qinghai-Tibet Plateau and provides important basic resources for human life on the plateau. Domestic yaks have been subjected to strong artificial selection and environmental pressures over the long-term. Understanding the molecular mechanisms of phenotypic differences in yak populations can reveal key functional genes involved in the domestication process and improve genetic breeding. MATERIAL AND METHOD: Here, we re-sequenced 80 yaks (Maiwa, Yushu, and Huanhu populations) to identify single-nucleotide polymorphisms (SNPs) as genetic variants. After filtering and quality control, remaining SNPs were kept to identify the genome-wide regions of selective sweeps associated with domestic traits. The four methods (π, XPEHH, iHS, and XP-nSL) were used to detect the population genetic separation. RESULTS: By comparing the differences in the population stratification, linkage disequilibrium decay rate, and characteristic selective sweep signals, we identified 203 putative selective regions of domestic traits, 45 of which were mapped to 27 known genes. They were clustered into 4 major GO biological process terms. All known genes were associated with seven major domestication traits, such as dwarfism (ANKRD28), milk (HECW1, HECW2, and OSBPL2), meat (SPATA5 and GRHL2), fertility (BTBD11 and ARFIP1), adaptation (NCKAP5, ANTXR1, LAMA5, OSBPL2, AOC2, and RYR2), growth (GRHL2, GRID2, SMARCAL1, and EPHB2), and the immune system (INPP5D and ADCYAP1R1). CONCLUSIONS: We provided there is an obvious genetic different among domestic progress in these three yak populations. Our findings improve the understanding of the major genetic switches and domestic processes among yak populations.

Shared evolutionary processes shape landscapes of genomic variation in the great apes.

For at least the past 5 decades, population genetics, as a field, has worked to describe the precise balance of forces that shape patterns of variation in genomes. The problem is challenging because modeling the interactions between evolutionary processes is difficult, and different processes can impact genetic variation in similar ways. In this paper, we describe how diversity and divergence between closely related species change with time, using correlations between landscapes of genetic variation as a tool to understand the interplay between evolutionary processes. We find strong correlations between landscapes of diversity and divergence in a well-sampled set of great ape genomes, and explore how various processes such as incomplete lineage sorting, mutation rate variation, GC-biased gene conversion and selection contribute to these correlations. Through highly realistic, chromosome-scale, forward-in-time simulations, we show that the landscapes of diversity and divergence in the great apes are too well correlated to be explained via strictly neutral processes alone. Our best fitting simulation includes both deleterious and beneficial mutations in functional portions of the genome, in which 9% of fixations within those regions is driven by positive selection. This study provides a framework for modeling genetic variation in closely related species, an approach which can shed light on the complex balance of forces that have shaped genetic variation.

Two-speed genomes of Epichloe fungal pathogens show contrasting signatures of selection between species and across populations.

Antagonistic selection between pathogens and their hosts can drive rapid evolutionary change and leave distinct molecular footprints of past and ongoing selection in the genomes of the interacting species. Despite an increasing availability of tools able to identify signatures of selection, the genetic mechanisms underlying coevolutionary interactions and the specific genes involved are still poorly understood, especially in heterogeneous natural environments. We searched the genomes of two species of Epichloe plant pathogen for evidence of recent selection. The Epichloe genus includes highly host-specific species that can sterilize their grass hosts. We performed selection scans using genome-wide SNP data from seven natural populations of two co-occurring Epichloe sibling species specialized on different hosts. We found evidence of recent (and ongoing) selective sweeps across the genome in both species. However, selective sweeps were more abundant in the species with a larger effective population size. Sweep regions often overlapped with highly polymorphic AT-rich regions supporting the role of these genome compartments in adaptive evolution. Although most loci under selection were specific to individual populations, we could also identify several candidate genes targeted by selection in sweep regions shared among populations. The genes encoded small secreted proteins typical of fungal effectors and cell wall-degrading enzymes. By investigating the genomic signatures of selection across multiple populations and species, this study contributes to our understanding of complex adaptive processes in natural plant pathogen systems.

Research Note: Development and application of specific molecular identity cards for "Yufen 1" H line chickens.

The "Yufen 1" H line chicken (YF) has excellent characteristics including early sexual maturity and high egg production, and the conservation of its genetic diversity is the core of the breeding activity. To overcome misrepresented breeds and protect the integrity of the germplasm genetic resources, it is important to develop accurate and convenient methods to identify YF. In this study, whole genome resequencing was performed on the YF population, and bioinformatics analysis was conducted by combining the data from different breeds. Linkage disequilibrium (LD) analysis revealed that YF had the slowest LD-decay rate, suggesting strong natural and artificial selection in its history. Through selective sweep analysis, 1,126 selected regions in YF were identified, which contained 163,661 single nucleotide polymorphisms (SNPs). In particular, 5 specific SNPs (SNP1: Chr2:45509616, SNP2: Chr2:45510792, SNP3: Chr9:13788193, SNP4: Chr9:13795646, SNP5: Chr9:13798154) were found exclusively in the YF population. Subsequently, PCR amplification and Sanger sequencing confirmed the presence of these 5 SNPs in YF. Finally, 4 SNPs (SNP1, SNP2, SNP4, SNP5) were screened and verified using the Kompetitive Allele Specific PCR (KASP) typing technique. These SNPs can be used as specific molecular identity cards (IDs) for YF authentication. The present study is of great significance to ensure sustainable conservation and promotion of YF germplasm resources.