Boris A E S Savassi

Gauthier G Dobigny

Jonas R Etougbétché

Thalasse T Avocegan

François T Quinsou

Philippe Gauthier

Moudachirou Ibikounlé

Hélène Moné

Gabriel Mouahid
email: mouahid@univ-perp.fr

Sonsino

x Schistosoma bovis Sonsino, 1876 introgressive hybrids

Keywords: Schistosoma haematobium x Schistosoma bovis, Cercarial chronobiology, Mastomys natalensis, Rattus rattus, Schistosome transmission

établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

Schistosomiasis is an acute and chronic disease caused by flatworms of the genus Schistosoma and estimates have shown that at least 229 million people required preventive treatment in 2018 (WHO 2020). The lifecycle of schistosomes is complex and includes a vertebrate definitive host, where the adult schistosomes pair and reproduce, and a gastropod snail intermediate host where the larvae multiply asexually. Transmission from the vertebrate host to the snail host is ensured by the miracidium larva. Transmission from the snail host to the vertebrate host is ensured by the cercaria larva which actively penetrates the vertebrate host. The timing of Schistosoma cercarial shedding from the intermediate snail hosts is synchronized by exogenous factors (the major synchronizer is the photoperiod) and is highly correlated to the time interval when the definitive host goes to water [START_REF] Combes | Behaviors in trematode cercariae that enhance parasite transmission: patterns and processes[END_REF]). As such, cercarial emission of schistosomes is a primordial parameter that determines successful transmission to the definitive mammalian host. As examples, the cercariae of Schistosoma bovis Sonsino, 1876 are shed during the early morning when cattle go to watering [START_REF] Mouahid | Schistosoma bovis: patterns of cercarial emergence from snails of the genera Bulinus and Planorbarius[END_REF][START_REF] Mouahid | Analyse comparative du rythme d'émission des cercaires de trois souches de Schistosoma bovis[END_REF][START_REF] Mouahid | Cercarial shedding patterns of Schistosoma bovis and Schistosoma haematobium from single and mixed infections of Bulinus truncatus[END_REF]; the cercariae of Schistosoma haematobium (Bilharz, 1852) Weinland, 1858 andS. mansoni Sambon, 1907 are shed around mid-day when humans go to water for bathing, washing, or fishing [START_REF] Mouahid | Cercarial shedding patterns of Schistosoma bovis and Schistosoma haematobium from single and mixed infections of Bulinus truncatus[END_REF][START_REF] Mouahid | A new chronotype of Schistosoma mansoni: adaptive significance[END_REF]; finally, the cercariae of S. rodhaini Brumpt 1931 and of some chronotypes of S. mansoni and of S. japonicum Katsurada 1904 are shed during the early night, when nocturnal rodents actively use water [START_REF] Théron | Hybrids between Schistosoma mansoni and S. rodhaini: characterization by cercarial emergence rhythms[END_REF][START_REF] Mouahid | A new chronotype of Schistosoma mansoni: adaptive significance[END_REF][START_REF] Su | A circular analysis of chronobiology of Schistosoma japonicum cercarial emergence from hilly areas of Anhui, China[END_REF]. This makes cercarial pattern a good marker for identifying which definitive host is responsible for Schistosoma transmission, mainly in the introgressive hybridization situations between different species of schistosomes that have increasingly been described in the literature over the last decade [START_REF] Boissier | Schistosoma spp[END_REF]. [START_REF] Savassi | Cattle as natural host for Schistosoma haematobium (Bilharz, 1852) Weinland, 1858 x Schistosoma bovis Sonsino, 1876 interactions, with new cercarial emergence and genetic patterns[END_REF] showed, for the first time, that cattle could be responsible for the transmission of S. haematobium x S. bovis populations in Benin, West Africa. The authors found atypical cercarial chronotypes: (i) early shedding that may correspond to S. bovis, (ii) mid-day shedding that may correspond to S. haematobium, and (iii) nocturnal cercarial emission unknown for the two species. The latter nocturnal chronobiological profile for S. haematobium x S. bovis raises suspicion that nocturnal mammals, such as rodents, could play a role in the transmission. The goal of this study is to test the hypothesis that micromammals play a role in schistosome transmission.

Materials and methods

Sites of collection

Small mammal sampling was performed in December 2019 in seven villages of southern Benin (Fig. 1): four villages within the Sô River hydrosystem (Ahomey Gblon, Gbessou, Sô-Ava Centre, and Vekky, all belonging to the Atlantic Department), and three along the Ouémé River (Kessounou, Avagbodji Bembè, and Avagbodji Diékomey, all belonging to the Ouémé Department). Two different areas of Kessounou were sampled: Kessounou Glo and Kessounou Bahoué located on the west and east banks of the Ouémé River, respectively.

Rodent trapping and identification

Since trapping success of different rodent species may depend on the types of traps used [START_REF] Garba | Spatial segregation between invasive and native commensal rodents in an urban environment: a case study in Niamey, Niger[END_REF], both Sherman traps (B. B. Sherman Traps, Inc., Tallahassee, FL, USA) and locally made wire-mesh live traps were systematically used together for a total of 498 night traps. Each trap was baited with the same mixture of peanut butter and fish and then set within houses over two successive nights. Captured animals were kept alive until their processing, which occurred within a week from the time of collection.

