Search

Antimicrobial resistance (AMR) is a major public health problem of the 21st century. The ability of some bacteria to develop resistance to specific antibiotics is the cause of an increased morbidity, mortality and health expenditure. Several surveillance programms have been introduced in the last decades to monitor the spread of antimicrobial resistance. The present work has been conducted within the JPIAMR-project MAGIcIAN whose aim is to support the sustainable introduction of novel class and last-resort antimicrobial drugs minimising the emergence of AMR. Within this project we developed a multi-level model to describe the spread of the sexually transmitted disease of gonorrhoea, caused by the Neisseria gonorrhoeae bacterium, a multidrug resistant bacteria who has progressively developed resistance to many treatment options. The multi-level model includes a dynamic sexual contact network, that describes the dynamic of sexual partnerships, a transmission model that describes the probability of infection during intercourse, and a within-host model, that describes the dynamic of gonorrhoea infection within an individual. The novelty of the proposed model is in including communities having different sexual orientations and behaviour and the possibility of these communities to interact in a dynamic framework. In this work, we calibrate the model using data coming from several clinics located in Amsterdam.

Dysregulated inflammation underlies various diseases. Specialized pro-resolving mediators (SPMs) like Resolvin D1 (RvD1) have been shown to resolve inflammation and halt disease progression. Macrophages, key immune cells that drive inflammation, respond to the presence of RvD1 by polarizing to an anti-inflammatory type ( ). However, RvD1’s mechanisms, roles, and utility are not fully understood. This paper introduces a gene-regulatory network (GRN) model that contains pathways for RvD1 and other SPMs and proinflammatory molecules like lipopolysaccharides. We couple this GRN model to a partial differential equation–agent-based hybrid model using a multiscale framework to simulate an acute inflammatory response with and without the presence of RvD1. We calibrate and validate the model using experimental data from two animal models. The model reproduces the dynamics of key immune components and the effects of RvD1 during acute inflammation. Our results suggest RvD1 can drive macrophage polarization through the G protein-coupled receptor 32 (GRP32) pathway. The presence of RvD1 leads to an earlier and increased polarization, reduced neutrophil recruitment, and faster apoptotic neutrophil clearance. These results support a body of literature that suggests that RvD1 is a promising candidate for promoting the resolution of acute inflammation. We conclude that once calibrated and validated on human data, the model can identify critical sources of uncertainty, which could be further elucidated in biological experiments and assessed for clinical use., Dysregulated inflammation underlies various diseases. Specialized pro-resolving mediators (SPMs) like Resolvin D1 (RvD1) have been shown to resolve inflammation and halt disease progression. Macrophages, key immune cells that drive inflammation, respond to the presence of RvD1 by polarizing to an anti-inflammatory type (M2). However, RvD1's mechanisms, roles, and utility are not fully understood. This paper introduces a gene-regulatory network (GRN) model that contains pathways for RvD1 and other SPMs and proinflammatory molecules like lipopolysaccharides. We couple this GRN model to a partial differential equation–agent-based hybrid model using a multiscale framework to simulate an acute inflammatory response with and without the presence of RvD1. We calibrate and validate the model using experimental data from two animal models. The model reproduces the dynamics of key immune components and the effects of RvD1 during acute inflammation. Our results suggest RvD1 can drive macrophage polarization through the G protein-coupled receptor 32 (GRP32) pathway. The presence of RvD1 leads to an earlier and increased M2 polarization, reduced neutrophil recruitment, and faster apoptotic neutrophil clearance. These results support a body of literature that suggests that RvD1 is a promising candidate for promoting the resolution of acute inflammation. We conclude that once calibrated and validated on human data, the model can identify critical sources of uncertainty, which could be further elucidated in biological experiments and assessed for clinical use.

