Chromosome-length genome assemblies and cytogenomic analyses of pangolins reveal remarkable chromosome counts and plasticity.

We report the first chromosome-length genome assemblies for three species in the mammalian order Pholidota: the white-bellied, Chinese, and Sunda pangolins. Surprisingly, we observe extraordinary karyotypic plasticity within this order and, in female white-bellied pangolins, the largest number of chromosomes reported in a Laurasiatherian mammal: 2n = 114. We perform the first karyotype analysis of an African pangolin and report a Y-autosome fusion in white-bellied pangolins, resulting in 2n = 113 for males. We employ a novel strategy to confirm the fusion and identify the autosome involved by finding the pseudoautosomal region (PAR) in the female genome assembly and analyzing the 3D contact frequency between PAR sequences and the rest of the genome in male and female white-bellied pangolins. Analyses of genetic variability show that white-bellied pangolins have intermediate levels of genome-wide heterozygosity relative to Chinese and Sunda pangolins, consistent with two moderate declines of historical effective population size. Our results reveal a remarkable feature of pangolin genome biology and highlight the need for further studies of these unique and endangered mammals.

A highly contiguous genome assembly for the Yellow Warbler (Setophaga petechia).

The Yellow Warbler (Setophaga petechia) is a small songbird in the wood-warbler family (Parulidae) that exhibits phenotypic and ecological differences across a widespread distribution and is important to California's riparian habitat conservation. Here, we present a high-quality de novo genome assembly of a vouchered female Yellow Warbler from southern California. Using HiFi long-read and Omni-C proximity sequencing technologies, we generated a 1.22 Gb assembly including 687 scaffolds with a contig N50 of 6.80 Mb, scaffold N50 of 21.18 Mb, and a BUSCO completeness score of 96.0%. This highly contiguous genome assembly provides an essential resource for understanding the history of gene flow, divergence, and local adaptation in Yellow Warblers and can inform conservation management of this charismatic bird species.

Evidence for ecotone speciation across an African rainforest-savanna gradient.

Accelerating climate change and habitat loss make it imperative that plans to conserve biodiversity consider species' ability to adapt to changing environments. However, in biomes where biodiversity is highest, the evolutionary mechanisms responsible for generating adaptative variation and, ultimately, new species are frequently poorly understood. African rainforests represent one such biome, as decadal debates continue concerning the mechanisms generating African rainforest biodiversity. These debates hinge on the relative importance of geographic isolation versus divergent natural selection across environmental gradients. Hindering progress is a lack of robust tests of these competing hypotheses. Because African rainforests are severely at-risk due to climate change and other anthropogenic activities, addressing this long-standing debate is critical for making informed conservation decisions. We use demographic inference and allele frequency-environment relationships to investigate mechanisms of diversification in an African rainforest skink, Trachylepis affinis, a species inhabiting the gradient between rainforest and rainforest-savanna mosaic (ecotone). We provide compelling evidence of ecotone speciation, in which gene flow has all but ceased between rainforest and ecotone populations, at a level consistent with infrequent hybridization between sister species. Parallel patterns of genomic, morphological, and physiological divergence across this environmental gradient and pronounced allele frequency-environment correlation indicate speciation is mostly probably driven by ecological divergence, supporting a central role for divergent natural selection. Our results provide strong evidence for the importance of ecological gradients in African rainforest speciation and inform conservation strategies that preserve the processes that produce and maintain biodiversity.

Evolutionary and ecological drivers of local adaptation and speciation in a North American avian species complex.

Throughout the speciation process, genomic divergence can be differentially impacted by selective pressures, as well as gene flow and genetic drift. Disentangling the effects of these evolutionary mechanisms remains challenging, especially for nonmodel organisms. Accounting for complex evolutionary histories and contemporary population structure often requires sufficient sample sizes, for which the expense of full genomes remains prohibitive. Here, we demonstrate the utility of partial-genome sequence data for range-wide samples to shed light into the divergence process of two closely related ducks, the Mexican duck (Anas diazi) and mallard (A. platyrhynchos). We determine the role of selective and neutral processes during speciation of Mexican ducks by integrating evolutionary and demographic modelling with genotype-environment and genotype-phenotype association testing. First, evolutionary models and demographic analyses support the hypothesis that Mexican ducks originally diverged ~300,000 years ago in climate refugia arising during a glacial period in southwest North America, and that subsequent environmental selective pressures played a key role in divergence. Mexican ducks then showed cyclical demographic patterns that probably reflected repeated range expansions and contractions, along with bouts of gene flow with mallards during glacial cycles. Finally, we provide evidence that sexual selection acted on several phenotypic traits as a co-evolutionary process, facilitating the development of reproductive barriers that initially arose due to strong ecological selection. More broadly, this work reveals that the genomic and phenotypic patterns observed across species complexes are the result of myriad factors that contribute in dynamic ways to the evolutionary trajectories of a lineage.

