Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Bayesian modelling”

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 1,513 records · Page 84Linked to original sources
PubMed2006

Gestational pyelonephritis as an indicator of the quality of ambulatory maternal health care services.

OBJECTIVE: Inpatient conditions that might be avoided through improved outpatient services are called Ambulatory Care Sensitive Indicators, and they include pyelonephritis in nonpregnant adults. No such indicators have been developed for pregnant women. We examine whether hospital-specific rates of gestational pyelonephritis may serve as a measure of the quality of ambulatory maternal care. METHODS: The California Department of Health Services provided an administrative data set linking maternal and newborn delivery records for 1997 with antepartum hospital admissions. We created a "low-risk" study population by largely excluding women with maternal, fetal, and placental morbidities and those with no first-trimester prenatal care. We generated hospital-specific infection rates using a Bayesian hierarchical logistic regression model. RESULTS: We identified 280,816 low-risk women, of whom 1,853 (0.66%) had at least 1 inpatient admission for gestational pyelonephritis. The model suggested only 2 risk factors: MediCal as a payer (odds ratio 1.60, 95% confidence interval 1.46-1.80 compared with all other payers), and African-American race (odds ratio 1.24, 95% confidence interval 1.10-1.41 compared with white race). Women with pyelonephritis were more than twice as likely to deliver preterm. Adjusted rates of gestational pyelonephritis for the 291 hospitals in the sample ranged from 0.22% to 2.64%. CONCLUSION: These findings suggest that because of its preventability, its consequent related morbidity, and the variation in hospital-specific rates, gestational pyelonephritis meets both clinical and technical requirements as a quality indicator for ambulatory maternal care. The use of such rates would provide an opportunity for hospitals to improve patient outcomes through partnership with obstetricians in the management of women at risk. LEVEL OF EVIDENCE: II-3.

California↗
PubMed2004

Comparing alternative measures of functional limitation.

BACKGROUND: Functional limitation is a central concept in the disability process, but its operationalization has varied widely, making it impossible to compare results across studies. OBJECTIVE: Our goal is to systematically compare the effectiveness of alternative measures of functional limitations to predict disability and to provide guidelines for their standardization. DESIGN: Over 100 alternative scales of functional limitations are compared in regression models of disability. The Bayesian information criterion is used to compare the performance of measures. SUBJECTS: The subjects were 5,764 women, aged 37 to 68, from the Mature and Young Women's cohorts of the National Longitudinal Surveys. MEASURES: Scales are constructed from 9 indicators of activity limitations: (1) grasping, (2) lifting/carrying up to 10 pounds, (3) lifting/carrying heavy weights, (4) reaching, (5) sitting, (6) stairs, (7) standing, (8) stooping, and (9) walking. Respondents were asked whether they could perform the activity without difficulty, with difficulty, or could not perform it at all. The 4 disability outcomes include binary and ordinal measures of work limitation, a single ADL measure (bathing), and total number of ADLs. RESULTS: The 3 best scales used all 9 indicators. The top scale truncated values above 5, while the others logged the sum. Two scales were based on sums of binary indicators, while the third used 3 category indicators. CONCLUSIONS: Simple scales perform better than complex scales that separate upper and lower body limitations or weight some limitations more than others. Scales that transform higher values are preferable.

Activities of Daily Living↗
PubMed2026

Bayesian inference of lineage trees by joint analysis of single-cell multimodal lineage-tracing data with BiLinT.

The advent of single-cell lineage-tracing technologies has enabled the simultaneous profiling of gene expression and lineage barcodes. However, accurate, high-resolution reconstruction of cell lineage trees remains challenging because most existing approaches treat these modalities separately and therefore fail to fully exploit their complementary information. Here we present BiLinT, a Bayesian framework that jointly models multimodal single-cell lineage-tracing data for lineage tree reconstruction. BiLinT integrates barcode evolution (a continuous-time Markov chain) with gene expression dynamics (an Ornstein-Uhlenbeck process) within a unified probabilistic model. Across synthetic and real data sets, BiLinT provides accurate lineage-tree reconstruction and reveals differentiation-associated clonal structure and developmental fate biases.

