SEARCH · PubMed Health

Results for “Bayesian estimation and selection”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

16 recordsLinked to original sources

PubMed2025

Bayesian Modeling of Cancer Outcomes Using Genetic Variables Assisted by Pathological Imaging Data.

With the increasing maturity of genetic profiling, an essential and routine task in cancer research is to model disease outcomes/phenotypes using genetic variables. Many methods have been successfully developed. However, oftentimes, empirical performance is unsatisfactory because of a "lack of information." In cancer research and clinical practice, a source of information that is broadly available and highly cost-effective comes from pathological images, which are routinely collected for definitive diagnosis and staging. In this article, we consider a Bayesian approach for selecting relevant genetic variables and modeling their relationships with a cancer outcome/phenotype. We propose borrowing information from (manually curated, low-dimensional) pathological imaging features via reinforcing the same selection results for the cancer outcome and imaging features. We further develop a weighting strategy to accommodate the scenario where information borrowing may not be equally effective for all subjects. Computation is carefully examined. Simulations demonstrate competitive performance of the proposed approach. We analyze TCGA (The Cancer Genome Atlas) LUAD (lung adenocarcinoma) data, with overall survival and gene expressions being the outcome and genetic variables, respectively. Findings different from the alternatives and with sound properties are made.

Humans↗

PubMed2026

Spatially Contextualized Integrative Genomics Highlights Neuronal and Glial Regulatory Programs in Low Back Pain.

PURPOSE: Low back pain (LBP) is a heterogeneous pain condition with a measurable genetic contribution, but the genes, brain cell types, and spatial tissue contexts through which inherited risk is expressed remain unclear. We aimed to define cell-type-specific and spatially contextualized genetic mechanisms underlying LBP. METHODS: FinnGen R12 LBP GWAS summary statistics (42,521 cases and 353,224 controls) were integrated with brain single-nuclei eQTL data across eight major brain cell classes. We evaluated genome-wide polygenic signal using LDSC, prioritized genes using MAGMA and PoPS, and performed brain cell-type-specific eQTL-anchored Mendelian randomization, primarily based on single-instrument Wald ratio estimates, followed by Bayesian colocalization. Spatial genetic mapping was conducted using gsMap in an E16.5 murine embryonic atlas and two adult human lumbar spinal cord Visium sections. Selected candidates were assessed by RT-qPCR in neuronal-like and astroglial-like inflammatory cell models. RESULTS: LDSC supported interpretable polygenic signal for LBP. MAGMA and PoPS showed partial gene-level convergence, with TCF4 and TMEFF2 supported by both approaches. Across 1641 tested gene-cell type exposures, significant eQTL-anchored MR associations were concentrated in excitatory neurons, oligodendrocytes, inhibitory neurons, and astrocytes. Integrated eQTL-anchored MR, colocalization, and gene-prioritization evidence highlighted CLEC18A, QPRT, and GMPPB as higher-priority non-MHC candidates with moderate, but not strong, colocalization support. gsMap localized LBP-associated enrichment to neuroaxis-related embryonic regions, including brain, spinal cord, sympathetic nerve, and dorsal root ganglion, and to neuronal-like niches in adult lumbar spinal cord. RT-qPCR showed model-dependent expression changes, with QPRT and LGI4 preferentially responsive in neuronal-like SH-SY5Y cells and GMPPB and DPYSL5 responsive in astroglial-like U251 cells. CONCLUSION: These findings support neuronal and glial regulatory programs as plausible contributors to LBP genetic susceptibility and highlight CLEC18A, QPRT, and GMPPB as higher-priority non-MHC candidates with moderate colocalization support. The results provide a spatially contextualized framework for candidate prioritization in LBP, while emphasizing the need for larger cell-type-specific eQTL resources and functional validation before therapeutic or mechanistic conclusions can be drawn.

Mendelian randomization↗

PubMed2026

Introgression among maternal lineages inferred from complete mitogenomes and molecular dating helps resolve phylogeography of European roe deer.