Landscape Genomics to Enable Conservation Actions: The California Conservation Genomics Project.

The California Conservation Genomics Project (CCGP) is a unique, critically important step forward in the use of comprehensive landscape genetic data to modernize natural resource management at a regional scale. We describe the CCGP, including all aspects of project administration, data collection, current progress, and future challenges. The CCGP will generate, analyze, and curate a single high-quality reference genome and 100-150 resequenced genomes for each of 153 species projects (representing 235 individual species) that span the ecological and phylogenetic breadth of California's marine, freshwater, and terrestrial ecosystems. The resulting portfolio of roughly 20 000 resequenced genomes will be analyzed with identical informatic and landscape genomic pipelines, providing a comprehensive overview of hotspots of within-species genomic diversity, potential and realized corridors connecting these hotspots, regions of reduced diversity requiring genetic rescue, and the distribution of variation critical for rapid climate adaptation. After 2 years of concerted effort, full funding ($12M USD) has been secured, species identified, and funds distributed to 68 laboratories and 114 investigators drawn from all 10 University of California campuses. The remaining phases of the CCGP include completion of data collection and analyses, and delivery of the resulting genomic data and inferences to state and federal regulatory agencies to help stabilize species declines. The aspirational goals of the CCGP are to identify geographic regions that are critical to long-term preservation of California biodiversity, prioritize those regions based on defensible genomic criteria, and provide foundational knowledge that informs management strategies at both the individual species and ecosystem levels.

Landscape analyses using eDNA metabarcoding and Earth observation predict community biodiversity in California.

Ecosystems globally are under threat from ongoing anthropogenic environmental change. Effective conservation management requires more thorough biodiversity surveys that can reveal system-level patterns and that can be applied rapidly across space and time. Using modern ecological models and community science, we integrate environmental DNA and Earth observations to produce a time snapshot of regional biodiversity patterns and provide multi-scalar community-level characterization. We collected 278 samples in spring 2017 from coastal, shrub, and lowland forest sites in California, a complex ecosystem and biodiversity hotspot. We recovered 16,118 taxonomic entries from eDNA analyses and compiled associated traditional observations and environmental data to assess how well they predicted alpha, beta, and zeta diversity. We found that local habitat classification was diagnostic of community composition and distinct communities and organisms in different kingdoms are predicted by different environmental variables. Nonetheless, gradient forest models of 915 families recovered by eDNA analysis and using BIOCLIM variables, Sentinel-2 satellite data, human impact, and topographical features as predictors, explained 35% of the variance in community turnover. Elevation, sand percentage, and photosynthetic activities (NDVI32) were the top predictors. In addition to this signal of environmental filtering, we found a positive relationship between environmentally predicted families and their numbers of biotic interactions, suggesting environmental change could have a disproportionate effect on community networks. Together, these analyses show that coupling eDNA with environmental predictors including remote sensing data has capacity to test proposed Essential Biodiversity Variables and create new landscape biodiversity baselines that span the tree of life.

Precipitation and vegetation shape patterns of genomic and craniometric variation in the central African rodent Praomys misonnei.

Predicting species' capacity to respond to climate change is an essential first step in developing effective conservation strategies. However, conservation prioritization schemes rarely take evolutionary potential into account. Ecotones provide important opportunities for diversifying selection and may thus constitute reservoirs of standing variation, increasing the capacity for future adaptation. Here, we map patterns of environmentally associated genomic and craniometric variation in the central African rodent Praomys misonnei to identify areas with the greatest turnover in genomic composition. We also project patterns of environmentally associated genomic variation under future climate change scenarios to determine where populations may be under the greatest pressure to adapt. While precipitation gradients influence both genomic and craniometric variation, vegetation structure is also an important determinant of craniometric variation. Areas of elevated environmentally associated genomic and craniometric variation overlap with zones of rapid ecological transition underlining their importance as reservoirs of evolutionary potential. We also find that populations in the Sanaga river basin, central Cameroon and coastal Gabon are likely to be under the greatest pressure from climate change. Lastly, we make specific conservation recommendations on how to protect zones of high evolutionary potential and identify areas where populations may be the most susceptible to climate change.

