Dynamic of SARS-CoV-2 variants circulation in Tunisian pediatric population, during successive waves, from March 2020 to September 2022.

The emergence of SARS-CoV-2 variants has led to several cases among children. However, limited information is available from North African countries. This study describes the SARS-CoV-2 strains circulating in Tunisian pediatric population during successive waves. A total of 447 complete sequences were obtained from individuals aged from 13 days to 18 years, between March 2020 and September 2022: 369 sequences generated during this study and 78 ones, available in GISAID, previously obtained from Tunisian pediatric patients. These sequences were compared with 354 and 274 ones obtained from Tunisian adults and a global dataset, respectively. The variant circulation dynamics of predominant variants were investigated during the study period using maximum-likelihood phylogenetic analysis. Among the studied population, adolescents were the predominant age group, comprising 55.26% of cases. Twenty-three lineages were identified; seven of which were not previously reported in Tunisia. Phylogenetic analysis showed a close relationship between the sequences from Tunisian adults and children. The connections of sequences from other countries were variable according to variants: close relationships were observed for Alpha, B1.160 and Omicron variants, while independent Tunisian clusters were observed for Delta and B.1.177 lineages. These findings highlight the pivotal role of children in virus transmission and underscore the impact of vaccination on virus spread. Vaccination of children, with booster doses, may be considered for better management of future emergences.

Development and evaluation of an easy to use real-time reverse-transcription loop-mediated isothermal amplification assay for clinical diagnosis of West Nile virus.

West Nile Virus (WNV) causes a serious public health concern in many countries around the world. Virus detection in pathological samples is a key component of WNV infection diagnostic, classically performed by real-time PCR. In outbreak situation, rapid detection of the virus, in peripheral laboratories or at point of care, is crucial to guide decision makers and for the establishment of adequate action plans to prevent virus dissemination. Here, we evaluate a Loop-mediated isothermal amplification (LAMP) tool for WNV detection. Amplifications were performed comparatively on extracted viral RNA and on crude samples using a classical thermal cycler and a portable device (pebble device). qRT-PCR was used as gold standard and two sets of urine samples (n = 62 and n = 74) were used to evaluate the retained amplification protocols and assess their sensitivity and specificity. RT-LAMP on RNA extracts and crude samples showed a sensitivity of 90 % and 87 %, respectively. The specificity was 100 % for extracts and 97 % for crude samples. Using the device, the RT-LAMP on extracted RNA was comparable to the gold standard results (100 % sensitivity and specificity) and it was a bit lower on crude samples (65 % sensitivity and 94 % specificity). These results show that RT-LAMP is an efficient technique to detect WNV. RT-LAMP provides a rapid, sensitive, high-throughput and portable tool for accurate WNV detection and has potentials to facilitate diagnostic and surveillance efforts both in the laboratory and in the field, especially in developing countries.

Investigation of the Sandfly Fauna of Central Arid Areas and Northern Humid Regions of Tunisia, with Morphological and Molecular Identification of the Recently Established Population of Phlebotomus (Larroussius) perfiliewi.

Based on nucleotide sequences, we re-identified representative samples of Phlebotomus perfiliewi originating from two different biogeographical areas of Tunisia, whose populations had previously been identified based on morphological criteria. A partial region of the mitochondrial DNA cytochrome b gene was targeted, and sandfly species was determined by analogy with DNA sequences available in the GenBank database via a BLAST analysis, taking into account the query coverage and percentage identity. The recognized species presents the most substantial homology with the analyzed sequence. The results of the molecular identification showed complete agreement with the morphological identifications, and Phlebotomus perfiliewi is genetically variable.

The value of West Nile virus RNA detection by real-time RT-PCR in urine samples from patients with neuroinvasive forms.