Dysregulated inflammation underlies various diseases. Specialized pro-resolving mediators (SPMs) like Resolvin D1 (RvD1) have been shown to resolve inflammation and halt disease progression. Macrophages, key immune cells that drive inflammation, respond to the presence of RvD1 by polarizing to an anti-inflammatory type ( ). However, RvD1’s mechanisms, roles, and utility are not fully understood. This paper introduces a gene-regulatory network (GRN) model that contains pathways for RvD1 and other SPMs and proinflammatory molecules like lipopolysaccharides. We couple this GRN model to a partial differential equation–agent-based hybrid model using a multiscale framework to simulate an acute inflammatory response with and without the presence of RvD1. We calibrate and validate the model using experimental data from two animal models. The model reproduces the dynamics of key immune components and the effects of RvD1 during acute inflammation. Our results suggest RvD1 can drive macrophage polarization through the G protein-coupled receptor 32 (GRP32) pathway. The presence of RvD1 leads to an earlier and increased polarization, reduced neutrophil recruitment, and faster apoptotic neutrophil clearance. These results support a body of literature that suggests that RvD1 is a promising candidate for promoting the resolution of acute inflammation. We conclude that once calibrated and validated on human data, the model can identify critical sources of uncertainty, which could be further elucidated in biological experiments and assessed for clinical use., Dysregulated inflammation underlies various diseases. Specialized pro-resolving mediators (SPMs) like Resolvin D1 (RvD1) have been shown to resolve inflammation and halt disease progression. Macrophages, key immune cells that drive inflammation, respond to the presence of RvD1 by polarizing to an anti-inflammatory type (M2). However, RvD1's mechanisms, roles, and utility are not fully understood. This paper introduces a gene-regulatory network (GRN) model that contains pathways for RvD1 and other SPMs and proinflammatory molecules like lipopolysaccharides. We couple this GRN model to a partial differential equation–agent-based hybrid model using a multiscale framework to simulate an acute inflammatory response with and without the presence of RvD1. We calibrate and validate the model using experimental data from two animal models. The model reproduces the dynamics of key immune components and the effects of RvD1 during acute inflammation. Our results suggest RvD1 can drive macrophage polarization through the G protein-coupled receptor 32 (GRP32) pathway. The presence of RvD1 leads to an earlier and increased M2 polarization, reduced neutrophil recruitment, and faster apoptotic neutrophil clearance. These results support a body of literature that suggests that RvD1 is a promising candidate for promoting the resolution of acute inflammation. We conclude that once calibrated and validated on human data, the model can identify critical sources of uncertainty, which could be further elucidated in biological experiments and assessed for clinical use.

Background: Meloidae (blister beetles) are known to synthetize cantharidin (CA), a toxic and defensive terpene mainly stored in male accessory glands (MAG) and emitted outward through reflex-bleeding. Recent progresses in understanding CA biosynthesis and production organ(s) in Meloidae have been made, but the way in which self-protection is achieved from the hazardous accumulation and release of CA in blister beetles has been experimentally neglected. To provide hints on this pending question, a comparative de novo assembly transcriptomic approach was performed by targeting two tissues where CA is largely accumulated and regularly circulates in Meloidae: the male reproductive tract (MRT) and the haemolymph. Differential gene expression profiles in these tissues were examined in two blister beetle species, Lydus trimaculatus (Fabricius, 1775) (tribe Lyttini) and Mylabris variabilis (Pallas, 1781) (tribe Mylabrini). Upregulated transcripts were compared between the two species to identify conserved genes possibly involved in CA detoxification and transport. Results: Based on our results, we hypothesize that, to avoid auto-intoxication, ABC, MFS or other solute transporters might sequester purported glycosylated CA precursors into MAG, and lipocalins could bind CA and mitigate its reactivity when released into the haemolymph during the autohaemorrhaging response. We also found an over-representation in haemolymph of protein-domains related to coagulation and integument repairing mechanisms that likely reflects the need to limit fluid loss during reflex-bleeding. Conclusions: The de novo assembled transcriptomes of L. trimaculatus and M. variabilis here provided represent valuable genetic resources to further explore the mechanisms employed to cope with toxicity of CA in blister beetle tissues. These, if revealed, might help conceiving safe and effective drug-delivery approaches to enhance the use of CA in medicine.

