Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Single nucleotide polymorphisms (SNPs)”

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.

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

At least 19 recordsLinked to original sources
PubMed2024

Mendelian randomization study of lipid metabolism characteristics and migraine risk.

BACKGROUND: The association between serum lipids and migraine is controversial. However, randomized controlled trials have suggested that statins may be efficacious for the prevention of migraine. In this study, we aim to investigate the relationship between lipids metabolism and migraine risk. METHODS: Single-nucleotide polymorphisms (SNPs), relating to the serum lipid traits and the effect of lipid-lowering drugs that target APOB, CETP, HMGCR, NPC1L1, and PCSK9, were extracted from genome-wide association studies (GWAS) summary data. The GWAS summary data were obtained from the Global Lipids Genetic Consortium (GLGC), the UK Biobank, and the FinnGen study, respectively. Mendelian randomization (MR) analysis was performed to evaluate the association between serum lipid traits and lipid-lowering drugs with migraine risk. RESULTS: Regarding serum lipids, it was found that SNPs related to high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), non-high-density lipoprotein cholesterol (non-HDL-C), total cholesterol (TC), or triglycerides (TG) levels were not associated with migraine, migraine with aura (MA) or migraine without aura (MO). In addition, genotypes of HMGCR related to higher LDL-C levels were associated with increased risk of migraine (OR = 1.46, p = 0.035) and MA (OR = 2.03, p = 0.008); However, genotypes of PCSK9 related to higher LDL-C levels were associated with decreased risk of migraine (OR = 0.75, p = 0.001) and MA (OR = 0.69, p = 0.004); And genotypes of APOB related to higher LDL-C levels were associated with decreased risk of MO (OR = 0.62, p = 0.000). CONCLUSIONS: There is a relationship between lipid metabolism characteristics and migraine risk. SIGNIFICANCE: Based on the genome-wide association summary data, single-nucleotide polymorphisms (SNPs) related to high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), non-high-density lipoprotein cholesterol (non-HDL-C), total cholesterol (TC), or triglycerides (TG) level were not associated with risk of migraine, migraine with aura (MA) or migraine without aura (MO). However, genotypes of HMGCR related to higher LDL-C levels have shown an increased risk on migraine and MA. And genotypes of APOB or PCSK9 related to higher LDL-C levels have shown a decreased risk on MO, or migraine and MA, respectively. These results suggested that there may be a relationship between lipid metabolism characteristics and the risk for migraine development.

Humans↗
PubMed2026

Upper airway microbiome interacts with GSDMB and ORMDL3 asthma risk SNPs to influence early-life wheeze risk.

BACKGROUND: Single-nucleotide polymorphisms (SNPs) in the chromosome 17q12-q21 region and, independently, early-life nasal microbiota dominated by Moraxella, Streptococcus, or Haemophilus (MSH) increase risk of chronic wheeze and asthma development. OBJECTIVE: We sought to determine whether 17q12-q21 risk SNPs and nasal microbiota interact to modulate childhood wheeze risk. METHODS: Nasal wash samples from 12-month-old infants in 2 birth cohorts, COAST (Childhood Origins of Asthma; n = 180) and URECA (Urban Environment and Childhood Asthma; n = 139), underwent 16S ribosomal RNA variable region 4 sequencing. Nasal microbiota dominated by MSH or Corynebacterium, Dolosigranulum, Staphylococcus, or Bacillus (CDSB) were assessed. Paired blood was genotyped for 9 17q12-q21 risk SNPs. Logistic regression tested interactions between 17q12-q21 SNPs and MSH or CDSB on wheeze risk in the first 3 years of life. A549 lung epithelial cells, CRISPR-edited to encode the rs7216389 risk genotype (rs7216389TT) were compared to the heterozygous (rs7216389CT) line using bulk RNA sequencing. RESULTS: SNPs, particularly those in the ORMDL3 (rs8076131; odds ratio [OR]: 1.72; 95% CI: 1.09-2.71; Pint = .031) and GSDMB (rs2305480; OR: 1.72; 95% CI: 1.09-2.71; Pint = 0.042; and rs7216389; OR: 1.73; 95% CI: 1.09-2.70; Pint = .047) genes, interact with MSH microbiota to increase early-life wheeze risk (false discovery rate Pint = .016 for all), while interactions with CDSB reduce risk. A549 airway epithelial cells homozygous for rs7216389TT exhibited decreased expression of genes involved in antimicrobial responses and neutrophil recruitment and evidence increased microbial adherence compared with the heterozygous cell line. CONCLUSION: Airway microbiota interact with SNPs at the 17q12-q21 locus in genes involved in sphingolipid metabolism and intracellular antimicrobial responses, to modulate wheeze risk.

