Long-term Small Population Size, Deleterious Variation, and Altitude Adaptation in the Ethiopian Wolf, a Severely Endangered Canid.

Ethiopian wolves, a canid species endemic to the Ethiopian Highlands, have been steadily declining in numbers for decades. Currently, out of 35 extant species, it is now one of the world's most endangered canids. Most conservation efforts have focused on preventing disease, monitoring movements and behavior, and assessing the geographic ranges of sub-populations. Here, we add an essential layer by determining the Ethiopian wolf's demographic and evolutionary history using high-coverage (∼40×) whole-genome sequencing from 10 Ethiopian wolves from the Bale Mountains. We observe exceptionally low diversity and enrichment of weakly deleterious variants in the Ethiopian wolves in comparison with two North American gray wolf populations and four dog breeds. These patterns are consequences of long-term small population size, rather than recent inbreeding. We infer the demographic history of the Ethiopian wolf and find it to be concordant with historic records and previous genetic analyses, suggesting Ethiopian wolves experienced a series of both ancient and recent bottlenecks, resulting in a census population size of fewer than 500 individuals and an estimated effective population size of approximately 100 individuals. Additionally, long-term small population size may have limited the accumulation of strongly deleterious recessive mutations. Finally, as the Ethiopian wolves have inhabited high-altitude areas for thousands of years, we searched for evidence of high-altitude adaptation, finding evidence of positive selection at a transcription factor in a hypoxia-response pathway [CREB-binding protein (CREBBP)]. Our findings are pertinent to continuing conservation efforts and understanding how demography influences the persistence of deleterious variation in small populations.

Understanding the Hidden Complexity of Latin American Population Isolates.

Most population isolates examined to date were founded from a single ancestral population. Consequently, there is limited knowledge about the demographic history of admixed population isolates. Here we investigate genomic diversity of recently admixed population isolates from Costa Rica and Colombia and compare their diversity to a benchmark population isolate, the Finnish. These Latin American isolates originated during the 16th century from admixture between a few hundred European males and Amerindian females, with a limited contribution from African founders. We examine whole-genome sequence data from 449 individuals, ascertained as families to build mutigenerational pedigrees, with a mean sequencing depth of coverage of approximately 36×. We find that Latin American isolates have increased genetic diversity relative to the Finnish. However, there is an increase in the amount of identity by descent (IBD) segments in the Latin American isolates relative to the Finnish. The increase in IBD segments is likely a consequence of a very recent and severe population bottleneck during the founding of the admixed population isolates. Furthermore, the proportion of the genome that falls within a long run of homozygosity (ROH) in Costa Rican and Colombian individuals is significantly greater than that in the Finnish, suggesting more recent consanguinity in the Latin American isolates relative to that seen in the Finnish. Lastly, we find that recent consanguinity increased the number of deleterious variants found in the homozygous state, which is relevant if deleterious variants are recessive. Our study suggests that there is no single genetic signature of a population isolate.

Determining the factors driving selective effects of new nonsynonymous mutations.

The distribution of fitness effects (DFE) of new mutations plays a fundamental role in evolutionary genetics. However, the extent to which the DFE differs across species has yet to be systematically investigated. Furthermore, the biological mechanisms determining the DFE in natural populations remain unclear. Here, we show that theoretical models emphasizing different biological factors at determining the DFE, such as protein stability, back-mutations, species complexity, and mutational robustness make distinct predictions about how the DFE will differ between species. Analyzing amino acid-changing variants from natural populations in a comparative population genomic framework, we find that humans have a higher proportion of strongly deleterious mutations than Drosophila melanogaster. Furthermore, when comparing the DFE across yeast, Drosophila, mice, and humans, the average selection coefficient becomes more deleterious with increasing species complexity. Last, pleiotropic genes have a DFE that is less variable than that of nonpleiotropic genes. Comparing four categories of theoretical models, only Fisher's geometrical model (FGM) is consistent with our findings. FGM assumes that multiple phenotypes are under stabilizing selection, with the number of phenotypes defining the complexity of the organism. Our results suggest that long-term population size and cost of complexity drive the evolution of the DFE, with many implications for evolutionary and medical genomics.

PReFerSim: fast simulation of demography and selection under the Poisson Random Field model.

The Poisson Random Field (PRF) model has become an important tool in population genetics to study weakly deleterious genetic variation under complicated demographic scenarios. Currently, there are no freely available software applications that allow simulation of genetic variation data under this model. Here we present PReFerSim, an ANSI C program that performs forward simulations under the PRF model. PReFerSim models changes in population size, arbitrary amounts of inbreeding, dominance and distributions of selective effects. Users can track summaries of genetic variation over time and output trajectories of selected alleles. AVAILABILITY AND IMPLEMENTATION: PReFerSim is freely available at: https://github.com/LohmuellerLab/PReFerSim CONTACT: klohmueller@ucla.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Genomic Flatlining in the Endangered Island Fox.