Journal Article↗
PubMed2005

Active and dynamic information fusion for facial expression understanding from image sequences.

This paper explores the use of multisensory information fusion technique with Dynamic Bayesian networks (DBNs) for modeling and understanding the temporal behaviors of facial expressions in image sequences. Our facial feature detection and tracking based on active IR illumination provides reliable visual information under variable lighting and head motion. Our approach to facial expression recognition lies in the proposed dynamic and probabilistic framework based on combining DBNs with Ekman's Facial Action Coding System (FACS) for systematically modeling the dynamic and stochastic behaviors of spontaneous facial expressions. The framework not only provides a coherent and unified hierarchical probabilistic framework to represent spatial and temporal information related to facial expressions, but also allows us to actively select the most informative visual cues from the available information sources to minimize the ambiguity in recognition. The recognition of facial expressions is accomplished by fusing not only from the current visual observations, but also from the previous visual evidences. Consequently, the recognition becomes more robust and accurate through explicitly modeling temporal behavior of facial expression. In this paper, we present the theoretical foundation underlying the proposed probabilistic and dynamic framework for facial expression modeling and understanding. Experimental results demonstrate that our approach can accurately and robustly recognize spontaneous facial expressions from an image sequence under different conditions.

Algorithms↗
PubMed2006

Bayesian analyses of admixture in wild and domestic cats (Felis silvestris) using linked microsatellite loci.

Methods recently developed to infer population structure and admixture mostly use individual genotypes described by unlinked neutral markers. However, Hardy-Weinberg and linkage disequilibria among independent markers decline rapidly with admixture time, and the admixture signals could be lost in a few generations. In this study, we aimed to describe genetic admixture in 182 European wild and domestic cats (Felis silvestris), which hybridize sporadically in Italy and extensively in Hungary. Cats were genotyped at 27 microsatellites, including 21 linked loci mapping on five distinct feline linkage groups. Genotypes were analysed with structure 2.1, a Bayesian procedure designed to model admixture linkage disequilibrium, which promises to assess efficiently older admixture events using tightly linked markers. Results showed that domestic and wild cats sampled in Italy were split into two distinct clusters with average proportions of membership Q > 0.90, congruent with prior morphological identifications. In contrast, free-living cats sampled in Hungary were assigned partly to the domestic and the wild cat clusters, with Q < 0.50. Admixture analyses of individual genotypes identified, respectively, 5/61 (8%), and 16-20/65 (25-31%) hybrids among the Italian wildcats and Hungarian free-living cats. Similar results were obtained in the past using unlinked loci, although the new linked markers identified additional admixed wildcats in Italy. Linkage analyses confirm that hybridization is limited in Italian, but widespread in Hungarian wildcats, a population that is threatened by cross-breeding with free-ranging domestic cats. The total panel of 27 loci performed better than the linked loci alone in the identification of domestic and known hybrid cats, suggesting that a large number of linked plus unlinked markers can improve the results of admixture analyses. Inferred recombination events led to identify the population of origin of chromosomal segments, suggesting that admixture mapping experiments can be designed also in wild populations.

Animals↗
PubMed2005

Prevalence and associations between bacterial isolates from dry mammary glands of dairy cows.

To assess the prevalence and patterns of bacterial isolates, cultures were made from the dry mammary glands of dairy cows in six commercial dairy herds in the UK. Milk samples were taken from all four quarters of 480 cows at drying off and at weekly intervals from 14 days before to seven days after calving. A major mastitis pathogen was isolated from at least one quarter of 220 (45.8 per cent) of the cows and from more than one quarter of 90 (18.8 per cent) of them. During the late dry to calving period, of the 957 quarters with three culture results, a major mastitis pathogen was cultured from 236 (24.7 per cent) quarters of 186 (38.8 [corrected] per cent) cows. The most commonly isolated major pathogen was Escherichia coli, followed by Streptococcus uberis and coagulase-positive staphylococci. There were significant differences between the patterns of isolates from different farms and in different calving months, suggesting that the rate of infection was partially dependent on external conditions. The isolation of E. coli, S. uberis or coagulase-positive staphylococci from a cow during the late dry/periparturient period was associated with an increased risk of that cow being culled in the next lactation. Bayesian general linear mixed models were used to assess the associations between the different bacterial species. The probability of isolating either E. coli or S. uberis was significantly greater when the other organism was cultured in a milk sample; this was also true of coagulase-positive staphylococci and S. uberis. When Corynebacterium species were isolated from a milk sample, the probability of isolating coagulase-positive staphylococci or S. uberis decreased significantly, and when coagulase-negative staphylococci were isolated the probability of isolating coagulase-positive staphylococci was reduced.

