Amanda J. Lea

Genome-wide Evidence Reveals that African and Eurasian Golden Jackals Are Distinct Species

The golden jackal of Africa (Canis aureus) has long been considered a conspecific of jackals distributed throughout Eurasia, with the nearest source populations in the Middle East. However, two recent reports found that mitochondrial haplotypes of some African golden jackals aligned more closely to gray wolves (Canis lupus), which is surprising given the absence of gray wolves in Africa and the phenotypic divergence between the two species. Moreover, these results imply the existence of a previously unrecognized phylogenetically distinct species despite a long history of taxonomic work on African canids. To test the distinct-species hypothesis and understand the evolutionary history that would account for this puzzling result, we analyzed extensive genomic data including mitochondrial genome sequences, sequences from 20 autosomal loci (17 introns and 3 exon segments), microsatellite loci, X- and Y-linked zinc-finger protein gene (ZFX and ZFY) sequences, and whole-genome nuclear sequences in African and Eurasian golden jackals and gray wolves. Our results provide consistent and robust evidence that populations of golden jackals from Africa and Eurasia represent distinct monophyletic lineages separated for more than one million years, sufficient to merit formal recognition as different species: C. anthus (African golden wolf) and C. aureus (Eurasian golden jackal). Using morphologic data, we demonstrate a striking morphologic similarity between East African and Eurasian golden jackals, suggesting parallelism, which may have misled taxonomists and likely reflects uniquely intense interspecific competition in the East African carnivore guild. Our study shows how ecology can confound taxonomy if interspecific competition constrains size diversification.

Resource base influences genome-wide DNA methylation levels in wild baboons (Papio cynocephalus)

Variation in resource availability commonly exerts strong effects on fitness‐related traits in wild animals. However, we know little about the molecular mechanisms that mediate these effects, or about their persistence over time. To address these questions, we profiled genome‐wide whole‐blood DNA methylation levels in two sets of wild baboons: (i) ‘wild‐feeding’ baboons that foraged naturally in a savanna environment and (ii) ‘Lodge’ baboons that had ready access to spatially concentrated human food scraps, resulting in high feeding efficiency and low daily travel distances. We identified 1014 sites (0.20% of sites tested) that were differentially methylated between wild‐feeding and Lodge baboons, providing the first evidence that resource availability shapes the epigenome in a wild mammal. Differentially methylated sites tended to occur in contiguous stretches (i.e., in differentially methylated regions or DMRs), in promoters and enhancers, and near metabolism‐related genes, supporting their functional importance in gene regulation. In agreement, reporter assay experiments confirmed that methylation at the largest identified DMR, located in the promoter of a key glycolysis‐related gene, was sufficient to causally drive changes in gene expression. Intriguingly, all dispersing males carried a consistent epigenetic signature of their membership in a wild‐feeding group, regardless of whether males dispersed into or out of this group as adults. Together, our findings support a role for DNA methylation in mediating ecological effects on phenotypic traits in the wild and emphasize the dynamic environmental sensitivity of DNA methylation levels across the life course.

A Flexible, Efficient Binomial Mixed Model for Identifying Differential DNA Methylation in Bisulfite Sequencing Data

Identifying sources of variation in DNA methylation levels is important for understanding gene regulation. Recently, bisulfite sequencing has become a popular tool for investigating DNA methylation levels. However, modeling bisulfite sequencing data is complicated by dramatic variation in coverage across sites and individual samples, and because of the computational challenges of controlling for genetic covariance in count data. To address these challenges, we present a binomial mixed model and an efficient, sampling-based algorithm (MACAU: Mixed model association for count data via data augmentation) for approximate parameter estimation and p-value computation. This framework allows us to simultaneously account for both the over-dispersed, count-based nature of bisulfite sequencing data, as well as genetic relatedness among individuals. Using simulations and two real data sets (whole genome bisulfite sequencing (WGBS) data from Arabidopsis thaliana and reduced representation bisulfite sequencing (RRBS) data from baboons), we show that our method provides well-calibrated test statistics in the presence of population structure. Further, it improves power to detect differentially methylated sites: in the RRBS data set, MACAU detected 1.6-fold more age-associated CpG sites than a beta-binomial model (the next best approach). Changes in these sites are consistent with known age-related shifts in DNA methylation levels, and are enriched near genes that are differentially expressed with age in the same population. Taken together, our results indicate that MACAU is an efficient, effective tool for analyzing bisulfite sequencing data, with particular salience to analyses of structured populations. MACAU is freely available at www.xzlab.org/software.html.

