Signatures of adaptation at key insecticide resistance loci in Anopheles gambiae in Southern Ghana revealed by reduced-coverage WGS.

Resistance to insecticides and adaptation to a diverse range of environments present challenges to Anopheles gambiae s.l. mosquito control efforts in sub-Saharan Africa. Whole-genome-sequencing is often employed for identifying the genomic basis underlying adaptation in Anopheles, but remains expensive for large-scale surveys. Reduced coverage whole-genome-sequencing can identify regions of the genome involved in adaptation at a lower cost, but is currently untested in Anopheles mosquitoes. Here, we use reduced coverage WGS to investigate population genetic structure and identify signatures of local adaptation in Anopheles mosquitoes across southern Ghana. In contrast to previous analyses, we find no structuring by ecoregion, with Anopheles coluzzii and Anopheles gambiae populations largely displaying the hallmarks of large, unstructured populations. However, we find signatures of selection at insecticide resistance loci that appear ubiquitous across ecoregions in An. coluzzii, and strongest in forest ecoregions in An. gambiae. Our study highlights resistance candidate genes in this region, and validates reduced coverage WGS, potentially to very low coverage levels, for population genomics and exploratory surveys for adaptation in Anopheles taxa.

Long-Term Balancing Selection and the Genetic Load Linked to the Self-Incompatibility Locus in Arabidopsis halleri and A. lyrata.

Balancing selection is a form of natural selection maintaining diversity at the sites it targets and at linked nucleotide sites. Due to selection favoring heterozygosity, it has the potential to facilitate the accumulation of a "sheltered" load of tightly linked recessive deleterious mutations. However, precisely evaluating the extent of these effects has remained challenging. Taking advantage of plant self-incompatibility as one of the best-understood examples of long-term balancing selection, we provide a highly resolved picture of the genomic extent of balancing selection on the sheltered genetic load. We used targeted genome resequencing to reveal polymorphism of the genomic region flanking the self-incompatibility locus in three sample sets in each of the two closely related plant species Arabidopsis halleri and Arabidopsis lyrata, and used 100 control regions from throughout the genome to factor out differences in demographic histories and/or sample structure. Nucleotide polymorphism increased strongly around the S-locus in all sample sets, but only over a limited genomic region, as it became indistinguishable from the genomic background beyond the first 25-30 kb. Genes in this chromosomal interval exhibited no excess of mutations at 0-fold degenerated sites relative to putatively neutral sites, hence revealing no detectable weakening of the efficacy of purifying selection even for these most tightly linked genes. Overall, our results are consistent with the predictions of a narrow genomic influence of linkage to the S-locus and clarify how natural selection in one genomic region affects the evolution of the adjacent genomic regions.

Genomic landscape of drug response reveals mediators of anthelmintic resistance.

Like other pathogens, parasitic helminths can rapidly evolve resistance to drug treatment. Understanding the genetic basis of anthelmintic drug resistance in parasitic nematodes is key to tracking its spread and improving the efficacy and sustainability of parasite control. Here, we use an in vivo genetic cross between drug-susceptible and multi-drug-resistant strains of Haemonchus contortus in a natural host-parasite system to simultaneously map resistance loci for the three major classes of anthelmintics. This approach identifies new alleles for resistance to benzimidazoles and levamisole and implicates the transcription factor cky-1 in ivermectin resistance. This gene is within a locus under selection in ivermectin-resistant populations worldwide; expression analyses and functional validation using knockdown experiments support that cky-1 is associated with ivermectin survival. Our work demonstrates the feasibility of high-resolution forward genetics in a parasitic nematode and identifies variants for the development of molecular diagnostics to combat drug resistance in the field.

Target-enrichment sequencing yields valuable genomic data for challenging-to-culture bacteria of public health importance.

