Methicillin-resistant Staphylococcus aureus resensitization by protein synthesis inhibitors : proposal for a mechanism

RESUME Staphylococcus aureus est un important pathogène à gram-positif, à la fois responsable d'infections nosocomiales et communautaires. Le S. aureus résistant à la méthicilline est intrinsèquement résistant aux bêta-lactamines, inhibiteurs de la synthèse de la paroi bactérienne, grâce à une enzyme nouvellement acquise, la protéine liant la pénicilline 2A, caractérisée par une faible affinité pour ces agents et pouvant poursuivre la synthèse de la paroi, alors que les autres enzymes sont bloquées. Ce micro-organisme a également développé des résistances contre quasiment tous les antibiotiques couramment utilisés en clinique. Parallèlement au développement de molécules entièrement nouvelles, il peut être utile d'explorer d'éventuelles caractéristiques inattendues de médicaments déjà existants, par exemple en les combinant, dans l'espoir d'un potentiel effet synergique. Comprendre les mécanismes de tels effets synergiques pourrait contribuer à la justification de leur utilisation clinique potentielle. Récemment, un effet synergique contre le S. aureus résistant à la méthicilline a été décrit entre la streptogramine quinupristine-datfopristine et les bêta-lactamines, aussi bien in vitro qu'in vivo. Le présent travail a pour but de proposer un modèle pour le mécanisme de cette interaction positive et de l'étendre à d'autres classes d'antibiotiques. Premièrement, un certain nombre de méthodes microbiologiques ont permis de mieux cerner la nature de cette interaction, en montrant qu'elle agissait spécifiquement sur le S. aureus résistant à la méthicilline et qu'elle était restreinte à l'association entre inhibiteurs de la synthèse des protéines et bêta-lactamines. Deuxièmement, L'observation de l'influence des inhibiteurs de la synthèse des protéines sur la machinerie de la paroi bactérienne, c'est-à-dire sur l'expression des protéines liant la pénicilline, responsables de la synthèse du peptidoglycan, a montré une diminution de la quantité de ta protéine liant la pénicilline 2, connue pour posséder une activité de transglycosylation, indispensable au bon fonctionnement de la protéine liant la pénicilline 2A, responsable de la résistance à la méthicilline. Troisièmement, l'analyse fine de la composition du peptidoglycan extrait de bactéries, avant ou après traitement par des inhibiteurs de la synthèse des protéines, a montré des altérations corrélant avec leur capacité à agir en synergie avec les bêta-lactamines contre S. aureus résistant à ta méthicilline. Ces altérations dans les muropeptides pourraient représenter une signature de la diminution de la quantité de la protéine liant la pénicilline 2. Le modèle mécanistique retenu considère que les inhibiteurs de la synthèse des protéines pourraient diminuer l'expression de la protéine Liant la pénicilline 2, indispensable à la résistance à la méthiciltine, et que ce déséquilibre dans les enzymes synthétisant la paroi bactérienne pourrait générer une signature dans les muropeptides. SUMMARY Staphylococcus aureus is a major gram-positive pathogen causing both hospital-acquired and community-acquired infections. Methicillin- resistant Staphylococcus aureus is intrinsically resistant to the cell wall inhibitors beta-lactams by virtue of a newly acquired cell-wall-building enzyme, tow-affinity penicillin-binding protein 2A, which can build the wall when other penicillin-binding proteins are blocked. Moreover, the microorganism has developed resistance to virtually all non-experimental antibiotics. In addition of producing entirely new molecules, it is useful to explore unexpected features of existing drugs, for example by using them in combination, expecting drug synergisms. Understanding the mechanisms of such synergisms would help justify their putative clinical utilization. Recently, a synergism between the streptogramin quinupristin-dalfopristin and beta-lactams was reported against methicillin-resistant S. aureus, both in vitro and in vivo. The present work intends to propose a model for the mechanism of this positive interaction and to extend it to other drug classes. First, microbiological experimentation helped better defining the nature of this interaction, restricting it to methicillin-resistant S. aureus, and to the association of protein synthesis inhibitors with beta-lactams. Second, the observation of inhibitors of protein synthesis influence on the cell-wall-building machinery, i.e. on the expression of penicillin-binding proteins responsible for peptidoglycan synthesis, showed a decrease in the amount of penicillin-binding protein 2, known to provide a transglycosylase activity for glycan chain elongation, indispensable for the functionality of the low-affinity penicillin-binding protein 2A responsible for methicillin resistance. Third, the fine analysis of the peptidoglycan composition purified from bacteria before or after treatment with inhibitors of protein synthesis showed alterations that correlated with their ability to synergize with beta-lactams against methicillin-resistant S. aureus. These muropeptide alterations could be the signature of decrease in the amount of penicillin-binding protein 2. The retained mechanistic model is that inhibitors of protein synthesis could decrease the expression of penicillin-binding protein 2, wich is indispensable for methicillin-resistance, and that this imbalance in cell-wall-building enzymes could generate a muropeptide signature.

