SEARCH · Search PubMed

Results for “nucleotide diversity”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

PubMed2026

Genetic rescue stabilizes diversity in small isolated populations of Bonneville cutthroat trout.

Genetic diversity loss due to anthropogenic factors is occurring rapidly on a global scale, putting many species at risk of extirpation and extinction. Different management strategies have been developed to slow this loss; however, it is often unknown whether these strategies reach their intended goals. In this study, we evaluate population structure and changes in nucleotide diversity (π) in isolated populations of Bonneville cutthroat trout (Oncorhynchus clarkii utah) from the Snake Range (Nevada, USA). Starting in the 1990s, three of these populations were used to reestablish populations in the Snake Range because many of the historic populations were extirpated. Some populations were stocked using a single-source and others were stocked using multiple-sources. Using low-coverage whole-genome sequencing coupled with historic samples (2003-2010) and contemporary samples (2019-2022), we find that single-source populations lost nucleotide diversity while mixed-source populations maintained nucleotide diversity. Further, source populations used to restore populations throughout the Snake Range lost the most nucleotide diversity over the time span evaluated. Our findings provide insight into how small, isolated populations can be managed to maintain genetic diversity.

Animals↗

PubMed2026

Targeted population genomics uncovers demographic history and genetic divergence in north American wild cranberry.

Wild populations of North American cranberry (Vaccinium macrocarpon Aiton) are reservoirs of genetic variation that may contribute to the improvement of breeding-relevant traits. However, the extent to which wild genetic variation is geographically structured and represented in elite germplasm remains unclear. We analysed 179 wild cranberry accessions from the upper Midwest and Eastern North America to estimate nucleotide diversity (π), population structure, and loci associated with genetic differentiation and environmental variables using a genome-informed targeted genotyping panel. Additionally, 14 demographic scenarios were evaluated using site-frequency-spectrum-based inference to identify historical events that could explain current genetic diversity. We observed extremely low nucleotide diversity within the targeted panel (π = 5 × 10-6). Rare allele distributions strongly influenced π and Tajima's D values, suggesting constrained diversity in the genomic regions assayed that is not captured by heterozygosity-based estimates alone. However, we interpreted these results as conservative lower bounds on genome-wide neutral diversity because the targeted panel is enriched for genic and conserved regions. A clear separation between the Midwest and East populations was observed, with inbreeding coefficients ranging from -0.13 to 0.15. Furthermore, site frequency spectrum inference from the targeted panel supported a demographic scenario consistent with a significant population reduction ≈15-14 thousand years ago (kya), followed by a divergence between the two regions ≈12 kya, and an asymmetric gene flow ≈1.3 kya. We detected 254 candidate loci showing regional allele-frequency differentiation. Several of these loci colocalized with candidate genes linked to stress response, development, and metabolic processes. To evaluate the representation of geographically differentiated wild alleles in a breeding context, we analysed Rutgers breeding materials (n = 484) and found that this panel is enriched for common alleles in Eastern wild populations. These findings indicate regionally structured allele-frequency variation in wild cranberry, with potential relevance to environmental response and breeding. This study extends prior wild cranberry population-genetic research by providing targeted-panel estimates of diversity, comparisons of demographic models, and breeding insights on geographically differentiated alleles, while highlighting the importance of conserving wild cranberry germplasm for use in modern breeding programs.

Journal Article↗

PubMed2026

Whole-Genome Sequencing Reveals Population Structure, Genetic Diversity, and Selection Signatures in Kazakh Dromedary and Bactrian Camels.

