Quantitative estimates of glacial refugia for chimpanzees (Pan troglodytes) since the Last Interglacial (120,000 BP).

Paleoclimate reconstructions have enhanced our understanding of how past climates have shaped present-day biodiversity. We hypothesize that the geographic extent of Pleistocene forest refugia and suitable habitat fluctuated significantly in time during the late Quaternary for chimpanzees (Pan troglodytes). Using bioclimatic variables representing monthly temperature and precipitation estimates, past human population density data, and an extensive database of georeferenced presence points, we built a model of changing habitat suitability for chimpanzees at fine spatio-temporal scales dating back to the Last Interglacial (120,000 BP). Our models cover a spatial resolution of 0.0467° (approximately 5.19 km2 grid cells) and a temporal resolution of between 1000 and 4000 years. Using our model, we mapped habitat stability over time using three approaches, comparing our modeled stability estimates to existing knowledge of Afrotropical refugia, as well as contemporary patterns of major keystone tropical food resources used by chimpanzees, figs (Moraceae), and palms (Arecacae). Results show habitat stability congruent with known glacial refugia across Africa, suggesting their extents may have been underestimated for chimpanzees, with potentially up to approximately 60,000 km2 of previously unrecognized glacial refugia. The refugia we highlight coincide with higher species richness for figs and palms. Our results provide spatio-temporally explicit insights into the role of refugia across the chimpanzee range, forming the empirical foundation for developing and testing hypotheses about behavioral, ecological, and genetic diversity with additional data. This methodology can be applied to other species and geographic areas when sufficient data are available.

Using MALDI-TOF mass spectrometry to identify ticks collected on domestic and wild animals from the Democratic Republic of the Congo.

Matrix-assisted laser desorption/ionization mass spectrometry (MALDI-TOF MS) has recently emerged as an alternative to morphological and molecular tools to identify tick species. In this study, we set out to evaluate and confirm the ability of MALDI-TOF MS to identify different species of ticks collected in the Democratic Republic of the Congo and preserved in 70% ethanol. A total of 575 ticks, of which 530 were collected from domestic pigs and 45 from wild animals, were subjected to MALDI-TOF MS analysis to evaluate the intraspecies reproducibility and interspecies specificity of MS profiles obtained from the different species. Morphologically, the ticks belonged to seven different species, namely Rhipicephalus complanatus, Rhipicephalus congolensis, Haemaphysalis muhsamae, Ixodes cumulatimpunctatus, Amblyomma exornatum, Amblyomma compressum and an unidentified Rhipicephalus sp. A total of 535/575 (93%) of the spectra obtained were of good enough quality to be used for our analyses. Our home-made MALDI-TOF MS arthropod database was upgraded with spectra obtained from between one and five randomly selected specimens per species. For these reference specimens, molecular identification of the ticks was also made using 16S, 12S rDNA genes and the Cox1 mtDNA gene sequencing. The remaining good quality spectra were then queried against the upgraded MALDI-TOF MS database, showing that 100% were in agreement with the morphological identification, with logarithmic score values (LSVs) between 1.813 and 2.51. The consistency between our morphological, molecular and MALDI-TOF MS identification confirms the capability and precision of MALDI-TOF MS for tick identification.

First evidence of a new simian adenovirus clustering with Human mastadenovirus F viruses.

BACKGROUND: Adenoviruses play an important role as human pathogens, though most infections are believed to be asymptomatic. The over 100 human adenovirus types are classified into seven species (A-G), some of which include simian adenoviruses. Recent findings have highlighted that simian adenoviruses have a zoonotic potential and that some human adenoviruses are likely the result of relatively recent spillover events. METHODS: In order to evaluate the risks associated with primates hunted and sold as bushmeat, multiple samples from 24 freshly killed monkeys were collected in the Republic of the Congo and tested for adenovirus DNA by PCRs targeting the conserved DNA polymerase and hexon genes. RESULTS: The DNA of a novel simian adenovirus was detected in a moustached monkey (Cercopithecus cephus) by the DNA polymerase PCR, but not by the hexon PCR. The 275 nucleotide amplicon was most closely related to members of the Human mastadenovirus F species (93% HAdV-40 and 89% HAdV-41 amino acid identity), rather than to other known simian adenoviruses. CONCLUSIONS: The phylogenetic clustering with Human mastadenovirus F sequences suggests a common ancestor, more recent than the last common ancestor of humans and moustached monkeys. The findings increase concerns about the zoonotic potential of simian adenoviruses and highlight the need for more research and surveillance on the issue.