Rodents and shrews were identified using morphology, weight, and external measurements (head and body length, tail length, hindfoot length, and ear maximum length) according to standard taxonomy keys [START_REF] Granjon | Les rongeurs de l'Afrique sahélosoudanienne[END_REF]. However, coexistence of sibling species within African rodent genera is frequent. This is the reason why species-specific identification of all Schistosoma-carrying individuals was assessed following a complete cytochrome b mitochondrial gene sequencing as previously described in [START_REF] Dobigny | An update on the taxonomy and geographic distribution of the cryptic species Mastomys kollmannspergeri (Muridae, Murinae) using combined cytogenetic and molecular data[END_REF].

Adult schistosome recovery

The small mammals were euthanized by a lethal intraperitoneal injection of sodium pentobarbital. They were perfused [START_REF] Duvall | An improved perfusion technique for recovering adult schistosomes from laboratory animals[END_REF] using Masterflex L/S apparatus (Cole-Parmer Instrument Co., Cambridgeshire, UK), and adult schistosomes were collected. The trapped worms in the liver and mesenteric (Fig. 2) or vesical veins were collected after excision of these organs under dissecting microscope. Male and female worms were numbered and placed in 95° alcohol for molecular analyses.

Schistosoma egg morphology and miracidial hatching

Liver and duodenum of each rodent harboring Schistosoma adults were ground in saline solution (NaCl 0.9%) and filtered through a series of sieves (425-, 180-, 106-, and 45-μm, successively). The eggs were trapped on a 45-μm sieve. They were randomly sampled for both morphological and molecular studies and mounted individually in a 0.9% NaCl solution beneath glass cover slips on glass slides using a Pasteur pipette. They were then photographed and those which contained a living miracidium were measured (length, width, and spine length) under an optical microscope. In total, 24 schistosome eggs from Mastomys natalensis (Smith, 1834) and two eggs from Rattus rattus (Linnaeus, 1758) were measured. For molecular biology analyses, miracidial hatching was conducted by placing some eggs in a beaker containing well water and then several individual miracidia were plotted onto FTA® classic cards (Whatman, GE Healthcare companies, Little Chalfont, UK). The rest of the miracidia from the rodents were used for snail infections.

Snail infection

Seventy-two Bulinus globosus (Morelet, 1866) (6.4 ± 1 mm height) were individually exposed to a single miracidium. The presence of only one miracidium was verified for each snail. This was made because we wanted that all the cercariae coming from each monomiracidially infected snail had the same genotype (one miracidium gives rise, by asexual multiplication, to thousands of cercariae of only one genotype). This ensured us that all the cercariae that were counted for the cercarial emission pattern analyses were from the same genotype and that multi-peaks were coming without any doubt from only one miracidium. The snails were then fed with dry and fresh lettuce and maintained individually in a glass containing 40 ml of well water at 26 °C with a balanced photoperiod (12-h light, 12-h dark) and a photophase from 6:00 am to 6:00 pm. Once a day, each snail was transferred to a new glass containing well water, and the cercarial emission was checked.

Cercarial emergence pattern

The cercarial emergence pattern was analyzed every hour. For this analysis, each infected snail was transferred exactly every hour to a new glass containing 40 ml of well water at 26 °C; each transfer lasted approximately 3 s per snail. Regarding the specific transfers at 6 am and 6 pm, snails were transferred 1 min before light shifts. The water samples containing the emerged cercariae were filtered through polyamide nitrile filters (25-μm mesh sieve). The cercariae retained on the filter were stained with Lugol's solution and counted under a dissecting microscope.

The cercarial emission pattern was analyzed every hour over 6 consecutive days, i.e., over a period of 144 h (6 × 24). Several fresh cercariae were individually stored onto FTA® classic cards for further molecular investigations.

Schistosoma miracidium recoveries from urine and stool of schoolchildren

In order to get comparable data from the same period (December 2019) on schistosomiasis in humans, urine samples from 50 primary schoolchildren volunteers of Kessounou village (11 girls and 16 boys in the 4th grade; and 9 girls and 14 boys in the 5th grade) were collected individually and then mixed and passed through a 45-μm pore size sieve. The retained eggs were placed in a beaker containing well water for miracidial hatching. The stools from 46 of the primary schoolchildren were collected individually then mixed in a 0.9% NaCl saline solution. The mixture was passed through a series of sieves (425-, 180-, 106-, and 45-μm pore size, successively) and washed through with another 0.9% NaCl saline solution using a pressure pump. The eggs that were retained on the 45-μm pore size sieve were collected in a beaker containing well water for miracidial hatching. Using such a protocol, the infection status of each individual urine and feces sample could not be analyzed but all the schoolchildren were treated. Several miracidia retrieved from human urine and stool samples were collected individually and placed onto FTA® classic cards for further molecular analyses following the protocols detailed above.

Molecular analyses DNA extraction

Total DNA was extracted from each adult worm using the E.Z.N.A.® Tissue DNA Kit (Omega Bio-Tek, Norcross, GA, USA), as recommended by the manufacturer, and stored at -20 °C. DNA from miracidia and cercariae that were individually stored on FTA® classic cards were extracted as follows. 3-mm2 disks were cut and removed with a Craft Punch from the FTA® classic cards where the miracidium or cercaria was loaded. The disks were then deposited in Eppendorf tubes (1.5 ml) and an initial washing step with 100 μl of Milli-Q water was performed. After a 10-min incubation at room temperature, water was removed and replaced by 80 μl of 5% Chelex® 100 Molecular Biology Grade Resin solution (Bio-Rad Laboratories, Hercules, CA, USA). The samples were then heated to 65 °C for 30min at a stirring speed of 800 rpm, followed by a second heating phase at 99 °C for 8 min without stirring. DNA was finally collected by centrifugation at 14,000 rpm for 2 min and 50 μl of the supernatant (DNA) of each sample were stored at -20 °C for further molecular analyses.