BACKGROUND: The European roe deer (Capreolus capreolus) is one of the most widespread ungulates in Europe, with a phylogeographic structure mainly shaped by Pleistocene glacial cycles and secondary contacts with the Siberian roe deer (C. pygargus). METHODS: We sequenced 52 complete mitogenomes of C. capreolus from Slovenia, Poland and France, and combined them with 24 publicly available sequences of C. capreolus and C. pygargus, yielding an alignment of 76 genomes representing 59 haplotypes (42 from C. capreolus and 17 from C. pygargus). Phylogeographic structure was assessed using a median-joining network, and divergence times were estimated using a time-calibrated Bayesian phylogeny based on mitochondrial coding regions, incorporating published ancient C. pygargus mitogenomes. We additionally screened mitochondrial protein-coding genes for selection. RESULTS: The haplotype network recovered the three major European roe deer clades (Eastern, Central, and Western) and detected Central-clade haplotypes in France. Two Polish haplotypes (Cp9 and Cp10), detected in C. capreolus, clustered within the C. pygargus mitochondrial lineage, supporting mitochondrial introgression. Time-calibrated phylogenies placed introgressed haplotypes within established C. pygargus lineages. Selection analyses provided limited evidence for episodic positive selection restricted to a small number of codons. CONCLUSIONS: Whole mitogenomes improve resolution of roe deer phylogeography and reveal introgressed maternal lineages, while time-calibrated phylogenies and selection tests add evolutionary context for interpreting mtDNA diversity in genus Capreolus.

Animals↗

PubMed2026

Genomic background of gestation length and calving-related traits in Holstein cattle.

The reproductive success of cows directly influences the profitability of dairy farms. Reproductive traits, particularly calving-related traits, generally have low heritability but sufficient additive genetic variance to enable genetic progress through genomic selection. Thus, the primary objectives of this study were to estimate genetic parameters and perform single-step genome-wide association studies (ssGWAS) for calf size, calving ease, gestation length, and stillbirth in Holstein cattle. Variance components were estimated based on animal models and Bayesian inference using a data set containing 226,717 animals with phenotypic records, 15,761 animals genotyped with 45,101 SNP markers, and 461,819 animals in the pedigree. SNP effects were estimated using the single-step GBLUP method. For direct and maternal genetic effects, heritability estimates (posterior standard deviation) ranged from 0.001 (0.002) for gestation length in heifers to 0.16 (0.001) for gestation length in cows. Genetic correlations ranged from -0.57 (0.01) between calving ease and stillbirth in heifers to 0.74 (0.01) between gestation length evaluated in heifers and cows. The ssGWAS results supported a highly polygenic architecture for calving-related traits, with most genomic signals not reaching genome-wide significance. A genome-wide significant association was detected for calving ease in cows on BTA23, highlighting FARS2 as a positional candidate gene. The strongest GWAS signals for each trait harbored additional biologically important candidate genes, including NPPA, NPPB, BCHE, EPHA4, DLD, and GTF2I. Given the generally low heritability estimates and the predominantly polygenic architecture observed for these traits, genomic selection may contribute to the genetic improvement of calving-related traits in Holstein cattle, with potential benefits for cow welfare, calf survival, and overall dairy production efficiency.

dairy cattle↗

PubMed2026

SNP genotyping in Pseudotsuga menziesii and Pinus radiata using targeted genotyping-by-sequencing (GBS): improved Bayesian SNP calling using a beta-binomial distribution and other optimized input parameters.

