Filling the gaps in the global prevalence map of clinical antimicrobial resistance.

Surveillance is critical in containing globally increasing antimicrobial resistance (AMR). Affordable methodologies to prioritize AMR surveillance efforts are urgently needed, especially in low- and middle-income countries (LMICs), where resources are limited. While socioeconomic characteristics correlate with clinical AMR prevalence, this correlation has not yet been used to estimate AMR prevalence in countries lacking surveillance. We captured the statistical relationship between AMR prevalence and socioeconomic characteristics in a suite of beta-binomial principal component regression models for nine pathogens resistant to 19 (classes of) antibiotics. Prevalence data from ResistanceMap were combined with socioeconomic profiles constructed from 5,595 World Bank indicators. Cross-validated models were used to estimate clinical AMR prevalence and temporal trends for countries lacking data. Our approach provides robust estimates of clinical AMR prevalence in LMICs for most priority pathogens (cross-validated q2 > 0.78 for six out of nine pathogens). By supplementing surveillance data, 87% of all countries worldwide, which represent 99% of the global population, are now informed. Depending on priority pathogen, our estimates benefit 2.1 to 4.9 billion people living in countries with currently insufficient diagnostic capacity. By estimating AMR prevalence worldwide, our approach allows for a data-driven prioritization of surveillance efforts. For carbapenem-resistant Acinetobacter baumannii and third-generation cephalosporin-resistant Escherichia coli, specific countries of interest are located in the Middle East, based on the magnitude of estimates; sub-Saharan Africa, based on the relative prevalence increase over 1998 to 2017; and the Pacific Islands, based on improving overall model coverage and performance.

Uncovering Active Bacterial Symbionts in Three Species of Pollen-feeding Beetles (Nitidulidae: Meligethinae).

Microbial symbionts enable many phytophagous insects to specialize on plant-based diets through a range of metabolic services. Pollen comprises one-plant tissue consumed by such herbivores. While rich in lipids and proteins, its nutrient content is often imbalanced and difficult-to-access due to a digestibly recalcitrant cell wall. Pollen quality can be further degraded by harmful allelochemicals. To identify microbes that may aid in palynivory, we performed cDNA-based 16S rRNA metabarcoding on three related pollen beetles (Nitidulidae: Meligethinae) exhibiting different dietary breadths: Brassicogethes aeneus, B. matronalis, and Meligethes atratus. Nine bacterial symbionts (i.e., 97% OTUs) exhibited high metabolic activity during active feeding. Subsequent PCR surveys revealed varying prevalence of those from three Rickettsialles genera-Lariskella, Rickettsia, and Wolbachia-within beetle populations. Our findings lay the groundwork for future studies on the influence of phylogeny and diet on palynivorous insect microbiomes, and roles of symbionts in the use of challenging diets.

New Evidence of Canthariphily: Tilloidea transversalis (Coleoptera: Cleridae) Sequestering Cantharidin From Lydus trimaculatus (Coleoptera: Meloidae).

Cantharidin (CTD) is a defensive compound autogenously and exclusively produced by two phylogenetically related beetle families: Meloidae and Oedemeridae. Although this molecule usually acts as a strong deterrent against potential predators and parasites, some arthropod species, collectively named 'canthariphilous species', are attracted to CTD. Some species can sequester CTD from the CTD-producing species, using it as a chemical defense against enemies. The present paper focuses on the first-ever description of canthariphilous interactions between a checkered beetle species (Coleoptera: Cleridae) and a CTD -producing species. Field observations revealed individuals of the phytophagous beetle Tilloidea transversalis (Charpentier, 1825) (Coleoptera: Cleridae) biting individuals of the blister beetle Lydus trimaculatus (Fabricius, 1775) (Coleoptera: Meloidae). Laboratory behavioral experiments followed to verify if this peculiar behavior of T. transversalis also occurs on other co-occurring species. Moreover, chemical analyses were performed to assess whether T. transversalis can sequester CTD. Our results show that T. transversalis only attacks CTD-producing species. However, while chemical analyses prove that T. transversalis can sequester CTD from the hemolymph of L. trimaculatus, some clues (based on a CTD-baited traps sampling) suggest that this beetle, contrarily to other canthariphilous species, does not appear to show a high attraction to pure synthetic CTD. Thus, other unknown signals, alone or in combination with CTD, could be implicated in triggering the canthariphilous behaviors of T. transversalis.