INTRODUCTION: Routine laboratory screening is based on the detection of WNV specific IgM and IgG in blood and cerebrospinal fluid. Confirmation is then classically applied by real time RT-PCR (rRT-PCR) in Cerebrospinal fluid (CSF), which often gives negative results due to too short virorachia and late sampling. rRT-PCR was applied-for the first time for routine diagnosis purpose-on urine samples. METHODS: During 2018 outbreak in Tunisia, 107 patients presented WNV neurologic symptoms and were positive for WNV serology. Of them, 95 patients were sampled for urine and 35 were sampled for CSF. Qualitative rRT-PCR was performed on both type of samples. RESULTS: WNV RNA was detected in 50.5% of urine samples (48/95) and in 2.8% of CSF samples (1/35). WNV RNA was detectable from day 1 to day 41 from symptom onset, however, positive urine rate was 53.1% during the first 10 days from symptom onset. The proportions of urine-positive and urine-negative samples, based on day of collection, showed no statistical difference (p > 0.005). Cycle threshold (Ct) values ranged from 12 to 39, with no correlation with the day of collection. The lowest Ct value was detected for urine sampled on day 5 after symptom onset. A statistically significant difference was found between age groups of confirmed and non confirmed cases (p < 0.001). DISCUSSION/CONCLUSION: Our study reported the use of rRT-PCR on urine samples as a confirmatory diagnostic tool for WNV "probable cases" during an outbreak. Our findings underlined the reliability and the rapidity of this confirmatory tool, even late, and showed its superiority on CSF investigation.

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance.

Investment in Africa over the past year with regards to SARS-CoV-2 genotyping has led to a massive increase in the number of sequences, exceeding 100,000 genomes generated to track the pandemic on the continent. Our results show an increase in the number of African countries able to sequence within their own borders, coupled with a decrease in sequencing turnaround time. Findings from this genomic surveillance underscores the heterogeneous nature of the pandemic but we observe repeated dissemination of SARS-CoV-2 variants within the continent. Sustained investment for genomic surveillance in Africa is needed as the virus continues to evolve, particularly in the low vaccination landscape. These investments are very crucial for preparedness and response for future pathogen outbreaks. One-Sentence SummaryExpanding Africa SARS-CoV-2 sequencing capacity in a fast evolving pandemic.

Molecular Epidemiology of SARS-CoV-2 in Tunisia (North Africa) through Several Successive Waves of COVID-19.

Documenting the circulation dynamics of SARS-CoV-2 variants in different regions of the world is crucial for monitoring virus transmission worldwide and contributing to global efforts towards combating the pandemic. Tunisia has experienced several waves of COVID-19 with a significant number of infections and deaths. The present study provides genetic information on the different lineages of SARS-CoV-2 that circulated in Tunisia over 17 months. Lineages were assigned for 1359 samples using whole-genome sequencing, partial S gene sequencing and variant-specific real-time RT-PCR tests. Forty-eight different lineages of SARS-CoV-2 were identified, including variants of concern (VOCs), variants of interest (VOIs) and variants under monitoring (VUMs), particularly Alpha, Beta, Delta, A.27, Zeta and Eta. The first wave, limited to imported and import-related cases, was characterized by a small number of positive samples and lineages. During the second wave, a large number of lineages were detected; the third wave was marked by the predominance of the Alpha VOC, and the fourth wave was characterized by the predominance of the Delta VOC. This study adds new genomic data to the global context of COVID-19, particularly from the North African region, and highlights the importance of the timely molecular characterization of circulating strains.

SARS-CoV2 RT-PCR assays: In vitro comparison of 4 WHO approved protocols on clinical specimens and its implications for real laboratory practice through variant emergence.