Rapid diagnostic multiplex polymerase chain reaction

A rapid diagnostic multiplex polymerase chain reaction (PCR) targeting the cytochrome C oxidase subunit I mitochondrial DNA gene (COI mtDNA) was conducted as described in [START_REF] Savassi | Cattle as natural host for Schistosoma haematobium (Bilharz, 1852) Weinland, 1858 x Schistosoma bovis Sonsino, 1876 interactions, with new cercarial emergence and genetic patterns[END_REF]. The primers used were one universal reverse primer (Shmb.R: 5′-CAA GTA TCA TGA AAY ART ATR TCT AA-3′) and three species-specific forward primers: Sh.F: 5′-GGT CTC GTG TAT GAG ATC CTA TAG TTT G-3′ for S. haematobium (120 bp); Sb.F: 5′-GTT TAG GTA GTG TAG TTT GGG CTC AC-3′ for S. bovis (260 bp); and Sm.F: 5′-CTT TGA TTC GTT AAC TGG AGT G -3′ for S. mansoni (215 bp).

COI mtDNA PCR

Partial COI mtDNA amplification was performed by PCR using the Cox1_schist F forward primer (5′-TCTTTRGATCATAA GCG-3′), and the Cox1_schist-R reverse primer (5′-TAAT GCATMGGAAAAAAACA-3′) [START_REF] Lockyer | The phylogeny of the Schistosomatidae based on three genes with emphasis on the interrelationships of Schistosoma Weinland, 1858[END_REF]. PCR conditions used were similar to those described in [START_REF] Moné | Introgressive hybridizations of Schistosoma haematobium by Schistosoma bovis at the origin of the first case report of schistosomiasis in Corsica (France, Europe)[END_REF].

Internal transcribed spacer partial 18S, internal transcribed spacer 1, 5.8S, and internal transcribed spacer 2 ribosomal DNA PCR

PCR amplification of internal transcribed spacer (ITS) rDNA (partial 18S, internal transcribed spacer 1 (ITS1), 5.8S, and internal transcribed spacer 2 (ITS2)) was performed using the primers of [START_REF] Barber | PCR-RFLP analysis of the ITS2 region to identify Schistosoma haematobium and S. bovis from Kenya[END_REF] (ITS4F: 5′-TCCT CCGCTTATTGATATGC-3′, and ITS5R: 5′-GGAA GTAAAAGTCGTAACAAGG-3′) as well as PCR conditions described in [START_REF] Savassi | Cattle as natural host for Schistosoma haematobium (Bilharz, 1852) Weinland, 1858 x Schistosoma bovis Sonsino, 1876 interactions, with new cercarial emergence and genetic patterns[END_REF].

Sequencing

Partial COI mtDNA gene and ITS rDNA region were sequenced (Genoscreen; Lille, France) using the reverse primers; some regions were also sequenced using the forward primer in order to confirm the sequence. Sequences were all manually edited and verified using Sequencher 4.5 (Gene Codes Corporation, Ann Arbor, USA). The nucleotidic polymorphism of the nuclear sequences (ITS rDNA) was systematically investigated through the visualization of the raw chromatograms: since this gene has biparental transmission and segregation sites (polymorphic sites) between S. haematobium and S. bovis, the sequence chromatograms were carefully checked to identify the presence of possible heterozygous sites. At each polymorphic site (i.e., where two chromatogram peaks unambiguously overlapped, thus indicating the genetic signatures of both parent species), the IUPAC ambiguity symbols were used to indicate the individual nucleotide polymorphisms: thus, Y indicates the observed presence of both bases T and C, and not an ambiguous reading between T and C. Similarly, R indicates the observed presence of both bases A and G while W corresponds to the observed presence of both bases A and T.

Partial COI mtDNA phylogenetic analysis

DNA multiple sequence alignments were performed using Muscle program [START_REF] Edgar | MUSCLE: multiple sequence alignment with high accuracy and high throughput[END_REF]) implemented in MEGA 7.0 software [START_REF] Kumar | MEGA7:Molecular evolutionary genetics analysis Version 7.0[END_REF], and then refined using Gblocks 0.91b [START_REF] Castresana | Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis[END_REF][START_REF] Dereeper | Phylogeny.fr: robust phylogenetic analysis for the non-specialist[END_REF][START_REF] Dereeper | BLAST-EXPLORER helps you building datasets for phylogenetic analysis[END_REF]. The probabilistic model of sequence evolution [START_REF] Nei | Molecular evolution and phylogenetics[END_REF] and the gamma distribution (G) to approximate rate heterogeneity among haplotypes was applied using MEGA 7.0 software. The best model showing the lowest BIC score (Bayesian information criterion) was HKY + G (0.16) (Hasegawa-Kishino-Yano with the gamma distribution). A phylogenetic tree using the maximum likelihood method was constructed using this model of molecular evolution under MEGA 7.0 software. Robustness of internal nodes was systematically assessed using bootstrapping procedures (1000 replicates). The topology of the tree was rooted using S. intercalatum Fisher, 1934 as a closely related but external outgroup.

Statistical analyses

Means and standard errors were calculated. The Mann-Whitney test was used for morphological measure comparisons, and the Fisher's exact test was implemented for the comparisons of proportions using BiostaTGV (https://biostatgv.sentiweb.fr). A p value less than 0.05 (<0.05) was considered to be statistically significant.

Results

Trapping success and schistosome prevalence in small Mammals

