Search

Hybridization in nature offers unique insights into the process of natural selection in incipient species and their hybrids. In order to evaluate the patterns and targets of selection, we examine a recently discovered baboon hybrid zone in the Kafue River Valley of Zambia, where Kinda baboons (Papio kindae) and grey-footed chacma baboons (P. ursinus griseipes) coexist with hybridization. We genotyped baboons at 14,962 variable genome-wide autosomal markers using double-digest RADseq. We compared ancestry patterns from this genome-wide data set to previously reported ancestry from mitochondrial-DNA and Y-chromosome sources. We also fit a Bayesian genomic cline model to scan for genes with extreme patterns of introgression. We show that the Kinda baboon Y chromosome has penetrated the species boundary to a greater extent than either mitochondrial DNA or the autosomal chromosomes. We also find evidence for overall restricted introgression in the JAK/STAT signalling pathway. Echoing results in other species including humans, we find evidence for enhanced and/or directional introgression of immune-related genes or pathways including the toll-like receptor pathway, the blood coagulation pathway, and the LY96 gene. Finally we show enhanced introgression and excess chacma baboon ancestry in the sperm tail gene ODF2. Together, our results elucidate the dynamics of introgressive hybridization in a primate system while identifying genes and pathways possibly under selection.

Human culture, biology, and health were shaped dramatically by the onset of agriculture ∼12,000 y B.P. This shift is hypothesized to have resulted in increased individual fitness and population growth as evidenced by archaeological and population genomic data alongside a decline in physiological health as inferred from skeletal remains. Here, we consider osteological and ancient DNA data from the same prehistoric individuals to study human stature variation as a proxy for health across a transition to agriculture. Specifically, we compared “predicted” genetic contributions to height from paleogenomic data and “achieved” adult osteological height estimated from long bone measurements for 167 individuals across Europe spanning the Upper Paleolithic to Iron Age (∼38,000 to 2,400 B.P.). We found that individuals from the Neolithic were shorter than expected (given their individual polygenic height scores) by an average of −3.82 cm relative to individuals from the Upper Paleolithic and Mesolithic (P = 0.040) and −2.21 cm shorter relative to post-Neolithic individuals (P = 0.068), with osteological vs. expected stature steadily increasing across the Copper (+1.95 cm relative to the Neolithic), Bronze (+2.70 cm), and Iron (+3.27 cm) Ages. These results were attenuated when we additionally accounted for genome-wide genetic ancestry variation: for example, with Neolithic individuals −2.82 cm shorter than expected on average relative to pre-Neolithic individuals (P = 0.120). We also incorporated observations of paleopathological indicators of nonspecific stress that can persist from childhood to adulthood in skeletal remains into our model. Overall, our work highlights the potential of integrating disparate datasets to explore proxies of health in prehistory.

Dietary changes associated with industrialization increase the prevalence of chronic diseases, such as obesity, type II diabetes, and cardiovascular disease. This relationship is often attributed to an ‘evolutionary mismatch’ between human physiology and modern nutritional environments. Western diets enriched with foods that were scarce throughout human evolutionary history (e.g. simple sugars and saturated fats) promote inflammation and disease relative to diets more akin to ancestral human hunter-gatherer diets, such as a Mediterranean diet. Peripheral blood monocytes, precursors to macrophages and important mediators of innate immunity and inflammation, are sensitive to the environment and may represent a critical intermediate in the pathway linking diet to disease. We evaluated the effects of 15 months of whole diet manipulations mimicking Western or Mediterranean diet patterns on monocyte polarization in a well-established model of human health, the cynomolgus macaque (Macaca fascicularis). Monocyte transcriptional profiles differed markedly between diets, with 40% of transcripts showing differential expression (FDR < 0.05). Monocytes from Western diet consumers were polarized toward a more proinflammatory phenotype. The Western diet shifted the co-expression of 445 gene pairs, including small RNAs and transcription factors associated with metabolism and adiposity in humans, and dramatically altered behavior. For example, Western-fed individuals were more anxious and less socially integrated. These behavioral changes were also associated with some of the effects of diet on gene expression, suggesting an interaction between diet, central nervous system activity, and monocyte gene expression. This study provides new molecular insights into an evolutionary mismatch and uncovers new pathways through which Western diets alter monocyte polarization toward a proinflammatory phenotype.

Survival at high altitude requires adapting to extreme conditions such as environmental hypoxia. To understand high-altitude adaptations in a primate, we assembled the genome of the gelada (Theropithecus gelada), an endemic Ethiopian monkey, and complemented it with population resequencing, hematological, and morphometric data. Unexpectedly, we identified a novel karyotype that may contribute to reproductive isolation between gelada populations. We also identified genomic elements including protein-coding sequences and gene families that exhibit accelerated changes in geladas and may contribute to high-altitude adaptation. Our findings lend insight into mechanisms of speciation and adaptation while providing promising avenues for functional hypoxia research.

Anopheles mosquitoes are the principal vectors for malaria and lymphatic filariasis, and evidence for arboviral transmission under laboratory and natural contexts has been demonstrated. Vector management approaches require an understanding of the ecological, epidemiological, and biological context of the species in question, and increased interest in gene drive systems for vector control applications has resulted in an increased need for genome assemblies from understudied mosquito species. While the genomes for many Anopheles species have been sequenced, North American Anopheles taxa have been neglected. In this study we present novel genome assemblies for the North American species An. crucians, An. freeborni, An. albimanus, and An. quadrimaculatus, and examine the evolutionary relationship between these species. We identified 790 shared single copy orthologs between the newly sequenced genomes and created a phylogeny using 673 of the orthologs, identifying 525 orthologs with evidence for positive selection on at least one branch of the phylogeny. Gene ontology terms such as calcium ion signaling, histone binding, and protein acetylation were identified as being biased in the set of selected genes. These novel genome sequences will be useful in developing our understanding of the diverse biological traits that drive vectorial capacity in anophelines.