INTRODUCTION: RT-PCR testing on nasopharyngeal swabs is a key component in the COVID-19 fighting, provided to use sensitive and specific SARS-CoV2 genome targets. In this study, we aimed to evaluate and to compare 4 widely used WHO approved RT-PCR protocols on real clinical specimens, to decrypt the reasons of the diverging results and to propose recommendations for the choice of the genome targets. METHODS: 260 nasopharyngeal samples were randomly selected among the samples tested between Week-16, 2020 and week-16 2021, in the Institut Pasteur de Tunis, Tunisia, one of the referent laboratories of COVID-19 in Tunisia. All samples were tested by Charité, Berlin protocol (singleplex envelop (E) and singleplex RNA-dependent RNA polymerase (RdRp)), Hong Kong Universiy, China protocol (singleplex nucleoprotein (N) and singleplex Open reading frame Orf1b), commercial test DAAN Gene® (using the CDC China protocol), (triplex N, Orf1ab with internal control) and Institut Pasteur Paris protocol (IPP) (triplex IP2(nsp9) and IP4 (nsp12) with internal control). For IPP, a selection from samples positive by IP2 but negative with IP4 was re-tested by exactly the same protocol but this time in singleplex. New results were described and analyzed. RESULTS: In vitro analysis showed discordant results in 29.2% of cases (76 out of 260). The most discordant protocol is DAAN Gene® due to the false positive late signals with N target. Discordant results between the two protocol's targets are more frequent when viral load are low (high Ct values). Our results demonstrated that the multiplexing has worsened the sensitivity of the IP4 target. CONCLUSION: We provide concise recommendations for the choice of the genome targets, the interpretation of the results and the alarm signals which makes suspect a gene mutation.

The Delta variant wave in Tunisia: Genetic diversity, spatio-temporal distribution and evidence of the spread of a divergent AY.122 sub-lineage.

Introduction: The Delta variant posed an increased risk to global public health and rapidly replaced the pre-existent variants worldwide. In this study, the genetic diversity and the spatio-temporal dynamics of 662 SARS-CoV2 genomes obtained during the Delta wave across Tunisia were investigated. Methods: Viral whole genome and partial S-segment sequencing was performed using Illumina and Sanger platforms, respectively and lineage assignemnt was assessed using Pangolin version 1.2.4 and scorpio version 3.4.X. Phylogenetic and phylogeographic analyses were achieved using IQ-Tree and Beast programs. Results: The age distribution of the infected cases showed a large peak between 25 to 50 years. Twelve Delta sub-lineages were detected nation-wide with AY.122 being the predominant variant representing 94.6% of sequences. AY.122 sequences were highly related and shared the amino-acid change ORF1a:A498V, the synonymous mutations 2746T>C, 3037C>T, 8986C>T, 11332A>G in ORF1a and 23683C>T in the S gene with respect to the Wuhan reference genome (NC_045512.2). Spatio-temporal analysis indicates that the larger cities of Nabeul, Tunis and Kairouan constituted epicenters for the AY.122 sub-lineage and subsequent dispersion to the rest of the country. Discussion: This study adds more knowledge about the Delta variant and sub-variants distribution worldwide by documenting genomic and epidemiological data from Tunisia, a North African region. Such results may be helpful to the understanding of future COVID-19 waves and variants.

Sequencing Using a Two-Step Strategy Reveals High Genetic Diversity in the S Gene of SARS-CoV-2 after a High-Transmission Period in Tunis, Tunisia.

Recent efforts have reported numerous variants that influence severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral characteristics, including pathogenicity, transmission rate, and detectability by molecular tests. Whole-genome sequencing based on next-generation sequencing technologies is the method of choice to identify all viral variants; however, the resources needed to use these techniques for a representative number of specimens remain limited in many low- and middle-income countries. To decrease sequencing costs, we developed a primer set allowing partial sequences to be generated in the viral S gene, enabling rapid detection of numerous variants of concern (VOCs) and variants of interest (VOIs); whole-genome sequencing is then performed on a selection of viruses based on partial sequencing results. Two hundred one nasopharyngeal specimens collected during the decreasing phase of a high-transmission COVID-19 wave in Tunisia were analyzed. The results reveal high genetic variability within the sequenced fragment and allow the detection of first introductions in the country of already-known VOCs and VOIs, as well as other variants that have interesting genomic mutations and need to be kept under surveillance. IMPORTANCE The method of choice for SARS-CoV-2 variant detection is whole-genome sequencing using next-generation sequencing (NGS) technologies. Resources for this technology remain limited in many low- and middle-income countries, where it is not possible to perform whole-genome sequencing for representative numbers of SARS-CoV-2-positive cases. In the present work, we developed a novel strategy based on a first partial Sanger screening in the S gene, which includes key mutations of the already known VOCs and VOIs, for rapid identification of these VOCs and VOIs and to help better select specimens that need to be sequenced by NGS technologies. The second step consists of whole-genome sequencing to allow a holistic view of all variants within the selected viral strains and confirm the initial classification of the strains based on partial S gene sequencing.