Understanding the genomic basis of environmental adaptation is essential for the conservation and genetic improvement of domestic camels. In this study, we investigated the population structure, genetic diversity, and genomic variation potentially associated with environmental adaptation of Kazakh dromedary and Bactrian camels using whole-genome sequencing. Whole-genome sequencing data were generated for Kazakh camels (15 dromedaries and 16 Bactrian camels) and integrated with 131 publicly available genomes representing camel populations from the Arabian Peninsula, Iran, Xinjiang, Inner Mongolia, and Mongolian wild camels. Population structure, genetic diversity, and genome-wide selection were evaluated using principal component analysis, ADMIXTURE, nucleotide diversity, linkage disequilibrium, runs of homozygosity, genomic inbreeding (FROH), and selection scans based on FST, θπ ratio, and XP-EHH. Population genomic analyses revealed clear differentiation between dromedary and Bactrian camels, whereas Kazakh camel populations exhibited higher nucleotide diversity (θπ = 1.307-1.551 × 10-3), and lower genomic inbreeding (median FROH: 0.037-0.056) than Arabian populations. Genome-wide selection analyses identified MC4R as the prominent candidate gene in Kazakh dromedaries and RYR1 as a prominent candidate gene in Kazakh Bactrian camels. Functional enrichment analyses highlighted pathways related to energy metabolism, thermogenesis, calcium signaling, skeletal muscle function, mitochondrial activity, and oxidative stress response. These findings provide new insights into genomic variation potentially associated with environmental adaptation in Kazakh camels and offer valuable genomic resources for future conservation, breeding, and evolutionary studies.

MC4R↗

PubMed2024

A genetic atlas for the butterflies of continental Canada and United States.

Multi-locus genetic data for phylogeographic studies is generally limited in geographic and taxonomic scope as most studies only examine a few related species. The strong adoption of DNA barcoding has generated large datasets of mtDNA COI sequences. This work examines the butterfly fauna of Canada and United States based on 13,236 COI barcode records derived from 619 species. It compiles i) geographic maps depicting the spatial distribution of haplotypes, ii) haplotype networks (minimum spanning trees), and iii) standard indices of genetic diversity such as nucleotide diversity (π), haplotype richness (H), and a measure of spatial genetic structure (GST). High intraspecific genetic diversity and marked spatial structure were observed in the northwestern and southern North America, as well as in proximity to mountain chains. While species generally displayed concordance between genetic diversity and spatial structure, some revealed incongruence between these two metrics. Interestingly, most species falling in this category shared their barcode sequences with one at least other species. Aside from revealing large-scale phylogeographic patterns and shedding light on the processes underlying these patterns, this work also exposed cases of potential synonymy and hybridization.

Animals↗

PubMed2026

Paralog-aware assembly and filtering strategies reveal minimal nucleotide variation on the macro germline-restricted chromosome of the zebra finch.

The germline-restricted chromosome (GRC) of passerines is a remarkable tissue-specific chromosome that accumulated paralogs of genes from the regular "A chromosomes" over millions of years, often amplified into dozens of gene copies. In addition to its repetitive content, typically uniparental inheritance, and lack of recombination, the GRC resembles non-recombining sex chromosomes and some B chromosomes, for all of which assembly and single-nucleotide polymorphisms (SNPs) calling are difficult. Here, we first show that much of the Australian zebra finch macro-GRC can be assembled using accurate long reads. We then describe a paralog-aware Snakemake pipeline, ParaVar, to map short reads from the GRC to retrieve GRC regions suitable for haplotype-based analysis. ParaVar reliably calls hundreds of SNPs across the GRC, thereby providing an estimate of nucleotide diversity on the highly repetitive zebra finch macro-GRC. Our results show significantly lower nucleotide diversity (20- to 50-fold lower) on the GRC compared to the mitogenome and autosomes, and a strong phylogenetic discordance between the GRC and the mitochondrial genome. Beyond the contribution of background selection, our results suggest that a single GRC haplotype recently spread through the populations while jumping across matrilines via occasional paternal inheritance. We anticipate that our paralog-aware pipeline will be useful for SNP calling and population genetics analyses of repetitive GRCs, sex chromosomes, and B chromosomes.

Animals↗

PubMed2026

No receptor-binding domain adaptation detected in within-host H5N1 surveillance of 4,559 US dairy outbreak sequences.