Human culture, biology, and health were shaped dramatically by the onset of agriculture ~12,000 years before present (BP). Subsistence shifts from hunting and gathering to agriculture are hypothesized to have resulted in increased individual fitness and population growth as evidenced by archaeological and population genomic data alongside a simultaneous decline in physiological health as inferred from paleopathological analyses and stature reconstructions of skeletal remains. A key component of the health decline inference is that relatively shorter statures observed for early farmers may (at least partly) reflect higher childhood disease burdens and poorer nutrition. However, while such stresses can indeed result in growth stunting, height is also highly heritable, and substantial inter-individual variation in the height genetic component within a population is typical. Moreover, extensive migration and gene flow were characteristics of multiple agricultural transitions worldwide. Here, we consider both osteological and ancient DNA data from the same prehistoric individuals to comprehensively study the trajectory of human stature variation as a proxy for health across a transition to agriculture. Specifically, we compared ‘predicted’ genetic contributions to height from paleogenomic data and ‘achieved’ adult osteological height estimated from long bone measurements on a per-individual basis for n=160 ancient Europeans from sites spanning the Upper Paleolithic to the Iron Age (~38,000-2,400 BP). We found that individuals from the Neolithic were shorter than expected (given their individual polygenic height scores) by an average of −4.47 cm relative to individuals from the Upper Paleolithic and Mesolithic (P=0.016). The average osteological vs. expected stature then increased relative to the Neolithic over the Copper (+2.67 cm, P=0.052), Bronze (+3.33 cm, P=0.032), and Iron Ages (+3.95 cm, P=0.094). These results were partly attenuated when we accounted for genome-wide genetic ancestry variation in our sample (which we note is partly duplicative with the individual polygenic score information). For example, in this secondary analysis Neolithic individuals were −3.48 cm shorter than expected on average relative to individuals from the Upper Paleolithic and Mesolithic (P=0.056). We also incorporated observations of paleopathological indicators of non-specific stress that can persist from childhood to adulthood in skeletal remains (linear enamel hypoplasia, cribra orbitalia, and porotic hyperostosis) into our model. Overall, our work highlights the potential of integrating disparate datasets to explore proxies of health in prehistory.

Malaria is infamous for the massive toll it exacts on human life and health. In the face of this intense selection, many human populations have evolved mechanisms that confer some resistance to the disease, such as sickle-cell hemoglobin or the Duffy null blood group. Less understood are adaptations in other vertebrate hosts, many of which have a longer history of co-evolution with malaria parasites. By comparing malaria resistance adaptations across host species, we can gain fundamental insight into host-parasite co-evolution. In particular, understanding the mechanisms by which non-human primate immune systems combat malaria may be fruitful in uncovering transferable therapeutic targets for humans. However, most research on primate response to malaria has focused on a single or few loci, typically in experimentally-infected captive primates. Here, we report the first transcriptomic study of a wild primate response to a malaria-like parasite, investigating gene expression response of red colobus monkeys (Piliocolobus tephrosceles) to natural infection with the malaria-like parasite, Hepatocystis. We identified colobus genes with expression strongly correlated with parasitemia, including many implicated in human malaria and suggestive of common genetic architecture of disease response. For instance, the expression of ACKR1 (alias DARC) gene, previously linked to resistance in humans, was found to be positively correlated with parasitemia. Other similarities to human parasite response include induction of changes in immune cell type composition and, potentially, increased extramedullary hematopoiesis and altered biosynthesis of neutral lipids. Our results illustrate the utility of comparative immunogenetic investigation of malaria response in primates. Such inter-specific comparisons of transcriptional response to pathogens afford a unique opportunity to compare and contrast the adaptive genetic architecture of disease resistance, which may lead to the identification of novel intervention targets to improve human health.Author SummaryThe co-evolutionary arms race between humans and malaria parasites has been ongoing for millennia. Fully understanding the evolved human response to malaria is impossible without comparative study of parasites in our non-human primate relatives. Though laboratory primates are fruitful models, the complexity of wild primates infected in a natural transmission system may be a more suitable comparison for contextualizing malaria infections in human patients. Here, we investigate the genetic mechanisms underlying the immune response to Hepatocystis, a close relative of human-infective malaria, in a population of wild Ugandan red colobus monkeys. We find that the genes involved have considerable overlap with those active in human malaria patients. Like Plasmodium, Hepatocystis induces changes in blood cell type and may cause the host to produce blood components outside of the bone marrow or alter metabolism related to the production of lipids. Our work helps to identify the genetic mechanisms underlying the arms race between primates and malaria parasites, providing fundamental evolutionary insight. Such comparative work on the interaction between wild non-human primates and malaria parasites can identify ways in which primates have evolved resistance to malaria parasites, and further investigation of such implicated genes may lead to novel potential therapeutic and vaccine targets.