In total, 62 small mammals were captured: 13/195 (6.7%) in Sherman traps and 49/303 (16.2%) in locally made wire-mesh traps (Fisher Exact Test, p=0.0031). The global trapping success was 12.5%, with local success ranging from 1.8% (Avagbodji Bembè) to 21.1% (Sô-Ava Centre) (Table 1). The small mammals included 61 rodents and 1 shrew (Table 1): 52 Rattus rattus, 3 Rattus norvegicus (Berkenhout, 1769), 6 Mastomys natalensis, and 1 Crocidura cf. olivieri (Lesson, 1827).

Among the 62 animals trapped, 50 individuals could be perfused and investigated for Schistosoma adults and eggs: 37 R. rattus (22 females and 15 males); 3 R. norvegicus (all females); 5 M. natalensis (3 females and 2 males) and 1 C. olivieri (female). Details by locality and species are provided in Table 1. In total, six rodents were found carrying Schistosoma adults: four male R. rattus from Vekky (V42_Rr_GDOB2229, V20_Rr_GDOB2231, V132_Rr_GDOB2232, and V17_Rr_GDOB2233) and two female M. natalensis from Kessounou Glo (K14_Mn _GDOB2246 and K34_Mn_GDOB2248). When considering only the individuals that could be perfused, this makes a species-specific prevalence of 9.8% (4/41) and 40% (2/5) in R. rattus and M. natalensis, respectively. Local prevalence reached 80% (4/5) for R. rattus from Vekky, and 50% (2/4) for M. natalensis from Kessounou Glo (Table 1).

Schistosome burden and egg presence

Schistosome burden (number of adult schistosomes recovered from the rodents) was variable between the rodents. In the four positive R. rattus from Vekky, we found 4 worms (2 females and 2 males) and only two eggs in V42_Rr (Field code=GDOB2229), 8 male worms in V20_Rr (GDOB2231), 3 worms (2 females, and 1 male) in V132_Rr (GDOB2232) and 1 male worm in V17_Rr (GDOB2233). The two eggs were used for morphology and all the adult worms were used for molecular studies.

In the two M. natalensis from Kessounou Glo, we found 30 worms (16 males and 7 pairs) in K14_Mn (GDOB2246) and 3 males in K34_Mn (GDOB2248). All of the adult worms were used for molecular studies. Several eggs could be extracted from K14_Mn and were used for both morphology and molecular studies.

Schistosoma egg morphology

Two morphotypes were observed in the eggs isolated from M. natalensis: S. haematobium morphotype and an intermediate morphotype between S. haematobium and S. bovis. The percentages of each morphotype were not significantly different (41.7 and 58.3%, respectively; Fisher exact test, p>0.05). Regarding egg morphometry, mean egg length was significantly lower in the S. haematobium morphotype (164.00 μm ±3.89; N=10) compared to the intermediate morphotype group (177.14 μm±3.88; N=14; Mann-Whitney test, p=0.029); the same trend was observed for the mean length/width ratio with 1.93±0.03 and 2.22±0.04, respectively (Mann-Whitney test, p=0.0001). The two eggs from R. rattus displayed intermediate morphotypes, the egg lengths were 165 and 240 μm and the mean length/width ratio 2.12 and 2.34, respectively.

Snail infection

Among the 72 B. globosus that were exposed to one miracidium per snail, 52 survived and 3 shed cercariae (5.8%) (Bg1, Bg2, and Bg3) (the miracidia were all extracted from the Mastomys natalensis K14_Mn (GDOB2246) from Kessounou Glo). Prepatent periods lasted 39 days from Bg1, 45 days from Bg2, and 46 days from Bg3.

Cercarial emission patterns

The cercarial emission patterns were analyzed for Bg1 and Bg2 over 6 consecutive days, from the 39th to the 44th day after infection for Bg1 and from the 49th to the 54th day after infection for the Bg2. The analysis could not be done for the Bg3 which died. The results showed that the two snails harbored different patterns.

The Bg1 snail had a late diurnal and nocturnal pattern, with 72% of the cercariae that were shed during the day. Cercarial emission occurred mainly from 2 to 7 pm, with a high plateau of emission between 1 and 6 pm and an emission peak focalized at 7 pm (Fig. 3a).

The Bg2 snail differed by the simultaneous presence of three peaks: an early diurnal peak at 7 am, a late diurnal peak at 4 pm, and a nocturnal peak at 7 pm. The percentage of cercariae shed during the day was nearly 53%. The nocturnal cercarial emission occurred as a very sharp peak at 7 pm, the first hour of darkness, representing 47% of the total cercarial daily emission (Fig. 3b).

COI mtDNA (1027 bp)

One hundred and fifty-six sequences were analyzed and showed 30 mutation sites and 14 haplotypes (Table 2). From the 96 sequences of schistosomes extracted from the two M. natalensis, we found only 1 haplotype (K14M1_Mn_Hap1) shared by the 33 adult worms from the two M. natalensis (30 from K14_Mn and 3 from K34_Mn), by the 34 miracidia obtained from K14_Mn, and by the 29 cercariae emitted by the 3 snails infected with one miracidium (13, 14, and 2, respectively).

Among the 13 sequences of schistosomes that were obtained from the four R. rattus, four haplotypes were identified: V20M1_Rr_Hap2, V20M2_Rr_Hap3, V20M3_Rr_Hap4, and V132M1_Rr_Hap5, all different from the haplotype K14M1_Mn_Hap1 found from M. natalensis. One to three of these Rattus-originating haplotypes were found within the same R. rattus individual.