Animals↗
PubMed2006

Bayesian network to predict breast cancer risk of mammographic microcalcifications and reduce number of benign biopsy results: initial experience.

PURPOSE: To retrospectively determine whether a Bayesian network (BN) computer model can accurately predict the probability of breast cancer on the basis of risk factors and mammographic appearance of microcalcifications, to improve the positive predictive value (PPV) of biopsy, with pathologic examination and follow-up as reference standards. MATERIALS AND METHODS: The institutional review board approved this HIPAA-compliant study; informed consent was not required. Results of 111 consecutive image-guided breast biopsies performed for microcalcifications deemed suspicious by radiologists were analyzed. Mammograms obtained before biopsy were analyzed in a blinded manner by a breast imager who recorded Breast Imaging Reporting and Data System (BI-RADS) descriptors and provided a probability of malignancy. The BN uses probabilistic relationships between breast disease and mammography findings to estimate the risk of malignancy. Probability estimates from the radiologist and the BN were used to create receiver operating characteristic (ROC) curves, and area under the ROC curve (A(z)) values were compared. PPV of biopsy was also evaluated on the basis of these probability estimates. RESULTS: The BN and the radiologist achieved A(z) values of 0.919 and 0.916, respectively, which were not significantly different. If the 34 patients estimated by the BN to have less than a 10% probability of malignancy had not undergone biopsy, the PPV of biopsy would have increased from 21.6% to 31.2% without missing a breast cancer (P < .001). At this level, the radiologist's probability estimation improved the PPV to 30.0% (P < .001). CONCLUSION: A probabilistic model that includes BI-RADS descriptors for microcalcifications can distinguish between benign and malignant abnormalities at mammography as well as a breast imaging specialist can and may be able to improve the PPV of image-guided breast biopsy.

Adult↗
PubMedDate not supplied

Clinical outcomes of therapeutic agents that block the platelet glycoprotein IIb/IIIa integrin in ischemic heart disease.

BACKGROUND: Several platelet glycoprotein (GP) IIb/IIIa receptor antagonists have been evaluated in clinical trials. We conducted a systematic overview (meta-analysis) to assess the effect of these compounds on death, myocardial infarction (MI), and revascularization. METHODS AND RESULTS: ORs were calculated for 16 randomized, controlled trials of GP IIb/IIIa inhibitors. An empirical Bayesian random-effects model combined the outcomes of 32 135 patients. There was a significant mortality reduction by GP IIb/IIIa inhibitors at 48 to 96 hours, with an OR of 0.70 (95% CI, 0. 51 to 0.96; P<0.03), equivalent to a reduction of 1 death per 1000 patients treated. Mortality benefits at 30 days (OR, 0.87; 95% CI, 0. 74 to 1.02; P=0.08) and 6 months (OR, 0.97; 95% CI, 0.86 to 1.10; P=0.67) were not statistically significant. For the combined end point of death or MI, there was a highly significant (P<0.001) benefit for GP IIb/IIIa inhibitors at each time point. The 30-day OR was 0.76 (95% CI, 0.66 to 0.87), or 20 fewer events per 1000 patients treated. For the composite end point of death, MI, or revascularization, there was also a highly significant (P<0.001) benefit for GP IIb/IIIa inhibitors. At 30 days, the OR was 0.77 (95% CI, 0.68 to 0.86), or 30 fewer events per 1000 patients treated. The risk differences for death, death or MI, and composite outcomes were similar at 6 months, indicating a sustained absolute improvement. Similar benefit was seen when trials were subgrouped by therapeutic indication (percutaneous intervention versus acute coronary syndromes). CONCLUSIONS: Application of this new therapeutic class to clinical practice promises substantial benefit for both indications.