BACKGROUND: Single-nucleotide polymorphism markers (SNPs) have important applications in gene conservation, breeding, and fundamental genetics research. Our long-term goal is to develop routine approaches for SNP genotyping in forest trees. Ideally, these approaches would be inexpensive, able to accommodate a wide range of samples and SNPs, available through commercial providers, and produce high-quality SNP data. RESULTS: Using targeted genotyping-by-sequencing (GBS), we developed SNP assays for two highly heterozygous tree species, Douglas-fir (Pseudotsuga menziesii) and radiata pine (Pinus radiata). Using Douglas-fir haploid and diploid data, we optimized Bayesian SNP calling by testing four input parameters: (1) allele and genotype prior probabilities, (2) Rho, the beta-binomial dispersion parameter, (3) estimated read error (BayesReadError), and (4) the logPO cutoff used to filter low confidence SNP calls. logPO is the Bayesian posterior odds ratio for a called SNP. Compared to assuming a binomial distribution of read counts (Rho = 0), the beta-binomial distribution (Rho = 0.33) substantially reduced call error and heterozygote undercalling. Compared to the other Bayesian parameters, genotype priors had little effect on genotyping success. For Douglas-fir, we tested 5,360 SNP assays, and then studied the performance of the best 4,000. For radiata pine, we tested 6,000 SNP assays, and then studied the performance of the best 4,570. In Douglas-fir and radiata pine, our Bayesian approach resulted in median call rates of 95% to 98% for the top-ranked SNPs, with an estimated call error of 1.60% for known homozygous genotypes and 2.27% for known heterozygotes. In radiata pine, median and mean call rates were above 91% for GBS and SNP genotyping using an Axiom fixed genotyping array. Additionally, the median correspondence between the GBS and Axiom genotypes was about 98% overall (mean 96%). CONCLUSIONS: By optimizing Bayesian SNP calling, selecting the best 4-5 K SNPs, and excluding samples with low DNA amounts, we substantially reduced call error and heterozygote undercalling, resulting in SNP genotypes that were nearly identical to genotypes obtained using the Axiom array. Furthermore, genotyping performance should increase even further if our SNP rankings were used to develop less complex probe pools that target fewer SNPs.

Pinus↗

PubMed2026

Non-destructive prediction of lead content in oilseed rape leaves by fluorescence hyperspectral technology based on neural network.

Based on fluorescence hyperspectral imaging (FHSI), this study targeted rapid, non-destructive quantification of lead (Pb) content in oilseed rape leaves treated with varying silicon (Si) concentrations, acquiring fluorescence spectra over the 484.43-1001.61 nm wavelength range. To optimize spectral data quality, preprocessing methods (Savitzky-Golay smoothing, first derivative, detrending) were comprehensively compared. Characteristic wavelengths were then selected via interval variable iterative shrinkage, which effectively compressed data dimensionality and reduced computational load. A hybrid SE-CL1DA model, fusing a 1D convolutional neural network, a long short-term memory network and SE attention mechanism was constructed, with Bayesian optimization tuning hyperparameters to boost stability. The BO-SE-CL1DA outperformed both traditional machine learning and insufficiently optimized deep learning model (Rp2=0.9609, RMSE = 0.0377 mg/kg, RPD = 5.1736), thus enabling accurate Pb estimation, supporting Si-regulated heavy metal stress management and facilitating agricultural contamination monitoring.

Plant Leaves↗

PubMed2025

IsoBayes: a Bayesian approach for single-isoform proteomics inference.

MOTIVATION: Studying protein isoforms is an essential step in biomedical research; at present, the main approach for analyzing proteins is via bottom-up mass spectrometry proteomics, which return peptide identifications, that are indirectly used to infer the presence of protein isoforms. However, the detection and quantification processes are noisy; in particular, peptides may be erroneously detected, and most peptides, known as shared peptides, are associated to multiple protein isoforms. As a consequence, studying individual protein isoforms is challenging, and inferred protein results are often abstracted to the gene-level or to groups of protein isoforms. RESULTS: Here, we introduce IsoBayes, a novel statistical method to perform inference at the isoform level. Our method enhances the information available, by integrating mass spectrometry proteomics and transcriptomics data in a Bayesian probabilistic framework. To account for the uncertainty in the measurement process, we propose a two-layer latent variable approach: first, we sample if a peptide has been correctly detected (or, alternatively filter peptides); second, we allocate the abundance of such selected peptides across the protein(s) they are compatible with. This enables us, starting from peptide-level data, to recover protein-level data; in particular, we: (i) infer the presence/absence of each protein isoform (via a posterior probability), (ii) estimate its abundance (and credible interval), and (iii) target isoforms where transcript and protein relative abundances significantly differ. We benchmarked our approach in simulations, and in two multi-protease real datasets: our method displays good sensitivity and specificity when detecting protein isoforms, its estimated abundances highly correlate with the ground truth, and can detect changes between protein and transcript relative abundances. AVAILABILITY AND IMPLEMENTATION: IsoBayes is freely distributed as a Bioconductor R package, and is accompanied by an example usage vignette.