New Approach to Privacy-Preserving Clinical Decision Support Systems for HIV Treatment.

BACKGROUND: HIV treatment prescription is a complex process. Clinical decision support systems (CDSS) are a category of health information technologies that can assist clinicians to choose optimal treatments based on clinical trials and expert knowledge. The usability of some CDSSs for HIV treatment would be significantly improved by using the knowledge obtained by treating other patients. This knowledge, however, 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 HIV 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. FINDINGS: Our solution enables to compute an effectiveness measure of an HIV treatment, the average time-to-treatment-failure, while preserving privacy. Experimental results show that our solution, although at proof-of-concept stage, has good efficiency and provides a result to a query within 24 min for a dataset of realistic size. INTERPRETATION: This paper presents a novel and efficient approach HIV clinical decision support systems, that harnesses the potential and insights acquired from treatment data, while preserving the privacy of patient records and the confidentiality of clinician decisions.

Unraveling the role of male reproductive tract and haemolymph in cantharidin-exuding Lydus trimaculatus and Mylabris variabilis (Coleoptera: Meloidae): a comparative transcriptomics approach.

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.

Rosaceae, Brassicaceae and pollen beetles: exploring relationships and evolution in an anthophilous beetle lineage (Nitidulidae, Meligethes-complex of genera) using an integrative approach.

BACKGROUND: Meligethes are pollen-beetles associated with flowers of Rosaceae as larvae. This genus currently consists of 63 known species in two subgenera, Meligethes and Odonthogethes, predominantly occurring in the eastern Palaearctic. We analyzed 74 morphological and ecological characters (169 states) of all species, as well as of 11 outgroup species from 7 Meligethinae genera (including Brassicogethes), to investigate their phylogeny. We also conducted a parallel molecular analysis on 9 Meligethes, 9 Odonthogethes, 3 Brassicogethes and 2 Meligethinus species based on DNA sequence data from mitochondrial (COI, 16S) and nuclear (CAD) genes. RESULTS: Morphological phylogenetic reconstructions supported the monophyly of the whole genus and clades corresponding to purported subgenera Meligethes s.str. and Odonthogethes. Main species-groups were mostly confirmed, however some unresolved polytomies remained. Molecular data placed members of Brassicogethes (including 42 mostly W Palearctic species associated with Brassicaceae) as sister to Odonthogethes, with this clade being sister to Meligethes s.str. This phylogenetic scenario suggests that monophyletic Meligethes s.str., Odonthogethes and Brassicogethes should be regarded alternatively as three subgenera of a monophyletic Meligethes, or three genera in a monophyletic genus-complex, with mutually monophyletic Brassicogethes and Odonthogethes. Molecular analyses estimated the origin of this lineage at ca. 14-15 Mya from a common stem including Meligethinus. CONCLUSIONS: We hypothesize that the ancestor of Meligethes specialized on Rosaceae in the Middle Miocene (likely in Langhian Age) and subsequently radiated during Late Miocene and Plio-Pleistocene maintaining a trophic niche on this plant family. This radiation was primarily due to geographic isolation in E Asiatic mountain systems. Combined evidence from morphology, ancestral state parsimony reconstruction of host-plant associations and molecular evidence suggested that Rosoideae (Rosa spp.) represented the ancestral hosts of Meligethes s.str., followed by an independent shift of ancestral Odonthogethes (ca. 9-15 Mya) on Rubus (Rosoideae) and members of Rosaceae Spiraeoideae. Other ancestral Odonthogethes probably shifted again on the unrelated plant family Brassicaceae (maybe 8-14 Mya in S China), allowing a rapid westward radiation of the Brassicogethes clade.

Phylogeny, biogeography and systematics of the hyper-diverse blister beetle genus Hycleus (Coleoptera: Meloidae).