Humans↗
PubMed2026

Correlations between causal effect sizes of proximal SNPs vary with functional annotations and implicate stabilizing selection.

Causal disease effect sizes of proximal single-nucleotide polymorphisms (SNPs) are widely assumed to be independent but could be correlated. Here we introduce a new method, linkage disequilibrium SNP-pair effect correlation regression (LDSPEC), to estimate the correlation of causal disease effect sizes of derived alleles between proximal SNPs; LDSPEC produced robust estimates in simulations. Analyzing 70 UK Biobank diseases and traits (average N = 305,646), we detected significantly non-zero SNP-pair effect correlations (for example, -0.37 ± 0.09 for low-frequency positive linkage disequilibrium 0-100-bp SNP pairs) that decayed with distance and varied with allele frequency and linkage disequilibrium between SNPs. SNP pairs with shared functions had stronger effect correlations that spanned longer genomic distances. Consequently, SNP heritability estimates were smaller than estimates of the sum of causal effect size variances across SNPs, particularly for certain functional annotations. We recapitulated our findings via forward simulations involving stabilizing selection, implicating the action of linkage masking, whereby haplotypes containing linked SNPs with opposite effects on disease have reduced effects on fitness and escape negative selection.

Polymorphism, Single Nucleotide↗
PubMed2026

Emerging terbinafine-resistant Trichophyton indotineae between 2018 and 2023: a multinational genomic epidemiology study.

BACKGROUND: Trichophyton species cause the greatest burden of dermatophytosis worldwide, with the Trichophyton mentagrophytes species complex being particularly associated with the emergence of new aggressive infections. One emerging species, Trichophyton indotineae is notable for its clinical resistance to terbinafine antifungal treatment and rapid global spread. In this study we aim to characterise the epidemiology of this emerging pathogen using genomics. METHODS: In this retrospective genomic epidemiology study, to better understand the epidemiology of this disease, we sourced isolates collected from patients with severe cases of dermatophytosis (identified either by internal transcribed spacer sequencing or phenotypic characterisation) in the UK, Ireland, France, Canada, and India for the period 2014-23, including the T indotineae type strain from Japan. We used whole-genome sequencing to confirm 90 isolates were T indotineae, and antifungal susceptibility testing to assess susceptibility to terbinafine. FINDINGS: 103 cases of severe dermatophytosis caused by Trichophyton species collected between 2018 and 2023 in the UK, France, Canada, Ireland, and India were included in this study. Susceptibility testing indicated that 63 (70%) of 90 T indotineae isolates were resistant to terbinafine (minimum inhibitory concentration [MIC] ≥0·5 mg/L). Pairwise genetic distances showed very high identity with only 147 (range 1-414) single-nucleotide polymorphisms (SNPs) separating isolates that were nested within a monophyletic phylogeny, supporting a single evolutionary origin of T indotineae. Genome-wide analyses identified multiple non-synonymous SNPs in SQLE (ERG1), the squalene epoxidase target of terbinafine, that were associated with terbinafine in vitro resistance of 0·5 mg/L or higher. However, six isolates exhibited high MIC values without SQLE mutations, suggesting the presence of alternative resistance mechanisms. INTERPRETATION: That no clear geographical clustering of isolates was observed confirms the rapid transcontinental spread of T indotineae from its likely centre of diversity in Asia. Our findings highlight the importance of better genomic surveillance to understand and manage this severe and rapidly emerging terbinafine-resistant dermatophyte. FUNDING: None.

Terbinafine↗
PubMed2026