Detecting plague-host abundance from space: Using a spectral vegetation index to identify occupancy of great gerbil burrows.

In Kazakhstan, plague outbreaks occur when its main host, the great gerbil, exceeds an abundance threshold. These live in family groups in burrows, which can be mapped using remote sensing. Occupancy (percentage of burrows occupied) is a good proxy for abundance and hence the possibility of an outbreak. Here we use time series of satellite images to estimate occupancy remotely. In April and September 2013, 872 burrows were identified in the field as either occupied or empty. For satellite images acquired between April and August, 'burrow objects' were identified and matched to the field burrows. The burrow objects were represented by 25 different polygon types, then classified (using a majority vote from 10 Random Forests) as occupied or empty, using Normalized Difference Vegetation Indices (NDVI) calculated for all images. Throughout the season NDVI values were higher for empty than for occupied burrows. Occupancy status of individual burrows that were continuously occupied or empty, was classified with producer's and user's accuracy values of 63 and 64% for the optimum polygon. Occupancy level was predicted very well and differed 2% from the observed occupancy. This establishes firmly the principle that occupancy can be estimated using satellite images with the potential to predict plague outbreaks over extensive areas with much greater ease and accuracy than previously.

Variable effects of host characteristics on species richness of flea infracommunities in rodents from three continents.

We studied the effect of host gender and body mass on species richness of flea infracommunities in nine rodent host species from three biomes (temperate zone of central Europe, desert of the Middle East and the tropics of East Africa). Using season- and species-specific generalized linear mixed models and controlling for year-to-year variation, spatial clustering of rodent sampling and over-dispersion of the data, we found inconsistent associations between host characteristics and flea species richness. We found strong support for male-biased flea parasitism, especially during the reproductive period (higher species richness in male hosts than in females) in all considered European rodents (Apodemus agrarius, Myodes glareolus and Microtus arvalis) and in one rodent species from the Middle East (Dipodillus dasyurus). In contrast, two of three African rodent species (Lophuromys kilonzoi and Praomys delectorum) demonstrated a trend of female-biased flea species richness. Positive associations between body mass and the number of flea species were detected mainly in males (five of nine species: A. agrarius, M. glareolus, M. arvalis, D. dasyurus and Mastomys natalensis) and not in females (except for M. natalensis). The results of this study support earlier reports that gender-biased, in general, and male-biased, in particular, infestation by ectoparasites is not a universal rule. This suggests that mechanisms of parasite acquisition by an individual host are species-specific and have evolved independently in different rodent host-flea systems.

Local genetic adaptation to habitat in wild chimpanzees.

How populations adapt to their environment is a fundamental question in biology. Yet we know surprisingly little about this process, especially for endangered species such as non-human great apes. Chimpanzees, our closest living relatives, are particularly interesting because they inhabit diverse habitats, from rainforest to woodland-savannah. Whether genetic adaptation facilitates such habitat diversity remains unknown, despite having wide implications for evolutionary biology and conservation. Using 828 newly generated exomes from wild chimpanzees, we find evidence of fine-scale genetic adaptation to habitat. Notably, adaptation to malaria in forest chimpanzees is mediated by the same genes underlying adaptation to malaria in humans. This work demonstrates the power of non-invasive samples to reveal genetic adaptations in endangered populations and highlights the importance of adaptive genetic diversity for chimpanzees.

Plague and climate: scales matter.

Plague is enzootic in wildlife populations of small mammals in central and eastern Asia, Africa, South and North America, and has been recognized recently as a reemerging threat to humans. Its causative agent Yersinia pestis relies on wild rodent hosts and flea vectors for its maintenance in nature. Climate influences all three components (i.e., bacteria, vectors, and hosts) of the plague system and is a likely factor to explain some of plague's variability from small and regional to large scales. Here, we review effects of climate variables on plague hosts and vectors from individual or population scales to studies on the whole plague system at a large scale. Upscaled versions of small-scale processes are often invoked to explain plague variability in time and space at larger scales, presumably because similar scale-independent mechanisms underlie these relationships. This linearity assumption is discussed in the light of recent research that suggests some of its limitations.