A genoscape-network model for conservation prioritization in a migratory bird.

Migratory animals are declining worldwide and coordinated conservation efforts are needed to reverse current trends. We devised a novel genoscape-network model that combines genetic analyses with species distribution modeling and demographic data to overcome challenges with conceptualizing alternative risk factors in migratory species across their full annual cycle. We applied our method to the long distance, Neotropical migratory bird, Wilson's Warbler (Cardellina pusilla). Despite a lack of data from some wintering locations, we demonstrated how the results can be used to help prioritize conservation of breeding and wintering areas. For example, we showed that when genetic, demographic, and network modeling results were considered together it became clear that conservation recommendations will differ depending on whether the goal is to preserve unique genetic lineages or the largest number of birds per unit area. More specifically, if preservation of genetic lineages is the goal, then limited resources should be focused on preserving habitat in the California Sierra, Basin Rockies, or Coastal California, where the 3 most vulnerable genetic lineages breed, or in western Mexico, where 2 of the 3 most vulnerable lineages overwinter. Alternatively, if preservation of the largest number of individuals per unit area is the goal, then limited conservation dollars should be placed in the Pacific Northwest or Central America, where densities are estimated to be the highest. Overall, our results demonstrated the utility of adopting a genetically based network model for integrating multiple types of data across vast geographic scales and better inform conservation decision-making for migratory animals.

Un Modelo de Redes de Panorama Poblacional para la Priorización de la Conservación de un Ave Migratoria Resumen Los animales migratorios están pasando por una declinación mundial y se requieren esfuerzos coordinados de conservación para revertir las tendencias actuales. Diseñamos un modelo novedoso de redes de panorama poblacional que combina el análisis genético con el modelado de la distribución de especies y los datos demográficos para sobreponerse a los obstáculos con la conceptualización de los factores alternativos de riesgo en las especies migratorias durante su ciclo anual completo. Aplicamos nuestro método al chipe de corona negra (Cardellina pusilla), un ave migratoria neotropical que recorre largas distancias. A pesar de la falta de datos de algunas localidades de invernación, mostramos cómo pueden usarse los resultados para ayudar a priorizar la conservación de las áreas de reproducción y de invernación. Por ejemplo, mostramos que cuando se consideraron en conjunto los resultados del modelado genético, demográfico y de redes queda claro que las recomendaciones de conservación diferirán dependiendo de si el objetivo es preservar linajes genéticos únicos o el mayor número de aves por unidad de área. Más específicamente, si el objetivo es la conservación de los linajes genéticos, entonces los recursos limitados deberían enfocarse en preservar el hábitat en la Sierra de California, la Cuenca de las Rocallosas, la costa de California (lugares en donde se reproducen los tres linajes genéticos más vulnerables) o en el oeste de México (en donde dos de los tres linajes más vulnerables pasan el invierno). Alternativamente, si el objetivo es la conservación del mayor número de individuos por unidad de área, entonces el financiamiento limitado debería aplicarse en el noroeste del Pacífico o en América Central, en donde se estima que las densidades poblacionales son las más altas. En general, nuestros resultados demostraron la utilidad de adoptar un modelo de redes basadas en la genética para la integración de datos a lo largo de escalas geográficas amplias y para informar de mejor manera la toma de decisiones de conservación para los animales migratorios.

Natural re-colonization and admixture of wolves (Canis lupus) in the US Pacific Northwest: challenges for the protection and management of rare and endangered taxa.