BACKGROUND: The 2024-2026 US H5N1 clade 2.3.4.4b dairy cattle outbreak has been characterised primarily through consensus-level phylogenetics. Whether mammalian-adaptation variants are emerging at sub-consensus frequencies within infected hosts, particularly at the haemagglutinin receptor-binding domain (RBD), remains unknown because no systematic within-host variant analysis of the public sequencing corpus has been performed. METHODS: We conducted a pre-registered, corpus-wide intrahost single-nucleotide variant (iSNV) analysis of all publicly available H5N1 cattle, feline-spillover, and retail-milk sequences on the NCBI Sequence Read Archive (4559 samples across 7 BioProjects). A dual-caller concordance pipeline (iVar + LoFreq) with empirically determined allele frequency (AF) threshold (3%, set via four-criterion validation including synthetic spike-in controls) was applied to an 11-site Tier 1 mammalian-adaptation panel spanning the polymerase complex, haemagglutinin RBD, and accessory proteins. Within-host nucleotide diversity was compared across host categories. RESULTS: The HA RBD sites Q226L and G228S (H3 numbering) showed zero detections across >4300 adequately sequenced samples at all AF thresholds tested (1-5%), despite the pipeline detecting other non-synonymous variants at these exact codon positions (upper 95% CI for prevalence: 0.08%). Seven of eleven adaptation sites carried statistically significant iSNV signals after Bonferroni correction (corrected α = 0.00417), though all at low prevalence (≤2.95%). Genotype stratification showed that most polymerase-site detections reflected genotype structure rather than within-host emergence: the apparent PB2 631 L→M "reversion" was largely the ancestral avian state of the D1.1 genotype (20 of 23 detections), which never acquired the 631L mammalian adaptation, with only two genuine sub-consensus events in the B3.13 background, while consensus-level PB2 701N was a fixed feature of the D1.1 genotype (10 of 14 detections) rather than independent sub-consensus emergence. Cattle exhibited significantly higher within-host nucleotide diversity than feline-spillover samples (π = 1.59 × 10-4 vs 6.11 × 10-5; Kruskal-Wallis p = 6.6 × 10-15), a finding that persisted after depth-matching (p = 4.6 × 10-5); this may reflect prolonged mammary-gland infection, though sampling differences and host biology cannot be excluded. CONCLUSIONS: We did not detect HA receptor-switching adaptation (the acquisition of human-type α2,6 receptor binding via Q226L/G228S) at any tested allele frequency in the US dairy H5N1 outbreak. Sub-consensus mammalian-adaptation signals exist at polymerase-complex sites but at low prevalence, are genotype-structured rather than independently recurrent, and require functional characterisation before informing risk assessment.

Dairy cattle↗

PubMed2026

Comparative Analysis of Chloroplast Genomes Reveals Molecular Evolution and Phylogenetic Relationships in Fraxinus (Fraxinus mandshurica).

Fraxinus mandshurica (Manchurian ash) is an ecologically and economically valuable hardwood tree native to Northeast Asia, yet its genomic resources remain limited. We assembled its complete chloroplast (cp) genome (155,559 bp) using hybrid PacBio and Illumina sequencing and performed comparative, phylogenetic, and evolutionary analyses. The cp genome exhibits a typical quadripartite structure encoding 132 gene copies, comprising 114 unique genes (80 protein-coding, 30 tRNA, and 4 rRNA genes), with 18 genes duplicated in the inverted repeat (IR) regions. Simple sequence repeat analysis revealed dominance of mononucleotide A/T repeats. Phylogenetic analysis of 53 complete cp genomes strongly supported the monophyly of Oleaceae and resolved F. mandshurica as sister to the North American F. nigra, consistent with previously proposed Miocene intercontinental dispersal scenarios between East Asia and North America. Most protein-coding genes were under strong purifying selection (Ka/Ks << 1), whereas petB, rpl2, and several ndh genes showed elevated Ka/Ks values that are suggestive of altered selective constraint but are based on very few substitutions and are therefore not, on their own, evidence of positive selection. Nucleotide diversity (Pi) analysis identified 15 hypervariable intergenic spacers (mean Pi = 0.067), among which trnM-CAU-rps14, ndhJ-ndhK, and petL-petG represent promising candidate barcode regions requiring further validation. This study provides a high-quality, fully annotated cp genome of F. mandshurica and a valuable genomic resource for future phylogenetic, population genetic, and conservation studies of this important genus.