Adaptation to Plant Defence in an Agricultural Insect Pest: Integrating Genome Scans and Gene Expression in the Soybean Aphid Reveals Multi-Genic Pathways.

In agroecosystems, intense selection pressures cause species to adapt and spread, often leading to the evolution and persistence of pests. Understanding how pests rapidly adapt can help develop sustainable strategies for their management and improve agroecosystem health. Pest adaptation involves stable variations in DNA sequence, as well as dynamic shifts in gene expression, often mediated by non-coding regulatory elements. We examined adaptation to plant defences in the soybean aphid, Aphis glycines, in which virulent aphids have overcome plant defences and avirulent aphids have not. Previous data with laboratory colonies suggested that virulent aphids have higher overall gene expression, including transposable elements, some of which influence gene regulation. However, we lack information on how genetic variation in natural populations impacts adaptation and potentially gene regulation. We integrated population genome scans of field-collected, soybean aphid populations with gene expression profiles of virulent and avirulent laboratory colonies to uncover connections between genetic differentiation and gene regulation for virulence. Genome scan methods found 2144 single nucleotide polymorphisms (SNPs) with significant genetic differentiation (i.e., outliers) in field-collected populations. These SNPs were near 1004 genes, representing 5.16% of the effective number of genes. Based on previous RNA-Seq data with laboratory colonies, we found 3160 genes and 147 long non-coding RNAs (lncRNAs) with differential expression among virulent and avirulent biotypes. By integrating both data sets, we identified 16 genes and 5 long non-coding RNAs with differential expression and that were associated with an outlier SNP (within 10 kbp). We validated SNPs with additional field collected aphids and found an aphid clone with stronger virulence than our laboratory virulent colony, surviving on 2 different aphid-resistant soybean varieties. This new virulent clone had fixed allele differences at 9 SNPs compared to our avirulent and other virulent colony. Field collected soybean aphids matching the phenotype of this new virulent clone had significant genetic differentiation with 3 outlier SNPs near genes related to zinc transport and lachesin compared to field collected avirulent aphids. Our entire data reinforced the importance of a potential multi-genetic response to overcome plant defence and generates new insights into complex genetic and regulatory mechanisms involved in insect-plant interactions.

Animals↗
PubMed2025

Genetic Determinants of Leisure-Time Physical Activity in the Taiwanese Population: A Genome-Wide Association Study.