Admixture resulting from natural dispersal processes can potentially generate novel phenotypic variation that may facilitate persistence in changing environments or result in the loss of population-specific adaptations. Yet, under the US Endangered Species Act, policy is limited for management of individuals whose ancestry includes a protected taxon; therefore, they are generally not protected under the Act. This issue is exemplified by the recently re-established grey wolves of the Pacific Northwest states of Washington and Oregon, USA. This population was likely founded by two phenotypically and genetically distinct wolf ecotypes: Northern Rocky Mountain (NRM) forest and coastal rainforest. The latter is considered potentially threatened in southeast Alaska and thus the source of migrants may affect plans for their protection. To assess the genetic source of the re-established population, we sequenced a ~ 300 bp portion of the mitochondrial control region and ~ 5 Mbp of the nuclear genome. Genetic analysis revealed that the Washington wolves share ancestry with both wolf ecotypes, whereas the Oregon population shares ancestry with NRM forest wolves only. Using ecological niche modelling, we found that the Pacific Northwest states contain environments suitable for each ecotype, with wolf packs established in both environmental types. Continued migration from coastal rainforest and NRM forest source populations may increase the genetic diversity of the Pacific Northwest population. However, this admixed population challenges traditional management regimes given that admixture occurs between an adaptively distinct ecotype and a more abundant reintroduced interior form. Our results emphasize the need for a more precise US policy to address the general problem of admixture in the management of endangered species, subspecies, and distinct population segments.

Ecological genomics predicts climate vulnerability in an endangered southwestern songbird.

Few regions have been more severely impacted by climate change in the USA than the Desert Southwest. Here, we use ecological genomics to assess the potential for adaptation to rising global temperatures in a widespread songbird, the willow flycatcher (Empidonax traillii), and find the endangered desert southwestern subspecies (E. t. extimus) most vulnerable to future climate change. Highly significant correlations between present abundance and estimates of genomic vulnerability - the mismatch between current and predicted future genotype-environment relationships - indicate small, fragmented populations of the southwestern willow flycatcher will have to adapt most to keep pace with climate change. Links between climate-associated genotypes and genes important to thermal tolerance in birds provide a potential mechanism for adaptation to temperature extremes. Our results demonstrate that the incorporation of genotype-environment relationships into landscape-scale models of climate vulnerability can facilitate more precise predictions of climate impacts and help guide conservation in threatened and endangered groups.

Hierarchical decision-making balances current and future reproductive success.

Parental decisions in animals are often context-dependent and shaped by fitness trade-offs between parents and offspring. For example, the selection of breeding habitats can considerably impact the fitness of both offspring and parents, and therefore, parents should carefully weigh the costs and benefits of available options for their current and future reproductive success. Here, we show that resource-use preferences are shaped by a trade-off between parental effort and offspring safety in a tadpole-transporting frog. In a large-scale in situ experiment, we investigated decision strategies across an entire population of poison frogs that distribute their tadpoles across multiple water bodies. Pool use followed a dynamic and sequential selection process, and transportation became more efficient over time. Our results point to a complex suite of environmental variables that are considered during offspring deposition, which necessitates a highly dynamic and flexible decision-making process in tadpole-transporting frogs.

Determining the drivers of population structure in a highly urbanized landscape to inform conservation planning.

Understanding the environmental contributors to population structure is of paramount importance for conservation in urbanized environments. We used spatially explicit models to determine genetic population structure under current and future environmental conditions across a highly fragmented, human-dominated environment in Southern California to assess the effects of natural ecological variation and urbanization. We focused on 7 common species with diverse habitat requirements, home-range sizes, and dispersal abilities. We quantified the relative roles of potential barriers, including natural environmental characteristics and an anthropogenic barrier created by a major highway, in shaping genetic variation. The ability to predict genetic variation in our models differed among species: 11-81% of intraspecific genetic variation was explained by environmental variables. Although an anthropogenically induced barrier (a major highway) severely restricted gene flow and movement at broad scales for some species, genetic variation seemed to be primarily driven by natural environmental heterogeneity at a local level. Our results show how assessing environmentally associated variation for multiple species under current and future climate conditions can help identify priority regions for maximizing population persistence under environmental change in urbanized regions.

Persistent impacts of West Nile virus on North American bird populations.

Since its introduction to North America in 1999, West Nile virus (WNV) has had devastating impacts on native host populations, but to date these impacts have been difficult to measure. Using a continental-scale dataset comprised of a quarter-million birds captured over nearly two decades and a recently developed model of WNV risk, we estimated the impact of this emergent disease on the survival of avian populations. We find that populations were negatively affected by WNV in 23 of the 49 species studied (47%). We distinguished two groups of species: those for which WNV negatively impacted survival only during initial spread of the disease (n = 11), and those that show no signs of recovery since disease introduction (n = 12). Results provide a novel example of the taxonomic breadth and persistent impacts of this wildlife disease on a continental scale. Phylogenetic analyses further identify groups (New World sparrows, finches, and vireos) disproportionally affected by temporary or persistent WNV effects, suggesting an evolutionary dimension of disease risk. Identifying the factors affecting the persistence of a disease across host species is critical to mitigating its effects, particularly in a world marked by rapid anthropogenic change.