Genomic data contribute invaluable information to the epidemiological investigation of pathogens of public health importance. However, whole-genome sequencing (WGS) of bacteria typically relies on culture, which represents a major hurdle for generating such data for a wide range of species for which culture is challenging. In this study, we assessed the use of culture-free target-enrichment sequencing as a method for generating genomic data for two bacterial species: (1) Bacillus anthracis, which causes anthrax in both people and animals and whose culture requires high-level containment facilities; and (2) Mycoplasma amphoriforme, a fastidious emerging human respiratory pathogen. We obtained high-quality genomic data for both species directly from clinical samples, with sufficient coverage (>15×) for confident variant calling over at least 80% of the baited genomes for over two thirds of the samples tested. Higher qPCR cycle threshold (Ct) values (indicative of lower pathogen concentrations in the samples), pooling libraries prior to capture, and lower captured library concentration were all statistically associated with lower capture efficiency. The Ct value had the highest predictive value, explaining 52 % of the variation in capture efficiency. Samples with Ct values ≤30 were over six times more likely to achieve the threshold coverage than those with a Ct > 30. We conclude that target-enrichment sequencing provides a valuable alternative to standard WGS following bacterial culture and creates opportunities for an improved understanding of the epidemiology and evolution of many clinically important pathogens for which culture is challenging.

Blood meal analysis of tsetse flies ( Glossina pallidipes: Glossinidae) reveals higher host fidelity on wild compared with domestic hosts.

Background: Changes in climate and land use can alter risk of transmission of parasites between domestic hosts and wildlife, particularly when mediated by vectors that can travel between populations. Here we focused on tsetse flies (genus Glossina), the cyclical vectors for both Human African Trypanosomiasis (HAT) and Animal African Trypanosomiasis (AAT). The aims of this study were to investigate three issues related to G. palldipes from Kenya: 1) the diversity of vertebrate hosts that flies fed on; 2) whether host feeding patterns varied in relation to type of hosts, tsetse feeding behaviour, site or tsetse age and sex; and 3) if there was a relationship between trypanosome detection and host feeding behaviours or host types. Methods: Sources of blood meals of Glossina pallidipes were identified by sequencing of the mitochondrial cytochrome b gene and analyzed in relationship with previously determined trypanosome detection in the same flies. Results: In an area dominated by wildlife but with seasonal presence of livestock (Nguruman), 98% of tsetse fed on single wild host species, whereas in an area including a mixture of resident domesticated animals, humans and wildlife (Shimba Hills), 52% of flies fed on more than one host species. Multiple Correspondence Analysis revealed strong correlations between feeding pattern, host type and site but these were resolved along a different dimension than trypanosome status, sex and age of the flies. Conclusions: Our results suggest that individual G. pallidipes in interface areas may show higher feeding success on wild hosts when available but often feed on both wild and domesticated hosts. This illustrates the importance of G. pallidipes as a vector connecting the sylvatic and domestic cycles of African trypanosomes.

Phylogenomic inference of species and subspecies diversity in the Palearctic salamander genus Salamandra.

The salamander genus Salamandra is widespread across Europe, North Africa, and the Near East and is renowned for its conspicuous and polymorphic colouration and diversity of reproductive modes. The phylogenetic relationships within the genus, and especially in the highly polymorphic species S. salamandra, have been very challenging to elucidate, leaving its real evolutionary history and classification at species and subspecies levels a topic of debate and contention. However, the distribution of diversity and species delimitation within the genus are critically important for identifying evolutionarily significant units for conservation and management, especially in light of threats posed by the pathogenic chytrid fungus Batrachochytrium salamandrivorans that is causing massive declines of S. salamandra populations in central Europe. Here, we conducted a phylogenomic analysis from across the taxonomic and geographic breadth of the genus Salamandra in its entire range. Bayesian, maximum likelihood and network-based phylogenetic analyses of up to 4905 ddRADseq-loci (294,300 nucleotides of sequence) supported the distinctiveness of all currently recognised species (Salamandra algira, S. atra, S. corsica, S. infraimmaculata, S. lanzai, and S. salamandra), and all five species for which we have multiple exemplars were confirmed as monophyletic. Within S. salamandra, two main clades can be distinguished: one clade with the Apenninic subspecies S. s. gigliolii nested within the Iberian S. s. bernardezi/fastuosa; and a second clade comprising all other Iberian, Central and East European subspecies. Our analyses revealed that some of the currently recognized subspecies of S. salamandra are paraphyletic and may require taxonomic revision, with the Central- and Eastern-European subspecies all being poorly differentiated at the analysed genomic markers. Salamandra s. longirostris - sometimes considered a separate species - was nested within S. salamandra, consistent with its subspecies status. The relationships identified within and between Salamandra species provide valuable context for future systematic and biogeographic studies, and help elucidate critical evolutionary units for conservation and taxonomy.