Fraxinus↗

PubMed2026

Low-pass whole-genome sequencing reveals genomic diversity and ecotype-specific adaptation in indigenous Tigrayan chickens.

Indigenous chickens play a critical role in food security and climate resilience in smallholder systems, yet their genomic diversity and adaptive potential remain insufficiently characterised. This study employed low-pass whole-genome sequencing (LP-WGS; 0.2-1.99×) to investigate genomic diversity, population structure, inbreeding and candidate environment-associated genomic variation in 33 chickens from highland, midland, and lowland agroecologies in the Tigray region of northern Ethiopia. After imputation and stringent filtering, 23.4 million high-confidence SNPs were retained, including ~ 17% novel variants, indicating substantial uncharacterised genetic diversity in these populations. SNP density (13.8 ± 8.6 SNPs/kb) was comparable to values reported from high-coverage Ethiopian chicken datasets, demonstrating the suitability of LP-WGS for population genomics in resource-limited settings. Marked differences in genomic diversity were observed among ecotypes: midland chickens showed the highest nucleotide diversity (π = 0.00267), followed by lowland (π = 0.00233), whereas highland chickens showed the lowest diversity (π = 0.00203) and elevated genomic inbreeding (FROH and FHOM ≈ 0.18). Population structure analyses revealed clear genetic separation among ecotypes. PCA (13.91% variation explained) distinguished lowland chickens along PC1 and separated highland from midland along PC2, while ADMIXTURE and FST patterns supported three major ancestral genomic backgrounds. Functional annotation of private missense variants uncovered distinct adaptive signatures reflecting the contrasting agroecological conditions. Highland chickens showed enrichment of candidate genes potentially involved in physiological processes relevant to high-altitude environments, including cold response, angiogenesis, cardiovascular regulation and metabolic homeostasis (eg., PARP1, ACOX2, ITGB3, EDNRB, SOX8, and SOX10). Midland chickens exhibited candidate signals of selection in genes with known roles in innate antiviral immunity, bacterial defence and inflammatory regulation (eg., BAK1, CLSTN1, CYSLTR1, CYSLTR2, CXCR7, GIPR, DSCAM, GDAP1, TLR3, TLR4, TLR7, IFIH1, ADORA1, EPHB1, and TMPRSS2). Lowland chickens displayed candidate variants associated with heat-stress response, DNA damage repair, oxidative balance and cardiovascular support under extreme temperatures (e.g., MLH1, BDKRB1, GPR19, FLT1, CCL18, TGM2, and RAMP3). Overall, the results indicate substantial genomic differentiation among ecotypes and suggest candidate environment-associated genetic divergence across Tigray's diverse agroecological zones. These populations may represent important reservoirs of adaptive genetic variation for climate-resilient poultry breeding, warranting further functional validation and conservation-oriented management.

Animals↗

PubMed2026

Intra-colony divergence and global allele sharing reflect purifying selection and recombination at the Botryllus histocompatibility factor locus.

Urochordates, the closest relatives of vertebrates, lack adaptive immunity. However, some taxa, such as the colonial species Botryllus schlosseri, provide a unique model for studying innate self/non-self recognition through natural allogeneic transplantation responses. In this species, interactions between colonies are controlled by a highly polymorphic locus, with the Botryllus histocompatibility factor (BHF) being the only gene known to predict tissue fusion or rejection outcomes with complete accuracy. Here, we analyzed full-length BHF alleles from 19 laboratory-born and wild colonies and found that highly divergent alleles tend to coexist within individuals, whereas identical alleles can be shared across continental-scale distances. Despite extensive length variation, evidence of intragenic recombination, and pronounced nucleotide diversity, BHF exhibits limited protein divergence, with 33 alleles encoding only 17 distinct protein variants. Contrary to expectations for polymorphic recognition genes, no evidence of balancing or directional selection was detected. Instead, signatures of purifying selection were observed. We propose that this contrast between nucleotide and protein diversity arises from the combined effects of recombination, human-mediated gene flow, and linkage to nearby loci under balancing selection, while functional constraints maintain protein stability. These findings suggest that extensive protein diversification may not be a universal driver of allorecognition gene evolution.