Hycleus is a hyper-diverse genus of blister beetles including ~500 species widely distributed in the Old World, currently divided into three "sections" and into 45 "phenetic" species groups according to morphological characters. Recently the monophyly of Hycleus was questioned pointing out its paraphyly with respect to the genera Ceroctis and Paractenodia. In this study, we built a time-calibrated phylogenetic tree based on DNA sequence data from mitochondrial and nuclear genes obtained from 125 species, to understand the phylogenetic relationships among the species of this genus, to infer the biogeographic processes behind their diversification, and to assess their taxonomy and classification. Our results identified four main lineages one of which included the species belonging to Ceroctis and Paractenodia; therefore, both taxa are now referred to Hycleus as new synonyms. The three described sections of Hycleus resulted polyphyletic and are rejected, whereas several species groups represented well supported clades. Hycleus likely originated in Africa during the Early Miocene (~20 Mya), and subsequently spread in Europe and western Asia. Later, in the Late Miocene (~6 Mya) a Saharo-Sindian group branched off from the Palaearctic lineage, whereas the Oriental Region was colonized following a dispersal event through the Arabian Peninsula from the Afrotropical Region (~5 Mya).

Convection and the Extracellular Matrix Dictate Inter- and Intra-Biofilm Quorum Sensing Communication in Environmental Systems.

The mechanisms and impact of bacterial quorum sensing (QS) for the coordination of population-level behaviors are well studied under laboratory conditions. However, it is unclear how, in otherwise open environmental systems, QS signals accumulate to sufficient concentration to induce QS phenotypes, especially when quorum quenching (QQ) organisms are also present. We explore the impact of QQ activity on QS signaling in spatially organized biofilms in scenarios that mimic open systems of natural and engineered environments. Using a functionally differentiated biofilm system, we show that the extracellular matrix, local flow, and QQ interact to modulate communication. In still aqueous environments, convection facilitates signal dispersal while the matrix absorbs and relays signals to the cells. This process facilitates inter-biofilm communication even at low extracellular signal concentrations. Within the biofilm, the matrix further regulates the transport of the competing QS and QQ molecules, leading to heterogenous QS behavior. Importantly, only extracellular QQ enzymes can effectively control QS signaling, suggesting that the intracellular QQ enzymes may not have evolved to degrade environmental QS signals for competition.

Game of neutrophils: modeling the balance between apoptosis and necrosis.

BACKGROUND: Neutrophils are one of the key players in the human innate immune system (HIIS). In the event of an insult where the body is exposed to inflammation triggering moieties (ITMs), neutrophils are mobilized towards the site of insult and antagonize the inflammation. If the inflammation is cleared, neutrophils go into a programmed death called apoptosis. However, if the insult is intense or persistent, neutrophils take on a violent death pathway called necrosis, which involves the rupture of their cytoplasmic content into the surrounding tissue that causes local tissue damage, thus further aggravating inflammation. This seemingly paradoxical phenomenon fuels the inflammatory process by triggering the recruitment of additional neutrophils to the site of inflammation, aimed to contribute to the complete neutralization of severe inflammation. This delicate balance between the cost and benefit of the neutrophils' choice of death pathway has been optimized during the evolution of the innate immune system. The goal of our work is to understand how the tradeoff between the cost and benefit of the different death pathways of neutrophils, in response to various levels of insults, has been optimized over evolutionary time by using the concepts of evolutionary game theory. RESULTS: We show that by using evolutionary game theory, we are able to formulate a game that predicts the percentage of necrosis and apoptosis when exposed to various levels of insults. CONCLUSION: By adopting an evolutionary perspective, we identify the driving mechanisms leading to the delicate balance between apoptosis and necrosis in neutrophils' cell death in response to different insults. Using our simple model, we verify that indeed, the global cost of remaining ITMs is the driving mechanism that reproduces the percentage of necrosis and apoptosis observed in data and neutrophils need sufficient information of the overall inflammation to be able to pick a death pathway that presumably increases the survival of the organism.

A multilocus view on Mediterranean aeolid nudibranchs (Mollusca): Systematics and cryptic diversity of Flabellinidae and Piseinotecidae.