Risk Assessment of the Role of the Ecotones in the Transmission of Zoonotic Cutaneous Leishmaniasis in Central Tunisia.

Zoonotic cutaneous leishmaniasis (ZCL), endemic in Central and Southern Tunisia, is caused by Leishmania major (Kinetoplastida: Trypanosomatidae), which is transmitted by the sand fly Phlebotomus papatasi. In Tunisia, the fat sand rat Psammomys obesus and the desert jird Meriones shawi are the principal reservoir hosts of L. major. The presence of the P. papatasi vector of the L. major etiologic agent of ZCL was assessed in the vicinity of villages in endemic areas of Central Tunisia. The study was performed from September through October 2019, a period corresponding to the main peak of activity of P. papatasi. Sand flies were collected from rodent burrows located at the ecotone level, which is the transition zone between the natural environment and human settlement. Sand flies were identified to species level and tested for the presence of L. major by PCR. Our entomological survey showed that P. papatasi is the most abundant sand fly species associated with rodent burrows, and this abundance is even higher in ecotones primarily occupied by P. obesus in comparison to ecotones occupied by M. shawi. Infections with Leishmania major were detected only in P. papatasi, with an overall minimum infection rate (MIR) of 2.64%. No significant difference was observed between the MIRs in ecotones of P. obesus and of M. shawi. Incidence of ZCL in the studied areas ranged from 200 to 700 cases per 100,000 inhabitants, with a mean incidence of 385.41 per 100,000. Higher ZCL incidence was identified in ecotones of M. shawi compared to ecotones of P. obesus. ZCL cases are positively correlated with the MIRs. Considering the short flight range of P. papatasi, increases in its densities associated with burrows of P. obesus or M. shawi at the ecotone level expand the overlap of infected vectors with communities and subsequently increase ZCL incidence. Therefore, control measures should target P. papatasi populations at the ecotones.

Multiplexed Magnetofluorescent Bioplatform for the Sensitive Detection of SARS-CoV-2 Viral RNA without Nucleic Acid Amplification.

Rapid and sensitive detection of SARS-CoV-2 virus genetic material is of paramount importance to mitigate the COVID-19 pandemic outbreak and lower the death toll. Herein, we report the design of a magnetofluorescent bioplatform for the direct and specific detection of the viral RNA of SARS-CoV-2 in the total RNA extracted from nasopharyngeal swabs of COVID-19-positive patients. A higher fluorescence response was achieved using two capture probes tethered to magnetic beads using a biotin/streptavidin linkage, targeting two specific sites in the ORF1a and S genes. Two horseradish peroxidase (HRP)-conjugated reporter sequences, complementary to the loci of the S and N genes, were used to reveal the presence of the viral RNA through the oxidation of o-phenylenediamine to fluorescent 2,3-diaminophenazine. Under optimal conditions, the bioplatform showed high selectivity and sensitivity and was able to detect as low as 0.01 ng of viral RNA (1 × 103 copies/µL) with a linear dynamic range varying from 0.01 to 3.0 ng (1 × 103 to 9 × 107 copies/µL). The bioplatform was also able to discriminate the SARS-CoV-2 RNA from those of other related viruses such as hepatitis C, West Nile, measles, and non-polio viruses. Furthermore, the developed biosensor was validated in 46 clinical samples (36 COVID-19-positive patients and 10 COVID-19-negative subjects, as assessed with the gold standard RT-qPCR method). Both sensitivity and specificity of the developed method reached 100%. Finally, making such a simple and specific method available in the field, at a primary point of care, can better help the detection of SARS-CoV-2 infection in low-resource settings.