Sex-biased parasitism is not universal: evidence from rodent-flea associations from three biomes.

The distribution of parasites among individual hosts is characterised by high variability that is believed to be a result of variations in host traits. To find general patterns of host traits affecting parasite abundance, we studied flea infestation of nine rodent species from three different biomes (temperate zone of central Europe, desert of Middle East and tropics of East Africa). We tested for independent and interactive effects of host sex and body mass on the number of fleas harboured by an individual host while accounting for spatial clustering of host and parasite sampling and temporal variation. We found no consistent patterns of the effect of host sex and body mass on flea abundance either among species within a biome or among biomes. We found evidence for sex-biased flea infestation in just five host species (Apodemus agrarius, Myodes glareolus, Microtus arvalis, Gerbillus andersoni, Mastomys natalensis). In six rodent species, we found an effect of body mass on flea abundance (all species mentioned above and Meriones crassus). This effect was positive in five species and negative in one species (Microtus arvalis). In M. glareolus, G. andersoni, M. natalensis, and M. arvalis, the relationship between body mass and flea abundance was mediated by host sex. This was manifested in steeper change in flea abundance with increasing body mass in male than female individuals (M. glareolus, G. andersoni, M. natalensis), whereas the opposite pattern was found in M. arvalis. Our findings suggest that sex and body mass are common determinants of parasite infestation in mammalian hosts, but neither of them follows universal rules. This implies that the effect of host individual characteristics on mechanisms responsible for flea acquisition may be manifested differently in different host species.

Potential corridors and barriers for plague spread in Central Asia.

BACKGROUND: Plague (Yersinia pestis infection) is a vector-borne disease which caused millions of human deaths in the Middle Ages. The hosts of plague are mostly rodents, and the disease is spread by the fleas that feed on them. Currently, the disease still circulates amongst sylvatic rodent populations all over the world, including great gerbil (Rhombomys opimus) populations in Central Asia. Great gerbils are social desert rodents that live in family groups in burrows, which are visible on satellite images. In great gerbil populations an abundance threshold exists, above which plague can spread causing epizootics. The spatial distribution of the host species is thought to influence the plague dynamics, such as the direction of plague spread, however no detailed analysis exists on the possible functional or structural corridors and barriers that are present in this population and landscape. This study aims to fill that gap. METHODS: Three 20 by 20 km areas with known great gerbil burrow distributions were used to analyse the spatial distribution of the burrows. Object-based image analysis was used to map the landscape at several scales, and was linked to the burrow maps. A novel object-based method was developed - the mean neighbour absolute burrow density difference (MNABDD) - to identify the optimal scale and evaluate the efficacy of using landscape objects as opposed to square cells. Multiple regression using raster maps was used to identify the landscape-ecological variables that explain burrow density best. Functional corridors and barriers were mapped using burrow density thresholds. Cumulative resistance of the burrow distribution to potential disease spread was evaluated using cost distance analysis. A 46-year plague surveillance dataset was used to evaluate whether plague spread was radially symmetric. RESULTS: The burrow distribution was found to be non-random and negatively correlated with Greenness, especially in the floodplain areas. Corridors and barriers showed a mostly NWSE alignment, suggesting easier spreading along this axis. This was confirmed by the analysis of the plague data. CONCLUSIONS: Plague spread had a predominantly NWSE direction, which is likely due to the NWSE alignment of corridors and barriers in the burrow distribution and the landscape. This finding may improve predictions of plague in the future and emphasizes the importance of including landscape analysis in wildlife disease studies.

Human interactions with bats and bat coronaviruses in rural Côte d'Ivoire.