Recent molecular studies revealed high level of endemism and numerous cryptic species within opisthobranchs, with Mediterranean taxa clearly understudied. Here we used genetic data from both mitochondrial and nuclear gene fragments as well as morphological data from taxonomically relevant characters to investigate the phylogenetic relationships and systematics of Mediterranean taxa of the Flabellinidae and Piseinotecidae families. Phylogenetic analyses based on Bayesian and Maximum-Likelihood methods indicate that Flabellinidae and Pisenotecidae taxa and species within the genera Flabellina, Calmella and Piseinotecus do not form monophyletic clades. These results are supported by our morphological analyses which allowed the re-evaluation of the triseriate radula condition in Pisenotecidae and Calmella taxa and their inclusion in the genus Flabellina as Flabellina gaditanacomb. nov. (synonym of F. confusa), Flabellina gabiniereicomb. nov. and Flabellina cavolinicomb. nov. Species delimitation and barcoding gap analyses allowed uncovering cryptic species within Flabellina gracilis (Alder and Hancock, 1844), F. trophina (Bergh, 1890), F. verrucosa (M. Sars, 1829) and F. ischitana Hirano and Thompson, 1990, the latter with an Atlantic form which is under description. This study corroborates the relevance of combining molecular and morphological data from multiple populations and species in the assessment of nudibranch diversity and classification.

The last bastion? X chromosome genotyping of Anopheles gambiae species pair males from a hybrid zone reveals complex recombination within the major candidate 'genomic island of speciation'.

Speciation with gene flow may be aided by reduced recombination helping to build linkage between genes involved in the early stages of reproductive isolation. Reduced recombination on chromosome X has been implicated in speciation within the Anopheles gambiae complex, species of which represent the major Afrotropical malaria vectors. The most recently diverged, morphologically indistinguishable, species pair, A. gambiae and Anopheles coluzzii, ubiquitously displays a 'genomic island of divergence' spanning over 4 Mb from chromosome X centromere, which represents a particularly promising candidate region for reproductive isolation genes, in addition to containing the diagnostic markers used to distinguish the species. Very low recombination makes the island intractable for experimental recombination studies, but an extreme hybrid zone in Guinea Bissau offers the opportunity for natural investigation of X-island recombination. SNP analysis of chromosome X hemizygous males revealed: (i) strong divergence in the X-island despite a lack of autosomal divergence; (ii) individuals with multiple-recombinant genotypes, including likely double crossovers and localized gene conversion; (iii) recombination-driven discontinuity both within and between the molecular species markers, suggesting that the utility of the diagnostics is undermined under high hybridization. The largely, but incompletely protected nature of the X centromeric genomic island is consistent with a primary candidate area for accumulation of adaptive variants driving speciation with gene flow, while permitting some selective shuffling and removal of genetic variation.

Remarkable diversity of intron-1 of the para voltage-gated sodium channel gene in an Anopheles gambiae/Anopheles coluzzii hybrid zone.

BACKGROUND: Genomic differentiation between Anopheles gambiae and Anopheles coluzzii--the major malaria vectors in sub-Saharan Africa--is localized into large "islands" toward the centromeres of chromosome-X and the two autosomes. Linkage disequilibrium between these genomic islands was first detected between species-specific polymorphisms within ribosomal DNA genes (IGS-rDNA) on the X-chromosome and a single variant at position 702 of intron 1 (Int-1702) of the para Voltage-Gated Sodium Channel (VGSC) gene on chromosome arm 2 L. Intron-1 sequence data from West and Central Africa revealed two clearly distinct and species-specific haplogroups, each characterized by very low polymorphism, which has been attributed to a selective sweep. The aim of this study was to analyse Int-1 sequence diversity in A. gambiae and A. coluzzii populations from the Far-West of their range, in order to assess whether this selective-sweep signature could persist in a zone of high interspecific hybridization. METHODS: A 531 bp region of VGSC Int-1 was sequenced in 21 A. coluzzii, 31 A. gambiae, and 12 hybrids from The Gambia and Guinea Bissau, located within the Far-West geographical region, and in 53 A. gambiae s.l. samples from the rest of the range. RESULTS: Far-West samples exhibit dramatic Int-1 polymorphism, far higher within each country than observed throughout the rest of the species range. Moreover, patterning of haplotypes within A. coluzzii confirms previous evidence of a macro-geographic subdivision into a West and a Central African genetic cluster, and reveals a possible genetic distinction of A. coluzzii populations from the Far-West. CONCLUSIONS: The results suggest a relaxation of selective pressures acting across the VGSC gene region in the hybrid zone. Genetic differentiation in the Far-West could be attributable to a founder effect within A. coluzzii, with subsequent extensive gene flow with secondarily-colonizing A. gambiae, potentially yielding a novel insight on the dynamic processes impacting genetic divergence of these key malaria vectors.