Genomic divergence across ecological gradients in the Central African rainforest songbird (Andropadus virens).

The little greenbul, a common rainforest passerine from sub-Saharan Africa, has been the subject of long-term evolutionary studies to understand the mechanisms leading to rainforest speciation. Previous research found morphological and behavioural divergence across rainforest-savannah transition zones (ecotones), and a pattern of divergence with gene flow suggesting divergent natural selection has contributed to adaptive divergence and ecotones could be important areas for rainforests speciation. Recent advances in genomics and environmental modelling make it possible to examine patterns of genetic divergence in a more comprehensive fashion. To assess the extent to which natural selection may drive patterns of differentiation, here we investigate patterns of genomic differentiation among populations across environmental gradients and regions. We find compelling evidence that individuals form discrete genetic clusters corresponding to distinctive environmental characteristics and habitat types. Pairwise FST between populations in different habitats is significantly higher than within habitats, and this differentiation is greater than what is expected from geographic distance alone. Moreover, we identified 140 SNPs that showed extreme differentiation among populations through a genomewide selection scan. These outliers were significantly enriched in exonic and coding regions, suggesting their functional importance. Environmental association analysis of SNP variation indicates that several environmental variables, including temperature and elevation, play important roles in driving the pattern of genomic diversification. Results lend important new genomic evidence for environmental gradients being important in population differentiation.

A continental risk assessment of West Nile virus under climate change.

Since first introduced to North America in 1999, West Nile virus (WNV) has spread rapidly across the continent, threatening wildlife populations and posing serious health risks to humans. While WNV incidence has been linked to environmental factors, particularly temperature and rainfall, little is known about how future climate change may affect the spread of the disease. Using available data on WNV infections in vectors and hosts collected from 2003-2011 and using a suite of 10 species distribution models, weighted according to their predictive performance, we modeled the incidence of WNV under current climate conditions at a continental scale. Models were found to accurately predict spatial patterns of WNV that were then used to examine how future climate may affect the spread of the disease. Predictions were accurate for cases of human WNV infection in the following year (2012), with areas reporting infections having significantly higher probability of presence as predicted by our models. Projected geographic distributions of WNV in North America under future climate for 2050 and 2080 show an expansion of suitable climate for the disease, driven by warmer temperatures and lower annual precipitation that will result in the exposure of new and naïve host populations to the virus with potentially serious consequences. Our risk assessment identifies current and future hotspots of West Nile virus where mitigation efforts should be focused and presents an important new approach for monitoring vector-borne disease under climate change.

Spillover of pH1N1 to swine in Cameroon: an investigation of risk factors.

BACKGROUND: The 2009 pH1N1 influenza pandemic resulted in at least 18,500 deaths worldwide. While pH1N1 is now considered to be in a post-pandemic stage in humans it has nevertheless spilled back into swine in at least 20 countries. Understanding the factors that increase the risk of spillover events between swine and humans is essential to predicting and preventing future outbreaks. We assessed risk factors that may have led to spillover of pH1N1 from humans to swine in Cameroon, Central Africa. We sampled swine, domestic poultry and wild birds for influenza A virus at twelve sites in Cameroon from December 2009 while the pandemic was ongoing, to August 2012. At the same time we conducted point-count surveys to assess the abundance of domestic livestock and wild birds and assess interspecific contact rates. Random forest models were used to assess which variables were the best predictors of influenza in swine. RESULTS: We found swine with either active pH1N1 infections or positive for influenza A at four of our 12 sites. Only one swine tested positive by competitive ELISA in 2011-2012. To date we have found pH1N1 only in the North and Extreme North regions of Cameroon (regions in Cameroon are administrative units similar to provinces), though half of our sites are in the Central and Western regions. Swine husbandry practices differ between the North and Extreme North regions where it is common practice in to let swine roam freely, and the Central and Western regions where swine are typically confined to pens. Random forest analyses revealed that the three best predictors of the presence of pH1N1 in swine were contact rates between free-ranging swine and domestic ducks, contact rates between free-ranging swine and wild Columbiformes, and contact rates between humans and ducks. Sites in which swine were allowed to range freely had closer contact with other species than did sites in which swine were kept penned. CONCLUSIONS: Results suggest that the practice of allowing swine to roam freely is a significant risk factor for spillover of influenza from humans into swine populations.