Animals↗

PubMed2025

Whole-genome sequencing identifies genetic diversity and adaptive signatures of hypoxia and ultraviolet radiation in Chinese chickens.

INTRODUCTION: Domestic chickens primarily descended from the wild red junglefowl, play a crucial role in global egg and meat production. China hosts diverse indigenous chicken populations that have adapted to various environmental conditions, including high-altitude with hypoxic and ultraviolet radiation stress. METHOD: We analyzed whole-genome sequences of 118 birds from five Indigenous Chinese chicken populations and 295 chicken genomes from publicly available databases to identify genomic diversity, admixture, and selection signatures of chickens adapted to high-altitude environments. Selection signatures were identified using nucleotide diversity (π), Tajima's D, XPEHH, and XP-CLR, selection scan methods. RESULTS: We observed a reduction in genetic diversity and historical declines in effective population size in high-altitude chicken, suggesting ongoing selection pressures shaping these populations. Selection scans identified nine genomic regions under strong positive selection, enriched for genes associated with hypoxia and ultraviolet radiation. Notably, five genes (TPK1, BAZ2B, MARCHF7, LLGL2, and RCAN3) were repeatedly detected across multiple selection signature analyses. RNA-seq analysis further confirmed the differential expression of these genes in the lung and heart tissues of chickens adapted to high and low altitudes, reinforcing their role in physiological adaptation to hypoxic environments. Altitude adaptation is driven by the selection of genes involved in oxygen metabolism, cellular stress response, and energy regulation. CONCLUSION: Our study provides compelling genetic evidence for differentiation between high and low and high-altitude Chinese chicken populations. These findings also ensure our understanding of local adaptation in poultry and establish a genomic framework for breeding strategies to improve environmental resilience to altitude-related stressors.

Animals↗

PubMed2026

Sexual selection purges mutation load, but not overall genetic diversity, decreasing vulnerability to extinction.

Theory suggests sexual selection will enhance population viability by purging deleterious alleles. However, direct genomic evidence for this fundamental idea is scarce and contradictory. We combined long-term experimental evolution with whole-genome resequencing to directly test how sexual selection affects mutation load, genomic divergence, and extinction risk in small populations (maximum Ne = 40) of Tribolium castaneum. After 156 generations, populations evolving under strong sexual selection carried substantially fewer deleterious alleles than populations under weak sexual selection, based on both individual-level estimates of missense and nonsense variants and population-level Rxy analyses, indicating more efficient purging of deleterious alleles. In contrast, nucleotide diversity and runs of homozygosity were similar across treatments, indicating that purging acted most strongly on deleterious variation, and that reduced mutation load in these small populations under strong sexual selection was not explained by demographic effects. Importantly, population-level mutation load estimates best explained extinction risk under inbreeding, directly linking sexual selection to purging and population viability. Genome scans of high and low sexual selection populations revealed peaks of divergence, which included genes involved in courtship, sex discrimination, and seminal fluid proteins. Our results provide direct genomic evidence that sexual selection can reduce mutation load without eroding standing genetic diversity and thus adaptive potential, while driving adaptive divergence in reproductive traits. This beneficial purging may help explain the widespread prevalence of sexual reproduction in nature despite inherent costs and have important ramifications as to how we manage populations of conservation concern.

Animals↗

PubMed2025

Combined Evidence Reveals the Origin of a Rapid Range Expansion Despite Retained Genetic Diversity and a Weak Founder Effect.