Genetic studies of rare and endangered species often focus on defining and preserving genetically distinct populations, especially those having unique adaptations [1, 2]. Much less attention is directed at understanding the landscape of deleterious variation, an insidious consequence of geographic isolation and the inefficiency of natural selection to eliminate harmful variants in small populations [3-5]. With population sizes of many vertebrates decreasing and isolation increasing through habitat fragmentation and loss, understanding the extent and nature of deleterious variation in small populations is essential for predicting and enhancing population persistence. The Channel Island fox (Urocyon littoralis) is a dwarfed species that inhabits six of California's Channel Islands and is derived from the mainland gray fox (U. cinereoargenteus). These isolated island populations have persisted for thousands of years at extremely small population sizes [6, 7] and, consequently, are a model for testing ideas about the accumulation of deleterious variation in small populations under natural conditions. Analysis of complete genome sequence data from island foxes shows a dramatic decrease in genome-wide variation and a sharp increase in the homozygosity of deleterious variants. The San Nicolas Island population has a near absence of variation, demonstrating a unique genetic flatlining that is punctuated by heterozygosity hotspots, enriched for olfactory receptor genes and other genes with high levels of ancestral variation. These findings question the generality of the small-population paradigm that maintains substantial genetic variation is necessary for short- and long-term persistence.

An assessment of the information content of likelihood ratios derived from complex mixtures.

With the increasing sensitivity of DNA typing methodologies, as well as increasing awareness by law enforcement of the perceived capabilities of DNA typing, complex mixtures consisting of DNA from two or more contributors are increasingly being encountered. However, insufficient research has been conducted to characterize the ability to distinguish a true contributor (TC) from a known non-contributor (KNC) in these complex samples, and under what specific conditions. In order to investigate this question, sets of six 15-locus Caucasian genotype profiles were simulated and used to create mixtures containing 2-5 contributors. Likelihood ratios were computed for various situations, including varying numbers of contributors and unknowns in the evidence profile, as well as comparisons of the evidence profile to TCs and KNCs. This work was intended to illustrate the best-case scenario, in which all alleles from the TC were detected in the simulated evidence samples. Therefore the possibility of drop-out was not modeled in this study. The computer program DNAMIX was then used to compute LRs comparing the evidence profile to TCs and KNCs. This resulted in 140,000 LRs for each of the two scenarios. These complex mixture simulations show that, even when all alleles are detected (i.e. no drop-out), TCs can generate LRs less than 1 across a 15-locus profile. However, this outcome was rare, 7 of 140,000 replicates (0.005%), and associated only with mixtures comprising 5 contributors in which the numerator hypothesis includes one or more unknown contributors. For KNCs, LRs were found to be greater than 1 in a small number of replicates (75 of 140,000 replicates, or 0.05%). These replicates were limited to 4 and 5 person mixtures with 1 or more unknowns in the numerator. Only 5 of these 75 replicates (0.004%) yielded an LR greater than 1,000. Thus, overall, these results imply that the weight of evidence that can be derived from complex mixtures containing up to 5 contributors, under a scenario in which no drop-out is required to explain any of the contributors, is remarkably high. This is a useful benchmark result on top of which to layer the effects of additional factors, such as drop-out, peak height, and other variables.

Bottlenecks and selective sweeps during domestication have increased deleterious genetic variation in dogs.

Population bottlenecks, inbreeding, and artificial selection can all, in principle, influence levels of deleterious genetic variation. However, the relative importance of each of these effects on genome-wide patterns of deleterious variation remains controversial. Domestic and wild canids offer a powerful system to address the role of these factors in influencing deleterious variation because their history is dominated by known bottlenecks and intense artificial selection. Here, we assess genome-wide patterns of deleterious variation in 90 whole-genome sequences from breed dogs, village dogs, and gray wolves. We find that the ratio of amino acid changing heterozygosity to silent heterozygosity is higher in dogs than in wolves and, on average, dogs have 2-3% higher genetic load than gray wolves. Multiple lines of evidence indicate this pattern is driven by less efficient natural selection due to bottlenecks associated with domestication and breed formation, rather than recent inbreeding. Further, we find regions of the genome implicated in selective sweeps are enriched for amino acid changing variants and Mendelian disease genes. To our knowledge, these results provide the first quantitative estimates of the increased burden of deleterious variants directly associated with domestication and have important implications for selective breeding programs and the conservation of rare and endangered species. Specifically, they highlight the costs associated with selective breeding and question the practice favoring the breeding of individuals that best fit breed standards. Our results also suggest that maintaining a large population size, rather than just avoiding inbreeding, is a critical factor for preventing the accumulation of deleterious variants.