Bats are presumed reservoirs of diverse α- and ß- coronaviruses (CoVs) and understanding the diversity of bat-CoVs and the role bats play in CoV transmission is highly relevant in the context of the current COVID pandemic. We sampled bats in Côte d'Ivoire (2016-2018) living at ecotones between anthropogenic and wild habitats in the Marahoué National Park, a recently encroached protected area, to detect and characterize the CoVs circulating in bats and humans. A total of 314 bats were captured, mostly during the rainy season (78%), and CoV RNA was detected in three of the bats (0.96%). A CoV RNA sequence similar to Chaerephon bat coronavirus/Kenya/KY22/2006 (BtKY22) was found in a Chaerephon cf. pumilus and a Mops sp. fecal swab, while a CoV RNA sequence similar to the two almost identical Kenya bat coronaviruses BtKY55 and BtKY56 (BtKY55/56) was detected in an Epomops buettikoferi oral swab. Phylogenetic analyses indicated differences in the degree of evolutionary host-virus co-speciation for BtKY22 and BtKY55/56. To assess potential for human exposure to these viruses, we conducted human syndromic and community-based surveillance in clinics and high-risk communities. We collected data on participant characteristics, livelihoods, animal contact, and high-risk behaviors that may be associated with exposure to zoonotic diseases. We then collected biological samples for viral testing from 401 people. PCR testing of these biological samples revealed no evidence of CoV infection among the enrolled individuals. We identified higher levels of exposure to bats in people working in crop production and in hunting, trapping and fishing. Finally, we used the 'Spillover' risk-ranking tool to assess the potential for viral spillover and concluded that, while there is no evidence to suggest imminent risk of spillover for these CoVs, their host range and other traits suggest caution and vigilance are warranted in people with high exposure risk.

Orbivirus RNA in a Banana Serotine (Afronycteris nanus) Bat in the Republic of the Congo.

Orbiviruses are arthropod borne viruses of vertebrates, with some of them being important pathogens of veterinary, conservation and economic importance, while others are occasionally associated with human disease. Some apparently bat specific orbiviruses have been detected, but little is known about their distribution and diversity. We thus sampled and screened 52 bats living in the Congo Basin, and detected RNA indicative of a novel orbivirus in a single banana serotine (Afronycteris nanus) by PCR. The detected RNA clusters with epizootic haemorrhagic disease virus, bluetongue virus, and others. The findings highlight the need for more studies into arbovirus presence and diversity in bat species.

Identification of the onchocerciasis vector in the Kakoi-Koda focus of the Democratic Republic of Congo.

BACKGROUND: The objective of this study was to characterise the vector in a small hyper-endemic focus of onchocerciasis (the Kakoi-Koda focus) which has recently been discovered on the western slopes of the rift valley above Lake Albert. METHODOLOGY/PRINCIPAL FINDINGS: Aquatic stages of blackflies were collected by hand from streams and rivers, and anthropophilic adult females were collected by human landing catches. Using a combination of morphotaxonomy and DNA barcoding, the blackflies collected biting humans within the focus were identified as Simulium dentulosum and Simulium vorax, which were also found breeding in local streams and rivers. Simulium damnosum s.l., Simulium neavei and Simulium albivirgulatum were not found (except for a single site in 2009 where crabs were carrying S. neavei). Anthropophilic specimens from the focus were screened for Onchocerca DNA using discriminant qualitative real-time triplex PCR. One specimen of S. vorax was positive for Onchocerca volvulus in the body, and out of 155 S. dentulosum, 30% and 11% were infected and infective (respectively). CONCLUSIONS/SIGNIFICANCE: Simulium dentulosum currently appears to be the main vector of human onchocerciasis within the Kakoi-Koda focus, and S. vorax may be a secondary vector. It remains possible that S. neavei was the main (or only) vector in the past having now become rare as a result of the removal of tree-cover and land-use changes. Simulium vorax has previously been shown to support the development of O. volvulus in the laboratory, but this is the first time that S. dentulosum has been implicated as a probable vector of onchocerciasis, and this raises the possibility that other blackfly species which are not generally considered to be anthropophilic vectors might become vectors under suitable conditions. Because S. dentulosum is not a vector in endemic areas surrounding the Kakoi-Koda focus, it is probable that the Kakoi-Koda focus is significantly isolated.

Correction: Dusabimana et al. Surveillance for Onchocerciasis-Associated Epilepsy and Ov16 IgG4 Testing of Children 6-10 Years Old Should Be Used to Identify Areas Where Onchocerciasis Elimination Programs Need Strengthening. Pathogens 2022, 11, 281.

There was an error in the original publication [...].