Comparative genomics provides insights into molecular adaptation to hypermetamorphosis and cantharidin metabolism in blister beetles (Coleoptera: Meloidae).

Blister beetles (Coleoptera: Meloidae) are currently subdivided into three subfamilies: Eleticinae (a basal group), Nemognathinae, and Meloinae. These are all characterized by the endogenous production of the defensive terpene cantharidin (CA), whereas the two most derived subfamilies show a hypermetamorphic larval development. Here, we provide novel draft genome assemblies of five species sampled across the three blister beetle subfamilies (Iselma pallidipennis, Stenodera caucasica, Zonitis immaculata, Lydus trimaculatus, and Mylabris variabilis) and performed a comparative analysis with other available Meloidae genomes and the closely-related canthariphilous species (Pyrochroa serraticornis) to disclose adaptations at a molecular level. Our results highlighted the expansion and selection of genes potentially responsible for CA production and metabolism, as well as its mobilization and vesicular compartmentalization. Furthermore, we observed adaptive selection patterns and gain of genes devoted to epigenetic regulation, development, and morphogenesis, possibly related to hypermetamorphosis. We hypothesize that most genetic adaptations occurred to support both CA biosynthesis and hypermetamorphosis, two crucial aspects of Meloidae biology that likely contributed to their evolutionary success.

The cranial apparatus glands of the canthariphilous Pyrochroa coccinea (Coleoptera: Pyrochroidae: Pyrochroinae), and their implications in sexual behaviour.

Some Pyrochroidae species are known as "canthariphilous" for their attraction to cantharidin (CTD), a toxic terpene with anti-predatory effects, produced in nature by only two beetle families (Meloidae and Oedemeridae). It has been demonstrated that males of Neopyrochroa flabellata ingesting CTD are positively selected by females. Indeed, the compound is re-emitted from a glandular cranial apparatus as secretions that are licked up by females during courtship behaviour, inducing copulation. Herein, we provide the first description of the glands associated to the cranial apparatus of male Pyrochroinae using the European species Pyrochroa coccinea as a model. Morphological analyses show that the cranial apparatus consists of a concave pit lined with short setae retaining secretions emitted through numerous glandular pores. Ultrastructural investigations reveal the presence of two different class 3 glands (Gl.A and Gl.B), intermixed at the level of the pit but exhibiting distinct features. Gl.A are mainly characterised by short conducting canals, rounded nuclei and electrondense vesicles while Gl.B are characterised by long conducting canals, irregular nuclei, vesicles containing a particulate substance and a multifolded plasma membrane. Observations of sexual behaviour are also reported for P. coccinea and compared to N. flabellata, confirming the involvement of cranial apparatus secretions in courtship behaviour.

A multiscale hybrid model for exploring the effect of Resolvin D1 on macrophage polarization during acute inflammation.

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.

Male Accessory Glands of Blister Beetles and Cantharidin Release: A Comparative Ultrastructural Analysis.

Members of the family Meloidae are known to produce cantharidin, a highly toxic monoterpene found in their hemolymph and exuded as droplets capable of deterring many predators. As a nuptial gift, males transfer large amounts of cantharidin to females via a spermatophore, which is formed by specific accessory glands containing high concentrations of this terpene. Using light, electron and ion beam microscopy, the ultrastructural features of the three pairs of male accessory glands as well as the glandular part of the vasa deferentia were comparatively investigated in seven species of blister beetles belonging to five different tribes and two subfamilies. All gland pairs examined share common features such as mesodermal derivation, the presence of muscle sheath, a developed rough endoplasmic reticulum, abundant mitochondria, secretory vesicles, and microvillated apical membranes. Within the same species, glands exhibit distinctive features, suggesting that each pair is responsible for the formation of a specific substance. The vasa deferentia, while showing many similarities within the family, often exhibit features unique to each of the individual species investigated, whereas the accessory glands of the first and second pairs display the highest degree of ultrastructural variability. A comparison across the species shows an interesting constancy limited to ultrastructural features in the third pair of accessory glands. The similarities and differences among the species are discussed in the light of the available literature and in relation to the potential role that blister beetles' male accessory glands could play in the storage and management of cantharidin.