Saccharomyces cerevisiae, a tool to study the β-oxidation through the synthesis of polyhydroxyalkanoates

Summary Polyhydroxyalkanoates (PHAs) represent a family of polyesters naturally synthesized by a wide variety of bacteria. Through their thermoplastic and elastomeric qualities, together with their biodegradable and renewable properties, they are predicted to be a good alternative to the petroleum- derived plastics. Nevertheless, as PHA production costs using bacteria fermentation are still too high, PHA synthesis within eukaryotic systems, such as plants, has been elaborated. Although the costs were then efficiently lowered, the yield of PHAs produced remained low. In this study, Saccharomyces cerevisae has been used as another eukaryotic model in order to reveal the steps which limit PHA production. These cells express the PHA synthase of Pseudomonas aeruginosa and the PHAs obtained were analyzed to understand the flux of fatty acids towards and through the peroxisomal β-oxidation core cycle, generating the main substrate of the PHA synthase. When S. cerevisiae wild-type cells are grown in a media containing glucose as carbon source as well as fatty acids, the PHA monomer composition is largely influenced by the nature of the external fatty acid used. Thus, even-chain PHA monomers are generated from oleic acid (18:1Δ9cis) and odd- chain PHA monomers are generated from heptadecenoic acid (17:1Δ. 10 cis). Moreover, PHA synthesis is dependent on the first two enzymes of the 0-oxidation core cycle, the acyl-CoA oxidase and the multifunctional enzyme enoyl-CoA hydratase II / R-3-hydroxyacyl-CoA dehydrogenase. S. cerevisiae mutant cells growing on oleic or heptadecenoic acid and deficient in either the R-3- hydroxyacyl-CoA dehydrogenase or in the 3-ketothiolase activity, the last β-oxidation cycle steps, surprisingly contained PHAs of predominantly even-chain monomers. This is also noticed in wild- type and mutants grown on glucose or raffinose, indicating that the substrate used for PHA synthesis is generated from the degradation of intracellular short- and medium-chain fatty acids by the 3- oxidation cycle. Inhibition of fatty acid biosynthesis by cerulenin blocks the synthesis of PHAs from intracellular fatty acids but still enables the use of extracellular fatty acids for polymer production. Together, these results uncovered the existence of a substantial futile cycle whereby short- and medium-chain intermediates of the cytoplasmic fatty acid biosynthetic pathway are directed towards the peroxisomal β-oxidation pathway. In this thesis, no increase of the yield of PHA produced could be obtained. But the PHA synthesis confirmed the carbon flux into and through the β-oxidation core cycle and unveiled the existence of novel mechanisms. It is thus a good tool to study in vivo the flux of carbons in S. cerevisiae cells. Résumé Les polyhydroxyalkanoates (PHAs) sont une famille de polyesters naturellement synthétisés par un grand nombre de bactéries. Ayant des propriétés de thermoplastiques, d'élastomères et étant des ressources biodégradables et renouvelables, les PHAs représentent une bonne alternative aux plastiques dérivés du pétrole. Pour pallier aux coûts considérables de la production de PHAs par fermentation bactérienne, la synthèse de PHAs par des systèmes eucaryotes telles les plantes a été élaborée. Les coûts ont ainsi efficacement été diminués, mais le rendement de PHAs produits reste faible. Dans cette étude, Saccharomyces cerevisiae a été utilisé comme autre modèle eucaryote pour révéler les étapes limitantes de la production de PHAs. Les PHAs obtenus dans les cellules exprimant la F'HA synthase de Pseudomonas aeruginosa ont été analysés afin de comprendre le flux d'acides gras vers et à travers le cycle péroxisomal de la β-oxidation, principal producteur du substrat de la PHA synthase. Lorsque la souche S. cerevisiae de type sauvage se développe dans un milieu contenant du glucose et des acides gras, la composition des monomères de PHAs est influencée par la nature des acides gras extracellulaires. Ainsi, les monomères pairs sont générés par l'acide oléique (18:1Δ9cis), tandis que les impairs le sont par l'acide heptadécénoïque (17:1Δ10cis). La synthèse de PHAs est dépendante des deux premières enzymes de la β-oxidation; l'acyl-CoA oxidase et l'enzyme multifonctionnelle enoyl-CoA hydratase II / R-3-hydroxyacyl-CoA déshydrogénase. Les souches mutantes ne possédant pas les activités de la R-3-hydroxyacyl-CoA déshydrogénase ou de la 3- ketothiolase contiennent, en présence d'acide oléique ou heptadécénoïque, des PHAs composés essentiellement de monomères pairs. Cela a également été observé en présence de glucose ou de raffinose uniquement. Le substrat utilisé pour la synthèse de PHAs a ainsi été généré par la dégradation d'acides gras intracellulaires à chaîne courte et moyenne via le cycle de la β-oxidation. L'inhibition de la synthèse d'acides gras par la cérulénine a bloqué la synthèse de PHAs par les acides gras internes. Ces résultats ont révélés l'existence d'un cycle futile par lequel des intermédiaires à chaîne courte et moyenne de la synthèse cytoplasmique d'acides gras sont dirigés vers le cycle péroxisomal de la β-oxidation. Dans cette étude, le rendement de PHAs produits reste inchangé, mais l'analyse des PHAs permet de confirmer le flux de carbones vers et à travers le cycle péroxisomal de la β-oxidation et l'existence de nouveaux méchanismes a été dévoilée. Cette synthèse s'avère être un bon outil pour étudier in vivo le flux de carbones dans les cellules de S. cerevisiae.