Whole genome sequencing and phylogenetic analysis of six SARS-CoV-2 strains isolated during COVID-19 pandemic in Tunisia, North Africa.

BACKGROUND: In Tunisia a first SARS-CoV-2 confirmed case was reported in March 03, 2020. Since then, an increase of cases number was observed from either imported or local cases. The aim of this preliminary study was to better understand the molecular epidemiology and genetic variability of SARS-CoV-2 viruses circulating in Tunisia and worldwide. METHODS: Whole genome sequencing was performed using NGS approach on six SARS. CoV-2 highly positive samples detected during the early phase of the outbreak. RESULTS: Full genomes sequences of six Tunisian SARS-CoV-2 strains were obtained from imported and locally transmission cases during the COVID-19 outbreak. Reported sequences were non-identical with 0.1% nucleotide divergence rate and clustered into 6 different clades with worldwide sequences. SNPs results favor the distribution of the reported Tunisian sequences into 3 major genotypes. These SNP mutations are critical for diagnosis and vaccine development. CONCLUSIONS: These results indicate multiple introductions of the virus in Tunisia and add new genomic data on SARS-CoV-2 at the international level.

A year of genomic surveillance reveals how the SARS-CoV-2 pandemic unfolded in Africa

The progression of the SARS-CoV-2 pandemic in Africa has so far been heterogeneous and the full impact is not yet well understood. Here, we describe the genomic epidemiology using a dataset of 8746 genomes from 33 African countries and two overseas territories. We show that the epidemics in most countries were initiated by importations, predominantly from Europe, which diminished following the early introduction of international travel restrictions. As the pandemic progressed, ongoing transmission in many countries and increasing mobility led to the emergence and spread within the continent of many variants of concern and interest, such as B.1.351, B.1.525, A.23.1 and C.1.1. Although distorted by low sampling numbers and blind-spots, the findings highlight that Africa must not be left behind in the global pandemic response, otherwise it could become a breeding ground for new variants.

Tick-borne encephalitis virus in Ixodes ricinus (Acari: Ixodidae) ticks, Tunisia.

Tick-borne encephalitis virus (TBEV) is the etiologic agent of tick-borne encephalitis (TBE) and transmitted by Ixodes ricinus (Linnaeus, 1758) in Europe. The geographical distribution of I. ricinus in the Palearctic region covers also northern Africa, including northwestern Tunisia. While the eco-epidemiology of TBE in Europe is well documented, no data concerning TBEV from northern Africa are available. We investigated whether TBEV is circulating in Tunisia. A total of 877 adult I. ricinus collected from northwestern Tunisia were examined in pools for the presence of TBEV by nRT-PCR. Viral RNA was detected in one pool of three engorged ticks, yielding a minimum infection rate of 0.11 % (1/877). Phylogenetic analysis showed that the Tunisian TBEV strain belongs to the European lineage. We report for the first time the presence of TBEV in I. ricinus from northern Africa. Therefore, more studies are needed to assess the public health importance of TBEV in northern Africa.

Genetic characterization of West Nile Virus strains during neuroinvasives infection outbreak in Tunisia, 2018.

West Nile Virus (WNV) is an arbovirus transmitted by mosquito bite involving birds as reservoirs, humans and equines as accidental hosts. Eight distinct lineages (WNV-1 to WNV-8) have been identified: WNV-1 and WNV-2 infect humans and animals, and WNV-3 to WNV-8 have been identified only in vectors. WNV has been implicated in neuroinvasives infections, especially meningitis and encephalitis. Tunisia experienced three epidemics in 1997, 2003 and 2012. Serological studies on humans, equines and birds as well as the detection of the virus in the vector favour a fairly frequent circulation in the country. A new epidemic has been observed in Tunisia between August and November 2018. The obtained sequences of the VWN from Tunisia 2018 grouped in a distinct monophyletic group within the Mediterranean subtype in Cluster 1, with a maximum of 2% nucleotide divergence. These sequences were clearly distinct from the Tunisia 1997, which grouped with sequences mainly from USA in Cluster 2. This work reports the genetic characterization of the Tunisia 2018 strain in comparison with the previously identified strains in Tunisia and worldwide. The epidemic virus Tunisia 2018 was genetically close to the Mediterranean basin and Eastern Europe sequences but distinct from the Tunisia 1997 closely related to the American sequences.