Double-Blind Method↗
PubMed2003

Hierarchical Bayesian inference in the visual cortex.

Traditional views of visual processing suggest that early visual neurons in areas V1 and V2 are static spatiotemporal filters that extract local features from a visual scene. The extracted information is then channeled through a feedforward chain of modules in successively higher visual areas for further analysis. Recent electrophysiological recordings from early visual neurons in awake behaving monkeys reveal that there are many levels of complexity in the information processing of the early visual cortex, as seen in the long-latency responses of its neurons. These new findings suggest that activity in the early visual cortex is tightly coupled and highly interactive with the rest of the visual system. They lead us to propose a new theoretical setting based on the mathematical framework of hierarchical Bayesian inference for reasoning about the visual system. In this framework, the recurrent feedforward/feedback loops in the cortex serve to integrate top-down contextual priors and bottom-up observations so as to implement concurrent probabilistic inference along the visual hierarchy. We suggest that the algorithms of particle filtering and Bayesian-belief propagation might model these interactive cortical computations. We review some recent neurophysiological evidences that support the plausibility of these ideas.

Animals↗
PubMed2004

Evidence for genetic control of adult weight plasticity in the snail Helix aspersa.

Phenotypic plasticity and canalization are important topics in quantitative genetics and evolution. Both concepts are related to environmental sensitivity. The latter can be modeled using a model with genetically structured environmental variance. This work reports the results of a genetic analysis of adult weight in the snail Helix aspersa. Several models of heterogeneous variance are fitted using a Bayesian, MCMC approach. Exploratory analyses using posterior predictive model checking and model comparisons based on the deviance information criterion favor a model postulating a genetically structured heterogeneous environmental variance. Our analysis provides a strong indication of a positive genetic correlation between additive genetic values affecting the mean and those affecting environmental variation of adult body weight. The possibility of manipulating environmental variance by selection is illustrated numerically using estimates of parameters derived from the snail data set.

Animals↗
PubMed2006

A Bayesian heterogeneous analysis of variance approach to inferring recent selective sweeps.

The distribution of microsatellite allele sizes in populations aids in understanding the genetic diversity of species and the evolutionary history of recent selective sweeps. We propose a heterogeneous Bayesian analysis of variance model for inferring loci involved in recent selective sweeps by analyzing the distribution of allele sizes at multiple loci in multiple populations. Our model is shown to be consistent with a multilocus test statistic, ln RV, proposed for identifying microsatellite loci involved in recent selective sweeps. Our methodology differs in that it accepts original allele size data rather than summary statistics and allows the incorporation of prior knowledge about allele frequencies using a hierarchical prior distribution consisting of log normal and gamma probability distributions. Interesting features of the model are its ability to simultaneously analyze allele size data for any number of populations and to cope with the presence of any number of selected loci. The utility of the method is illustrated by application to two sets of microsatellite allele size data for a group of West African Anopheles gambiae populations. The results are consistent with the suppressed-recombination model of speciation, and additional candidate loci on chromosomes 2 (079 and 175) and 3 (088) are discovered that escaped former analysis.

Alleles↗
PubMed2006

Carryover aquatic effects on survival of metamorphic frogs during pond emigration.