Developmental Constraints in a Wild Primate

Early-life experiences can dramatically affect adult traits. However, the evolutionary origins of such early-life effects are debated. The predictive adaptive response hypothesis argues that adverse early environments prompt adaptive phenotypic adjustments that prepare animals for similar challenges in adulthood. In contrast, the developmental constraints hypothesis argues that early adversity is generally costly. To differentiate between these hypotheses, we studied two sets of wild female baboons: those born during low-rainfall, low-quality years and those born during normal-rainfall, high-quality years. For each female, we measured fertility-related fitness components during years in adulthood that matched and mismatched her early conditions. We found support for the developmental constraints hypothesis: females born in low-quality environments showed greater decreases in fertility during drought years than females born in high-quality environments, even though drought years matched the early conditions of females born in low-quality environments. Additionally, we found that females born in low-quality years to high-status mothers did not experience reduced fertility during drought years. These results indicate that early ecological adversity did not prepare individuals to cope with ecological challenges in later life. Instead, individuals that experienced at least one high-quality early environment—either ecological or social—were more resilient to ecological stress in later life. Together, these data suggest that early adversity carries lifelong costs, which is consistent with the developmental constraints hypothesis.

Disease and freeways drive genetic change in urban bobcat populations

Urbanization profoundly impacts animal populations by causing isolation, increased susceptibility to disease, and exposure to toxicants. Genetic effects include reduced effective population size, increased population substructure, and decreased adaptive potential. We investigated the influence that urbanization and a disease epizootic had on the population genetics of bobcats (Lynx rufus) distributed across a highly fragmented urban landscape. We genotyped more than 300 bobcats, sampled from 1996 to 2012, for variation at nine neutral and seven immune gene‐linked microsatellite loci. We found that two freeways are significant barriers to gene flow. Further, a 3‐year disease epizootic, associated with secondary anticoagulant rodenticide exposure, caused a population bottleneck that led to significant genetic differentiation between pre‐ and post‐disease populations that was greater than that between populations separated by major freeways for >60 years. However, balancing selection acted on immune‐linked loci during the epizootic, maintaining variation at functional regions. Conservation assessments need to assay loci that are potentially under selection to better preserve the adaptive potential of populations at the urban–wildland interface. Further, interconnected regions that contain appropriate habitat for wildlife will be critical to the long‐term viability of animal populations in urban landscapes.

Developmental plasticityBridging research in evolution and human health

Abstract Early life experiences can have profound and persistent effects on traits expressed throughout the life course, with consequences for later life behavior, disease risk, and mortality rates. The shaping of later life traits by early life environments, known as ‘developmental plasticity’, has been well-documented in humans and non-human animals, and has consequently captured the attention of both evolutionary biologists and researchers studying human health. Importantly, the parallel significance of developmental plasticity across multiple fields presents a timely opportunity to build a comprehensive understanding of this phenomenon. We aim to facilitate this goal by highlighting key outstanding questions shared by both evolutionary and health researchers, and by identifying theory and empirical work from both research traditions that is designed to address these questions. Specifically, we focus on: (i) evolutionary explanations for developmental plasticity, (ii) the genetics of developmental plasticity and (iii) the molecular mechanisms that mediate developmental plasticity. In each section, we emphasize the conceptual gains in human health and evolutionary biology that would follow from filling current knowledge gaps using interdisciplinary approaches. We encourage researchers interested in developmental plasticity to evaluate their own work in light of research from diverse fields, with the ultimate goal of establishing a cross-disciplinary understanding of developmental plasticity.

Urbanization and anticoagulant poisons promote immune dysfunction in bobcats