Mayaro virus (MAYV) is an arboviral pathogen in the genus Alphavirus that is circulating in South America with potential to spread to naïve regions. MAYV is also one of the few viruses with the ability to be transmitted by mosquitoes in the genus Anopheles, as well as the typical arboviral transmitting mosquitoes in the genus Aedes. Few studies have investigated the infection response of Anopheles mosquitoes. In this study we detail the transcriptomic and small RNA responses of An. stephensi to infection with MAYV via infectious bloodmeal at 2, 7, and 14 days post infection (dpi). 487 unique transcripts were significantly regulated, 78 putative novel miRNAs were identified, and an siRNA response is observed targeting the MAYV genome. Gene ontology analysis of transcripts regulated at each timepoint suggested activation of the Toll pathway at 7 dpi and repression of pathways related to autophagy and apoptosis at 14 dpi. These findings provide a basic understanding of the infection response of An. stephensi to MAYV and help to identify host factors which might be useful to target to inhibit viral replication in Anopheles mosquitoes.AUTHOR SUMMARYMayaro virus (MAYV) is a mosquito-borne Alphavirus responsible for outbreaks in South America and the Caribbean. In this study we infected Anopheles stephensi with MAYV and sequenced mRNA and small RNA to understand how MAYV infection impacts gene transcription and the expression of small RNAs in the mosquito vector. Genes involved with innate immunity and signaling pathways related to cell death are regulated in response to MAYV infection of An. stephensi, we also discover novel miRNAs and describe the expression patterns of miRNAs, siRNAs, and piRNAs following bloodmeal ingestion. These results suggest that MAYV does induce a molecular response to infection in its mosquito vector species.

Significance Based on “subfossil” skeletal remains it is known that multiple now-extinct giant lemur (primate) species with estimated body masses of up to ∼160 kg survived on Madagascar into the past millennium. In this study, we used ancient DNA methods to sequence the nuclear genome of one of these megafaunal lemurs, Megaladapis edwardsi (∼85 kg). With the power of the nuclear genome, we robustly resolved the phylogenetic relationship between Megaladapis and other lemurs, which had been a lingering uncertainty. We also identified multiple signatures of past positive natural selection across the M. edwardsi genome that support reconstructions of this taxon as a large-bodied, specialized folivore., No endemic Madagascar animal with body mass >10 kg survived a relatively recent wave of extinction on the island. From morphological and isotopic analyses of skeletal “subfossil” remains we can reconstruct some of the biology and behavioral ecology of giant lemurs (primates; up to ∼160 kg) and other extraordinary Malagasy megafauna that survived into the past millennium. Yet, much about the evolutionary biology of these now-extinct species remains unknown, along with persistent phylogenetic uncertainty in some cases. Thankfully, despite the challenges of DNA preservation in tropical and subtropical environments, technical advances have enabled the recovery of ancient DNA from some Malagasy subfossil specimens. Here, we present a nuclear genome sequence (∼2× coverage) for one of the largest extinct lemurs, the koala lemur Megaladapis edwardsi (∼85 kg). To support the testing of key phylogenetic and evolutionary hypotheses, we also generated high-coverage nuclear genomes for two extant lemurs, Eulemur rufifrons and Lepilemur mustelinus, and we aligned these sequences with previously published genomes for three other extant lemurs and 47 nonlemur vertebrates. Our phylogenetic results confirm that Megaladapis is most closely related to the extant Lemuridae (typified in our analysis by E. rufifrons) to the exclusion of L. mustelinus, which contradicts morphology-based phylogenies. Our evolutionary analyses identified significant convergent evolution between M. edwardsi and an extant folivore (a colobine monkey) and an herbivore (horse) in genes encoding proteins that function in plant toxin biodegradation and nutrient absorption. These results suggest that koala lemurs were highly adapted to a leaf-based diet, which may also explain their convergent craniodental morphology with the small-bodied folivore Lepilemur.

The COVID-19 pandemic, caused by the coronavirus SARS-CoV-2, has devastated health infrastructure around the world. Both ACE2 (an entry receptor) and TMPRSS2 (used by the virus for spike protein priming) are key proteins to SARS-CoV-2 cell entry, enabling progression to COVID-19 in humans. Comparative genomic research into critical ACE2 binding sites, associated with the spike receptor binding domain, has suggested that Old World primates may also be susceptible to disease from SARS-CoV-2 infection. Savanna monkeys (Chlorocebus spp.) are a widespread non-human primate with well-established potential as a bi-directional zoonotic/anthroponotic agent due to high levels of human interaction throughout their range in sub-Saharan Africa and the Caribbean. To characterize potential functional variation in savanna monkey ACE2 and TMPRSS2, we inspected recently published genomic data from 245 savanna monkeys, including 163 wild monkeys from Africa and the Caribbean and 82 captive monkeys from the Vervet Research Colony (VRC). We found several missense variants. One missense variant in ACE2 (X:14,077,550; Asp30Gly), common in Ch. sabaeus, causes a change in amino acid residue that has been inferred to reduce binding efficiency of SARS-CoV-2, suggesting potentially reduced susceptibility. The remaining populations appear as susceptible as humans, based on these criteria for receptor usage. All missense variants observed in wild Ch. sabaeus populations are also present in the VRC, along with two splice acceptor variants (at X:14,065,076) not observed in the wild sample that are potentially disruptive to ACE2 function. The presence of these variants in the VRC suggests a promising model for SARS-CoV-2 infection and vaccine and therapy development. In keeping with a One Health approach, characterizing actual susceptibility and potential for bi-directional zoonotic/anthroponotic transfer in savanna monkey populations may be an important consideration for controlling COVID-19 epidemics in communities with frequent human/non-human primate interactions that, in many cases, may have limited health infrastructure.