From the 23 sequences of schistosomes retrieved fromchildren's urine, we found 7 haplotypes (KU1_Hs_Hap6; KU3_Hs_Hap7; KU4_Hs_Hap8; KU13_Hs_Hap9; KU5_Hs_Hap10; KU10_Hs_Hap11; KU11_Hs_Hap12) all different from those extracted from rodents. From the 24 sequences of schistosomes extracted from children stool, we found 8 haplotypes (KU1_Hs_Hap6; KU3_Hs_Hap7; KU4_Hs_Hap8; KU13_Hs_Hap9; KU10_Hs_Hap11, KS5_Hs_Hap13, KS10_Hs_Hap14), among which five were shared with some haplotypes found in the urine, and one haplotype (K14M1_Mn_Hap1) shared with the schistosome haplotypes found in M. natalensis. Sequence data were deposited in the NCBI GenBank database under the accession numbers from MW022134 to MW022147.

COI mtDNA phylogeny (884 bp)

The 14 haplotypes obtained from the COI mtDNA sequences (MW022134 toMW022147) were compared to 18 nucleotide sequences of Schistosoma with terminal-spined eggs [START_REF] Rollinson | The biology of schistosomes. From genes to latrines[END_REF] available in GenBank: 8 sequences of S. haematobium, from Gambia (JQ397349), Mali (AY157209), Egypt (JQ397368), Toho-Todougba_Benin (KT354661) , Melen_Gabon (KT354660) , Kenya (JQ397378), South Africa (JQ397397), and Madagascar (JQ397399); 3 sequences of S. bovis, from Senegal (AJ519521), Kenya (FJ897160), and Tanzania (AY157212); 4 sequences of S. bovis x S. haematobium, from Corsica, France (KT354656, KT354657, KT354658), and Sô-Tchanhoué, Benin (KT354662); 1 sequence of S. guineensis Pagès, Jourdane, Southgate & Tchuem Tchenté, 2003 from Sao Tome and Principe (AJ519517); 1 sequence of S. curassoni Brumpt, 1931 from Senegal (AJ519516); and 1 sequence of S. intercalatum from Democratic Republic of the Congo (AJ519515). Maximum likelihood tree topology showed that all the 14 Beninese haplotypes, which came either from rodents or schoolchildren, belonged to the S. bovis and S. bovis x S. haematobium clade, but not to the S. haematobium or to other Schistosoma species clades (Fig. 4). Genetic differences ranged between 0.11 and 1.13% among the haplotypes obtained from rodents; they ranged between 0 and 1.36% when both rodent-and human-originating haplotypes were considered. The genetic difference between schistosome haplotypes obtained from Benin rodents and the S. bovis haplotypes from Senegal, Kenya, and Tanzania ranged between 0.23 and 1.47%; this difference ranged between 0.23 and 1.36% between the rodentoriginating haplotypes and the previously published S. bovis x S. haematobium haplotypes from Benin and Corsica. The five haplotypes recovered from rodents were similar to haplotypes from cows, and one of them (K14M1_Mn_Hap1) was identical to a haplotype from humans.

ITS rDNA (946 bp)

One hundred and sixty-one 18S, ITS1, 5.8S, and ITS2 sequences were analyzed. The first 22 bp of our ITS rDNA sequences belonged to the 18S gene and no variability was found between any of our samples. We thus no longer gave this gene further consideration. Eleven profiles were found with the ITS1, 5.8S, and ITS2 sequences, with 1 to 3 profiles per rodent, only 1 profile for all the miracidia, only one profile for all the cercariae, and 6 and 7 profiles for urine and stool of schoolchildren, respectively (Table 1). Sequence data were deposited in the NCBI GenBank database under the accession numbers from MW027648 to MW027658.

Table 3 shows that the 11 profiles differ by their ITS1 gene (458 bp: from position 23 to position 480), with 5 genetic variants at four positions (TGAT, TGAY, TAAT, TART, and YAAT); by their 5.8S gene (155 bp: from position 481 to position 635), with 2 genetic variants at 1 position (C and Y); and by their ITS2 gene (311 bp: from position 636 to position 946), with 5 genetic variants at 6 positions (AATGAT, WATRAT, AATRAT, AGCGGC, and ARYGRY). Among these 11 profiles, four were found only in R. rattus (profiles 1 to 4: accession numbersMW027648 to MW027651), one profile corresponded to the typical S. bovis ITS rDNA (profile 1: accession number MW027648) and three profiles showed variants of the S. bovis typical profile with double peaks either in the ITS1 or in the ITS2 gene (profiles 2 to 4: accession numbers MW027649 to MW027651). Two profiles were found in both M. natalensis and schoolchildren (profiles 5 and 6: accession numbers MW027652 and MW027653), with profile 5 harboring both a typical S. haematobium ITS1 gene together with a typical S. bovis ITS2 gene, and profile 6 harboring a new combination which included the ITS1 gene (TART, variant of S. haematobium) together with a typical S. bovis ITS2 gene. Profile 6 was observed in all the miracidia and all the cercariae collected from the Bg1, Bg2, and Bg3 snails experimentally infected with one miracidium extracted from naturally infected M. natalensis. The other five profiles were found only in schoolchildren (profiles 7 to 11: accession numbers MW027654 to MW027658) and showed variants of the S. haematobium typical profile, with double peaks in the chromatograms of at least one of the ITS1, 5.8S, or ITS2 genes.

Discussion

