Impact of mutation rate and selection at linked sites on DNA variation across the genomes of humans and other homininae

DNA diversity varies across the genome of many species. Variation in diversity across a genome might arise from regional variation in the mutation rate, variation in the intensity and mode of natural selection, and regional variation in the recombination rate. We show that both non-coding and non-synonymous diversity are positively correlated to a measure of the mutation rate and the recombination rate and negatively correlated to the density of conserved sequences in 50KB windows across the genomes of humans and non-human homininae. Interestingly, we find that while non-coding diversity is equally affected by these three genomic variables, non-synonymous diversity is mostly dominated by the density of conserved sequences. The positive correlation between diversity and our measure of the mutation rate seems to be largely a direct consequence of regions with higher mutation rates having more diversity. However, the positive correlation with recombination rate and the negative correlation with the density of conserved sequences suggests that selection at linked sites also affect levels of diversity. This is supported by the observation that the ratio of the number of non-synonymous to non-coding polymorphisms is negatively correlated to a measure of the effective population size across the genome. We show these patterns persist even when we restrict our analysis to GC-conservative mutations, demonstrating that the patterns are not driven by GC biased gene conversion. In conclusion, our comparative analyses describe how recombination rate, gene density, and mutation rate interact to produce the patterns of DNA diversity that we observe along the hominine genomes

A hidden Markov model approach for determining expression from genomic tiling micro arrays

BACKGROUND: Genomic tiling micro arrays have great potential for identifying previously undiscovered coding as well as non-coding transcription. To-date, however, analyses of these data have been performed in an ad hoc fashion. RESULTS: We present a probabilistic procedure, ExpressHMM, that adaptively models tiling data prior to predicting expression on genomic sequence. A hidden Markov model (HMM) is used to model the distributions of tiling array probe scores in expressed and non-expressed regions. The HMM is trained on sets of probes mapped to regions of annotated expression and non-expression. Subsequently, prediction of transcribed fragments is made on tiled genomic sequence. The prediction is accompanied by an expression probability curve for visual inspection of the supporting evidence. We test ExpressHMM on data from the Cheng et al. (2005) tiling array experiments on ten Human chromosomes [1]. Results can be downloaded and viewed from our web site [2]. CONCLUSION: The value of adaptive modelling of fluorescence scores prior to categorisation into expressed and non-expressed probes is demonstrated. Our results indicate that our adaptive approach is superior to the previous analysis in terms of nucleotide sensitivity and transfrag specificity

Genomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca

Background Despite their regional economic importance and being increasingly reared globally, the origins and evolution of the llama and alpaca remain poorly understood. Here we report reference genomes for the llama, and for the guanaco and vicuña (their putative wild progenitors), compare these with the published alpaca genome, and resequence seven individuals of all four species to better understand domestication and introgression between the llama and alpaca. Results Phylogenomic analysis confirms that the llama was domesticated from the guanaco and the alpaca from the vicuña. Introgression was much higher in the alpaca genome (36%) than the llama (5%) and could be dated close to the time of the Spanish conquest, approximately 500 years ago. Introgression patterns are at their most variable on the X-chromosome of the alpaca, featuring 53 genes known to have deleterious X-linked phenotypes in humans. Strong genome-wide introgression signatures include olfactory receptor complexes into both species, hypertension resistance into alpaca, and fleece/fiber traits into llama. Genomic signatures of domestication in the llama include male reproductive traits, while in alpaca feature fleece characteristics, olfaction-related and hypoxia adaptation traits. Expression analysis of the introgressed region that is syntenic to human HSA4q21, a gene cluster previously associated with hypertension in humans under hypoxic conditions, shows a previously undocumented role for PRDM8 downregulation as a potential transcriptional regulation mechanism, analogous to that previously reported at high altitude for hypoxia-inducible factor 1α. Conclusions The unprecedented introgression signatures within both domestic camelid genomes may reflect post-conquest changes in agriculture and the breakdown of traditional management practices

HIV-1 Nef Targets MHC-I and CD4 for Degradation Via a Final Common β-COP–Dependent Pathway in T Cells