Background Malaria is the leading cause of global paediatric mortality in children below 5 years of age. The number of fatalities has reduced significantly due to an expansion of control interventions but the development of new technologies remains necessary in order to achieve elimination. Recent attention has been focused on the release of genetically modified (GM) mosquitoes into natural vector populations as a mechanism of interrupting parasite transmission but despite successful in vivo laboratory studies, a detailed population genetic assessment, which must first precede any proposed field trial, has yet to be undertaken systematically. Here, the genetic structure of Anopheles gambiae populations in north-western Lake Victoria is explored to assess their suitability as candidates for a pilot field study release of GM mosquitoes. Methods 478 Anopheles gambiae mosquitoes were collected from six locations and a subset (N = 96) was selected for restriction site-associated DNA sequencing (RADseq). The resulting single nucleotide polymorphism (SNP) marker set was analysed for effective size (Ne), connectivity and population structure (PCA, FST). Results 5175 high-quality genome-wide SNPs were identified. A principal components analysis (PCA) of the collinear genomic regions illustrated that individuals clustered in concordance with geographic origin with some overlap between sites. Genetic differentiation between populations was varied with inter-island comparisons having the highest values (median FST 0.0480–0.0846). Ne estimates were generally small (124.2–1920.3). Conclusions A reduced-representation SNP marker set for genome-wide An. gambiae genetic analysis in the north-western Lake Victoria basin is reported. Island populations demonstrated low to moderate genetic differentiation and greater structure suggesting some limitation to migration. Smaller estimates of Ne indicate that an introduced effector transgene will be more susceptible to genetic drift but to ensure that it is driven to fixation a robust gene drive mechanism will likely be needed. These findings, together with their favourable location and suitability for frequent monitoring, indicate that the Ssese Islands contain several candidate field locations, which merit further evaluation as potential GM mosquito pilot release sites. Electronic supplementary material The online version of this article (10.1186/s12936-018-2432-0) contains supplementary material, which is available to authorized users.

Obtaining high-quality samples from wild animals is a major obstacle for genomic studies of many taxa, particularly at the population level, as collection methods for such samples are typically invasive. DNA from feces is easy to obtain noninvasively, but is dominated by bacterial and other non-host DNA. The high proportion of non-host DNA drastically reduces the efficiency of high-throughput sequencing for host animal genomics. To address this issue, we developed an inexpensive capture method for enriching host DNA from noninvasive fecal samples. Our method exploits natural differences in CpG-methylation density between vertebrate and bacterial genomes to preferentially bind and isolate host DNA from majority-bacterial samples. We demonstrate that the enrichment is robust, efficient, and compatible with downstream library preparation methods useful for population studies (e.g., RADseq). Compared to other enrichment strategies, our method is quick and inexpensive, adding only a negligible cost to sample preparation. In combination with downstream methods such as RADseq, our approach allows for cost-effective and customizable genomic-scale genotyping that was previously feasible in practice only with invasive samples. Because feces are widely available and convenient to collect, our method empowers researchers to explore genomic-scale population-level questions in organisms for which invasive sampling is challenging or undesirable.

The sustainability of malaria control in Africa is threatened by the rise of insecticide resistance in Anopheles mosquitoes, which transmit the disease(1). To gain a deeper understanding of how mosquito populations are evolving, here we sequenced the genomes of 765 specimens of Anopheles gambiae and Anopheles coluzzii sampled from 15 locations across Africa, and identified over 50 million single nucleotide polymorphisms within the accessible genome. These data revealed complex population structure and patterns of gene flow, with evidence of ancient expansions, recent bottlenecks, and local variation in effective population size. Strong signals of recent selection were observed in insecticide-resistance genes, with several sweeps spreading over large geographical distances and between species. The design of new tools for mosquito control using gene-drive systems will need to take account of high levels of genetic diversity in natural mosquito populations.

The shift from a hunter-gatherer (HG) to an agricultural (AG) mode of subsistence is believed to have been associated with profound changes in the burden and diversity of pathogens across human populations. Yet, the extent to which the advent of agriculture may have impacted the evolution of the human immune system remains unknown. Here we present a comparative study of variation in the transcriptional responses of peripheral blood mononuclear cells to bacterial and viral stimuli between Batwa rainforest hunter-gatherers and Bakiga agriculturalists from Uganda. We observed increased divergence between hunter-gatherers and agriculturalists in the transcriptional response to viruses compared to that for bacterial stimuli. We demonstrate that a significant fraction of these transcriptional differences are under genetic control, and we show that positive natural selection has helped to shape population differences in immune regulation. Across the set of genetic variants underlying inter-population immune response differences, however, the signatures of positive selection were disproportionately observed in the rainforest hunter-gatherers. This result is counter to expectations based on the popularized notion that shifts in pathogen exposure due to the advent of agriculture imposed radically heightened selective pressures in agriculturalist populations.

Kinda and gray-footed chacma baboons occupy opposite extremes of the body size distribution in extant baboons (genusPapio). In order to detect signatures of natural selection in these two species, we genotyped 24,790 genome-wide autosomal SNPs from populations of Zambian baboons using double digest RADseq. We scanned the genome for evidence of selection by identifying regions with extreme differentiation between populations. We find evidence of selection on body size influencing multiple genes in one or both species, includingFGF1, ATXN2, andPRKCE. We also find an enriched signal of selection associated with biological processes involved in multicellular organism growth and development, cell proliferation and cell growth, nutrient metabolism, and chondrocyte differentiation. Finally, we find that selection has impacted components of the CCKR signaling pathway, which regulates food intake and metabolism, and the JAK/STAT signaling pathway, which mediates the effect of cytokine signals on processes including epiphyseal chondrocyte proliferation essential for longitudinal bone growth. Our findings highlight promising avenues for future studies disentangling the genetic architecture of body size in primates including humans.