Genetic dysregulation of glutathione synthesis predicts alteration of plasma thiol redox status in schizophrenia.

Abstract Genetic studies have shown an association between schizophrenia and a GAG trinucleotide repeat (TNR) polymorphism in the catalytic subunit (GCLC) of the glutamate cysteine ligase (GCL), the key enzyme for glutathione (GSH) synthesis. The present study was aimed at analyzing the influence of a GSH dysregulation of genetic origin on plasma thiols (total cysteine, homocysteine, and cysteine-glycine) and other free amino acid levels as well as fibroblast cultures GSH levels. Plasma thiols levels were also compared between patients and controls. As compared with patients with a low-risk GCLC GAG TNR genotype, patients with a high-risk genotype, having an impaired GSH synthesis, displayed a decrease of fibroblast GSH and plasma total cysteine levels, and an increase of the oxidized form of cysteine (cystine) content. Increased levels of plasma free serine, glutamine, citrulline, and arginine were also observed in the high-risk genotype. Taken together, the high-risk genotypes were associated with a subgroup of schizophrenia characterized by altered plasma thiols and free amino acid levels that reflect a dysregulation of redox control and an increased susceptibility to oxidative stress. This altered pattern potentially contributes to the development of a biomarker profile useful for early diagnosis and monitoring the effectiveness of novel drugs targeting redox dysregulation in schizophrenia. Antioxid. Redox Signal. 15, 2003-2010.

Low transformation growth factor-β1 production and collagen synthesis correlate with the lack of hepatic periportal fibrosis development in undernourished mice infected with Schistosoma mansoni

Undernourished mice infected (UI) submitted to low and long-lasting infections by Schistosoma mansoni are unable to develop the hepatic periportal fibrosis that is equivalent to Symmers’ fibrosis in humans. In this report, the effects of the host’s nutritional status on parasite (worm load, egg viability and maturation) and host (growth curves, biology, collagen synthesis and characteristics of the immunological response) were studied and these are considered as interdependent factors influencing the amount and distribution of fibrous tissue in hepatic periovular granulomas and portal spaces. The nutritional status of the host influenced the low body weight and low parasite burden detected in UI mice as well as the number, viability and maturation of released eggs. The reduced oviposition and increased number of degenerated or dead eggs were associated with low protein synthesis detected in deficient hosts, which likely induced the observed decrease in transformation growth factor (TGF)-β1 and liver collagen. Despite the reduced number of mature eggs in UI mice, the activation of TGF-β1 and hepatic stellate cells occurred regardless of the unviability of most miracidia, due to stimulation by fibrogenic proteins and eggshell glycoproteins. However, changes in the repair mechanisms influenced by the nutritional status in deficient animals may account for the decreased liver collagen detected in the present study.