The Impact of Illegal Waste Sites on the Transmission of Zoonotic Cutaneous Leishmaniasis in Central Tunisia.

Illegal waste disposal represents a risk health factor for vector-borne diseases by providing shelter for rodents and their ectoparasites. The presence of the Phlebotomus papatasi vector of Leishmania major, an etiologic agent of zoonotic cutaneous leishmaniasis (ZCL), was assessed at illegal waste sites located at the vicinity of villages in endemic areas of Central Tunisia. The study was performed over a two-year period over three nights from July to September 2017, and over three nights in September 2018. Household waste is deposited illegally forming dumpsites at the vicinity of each village and contains several rodent burrows of Psammomys obesus, the main reservoir host of L. major. Sandflies were collected from rodent burrows in the natural environment and in dumpsites using sticky traps and were identified at species level. Female sandflies were tested for the presence of L. major by PCR. Our entomological survey showed that Phlebotomus papatasi is the most abundant sandfly species associated with rodent burrows in these waste sites. The densities of P. papatasi in dumpsites are significantly higher compared to the natural environment. The minimum infection rate of P. papatasi with L. major in these illegal waste sites is not significantly different compared to the natural environment. Considering the short flight range of P. papatasi, increases in its densities, associated with burrows of P. obesus in illegal waste sites located at the edge of villages, expands the overlap of infected ZCL vectors with communities. Thus, illegal waste sites pose a high risk of spreading ZCL to neighboring home ranges. Waste management is an environmentally friendly method of controlling sandfly populations and should be included in an integrated management program for controlling ZCL in endemic countries.

Ixodes inopinatus and Ixodes ricinus (Acari: Ixodidae) Are Sympatric Ticks in North Africa.

In the present study, we report the sympatric occurrence of Ixodes ricinus (Linnaeus, 1758) and Ixodes inopinatus (Estrada-Peña, Nava, and Petney, 2014) in Tunisia. In total, 173 adult Ixodes ticks were collected from four sites (El Jouza, Tamra, Aïn Soltan, and Jbel Zaghouan) between February and April 2017, a period corresponding to the peak of activity of I. ricinus in North Africa. The morphological characters corresponded to both species; thus, we generated a total of 28 16S rRNA sequences and compared them with previously published data in GenBank. The two species were sympatric in Tamra, Aïn Soltan, and El Jouza, whereas collections in Jbel Zaghouan only yielded I. inopinatus. These results indicate that the two taxa are widespread in the humid area of northern Tunisia. The one tick collected in Jbel Zaghouan suggests that the distribution of at least I. inopinatus might extend to the sub-humid area. More studies are needed to fully comprehend the systematic status of the two taxonomic entities using multiple molecular markers and morphological characters; integrating these two identification methods are a necessary step toward a better understanding of the ecology and epidemiology of tick-borne diseases in Tunisia.

Absence of Crimean-Congo haemorrhagic fever virus in the tick Hyalomma aegyptium parasitizing the spur-thighed tortoise (Testudo graeca) in Tunisia.

Free-ranging spur-thighed tortoises Testudo graeca, captured in different habitat types of Northern Tunisia from March to April 2017, were examined for tick infestation: 134/147 (91%) were infested. The overall infestation intensity and abundance was 8.5 and 7.8, respectively. From these tortoises, 1174 ticks were collected, of which 10% (n = 120) taken from 18 randomly-selected tortoises were identified at the species level; the remaining ticks were examined for the presence of Crimean-Congo haemorrhagic fever virus (CCHFv) by real time RT-PCR. Only adult Hyalomma aegyptium were found, suggesting a high degree of host specificity to tortoises. No CCHFv was detected in ticks. Considering the absence of CCHFv in Hyalomma aegyptium infesting its main host, the spur-thighed tortoise, this tick species is unlikely to play a major role in the epidemiology of CCHF. Therefore, more studies are needed to investigate the circulation of this arbovirus between livestock and other tick species from North Africa.