Hybridization in nature offers unique insights into the process of natural selection in incipient species and their hybrids. In order to evaluate the patterns and targets of selection, we examine a recently discovered baboon hybrid zone in the Kafue River valley of Zambia, where Kinda baboons (Papio kindae) and gray-footed chacma baboons (P. ursinus griseipes) coexist with hybridization. We genotyped baboons at 14,962 variable genome-wide autosomal markers using double-digest RADseq. We compare ancestry patterns from this genome-wide dataset to previously reported ancestry from mitochondrial-DNA and Y-chromosome sources. We also fit a Bayesian genomic cline model to scan for genes with extreme patterns of introgression. We show that the Kinda baboon Y chromosome has penetrated the species boundary to a greater extent than either mitochondrial DNA or the autosomal chromosomes. We also find evidence for overall restricted introgression in the JAK/STAT signaling pathway. Echoing results in other species including humans, we find evidence for enhanced and/or directional introgression of immune-related genes or pathways including the toll-like receptor pathway, the blood coagulation pathway, and the LY96 gene. Finally we show enhanced introgression and excess chacma baboon ancestry in the sperm tail gene ODF2. Together, our results elucidate the dynamics of introgressive hybridization in a primate system while highlighting genes and pathways under selection.

Recent studies suggest that closely related species can accumulate substantial genetic and phenotypic differences despite ongoing gene flow, thus challenging traditional ideas regarding the genetics of speciation. Baboons (genus Papio) are Old World monkeys consisting of six readily distinguishable species. Baboon species hybridize in the wild, and prior data imply a complex history of differentiation and introgression. We produced a reference genome assembly for the olive baboon (Papio anubis) and whole-genome sequence data for all six extant species. We document multiple episodes of admixture and introgression during the radiation of Papio baboons, thus demonstrating their value as a model of complex evolutionary divergence, hybridization, and reticulation. These results help inform our understanding of similar cases, including modern humans, Neanderthals, Denisovans, and other ancient hominins., The sequencing and analysis activities at the Human Genome Sequencing Center, Baylor College of Medicine, were supported by NIH (NHGRI) grants U54-HG003273 and U54-HG006484 to R.A.G. and GAC grant 1 S10 RR026605 to J. G. Reid. This research was also supported by NIH grant R01-GM59290 to M.A.B.; grants from the Austrian Science Fund (FWF-P24551 and FWF-W1225) and Vienna Science and Technology Fund (WWTF-MA16-061) to C.K.; grants from the Wellcome Trust (WT108749/Z/15/Z) and EMBL to B.A., F.J.M., and M.M.; grants VEGA 1/0719/14 and APVV-14-0253 to T. Vinar (Consortium Member); MINECO/FEDER grant, NIH U01-MH106874 grant, Howard Hughes International Early Career award, and Obra Social “La Caixa” award to T.M.-B.; NSF grants BNS83-03506 to J.P.-C.; NSF1029302 to J.P.-C., J.R., and C.J.J.; BNS96-15150 to J.P.-C., C.J.J., and T.D.; and National Geographic Society and Leakey Foundation grants to J.P.-C. and C.J.J. E.E.E. is an investigator of the Howard Hughes Medical Institute. This work was supported, in part, by U.S. NIH grant HG002385 to E.E.E.

The 1M SNP genotype data are available at the European Genome–Phenome archive, www.ebi.ac.uk/ega/ (accession numbers EGAS00001000605 and EGAS00001000908).; International audience; The shift from a hunter-gatherer to an agricultural mode of subsistence is believed to have been associated with profound changes in the burden and diversity of pathogens across human populations. Yet, the extent to which the advent of agriculture affected the evolution of the human immune system remains unknown. Here we present a comparative study of variation in the transcriptional responses of peripheral blood mononuclear cells to bacterial and viral stimuli between Batwa rainforest hunter-gatherers and Bakiga agriculturalists from Uganda. We observed increased divergence between hunter-gatherers and agriculturalists in the early transcriptional response to viruses compared with that for bacterial stimuli. We demonstrate that a significant fraction of these transcriptional differences are under genetic control and we show that positive natural selection has helped to shape population differences in immune regulation. Across the set of genetic variants underlying inter-population immune-response differences, however, the signatures of positive selection were disproportionately observed in the rainforest hunter-gatherers. This result is counter to expectations on the basis of the popularized notion that shifts in pathogen exposure due to the advent of agriculture imposed radically heightened selective pressures in agriculturalist populations.