The rodent species assemblage identified here in the commensal compartment of the Lake Nokoué ecosystem is quite similar to those observed in the urban and peri-urban areas of south Benin, including the lacustrine city of Ganvié [START_REF] Hima | Native and invasive small mammals in urban habitats along the commercial axis connecting Benin and Niger, West Africa[END_REF][START_REF] Houéménou | Pathogenic Leptospira in commensal small mammals from the extensively urbanized coastal Benin[END_REF]. Such a small mammal community has already been associated with the circulation of various zoonotic agents in this area: Bartonella [START_REF] Leulmi | Detection of Rickettsia felis, Rickettsia typhi, Bartonella species and Yersinia pestis in fleas (Siphonaptera) from Africa[END_REF][START_REF] Houéménou | Bartonella spp. in smallmammals, Benin[END_REF], waterassociated Leptospira [START_REF] Houéménou | Pathogenic Leptospira in commensal small mammals from the extensively urbanized coastal Benin[END_REF], Rickettsia [START_REF] Leulmi | Detection of Rickettsia felis, Rickettsia typhi, Bartonella species and Yersinia pestis in fleas (Siphonaptera) from Africa[END_REF]Trypanosoma (Dobigny et al. 2019). The results obtained in the present manuscript showed that two species of rodents were found to harbor schistosomes: Rattus rattus and Mastomys natalensis. Although sampling size was quite low, the schistosome host-specific prevalence observed here was quite high: 80% in R. rattus from Vekky and 50% in M. natalensis from Kessounou Glo. For R. rattus, the four COI schistosome haplotypes and the four ITS profiles belong to S. bovis species. To our knowledge, this is the first time that R. rattus has been found to be naturally infected by S. bovis. This is surprising for two reasons: (i) S. bovis is known to have a large vertebrate definitive host spectrum, including rodent species, but not R. rattus (Mastomys natalensis and Lophuromys flavopunctatus Thomas 1888 by [START_REF] Pitchford | A check list of definitive hosts exhibiting evidence of the genus SchistosomaWeinland, 1858 acquired naturally in Africa and the Middle East[END_REF] and by [START_REF] Hanelt | Schistosomes of small mammals from the Lake Victoria Basin, Kenya: new species, familiar species, and implications for schistosomiasis control[END_REF]; Arvicanthis niloticus (Desmarest, 1822) by [START_REF] Catalano | Rodents as natural hosts of zoonotic Schistosoma species and hybrids: an epidemiological and evolutionary perspective from West Africa[END_REF]); (ii) R. rattus is widespread in the lacustrine and semi-lacustrine habitats of south Benin [START_REF] Hima | Native and invasive small mammals in urban habitats along the commercial axis connecting Benin and Niger, West Africa[END_REF][START_REF] Houéménou | Pathogenic Leptospira in commensal small mammals from the extensively urbanized coastal Benin[END_REF] and is known to act as a natural host for other species of schistosomes such as S. japonicum in the Philippines [START_REF] Fedorko | Schistosoma japonicum in the black rat, Rattus rattus mindanensis, from Leyte, Philippines in relation to Oncomelania snail colonies with reference to other endoparasites[END_REF] and S. mansoni in the West Indies [START_REF] Combes | Le rôle du Rat dans la dynamique de l'endémie schistosomienne en Guadeloupe[END_REF]) and in Oman [START_REF] Mouahid | A new chronotype of Schistosoma mansoni: adaptive significance[END_REF]. Limited by the fact that only two dead eggs were found, we were unable to prove that R. rattus plays an effective role in the transmission of S. bovis. Further work should be done on the role of R. rattus in the transmission of schistosomes in the researched zone. For M. natalensis, the molecular analysis of the different stages of the collected schistosomes showed a COI of the S. bovis type and an ITS of the S. haematobium x S. bovis type. Previously, in Senegal, [START_REF] Catalano | Rodents as natural hosts of zoonotic Schistosoma species and hybrids: an epidemiological and evolutionary perspective from West Africa[END_REF] found one schistosome pair composed of a S. mansoni male with a hybrid female presenting the profile of S. haematobium for ITS and S. bovis for COI, in Mastomys huberti Wroughton, 1909; however, they did not provide any evidence of the presence of eggs or miracidia that could infect snails and shed cercariae. Our results make thus M. natalensis a natural host of an introgressive hybrid population between S. haematobium and S. bovis in the studied area. We provide the first example of a rodent being involved in the transmission of introgressive hybrid populations between S. haematobium and S. bovis. Indeed, the miracidia of these introgressive hybrids succeeded to infect snails and to give rise to cercariae. Such introgressive hybrid populations were also found in humans (this paper) and in bovines [START_REF] Savassi | Cattle as natural host for Schistosoma haematobium (Bilharz, 1852) Weinland, 1858 x Schistosoma bovis Sonsino, 1876 interactions, with new cercarial emergence and genetic patterns[END_REF] in the same geographical area; in humans, neither pure S. haematobium nor S. bovis species were found, while in bovines, S. bovis was present. Together, the three hosts, humans, rodents, and bovines, were found to share the same habitat as well as an introgressive hybrid schistosome population. The absence of pure S. haematobium may be surprising; however, the intensity of the interactions between the two species, S. haematobium and S. bovis, is very high, particularly in the studied area where bovines, humans, and rodents are livingwith high proximity (we collected the rodents in the human stilt houses) and can explain in this area that all the schistosomes that we collected were introgressive hybrids between the two species. Nothing excludes the effective presence of S. haematobium in other parts of Benin but this has to be investigated.