Background: Meloidae (blister beetles) are known to synthetize cantharidin (CA), a toxic and defensive terpene mainly stored in male accessory glands (MAG) and emitted outward through reflex-bleeding. Recent progresses in understanding CA biosynthesis and production organ(s) in Meloidae have been made, but the way in which self-protection is achieved from the hazardous accumulation and release of CA in blister beetles has been experimentally neglected. To provide hints on this pending question, a comparative de novo assembly transcriptomic approach was performed by targeting two tissues where CA is largely accumulated and regularly circulates in Meloidae: the male reproductive tract (MRT) and the haemolymph. Differential gene expression profiles in these tissues were examined in two blister beetle species, Lydus trimaculatus (Fabricius, 1775) (tribe Lyttini) and Mylabris variabilis (Pallas, 1781) (tribe Mylabrini). Upregulated transcripts were compared between the two species to identify conserved genes possibly involved in CA detoxification and transport. Results: Based on our results, we hypothesize that, to avoid auto-intoxication, ABC, MFS or other solute transporters might sequester purported glycosylated CA precursors into MAG, and lipocalins could bind CA and mitigate its reactivity when released into the haemolymph during the autohaemorrhaging response. We also found an over-representation in haemolymph of protein-domains related to coagulation and integument repairing mechanisms that likely reflects the need to limit fluid loss during reflex-bleeding. Conclusions: The de novo assembled transcriptomes of L. trimaculatus and M. variabilis here provided represent valuable genetic resources to further explore the mechanisms employed to cope with toxicity of CA in blister beetle tissues. These, if revealed, might help conceiving safe and effective drug-delivery approaches to enhance the use of CA in medicine.

Cantharidin (CA), a toxic terpene produced by blister beetles (Coleoptera: Meloidae), attracts the interest of many researchers for its renowned medical properties. The CA content in blister beetles has been mainly quantified in some Oriental species, due to their use in traditional Chinese medicine, or in few other species of toxicological importance. As CA quantification has been largely ignored in many other members of this family, we aimed at estimating CA content in natural populations of two of the most abundant species in central Italy, Mylabris variabilis and Lydus trimaculatus. Nearly 100 individuals for each species were collected in the field, identified and reared in fauna-boxes until CA collection. Available protocols were optimized for CA extraction from both dried-body tissues and exuded hemolymph collected from leg joints using capillary tubes. The CA content was quantified in a gas chromatography system coupled to a mass spectrometer using a calibration curve with diethyl-ester of norcantharidin as the internal standard. We observed: (i) high variability in CA content among specimens, which was positively related to individual dry-weight; and (ii) slightly larger (though non-significant) amount of CA in males than in females for both species. Our data are consistent with the available published reports about CA content, transfer and distribution in blister beetles and suggest a conserved biological role of this terpene in Meloidae.

Background: Clinical decision support systems (CDSS) are a category of health information technologies that can assist clinicians to choose optimal treatments. These support systems are based on clinical trials and expert knowledge; however, the amount of data available to these systems is limited. For this reason, CDSSs could be significantly improved by using the knowledge obtained by treating patients. This knowledge is mainly contained in patient records, whose usage is restricted due to privacy and confidentiality constraints. Methods: A treatment effectiveness measure, containing valuable information for treatment prescription, was defined and a method to extract this measure from patient records was developed. This method uses an advanced cryptographic technology, known as secure Multiparty Computation (henceforth referred to as MPC), to preserve the privacy of the patient records and the confidentiality of the clinicians' decisions. Results: Our solution enables to compute the effectiveness measure of a treatment based on patient records, while preserving privacy. Moreover, clinicians are not burdened with the computational and communication costs introduced by the privacy-preserving techniques that are used. Our system is able to compute the effectiveness of 100 treatments for a specific patient in less than 24 minutes, querying a database containing 20,000 patient records. Conclusion: This paper presents a novel and efficient clinical decision support system, that harnesses the potential and insights acquired from treatment data, while preserving the privacy of patient records and the confidentiality of clinician decisions., Comment: 15 pages, 4 figures

Resolving complexes of closely related and cryptic insect species can be challenging, especially when dealing with rare and protected taxa that are difficult to collect for genetic and morphological analyses. Until recently, populations of the genus Osmoderma (Scarabaeidae), widespread in Europe, were treated as a single species O.eremita (Scopoli, 1763) in spite of observed geographic variation in morphology. A previous survey using sequence data from the mtDNA cytochrome C oxidase I gene (COI) revealed the occurrence of at least two distinct lineages within this species complex: O.eremita in the west and O.barnabita Motschulsky, 1845, in the east. Interestingly, beetles confined to Sicily have been described as a distinct species, O.cristinae Sparacio, 1994, based on morphological traits. Only few Sicilian specimens were included in the former genetic analysis, and the results led to a still questionable taxonomic rank for these populations. To explore the robustness of the previous taxonomic arrangement for O.cristinae, a combination of genetic, morphological and pheromonal analyses was used. A 617-bp fragment of the COI gene, aligned with O.cristinae and O.eremita sequences already available in GenBank, showed a clear genetic divergence between the two species (interspecific mean distance=6.6%). Moreover, results from AFLP markers sustained the distinction of the two species. In addition, geometric morphometric analyses of the shape of male genitalia revealed a clear differentiation between the two species. Via scent analysis and field trapping, we demonstrated the production of the sex pheromone (R)-(+)--decalactone by males of O.cristinae, the attraction by conspecific individuals (mostly females) to this compound, and a lack of antagonistic effect of (S)-(-)--decalactone. The fact that O.eremita and O.eremita use the same compound for mate finding suggests that this sex pheromone has not undergone a differentiation and probably the allopatry of these two spec