Feral populations of Brassica oleracea along Atlantic coasts in western Europe.

There has been growing emphasis on the role that crop wild relatives might play in supporting highly selected agriculturally valuable species in the face of climate change. In species that were domesticated many thousands of years ago, distinguishing wild populations from escaped feral forms can be challenging, but reintroducing variation from either source could supplement current cultivated forms. For economically important cabbages (Brassicaceae: Brassica oleracea), "wild" populations occur throughout Europe but little is known about their genetic variation or potential as resources for breeding more resilient crop varieties. The main aim of this study was to characterize the population structure of geographically isolated wild cabbage populations along the coasts of the UK and Spain, including the Atlantic range edges. Double-digest restriction-site-associated DNA sequencing was used to sample individual cabbage genomes, assess the similarity of plants from 20 populations, and explore environment-genotype associations across varying climatic conditions. Interestingly, there were no indications of isolation by distance; several geographically close populations were genetically more distinct from each other than to distant populations. Furthermore, several distant populations shared genetic ancestry, which could indicate that they were established by escapees of similar source cultivars. However, there were signals of local adaptation to different environments, including a possible relationship between genetic diversity and soil pH. Overall, these results highlight wild cabbages in the Atlantic region as an important genetic resource worthy of further research into their relationship with existing crop varieties.

Functional colour genes and signals of selection in colour-polymorphic salamanders.

Coloration has been associated with multiple biologically relevant traits that drive adaptation and diversification in many taxa. However, despite the great diversity of colour patterns present in amphibians the underlying molecular basis is largely unknown. Here, we use insight from a highly colour-variable lineage of the European fire salamander (Salamandra salamandra bernardezi) to identify functional associations with striking variation in colour morph and pattern. The three focal colour morphs-ancestral black-yellow striped, fully yellow and fully brown-differed in pattern, visible coloration and cellular composition. From population genomic analyses of up to 4,702 loci, we found no correlations of neutral population genetic structure with colour morph. However, we identified 21 loci with genotype-phenotype associations, several of which relate to known colour genes. Furthermore, we inferred response to selection at up to 142 loci between the colour morphs, again including several that relate to coloration genes. By transcriptomic analysis across all different combinations, we found 196 differentially expressed genes between yellow, brown and black skin, 63 of which are candidate genes involved in animal coloration. The concordance across different statistical approaches and 'omic data sets provide several lines of evidence for loci linked to functional differences between colour morphs, including TYR, CAMK1 and PMEL. We found little association between colour morph and the metabolomic profile of its toxic compounds from the skin secretions. Our research suggests that current ecological and evolutionary hypotheses for the origins and maintenance of these striking colour morphs may need to be revisited.

Changing environments and genetic variation: natural variation in inbreeding does not compromise short-term physiological responses.

Selfing plant lineages are surprisingly widespread and successful in a broad range of environments, despite showing reduced genetic diversity, which is predicted to reduce their long-term evolutionary potential. However, appropriate short-term plastic responses to new environmental conditions might not require high levels of standing genetic variation. In this study, we tested whether mating system variation among populations, and associated changes in genetic variability, affected short-term responses to environmental challenges. We compared relative fitness and metabolome profiles of naturally outbreeding (genetically diverse) and inbreeding (genetically depauperate) populations of a perennial plant, Arabidopsis lyrata, under constant growth chamber conditions and an outdoor common garden environment outside its native range. We found no effect of inbreeding on survival, flowering phenology or short-term physiological responses. Specifically, naturally occurring inbreeding had no significant effects on the plasticity of metabolome profiles, using either multivariate approaches or analysis of variation in individual metabolites, with inbreeding populations showing similar physiological responses to outbreeding populations over time in both growing environments. We conclude that low genetic diversity in naturally inbred populations may not always compromise fitness or short-term physiological capacity to respond to environmental change, which could help to explain the global success of selfing mating strategies.

Molecular detection and genetic characterization of Bartonella species from rodents and their associated ectoparasites from northern Tanzania.