Many species are currently experiencing range shifts in response to changing environmental conditions with potentially serious genetic consequences. Repeated founder events and strong genetic drift are expected to erode genetic variation at the range front, reducing adaptive potential and slowing or even halting the expansion. However, the severity of these consequences for common and highly mobile species undergoing environment-driven range shifts (c.f. invasions) is less clear. Here, we combined historical observations and contemporary movement data of the common reed warbler (Acrocephalus scirpaceus) with genomic evidence from across its European breeding range to (1) infer the origin and (2) quantify the genetic consequences of a recent and rapid northward range expansion. Although there were no reductions in levels of nucleotide diversity or allelic richness, nor a signal of founder effect in the directionality index (ψ), our combined dataset approach was able to infer an expansion origin from the southwest. Furthermore, we found that private allelic richness retained a slight but significant linear decline along the colonisation route. These results suggest that high dispersal capabilities can allow even philopatric species to avoid the loss of genetic diversity during rapid range expansions. Nevertheless, if multiple lines of evidence enable identification of an expansion pathway, we may still detect genetic signals of expansion.

Founder Effect↗

PubMed2026

Genome-wide variation analysis of two Salvia hispanica L. genotypes and implication for associations with metabolic and adaptive traits.

BACKGROUND: Advances in next-generation sequencing have accelerated genome-wide exploration of genetic diversity in underutilized oilseed crops. Salvia hispanica L. (chia), a high-nutrient pseudocereal rich in omega-3 fatty acids, is increasingly valued for its health benefits and commercial potential, yet it remains poorly characterized at the genomic level. Understanding the scale and nature of genomic variation is essential for improving complex traits such as oil yield, stress tolerance, and seed quality. METHODS: Two contrasting chia genotypes, Black-chia (CACH-B) and White- chia (CACH-W), were resequenced using the Bio-Resequencing Toolkit (BRT) pipeline. High-coverage sequencing, with a mapping rate exceeding 99% and an average depth of approximately 28×, facilitated the detection and annotation of single-nucleotide polymorphisms (SNPs), insertions and deletions (InDels), copy-number variations (CNVs), and structural variants (SVs). The functional classification of variant impacts enabled the identification of genes potentially linked to metabolic and adaptive traits. RESULTS: A total of 1.97 million SNPs, 401,493 InDels, 836 CNVs, and 15,288 SVs were identified across the chia genome. Notably, approximately 53% of exonic SNPs were non-synonymous (dN/dS ≈ 1.28), predominantly affecting lipid metabolism, transcriptional regulation, and stress response pathways, potentially altering key agronomic traits. In addition, CNV hotspots were concentrated in chromosomes 3 and 6, overlapping MYB, WRKY, and bZIP transcription factor loci, may potentially be involved in stress tolerance and yield. Furthermore, structural rearrangements, including inversions and duplications within the FAD2, FAD3, and CYP450 gene clusters, were potentially associated with seed pigmentation and omega-3 biosynthesis, pointing to their potential breeding relevance. Observed heterozygosity (Hₒ ≈ 0.71) and nucleotide diversity (π ≈ 7 × 10-3) indicated moderate to high allelic richness. In addition, the low FST value (0.038) indicates substantial genomic similarity between the two genotypes. CONCLUSION: This study presents the first comprehensive map integrating SNPs, CNVs, and SVs in S. hispanica L. The results reveal a structurally dynamic genome characterized by substantial sequence and structural variation, providing valuable insights into genomic diversity and potential adaptive mechanisms in chia. The coexistence of high SNP diversity and abundant structural variation underpins chia's nutritional specialization and environmental resilience. These results deliver a foundational genomic resource for marker-assisted breeding, genome-wide association studies, and the development of climate-resilient chia cultivars.

Copy-number variation, structural variation↗

PubMed2025

Host-Associated Genetic Differentiation in the Face of Ongoing Gene Flow: Ecological Speciation in a Pathogenic Parasite of Freshwater Fish.