A new species of Epuraea, Epuraea sutcuimamun Avgin, Lason andamp; Audisio sp. n., from southern Turkey (Taurus Chain) was identified using morphological analyses. This species is a member of a circum-Mediterranean endemic group of sap beetles, the Epuraea latipes species group of the Dadopora lineage, which are mostly associated with cedar (Cedrus spp.) forests. Herein, the new species is described, the previously unknown association of the little-known Epuraea subparallela Grouvelle 1896 with meso-xerophilous oaks is reported, and some new records from southern Turkey (Taurus Chain) are listed for the latter species (previously only known from the Nur Mts. in SE Turkey, Osmaniye and Iskenderun provinces). An identification key to species of Dadopora lineage is also provided. Finally, the main palaeogeographic events affecting the Miocene to Holocene dispersal and evolution of species/populations of the Dadopora lineage on cedars and oaks, and the coincident distributional dynamics of Western Palaearctic Cedrus populations throughout the Mediterranean Basin, are discussed., Funding Agencies|M.I.U.R. - "Sapienza, University of Rome 60% funds"||Centro Nazionale per la Biodiversita Forestale, Verona, Italy

The main Afrotropical malaria vector, Anopheles gambiae sensu stricto, is undergoing a process of sympatric ecological diversification leading to at least two incipient species (the M and S molecular forms) showing heterogeneous levels of divergence across the genome. The physically unlinked centromeric regions on all three chromosomes of these closely related taxa contain fixed nucleotide differences which have been found in nearly complete linkage disequilibrium in geographic areas of no or low M-S hybridization. Assays diagnostic for SNP and structural differences between M and S forms in the three centromeric regions were applied in samples from the western extreme of their range of sympatry, the only area where high frequencies of putative M/S hybrids have been reported. The results reveal a level of admixture not observed in the rest of the range. In particular, we found: i) heterozygous genotypes at each marker, although at frequencies lower than expected under panmixia; ii) virtually all possible genotypic combinations between markers on different chromosomes, although genetic association was nevertheless detected; iii) discordant M and S genotypes at two X-linked markers near the centromere, suggestive of introgression and inter-locus recombination. These results could be indicative either of a secondary contact zone between M and S, or of the maintenance of ancestral polymorphisms. This issue and the perspectives opened by these results in the study of the M and S incipient speciation process are discussed.

Background: During copulation, the major Afro-tropical malaria vector Anopheles gambiae s.s. transfers male accessory gland (MAG) proteins to females as a solid mass (i.e. the “mating plug”). These proteins are postulated to function as important modulators of female post-mating responses. To understand the role of selective forces underlying the evolution of these proteins in the A. gambiae complex, we carried out an evolutionary analysis of gene sequence and expression divergence on a pair of paralog genes called AgAcp34A-1 and AgAcp34A-2. These encode MAG-specific proteins which, based on homology with Drosophila, have been hypothesized to play a role in sperm viability and function. Results: Genetic analysis of 6 species of the A. gambiae complex revealed the existence of a third paralog (68-78% of identity), that we named AgAcp34A-3. FISH assays showed that this gene maps in the same division (34A) of chromosome-3R as the other two paralogs. In particular, immuno-fluorescence assays targeting the C-terminals of AgAcp34A-2 and AgAcp34A-3 revealed that these two proteins are localized in the posterior part of the MAG and concentrated at the apical portion of the mating plug. When transferred to females, this part of the plug lies in proximity to the duct connecting the spermatheca to the uterus, suggesting a potential role for these proteins in regulating sperm motility. AgAcp34A-3 is more polymorphic than the other two paralogs, possibly because of relaxation of purifying selection. Since both unequal crossing-over and gene conversion likely homogenized the members of this gene family, the interpretation of the evolutionary patterns is not straightforward. Although several haplotypes of the three paralogs are shared by most A. gambiae s.l. species, some fixed species-specific replacements (mainly placed in the N- and C-terminal portions of the secreted peptides) were also observed, suggesting some lineage-specific adaptation. Conclusions: Progress in understandi