Understanding how human activities influence immune response to environmental stressors can support biodiversity conservation across increasingly urbanizing landscapes. We studied a bobcat (Lynx rufus) population in urban southern California that experienced a rapid population decline from 2002–2005 due to notoedric mange. Because anticoagulant rodenticide (AR) exposure was an underlying complication in mange deaths, we aimed to understand sublethal contributions of urbanization and ARs on 65 biochemical markers of immune and organ function. Variance in immunological variables was primarily associated with AR exposure and secondarily with urbanization. Use of urban habitat and AR exposure has pervasive, complex and predictable effects on biochemical markers of immune and organ function in free-ranging bobcats that include impacts on neutrophil, lymphocyte and cytokine populations, total bilirubin and phosphorus. We find evidence of both inflammatory response and immune suppression associated with urban land use and rat poison exposure that could influence susceptibility to opportunistic infections. Consequently, AR exposure may influence mortality and has population-level effects, as previous work in the focal population has revealed substantial mortality caused by mange infection. The secondary effects of anticoagulant exposure may be a worldwide, largely unrecognized problem affecting a variety of vertebrate species in human-dominated environments.

Maximizing ecological and evolutionary insight in bisulfite sequencing data sets

Genome-scale bisulfite sequencing approaches have opened the door to ecological and evolutionary studies of DNA methylation in many organisms. These approaches can be powerful. However, they introduce new methodological and statistical considerations, some of which are particularly relevant to non-model systems. Here, we highlight how these considerations influence a study’s power to link methylation variation with a predictor variable of interest. Relative to current practice, we argue that sample sizes will need to increase to provide robust insights. We also provide recommendations for overcoming common challenges and an R Shiny app to aid in study design.Bisulfite sequencing is widely used to study genome-scale DNA methylation. In this Review, the authors discuss methodological and statistical considerations related to bisulfite sequencing that are particularly relevant when studying non-model organisms.

Developmental plasticity research in evolution and human healthResponse to commentaries

Department of Biology, Duke University, Durham, NC 27708, USA, Institute of Primate Research, National Museums of Kenya, Karen, Nairobi, Kenya, Duke University Population Research Institute, Duke University, Durham, NC 27708, USA, Department of Evolutionary Anthropology, Duke University, Durham, NC 27708, USA, Department of Biological Sciences, University of Notre Dame, Notre Dame, IN 46556, USA, and 6 Present address: Lewis-Sigler Institute for Integrative Genomics, Carl Icahn Laboratory, Washington Road, Princeton University, Princeton, NJ 08544, USA *Corresponding author. Departments of Biology and Evolutionary Anthropology, Duke University, Box 90338, Durham, NC 27708, USA. Tel: 919-660-7272; Fax: 919-660-7293; E-mail: alberts@duke.edu Received 21 February 2018; revised version accepted 21 February 2018

Dominance rank-associated immune gene expression is widespread, sex-specific, and a precursor to high social status in wild male baboons

In humans and other hierarchical species, social status is tightly linked to variation in health and fitness-related traits. Experimental manipulations of social status in female rhesus macaques suggest that this relationship is partially explained by status effects on immune gene regulation. However, social hierarchies are established and maintained in different ways across species: while some are based on kin-directed nepotism, others emerge from direct physical competition. We investigated how this variation influences the relationship between social status and immune gene regulation in wild baboons, where hierarchies in males are based on fighting ability but female hierarchies are nepotistic. We measured rank-related variation in gene expression levels in adult baboons of both sexes at baseline and in response to ex vivo stimulation with the bacterial endotoxin lipopolysaccharide (LPS). We identified >2000 rank- associated genes in males, an order of magnitude more than in females. In males, high status predicted increased expression of genes involved in innate immunity and preferential activation of the NFkB-mediated pro-inflammatory pathway, a pattern previously associated with low status in female rhesus macaques. Using Mendelian randomization, we reconcile these observations by demonstrating that high status-associated gene expression patterns are precursors, not consequences, of high social status in males, in support of the idea that physiological condition determines who attains high rank. Together, our work provides the first test of the relationship between social status and immune gene regulation in wild primates. It also emphasizes the importance of social context in shaping the relationship between social status and immune function. SIGNIFICANCE Social status predicts fitness outcomes in social animals, motivating efforts to understand its physiological causes and consequences. We investigated the relationship between social status and immune gene expression in wild baboons, where female status is determined by kinship but male status is determined by fighting ability. We uncover pervasive status-gene expression associations in males, but not females. High status males exhibit high levels of pro-inflammatory gene expression, in contrast to previous findings in hierarchies that are not competitively determined. Using Mendelian randomization, we show that this status-associated variation precedes dominance rank attainment: males who compete successfully for high status are already immunologically distinct. The nature of social hierarchies thus fundamentally shapes the relationship between social status and immune function.