Adaptive evolution in response to varying environments, leading to population divergence, is among the most intriguing processes of speciation. However, the extent to which these adaptive processes effectively drive population divergence amidst ongoing gene flow remains controversial. Our study addresses this by analyzing population genetic structure, gene flow, and genomic divergence between lineages of a tapeworm parasite (Ligula intestinalis) isolated from sympatric fish hosts. This parasite, which must overcome host immunological defenses for successful infection, significantly impacts host health. Utilizing genome-wide Single Nucleotide Polymorphisms (SNPs) and transcriptome data, we investigated whether host species impose distinct selection pressures on parasite populations. Genetic clustering analyses revealed clear divergence, with parasites from bream (Abramis brama) forming a distinct genetic cluster separate from those infecting roach (Rutilus rutilus), rudd (Scardinius erythrophthalmus), and bleak (Alburnus alburnus). Demographic modeling indicated isolation with continuous gene flow as the most plausible scenario for this divergence. Selection analyses identified 896 SNPs under selection, displaying low to moderate nucleotide diversity and genetic divergence compared with neutral loci. Transcriptome profiling supported these findings, revealing distinct gene expression profiles between parasite populations. Examination of selected SNPs and differentially expressed genes identified candidate genes linked to immune evasion mechanisms, potentially driving ecological speciation. This research highlights the interplay of host specificity, population demography, and disruptive selection in ecological speciation. By dissecting genomic factors, our study improves the understanding of mechanisms facilitating population divergence despite ongoing gene flow.

Animals↗

PubMed2026

Genome-wide scans reveal candidate genes associated with wing morph differentiation in Tetrix japonica.

Wing dimorphism is an important dispersal-related trait in insects, but its genomic basis remains poorly understood in pygmy grasshoppers. Here, we integrated genome-wide single-nucleotide polymorphism (SNP) analyses, population structure inference, selection scans, and functional annotation to investigate genomic differentiation between long- and short-winged Tetrix japonica. Principal component analysis (PCA), ADMIXTURE, and phylogenetic analyses revealed weak genome-wide separation between morphs, indicating differentiation on a largely shared genetic background. Genome-wide scans based on the fixation index (FST), nucleotide diversity ratios, and Tajima's D, using 50-kb non-overlapping windows and empirical top-5% outlier thresholds, identified multiple candidate regions across seven chromosomes. The broader long- and short-winged candidate sets spanned 9.35 Mb and 9.37 Mb and directly overlapped 82 and 77 genes, respectively. Candidate genes were associated with signaling/hormone regulation, membrane transport, metabolism, cytoskeletal organization, extracellular matrix structure, and development. Short-winged candidate genes were significantly enriched for ABC-type transporter activity and ATP hydrolysis activity. Because all individuals originated from a single laboratory-maintained population with weak genome-wide structure, these regions should be regarded as candidate loci from a screening-stage analysis that require validation in independent populations and by functional assays, rather than as confirmed targets of selection.

Animals↗

PubMed2024

Genome-wide scan for selection signatures in Mexican Sardo Negro Zebu cattle.

The Sardo Negro cattle (SN) is the only zebu cattle breed developed in Mexico. Since its development, the selection could have led to an increase in the homozygosity level in some regions of the genome and made differentiation with other cattle populations. We aimed to identify and characterize selection signatures in SN using medium-density SNP data using four approaches: 1) Runs of homozygosity (ROH) 2) Nucleotide Diversity 3) Tajima's D and 4) the Wright's fixation index (FST). A sample of 555 SN animals genotyped for 65k SNPs was used to obtain ROH segments considered regions under selection. The FST values were estimated by comparing the sample of genotyped SN animals with samples of genotyped animals from the Gir, Brahman, and Ongole breeds. Only one region mapped to 35.78-42.51 Mb on BTA6 was considered a selection signature by the ROH method. This selection signature overlapped with the lowest diversity, negative values of Tajima's D and a diversification region between SN and the other Zebu breeds by FST. We found several candidate genes (LCORL, NCAPG, and SLIT2) related to growth and other economically important productive traits in this common region. Using the FST method, different regions, such as regions on BTA8 (8:93.4-93.9 Mb), BTA11 (11:99.2-99.7), and BTA14 (14: 26.1-26.8) related to growth and milk traits also were defined as candidate selection signatures. The selective signals identified in this study reflected the direction of the selection pressure that primarily involves the increase of live weight traits in the Sardo Negro cattle breeding program.

