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.

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.

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.

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.

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.

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.

R-gene variation across Arabidopsis lyrata subspecies: effects of population structure, selection and mating system.

BACKGROUND: Examining allelic variation of R-genes in closely related perennial species of Arabidopsis thaliana is critical to understanding how population structure and ecology interact with selection to shape the evolution of innate immunity in plants. We finely sampled natural populations of Arabidopsis lyrata from the Great Lakes region of North America (A. l. lyrata) and broadly sampled six European countries (A. l. petraea) to investigate allelic variation of two R-genes (RPM1 and WRR4) and neutral genetic markers (Restriction Associated DNA sequences and microsatellites) in relation to mating system, phylogeographic structure and subspecies divergence. RESULTS: Fine-scale sampling of populations revealed strong effects of mating system and population structure on patterns of polymorphism for both neutral loci and R-genes, with no strong evidence for selection. Broad geographic sampling revealed evidence of balancing selection maintaining polymorphism in R-genes, with elevated heterozygosity and diversity compared to neutral expectations and sharing of alleles among diverged subspecies. Codon-based tests detected both positive and purifying selection for both R-genes, as commonly found for animal immune genes. CONCLUSIONS: Our results highlight that combining fine and broad-scale sampling strategies can reveal the multiple factors influencing polymorphism and divergence at potentially adaptive genes such as R-genes.

Copy number variation and genetic diversity of MHC Class IIb alleles in an alien population of Xenopus laevis.

Xenopus laevis (the African clawed frog), which originated through hybridisation and whole genome duplication, has been used as a model for genetics and development for many years, but surprisingly little is known about immune gene variation in natural populations. The purpose of this study was to use an isolated population of X. laevis that was introduced to Wales, UK in the past 50 years to investigate how variation at the MHC compares to that at other loci, following a severe population bottleneck. Among 18 individuals, we found nine alleles based on exon 2 sequences of the Class IIb region (which includes the peptide binding region). Individuals carried from one to three of the loci identified from previous laboratory studies. Genetic variation was an order of magnitude higher at the MHC compared with three single-copy nuclear genes, but all loci showed high levels of heterozygosity and nucleotide diversity and there was not an excess of homozygosity or decrease in diversity over time that would suggest extensive inbreeding in the introduced population. Tajima's D was positive for all loci, which is consistent with a bottleneck. Moreover, comparison with published sequences identified the source of the introduced population as the Western Cape region of South Africa, where most commercial suppliers have obtained their stocks. These factors suggest that despite founding by potentially already inbred individuals, the alien population in Wales has maintained substantial genetic variation at both adaptively important and neutral genes.

Behavioural and physiological adaptations to low-temperature environments in the common frog, Rana temporaria.

BACKGROUND: Extreme environments can impose strong ecological and evolutionary pressures at a local level. Ectotherms are particularly sensitive to low-temperature environments, which can result in a reduced activity period, slowed physiological processes and increased exposure to sub-zero temperatures. The aim of this study was to assess the behavioural and physiological responses that facilitate survival in low-temperature environments. In particular, we asked: 1) do high-altitude common frog (Rana temporaria) adults extend the time available for larval growth by breeding at lower temperatures than low-altitude individuals?; and 2) do tadpoles sampled from high-altitude sites differ physiologically from those from low-altitude sites, in terms of routine metabolic rate (RMR) and freeze tolerance? Breeding date was assessed as the first day of spawn observation and local temperature recorded for five, paired high- and low-altitude R. temporaria breeding sites in Scotland. Spawn was collected and tadpoles raised in a common laboratory environment, where RMR was measured as oxygen consumed using a closed respiratory tube system. Freeze tolerance was measured as survival following slow cooling to the point when all container water had frozen. RESULTS: We found that breeding did not occur below 5°C at any site and there was no significant relationship between breeding temperature and altitude, leading to a delay in spawning of five days for every 100 m increase in altitude. The relationship between altitude and RMR varied by mountain but was lower for individuals sampled from high- than low-altitude sites within the three mountains with the highest high-altitude sites (≥900 m). In contrast, individuals sampled from low-altitudes survived freezing significantly better than those from high-altitudes, across all mountains. CONCLUSIONS: Our results suggest that adults at high-altitude do not show behavioural adaptations in terms of breeding at lower temperatures. However, tadpoles appear to have the potential to adapt physiologically to surviving at high-altitude via reduced RMR but without an increase in freeze tolerance. Therefore, survival at high-altitude may be facilitated by physiological mechanisms that permit faster growth rates, allowing completion of larval development within a shorter time period, alleviating the need for adaptations that extend the time available for larval growth.

