Broad-scale adaptive genetic variation in alpine plants is driven by temperature and precipitation.

Identifying adaptive genetic variation is a challenging task, in particular in non-model species for which genomic information is still limited or absent. Here, we studied distribution patterns of amplified fragment length polymorphisms (AFLPs) in response to environmental variation, in 13 alpine plant species consistently sampled across the entire European Alps. Multiple linear regressions were performed between AFLP allele frequencies per site as dependent variables and two categories of independent variables, namely Moran's eigenvector map MEM variables (to account for spatial and unaccounted environmental variation, and historical demographic processes) and environmental variables. These associations allowed the identification of 153 loci of ecological relevance. Univariate regressions between allele frequency and each environmental factor further showed that loci of ecological relevance were mainly correlated with MEM variables. We found that precipitation and temperature were the best environmental predictors, whereas topographic factors were rarely involved in environmental associations. Climatic factors, subject to rapid variation as a result of the current global warming, are known to strongly influence the fate of alpine plants. Our study shows, for the first time for a large number of species, that the same environmental variables are drivers of plant adaptation at the scale of a whole biome, here the European Alps.

Predicting stroke through genetic risk functions: the CHARGE Risk Score Project.

BACKGROUND AND PURPOSE: Beyond the Framingham Stroke Risk Score, prediction of future stroke may improve with a genetic risk score (GRS) based on single-nucleotide polymorphisms associated with stroke and its risk factors. METHODS: The study includes 4 population-based cohorts with 2047 first incident strokes from 22,720 initially stroke-free European origin participants aged ≥55 years, who were followed for up to 20 years. GRSs were constructed with 324 single-nucleotide polymorphisms implicated in stroke and 9 risk factors. The association of the GRS to first incident stroke was tested using Cox regression; the GRS predictive properties were assessed with area under the curve statistics comparing the GRS with age and sex, Framingham Stroke Risk Score models, and reclassification statistics. These analyses were performed per cohort and in a meta-analysis of pooled data. Replication was sought in a case-control study of ischemic stroke. RESULTS: In the meta-analysis, adding the GRS to the Framingham Stroke Risk Score, age and sex model resulted in a significant improvement in discrimination (all stroke: Δjoint area under the curve=0.016, P=2.3×10(-6); ischemic stroke: Δjoint area under the curve=0.021, P=3.7×10(-7)), although the overall area under the curve remained low. In all the studies, there was a highly significantly improved net reclassification index (P<10(-4)). CONCLUSIONS: The single-nucleotide polymorphisms associated with stroke and its risk factors result only in a small improvement in prediction of future stroke compared with the classical epidemiological risk factors for stroke.

Crocidura cossyrensis Contoli, 1989 (Mammalia, Soridiae): karyotype, biochemical genetics and hybridization experiments

Crocidura cossyrensis Contoli, 1989 (Mammalia, Soricidae): karyotype, biochemical genetics and hybridization experiments. - The shrew Crocidura cossyrensis Contoli, 1989 from Pantelleria (I), a Mediterranean island 100 km south of Sicily and 70 km west from Tunisia, was investigated in order to understand its origin and its relationship with C. russula from Tunisia, Morocco and Switzerland. With the exception of a single heterozygote centric fusion, C. cossyrensis had a karyotype identical with that of C russula from Tunisia (2N = 42, NF = 70 to 72), but it differed from C russula from Morocco and Switzerland (2N = 42, NF = 60). The former have 5-6 pairs of chromosomes with small arms that are acrocentric in the latter. Genetic comparisons with allozyme data revealed small genetic distance (0.04) between C cossyrensis and C russula from Tunisia. In contrast, this eastern clade (Tunisia and Pantelleria) is separated from the western clade (Switzerland and Morocco) by a genetic distance of 0.14. A hybridization experiment between shrews from Pantelleria and Switzerland lead rapidly to an F1 generation. From 12 F1 hybrids that were backcrossed, females reproduced normally, but none of the males did so. Concluding from the results, C. cossyrensis from Pantelleria and C. russula cf. agilis from Tunisia belong to the same taxon that may have reached the differentiation of a biological species within the C. russula group. More geographic samples are needed to determine the definitive taxonomic positions of these shrews.

Comparative protein profiling identifies elongation factor-1beta and tryparedoxin peroxidase as factors associated with metastasis in Leishmania guyanensis.