In organisms with complex life cycles, physiological stressors during early life stages may have fitness-level impacts that are delayed into later stages or habitats. We tested the hypothesis that body size and date of metamorphosis, which are highly responsive to aquatic stressors, influence post-metamorphic survival and movement patterns in the terrestrial phase of an ephemeral pond-breeding frog by examining these traits in two populations of northern red-legged frogs (Rana aurora aurora). To increase variation of body size at metamorphosis, we manipulated food availability for 314 of 1045 uniquely marked tadpoles and estimated the probability that frogs survived and emigrated using concentric rings of drift fencing surrounding ponds and Bayesian capture-recapture modeling. The odds of surviving and emigrating from the ponds to the innermost drift fences, approximately 12 m, increased by factors of 2.20 (95% credibility intervals 1.39-4.23) and 2.54 (0.94-4.91) with each millimeter increase in snout-vent length and decreased by factors of 0.91 (0.85-0.96) and 0.89 (0.80-1.00) with each day's delay in metamorphosis for the two ponds. The odds of surviving and moving to the next ring of fencing, 12 m to approximately 40 m from the ponds, increased by a factor of 1.20 (0.45-4.06) with each millimeter increase in size. Our results demonstrated that body size and timing of metamorphosis relate strongly to the performance of newly metamorphosed frogs during their initial transition into terrestrial habitat. Carryover effects of aquatic stressors that reduce size and delay metamorphosis may have population-level impacts that are not expressed until terrestrial stages. Since changes in both aquatic and terrestrial systems are implicated in many amphibian declines, quantifying both immediate and delayed effects of stressors on demographic rates is critical to sound management.

Age Factors↗
PubMed2004

Comparing alternative definitions of the contemporary group effect in Avilena Negra Iberica beef cattle using classical and Bayesian criteria.

Data on weaning weight from 12,740 animals were used to compare different definitions of contemporary groups (CG) for the genetic evaluation of the Avilena Negra Iberica beef cattle breed. Six alternative definitions for the CG effect were considered: herd-year-season of calving (HYS), with seasons defined according to the four natural seasons; herd-year-month of calving (HYM); herd clusters of 30 d (HC30-30) or 90 d (HC90-90); and adaptive herd clusters with two time limits, 30 and 90 d (HC30-90), and 30 and 180 d (HC30-180). A minimum of five observations in each CG class was required. This rendered substantial differences in loss of information, ranging from 0.7% of the total number of records for HC30-180 to 14% for HYM. Several classical statistics and Bayesian criteria for statistical model comparison were used. The use of classical criteria, such as the between- and within-CG variation and the accuracy of prediction, can be controversial because of their dependency on the unknown variance components. Residual variance decreased with the decrease in time span associated with the definition of CG. This was expected in this population because environmental conditions are highly variable throughout the year. However, estimates of the additive genetic variance for direct effects, which should not be affected by the definition of CG, were substantially larger for definitions involving larger time periods (HYS, HC90-90). When parameters used in the current evaluation procedure were used with all data sets, CG involving 30 d (HYM and HC30-30) were optimal in terms of providing the lowest/largest within-/between-CG variation. On the other hand, CG involving 90 d (HYS and HC90-90) yielded the poorest within-/between CG variation, with only a slight improvement of accuracy of prediction of direct genetic values over the other definitions. Bayes factors and cross-validation predictive densities allowed for improved discrimination among models. Models including CG spanning 30 d were more plausible and showed better predicting ability than models spanning 90 d. Adaptive CG showed intermediate results. Overall, it seems that average time span rendered by the different definitions had a major effect on the ranking of models. However, from the breeder's point of view, the loss of information associated with definitions involving shorter periods of time, such as HYM or HC30-30, might be unacceptable.

Animals↗
PubMed2006

Analysis of litter size and days to lambing in the Ripollesa ewe. II. Estimation of variance components and response to phenotypic selection on litter size.

A performance data set from 376 Ripollesa purebred ewes of the experimental flock of the Universitat Autònoma of Barcelona was analyzed using a bivariate Bayesian threshold-linear model. The data set contained 1,598 litter size records and 1,699 days-to-lambing records. The model included the additive genetic effect of each animal and 3 nongenetic sources of variation: ewe age, year of lambing, and the permanent environmental effect characterized by the ewe. The flock was phenotypically selected for litter size since 1986, and replacement ewes and rams were selected from the progeny of the more prolific ewes, which had at least 3 deliveries recorded. The phenotypic trend for litter size was positive, whereas days to lambing followed an unclear pattern. Both traits had low heritabilities; 0.13 for litter size and 0.11 for days to lambing. Response to selection was evaluated through (a) the average breeding value of the ewe lambs chosen annually, and (b) the average breeding value of the overall flock. The first measurement suggested a positive trend for litter size, although it showed important oscillations. On the other hand, the average breeding value for the overall flock showed a stable positive tendency after yr 4 of selection, with estimates clearly different from zero after yr 11 of selection. A significant increase in the incidence of multiple births was observed, with a mode of approximately 10%. The correlated response in days to lambing did not show a significant trend. The effect of year of lambing also positively influenced both litter size and days to lambing, although important oscillations were observed between years. Results indicated that litter size in sheep can be effectively improved through phenotypic selection, even in small flocks; moreover, days to lambing could also be genetically improved, given the estimate obtained for its heritability.