To facilitate viral infection and spread, HIV-1 Nef disrupts the surface expression of the viral receptor (CD4) and molecules capable of presenting HIV antigens to the immune system (MHC-I). To accomplish this, Nef binds to the cytoplasmic tails of both molecules and then, by mechanisms that are not well understood, disrupts the trafficking of each molecule in different ways. Specifically, Nef promotes CD4 internalization after it has been transported to the cell surface, whereas Nef uses the clathrin adaptor, AP-1, to disrupt normal transport of MHC-I from the TGN to the cell surface. Despite these differences in initial intracellular trafficking, we demonstrate that MHC-I and CD4 are ultimately found in the same Rab7+ vesicles and are both targeted for degradation via the activity of the Nef-interacting protein, β-COP. Moreover, we demonstrate that Nef contains two separable β-COP binding sites. One site, an arginine (RXR) motif in the N-terminal α helical domain of Nef, is necessary for maximal MHC-I degradation. The second site, composed of a di-acidic motif located in the C-terminal loop domain of Nef, is needed for efficient CD4 degradation. The requirement for redundant motifs with distinct roles supports a model in which Nef exists in multiple conformational states that allow access to different motifs, depending upon which cellular target is bound by Nef

Species-specific responses of Late Quaternary megafauna to climate and humans

Despite decades of research, the roles of climate and humans in driving the dramatic extinctions of large-bodied mammals during the Late Quaternary remain contentious. We use ancient DNA, species distribution models and the human fossil record to elucidate how climate and humans shaped the demographic history of woolly rhinoceros, woolly mammoth, wild horse, reindeer, bison and musk ox. We show that climate has been a major driver of population change over the past 50,000 years. However, each species responds differently to the effects of climatic shifts, habitat redistribution and human encroachment. Although climate change alone can explain the extinction of some species, such as Eurasian musk ox and woolly rhinoceros, a combination of climatic and anthropogenic effects appears to be responsible for the extinction of others, including Eurasian steppe bison and wild horse. We find no genetic signature or any distinctive range dynamics distinguishing extinct from surviving species, underscoring the challenges associated with predicting future responses of extant mammals to climate and human-mediated habitat change.This paper is in the memory of our friend and colleague Dr. Andrei Sher, who was a major contributor of this study. Dr Sher died unexpectedly, but his major contributions to the field of Quaternary science will be remembered and appreciated for many years to come. We are grateful to Dr. Adrian Lister and Dr. Tony Stuart for guides and discussions. Thanks to Tina B. Brandt, Dr. Bryan Hockett and Alice Telka for laboratory help and samples and to L. Malik R. Thrane for his work on the megafauna locality database. Data taken from the Stage 3 project was partly funded by Grant #F/757/A from the Leverhulme Trust, together with a grant from the McDonald Grants and Awards Fund. We acknowledge the Danish National Research Foundation, the Lundbeck Foundation, the Danish Council for Independent Research and the US National Science Foundation for financial suppor

The time scale of recombination rate evolution in great apes

We present three linkage-disequilibrium (LD)-based recombination maps generated using whole-genome sequence data from 10 Nigerian chimpanzees, 13 bonobos, and 15 western gorillas, collected as part of the Great Ape Genome Project (Prado-Martinez J, et al. 2013. Great ape genetic diversity and population history. Nature 499:471-475). We also identified species-specific recombination hotspots in each group using a modified LDhot framework, which greatly improves statistical power to detect hotspots at varying strengths. We show that fewer hotspots are shared among chimpanzee subspecies than within human populations, further narrowing the time scale of complete hotspot turnover. Further, using species-specific PRDM9 sequences to predict potential binding sites (PBS), we show higher predicted PRDM9 binding in recombination hotspots as compared to matched cold spot regions in multiple great ape species, including at least one chimpanzee subspecies. We found that correlations between broad-scale recombination rates decline more rapidly than nucleotide divergence between species. We also compared the skew of recombination rates at centromeres and telomeres between species and show a skew from chromosome means extending as far as 10-15Mb from chromosome ends. Further, we examined broad-scale recombination rate changes near a translocation in gorillas and found minimal differences as compared to other great ape species perhaps because the coordinates relative to the chromosome ends were unaffected. Finally, on the basis of multiple linear regression analysis, we found that various correlates of recombination rate persist throughout the African great apes including repeats, diversity, and divergence. Our study is the first to analyze within- And between-species genome-wide recombination rate variation in several close relatives