University and Network Society : research report (synthesis document)

Document de síntesi d'aquest estudi que analitza -seguint una metodologia quantitativa basada en una mostra representativa de 2.093 professors i 23.864 estudiants i reforçada amb elements qualitatius- la transició que es produeix en el sistema universitari públic català cap a un model més adaptat a les noves necessitats de la societat xarxa. Per a això, es posa especial èmfasi en l'anàlisi dels usos que es fa d'Internet (l'eina clau de la societat xarxa) en el món universitari i en les transformacions que es donen o es donaran com a conseqüència d'aquests usos.

Rapid profiling of intact glucosinolates in Arabidopsis leaves by UHPLC-QTOFMS using a charged surface hybrid column.

INTRODUCTION: The analysis of glucosinolates (GS) is traditionally performed by reverse-phase liquid chromatography coupled to ultraviolet detection after a time-consuming desulphation step, which is required for increased retention. Simpler and more efficient alternative methods that can shorten both sample preparation and analysis are much needed. OBJECTIVE: To evaluate the feasibility of using ultrahigh-pressure liquid chromatography coupled to quadrupole time-of-flight mass spectrometry (UHPLC-QTOFMS) for the rapid profiling of intact GS. METHODOLOGY: A simple and short extraction of GS from Arabidopsis thaliana leaves was developed. Four sub-2 µm reverse-phase columns were tested for the rapid separation of these polar compounds using formic acid as the chromatographic additive. High-resolution QTOFMS was used to detect and identify GS. RESULTS: A novel charged surface hybrid (CSH) column was found to provide excellent retention and separation of GS within a total running time of 11 min. Twenty-one GS could be identified based on their accurate mass as well as isotopic and fragmentation patterns. The method was applied to determine the changes in GS content that occur after herbivory in Arabidopsis. In addition, we evaluated its applicability to the profiling of other Brassicaceae species. CONCLUSION: The method developed can profile the full range of GS, including the most polar ones, in a shorter time than previous methods, and is highly compatible with mass spectrometric detection.

A hybrid method for sag source location in power network

The work presented in this paper belongs to the power quality knowledge area and deals with the voltage sags in power transmission and distribution systems. Propagating throughout the power network, voltage sags can cause plenty of problems for domestic and industrial loads that can financially cost a lot. To impose penalties to responsible party and to improve monitoring and mitigation strategies, sags must be located in the power network. With such a worthwhile objective, this paper comes up with a new method for associating a sag waveform with its origin in transmission and distribution networks. It solves this problem through developing hybrid methods which hire multiway principal component analysis (MPCA) as a dimension reduction tool. MPCA reexpresses sag waveforms in a new subspace just in a few scores. We train some well-known classifiers with these scores and exploit them for classification of future sags. The capabilities of the proposed method for dimension reduction and classification are examined using the real data gathered from three substations in Catalonia, Spain. The obtained classification rates certify the goodness and powerfulness of the developed hybrid methods as brand-new tools for sag classification

Synthesis of bone-like structured foams

Here we present a processing route to produce multi-structured ceramic foams based on the combination of particle-stabilized foams with polymeric sponges to produce positive and negative templating structures. Polyester sponges are infiltrated with freshly produced calcium aluminate alumina foams and upon sintering either positive templating structures are produced when wetting the sponges, or negative templating foams with a percolating pore network are obtained when completely filling the sponges. Additionally, by combining different layers of these particle-stabilized foam infiltrated sponges, various different structures can be produced, including sandwich structures, pore size gradients, and ceramic bone-like structures applying to different types of bone. The particle-stabilized foams used were in situ self-hardening calcium aluminate cement enriched alumina foams to obtain crack-free samples with pore interconnections and tailorable pore sizes.