BACKGROUND: Bartonellae are intracellular bacteria, which can cause persistent bacteraemia in humans and a variety of animals. Several rodent-associated Bartonella species are human pathogens but data on their global distribution and epidemiology are limited. The aims of the study were to: 1) determine the prevalence of Bartonella infection in rodents and fleas; 2) identify risk factors for Bartonella infection in rodents; and 3) characterize the Bartonella genotypes present in these rodent and flea populations. METHODS AND RESULTS: Spleen samples collected from 381 rodents representing six different species were tested for the presence of Bartonella DNA, which was detected in 57 individuals (15.0%; 95% CI 11.3-18.5), of three rodent species (Rattus rattus n = 54, Mastomys natalensis n = 2 and Paraxerus flavovottis n = 1) using a qPCR targeting the ssrA gene. Considering R. rattus individuals only, risk factor analysis indicated that Bartonella infection was more likely in reproductively mature as compared to immature individuals (OR = 3.42, p <0.001). Bartonella DNA was also detected in 53 of 193 Xenopsylla cheopis fleas (27.5%: 95% CI 21.3-34.3) collected from R.rattus individuals. Analysis of ssrA and gltA sequences from rodent spleens and ssrA sequences from fleas identified multiple genotypes closely related (≥ 97% similar) to several known or suspected zoonotic Bartonella species, including B. tribocorum, B. rochalimae, B. elizabethae and B. quintana. CONCLUSIONS: The ssrA and gltA sequences obtained from rodent spleens and ssrA sequences obtained from fleas reveal the presence of a diverse set of Bartonella genotypes and increase our understanding of the bartonellae present in Tanzanian. Further studies are needed to fully characterise the prevalence, genotypes and diversity of Bartonella in different host populations and their potential impacts on human health.

Refugia and anthelmintic resistance: Concepts and challenges.

Anthelmintic resistance is a threat to global food security. In order to alleviate the selection pressure for resistance and maintain drug efficacy, management strategies increasingly aim to preserve a proportion of the parasite population in 'refugia', unexposed to treatment. While persuasive in its logic, and widely advocated as best practice, evidence for the ability of refugia-based approaches to slow the development of drug resistance in parasitic helminths is currently limited. Moreover, the conditions needed for refugia to work, or how transferable those are between parasite-host systems, are not known. This review, born of an international workshop, seeks to deconstruct the concept of refugia and examine its assumptions and applicability in different situations. We conclude that factors potentially important to refugia, such as the fitness cost of drug resistance, the degree of mixing between parasite sub-populations selected through treatment or not, and the impact of parasite life-history, genetics and environment on the population dynamics of resistance, vary widely between systems. The success of attempts to generate refugia to limit anthelmintic drug resistance are therefore likely to be highly dependent on the system in hand. Additional research is needed on the concept of refugia and the underlying principles for its application across systems, as well as empirical studies within systems that prove and optimise its usefulness.

Sodalis glossinidius presence in wild tsetse is only associated with presence of trypanosomes in complex interactions with other tsetse-specific factors.

BACKGROUND: Susceptibility of tsetse flies (Glossina spp.) to trypanosomes of both humans and animals has been associated with the presence of the endosymbiont Sodalis glossinidius. However, intrinsic biological characteristics of the flies and environmental factors can influence the presence of both S. glossinidius and the parasites. It thus remains unclear whether it is the S. glossinidius or other attributes of the flies that explains the apparent association. The objective of this study was to test whether the presence of Trypanosoma vivax, T. congolense and T. brucei are related to the presence of S. glossinidius in tsetse flies when other factors are accounted for: geographic location, species of Glossina, sex or age of the host flies. RESULTS: Flies (n = 1090) were trapped from four sites in the Shimba Hills and Nguruman regions in Kenya. Sex and species of tsetse (G. austeni, G. brevipalpis, G. longipennis and G. pallidipes) were determined based on external morphological characters and age was estimated by a wing fray score method. The presence of trypanosomes and S. glossinidius was detected using PCR targeting the internal transcribed spacer region 1 and the haemolysin gene, respectively. Sequencing was used to confirm species identification. Generalised Linear Models (GLMs) and Multiple Correspondence Analysis (MCA) were applied to investigate multivariable associations. The overall prevalence of trypanosomes was 42.1%, but GLMs revealed complex patterns of associations: the presence of S. glossinidius was associated with trypanosome presence but only in interactions with other factors and only in some species of trypanosomes. The strongest association was found for T. congolense, and no association was found for T. vivax. The MCA also suggested only a weak association between the presence of trypanosomes and S. glossinidius. Trypanosome-positive status showed strong associations with sex and age while S. glossinidius-positive status showed a strong association with geographic location and species of fly. CONCLUSIONS: We suggest that previous conclusions about the presence of endosymbionts increasing probability of trypanosome presence in tsetse flies may have been confounded by other factors, such as community composition of the tsetse flies and the specific trypanosomes found in different regions.