Parasites of the Leishmania Viannia subgenus are major causative agents of mucocutaneous leishmaniasis (MCL), a disease characterised by parasite dissemination (metastasis) from the original cutaneous lesion to form debilitating secondary lesions in the nasopharyngeal mucosa. We employed a protein profiling approach to identify potential metastasis factors in laboratory clones of L. (V.) guyanensis with stable phenotypes ranging from highly metastatic (M+) through infrequently metastatic (M+/M-) to non-metastatic (M-). Comparison of the soluble proteomes of promastigotes by two-dimensional electrophoresis revealed two abundant protein spots specifically associated with M+ and M+/M- clones (Met2 and Met3) and two others exclusively expressed in M- parasites (Met1 and Met4). The association between clinical disease phenotype and differential expression of Met1-Met4 was less clear in L. Viannia strains from mucosal (M+) or cutaneous (M-) lesions of patients. Identification of Met1-Met4 by biological mass spectrometry (LC-ES-MS/MS) and bioinformatics revealed that M+ and M- clones express distinct acidic and neutral isoforms of both elongation factor-1 subunit beta (EF-1beta) and cytosolic tryparedoxin peroxidase (TXNPx). This interchange of isoforms may relate to the mechanisms by which the activities of EF-1beta and TXNPx are modulated, and/or differential post-translational modification of the gene product(s). The multiple metabolic functions of EF-1 and TXNPx support the plausibility of their participation in parasite survival and persistence and thereby, metastatic disease. Both polypeptides are active in resistance to chemical and oxidant stress, providing a basis for further elucidation of the importance of antioxidant defence in the pathogenesis underlying MCL.

Chromosomal rearrangement and genetic structure in the shrews of the Sorex araneus group