Two different Mastomys-related chronotypes were found and came fromtwo different snails, each infected with a single miracidium. The first chronotype, obtained from the Bg1 snail, exhibited a late diurnal and nocturnal pattern (Fig. 3a). We recently identified a similar pattern for schistosomes collected from cows in the very same area [START_REF] Savassi | Cattle as natural host for Schistosoma haematobium (Bilharz, 1852) Weinland, 1858 x Schistosoma bovis Sonsino, 1876 interactions, with new cercarial emergence and genetic patterns[END_REF]): it was the first time that a nocturnal shedding pattern was observed for schistosomes from cattle and led us to question the role of rodents as nocturnal animals. The second chronotype, obtained from the Bg2 snail (Fig. 3b), corresponds to a new chronotype with three peaks: an early diurnal, a late diurnal, and a nocturnal pattern. This chronotype cumulates all of the time periods associated with cercarial shedding patterns known to date: S. bovis typical pattern (early morning), S. haematobium typical pattern (diurnal or late diurnal), and a yet undescribed nocturnal pattern within S. haematobium and S. bovis species. The nocturnal period of cercarial emergence observed in both chronotypes can be linked to the nocturnal aquatic activity of M. natalensis. [START_REF] Duplantier | Rythmes d'activité chez six espèces de Muridés du Sénégal appartenant aux genres Mastomys, Arvicanthis[END_REF] showed via experimentation that the activity pattern of M. natalensis is primarily nocturnal. M. natalensis is the most widespread and common rodent in Africa south of the Sahara; it has a semicommensal to commensal habitat, where it plays primordial role in the spread of zoonotic diseases and is dependent on the presence of water [START_REF] Coetzee | The biology, behaviour, and ecology of Mastomys natalensis in southern Africa[END_REF][START_REF] Christensen | Home range and abundance of Mastomys natalensis (Smith, 1834) in habitats affected by cultivation[END_REF][START_REF] Granjon | Les rongeurs de l'Afrique sahélosoudanienne[END_REF]. Nocturnal cercarial emergences have already been observed in other species of schistosomes, namely, S. rodhaini [START_REF] Théron | Hybrids between Schistosoma mansoni and S. rodhaini: characterization by cercarial emergence rhythms[END_REF], S. mansoni [START_REF] Mouahid | A new chronotype of Schistosoma mansoni: adaptive significance[END_REF], and S. japonicum [START_REF] Su | A circular analysis of chronobiology of Schistosoma japonicum cercarial emergence from hilly areas of Anhui, China[END_REF], and were all related to the role of rodent hosts in the transmission. The fact that both of the Mastomysderived chronotypes, the one with both late diurnal and nocturnal patterns and the one with the three peaks (early diurnal, late diurnal, and nocturnal), came from snails infected with only onemiracidiumputs into question the circadian nature of all of the previously published results on schistosome chronobiology. The situation in which the same cercariae "clone" (since issued from only one miracidium) exhibits distinct periods of cercarial emission during the same day was demonstrated for the first generation of offspring, during laboratory crossbreeding experiments, between two S. mansoni chronotypes with two different cercarial emission patterns: early diurnal and late diurnal [START_REF] Théron | Analyse génétique du rythme d'émergence des cercaires de Schistosoma mansoni par croisement de souches à pics d'émission précoces ou tardifs[END_REF]. These authors hypothesized genetic support for the cercarial emission with genetic codominance (the two phenotypes are expressed). More data are now needed regarding the frequency of the M. natalensis-derived chronotypes observed in the studied area and on their genetic and transcriptomic determinants. With the specific mechanisms still unknown, we can still propose a scenario for the various events which led to the establishment of the various chronotypes described in Kessounou (Benin), including those of [START_REF] Savassi | Cattle as natural host for Schistosoma haematobium (Bilharz, 1852) Weinland, 1858 x Schistosoma bovis Sonsino, 1876 interactions, with new cercarial emergence and genetic patterns[END_REF]. Five chronotypes were identified (Fig. 5): onepeak S. bovis typical chronotype (Fig. 5a), one-peak S. haematobium typical chronotype (Fig. 5b), two-peak S. bovis and nocturnal chronotype (Fig. 5c), two-peak S. haematobium and nocturnal chronotype (Fig. 5d), and three-peak S. bovis, S. haematobium, and nocturnal chronotype (Fig. 5e). It can be noted that two other possible chronotypes are absent: the chronotype with the nocturnal emission alone and the two-peak S. bovis and S. haematobium chronotype. Molecular analysis based on two genes (COI and ITS) has shown that the species S. bovis is only present in cattle, the species S. haematobium is absent, and the introgressive hybrid S. haematobium x S. bovis is present in cattle, schoolchildren, and M. natalensis (this paper and [START_REF] Savassi | Cattle as natural host for Schistosoma haematobium (Bilharz, 1852) Weinland, 1858 x Schistosoma bovis Sonsino, 1876 interactions, with new cercarial emergence and genetic patterns[END_REF]). The species S. haematobium seems to give way to its descendants only through introgressive hybridization with S. bovis, descendants which seem adapted to the three hosts. However, it can be noted that the cercarial emission pattern of the S. haematobium is maintained in the introgressive hybrid descendants (see Fig. 5b). Concerning the time position of the different peaks, themorning peak is comparable to a typical S. bovis peak [START_REF] Mouahid | Schistosoma bovis: patterns of cercarial emergence from snails of the genera Bulinus and Planorbarius[END_REF][START_REF] Mouahid | Analyse comparative du rythme d'émission des cercaires de trois souches de Schistosoma bovis[END_REF][START_REF] Mouahid | Cercarial shedding patterns of Schistosoma bovis and Schistosoma haematobium from single and mixed infections of Bulinus truncatus[END_REF]; in this case, the cercariae are shed when cattle go to watering. The afternoon peak is comparable to a typical S. haematobium peak [START_REF] Mouahid | Cercarial shedding patterns of Schistosoma bovis and Schistosoma haematobium from single and mixed infections of Bulinus truncatus[END_REF]; in this case, the cercariae are shed when humans go to water for bathing, washing, or fishing. The nocturnal peak is comparable to a typical S. rodhaini Brumpt 1931 peak and to the nocturnal chronotype of S. mansoni [START_REF] Théron | Hybrids between Schistosoma mansoni and S. rodhaini: characterization by cercarial emergence rhythms[END_REF][START_REF] Mouahid | A new chronotype of Schistosoma mansoni: adaptive significance[END_REF]; in this case, the cercariae are shed when nocturnal rodents actively use water.