Proteomics↗

PubMed2026

Multi-omics integration and colocalization analyses prioritize candidate molecular loci associated with hypothermia.

BACKGROUND: Hypothermia is a life-threatening condition lacking specific pharmacological treatments. This study aimed to prioritize genetically supported molecular loci associated with hypothermia and to explore their pharmacological tractability using multi-omics data. METHODS: Initially, 2532 druggable genes were curated from the Drug-Gene Interaction Database and established literature. These were cross-referenced with cis-eQTL and cis-pQTL datasets, encompassing 870,655 and 114,281 SNPs for blood, respectively, alongside 2379 shared SNPs across adipose, skeletal muscle, and heart tissues. Matched instrumental variables were integrated with hypothermia GWAS summary statistics for two-sample Mendelian randomization (MR) and Bayesian colocalization. Transcriptomic differential expression analysis (DEA) was subsequently conducted as an exploratory analysis of cold-exposure-associated expression changes. Database-derived compound annotations were systematically re-evaluated according to target specificity, established pharmacological mechanism, and concordance with the direction of the MR estimates. RESULTS: Among 671 gene-level MR tests, 36 genes reached nominal significance, whereas only ABCC8 remained significant after FDR correction. Colocalization was evaluable for 8 of these 36 genes, and 4 loci (COL18A1, SLC1A7, ADIPOQ, and MERTK) met the prespecified PP.H4>0.90 threshold. The remaining 28 loci were not evaluable because sufficient overlapping regional variants were unavailable after harmonization. Transcriptomic analysis identified altered expression of SLC1A3 and SLCO4A1 under cold exposure, although these findings did not directly validate the colocalization-supported loci. Re-evaluation of database-derived compound annotations did not identify any direct, selective, and directionally concordant drug-repurposing candidate for hypothermia. CONCLUSIONS: COL18A1, SLC1A7, ADIPOQ, and MERTK showed colocalization support among the 8 evaluable nominal MR-associated loci. Because colocalization coverage was limited, these genes should be regarded as preliminary candidate loci rather than established therapeutic targets. The pharmacological annotations were indirect, non-selective, unsupported, or directionally inconsistent and should be interpreted solely as hypothesis-generating information.

Bayesian colocalization↗

PubMed2026

Bayesian Mendelian randomization reveals a protective effect of later age at first sexual intercourse against erectile dysfunction.

Erectile dysfunction (ED) is a prevalent health condition with significant psychosocial impacts, yet the causal role of age at first sexual intercourse (AFS) remains unclear. This study investigated the causal effect of AFS on the risk of ED using Mendelian randomization (MR) and Bayesian methods. Five traditional 2-sample MR analyses and 5 Bayesian MR analyses were performed using genome-wide association studies summary statistics from European populations. Sensitivity analyses included MR Egger regression, MR-pleiotropy residual sum and outlier, and Cochran Q-test. In mixed-sex cohorts (Groups 1 and 2), inverse variance weighted results demonstrated significant protective effects: odds ratio (OR) = 0.626, θ = -0.469, P = 2.73 × 10-6 for Group 1 and OR = 0.617, θ = -0.483, P = 3.56 × 10-5 for Group 2. The analyses for male-specific cohorts (Groups 3-10) showed weaker but consistent effects. For Group 3, OR = 0.643, θ = -0.442, P = .010. For Group 4, some instrumental variables associated with confounders were removed. The result became statistically insignificant: OR = 0.680, θ = -0.385, P = .064. For Group 5, the instrument selection criteria were relaxed and significance was retained: OR = 0.695, θ = -0.364, P = .016. For Groups 6 to 10, Bayesian MR was used to strengthen the inferences. In particular, for Group 8, which has a strongly informed prior, a posterior mean θ = -0.358 and a 95% credible interval (-0.575, -0.136) were obtained. This study provides evidence supporting a causal protective effect of later AFS on ED risk. While traditional MR analyses in male-specific cohorts yielded suggestive results, Bayesian MR analyses, which allow for the integration of prior evidence, provided more precise estimates and strengthened the causal inference. These findings may inform future sexual health policies. Strengths include the use of male-specific cohorts and Bayesian enhancement for weak instruments. Limitations include reliance on European-ancestry data and inability to stratify ED subtypes.