ABSTRACT The role of chromosomal rearrangements in the speciation process is much debated and many theoretical models have been developed. The shrews of the Sorex araneus group offer extraordinary opportunities to study the relationship between chromosomal variation and speciation. Indeed, this group of morphologically very similar species received a great deal of attention due to its karyotypic variability, which is mainly attributed to Robertsonian fusions. To explore the impact of karyotypic changes on genetic differentiation, we first studied the relationship between genetic and karyotypic structure among Alpine species and among chromosome races of the S. araneus group using Bayesian admixture analyses. The results of these analyses confirmed the taxonomic status of the studied species even though introgression can still be detected between species. Moreover, the strong spatial sub-structure highlighted the role of historical factors (e.g. geographical isolation) on genetic structure. Next, we studied gene flow at the chromosome level to address the question of the impact of chromosomal rearrangements on genetic differentiation. We used flow sorted chromosomes from three different karyotypic taxa of the S. araneus group to map microsatellite markers at the chromosóme arm level. We have been able to map 24 markers and to show that the karyotypic organisation of these taxa is well conserved, which suggests that these markers can be used for further inter-taxa studies. A general prediction of chromosomal speciation models is that genetic differentiation between two taxa should be larger across rearranged chromosomes than across chromosomes common to both taxa. We combined two approaches using mapped microsatellites to test this prediction. First, we studied the genetic differentiation among five shrew taxa placed at different evolutionary levels (i.e. within and among species). In this large scale study, we detected an overall significant difference in genetic structure between rearranged vs. common chromosomes. Moreover, this effect varied among pairwise comparisons, which allowed us to differentiate the role of the karyotypic complexity of hybrids and of the evolutionary divergence between taxa. Secondly, we compared the levels of gene flow measured across common vs. rearranged chromosomes in two karyotypically different hybrid zones (strong vs. low complexity of hybrids), which show similar levels of genetic structure. We detected a significantly stronger genetic structure across rearranged chromosomes in the hybrid zone showing the highest level of hybrid complexity. The large variance observed among loci suggested that other factors, such as the position of markers within the chromosome, also certainly affects genetic structure. In conclusion, our results strongly support the role of chromosomal rearrangements in the reproductive barrier and suggest their importance in speciation process of the S. araneus group. RESUME Le rôle des réarrangements chromosomiques dans les processus de spéciation est fortement débattu et de nombreux modèles théoriques ont été développés sur le sujet. Les musaraignes du groupe Sorex araneus présentent de nombreuses opportunités pour étudier les relations entre les variations chromosomiques et la spéciation. En effet, ce groupe d'espèces morphologiquement très proches a attiré l'attention des chercheurs en raison de sa variabilité caryotypique principalement attribuée à des fusions Robertsoniennes. Pour explorer l'impact des changements caryotypiques sur la différenciation génétique, nous avons tout d'abord étudié les relations entre la structure génétique et caryotypique de races chromosomiques et d'espèces alpine du groupe S. araneus en utilisant des analyses Bayesiennes d' « admixture ». Les résultats de ces analyses ont confirmé le statut taxonomique des espèces étudiées bien que nous ayons détecté de l'introgression entre espèces. L'observation d'une sous structure spatiale relativement forte souligne l'importance des facteurs historiques (telle que l'isolation géographique) sur la structure génétique de ce groupe. Ensuite, nous avons étudié le flux de gène au niveau des chromosomes pour aborder de manière directe la question de l'impact des réarrangements chromosomiques sur la différenciation génétique. En conséquence, nous avons utilisé des tris de chromosomes de trois taxons du groupe S. araneus pour localiser des marqueurs microsatellites au niveau du bras chromosomique. Au cours de cette étude, nous avons pu localiser 24 marqueurs et montrer une forte conservation dans l'organisation du caryotype de ces taxa. Ce résultat suggère que leur utilisation est appropriée pour des études entre taxa. Une prédiction générale à tous les modèles de spéciation chromosomique correspond à la plus grande différenciation génétique des chromosomes réarrangés que des chromosomes communs. Nous avons combiné deux approches utilisant des microsatellites localisés au niveau du bras chromosomique pour tester cette prédiction. Premièrement, nous avons étudié la différenciation génétique entre cinq taxa du groupe S. araneus se trouvant à des niveaux évolutifs différents (i.e. à l'intérieur et entre espèce). Au cours de cette étude, nous avons détecté une différenciation globale significativement plus élevée sur les chromosomes réarrangés. Cet effet varie entre les comparaisons, ce qui nous a permis de souligner le rôle de la complexité caryotypique des hybrides et du niveau de divergence évolutive entre taxa. Deuxièmement, nous avons comparé le flux de gènes des chromosomes communs et réarrangés dans deux zones d'hybridation caryotypiquement différentes (forte vs. Faible complexité des hybrides) mais présentant un niveau de différenciation génétique similaire. Ceci nous a permis de détecter une structure génétique significativement plus élevée sur les chromosomes réarrangés au centre de la zone d'hybridation présentant la plus grande complexité caryotypic. La forte variance observée entre loci souligne en outre le fait que d'autres facteurs, tel que la position du marqueur sur le chromosome, affectent probablement aussi la structure génétique mesurée. En conclusion, nos résultats supportent fortement le rôle des réarrangements chromosomiques dans la barrière reproductive entre espèces ainsi que leur importance dans les processus de spéciation des musaraignes du groupe S. araneus.

The biochemistry of drug metabolism : an introduction : part 6, Inter-individual factors affecting drug metabolism.

This review is part of a series of review articles on the metabolism of drugs and other xenobiotics published in Chemistry & Biodiversity. After a thorough discussion of metabolic reactions and their enzymes, this article focuses on genetically determined differences in drug and xenobiotic metabolism. After a short introduction on the causes for genetic differences, the first focus is on species differences in drug and xenobiotic metabolism. A major chapter is then dedicated to clinically relevant genetic polymorphisms in human drug metabolism and resultant ethnic differences. The last two chapters deal with sex-dependent differences in drug metabolism and personalized pharmacotherapy related to inter-individual differences in drug metabolism.

Neurobiology of suicide: do biomarkers exist?

Clinical risk factors have a low predictive value on suicide. This may explain the increasing interest in potential neurobiological correlates and specific heritable markers of suicide vulnerability. This review aims to present the current neurobiological findings that have been shown to be implicated in suicide completers and to discuss how postmortem studies may be useful in characterizing these individuals. Data on the role of the main neurobiological systems in suicidality, such as the neurotransmitter families, hypothalamic-pituitary-adrenal axis, neurotrophic factors, and polyamines, are exposed at the different biochemical, genetic, and epigenetic levels. Some neuroanatomic and neuropathological aspects as well as their in vivo morphological and functional neuroimaging correlates are also described. Except for the serotoninergic system, particularly with respect to the polymorphism of the gene coding for the serotonin transporter (5-HTTLPR) and brain-derived neurotrophic factor, data did not converge to produce a univocal consensus. The possible limitations of currently published studies are discussed, as well as the scope for long-term prospective studies.