Background: During copulation, the major Afro-tropical malaria vector Anopheles gambiae s.s. transfers male accessory gland (MAG) proteins to females as a solid mass (i.e. the “mating plug”). These proteins are postulated to function as important modulators of female post-mating responses. To understand the role of selective forces underlying the evolution of these proteins in the A. gambiae complex, we carried out an evolutionary analysis of gene sequence and expression divergence on a pair of paralog genes called AgAcp34A-1 and AgAcp34A-2. These encode MAG-specific proteins which, based on homology with Drosophila, have been hypothesized to play a role in sperm viability and function. Results: Genetic analysis of 6 species of the A. gambiae complex revealed the existence of a third paralog (68-78% of identity), that we named AgAcp34A-3. FISH assays showed that this gene maps in the same division (34A) of chromosome-3R as the other two paralogs. In particular, immuno-fluorescence assays targeting the C-terminals of AgAcp34A-2 and AgAcp34A-3 revealed that these two proteins are localized in the posterior part of the MAG and concentrated at the apical portion of the mating plug. When transferred to females, this part of the plug lies in proximity to the duct connecting the spermatheca to the uterus, suggesting a potential role for these proteins in regulating sperm motility. AgAcp34A-3 is more polymorphic than the other two paralogs, possibly because of relaxation of purifying selection. Since both unequal crossing-over and gene conversion likely homogenized the members of this gene family, the interpretation of the evolutionary patterns is not straightforward. Although several haplotypes of the three paralogs are shared by most A. gambiae s.l. species, some fixed species-specific replacements (mainly placed in the N- and C-terminal portions of the secreted peptides) were also observed, suggesting some lineage-specific adaptation. Conclusions: Progress in understandi

Hydraena is the largest genus within the water beetle family Hydraenidae, with ca. 1000 species distributed worldwide. Within this large genus some monophyletic groups of species are recognised, among them the “Haenydra” lineage, including ca. 90 species distributed in the western Palaearctic from the Iberian peninsula to Iran. Species of “Haenydra” have often very restricted distributions, and are typical of clean small rivers and streams. We obtained ca. 2.5 Kb of mitochondrial and nuclear protein-code and ribosomal markers of 101 specimens of 69 species of “Haenydra”, and used Bayesian and Maximum Likelihood phylogenetic methods to reconstruct their phylogeny and diversification history. We found a derived phylogenetic position of the “Haenydra” lineage within the genus Hydraena, as sister to the species of the Hydraenabisulcata group., Within “Haenydra” three main lineages were recognised, with poorly resolved relationships among them: the Hydraena iberica, Hydraena gracilis and Hydraena dentipes lineages, the former restricted to the Iberian peninsula but the latter two distributed through the whole north-Mediterranean area. A Bayesian relaxed molecular clock approach using a combined mitochondrial rate of 2% divergence per MY estimated the origin of “Haenydra” in the Tortonian, ca. 8 Mya, and the main diversification and the origin of most extant species in the Pliocene and Pleistocene. We did not found evidence of a phylogenetic connection between the western and eastern species that could be traced to the Messinian salinity crisis, with dispersal only at small geographical scales (e.g. the colonisation of Corsica and Sardinia from NW Italy and SW France). The H. gracilis and H. iberica lineages were estimated to have diversified under a pure birth model with a speciation rate of 0.64 and 0.23 species/MY respectively, while the H. dentipes lineage was estimated to have a decreasing diversification rate with time, with an average rate of 0.29 sp/MY.