Plasmodium falciparum infection rates for some Anopheles spp. from Guinea-Bissau, West Africa.

Presence of Plasmodium falciparum circumsporozoite protein (CSP) was detected by enzyme linked immunosorbent assay (ELISA) in a sample of Anopheles gambiae s.s., A. melas and A. pharoensis collected in Guinea-Bissau during October and November 2009. The percentage of P. falciparum infected samples (10.2% overall) was comparable to earlier studies from other sites in Guinea-Bissau (9.6-12.4%). The majority of the specimens collected were identified as A. gambiae which had an individual infection rate of 12.6 % across collection sites. A small number of specimens of A. coluzzii, A. coluzzii x A. gambiae hybrids, A. melas and A. pharoensis were collected and had infection rates of 4.3%, 4.1%, 11.1% and 33.3% respectively. Despite being present in low numbers in indoor collections, the exophilic feeding behaviors of A. melas (N=18) and A. pharoensis (N=6) and high infection rates observed in this survey suggest falciparum-malaria transmission potential outside of the protection of bed nets.

A new multiplex SNP genotyping assay for detecting hybridization and introgression between the M and S molecular forms of Anopheles gambiae.

The M and S forms of Anopheles gambiae have been the subject of intense study, but are morphologically indistinguishable and can only be identified using molecular techniques. PCR-based assays to distinguish the two forms have been designed and applied widely. However, the application of these assays towards identifying hybrids between the two forms, and backcrossed hybrids in particular, has been problematic as the currently available diagnostic assays are based on single locus and/or are located within a multicopy gene. Here, we present an alternative genotyping method for detecting hybridization and introgression between M and S molecular forms based on a multilocus panel of single-nucleotide polymorphisms (SNPs) fixed between the M and S forms. The panel of SNPs employed is located in so-called islands of divergence leading us to describe this method as the 'Divergence Island SNP' (DIS) assay. We show this multilocus SNP genotyping approach can robustly and accurately detect F1 hybrids as well as backcrossed individuals.

Spatiotemporal dynamics of gene flow and hybrid fitness between the M and S forms of the malaria mosquito, Anopheles gambiae.

The M and S forms of Anopheles gambiae have been the focus of intense study by malaria researchers and evolutionary biologists interested in ecological speciation. Divergence occurs at three discrete islands in genomes that are otherwise nearly identical. An "islands of speciation" model proposes that diverged regions contain genes that are maintained by selection in the face of gene flow. An alternative "incidental island" model maintains that gene flow between M and S is effectively zero and that divergence islands are unrelated to speciation. A "divergence island SNP" assay was used to explore the spatial and temporal distributions of hybrid genotypes. Results revealed that hybrid individuals occur at frequencies ranging between 5% and 97% in every population examined. A temporal analysis revealed that assortative mating is unstable and periodically breaks down, resulting in extensive hybridization. Results suggest that hybrids suffer a fitness disadvantage, but at least some hybrid genotypes are viable. Stable introgression of the 2L speciation island occurred at one site following a hybridization event.

A preliminary investigation of the relationship between water quality and Anopheles gambiae larval habitats in Western Cameroon.

BACKGROUND: Water quality and anopheline habitat have received increasing attention due to the possibility that challenges during larval life may translate into adult susceptibility to malaria parasite infection and/or insecticide resistance. METHODS: A preliminary study of Anopheles gambiae s.s. larval habitats in the north-west and south-west regions of Cameroon was conducted in order to detect associations between An. gambiae s.s. molecular form and 2La inversion distributions with basic water quality parameters. Water quality was measured by temperature, pH, conductivity, total dissolved solids (TDS) at seven sites in Cameroon and one site in Selinkenyi, Mali. RESULTS: Principal components and correlation analyses indicated a complex relationship between 2La polymorphism, temperature, conductivity and TDS. Cooler water sites at more inland locations yielded more S form larvae with higher 2La inversion polymorphism while warmer water sites yielded more M form larvae with rare observations of the 2La inversion. DISCUSSION: More detailed studies that take into account the population genetics but also multiple life stages, environmental data relative to these life stages and interactions with both humans and the malaria parasite may help us to understand more about how and why this successful mosquito is able to adapt and diverge, and how it can be successfully managed.

An analysis of two island groups as potential sites for trials of transgenic mosquitoes for malaria control.