Early onset collagen VI myopathies: Genetic and clinical correlations.

OBJECTIVE: Mutations in the genes encoding the extracellular matrix protein collagen VI (ColVI) cause a spectrum of disorders with variable inheritance including Ullrich congenital muscular dystrophy, Bethlem myopathy, and intermediate phenotypes. We extensively characterized, at the clinical, cellular, and molecular levels, 49 patients with onset in the first 2 years of life to investigate genotype-phenotype correlations. METHODS: Patients were classified into 3 groups: early-severe (18%), moderate-progressive (53%), and mild (29%). ColVI secretion was analyzed in patient-derived skin fibroblasts. Chain-specific transcript levels were quantified by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR), and mutation identification was performed by sequencing of complementary DNA. RESULTS: ColVI secretion was altered in all fibroblast cultures studied. We identified 56 mutations, mostly novel and private. Dominant de novo mutations were detected in 61% of the cases. Importantly, mutations causing premature termination codons (PTCs) or in-frame insertions strikingly destabilized the corresponding transcripts. Homozygous PTC-causing mutations in the triple helix domains led to the most severe phenotypes (ambulation never achieved), whereas dominant de novo in-frame exon skipping and glycine missense mutations were identified in patients of the moderate-progressive group (loss of ambulation). INTERPRETATION: This work emphasizes that the diagnosis of early onset ColVI myopathies is arduous and time-consuming, and demonstrates that quantitative RT-PCR is a helpful tool for the identification of some mutation-bearing genes. Moreover, the clinical classification proposed allowed genotype-phenotype relationships to be explored, and may be useful in the design of future clinical trials.

Phytochrome interacting factors 4 and 5 redundantly limit seedling de-etiolation in continuous far-red light.

Phytochromes are red/far-red photosensors that regulate numerous developmental programs in plants. Among them, phytochrome A (phyA) is essential to enable seedling de-etiolation under continuous far-red (FR) light, a condition that mimics the environment under a dense canopy. The ecological relevance of this response is demonstrated by the high mortality rate of phyA mutant plants that germinate in deep vegetational shade. phyA signaling involves direct interaction of the photoreceptor with phytochrome-interacting factors PIF1 and PIF3, members of the bHLH transcription factor family. Here we investigated the involvement of PIF4 and PIF5 in phyA signaling, and found that they redundantly control de-etiolation in FR light. The pif4 pif5 double mutant is hypersensitive to low fluence rates of FR light. This phenotype is dependent on FR light perception by phyA, but does not rely on alterations in the phyA level. Our microarray analysis shows that PIF4 and PIF5 are part of an inhibitory mechanism that represses the expression of some light-responsive genes in the dark, and that they are also needed for full expression of several growth-related genes in the light. Unlike PIF1 and PIF3, PIF4 and PIF5 are not degraded in response to FR light, indicating that they are light-regulated by a different mechanism. Our genetic analysis suggests that this is achieved through sequestration of these PIFs by the closely related bHLH transcription factor HFR1 (long hypocotyl in FR light).

The genetic basis of sleep disorders.

The contribution of genes, environment and gene-environment interactions to sleep disorders is increasingly recognized. Well-documented familial and twin sleep disorder studies suggest an important influence of genetic factors. However, only few sleep disorders have an established genetic basis including four rare diseases that may result from a single gene mutation: fatal familial insomnia, familial advanced sleep-phase syndrome, chronic primary insomnia, and narcolepsy with cataplexy. However, most sleep disorders are complex in terms of their genetic susceptibility together with the variable expressivity of the phenotype even within a same family. Recent linkage, genome-wide and candidate gene association studies resulted in the identification of gene mutations, gene localizations, or evidence for susceptibility genes and/or loci in several sleep disorders. Molecular techniques including mainly genome-wide linkage and association studies are further required to identify the contribution of new genes. These identified susceptibility genetic determinants will provide clues to better understand pathogenesis of sleep disorders, to assess the risk for diseases and also to find new drug targets to treat and to prevent the underlying conditions. We reviewed here the role of genetic basis in most of key sleep disorders.