Background SINEs (Short INterspersed Elements) are homoplasy-free and co-dominant genetic markers which are considered to represent useful tools for population genetic studies, and could help clarifying the speciation processes ongoing within the major malaria vector in Africa, Anopheles gambiae s.s. Here, we report the results of the analysis of the insertion polymorphism of a nearly 200 bp-long SINE (SINE200) within genome areas of high differentiation (i.e. "speciation islands") of M and S A. gambiae molecular forms. Methods A SINE-PCR approach was carried out on thirteen SINE200 insertions in M and S females collected along the whole range of distribution of A. gambiae s.s. in sub-Saharan Africa. Ten specimens each for Anopheles arabiensis, Anopheles melas, Anopheles quadriannulatus A and 15 M/S hybrids from laboratory crosses were also analysed. Results Eight loci were successfully amplified and were found to be specific for A. gambiae s.s.: 5 on 2L chromosome and one on X chromosome resulted monomorphic, while two loci positioned respectively on 2R (i.e. S200 2R12D) and X (i.e. S200 X6.1) chromosomes were found to be polymorphic. S200 2R12D was homozygote for the insertion in most S-form samples, while intermediate levels of polymorphism were shown in M-form, resulting in an overall high degree of genetic differentiation between molecular forms (Fst = 0.46 p < 0.001) and within M-form (Fst = 0.46 p < 0.001). The insertion of S200 X6.1 was found to be fixed in all M- and absent in all S-specimens. This led to develop a novel easy-to-use PCR approach to straightforwardly identify A. gambiae molecular forms. This novel approach allows to overcome the constraints associated with markers on the rDNA region commonly used for M and S identification. In fact, it is based on a single copy and irreversible SINE200 insertion and, thus, is not subjected to peculiar evolutionary patterns affecting rDNA markers, e.g. incomplete homogenization of the arrays through concerted

Background SINEs (Short INterspersed Elements) are homoplasy-free and co-dominant genetic markers which are considered to represent useful tools for population genetic studies, and could help clarifying the speciation processes ongoing within the major malaria vector in Africa, Anopheles gambiae s.s. Here, we report the results of the analysis of the insertion polymorphism of a nearly 200 bp-long SINE (SINE200) within genome areas of high differentiation (i.e. "speciation islands") of M and S A. gambiae molecular forms. Methods A SINE-PCR approach was carried out on thirteen SINE200 insertions in M and S females collected along the whole range of distribution of A. gambiae s.s. in sub-Saharan Africa. Ten specimens each for Anopheles arabiensis, Anopheles melas, Anopheles quadriannulatus A and 15 M/S hybrids from laboratory crosses were also analysed. Results Eight loci were successfully amplified and were found to be specific for A. gambiae s.s.: 5 on 2L chromosome and one on X chromosome resulted monomorphic, while two loci positioned respectively on 2R (i.e. S200 2R12D) and X (i.e. S200 X6.1) chromosomes were found to be polymorphic. S200 2R12D was homozygote for the insertion in most S-form samples, while intermediate levels of polymorphism were shown in M-form, resulting in an overall high degree of genetic differentiation between molecular forms (Fst = 0.46 p < 0.001) and within M-form (Fst = 0.46 p < 0.001). The insertion of S200 X6.1 was found to be fixed in all M- and absent in all S-specimens. This led to develop a novel easy-to-use PCR approach to straightforwardly identify A. gambiae molecular forms. This novel approach allows to overcome the constraints associated with markers on the rDNA region commonly used for M and S identification. In fact, it is based on a single copy and irreversible SINE200 insertion and, thus, is not subjected to peculiar evolutionary patterns affecting rDNA markers, e.g. incomplete homogenization of the arrays through concerted

Background SINEs (Short INterspersed Elements) are homoplasy-free and co-dominant genetic markers which are considered to represent useful tools for population genetic studies, and could help clarifying the speciation processes ongoing within the major malaria vector in Africa, Anopheles gambiae s.s. Here, we report the results of the analysis of the insertion polymorphism of a nearly 200 bp-long SINE (SINE200) within genome areas of high differentiation (i.e. "speciation islands") of M and S A. gambiae molecular forms. Methods A SINE-PCR approach was carried out on thirteen SINE200 insertions in M and S females collected along the whole range of distribution of A. gambiae s.s. in sub-Saharan Africa. Ten specimens each for Anopheles arabiensis, Anopheles melas, Anopheles quadriannulatus A and 15 M/S hybrids from laboratory crosses were also analysed. Results Eight loci were successfully amplified and were found to be specific for A. gambiae s.s.: 5 on 2L chromosome and one on X chromosome resulted monomorphic, while two loci positioned respectively on 2R (i.e. S200 2R12D) and X (i.e. S200 X6.1) chromosomes were found to be polymorphic. S200 2R12D was homozygote for the insertion in most S-form samples, while intermediate levels of polymorphism were shown in M-form, resulting in an overall high degree of genetic differentiation between molecular forms (Fst = 0.46 p < 0.001) and within M-form (Fst = 0.46 p < 0.001). The insertion of S200 X6.1 was found to be fixed in all M- and absent in all S-specimens. This led to develop a novel easy-to-use PCR approach to straightforwardly identify A. gambiae molecular forms. This novel approach allows to overcome the constraints associated with markers on the rDNA region commonly used for M and S identification. In fact, it is based on a single copy and irreversible SINE200 insertion and, thus, is not subjected to peculiar evolutionary patterns affecting rDNA markers, e.g. incomplete homogenization of the arrays through concerted