RADseq and mate choice assays reveal unidirectional gene flow among three lamprey ecotypes despite weak assortative mating: Insights into the formation and stability of multiple ecotypes in sympatry.

Adaptive divergence with gene flow often results in complex patterns of variation within taxa exhibiting substantial ecological differences among populations. One example where this may have occurred is the parallel evolution of freshwater-resident nonparasitic lampreys from anadromous-parasitic ancestors. Previous studies have focused on transitions between these two phenotypic extremes, but here, we considered more complex evolutionary scenarios where an intermediate freshwater form that remains parasitic is found sympatrically with the other two ecotypes. Using population genomic analysis (restriction-associated DNA sequencing), we found that a freshwater-parasitic ecotype was highly distinct from an anadromous-parasitic form (Qlake-P  = 96.8%, Fst  = 0.154), but that a freshwater-nonparasitic form was almost completely admixed in Loch Lomond, Scotland. Demographic reconstructions indicated that both freshwater populations likely derived from a common freshwater ancestor. However, while the nonparasitic ecotype has experienced high levels of introgression from the anadromous-parasitic ecotype (Qanad-P  = 37.7%), there is no evidence of introgression into the freshwater-parasitic ecotype. Paradoxically, mate choice experiments predicted high potential for gene flow: Males from all ecotypes were stimulated to spawn with freshwater-parasitic females, which released gametes in response to all ecotypes. Differentially fixed single nucleotide polymorphisms identified genes associated with growth and development, which could possibly influence the timing of metamorphosis, resulting in significant ecological differences between forms. This suggests that multiple lamprey ecotypes can persist in sympatry following shifts in adaptive peaks, due to environmental change during their repeated colonization of post-glacial regions, followed by periods of extensive gene flow among such diverging populations.

Restriction associated DNA-genotyping at multiple spatial scales in Arabidopsis lyrata reveals signatures of pathogen-mediated selection.

BACKGROUND: Genome scans based on outlier analyses have revolutionized detection of genes involved in adaptive processes, but reports of some forms of selection, such as balancing selection, are still limited. It is unclear whether high throughput genotyping approaches for identification of single nucleotide polymorphisms have sufficient power to detect modes of selection expected to result in reduced genetic differentiation among populations. In this study, we used Arabidopsis lyrata to investigate whether signatures of balancing selection can be detected based on genomic smoothing of Restriction Associated DNA sequencing (RAD-seq) data. We compared how different sampling approaches (both within and between subspecies) and different background levels of polymorphism (inbreeding or outcrossing populations) affected the ability to detect genomic regions showing key signatures of balancing selection, specifically elevated polymorphism, reduced differentiation and shifts towards intermediate allele frequencies. We then tested whether candidate genes associated with disease resistance (R-gene analogs) were detected more frequently in these regions compared to other regions of the genome. RESULTS: We found that genomic regions showing elevated polymorphism contained a significantly higher density of R-gene analogs predicted to be under pathogen-mediated selection than regions of non-elevated polymorphism, and that many of these also showed evidence for an intermediate site-frequency spectrum based on Tajima's D. However, we found few genomic regions that showed both elevated polymorphism and reduced FST among populations, despite strong background levels of genetic differentiation among populations. This suggests either insufficient power to detect the reduced population structure predicted for genes under balancing selection using sparsely distributed RAD markers, or that other forms of diversifying selection are more common for the R-gene analogs tested. CONCLUSIONS: Genome scans based on a small number of individuals sampled from a wide range of populations were sufficient to confirm the relative scarcity of signatures of balancing selection across the genome, but also identified new potential disease resistance candidates within genomic regions showing signatures of balancing selection that would be strong candidates for further sequencing efforts.