Genetic loci for retinal arteriolar microcirculation.

Narrow arterioles in the retina have been shown to predict hypertension as well as other vascular diseases, likely through an increase in the peripheral resistance of the microcirculatory flow. In this study, we performed a genome-wide association study in 18,722 unrelated individuals of European ancestry from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium and the Blue Mountain Eye Study, to identify genetic determinants associated with variations in retinal arteriolar caliber. Retinal vascular calibers were measured on digitized retinal photographs using a standardized protocol. One variant (rs2194025 on chromosome 5q14 near the myocyte enhancer factor 2C MEF2C gene) was associated with retinal arteriolar caliber in the meta-analysis of the discovery cohorts at genome-wide significance of P-value <5×10(-8). This variant was replicated in an additional 3,939 individuals of European ancestry from the Australian Twins Study and Multi-Ethnic Study of Atherosclerosis (rs2194025, P-value = 2.11×10(-12) in combined meta-analysis of discovery and replication cohorts). In independent studies of modest sample sizes, no significant association was found between this variant and clinical outcomes including coronary artery disease, stroke, myocardial infarction or hypertension. In conclusion, we found one novel loci which underlie genetic variation in microvasculature which may be relevant to vascular disease. The relevance of these findings to clinical outcomes remains to be determined.

Genetic Variants of Serum Uric Acid and Gout: An Analysis of > 170,000 Individuals

Background/Purpose: Gout is a common and excruciatingly painful inflammatory arthritis caused by hyperuricemia. In addition to various lifestyle risk factors, a substantial genetic predisposition to gout has long been recognized. The Global Urate Genetics Consortium (GUGC) has aimed to comprehensively investigate the genetics of serum uric acid and gout using data from _ 140,000 individuals of European-ancestry, 8,340 individuals of Indian ancestry, 5,820 African-Americans, and 15,286 Japanese. Methods: We performed discovery GWAS meta-analyses of serum urate levels (n_110,347 individuals) followed by replication analyses (n_32,813 different individuals). Our gout analysis involved 3,151 cases and 68,350 controls, including 1,036 incident gout cases that met the American College of Rheumatology Criteria. We also examined the association of gout with fractional excretion of uric acid (n_6,799). A weighted genetic urate score was constructed based on the number of risk alleles across urate-associated loci, and their association with the risk of gout was evaluated. Furthermore, we examined implicated transcript expression in cis (expression quantitative trait loci databases) for potential insights into the gene underlying the association signal. Finally, in order to further identify urate-associated genomic regions, we performed functional network analyses that incorporated prior knowledge on molecular interactions in which the gene products of implicated genes operate. Results: We identified and replicated 28 genome-wide significant loci in association with serum urate (P 5_10_8), including all previously-reported loci as well as 18 novel genetic loci. Unlike the majority of previouslyidentified loci, none of the novel loci appeared to be obvious candidates for urate transport. Rather, they were mapped to genes that encode for purine production, transcription, or growth factors with broad downstream responses. Besides SLC2A9 and ABCG2, no additional regions contained SNPs that differed significantly (P _ 5_10_8) between sexes. Urateincreasing alleles were associated with an increased risk of gout for all loci. The urate genetic risk score (ranging from 10 to 45) was significantly associated with an increased odds of prevalent gout (OR per unit increase, 1.11; 95% CI, 1.09-1.14) and incident gout (OR, 1.10; 95% CI, 1.08-1.13). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. Detailed characterization of the loci revealed associations with transcript expression and the fractional excretion of urate. Network analyses implicated the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. Conclusion: The novel genetic candidates identified in this urate/gout consortium study, the largest to date, highlight the importance of metabolic control of urate production and urate excretion. The modulation by signaling processes that influence metabolic pathways such as glycolysis and the pentose phosphate pathway appear to be central mechanisms underpinned by the novel GWAS candidates. These findings may have implications for further research into urate-lowering drugs to treat and prevent gout.

Childhood adversity and psychosis: examinig whether the association is due to genetic confounding using a monozygotic twin differences approach