An integrated overview of the midgut bacterial flora composition of Phlebotomus perniciosus, a vector of zoonotic visceral leishmaniasis in the Western Mediterranean Basin.

BACKGROUND: The Leishmania developmental life cycle within its sand fly vector occurs exclusively in the lumen of the insect's digestive tract in the presence of symbiotic bacteria. The composition of the gut microbiota and the factors that influence its composition are currently poorly understood. A set of factors, including the host and its environment, may influence this composition. It has been demonstrated that the insect gut microbiota influences the development of several human pathogens, such as Plasmodium falciparum. For sand flies and Leishmania, understanding the interactions between the parasite and the microbial environment of the vector midgut can provide new tools to control Leishmania transmission. METHODOLOGY/PRINCIPAL FINDINGS: The midguts of female Phlebotomus perniciosus from laboratory colonies or from the field were collected during the months of July, September and October 2011 and dissected. The midguts were analyzed by culture-dependent and culture-independent methods. A total of 441 and 115 cultivable isolates were assigned to 30 and 11 phylotypes from field-collected and colonized P. perniciosus, respectively. Analysis of monthly variations in microbiota composition shows a species diversity decline in October, which is to the end of the Leishmania infantum transmission period. In parallel, a compilation and a meta-analysis of all available data concerning the microbiota of two Psychodidae genera, namely Phlebotomus and Lutzomyia, was performed and compared to P. perniciosus, data obtained herein. This integrated analysis did not reveal any substantial divergences between Old and New world sand flies with regards to the midgut bacterial phyla and genera diversity. But clearly, most bacterial species (>76%) are sparsely distributed between Phlebotominae species. CONCLUSION/SIGNIFICANCE: Our results pinpoint the need for a more exhaustive understanding of the bacterial richness and abundance at the species level in Phlebotominae sand flies in order to capture the role of midgut bacteria during Leishmania development and transmission. The occurrence of Bacillus subtilis in P. perniciosus and at least two other sand fly species studied so far suggests that this bacterial species is a potential candidate for paratransgenic or biolological approaches for the control of sand fly populations in order to prevent Leishmania transmission.

Changes of Sand Fly Populations and Leishmania infantum Infection Rates in an Irrigated Village Located in Arid Central Tunisia.

The current spread of zoonotic visceral leishmaniasis (ZVL) throughout arid areas of Central Tunisia is a major public health concern. The main objective of this study is to investigate whether the development of irrigation in arid bio-geographical areas in Central Tunisia have led to the establishment of a stable cycle involving sand flies of the subgenus Larroussius and Leishmania infantum, and subsequently to the emergence of ZVL. Sand flies were collected from the village of Saddaguia, a highly irrigated zone located within an arid bio-geographical area of Central Tunisia by using modified Centers for Diseases Control (CDC) light traps. Morphological keys were used to identify sand flies. Collected sand flies were pooled with up to 30 specimens per pool according to date and tested by nested Polymerase Chain Reaction (PCR) DNA sequencing from positive pools was used to identify Leishmania spp. A total of 4915 sand flies (2422 females and 2493 males) were collected from Saddaguia in September and in October 2014. Morphological identification confirmed sand flies of the subgenus Larroussius to be predominant. PCR analysis followed by DNA sequencing indicated that 15 pools were infected with L. infantum yielding an overall infection rate of 0.6%. The majority of the infected pools were of sand fly species belonging to subgenus Larroussius. Intense irrigation applied to the arid bio-geographical areas in Central Tunisia is at the origin of the development of an environment capable of sustaining important populations of sand flies of the subgenus Larroussius. This has led to the establishment of stable transmission cycles of L. infantum and subsequently to the emergence of ZVL.