BACKGROUND: Physical inactivity contributes to systemic disease burden and premature mortality worldwide. Leisure-time physical activity (LTPA) improves health outcomes; however, its genetic determinants, particularly in Asian populations, remain unclear. This study aimed to identify genetic loci associated with LTPA in the Taiwanese population. METHODS: We conducted genome-wide association studies in 122,258 Taiwan Biobank participants. LTPA was assessed both as a binary trait (regular exerciser vs non-exerciser) and an ordinal trait (categorized by MET-hours per week into low, moderate, and high physical activity levels). Logistic and ordinal logistic regression models were used under an additive genetic model, adjusting for age, age 2 , sex, body mass index, smoking, and the first 10 genetic principal components. Candidate nonsynonymous mutations were further examined in 1494 whole-genome sequenced participants. RESULTS: Binary trait genome-wide association studies identified genome-wide significant (GWS) loci at ATXN2 (12q24.12), FTO (16q12.2), and NOTCH4 (6p21.32), with associations for FTO and NOTCH4 only observed in body mass index (BMI)-adjusted models. Ordinal trait analysis (<10, 10-<20, &#x2265;20 MET&#xb7;h&#xb7;wk -1 ) identified a single GWS locus at BRAP (12q24.12). Fine-mapping of 12q24.12 revealed multiple GWS single-nucleotide polymorphisms (SNPs) in strong linkage disequilibrium with lead variants; these signals largely disappeared after conditional analysis, consistent with a single underlying association. Whole-genome sequencing and linkage disequilibrium analysis identified three GWS nonsynonymous mutations, with ALDH2 rs671 emerging as the most likely causal variant. CONCLUSIONS: ATXN2-ALDH2 region on chromosome 12q24.12 was identified as a key locus for LTPA in Taiwanese individuals. These findings enhance our understanding of the genetic basis of physical activity and may inform future precision medicine and public health strategies.

Adult↗
PubMed2025

Association of CYP19 gene SNPs (rs7176005 and rs6493497) with polycystic ovary syndrome susceptibility in Northern Chinese women.

PURPOSE: The objective of this study was to elucidate the relationship between two single nucleotide polymorphisms (SNPs) rs7176005 and rs6493497 in CYP19 gene and the risk of polycystic ovary syndrome (PCOS) in Northern Chinese women. METHODS: In this case-control study, a total of 340 women with PCOS and 340 matched healthy controls were recruited. Polymerase chain reaction ligase detection reaction (PCR-LDR) method was used to investigate two SNPs (rs7176005 and rs6493497) in the 5'-flanking region of CYP19 gene exon 1. RESULTS: We observed a significant association of rs7176005 and rs6493497 with reduced risk of PCOS. Compared with CC genotype, a significant association of CT genotype (p&#x2009;=&#x2009;0.019), TT genotype (p&#x2009;<&#x2009;0.001) and combined CT&#x2009;+&#x2009;TT genotype (p&#x2009;<&#x2009;0.001) with reduced risk of PCOS was observed. The result of linkage disequilibrium analysis showed that these two SNPs are in complete linkage disequilibrium (r2 = 1). For rs7176005 SNP, compared with CC genotype, CT, TT and CT&#x2009;+&#x2009;TT genotypes reduced the risk of PCOS. The age, BMI-adjusted OR were 0.650 (95% CI&#x2009;=&#x2009;0.460-0.917), 0.158 (95% CI&#x2009;=&#x2009;0.066-0.376) and 0.545(95% CI&#x2009;=&#x2009;0.391-0.759), respectively. CONCLUSIONS: These findings highlight a significant association between CYP19 gene polymorphisms and PCOS susceptibility, implying potential protective effects of T and A alleles. Of course, the major limitation of this study is the sample size of the case-control study. Larger cohort studies are needed to confirm these findings and investigate the underlying causes.

Adult↗
PubMed2026

Comparisons Between Large-Scale Genomic Variants and SNPs in Driving Population Divergence and Local Adaptation.

Genomic variations, such as indels (2-49 bp) and structural variants (SVs, &#x2265;50 bp), are larger-scale mutations than single nucleotide polymorphisms (SNPs) and can substantially impact evolutionary processes, including speciation, adaptation, and phenotypes. Despite their functional importance, integrative population genetic analyses that jointly consider genome-wide SNPs, indels, and SVs remain under-explored. The ground tit (Pseudopodoces humilis), an endemic species to the Qinghai-Tibet Plateau (QTP), exhibits divergence across distinct glacial refugia, accompanied by habitat and morphological divergence, making it an excellent example for investigating how different types of genomic variants contribute to population divergence and local adaptation. Here, by retrieving 81 whole-genome sequence data, over 13 million SNPs, 2 million indels, and 22,101 SVs were identified. Variants were unevenly distributed across the genome, characterized by distinct hotspot regions. Indels and SVs revealed four genetic clusters consistent with previous SNP-based results, thereby validating the reliability of our variant datasets. FST and genotype-environment association (GEA) analyses independently revealed numerous candidate indels and SVs; each showed minimal overlap with previously identified SNPs, and were enriched in similar functional pathways such as signal transduction, skeletal muscle development, water transport, DNA repair, reproduction, nervous system development, and immunity. Collectively, our results demonstrated that indels and SVs could capture additional signatures besides SNPs. Furthermore, similar but distinct gene functions among different types of genomic variants collectively and complementarily drive genomic divergence across environmental gradients in such a high-elevation endemic species, underscoring its evolutionary relevance in local adaptation.

indels↗
PubMed2025

Genotypes of SNPs of key genes regulate susceptibility and drug sensitivity to neovascular AMD in the human population.

OBJECTIVE: To compare the genetic characteristics of the normal control group to those of neovascular age-related macular degeneration (AMD) patients and to detect single-nucleotide polymorphisms (SNPs) related to the pathogenesis of neovascular AMD and the sensitivity to anti-VEGF drug, combercept. METHOD: This is a prospective case-controlled study. A total of 104 neovascular AMD patients were treated with combercept and 106 normal subjects were served as the control group. SNPs associated with neovascular AMD and disease susceptibility and drug sensitivity were analysed. RESULTS: Significant differences existed between neovascular AMD patients and normal subjects among genotypes of the SNPs of two genes, ARMS2 (rs10490924 T) and HTRA 1 (rs11200638 A). The T alleles in rs1065489 of CFH and the rs2230205 of C3 significantly promoted neovascular AMD in males while having no significant effect in females. Six SNPs of five genes, including C3 (rs2250656 G), CFB (rs2072633 G), CFH (rs2274700 A, rs3766405 T), KDR (rs6828477 A) and FZD 4 (rs10898563 T), had significant impact in reducing neovascular AMD. Two SNPs of the CFH gene (rs2274700 A and rs3766405 T) and one SNP of the CFB gene, rs2072633 G, were statistically significantly associated with good response to combercept. Conversely, the other two SNPs of the CFH gene, rs1065489 T and rs3753396 G, and the rs7412 T of the APOE gene were associated with a relatively poor patient response to drug action. Two sets of SNPs of CFB have a combined positive effect on disease. The two SNPs of CFH (rs1065489 T and rs3753396 G) and the combination of the two SNPs of CFH and rs7412T of APOE have negative effects on the drug effectiveness. CONCLUSIONS: These genotype differences facilitate the selection of individualised treatment options towards obtaining the most efficacious clinical treatment. These findings need to be validated by studies with different ethnic populations and/or larger samples.

Humans↗
PubMed2025

Replication and Functional Prediction of Two GWAS-Reported SNPs Located on RAD50 Gene Associated with Asthma in Pakistani Children.

BACKGROUND: Genome-wide association studies (GWAS) have indicated that several single nucleotide variants (SNVs) of the RAD50 gene are significantly associated with childhood-onset asthma. However, the biological role of RAD50, and its genomic variants that predispose individuals to asthma, remains unclear. This case-control study aimed to investigate the association of two Single nucleotide polymorphisms (SNPs) rs2244012, and rs6871536 of RAD50 with asthma susceptibility using experimental and computational tools. METHODS: The case-control study involved 355 participants: "176 asthma cases [mean age (sd) = 8.91 &#xb1;3.05] and 179 healthy controls [mean age (sd) = 11.10 &#xb1;8.86] from local Punjabi population of Pakistan. The SNPs were analyzed using a modified single base extension method. The allelic association with asthma and linkage disequilibrium (LD) between the two main SNPs were performed using the SHEsis tool. SNPStats was used to assess the association of SNPs under genotypic models and interaction with non-genetic factors. The LD calculator of ENSEMBL employed for the identification of proxy SNPs in high LD (r^2 > 0.97) to main SNPs. Additionally, HaploReg(v4.1) was utilized to gauge the impact of SNPs on genomic regulations. RESULTS: In current study, both SNPs were found to have a significant association (p-value <0.05) with childhood-onset asthma development under allelic and genotypic models. The alternative "G" allele of rs2244012 is shown to modify two regulatory motifs: Nrf-2 and Zbtb12, while the alternative "C" allele of rs6871536 is predicted to alter the OSF-2 motif. Moreover, 10 SNVs proximal to rs2244012 and 21 SNVs near rs6871536 are in high LD in the Punjabi population of Lahore, Pakistan (PJL). These proxy/high-LD SNVs also displayed the potential to change DNA regulatory motifs. CONCLUSION: the rs2244012, and rs6871536 variants of RAD50 gene are significantly association with childhood asthma in Pakistan. Despite being intronic variants, it is our inference that these two SNPs have the potential to either independently or synergistically regulate inflammatory responses via nearby SNVs.

Asthma↗
PubMed2026

A systematic review and network meta-analysis of single nucleotide polymorphisms associated with oral submucous fibrosis risk.

BACKGROUND: Oral submucous fibrosis (OSF) is a chronic and insidious oral disease characterized by hyalinization of the subepithelial connective tissue and progressive fibrosis of the oral submucosa. It is a precancerous condition of oral squamous cell carcinoma. Studies have demonstrated that single nucleotide polymorphisms (SNPs) are closely associated with susceptibility to OSF. This study aims to comprehensively evaluate the association between SNPs and OSF risk and to rank the strength of the association between different genetic models and OSF susceptibility. METHODS: Literature related to OSF was comprehensively searched from PubMed, Web of Science, Embase, Cochrane Library, CNKI, and Wangfang databases up to July 2025. Full-text case-control studies with patients diagnosed with OSF were included. Quality assessment was performed to evaluate the risk of bias. RevMan 5.4, GeMTC 0.14.3, and STATA 17.0 were used for the pairwise and Bayesian network meta-analysis. RESULTS: A total of 24 studies with 2545 cases and 3772 controls, covering 13 SNPs in 11 genes, were included in our meta-analysis. We found that CYP1A1 rs4646903:T>C, CYP1A1 rs1048943:A>G, GSTT1 null genotype, GSTM1 null genotype, and XRCC3 rs861539:C>T were associated with an increased risk of OSF, while MMP2 rs243865:C>T and MMP3 rs3025058: 5A>6A were associated with a decreased risk of OSF. Further Bayesian network meta-analysis indicated the top 5 genetic models with the highest association with OSF risk in network group 1 were the dominant model, homozygous model, allelic model, and recessive model of CYP1A1 rs1048943:A>G (ranked 1-4), and the heterozygous/dominant model of CYP1A1 rs4646903:T>C (both ranked 5). While the allelic models of XRCC3 rs861539:C>T and MMP3 rs3025058: 5A>6A ranked first for predicting OSF in group 2 and group 3, respectively. CONCLUSION: Some specific SNPs are significantly related to the risk of OSF. Among them, the dominant model of CYP1A1 rs1048943:A>G may be the most strongly associated genetic model with OSF risk. Future large-sample, well-designed studies with detailed genotype data are needed to validate the roles of these SNPs in OSF risk.

Humans↗
PubMed2025

Construction of a Core Germplasm and Identification of Candidate SNPs Associated with Growth Performance of Epinephelus tukula by Whole-Genome Resequencing.

Epinephelus tukula is an economically important aquaculture animal, and a major parent in grouper crossbreeding. To better preserve and exploit E. tukula germplasm resources, a core collection (containing 34 individuals derived from 10 genetic groups) was first constructed based on phenotypic growth traits and whole-genome resequencing (WGS) data. The phenotypic traits of the individuals within the core collection were not significantly different from those in the original collection, suggesting effective representativeness of the core collection. Additionally, we performed genome-wide association study (GWAS) of E. tukula to identify candidate single nucleotide polymorphisms (SNPs) and genes associated with growth traits, to facilitate the improvements in the growth performance of this species. Twenty-six significant SNPs were identified, scattered among multiple chromosomes. Five SNPs were confirmed to be correlated with growth in another new group of 101 individuals. Based on the annotation results, these five SNPs were located in CCDC102A, NTRK2, CTSL, OTOF, and nestin, and were involved in cell development, differentiation and proliferation, glycolytic metabolism, neurological development, and myoblast differentiation. Our findings not only provide an effective basis for the conservation and utilization of E. tukula germplasm resources, but also promote the development of marker-assisted selection of E. tukula.

Polymorphism, Single Nucleotide↗
PubMed2026

European ash pangenome reveals widespread structural variation and genetic basis of low ash dieback susceptibility.

European Ash (Fraxinus excelsior) is a keystone tree species, whose populations are being decimated by ash dieback disease (ADB) - better characterisation of genetic variants associated with low susceptibility to the disease is needed. Here, we develop a F. excelsior pangenome to more fully capture sequence variability within this species compared with a linear reference genome, using a geographically diverse set of fifty F. excelsior samples. We identify 362,965 structural variants (SVs), including 174&#x2009;Mb of sequence absent from the linear reference genome (22% of the linear reference size), and identify 3,412 high-confidence dispensable genes (those present only in some individuals). We use the pangenome to analyse existing genomic data from over 1,200 individuals, revealing 220 single nucleotide polymorphisms (SNPs) showing consistent allele frequency shifts between healthy individuals and those highly damaged by ADB, across UK seed sources, explicitly demonstrating the existence of a shared genetic component to low ADB susceptibility.

Polymorphism, Single Nucleotide↗
PubMed2026

Phylogenetic inconsistency of pairwise SNP clustering for inferring tuberculosis transmission in a high-burden, endemic setting: a case study from Thailand.

Whole-genome sequence analysis is now widely used to delineate tuberculosis transmission clusters. A standard practice is to cluster bacterial isolates based on a fixed maximum genome-wide pairwise single nucleotide polymorphism (pwSNP) distance threshold. In this study, we evaluated the phylogenetic consistency of pwSNP-distance clustering with thresholds ranging between 1 and 25 single nucleotide polymorphisms (SNPs) using two contrasting data sets: (i) a data set from the UK (N = 390) published by T. M. Walker, C. L. C. Ip, R. H. Harrell, J. T. Evans, et al. (Lancet Infect Dis 13:137-146, 2013, https://doi.org/10.1016/S1473-3099(12)70277-3), which was foundational to the establishment of this method, and (ii) a data set from Thailand (N = 3,341), characterized by persistent transmission and sparse, non-systematic sampling. For the UK data set, the standard pwSNP-distance clustering using thresholds of &#x2265;12 SNPs yielded entirely monophyletic clusters and showed high concordance with a comparative monophyly constrained, tree-based method. In contrast, for the Thai data set, pwSNP-distance clustering often generated non-monophyletic clusters, even by the 25-SNP threshold. The pwSNP-distance and comparative tree-based clustering methods only showed large consistency at thresholds of &#x2265;22 SNPs. This suggests that SNP clusters defined by low distance thresholds (i.e., <12 SNPs for the UK data set, and <22 SNPs for the Thai data set) may lack robustness, and the problem is particularly severe for data sets characterized by persistent transmission, likely due to poorer cluster separation. Moreover, our findings indicate that large cluster sizes, high maximum intra-cluster genetic distances, and broad sample collection time spans may serve as useful indicators of potentially non-monophyletic clusters. We also demonstrate that mixed infections can produce spurious, phylogenetically long-range SNP linkages, underscoring the necessity of strict sequence quality control.IMPORTANCEFixed-threshold pairwise single nucleotide polymorphism (pwSNP)-distance clustering is commonly used to delineate tuberculosis transmission clusters. From an epidemiological perspective, a genuine transmission cluster must be monophyletic, originating from a single source. However, pwSNP-distance clustering is inherently simplistic and can therefore violate this principle, making the assessment of its phylogenetic consistency critical. Our results demonstrate that while this method effectively delineated complete transmission clusters for the data set from the UK, a low-burden and non-persistent transmission setting, it frequently generated non-monophyletic clusters when applied to the Thai data set, characterized by persistent transmission alongside sparse and non-systematic sampling. Furthermore, we found that clusters derived using low distance thresholds could notably vary between the pwSNP-distance and comparative tree-based clustering methods, suggesting limited reliability and robustness. To accurately delineate tuberculosis transmission clusters, especially for complex data from high-burden, endemic settings, we recommend transitioning from pwSNP-distance clustering toward more robust, phylogenetic clustering that respects evolutionary descent.

Mycobacterium tuberculosis↗
PubMed2025

Genomic prediction and genome-wide association study for liver abscesses in crossbred beef cattle.

Liver abscesses are a concern in feedlot cattle, and little is known about the role of genetics in their development. This study aimed to estimate genetic parameters and to identify single-nucleotide polymorphisms (SNPs) associated with liver abscesses. Crossbred cattle representing 18 breeds in the U.S. Meat Animal Research Center Germplasm Evaluation Program were phenotyped for liver abscesses at slaughter (n&#x2005;=&#x2005;9,044). Seventeen percent of cattle had liver abscesses. These cattle had genotypes that were imputed to sequence variant genotypes. After filtering and quality control, 340,723 SNPs were used in the analysis. Liver abscess prevalence was modeled with a single-step genomic best linear unbiased prediction (ssGBLUP) threshold model using a Bayesian framework. The model included contemporary group (sex, treatment group, and slaughter date), additive genomic, and residual effects. Genomic heritability was 0.039 (95% highest posterior density&#x2005;=&#x2005;0.005, 0.081), which was very small. To assess prediction quality, a 5-fold random cross-validation structure was used. Method Linear Regression was used to assess accuracy, bias, and dispersion by comparing estimated breeding values (EBV) from full and reduced analyses. Cross-validation metrics showed EBV based on genotypes had 0.05 reliability (SD&#x2005;<&#x2005;0.01) with no bias relative to EBV based on genotypes and phenotypes. For the genome-wide association study, SNP effects were back calculated from the EBV solutions from ssGBLUP. No SNPs were associated with liver abscesses at a Benjamini-Hochberg adjusted 0.05 significance level. Although a large dataset was used, this result was because of the low genomic heritability and imprecise EBV used to calculate SNP effects. Based on these results, environmental factors contribute to most of the variation in liver abscesses. Genetic selection to reduce liver abscesses would be slow because of the low genomic heritability, measurement late in life, and inability to measure breeding animals. A faster approach would be finding additional environmental interventions that maintain animal performance.

Animals↗
PubMed2026

Integrative Analysis Uncover the Effects and Multi-Omics Features of Thigh Muscle Fat Infiltration.

The health impacts and underlying biological pathways of thigh muscle fat infiltration (TMFI) remain incompletely understood. In this study, we analyzed TMFI measured by magnetic resonance imaging in 55,120 UK Biobank participants and found that higher TMFI was significantly associated with all-cause mortality as well as with all major system-specific diseases examined (p values ranged from 2.50&#x2009;&#xd7;&#x2009;10-88 to 9.97&#x2009;&#xd7;&#x2009;10-04). TMFI also mediated the effects of lifestyle factors on health-related outcomes, with mediation proportions ranging from 6.7% to 71.7%. A genome-wide association study (GWAS) identified 79 lead single nucleotide polymorphisms (SNPs) linked to TMFI, and the polygenic risk score for TMFI was significantly associated with mortality and all incident diseases across examined organ systems in an independent subset of UK Biobank participants of European ancestry who were not included in the TMFI GWAS (n&#x2009;=&#x2009;362,286, all p&#x2009;<&#x2009;0.05). Gene-drug interactions identified multiple drugs that could potentially modulate TMFI. Analysis of single-cell transcriptomic data indicated that myogenic cells were strongly linked to TMFI (p&#x2009;=&#x2009;7.08&#x2009;&#xd7;&#x2009;10-08). Summary-data-based Mendelian randomization and Transcriptome-Wide Association Study analyses revealed numerous genes whose expression in specific tissues was associated with TMFI. Proteomic and metabolomic profiling uncovered a broad array of circulating biomarkers associated with TMFI, many of which mediated the effects of modifiable factors and genetic risk on TMFI. Overall, our results highlight the biological relevance of TMFI to human health and provide insights into the multi-omics mechanisms underlying TMFI, identifying potential targets for interventions.

Humans↗
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↗
PubMed2024

Estimation of SARS-CoV-2 fitness gains from genomic surveillance data without prior lineage classification.

The emergence of SARS-CoV-2 variants with increased fitness has had a strong impact on the epidemiology of COVID-19, with the higher effective reproduction number of the viral variants leading to new epidemic waves. Tracking such variants and their genetic signatures, using data collected through genomic surveillance, is therefore crucial for forecasting likely surges in incidence. Current methods of estimating fitness advantages of variants rely on tracking the changing proportion of a particular lineage over time, but describing successful lineages in a rapidly evolving viral population is a difficult task. We propose a method of estimating fitness gains directly from nucleotide information generated by genomic surveillance, without a priori assigning isolates to lineages from phylogenies, based solely on the abundance of single nucleotide polymorphisms (SNPs). The method is based on mapping changes in the genetic population structure over time. Changes in the abundance of SNPs associated with periods of increasing fitness allow for the unbiased discovery of new variants, thereby obviating a deliberate lineage assignment and phylogenetic inference. We conclude that the method provides a fast and reliable way to estimate fitness advantages of variants without the need for a priori assigning isolates to lineages.

COVID-19↗