In conclusion, our results showed for the first time that a rodent, precisely Mastomys natalensis, can act as a definitive host for transmission of introgressive hybrids between S. haematobium and S. bovis. This result was obtained using different complementary approaches: egg morphology, cercarial emergence patterns, and molecular biology. The morphological parameter did not constitute a precise tool for diagnose but was useful to alert us on the presence of unusual schistosome egg morphology in rodents since around 50% of eggs were typical S. haematobium eggs and the other 50% were of intermediate morphotype between S. haematobium and S. bovis eggs. The cercarial emission pattern, even though obtained from only two snails infected with one miracidium of M. natalensis-derived schistosomes, showed that the typical cercarial emergence patterns of S. haematobium and S. bovis could be found but associated with other periods of emergence, notably a nocturnal one; this abled us to make the link with the previous nocturnal cercarial shedding that we could observe in bovines [START_REF] Savassi | Cattle as natural host for Schistosoma haematobium (Bilharz, 1852) Weinland, 1858 x Schistosoma bovis Sonsino, 1876 interactions, with new cercarial emergence and genetic patterns[END_REF]. Finally, the molecular parameter, based on both mitochondrial and nuclear genes, showed that those M. natalensis-derived schistosomes were introgressive hybrids between S. haematobium and S. bovis. More generally, schistosomes' adaptive strategy towards their different hosts aims to optimize host ecological niche, host spectrum, and cercarial emission time period. In Benin, intra-host ecological niche is optimally used by the S. haematobium x S. bovis individuals since, in children, both the urine and feces are used for parasite transmission. The adaptive use of a multiple-host spectrum also seems optimized since humans as well as Mastomys (this manuscript) and bovines [START_REF] Savassi | Cattle as natural host for Schistosoma haematobium (Bilharz, 1852) Weinland, 1858 x Schistosoma bovis Sonsino, 1876 interactions, with new cercarial emergence and genetic patterns[END_REF] share the same genetic variants (see Fig. 4). Finally, we showed that parasite time period transmission within a single day is also maximized since daily (morning to late afternoon) and nocturnal cercarial shedding is observed for the same single schistosome genotype. Table 1 Visited sites, micro-mammals trapped (Crocidura olivieri, Mastomys natalensis (Mn), Rattus norvegicus, Rattus rattus (Rr)), and infection with schistosomes Table 2 The 14 haplotypes of COI mtDNA gene harbored by schistosomes from rodents and human Table 3 The 11 profiles of ITS rDNA region harbored by the schistosome from rodents and humans

Figure 1 Figure 2 Figure 3 Figure 4 Figure 5

 12345 Figure 1 Geographical positions of the villages where samples were collected. a Position of the area in South Benin; b zoom view and position of the seven villages (white circles)

Acknowledgments

Thanks to Dr Gualbert Houéménou (Laboratoire de Recherche en Biologie Appliquée, Unité de Recherche sur les Invasions Biologiques; Ecole d'Abomey-Calavi, Université d'Abomey-Calavi, Benin) for welcoming us into his Lab and providing us with excellent working conditions. We thank all the councilors (Marius Codjo, Rémi Hounkanni, Nathael Houssou, Omer Kouton, and Mathieu Zingbé), the boatmen (Boniface Ahissouvou, and Moïse Hounvi-Edjê), and local volunteers for providing us with all the facilities during our study, as well as the director of the Kessounou primary public school (Pr Jean-Pierre Tokpo) and his schoolchildren for their collaboration. A big thanks to Jeanine Almany for the linguistic revision of the paper.

Funding

This research was funded by the French Ministry of Foreign Affairs, the International Joint Laboratory through the CONS-HELM project, the French Institute for Sustainable Development (IRD), the French National Agency for Research (ANR) [grant ANR-17-CE12-0005-01] CHRONOGET and CNRS. BS is a student fellow from the SCAC (Service de Coopération et d'Action Culturelle), French Embassy, in Benin.

Data availability

Sequence data were deposited in the NCBI GenBank database under the accession numbers MW022134 to MW022147 for COI mtDNA and MW027648 to MW027658 for ITS rDNA.

Author contribution

Conceptualization: Gabriel Mouahid, Gauthier Dobigny, Hélène Moné Field work: Boris A.E.S. Savassi, Gabriel Mouahid, Gauthier Dobigny, Jonas R. Etougbétché Material preparation, data collection, and analysis: Boris A.E.S.

Informed consent

The Director of Kessounou school (Pr Jean-Pierre Tokpo) obtained informed consent from the parents of all the schoolchildren who participated in the study and all the schoolchildren were done chemotherapy after the study. Study objectives and protocols were systematically explained to the councilors in order to obtain their formal oral agreement to work within the areas under their responsibility. Oral agreements to enter home were also formally obtained by all inhabitants before setting traps inside their home.

Conflict of interest

The authors declare no competing interests.