Strong habitat-specific phenotypic plasticity but no genome-wide differentiation across a rainforest gradient in an African butterfly.

Habitat-specific thermal responses are well documented in various organisms and likely determine the vulnerability of populations to climate change. However, the underlying roles of genetics and plasticity that shape such habitat-specific patterns are rarely investigated together. Here we examined the thermal plasticity of the butterfly Bicyclus dorothea originating from rainforest and ecotone habitats in Cameroon under common garden conditions. We also sampled wild-caught butterflies from forest and ecotone sites and used RADseq to explore genome-wide population differentiation. We found differences in the level of phenotypic plasticity across habitats. Specifically, ecotone populations exhibited greater sensitivity in wing eyespot features with variable development temperatures relative to rainforest populations. Known adaptive roles of wing eyespots in Bicyclus species suggest that this morphological plasticity is likely under divergent selection across environmental gradients. However, we found no distinct population structure of genome-wide variation between habitats, suggesting high level of ongoing gene flow between habitats is homogenizing most parts of the genome.

Genetic identification of avian samples recovered from solar energy installations.

Renewable energy production and development will drastically affect how we meet global energy demands, while simultaneously reducing the impact of climate change. Although the possible effects of renewable energy production (mainly from solar- and wind-energy facilities) on wildlife have been explored, knowledge gaps still exist, and collecting data from wildlife remains (when negative interactions occur) at energy installations can act as a first step regarding the study of species and communities interacting with facilities. In the case of avian species, samples can be collected relatively easily (as compared to other sampling methods), but may only be able to be identified when morphological characteristics are diagnostic for a species. Therefore, many samples that appear as partial remains, or "feather spots"-known to be of avian origin but not readily assignable to species via morphology-may remain unidentified, reducing the efficiency of sample collection and the accuracy of patterns observed. To obtain data from these samples and ensure their identification and inclusion in subsequent analyses, we applied, for the first time, a DNA barcoding approach that uses mitochondrial genetic data to identify unknown avian samples collected at solar facilities to species. We also verified and compared identifications obtained by our genetic method to traditional morphological identifications using a blind test, and discuss discrepancies observed. Our results suggest that this genetic tool can be used to verify, correct, and supplement identifications made in the field and can produce data that allow accurate comparisons of avian interactions across facilities, locations, or technology types. We recommend implementing this genetic approach to ensure that unknown samples collected are efficiently identified and contribute to a better understanding of wildlife impacts at renewable energy projects.

Genomic analyses reveal poaching hotspots and illegal trade in pangolins from Africa to Asia.

The white-bellied pangolin (Phataginus tricuspis) is the world's most trafficked mammal and is at risk of extinction. Reducing the illegal wildlife trade requires an understanding of its origins. Using a genomic approach for tracing confiscations and analyzing 111 samples collected from known geographic localities in Africa and 643 seized scales from Asia between 2012 and 2018, we found that poaching pressures shifted over time from West to Central Africa. Recently, Cameroon's southern border has emerged as a site of intense poaching. Using data from seizures representing nearly 1 million African pangolins, we identified Nigeria as one important hub for trafficking, where scales are amassed and transshipped to markets in Asia. This origin-to-destination approach offers new opportunities to disrupt the illegal wildlife trade and to guide anti-trafficking measures.

Host and habitat specialization of avian malaria in Africa.

Studies of both vertebrates and invertebrates have suggested that specialists, as compared to generalists, are likely to suffer more serious declines in response to environmental change. Less is known about the effects of environmental conditions on specialist versus generalist parasites. Here, we study the evolutionary strategies of malaria parasites (Plasmodium spp.) among different bird host communities. We determined the parasite diversity and prevalence of avian malaria in three bird communities in the lowland forests in Cameroon, highland forests in East Africa and fynbos in South Africa. We calculated the host specificity index of parasites to examine the range of hosts parasitized as a function of the habitat and investigated the phylogenetic relationships of parasites. First, using phylogenetic and ancestral reconstruction analyses, we found an evolutionary tendency for generalist malaria parasites to become specialists. The transition rate at which generalists become specialists was nearly four times as great as the rate at which specialists become generalists. We also found more specialist parasites and greater parasite diversity in African lowland rainforests as compared to the more climatically variable habitats of the fynbos and the highland forests. Thus, with environmental changes, we anticipate a change in the distribution of both specialist and generalist parasites with potential impacts on bird communities.