Considerable technological advances have been made towards the generation of genetically modified mosquitoes for vector control. In contrast, less progress has been made towards field evaluations of transformed mosquitoes which are critical for evaluating the success of, and hazards associated with, genetic modification. Oceanic islands have been highlighted as potentially the best locations for such trials. However, population genetic studies are necessary to verify isolation. Here, we used a panel of genetic markers to assess for evidence of genetic isolation of two oceanic island populations of the African malaria vector, Anopheles gambiae s.s. We found no evidence of isolation between the Bijagós archipelago and mainland Guinea-Bissau, despite separation by distances beyond the known dispersal capabilities of this taxon. Conversely, the Comoros Islands appear to be genetically isolated from the East African mainland, and thus represent a location worthy of further investigation for field trials. Based on assessments of gene flow within and between the Comoros islands, the island of Grande Comore was found to be genetically isolated from adjacent islands and also exhibited local population structure, indicating that it may be the most suitable site for trials with existing genetic modification technologies.

Chromosome inversions, genomic differentiation and speciation in the African malaria mosquito Anopheles gambiae.

The African malaria vector, Anopheles gambiae, is characterized by multiple polymorphic chromosomal inversions and has become widely studied as a system for exploring models of speciation. Near complete reproductive isolation between different inversion types, known as chromosomal forms, has led to the suggestion that A. gambiae is in early stages of speciation, with divergence evolving in the face of considerable gene flow. We compared the standard chromosomal arrangement (Savanna form) with genomes homozygous for j, b, c, and u inversions (Bamako form) in order to identify regions of genomic divergence with respect to inversion polymorphism. We found levels of divergence between the two sub-taxa within some of these inversions (2Rj and 2Rb), but at a level lower than expected and confined near the inversion breakpoints, consistent with a gene flux model. Unexpectedly, we found that the majority of diverged regions were located on the X chromosome, which contained half of all significantly diverged regions, with much of this divergence located within exons. This is surprising given that the Bamako and Savanna chromosomal forms are both within the S molecular form that is defined by a locus near centromere of X chromosome. Two X-linked genes (a heat shock protein and P450 encoding genes) involved in reproductive isolation between the M and S molecular forms of A. gambiae were also significantly diverged between the two chromosomal forms. These results suggest that genes mediating reproductive isolation are likely located on the X chromosome, as is thought to be the case for the M and S molecular forms. We conclude that genes located on the sex chromosome may be the major force driving speciation between these chromosomal forms of A. gambiae.

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.

Morphological differentiation may mediate mate-choice between incipient species of Anopheles gambiae s.s.

The M and S molecular forms of Anopheles gambiae s.s. have been considered incipient species for more than ten years, yet the mechanism underlying assortative mating of these incipient species has remained elusive. The discovery of the importance of harmonic convergence of wing beat frequency in mosquito mating and its relation to wing size have laid the foundation for exploring phenotypic divergence in wing size of wild populations of the two forms. In this study, wings from field collected mosquitoes were measured for wing length and wing width from two parts of the sympatric distribution, which differ with respect to the strength of assortative mating. In Mali, where assortative mating is strong, as evidenced by low rates of hybridization, mean wing lengths and wing widths were significantly larger than those from Guinea-Bissau. In addition, mean wing widths in Mali were significantly different between molecular forms. In Guinea-Bissau, assortative mating appears comparatively reduced and wing lengths and widths did not differ significantly between molecular forms. The data presented in this study support the hypothesis that wing beat frequency may mediate assortative mating in the incipient species of A. gambiae and represent the first documentation of a morphological difference between the M and S molecular forms.

Asymmetric introgression between the M and S forms of the malaria vector, Anopheles gambiae, maintains divergence despite extensive hybridization.

The suggestion that genetic divergence can arise and/or be maintained in the face of gene flow has been contentious since first proposed. This controversy and a rarity of good examples have limited our understanding of this process. Partially reproductively isolated taxa have been highlighted as offering unique opportunities for identifying the mechanisms underlying divergence with gene flow. The African malaria vector, Anopheles gambiae s.s., is widely regarded as consisting of two sympatric forms, thought by many to represent incipient species, the M and S molecular forms. However, there has been much debate about the extent of reproductive isolation between M and S, with one view positing that divergence may have arisen and is being maintained in the presence of gene flow, and the other proposing a more advanced speciation process with little realized gene flow because of low hybrid fitness. These hypotheses have been difficult to address because hybrids are typically rare (<1%). Here, we assess samples from an area of high hybridization and demonstrate that hybrids are fit and responsible for extensive introgression. Nonetheless, we show that strong divergent selection at a subset of loci combined with highly asymmetric introgression has enabled M and S to remain genetically differentiated despite extensive gene flow. We propose that the extent of reproductive isolation between M and S varies across West Africa resulting in a 'geographic mosaic of reproductive isolation'; a finding which adds further complexity to our understanding of divergence in this taxon and which has considerable implications for transgenic control strategies.