Male↗

PubMed2026

Polygenic Risk Scores for Incident Dementia in the Multi-Ethnic Study of Atherosclerosis.

Over 75 Alzheimer's disease (AD) and dementia-associated variants have been identified through genome-wide association studies, but the utility of polygenic risk scores (PRS) for predicting AD and dementia in diverse and admixed populations remains unclear. We compared how PRS approaches differing in p-value thresholds, variant weights, and source ancestry perform in predicting dementia in 6338 African American, Chinese, Hispanic, and White individuals from the Multi-Ethnic Study of Atherosclerosis. We tested clumping and thresholding (C+T) methods with varying parameters against Bayesian approaches (PRS-CS, PRS-CSx). We compared the ability of each method to predict incident dementia in all participants and in groups stratified by self-reported race/ethnicity. We additionally analyzed performance across groups stratified by estimated proportion of non-Finnish European (NFE)-like ancestry. Including more variants does not improve performance. We found comparable associations between dementia and PRS when comparing a C+T method with only 15 SNPs and PRS derived from Bayesian models that include > 800,000 SNPs (HR5e-08 = 1.18, 95% CI: 1.08-1.28; HRCSx = 1.17, 95% CI: 1.07-1.27). The p < 5e-08 C+T method was more strongly associated with incident dementia in populations genetically dissimilar from the source data (HRlowNFE_5e-08 = 1.27, 95% CI: 1.08-1.50; HRlowNFE_CSx = 1.12, 95% CI: 0.94-1.33). More selective PRS models using genome-wide significant SNPs may be preferable for dementia prediction in diverse populations.

Aged↗

PubMed2026

Comparative effectiveness of game-based learning modalities in nursing and medical education: a systematic review and Bayesian network meta-analysis.

BACKGROUND: Game-based learning (GBL) is increasingly used in healthcare education, but educators must choose among diverse modalities (e.g., quiz platforms, apps, serious games and metaverse environments). Comparative evidence on which modalities perform best across learning domains (knowledge, attitudes, and practice) remains limited. AIM: To compare the effects of distinct GBL modalities on knowledge, attitudes, and practice outcomes in nursing and medical education and to explore whether comparative effects differ by learner group (pre-licensure students and in-service professionals). DESIGN: PRISMA-NMA-aligned systematic review and Bayesian network meta-analysis. METHODS: We searched eight databases and trial registries through September 2, 2024, for randomized controlled trials comparing GBL with traditional teaching (TT). Outcomes were transformed to a 0-100 scale and analysed as change from baseline in Bayesian consistency models; random-effects models were selected using deviance information criterion (DIC). Risk of bias was assessed using RoB 2. We report mean differences (MDs) with 95% credible intervals (CrIs) versus TT, ranking probabilities, and subgroup NMAs by learner group. RESULTS: Thirty-one RCTs (n = 3439) were included; 15 contributed complete data to the network. Risk of bias was low in 15 trials and raised some concerns in 16. The network was modest for knowledge (11 trials) and sparse for attitudes (3) and practice (4). Compared with TT, metaverse-based learning showed improved attitudes (MD 15; 95% CrI 12 to 18), based on a single trial. For knowledge and practice, Kahoot-based quizzes (MD 9.1; 95% CrI -8.9 to 27) and app-based learning (MD 4.6; 95% CrI -4.4 to 14) had the highest estimated mean improvements, but credible intervals were wide and included the null for most comparisons. Subgroup rankings differed by learner group, but several comparisons were imprecise and uncertainty was substantial, particularly in sparse networks. CONCLUSIONS: GBL modalities may improve learning outcomes compared with TT, but relative effects appear domain-specific and the certainty of rankings is limited by sparse evidence and imprecision. Future trials should prioritise head-to-head comparisons, robust outcome measurement, and longer-term retention and transfer outcomes in both student and in-service populations.

Humans↗

PubMed2026