Purpose: To test whether the association between childhood adversity and positive and negative psychotic experiences is due to genetic confounding. Method: Childhood adversity and psychotic experiences were assessed in a sample of 226 twins from the general population. A monozygotic (MZ) twin differences approach was used to assess possible genetic confounding. Results: In the whole sample, childhood adversity was significantly associated with positive (β =.45; SE=.16; p=.008) and negative psychotic experiences (β=.77; SE=.18; p<.01). Within-pair MZ twin differences in exposure to childhood adversity were significantly associated with differences in positive (β =.71; SE=.29; p=.016) and negative psychotic experiences (β =.98; SE=.38; p=.014) in a subsample of 86 MZ twin pairs. Conclusions: Individuals exposed to childhood adversity are more likely to report psychotic experiences. Furthermore, our findings indicate that unique environmental effects of childhood adversity contribute to the development of psychotic experiences.

Prevalence of cardiovascular risk factors in a middle-income country and estimated cost of a treatment strategy.

BACKGROUND: We assessed the prevalence of risk factors for cardiovascular disease (CVD) in a middle-income country in rapid epidemiological transition and estimated direct costs for treating all individuals at increased cardiovascular risk, i.e. following the so-called "high risk strategy". METHODS: Survey of risk factors using an age- and sex-stratified random sample of the population of Seychelles aged 25-64 in 2004. Assessment of CVD risk and treatment modalities were in line with international guidelines. Costs are expressed as USD per capita per year. RESULTS: 1255 persons took part in the survey (participation rate of 80.2%). Prevalence of main risk factors was: 39.6% for high blood pressure (> or =140/90 mmHg or treatment) of which 59% were under treatment; 24.2% for high cholesterol (> or =6.2 mmol/l); 20.8% for low HDL-cholesterol (<1.0 mmol/l); 9.3% for diabetes (fasting glucose > or =7.0 mmol/l); 17.5% for smoking; 25.1% for obesity (body mass index > or =30 kg/m2) and 22.1% for the metabolic syndrome. Overall, 43% had HBP, high cholesterol or diabetes and substantially increased CVD risk. The cost for medications needed to treat all high-risk individuals amounted to USD 45.6, i.e. 11.2 dollars for high blood pressure, 3.8 dollars for diabetes, and 30.6 dollars for dyslipidemia (using generic drugs except for hypercholesterolemia). Cost for minimal follow-up medical care and laboratory tests amounted to 22.6 dollars. CONCLUSION: High prevalence of major risk factors was found in a rapidly developing country and costs for treatment needed to reduce risk factors in all high-risk individuals exceeded resources generally available in low or middle income countries. Our findings emphasize the need for affordable cost-effective treatment strategies and the critical importance of population strategies aimed at reducing risk factors in the entire population.

Effect of mutation and genetic background on drug resistance in Mycobacterium tuberculosis.

Bacterial factors may contribute to the global emergence and spread of drug-resistant tuberculosis (TB). Only a few studies have reported on the interactions between different bacterial factors. We studied drug-resistant Mycobacterium tuberculosis isolates from a nationwide study conducted from 2000 to 2008 in Switzerland. We determined quantitative drug resistance levels of first-line drugs by using Bactec MGIT-960 and drug resistance genotypes by sequencing the hot-spot regions of the relevant genes. We determined recent transmission by molecular methods and collected clinical data. Overall, we analyzed 158 isolates that were resistant to isoniazid, rifampin, or ethambutol, 48 (30.4%) of which were multidrug resistant. Among 154 isoniazid-resistant strains, katG mutations were associated with high-level and inhA promoter mutations with low-level drug resistance. Only katG(S315T) (65.6% of all isoniazid-resistant strains) and inhA promoter -15C/T (22.7%) were found in molecular clusters. M. tuberculosis lineage 2 (includes Beijing genotype) was associated with any drug resistance (adjusted odds ratio [OR], 3.0; 95% confidence interval [CI], 1.7 to 5.6; P < 0.0001). Lineage 1 was associated with inhA promoter -15C/T mutations (OR, 6.4; 95% CI, 2.0 to 20.7; P = 0.002). We found that the genetic strain background influences the level of isoniazid resistance conveyed by particular mutations (interaction tests of drug resistance mutations across all lineages; P < 0.0001). In conclusion, M. tuberculosis drug resistance mutations were associated with various levels of drug resistance and transmission, and M. tuberculosis lineages were associated with particular drug resistance-conferring mutations and phenotypic drug resistance. Our study also supports a role for epistatic interactions between different drug resistance mutations and strain genetic backgrounds in M. tuberculosis drug resistance.