Enhanced surveillance of monkeypox in Bas-Uélé, Democratic Republic of Congo: the limitations of symptom-based case definitions.

BACKGROUND: Following an outbreak of cases of vesicular-pustular rash with fever, evocative of human monkeypox, in Bas-Uélé province, Democratic Republic of Congo, surveillance was strengthened. METHODS: Households with at least one active generalized vesicular-pustular rash case were visited, and contact and clinical history information were collected from all household members. Whenever possible, skin lesions were screened by polymerase chain reaction for the monkeypox virus, followed by the varicella-zoster virus, when negative for the former. RESULTS: Polymerase chain reaction results were obtained for 77 suspected cases, distributed in 138 households, of which 27.3% were positive for monkeypox, 58.4% positive for chickenpox, and 14.3% negative for both. Confirmed monkeypox cases presented more often with monomorphic skin lesions on the palms of the hands and on the soles of the feet. Integrating these three features into the case definition raised the specificity to 85% but would miss 50% of true monkeypox cases. A predictive model fit on patient demographics and symptoms had 97% specificity and 80% sensitivity but only 80% and 33% in predicting out-of-sample cases. CONCLUSION: Few discriminating features were identified and the performance of clinical case definitions was suboptimal. Rapid field diagnostics are needed to optimize worldwide early detection and surveillance of monkeypox.

Surveillance for Onchocerciasis-Associated Epilepsy and OV16 IgG4 Testing of Children 6-10 Years Old Should Be Used to Identify Areas Where Onchocerciasis Elimination Programs Need Strengthening.

To eliminate onchocerciasis-associated morbidity, it is important to identify areas where there is still high ongoing Onchocerca volvulus transmission. Between 2015 and 2021, door-to-door surveys were conducted in onchocerciasis-endemic villages in Cameroon, the Democratic Republic of Congo (DRC), Nigeria, South Sudan, and Tanzania to determine epilepsy prevalence and incidence, type of epilepsy and ivermectin therapeutic coverage. Moreover, children aged between six and 10 years were tested for anti-Onchocerca antibodies using the Ov16 IgG4 rapid diagnostic test (RDT). A mixed-effect binary logistic regression analysis was used to assess significantly associated variables of Ov16 antibody seroprevalence. A high prevalence and incidence of epilepsy was found to be associated with a high Ov16 antibody seroprevalence among 6-10-year-old children, except in the Logo health zone, DRC. The low Ov16 antibody seroprevalence among young children in the Logo health zone, despite a high prevalence of epilepsy, may be explained by a recent decrease in O. volvulus transmission because of a decline in the Simulium vector population as a result of deforestation. In the Central African Republic, a new focus of O. volvulus transmission was detected based on the high Ov16 IgG4 seropositivity among children and the detecting of nodding syndrome cases, a phenotypic form of onchocerciasis-associated epilepsy (OAE). In conclusion, Ov16 IgG4 RDT testing of 6-10-year-old children is a cheap and rapid method to determine the level of ongoing O. volvulus transmission and to assess, together with surveillance for OAE, the performance of onchocerciasis elimination programs.

Population dynamics and genetic connectivity in recent chimpanzee history.

Knowledge on the population history of endangered species is critical for conservation, but whole-genome data on chimpanzees (Pan troglodytes) is geographically sparse. Here, we produced the first non-invasive geolocalized catalog of genomic diversity by capturing chromosome 21 from 828 non-invasive samples collected at 48 sampling sites across Africa. The four recognized subspecies show clear genetic differentiation correlating with known barriers, while previously undescribed genetic exchange suggests that these have been permeable on a local scale. We obtained a detailed reconstruction of population stratification and fine-scale patterns of isolation, migration, and connectivity, including a comprehensive picture of admixture with bonobos (Pan paniscus). Unlike humans, chimpanzees did not experience extended episodes of long-distance migrations, which might have limited cultural transmission. Finally, based on local rare variation, we implement a fine-grained geolocalization approach demonstrating improved precision in determining the origin of confiscated chimpanzees.

Male hosts drive infracommunity structure of ectoparasites.