Multistage Genetic, Transcriptomic, and Single-Cell Evidence Prioritizes MAP1LC3A among Ferroptosis-Related Genes in Glioblastoma.

Glioblastoma (GBM) remains a highly aggressive malignancy, and the contribution of ferroptosis-related genes to disease susceptibility remains incompletely understood. A genetically anchored, multistage framework was applied to prioritize ferroptosis-related genes associated with GBM. Among 483 genes curated from FerrDb V2, 315 had candidate cis-expression quantitative trait loci (cis-eQTLs) in eQTLGen, 250 retained at least three independent instruments after linkage disequilibrium clumping, and 226 yielded valid inverse-variance weighted (IVW) Mendelian randomization estimates using a GBM genome-wide association study comprising 6,183 cases and 18,169 controls. Thirty-four genes met the exploratory discovery criteria of P < 0.05 and a Benjamini-Hochberg false discovery rate (BH-FDR) < 0.20, with directionally concordant Bayesian weighted Mendelian randomization (BWMR) estimates. Replication-stage Mendelian randomization using GTEx V10 whole-blood cis-eQTLs supported four genes: ATG7, RPTOR, MAP1LC3A, and CHMP6. Evaluation across three independent tumor-control transcriptomic cohorts demonstrated that MAP1LC3A was consistently downregulated in tumor tissue and showed a significant random-effects pooled estimate (log₂ fold change, -1.273; 95% confidence interval, -1.625 to -0.920; false discovery rate = 0.016), whereas the other three genes lacked comparable cross-cohort statistical support. Single-cell virtual knockout analysis was subsequently performed in a patient-balanced subset of 2,400 malignant cells selected from 4,916 eligible cells across 20 adult IDH-wild-type GBM tumors. Across five independently seeded runs, 3, 15, 4, and 7 robust downstream genes were identified for ATG7, RPTOR, MAP1LC3A, and CHMP6, respectively. The resulting consensus sets comprised 17 unique genes, with RND3 shared across all four targets. Gene Ontology analysis indicated enrichment of cell-adhesion and cell-surface processes, whereas no KEGG or Reactome pathways remained significant after multiple-testing correction. Collectively, these findings prioritize MAP1LC3A for future experimental investigation while distinguishing genetic association, tumor-expression concordance, and computational perturbation from definitive evidence of causality or mechanism.

Humans↗

PubMed2026

Precision targeting of teacher burnout using network-informed ecological momentary interventions.

Teacher well-being affects classroom functioning and workforce stability, yet generic digital programs rarely use person-specific affect dynamics to select support. This cluster-randomised trial evaluated whether micro-interventions selected from high expected influence (EI) nodes in teachers' contemporaneous affect networks produced larger changes in burnout-related EI and everyday happiness than content-matched random allocation. The objectives were to estimate allocation effects on changes in estimated network summaries and happiness, evaluate network change as a statistical mediator, examine personality moderation, and benchmark simpler allocation rules. A two-arm cluster randomised platform trial was conducted in 84 public schools across four urban districts in H Province. After a 14 day baseline of ecological momentary assessment (EMA), person specific partial correlation networks were estimated for happiness, exhaustion, detachment, efficacy and rumination. An optimisation engine prioritised three brief micro-intervention types per teacher according to baseline EI, while the active control received the same library without network information. EMA continued for 8 weeks; Bayesian multilevel models, permutation-based mediation, and benchmarking analyses were applied. EI-based targeting produced larger reductions in the composite EI-change index than active control (mean difference 0.11, 95% credible interval 0.08 to 0.14) and higher week 7 EMA happiness (4.4 points on a 0 to 100 scale, 95% credible interval 2.7 to 6.0), with a positive arm by week slope difference of 0.62 points per week (95% credible interval 0.39 to 0.85). Model-based mediation estimates were consistent with approximately one half of the happiness difference being statistically associated with change in the composite EI-change index (average conditional mediation estimate 3.5 points, 95% credible interval 2.0 to 5.2). Benchmarking showed smaller gains under severity, threshold, or group-level centrality rules. Effects were stronger among teachers higher in conscientiousness. The findings indicate that integrating EMA, network modelling, and EI-driven optimisation yields measurable gains beyond content-matched exposure, providing a proof of concept for district-scale precision mental health that requires prospective implementation testing. Replication in additional regions, expanded node sets, and longer follow up are warranted to assess durability and generalisability.