Recent progress and challenges in population genetics of polyploid organisms: an overview of current state-of-the-art molecular and statistical tools.

Despite the importance of polyploidy and the increasing availability of new genomic data, there remain important gaps in our knowledge of polyploid population genetics. These gaps arise from the complex nature of polyploid data (e.g. multiple alleles and loci, mixed inheritance patterns, association between ploidy and mating system variation). Furthermore, many of the standard tools for population genetics that have been developed for diploids are often not feasible for polyploids. This review aims to provide an overview of the state-of-the-art in polyploid population genetics and to identify the main areas where further development of molecular techniques and statistical theory is required. We review commonly used molecular tools (amplified fragment length polymorphism, microsatellites, Sanger sequencing, next-generation sequencing and derived technologies) and their challenges associated with their use in polyploid populations: that is, allele dosage determination, null alleles, difficulty of distinguishing orthologues from paralogues and copy number variation. In addition, we review the approaches that have been used for population genetic analysis in polyploids and their specific problems. These problems are in most cases directly associated with dosage uncertainty and the problem of inferring allele frequencies and assumptions regarding inheritance. This leads us to conclude that for advancing the field of polyploid population genetics, most priority should be given to development of new molecular approaches that allow efficient dosage determination, and to further development of analytical approaches to circumvent dosage uncertainty and to accommodate 'flexible' modes of inheritance. In addition, there is a need for more simulation-based studies that test what kinds of biases could result from both existing and novel approaches.

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.

Disentangling the effects of breakdown of self-incompatibility and transition to selfing in North American Arabidopsis lyrata.

A breakdown of self-incompatibility (SI) followed by a shift to selfing is commonly observed in the evolution of flowering plants. Both are expected to reduce the levels of heterozygosity and genetic diversity. However, breakdown of SI should most strongly affect the region of the SI locus (S-locus) because of the relaxation of balancing selection that operates on a functional S-locus, and a potential selective sweep. In contrast, a transition to selfing should affect the whole genome. We set out to disentangle the effects of breakdown of SI and transition to selfing on the level and distribution of genetic diversity in North American populations of Arabidopsis lyrata. Specifically, we compared sequence diversity of loci linked and unlinked to the S-locus for populations ranging from complete selfing to fully outcrossing. Regardless of linkage to the S-locus, heterozygosity and genetic diversity increased with population outcrossing rate. High heterozygosity of self-compatible individuals in outcrossing populations suggests that SI is not the only factor preventing the evolution of self-fertilization in those populations. There was a strong loss of diversity in selfing populations, which was more pronounced at the S-locus. In addition, selfing populations showed an accumulation of derived mutations at the S-locus. Our results provide evidence that beyond the genome-wide consequences of the population bottleneck associated with the shift to selfing, the S-locus of A. lyrata shows a specific signal either reflecting the relaxation of balancing selection or positive selection.

Spatial and temporal patterns of neutral and adaptive genetic variation in the endangered African wild dog (Lycaon pictus).

Deciphering patterns of genetic variation within a species is essential for understanding population structure, local adaptation and differences in diversity between populations. Whilst neutrally evolving genetic markers can be used to elucidate demographic processes and genetic structure, they are not subject to selection and therefore are not informative about patterns of adaptive variation. As such, assessments of pertinent adaptive loci, such as the immunity genes of the major histocompatibility complex (MHC), are increasingly being incorporated into genetic studies. In this study, we combined neutral (microsatellite, mtDNA) and adaptive (MHC class II DLA-DRB1 locus) markers to elucidate the factors influencing patterns of genetic variation in the African wild dog (Lycaon pictus); an endangered canid that has suffered extensive declines in distribution and abundance. Our genetic analyses found all extant wild dog populations to be relatively small (N(e) < 30). Furthermore, through coalescent modelling, we detected a genetic signature of a recent and substantial demographic decline, which correlates with human expansion, but contrasts with findings in some other African mammals. We found strong structuring of wild dog populations, indicating the negative influence of extensive habitat fragmentation and loss of gene flow between habitat patches. Across populations, we found that the spatial and temporal structure of microsatellite diversity and MHC diversity were correlated and strongly influenced by demographic stability and population size, indicating the effects of genetic drift in these small populations. Despite this correlation, we detected signatures of selection at the MHC, implying that selection has not been completely overwhelmed by genetic drift.

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.

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.