Different human populations facing similar environmental challenges have sometimes evolved convergent biological adaptations, for example, hypoxia resistance at high altitudes and depigmented skin in northern latitudes on separate continents. The “pygmy” phenotype (small adult body size), characteristic of hunter-gatherer populations inhabiting both African and Asian tropical rainforests, is often highlighted as another case of convergent adaptation in humans. However, the degree to which phenotypic convergence in this polygenic trait is due to convergent versus population-specific genetic changes is unknown. To address this question, we analyzed high-coverage sequence data from the protein-coding portion of the genomes of two pairs of populations: Batwa rainforest hunter-gatherers and neighboring Bakiga agriculturalists from Uganda and Andamanese rainforest hunter-gatherers and Brahmin agriculturalists from India. We observed signatures of convergent positive selection between the rainforest hunter-gatherers across the set of genes with “growth factor binding” functions ( P < 0.001 ). Unexpectedly, for the rainforest groups, we also observed convergent and population-specific signatures of positive selection in pathways related to cardiac development (e.g., “cardiac muscle tissue development”; P = 0.001 ). We hypothesize that the growth hormone subresponsiveness likely underlying the adult small body-size phenotype may have led to compensatory changes in cardiac pathways, in which this hormone also plays an essential role. Importantly, in the agriculturalist populations, we did not observe similar patterns of positive selection on sets of genes associated with growth or cardiac development, indicating our results most likely reflect a history of convergent adaptation to the similar ecology of rainforests rather than a more general evolutionary pattern.

Different human populations facing similar environmental challenges have sometimes evolved convergent biological adaptations, for example hypoxia resistance at high altitudes and depigmented skin in northern latitudes on separate continents. The pygmy phenotype (small adult body size), a characteristic of hunter-gatherer populations inhabiting both African and Asian tropical rainforests, is often highlighted as another case of convergent adaptation in humans. However, the degree to which phenotypic convergence in this polygenic trait is due to convergent vs. population-specific genetic changes is unknown. To address this question, we analyzed high-coverage sequence data from the protein-coding portion of the genomes (exomes) of two pairs of populations, Batwa rainforest hunter-gatherers and neighboring Bakiga agriculturalists from Uganda, and Andamanese rainforest hunter-gatherers (Jarawa and Onge) and Brahmin agriculturalists from India. We observed signatures of convergent positive selection between the Batwa and Andamanese rainforest hunter-gatherers across the set of genes with annotated ‘growth factor binding’ functions (p< 0.001). Unexpectedly, for the rainforest groups we also observed convergent and population-specific signatures of positive selection in pathways related to cardiac development (e.g. ‘cardiac muscle tissue development’;p= 0.001). We hypothesize that the growth hormone sub-responsiveness likely underlying the pygmy phenotype may have led to compensatory changes in cardiac pathways, in which this hormone also plays an essential role. Importantly, in the agriculturalist populations we did not observe similar patterns of positive selection on sets of genes associated with either growth or cardiac development, indicating that our results most likely reflect a history of convergent adaptation to the similar ecology of rainforest hunter-gatherers rather than a more common or general evolutionary pattern for human populations.

Background Understanding population genetic structure in the malaria vector Anopheles gambiae (s.s.) is crucial to inform genetic control and manage insecticide resistance. Unfortunately, species characteristics such as high nucleotide diversity, large effective population size, recent range expansion, and high dispersal ability complicate the inference of genetic structure across its range in sub-Saharan Africa. The ocean, along with the Great Rift Valley, is one of the few recognized barriers to gene flow in this species, but the effect of inland lakes, which could be useful sites for initial testing of genetic control strategies, is relatively understudied. Here we examine Lake Victoria as a barrier between the Ugandan mainland and the Ssese Islands, which lie up to 60 km offshore. We use mitochondrial DNA (mtDNA) from populations sampled in 2002, 2012 and 2015, and perform Bayesian cluster analysis on mtDNA combined with microsatellite data previously generated from the same 2002 mosquito DNA samples. Results Hierarchical analysis of molecular variance and Bayesian clustering support significant differentiation between the mainland and lacustrine islands. In an mtDNA haplotype network constructed from this and previous data, haplotypes are shared even between localities separated by the Rift Valley, a result that more likely reflects retention of shared ancestral polymorphism than contemporary gene flow. Conclusions The relative genetic isolation of An. gambiae on the Ssese Islands, their small size, level terrain and ease of access from the mainland, the relative simplicity of the vectorial system, and the prevalence of malaria, are all attributes that recommend these islands as possible sites for the testing of genetic control strategies.

Natural history collections have long been used by morphologists, anatomists, and taxonomists to probe the evolutionary process and describe biological diversity. These biological archives also offer great opportunities for genetic research in taxonomy, conservation, systematics, and population biology. They allow assays of past populations, including those of extinct species, giving context to present patterns of genetic variation and direct measures of evolutionary processes. Despite this potential, museum specimens are difficult to work with because natural postmortem processes and preservation methods fragment and damage DNA. These problems have restricted geneticists’ ability to use natural history collections primarily by limiting how much of the genome can be surveyed. Recent advances in DNA sequencing technology, however, have radically changed this, making truly genomic studies from museum specimens possible. We review the opportunities and drawbacks of the use of museum specimens, and suggest how to best execute projects when incorporating such samples. Several high-throughput (HT) sequencing methodologies, including whole genome shotgun sequencing, sequence capture, and restriction digests (demonstrated here), can be used with archived biomaterials.

The sustainability of malaria control in Africa is threatened by rising levels of insecticide resistance, and new tools to prevent malaria transmission are urgently needed. To gain a better understanding of the mosquito populations that transmit malaria, we sequenced the genomes of 765 wild specimens ofAnopheles gambiaeandAnopheles coluzziisampled from 15 locations across Africa. The data reveal high levels of genetic diversity, with over 50 million single nucleotide polymorphisms across the 230 Mbp genome. We observe complex patterns of population structure and marked variations in local population size, some of which may be due at least in part to malaria control interventions. Insecticide resistance genes show strong signatures of recent selection associated with multiple independent mutations spreading over large geographical distances and between species. The genetic variability of natural populations substantially reduces the target space for novel gene-drive strategies for mosquito control. This large dataset provides a foundation for tracking the emergence and spread of insecticide resistance and developing new vector control tools.