Female↗

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 = 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 = 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 < 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↗

PubMed2022

Comparison on Major Gene Mutations Related to Rifampicin and Isoniazid Resistance between Beijing and Non-Beijing Strains of Mycobacterium tuberculosis: A Systematic Review and Bayesian Meta-Analysis.

Objective: The Beijing strain of Mycobacterium tuberculosis (MTB) is controversially presented as the predominant genotype and is more drug resistant to rifampicin and isoniazid compared to the non-Beijing strain. We aimed to compare the major gene mutations related to rifampicin and isoniazid drug resistance between Beijing and non-Beijing genotypes, and to extract the best evidence using the evidence-based methods for improving the service of TB control programs based on genetics of MTB. Method: Literature was searched in Google Scholar, PubMed and CNKI Database. Data analysis was conducted in R software. The conventional and Bayesian random-effects models were employed for meta-analysis, combining the examinations of publication bias and sensitivity. Results: Of the 8785 strains in the pooled studies, 5225 were identified as Beijing strains and 3560 as non-Beijing strains. The maximum and minimum strain sizes were 876 and 55, respectively. The mutations prevalence of rpoB, katG, inhA and oxyR-ahpC in Beijing strains was 52.40% (2738/5225), 57.88% (2781/4805), 12.75% (454/3562) and 6.26% (108/1724), respectively, and that in non-Beijing strains was 26.12% (930/3560), 28.65% (834/2911), 10.67% (157/1472) and 7.21% (33/458), separately. The pooled posterior value of OR for the mutations of rpoB was 2.72 ((95% confidence interval (CI): 1.90, 3.94) times higher in Beijing than in non-Beijing strains. That value for katG was 3.22 (95% CI: 2.12, 4.90) times. The estimate for inhA was 1.41 (95% CI: 0.97, 2.08) times higher in the non-Beijing than in Beijing strains. That for oxyR-ahpC was 1.46 (95% CI: 0.87, 2.48) times. The principal patterns of the variants for the mutations of the four genes were rpoB S531L, katG S315T, inhA-15C > T and oxyR-ahpC intergenic region. Conclusion: The mutations in rpoB and katG genes in Beijing are significantly more common than that in non-Beijing strains of MTB. We do not have sufficient evidence to support that the prevalence of mutations of inhA and oxyR-ahpC is higher in non-Beijing than in Beijing strains, which provides a reference basis for clinical medication selection.

Isoniazid↗

PubMed2025

Navigating Sampling Bias in Discrete Phylogeographic Analysis: Assessing the Performance of an Adjusted Bayes Factor.

Bayesian phylogeographic inference is widely used in molecular epidemiological studies to reconstruct the dispersal history of pathogens. Discrete phylogeographic analysis treats geographic locations as discrete traits and infers lineage transition events among them, and is typically followed by a Bayes factor (BF) test to assess the statistical support. In the standard BF (BFstd) test, the relative abundance of the involved trait states is not considered, which can be problematic in the case of unbalanced sampling. Existing methods to correct sampling bias in discrete phylogeographic analyses using continuous-time Markov chain (CTMC) model, often require additional epidemiological information to balance the sampling effort among locations. As such data is not necessarily available, alternative approaches that rely solely on available genomic data are needed. In this perspective, we assess the performance of a modification of the BFstd, the adjusted Bayes factor (BFadj), which incorporates information on the relative abundance of samples by location when inferring support for transition events and root location inference without requiring additional data. Using a simulation framework, we assess the statistical performance of BFstd and BFadj under varying levels of sampling bias, estimating their type I and type II error rates. Our results show that BFadj complements the BFstd by reducing type I errors at the cost increasing type II errors for inferred transition events, while improving type I and type II errors in root location inference. Our findings provide guidelines for implementing the complementary BFadj to detect and mitigate sampling bias in discrete phylogeographic inference using CTMC modeling.

Bayes Theorem↗