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.

Sex-dependent effects of acute stress on amyloid-ß in male and female mice.

The risk of developing Alzheimer's disease is mediated by a combination of genetics and environmental factors, such as stress, sleep abnormalities and traumatic brain injury. Women are at a higher risk of developing Alzheimer's disease than men, even when controlling for differences in lifespan. Women are also more likely to report high levels of stress than men. Sex differences in response to stress may play a role in the increased risk of Alzheimer's disease in women. In this study, we use in vivo microdialysis to measure levels of Aß in response to acute stress in male and female mice. We show that Aß levels are altered differently between female and male mice (APP/PS1 and wild-type) in response to stress, with females showing significantly increased levels of Aß while most males do not show a significant change. This response is mediated through ß-arrestin involvement in Corticotrophin Releasing Factor receptor signalling pathway differences in male and female mice as male mice lacking ß-arrestin show increase in Aß in response to stress similar to females.

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.

Population structure, inbreeding and stripe pattern abnormalities in plains zebras.

One of the most iconic wild equids, the plains zebra occupies a broad region of sub-Saharan Africa and exhibits a wide range of phenotypic diversity in stripe patterns that have been used to classify multiple subspecies. After decades of relative stability, albeit with a loss of at least one recognized subspecies, the total population of plains zebras has undergone an approximate 25% decline since 2002. Individuals with abnormal stripe patterns have been recognized in recent years but the extent to which their appearance is related to demography and/or genetics is unclear. Investigating population genetic health and genetic structure are essential for developing effective strategies for plains zebra conservation. We collected DNA from 140 plains zebra, including seven with abnormal stripe patterns, from nine locations across the range of plains zebra, and analyzed data from restriction site-associated and whole genome sequencing (RAD-seq, WGS) libraries to better understand the relationships between population structure, genetic diversity, inbreeding, and abnormal phenotypes. We found that genetic structure did not coincide with described subspecific variation, but did distinguish geographic regions in which anthropogenic habitat fragmentation is associated with reduced gene flow and increased evidence of inbreeding, especially in certain parts of East Africa. Further, zebras with abnormal striping exhibited increased levels of inbreeding relative to normally striped individuals from the same populations. Our results point to a genetic cause of stripe pattern abnormalities, and dramatic evidence of the consequences of habitat fragmentation.

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.

Targeted capture and resequencing of 1040 genes reveal environmentally driven functional variation in grey wolves.

In an era of ever-increasing amounts of whole-genome sequence data for individuals and populations, the utility of traditional single nucleotide polymorphisms (SNPs) array-based genome scans is uncertain. We previously performed a SNP array-based genome scan to identify candidate genes under selection in six distinct grey wolf (Canis lupus) ecotypes. Using this information, we designed a targeted capture array for 1040 genes, including all exons and flanking regions, as well as 5000 1-kb nongenic neutral regions, and resequenced these regions in 107 wolves. Selection tests revealed striking patterns of variation within candidate genes relative to noncandidate regions and identified potentially functional variants related to local adaptation. We found 27% and 47% of candidate genes from the previous SNP array study had functional changes that were outliers in sweed and bayenv analyses, respectively. This result verifies the use of genomewide SNP surveys to tag genes that contain functional variants between populations. We highlight nonsynonymous variants in APOB, LIPG and USH2A that occur in functional domains of these proteins, and that demonstrate high correlation with precipitation seasonality and vegetation. We find Arctic and High Arctic wolf ecotypes have higher numbers of genes under selection, which highlight their conservation value and heightened threat due to climate change. This study demonstrates that combining genomewide genotyping arrays with large-scale resequencing and environmental data provides a powerful approach to discern candidate functional variants in natural populations.

Genetic subdivision and candidate genes under selection in North American grey wolves.

Previous genetic studies of the highly mobile grey wolf (Canis lupus) found population structure that coincides with habitat and phenotype differences. We hypothesized that these ecologically distinct populations (ecotypes) should exhibit signatures of selection in genes related to morphology, coat colour and metabolism. To test these predictions, we quantified population structure related to habitat using a genotyping array to assess variation in 42 036 single-nucleotide polymorphisms (SNPs) in 111 North American grey wolves. Using these SNP data and individual-level measurements of 12 environmental variables, we identified six ecotypes: West Forest, Boreal Forest, Arctic, High Arctic, British Columbia and Atlantic Forest. Next, we explored signals of selection across these wolf ecotypes through the use of three complementary methods to detect selection: FST /haplotype homozygosity bivariate percentilae, bayescan, and environmentally correlated directional selection with bayenv. Across all methods, we found consistent signals of selection on genes related to morphology, coat coloration, metabolism, as predicted, as well as vision and hearing. In several high-ranking candidate genes, including LEPR, TYR and SLC14A2, we found variation in allele frequencies that follow environmental changes in temperature and precipitation, a result that is consistent with local adaptation rather than genetic drift. Our findings show that local adaptation can occur despite gene flow in a highly mobile species and can be detected through a moderately dense genomic scan. These patterns of local adaptation revealed by SNP genotyping likely reflect high fidelity to natal habitats of dispersing wolves, strong ecological divergence among habitats, and moderate levels of linkage in the wolf genome.

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.