Modeling of hybrid systems by piecewise decomposition

Piecewise linear models systems arise as mathematical models of systems in many practical applications, often from linearization for nonlinear systems. There are two main approaches of dealing with these systems according to their continuous or discrete-time aspects. We propose an approach which is based on the state transformation, more particularly the partition of the phase portrait in different regions where each subregion is modeled as a two-dimensional linear time invariant system. Then the Takagi-Sugeno model, which is a combination of local model is calculated. The simulation results show that the Alpha partition is well-suited for dealing with such a system

Pathways for the synthesis of polyesters in plants; cutin, suberin, and polyhydroxyalkanoates

Plants naturally produce the lipid-derived polyester cutin, which is found in the plant cuticle that is deposited at the outermost extracellular matrix of the epidermis covering nearly all aboveground tissues. Being at the interface between the cell and the external environment, cutin and the cuticle play important roles in the protection of plants from several stresses. A number of enzymes involved in the synthesis of cutin monomers have recently been identified, including several P450s and one acyl-CoA synthetase, thus representing the first steps toward the understanding of polyester formation and, potentially, polyester engineering to improve the tolerance of plants to stresses, such as drought, and for industrial applications. However, numerous processes underlying cutin synthesis, such as a controlled polymerization, still remain elusive. Suberin is a second polyester found in the extracellular matrix, most often synthesized in root tissues and during secondary growth. Similar to cutin, the function of suberin is to seal off the respective tissue to inhibit water loss and contribute to resistance to pathogen attack. Being the main constituent of cork, suberin is a plant polyester that has already been industrially exploited. Genetic engineering may be worth exploring in order to change the polyester properties for either different applications or to increase cork production in other species. Polyhydroxyalkanoates (PHAs) are attractive polyesters of 3-hydroxyacids because of their properties as bioplastics and elastomers. Although PHAs are naturally found in a wide variety of bacteria, biotechnology has aimed at producing these polymers in plants as a source of cheap and renewable biodegradable plastics. Synthesis of PHA containing various monomers has been demonstrated in the cytosol, plastids, and peroxisomes of plants. Several biochemical pathways have been modified in order to achieve this, including the isoprenoid pathway, the fatty acid biosynthetic pathway, and the fatty acid β-oxidation pathway. PHA synthesis has been demonstrated in a number of plants, including monocots and dicots, and up to 40% PHA per gram dry weight has been demonstrated in Arabidopsis thaliana. Despite some successes, production of PHA in crop plants remains a challenging project. PHA synthesis at high level in vegetative tissues, such as leaves, is associated with chlorosis and reduced growth. The challenge for the future is to succeed in synthesis of PHA copolymers with a narrow range of monomer compositions, at levels that do not compromise plant productivity. This goal will undoubtedly require a deeper understanding of plant biochemical pathways and how carbon fluxes through these pathways can be manipulated, areas where plant "omics" can bring very valuable contributions.

Strong reproductive barriers in a narrow hybrid zone of West-Mediterranean green toads (Bufo viridis subgroup) with Plio-Pleistocene divergence

Background One key question in evolutionary biology deals with the mode and rate at which reproductive isolation accumulates during allopatric speciation. Little is known about secondary contacts of recently diverged anuran species. Here we conduct a multi-locus field study to investigate a contact zone between two lineages of green toads with an estimated divergence time of 2.7 My, and report results from preliminary experimental crosses. Results The Sicilian endemic Bufo siculus and the Italian mainland-origin B. balearicus form a narrow hybrid zone east of Mt. Etna. Despite bidirectional mtDNA introgression over a ca. 40 km North-South cline, no F1 hybrids could be found, and nuclear genomes display almost no admixture. Populations from each side of the contact zone showed depressed genetic diversity and very strong differentiation (FST = 0.52). Preliminary experimental crosses point to a slightly reduced fitness in F1 hybrids, a strong hybrid breakdown in backcrossed offspring (F1 x parental, with very few reaching metamorphosis) and a complete and early mortality in F2 (F1 x F1). Conclusion Genetic patterns at the contact zone are molded by drift and selection. Local effective sizes are reduced by the geography and history of the contact zone, B. balearicus populations being at the front wave of a recent expansion (late Pleistocene). Selection against hybrids likely results from intrinsic genomic causes (disruption of coadapted sets of genes in backcrosses and F2-hybrids), possibly reinforced by local adaptation (the ranges of the two taxa roughly coincide with the borders of semiarid and arid climates). The absence of F1 in the field might be due to premating isolation mechanisms. Our results, show that these lineages have evolved almost complete reproductive isolation after some 2.7 My of divergence, contrasting sharply with evidence from laboratory experiments that some anuran species may still produce viable F1 offspring after > 20 My of divergence.