Animals↗
PubMedDate not supplied

Computer-based testing in family practice certification and recertification.

BACKGROUND: The member boards of the American Board of Medical Specialties have agreed to expand the scope of certification to include assessment of medical knowledge, practice-based learning and improvement, patient care, interpersonal and communication skills, systems-based practice, and professionalism. Multiple-choice examinations provide limited ability to assess these dimensions. METHODS: The American Board of Family Practice (ABFP) has developed a computer simulation system to facilitate more comprehensive candidate evaluation. The system consists of a knowledge base, a simulation program to create patient scenarios, an interface for presenting simulations to users, and an administrative database to track candidate performance and interactions with the system. The system uses population distributions for disease states to produce cases and evolves patients in response to candidate interventions, such as pharmacological and nonpharmacological therapies. We use Bayesian belief networks to model patient characteristics and comorbid condition interactions. RESULTS: Simulations have been created for 7 disease states; ultimately simulations will be available for 25 to 30 disease states. Initial testing will take place in regional examination centers but will ultimately use the Internet for convenient access for certification and recertification candidates. CONCLUSION: The ABFP will begin field-testing the system in early 2003 and will include simulations in the certification and recertification examination process in 2004.

Bayes Theorem↗
PubMed2006

Short communication: Genetic analysis of nonreturn rate and mastitis in first-lactation Norwegian Red cows.

Associations between clinical mastitis (CM) and nonreturn rate within 56 d after first insemination (NR56) were examined in Norwegian Red (NRF) cows. Records on absence or presence of CM within each of the intervals, -30 to 30, 31 to 150, and 151 to 300 d after first calving, and records on NR56 for 620,492 first-lactation daughters of 3,064 NRF sires were analyzed with a Bayesian multivariate threshold liability model. Point estimates of genetic correlations between NR56 and the 3 CM traits were between -0.05 and -0.02. Residual correlations were close to zero, and correlations between herd-5-yr effects on NR56 and CM in the 3 lactation intervals ranged from -0.15 to -0.17. It appears that CM and NR56 in first lactation are independent traits.

Animals↗
PubMed2005

Validation of in situ hybridisation and histology assays for the detection of the oyster parasite Marteilia refringens.

An in situ hybridisation technique has been developed for the detection of infection in oysters with Marteilia refringens with particular emphasis on light infections or confirmation of suspected cases by means of histology. Although validation of new diagnostic methods is usually achieved by comparison with standard techniques, in our case the sensitivity and specificity of the standard (histology) had not previously been established. Another point to consider is that surveillance and monitoring frequently target populations displaying different levels of prevalence under different field conditions. The objective of our study was to evaluate the sensitivity and specificity values of in situ hybridisation and histology for the detection of M. refringens, based on 3 populations of flat oysters, free of the disease and with mild and high levels of prevalence. A blind assay of 200 individuals from each population was performed using both techniques. Results were analysed by means of the classical approach and latent models (maximum likehood and Bayesian approach). Assumptions and results were found to vary slightly with the different statistical approaches. The more realistic estimate by the Bayesian approach shows a link between the level of prevalence and the sensitivity of the techniques. Values of sensitivity and specificity for histology were 0.7 and 0.99 respectively, and 0.9 and 0.99 respectively in the case of in situ hybridisation. Some uncertainty remains regarding these values because the study does not take into account the severity of infection or the developmental stages of the parasite actually present in each individual. This work provides valuable information with regard to the choice and potential use of those 2 diagnostic methods currently recommended by international standards.

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↗