Background SINEs (Short INterspersed Elements) are homoplasy-free and co-dominant genetic markers which are considered to represent useful tools for population genetic studies, and could help clarifying the speciation processes ongoing within the major malaria vector in Africa, Anopheles gambiae s.s. Here, we report the results of the analysis of the insertion polymorphism of a nearly 200 bp-long SINE (SINE200) within genome areas of high differentiation (i.e. "speciation islands") of M and S A. gambiae molecular forms. Methods A SINE-PCR approach was carried out on thirteen SINE200 insertions in M and S females collected along the whole range of distribution of A. gambiae s.s. in sub-Saharan Africa. Ten specimens each for Anopheles arabiensis, Anopheles melas, Anopheles quadriannulatus A and 15 M/S hybrids from laboratory crosses were also analysed. Results Eight loci were successfully amplified and were found to be specific for A. gambiae s.s.: 5 on 2L chromosome and one on X chromosome resulted monomorphic, while two loci positioned respectively on 2R (i.e. S200 2R12D) and X (i.e. S200 X6.1) chromosomes were found to be polymorphic. S200 2R12D was homozygote for the insertion in most S-form samples, while intermediate levels of polymorphism were shown in M-form, resulting in an overall high degree of genetic differentiation between molecular forms (Fst = 0.46 p < 0.001) and within M-form (Fst = 0.46 p < 0.001). The insertion of S200 X6.1 was found to be fixed in all M- and absent in all S-specimens. This led to develop a novel easy-to-use PCR approach to straightforwardly identify A. gambiae molecular forms. This novel approach allows to overcome the constraints associated with markers on the rDNA region commonly used for M and S identification. In fact, it is based on a single copy and irreversible SINE200 insertion and, thus, is not subjected to peculiar evolutionary patterns affecting rDNA markers, e.g. incomplete homogenization of the arrays through concerted

Background SINEs (Short INterspersed Elements) are homoplasy-free and co-dominant genetic markers which are considered to represent useful tools for population genetic studies, and could help clarifying the speciation processes ongoing within the major malaria vector in Africa, Anopheles gambiae s.s. Here, we report the results of the analysis of the insertion polymorphism of a nearly 200 bp-long SINE (SINE200) within genome areas of high differentiation (i.e. "speciation islands") of M and S A. gambiae molecular forms. Methods A SINE-PCR approach was carried out on thirteen SINE200 insertions in M and S females collected along the whole range of distribution of A. gambiae s.s. in sub-Saharan Africa. Ten specimens each for Anopheles arabiensis, Anopheles melas, Anopheles quadriannulatus A and 15 M/S hybrids from laboratory crosses were also analysed. Results Eight loci were successfully amplified and were found to be specific for A. gambiae s.s.: 5 on 2L chromosome and one on X chromosome resulted monomorphic, while two loci positioned respectively on 2R (i.e. S200 2R12D) and X (i.e. S200 X6.1) chromosomes were found to be polymorphic. S200 2R12D was homozygote for the insertion in most S-form samples, while intermediate levels of polymorphism were shown in M-form, resulting in an overall high degree of genetic differentiation between molecular forms (Fst = 0.46 p < 0.001) and within M-form (Fst = 0.46 p < 0.001). The insertion of S200 X6.1 was found to be fixed in all M- and absent in all S-specimens. This led to develop a novel easy-to-use PCR approach to straightforwardly identify A. gambiae molecular forms. This novel approach allows to overcome the constraints associated with markers on the rDNA region commonly used for M and S identification. In fact, it is based on a single copy and irreversible SINE200 insertion and, thus, is not subjected to peculiar evolutionary patterns affecting rDNA markers, e.g. incomplete homogenization of the arrays through concerted