Animals↗

PubMed2026

Integrative Genomic and Transcriptomic Insights into High-Altitude Adaptation in Changthangi Goats.

The Changthangi goat, native to the high-altitude Ladakh Plateau in northern India, thrives in oxygen-deficient environments above 4,000 m. This study investigated the genetic basis of high-altitude adaptation in Changthangi goats by integrating comparative genomics and transcriptomics, using the tropical lowland Jamunapari goat as a comparative model. Whole-genome sequence data from 15 individuals per breed were analyzed using complementary selection sweep metrics, including nucleotide diversity, Tajima's D, iHS, CLR, XP-EHH, and FST. These analyses identified candidate genomic regions under strong selective pressure, encompassing genes involved in hypoxia sensing (HIF-1α, HIF-2α/EPAS1, EGLN1), angiogenesis (VEGFA, AGGF1, ZEB1), cardiovascular regulation (PRKCB, ESR1, RYR2), mitochondrial and energy metabolism (ACADSB, ACSS3, ACSL1), cellular stress tolerance (BCL2, ATM), and thermogenesis (UCP1, FGF21). Unlike previous caprine studies that primarily infer hypoxia adaptation from genomic signals alone, our study integrates cardiac transcriptomics to demonstrate that genomic selection in Changthangi goats is accompanied by coordinated transcriptional remodeling across interconnected physiological systems in a physiologically relevant tissue. Comparative cardiac transcriptomic profiling revealed concordant expression divergence in genes associated with oxygen transport, vascular remodeling, mitochondrial function, substrate utilization, redox balance, and genome maintenance. This integrative multi-omics framework provides a mechanistic view of caprine high-altitude adaptation and highlights the value of combining genomic selection analyses with tissue-specific transcriptional profiling to resolve complex adaptive traits.

Animals↗

PubMed2025

Genomic analysis of differentiation and demography of the formerly conspecific agile (Dipodomys agilis) and Dulzura (D. simulans) kangaroo rats.

Karyotype variation within Pacific kangaroo rat Dipodomys agilis motivated its division in 1997 into the agile kangaroo rat (AKR, D. agilis, 2N = 62) in the north of its range in California, and Dulzura kangaroo rat (DKR, D. simulans, 2N = 60) to the south, with a suspected sympatric zone south of the San Gabriel and San Bernardino Mountains. This division was supported by our whole genome sequencing that sampled a ~120 km transect from north of the mountains to SW Riverside County. The taxa showed marked genetic differentiation, with no evidence of hybridization or sympatry. AKR was found at the southern edge of the mountains, precluding the mountain barrier driving isolation, suggesting ecological separation linked to habitat differences between the mountains and the arid area to the south. Adding four additional Dipodomys species, we estimated genetic divergence times in the genus back to ∼3.5 mya. AKR and DKR diverged from D. stephensi ∼1.7 mya, and from each other ∼0.5 mya, when their joint effective population size (Ne) was ~100,000. After separation, DKR's Ne declined to ~20,000, while AKR's was little changed. More recently their Ne converged at ~50,000. Runs of homozygosity were longer in AKR, indicating a smaller neighborhood size, which may have promoted the karyotype change; however, nucleotide diversity was higher in AKR, but both had levels typical for rodents, indicating neither experienced recent bottlenecks. These patterns provide a baseline for any future conservation efforts. More generally, this study shows how a detailed genomic study can resolve taxonomic and demographic questions among morphologically indistinguishable taxa.

Animals↗