Phylogeneticcomparativemethodsplayacriticalroleinourunderstandingof the adaptive origin of primate behaviors. To incorporate evolutionary history directly into comparative behavioral research, behavioral ecologists rely on strong, well- resolved phylogenetic trees. Phylogenies provide the framework on which behaviors can be compared and homologies can be distinguished from similarities due to conver- gent or parallel evolution. Phylogenetic reconstructions are also of critical importance when inferring the ancestral state of behavioral patterns and when suggesting the evolutionary changes that behavior has undergone. Improvements in genome sequenc- ing technologies have increased the amount of data available to researchers. Recently, several primate phylogenetic studies have used multiple loci to produce robust phylo- genetic trees that include hundreds of primate species. These trees are now commonly usedincomparativeanalysesandthereisaperceptionthatwehaveacompletepictureof the primate tree. But how confident can we be in those phylogenies? And how reliable are comparative analyses based on such trees? Herein, we argue that even recent molecular phylogenies should be treated cautiously because they rely on many assump- tions and have many shortcomings. Most phylogenetic studies do not model gene tree diversity and can produce misleading results, such as strong support for an incorrect speciestree,especiallyinthecaseofrapidandrecentradiations.Wediscussimplications thatincorrectphylogeniescanhaveforreconstructingtheevolutionofprimatebehaviors

Over the last two decades primatologists have benefited from the use of numerous molecular markers to study various aspects of primate behavior and evolutionary history. However, most of the studies to date have been based on a single locus, usually mitochondrial DNA, or a few nuclear markers, e.g., microsatellites. Unfortunately, the use of such markers not only is unable to address successfully important questions in primate population genetics and phylogenetics (mainly because of the discordance between gene tree and species tree), but also their development is often a time-consuming and expensive task. The advent of next-generation sequencing allows researchers to generate large amounts of genomic data for nonmodel organisms. However, whole genome sequencing is still cost prohibitive for most primate species. We here introduce a second-generation sequencing technique for genotyping thousands of genome-wide markers for nonmodel organisms. Restriction site–associated DNA sequencing (RAD-seq) reduces the complexity of the genome and allows inexpensive and fast discovery of thousands of markers in many individuals. Here, we describe the principles of this technique and we demonstrate its application in five primates, Microcebus sp., Cebus sp., Theropithecus gelada, Pan troglodytes, and Homo sapiens, representing some of the major lineages within the order. Despite technical and bioinformatic challenges, RAD-seq is a promising method for multilocus phylogenetic and population genetic studies in primates, particularly in young clades in which a high number of orthologous regions are likely to be found across populations or species.

In the endeavor to associate genetic variation with complex traits, closely related taxa are particularly fruitful for understanding the neurophysiological and genetic underpinnings of species-specific attributes. Similarity to humans has motivated research into nonhuman primate models, yet few studies of wild primates have investigated immediate causal factors of evolutionarily diverged social behaviors. Neurotransmitter differences have been invoked to explain the distinct behavioral suites of two baboon species in Awash, Ethiopia, which differ markedly in social behavior despite evolutionary propinquity. With this natural experiment, we test the hypothesis that genomic regions associated with monoamine neurotransmitters would be highly differentiated, and we identify a dopamine pathway as an outlier, highlighting the system as a potential cause of species-specific social behaviors. Dopamine levels and resultant variation in impulsivity were likely under differential selection in the species due to social system structure differences, with either brash or circumspect social behavior advantageous to secure mating opportunities depending on the social backdrop. Such comparative studies into the causes of the behavioral agendas that create and interact with social systems are of particular interest, and differences in temperament related to boldness and associated with dopamine variation likely played important roles in the evolution of all social, behaviorally complex animals, including baboons and humans.

Descriptions of primate diets are generally based on either direct observation of foraging behavior, morphological classification of food remains from feces, or analysis of the stomach contents of deceased individuals. Some diet items (e.g. insect prey), however, are difficult to identify visually, and observation conditions often do not permit adequate quantitative sampling of feeding behavior. Moreover, the taxonomically informative morphology of some food species (e.g. swallowed seeds, insect exoskeletons) may be destroyed by the digestive process. Because of these limitations, we used a metagenomic approach to conduct a preliminary, “proof of concept” study of interspecific variation in the insect component of the diets of six sympatric NewWorld monkeys known, based on observational field studies, to differ markedly in their feeding ecology. We used generalized arthropod polymerase chain reaction (PCR) primers and cloning to sequence mitochondrial DNA (mtDNA) sequences of the arthropod cytochrome b (CYT B) gene from fecal samples of wild woolly, titi, saki, capuchin, squirrel, and spider monkeys collected from a single sampling site in western Amazonia where these genera occur sympatrically. We then assigned preliminary taxonomic identifications to the sequences by basic local alignment search tool (BLAST) comparison to arthropod CYT Bsequences present in GenBank. This study is the first to use molecular techniques to identify insect prey in primate diets. The results suggest that a metagenomic approach may prove valuable in augmenting and corroborating observational data and increasing the resolution of primate diet studies, although the lack of comparative reference sequences for many South American insects limits the approach at present. As such reference data become available for more animal and plant taxa, this approach also holds promise for studying additional components of primate diets. Am. J. Primatol. 74:622–631, 2012. © 2012 Wiley Periodicals, Inc.