Redox dysregulation and oxidative stress in schizophrenia: genetic and functional anomalies of glutathione synthesis in patients and experimental models involving GABA interneurons and NMDA receptors

Objective: Converging evidence speak in favor of an abnormal susceptibility to oxidative stress in schizophrenia. A decreased level of glutathione (GSH), the principal non-protein antioxidant and redox regulator, was observed both in cerebrospinal-fluid and prefrontal cortex of schizophrenia patients (Do et al., 2000). Results: Schizophrenia patients have an abnormal GSH synthesis most likely of genetic origin: Two independent case-control studies showed a significant association between schizophrenia and a GAG trinucleotide repeat (TNR) polymorphism in the GSH key synthesizing enzyme glutamate-cysteine-ligase (GCL) catalytic subunit (GCLC) gene. The most common TNR genotype 7/7 was more frequent in controls, whereas the rarest TNR genotype 8/8 was three times more frequent in patients. The disease-associated genotypes correlated with a decrease in GCLC protein expression, GCL activity and GSH content. Such a redox dysregulation during development could underlie the structural and functional anomalies in connectivity: In experimental models, GSH deficit induced anomalies similar to those observed in patients. (a) morphology: In animal models with GSH deficit during the development we observed in prefrontal cortex a decreased dendritic spines density in pyramidal cells and an abnormal development of parvalbumine (but not of calretinine) immunoreactive GABA interneurones in anterior cingulate cortex. (b) physiology: GSH depletion in hippocampal slices induces NMDA receptors hypofunction and an impairment of long term potentiation. In addition, GSH deficit affected the modulation of dopamine on NMDA-induced Ca 2+ response in cultured cortical neurons. While dopamine enhanced NMDA responses in control neurons, it depressed NMDA responses in GSH-depleted neurons. Antagonist of D2-, but not D1-receptors, prevented this depression, a mechanism contributing to the efficacy of antipsychotics. The redox sensitive ryanodine receptors and L-type calcium channels underlie these observations. (c) cognition: Developing rats with low [GSH] and high dopamine lead deficit in olfactory integration and in object recognition which appears earlier in males that females, in analogy to the delay of the psychosis onset between man and woman. Conclusion: These clinical and experimental evidence, combined with the favorable outcome of a clinical trial with N-Acetyl Cysteine, a GSH precursor, on both the negative symptoms (Berk et al., submitted) and the mismatch negativity in an auditory oddball paradigm supported the proposal that a GSH synthesis impairment of genetic origin represent, among other factors, one major risk factor in schizophrenia.

Synthesis and characterization of kraft lignin-based epoxy resins

Epoxidization is an interesting way to develop a new application of lignin and therefore to improve its application potential. In this work, kraft lignin-based epoxy resins were obtained by the epoxidization reaction, using the kraft lignin recovered directly from pulping liquor and modified by a methylolation reaction. The methylolated lignins were obtained by the reaction of original kraft lignin with formaldehyde and glyoxal, which is a less volatile and less toxic aldehyde. 1H-NMR spectroscopy showed that methylolated kraft lignin has more hydroxymethyl groups than glyoxalated kraft lignin. For the epoxidization reaction we studied the influence of the lignin:NaOH (w/w) ratio, temperature, and time of the reaction on the properties of the prepared epoxidized lignins. The structures of lignin-based epoxy resins were followed by epoxy index test and FTIR spectroscopy. Optimal conditions were obtained for lignin-based epoxy resin produced at lignin/NaOH = 1/3 at 70 ºC for 3h. Thermogravimetry analysis (TGA) revealed that the epoxidization enhances the thermal stability of lignins and may allow a wider temperature range for applications with lignin epoxy-PF blends

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone.

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long-standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well-studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best-supported scenario was a long period of allopatric isolation during the first three-quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.