A Genome Resequencing-Based Genetic Map Reveals the Recombination Landscape of an Outbred Parasitic Nematode in the Presence of Polyploidy and Polyandry.

The parasitic nematode Haemonchus contortus is an economically and clinically important pathogen of small ruminants, and a model system for understanding the mechanisms and evolution of traits such as anthelmintic resistance. Anthelmintic resistance is widespread and is a major threat to the sustainability of livestock agriculture globally; however, little is known about the genome architecture and parameters such as recombination that will ultimately influence the rate at which resistance may evolve and spread. Here, we performed a genetic cross between two divergent strains of H. contortus, and subsequently used whole-genome resequencing of a female worm and her brood to identify the distribution of genome-wide variation that characterizes these strains. Using a novel bioinformatic approach to identify variants that segregate as expected in a pseudotestcross, we characterized linkage groups and estimated genetic distances between markers to generate a chromosome-scale F1 genetic map. We exploited this map to reveal the recombination landscape, the first for any helminth species, demonstrating extensive variation in recombination rate within and between chromosomes. Analyses of these data also revealed the extent of polyandry, whereby at least eight males were found to have contributed to the genetic variation of the progeny analyzed. Triploid offspring were also identified, which we hypothesize are the result of nondisjunction during female meiosis or polyspermy. These results expand our knowledge of the genetics of parasitic helminths and the unusual life-history of H. contortus, and enhance ongoing efforts to understand the genetic basis of resistance to the drugs used to control these worms and for related species that infect livestock and humans throughout the world. This study also demonstrates the feasibility of using whole-genome resequencing data to directly construct a genetic map in a single generation cross from a noninbred nonmodel organism with a complex lifecycle.

Assessing S. mansoni prevalence in Biomphalaria snails in the Gombe ecosystem of western Tanzania: the importance of DNA sequence data for clarifying species identification.

BACKGROUND: Snails are essential for the transmission and maintenance of schistosomiasis in endemic areas, as they serve as intermediate hosts for schistosome parasites. A clear understanding of the snail species present, their local distribution and infection status is therefore a prerequisite for effective control of schistosomiasis. The purpose of this study was to establish the infection status and distribution of Schistosoma mansoni in snails in the Gombe area along the shores of Lake Tanganyika in western Tanzania, using both detection of cercarial shedding and molecular approaches. METHODS: Snails were collected from streams located close to human settlements in Gombe National Park, as well as from nearby villages (Kiziba, Mtanga, Mwamgongo and Bugamba) and the largest town in the region (Kigoma). Snails were individually exposed to light to induce shedding of schistosome larvae, which were examined using a compound light microscope. Additionally, the internal transcribed spacer (ITS) region of the ribosomal RNA gene cluster was simultaneously amplified in both snails and their trematodes using a single polymerase chain reaction (PCR) and sequenced to confirm species identification. RESULTS: Snails morphologically identified as Biomphalaria pfeifferi were present in all streams except at Mtanga but their distribution was patchy in both time and space. Sequencing of PCR products indicated that not all snails were B. pfeifferi. None of the snails from Gombe or Bugamba shed schistosome larvae, while larvae were shed at all other sites. Overall, an infection prevalence of only 12% was observed in snails based on cercarial shedding. While 47% of the snails were PCR-positive for the 500 bp ITS fragment, which was predicted to indicate infection with S. mansoni, sequence data demonstrated that these bands are not species-specific and can be amplified from other trematode infections. In addition, a 1000 bp band was amplified in 14% of samples, which was identified as a trematode in the family Derogenidae. CONCLUSIONS: The results support the previous assumption that B. pfeifferi snails may be involved in transmitting schistosomiasis in the area but suggest that the community structure of both snails and trematodes may be more complicated than previously thought. This emphasises the importance of confirming species identifications using sequencing, rather than relying only on PCR-based diagnostics or cercarial shedding.