Molecular Characterization of Kunitz-Type Protease Inhibitors from Blister Beetles (Coleoptera, Meloidae).

Protease inhibitors are widely studied since the unrestricted activity of proteases can cause extensive organ lesions. In particular, elastase activity is involved in the pathophysiology of acute lung injury, for example during SARS-CoV-2 infection, while serine proteases and thrombin-like proteases are involved in the development and/or pathology of the nervous system. Natural protease inhibitors have the advantage to be reversible and with few side effects and thus are increasingly considered as new drugs. Kunitz-type protease inhibitors (KTPIs), reported in the venom of various organisms, such as wasps, spiders, scorpions, and snakes, have been studied for their potent anticoagulant activity and widespread protease inhibitor activity. Putative KTPI anticoagulants have been identified in transcriptomic resources obtained for two blister beetle species, Lydus trimaculatus and Mylabris variabilis. The KTPIs of L. trimaculatus and M. variabilis were characterized by combined transcriptomic and bioinformatics methodologies. The full-length mRNA sequences were divided on the base of the sequence of the active sites of the putative proteins. In silico protein structure analyses of each group of translational products show the biochemical features of the active sites and the potential protease targets. Validation of these genes is the first step for considering these molecules as new drugs for use in medicine.

Molecular evolution of a gene cluster of serine proteases expressed in the Anopheles gambiae female reproductive tract.

BACKGROUND: Genes involved in post-mating processes of multiple mating organisms are known to evolve rapidly due to coevolution driven by sexual conflict among male-female interacting proteins. In the malaria mosquito Anopheles gambiae - a monandrous species in which sexual conflict is expected to be absent or minimal - recent data strongly suggest that proteolytic enzymes specifically expressed in the female lower reproductive tissues are involved in the processing of male products transferred to females during mating. In order to better understand the role of selective forces underlying the evolution of proteins involved in post-mating responses, we analysed a cluster of genes encoding for three serine proteases that are down-regulated after mating, two of which specifically expressed in the atrium and one in the spermatheca of A. gambiae females. RESULTS: The analysis of polymorphisms and divergence of these female-expressed proteases in closely related species of the A. gambiae complex revealed a high level of replacement polymorphisms consistent with relaxed evolutionary constraints of duplicated genes, allowing to rapidly fix novel replacements to perform new or more specific functions. Adaptive evolution was detected in several codons of the 3 genes and hints of episodic selection were also found. In addition, the structural modelling of these proteases highlighted some important differences in their substrate specificity, and provided evidence that a number of sites evolving under selective pressures lie relatively close to the catalytic triad and/or on the edge of the specificity pocket, known to be involved in substrate recognition or binding. The observed patterns suggest that these proteases may interact with factors transferred by males during mating (e.g. substrates, inhibitors or pathogens) and that they may have differently evolved in independent A. gambiae lineages. CONCLUSIONS: Our results - also examined in light of constraints in the application of selection-inference methods to the closely related species of the A. gambiae complex - reveal an unexpectedly intricate evolutionary scenario. Further experimental analyses are needed to investigate the biological functions of these genes in order to better interpret their molecular evolution and to assess whether they represent possible targets for limiting the fertility of Anopheles mosquitoes in malaria vector control strategies.

Molecular characterization and evolution of a gene family encoding male-specific reproductive proteins in the African malaria vector Anopheles gambiae.

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 understanding the signaling cascade in the A. gambiae reproductive pathway will elucidate the interaction of this MAG-specific protein family with their female counterparts. This knowledge will allow a better evaluation of the relative importance of genes involved in the reproductive isolation and fertility of A. gambiae species and could help the interpretation of the observed evolutionary patterns.

Molecular phylogeny and diversification of the "Haenydra" lineage (Hydraenidae, genus Hydraena), a north-Mediterranean endemic-rich group of rheophilic Coleoptera.

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.5Kb 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. 8Mya, 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.