The ranges of small kinda (Papio kindae) and much larger grayfooted chacma (P. ursinus griseipes) baboons adjoin in the Kafue National Park, Zambia. In a visual survey of baboons at 48 sites in the Kafue River drainage we found that, contrary to previous reports, groups at the species interface near the town of Ngoma are phenotypically diverse and presumably formed by multigenerational hybridization. Mitochondrial and/or Y-chromosome genetic markers from fecal samples (N=164) collected at 29 sites support this conclusion. Groups with phenotypic signs of a history of hybridization also had taxon-specific mitochondria and Y-haplotypes from both parental species. Although the distribution of mitochondrial haplotypes largely mirrored that of external phenotypes, a significant proportion of male specimens from grayfoot as well as hybrid groups carried kinda Y-chromosomes, and kinda Y-chromosomes were involved in all observed cases of mitochondrial/Y-chromosome discordance. These observations are consistent with, though they do not prove, a population history in which the range of chacmas and the hybrid zone have advanced at the expense of the kinda range. They also suggest that, unexpectedly, kinda male×chacma female matings are much more common than the reciprocal cross in the ancestry of hybrids. We suggest that distinctive male kinda behavior and the "juvenile" appearance of kinda baboons of both sexes, perhaps combined with obstetric difficulties of a small kinda female carrying the large offspring of a chacma male, may account for this bias.

Obtaining high-quality samples from wild animals is a major obstacle for genomic studies of many taxa, particular at the population level, as collection methods for such samples are typically invasive. DNA from feces is easy to obtain noninvasively, but is dominated by a preponderance of bacterial and other non-host DNA. Because next-generation sequencing technology sequences DNA largely indiscriminately, the high proportion of exogenous DNA drastically reduces the efficiency of high-throughput sequencing for host animal genomics. In order to address this issue, we developed an inexpensive methylation-based capture method for enriching host DNA from noninvasively obtained fecal DNA samples. Our method exploits natural differences in CpG-methylation density between vertebrate and bacterial genomes to preferentially bind and isolate host DNA from majority-bacterial fecal DNA samples. We demonstrate that the enrichment is robust, efficient, and compatible with downstream library preparation methods useful for population studies (e.g., RADseq). Compared to other enrichment strategies, our method is quick and inexpensive, adding only a negligible cost to sample preparation for research that is often severely constrained by budgetary limitations. In combination with downstream methods such as RADseq, our approach allows for cost-effective and customizable genomic-scale genotyping that was previously feasible in practice only with invasive samples. Because feces are widely available and convenient to collect, our method empowers researchers to explore genomic-scale population-level questions in organisms for which invasive sampling is challenging or undesirable.

SummaryAnalysis of the Neandertal genome indicates gene flow between Neandertals and modern humans of Eurasia but not Africa. This surprising result is difficult to reconcile with current models of human origins and might have to do with insufficient African sampling.

Summary: HippDB catalogs every protein–protein interaction whose structure is available in the Protein Data Bank and which exhibits one or more helices at the interface. The Web site accepts queries on variables such as helix length and sequence, and it provides computational alanine scanning and change in solvent-accessible surface area values for every interfacial residue. HippDB is intended to serve as a starting point for structure-based small molecule and peptidomimetic drug development. Availability and implementation: HippDB is freely available on the web at http://www.nyu.edu/projects/arora/hippdb. The Web site is implemented in PHP, MySQL and Apache. Source code freely available for download at http://code.google.com/p/helidb, implemented in Perl and supported on Linux. Contact: ude.uyn@arora

Descriptions of primate diets are generally based on either direct observation of foraging behavior, morphological classification of food remains from feces, or analysis of the stomach contents of deceased individuals. Some diet items (e.g. insect prey), however, are difficult to identify visually, and observation conditions often do not permit adequate quantitative sampling of feeding behavior. Moreover, the taxonomically informative morphology of some food species (e.g. swallowed seeds, insect exoskeletons) may be destroyed by the digestive process. Because of these limitations, we used a metagenomic approach to conduct a preliminary, "proof of concept" study of interspecific variation in the insect component of the diets of six sympatric New World monkeys known, based on observational field studies, to differ markedly in their feeding ecology. We used generalized arthropod polymerase chain reaction (PCR) primers and cloning to sequence mitochondrial DNA (mtDNA) sequences of the arthropod cytochrome b (CYT B) gene from fecal samples of wild woolly, titi, saki, capuchin, squirrel, and spider monkeys collected from a single sampling site in western Amazonia where these genera occur sympatrically. We then assigned preliminary taxonomic identifications to the sequences by basic local alignment search tool (BLAST) comparison to arthropod CYT B sequences present in GenBank. This study is the first to use molecular techniques to identify insect prey in primate diets. The results suggest that a metagenomic approach may prove valuable in augmenting and corroborating observational data and increasing the resolution of primate diet studies, although the lack of comparative reference sequences for many South American insects limits the approach at present. As such reference data become available for more animal and plant taxa, this approach also holds promise for studying additional components of primate diets.

Summary: AluHunter is a database of taxon-specific primate Alu elements for use in phylogeny and population genetics. The software automatically isolates potentially polymorphic Alu insertions in sequences submitted to GenBank by screening the elements against reference genomes. The resultant database of variable markers is a valuable resource for researchers interested in characterizing Alu elements in their primate taxon of interest. Availability and Implementation: The AluHunter database can be accessed at http://www.aluhunter.com. Contact: cmb433@nyu.edu