Estimating the prevalence and intensity of Schistosoma mansoni infection among rural communities in Western Tanzania: The influence of sampling strategy and statistical approach.

BACKGROUND: Schistosoma mansoni is a parasite of major public health importance in developing countries, where it causes a neglected tropical disease known as intestinal schistosomiasis. However, the distribution of the parasite within many endemic regions is currently unknown, which hinders effective control. The purpose of this study was to characterize the prevalence and intensity of infection of S. mansoni in a remote area of western Tanzania. METHODOLOGY/PRINCIPAL FINDINGS: Stool samples were collected from 192 children and 147 adults residing in Gombe National Park and four nearby villages. Children were actively sampled in local schools, and adults were sampled passively by voluntary presentation at the local health clinics. The two datasets were therefore analysed separately. Faecal worm egg count (FWEC) data were analysed using negative binomial and zero-inflated negative binomial (ZINB) models with explanatory variables of site, sex, and age. The ZINB models indicated that a substantial proportion of the observed zero FWEC reflected a failure to detect eggs in truly infected individuals, meaning that the estimated true prevalence was much higher than the apparent prevalence as calculated based on the simple proportion of non-zero FWEC. For the passively sampled data from adults, the data were consistent with close to 100% true prevalence of infection. Both the prevalence and intensity of infection differed significantly between sites, but there were no significant associations with sex or age. CONCLUSIONS/SIGNIFICANCE: Overall, our data suggest a more widespread distribution of S. mansoni in this part of Tanzania than was previously thought. The apparent prevalence estimates substantially under-estimated the true prevalence as determined by the ZINB models, and the two types of sampling strategies also resulted in differing conclusions regarding prevalence of infection. We therefore recommend that future surveillance programmes designed to assess risk factors should use active sampling whenever possible, in order to avoid the self-selection bias associated with passive sampling.

Inferring the shallow phylogeny of true salamanders (Salamandra) by multiple phylogenomic approaches.

The rise of high-throughput sequencing techniques provides the unprecedented opportunity to analyse controversial phylogenetic relationships in great depth, but also introduces a risk of being misinterpreted by high node support values influenced by unevenly distributed missing data or unrealistic model assumptions. Here, we use three largely independent phylogenomic data sets to reconstruct the controversial phylogeny of true salamanders of the genus Salamandra, a group of amphibians providing an intriguing model to study the evolution of aposematism and viviparity. For all six species of the genus Salamandra, and two outgroup species from its sister genus Lyciasalamandra, we used RNA sequencing (RNAseq) and restriction site associated DNA sequencing (RADseq) to obtain data for: (1) 3070 nuclear protein-coding genes from RNAseq; (2) 7440 loci obtained by RADseq; and (3) full mitochondrial genomes. The RNAseq and RADseq data sets retrieved fully congruent topologies when each of them was analyzed in a concatenation approach, with high support for: (1) S. infraimmaculata being sister group to all other Salamandra species; (2) S. algira being sister to S. salamandra; (3) these two species being the sister group to a clade containing S. atra, S. corsica and S. lanzai; and (4) the alpine species S. atra and S. lanzai being sister taxa. The phylogeny inferred from the mitochondrial genome sequences differed from these results, most notably by strongly supporting a clade containing S. atra and S. corsica as sister taxa. A different placement of S. corsica was also retrieved when analysing the RNAseq and RADseq data under species tree approaches. Closer examination of gene trees derived from RNAseq revealed that only a low number of them supported each of the alternative placements of S. atra. Furthermore, gene jackknife support for the S. atra - S. lanzai node stabilized only with very large concatenated data sets. The phylogeny of true salamanders thus provides a compelling example of how classical node support metrics such as bootstrap and Bayesian posterior probability can provide high confidence values in a phylogenomic topology even if the phylogenetic signal for some nodes is spurious, highlighting the importance of complementary approaches such as gene jackknifing. Yet, the general congruence among the topologies recovered from the RNAseq and RADseq data sets increases our confidence in the results, and validates the use of phylotranscriptomic approaches for reconstructing shallow relationships among closely related taxa. We hypothesize that the evolution of Salamandra has been characterized by episodes of introgressive hybridization, which would explain the difficulties of fully reconstructing their evolutionary relationships.