We studied the co-occurrence of flea species in infracommunities of 16 rodents from four regions (South Africa, Tanzania, central Europe and western Siberia) using null models, and predicted that flea co-occurrences will be expressed more strongly in male than in female hosts. We examined patterns of co-occurrence (measured as the C score) in infracommunities of fleas that are parasitic on male and female hosts by comparing co-occurrence frequencies with those expected by chance. When a significant degree of nonrandomness in flea co-occurrences was detected, it indicated aggregative infracommunity structure. In Tanzanian rodents, no significant flea co-occurrences were detected in either male or female hosts. In a South African rodent, significant flea co-occurrences were not detected in males, but were found in females in some localities. In Palaearctic rodents, significant nonrandomness was detected either equally for males and females or more frequently in males than in females. Meta-analyses demonstrated that the frequency of the detection of nonrandomness in flea co-occurrences was significantly higher in male than in female hosts. The values of the standardized effect size (SES) for the C score differed significantly among host species, but not between host genders. When the Palaearctic hosts were analyzed separately, the effects of both host gender and species appeared to be significant, with the SES values for the C score in males being smaller than those in females. The strength of the gender difference in the manifestation of flea community structure increased with increasing gender difference in flea species richness, and with decreasing gender difference in flea prevalence for the Palaearctic hosts. We conclude that male hosts are the main drivers of flea infracommunity structure. However, the manifestation of gender bias in flea community structure varies among host species, and is likely determined by the pattern of species-specific spatial behavior.

No Evidence Known Viruses Play a Role in the Pathogenesis of Onchocerciasis-Associated Epilepsy. An Explorative Metagenomic Case-Control Study.

Despite the increasing epidemiological evidence that the Onchocerca volvulus parasite is strongly associated with epilepsy in children, hence the name onchocerciasis-associated epilepsy (OAE), the pathophysiological mechanism of OAE remains to be elucidated. In June 2014, children with unprovoked convulsive epilepsy and healthy controls were enrolled in a case control study in Titule, Bas-Uélé Province in the Democratic Republic of the Congo (DRC) to identify risk factors for epilepsy. Using a subset of samples collected from individuals enrolled in this study (16 persons with OAE and 9 controls) plasma, buffy coat, and cerebrospinal fluid (CSF) were subjected to random-primed next-generation sequencing. The resulting sequences were analyzed using sensitive computational methods to identify viral DNA and RNA sequences. Anneloviridae, Flaviviridae, Hepadnaviridae (Hepatitis B virus), Herpesviridae, Papillomaviridae, Polyomaviridae (Human polyomavirus), and Virgaviridae were identified in cases and in controls. Not unexpectedly, a variety of bacteriophages were also detected in all cases and controls. However, none of the identified viral sequences were found enriched in OAE cases, which was our criteria for agents that might play a role in the etiology or pathogenesis of OAE.

The evolutionary history of manatees told by their mitogenomes.

The manatee family encompasses three extant congeneric species: Trichechus senegalensis (African manatee), T. inunguis (Amazonian manatee), and T. manatus (West Indian manatee). The fossil record for manatees is scant, and few phylogenetic studies have focused on their evolutionary history. We use full mitogenomes of all extant manatee species to infer the divergence dates and biogeographical histories of these species and the effect of natural selection on their mitogenomes. The complete mitochondrial genomes of T. inunguis (16,851 bp), T. senegalensis (16,882 bp), and T. manatus (16,882 bp), comprise 13 protein-coding genes, 2 ribosomal RNA genes (rRNA - 12S and 16S), and 22 transfer RNA genes (tRNA), and (D-loop/CR). Our analyses show that the first split within Trichechus occurred during the Late Miocene (posterior mean 6.56 Ma and 95% HPD 3.81-10.66 Ma), followed by a diversification event in the Plio-Pleistocene (posterior mean 1.34 Ma, 95% HPD 0.1-4.23) in the clade composed by T. inunguis and T. manatus; T. senegalensis is the sister group of this clade with higher support values (pp > 0.90). The branch-site test identified positive selection on T. inunguis in the 181st position of the ND4 amino acid gene (LRT = 6.06, p = 0.0069, BEB posterior probability = 0.96). The ND4 gene encodes one subunit of the NADH dehydrogenase complex, part of the oxidative phosphorylation machinery. In conclusion, our results provide novel insight into the evolutionary history of the Trichechidae during the Late Miocene, which was influenced by geological events, such as Amazon Basin formation.