Detecting Selection in Multiple Populations by Modeling Ancestral Admixture Components.

One of the most powerful and commonly used approaches for detecting local adaptation in the genome is the identification of extreme allele frequency differences between populations. In this article, we present a new maximum likelihood method for finding regions under positive selection. It is based on a Gaussian approximation to allele frequency changes and it incorporates admixture between populations. The method can analyze multiple populations simultaneously and retains power to detect selection signatures specific to ancestry components that are not representative of any extant populations. Using simulated data, we compare our method to related approaches, and show that it is orders of magnitude faster than the state-of-the-art, while retaining similar or higher power for most simulation scenarios. We also apply it to human genomic data and identify loci with extreme genetic differentiation between major geographic groups. Many of the genes identified are previously known selected loci relating to hair pigmentation and morphology, skin, and eye pigmentation. We also identify new candidate regions, including various selected loci in the Native American component of admixed Mexican-Americans. These involve diverse biological functions, such as immunity, fat distribution, food intake, vision, and hair development.

Population genomics of the Viking world.

The maritime expansion of Scandinavian populations during the Viking Age (about AD 750-1050) was a far-flung transformation in world history1,2. Here we sequenced the genomes of 442 humans from archaeological sites across Europe and Greenland (to a median depth of about 1×) to understand the global influence of this expansion. We find the Viking period involved gene flow into Scandinavia from the south and east. We observe genetic structure within Scandinavia, with diversity hotspots in the south and restricted gene flow within Scandinavia. We find evidence for a major influx of Danish ancestry into England; a Swedish influx into the Baltic; and Norwegian influx into Ireland, Iceland and Greenland. Additionally, we see substantial ancestry from elsewhere in Europe entering Scandinavia during the Viking Age. Our ancient DNA analysis also revealed that a Viking expedition included close family members. By comparing with modern populations, we find that pigmentation-associated loci have undergone strong population differentiation during the past millennium, and trace positively selected loci-including the lactase-persistence allele of LCT and alleles of ANKA that are associated with the immune response-in detail. We conclude that the Viking diaspora was characterized by substantial transregional engagement: distinct populations influenced the genomic makeup of different regions of Europe, and Scandinavia experienced increased contact with the rest of the continent.

Ancestral Population Genomics with Jocx, a Coalescent Hidden Markov Model.

Coalescence theory lets us probe the past demographics of present-day genetic samples and much information about the past can be gleaned from variation in rates of coalescence event as we trace genetic lineages back in time. Fewer and fewer lineages will remain, however, so there is a limit to how far back we can explore. Without recombination, we would not be able to explore ancient speciation events because of this-any meaningful species concept would require that individuals of one species are closer related than they are to individuals of another species, once speciation is complete. Recombination, however, opens a window to the deeper past. By scanning along a genomic alignment, we get a sequential variant of the coalescence process as it looked at the time of the speciation. This pattern of coalescence times is fixed at speciation time and does not erode with time; although accumulated mutations and genomic rearrangements will eventually hide the signal, it enables us to glance at events in the past that would not be observable without recombination. So-called coalescence hidden Markov models allow us to exploit this, and in this chapter, we present the tool Jocx that uses a framework of these models to infer demographic parameters in ancient speciation events.

The comparative genomics and complex population history of Papio baboons.

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.

Parallel adaptation of rabbit populations to myxoma virus.

In the 1950s the myxoma virus was released into European rabbit populations in Australia and Europe, decimating populations and resulting in the rapid evolution of resistance. We investigated the genetic basis of resistance by comparing the exomes of rabbits collected before and after the pandemic. We found a strong pattern of parallel evolution, with selection on standing genetic variation favoring the same alleles in Australia, France, and the United Kingdom. Many of these changes occurred in immunity-related genes, supporting a polygenic basis of resistance. We experimentally validated the role of several genes in viral replication and showed that selection acting on an interferon protein has increased the protein's antiviral effect.

Physiological and Genetic Adaptations to Diving in Sea Nomads.

Understanding the physiology and genetics of human hypoxia tolerance has important medical implications, but this phenomenon has thus far only been investigated in high-altitude human populations. Another system, yet to be explored, is humans who engage in breath-hold diving. The indigenous Bajau people ("Sea Nomads") of Southeast Asia live a subsistence lifestyle based on breath-hold diving and are renowned for their extraordinary breath-holding abilities. However, it is unknown whether this has a genetic basis. Using a comparative genomic study, we show that natural selection on genetic variants in the PDE10A gene have increased spleen size in the Bajau, providing them with a larger reservoir of oxygenated red blood cells. We also find evidence of strong selection specific to the Bajau on BDKRB2, a gene affecting the human diving reflex. Thus, the Bajau, and possibly other diving populations, provide a new opportunity to study human adaptation to hypoxia tolerance. VIDEO ABSTRACT.

Natural Selection on Genes Related to Cardiovascular Health in High-Altitude Adapted Andeans.

The increase in red blood cell mass (polycythemia) due to the reduced oxygen availability (hypoxia) of residence at high altitude or other conditions is generally thought to be beneficial in terms of increasing tissue oxygen supply. However, the extreme polycythemia and accompanying increased mortality due to heart failure in chronic mountain sickness most likely reduces fitness. Tibetan highlanders have adapted to high altitude, possibly in part via the selection of genetic variants associated with reduced polycythemic response to hypoxia. In contrast, high-altitude-adapted Quechua- and Aymara-speaking inhabitants of the Andean Altiplano are not protected from high-altitude polycythemia in the same way, yet they exhibit other adaptive features for which the genetic underpinnings remain obscure. Here, we used whole-genome sequencing to scan high-altitude Andeans for signals of selection. The genes showing the strongest evidence of selection-including BRINP3, NOS2, and TBX5-are associated with cardiovascular development and function but are not in the response-to-hypoxia pathway. Using association mapping, we demonstrated that the haplotypes under selection are associated with phenotypic variations related to cardiovascular health. We hypothesize that selection in response to hypoxia in Andeans could have vascular effects and could serve to mitigate the deleterious effects of polycythemia rather than reduce polycythemia itself.

Fast admixture analysis and population tree estimation for SNP and NGS data.

MOTIVATION: Structure methods are highly used population genetic methods for classifying individuals in a sample fractionally into discrete ancestry components. CONTRIBUTION: We introduce a new optimization algorithm for the classical STRUCTURE model in a maximum likelihood framework. Using analyses of real data we show that the new method finds solutions with higher likelihoods than the state-of-the-art method in the same computational time. The optimization algorithm is also applicable to models based on genotype likelihoods, that can account for the uncertainty in genotype-calling associated with Next Generation Sequencing (NGS) data. We also present a new method for estimating population trees from ancestry components using a Gaussian approximation. Using coalescence simulations of diverging populations, we explore the adequacy of the STRUCTURE-style models and the Gaussian assumption for identifying ancestry components correctly and for inferring the correct tree. In most cases, ancestry components are inferred correctly, although sample sizes and times since admixture can influence the results. We show that the popular Gaussian approximation tends to perform poorly under extreme divergence scenarios e.g. with very long branch lengths, but the topologies of the population trees are accurately inferred in all scenarios explored. The new methods are implemented together with appropriate visualization tools in the software package Ohana. AVAILABILITY AND IMPLEMENTATION: Ohana is publicly available at https://github.com/jade-cheng/ohana . In addition to source code and installation instructions, we also provide example work-flows in the project wiki site. CONTACT: jade.cheng@birc.au.dk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx.

BACKGROUND: Genomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction. RESULTS: We generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species. CONCLUSIONS: The genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.

A genomic history of Aboriginal Australia.

The population history of Aboriginal Australians remains largely uncharacterized. Here we generate high-coverage genomes for 83 Aboriginal Australians (speakers of Pama-Nyungan languages) and 25 Papuans from the New Guinea Highlands. We find that Papuan and Aboriginal Australian ancestors diversified 25-40 thousand years ago (kya), suggesting pre-Holocene population structure in the ancient continent of Sahul (Australia, New Guinea and Tasmania). However, all of the studied Aboriginal Australians descend from a single founding population that differentiated ~10-32 kya. We infer a population expansion in northeast Australia during the Holocene epoch (past 10,000 years) associated with limited gene flow from this region to the rest of Australia, consistent with the spread of the Pama-Nyungan languages. We estimate that Aboriginal Australians and Papuans diverged from Eurasians 51-72 kya, following a single out-of-Africa dispersal, and subsequently admixed with archaic populations. Finally, we report evidence of selection in Aboriginal Australians potentially associated with living in the desert.

Nationwide Genomic Study in Denmark Reveals Remarkable Population Homogeneity.

Denmark has played a substantial role in the history of Northern Europe. Through a nationwide scientific outreach initiative, we collected genetic and anthropometrical data from ∼800 high school students and used them to elucidate the genetic makeup of the Danish population, as well as to assess polygenic predictions of phenotypic traits in adolescents. We observed remarkable homogeneity across different geographic regions, although we could still detect weak signals of genetic structure reflecting the history of the country. Denmark presented genomic affinity with primarily neighboring countries with overall resemblance of decreasing weight from Britain, Sweden, Norway, Germany, and France. A Polish admixture signal was detected in Zealand and Funen, and our date estimates coincided with historical evidence of Wend settlements in the south of Denmark. We also observed considerably diverse demographic histories among Scandinavian countries, with Denmark having the smallest current effective population size compared to Norway and Sweden. Finally, we found that polygenic prediction of self-reported adolescent height in the population was remarkably accurate (R2 = 0.639 ± 0.015). The high homogeneity of the Danish population could render population structure a lesser concern for the upcoming large-scale gene-mapping studies in the country.

Spitting for Science: Danish High School Students Commit to a Large-Scale Self-Reported Genetic Study.

Scientific outreach delivers science to the people. But it can also deliver people to the science. In this work, we report our experience from a large-scale public engagement project promoting genomic literacy among Danish high school students with the additional benefit of collecting data for studying the genetic makeup of the Danish population. Not only did we confirm that students have a great interest in their genetic past, but we were also gratified to see that, with the right motivation, adolescents can provide high-quality data for genetic studies.

Caveolae: One Function or Many?

Caveolae are small, bulb-shaped plasma membrane invaginations. Mutations that ablate caveolae lead to diverse phenotypes in mice and humans, making it challenging to uncover their molecular mechanisms. Caveolae have been described to function in endocytosis and transcytosis (a specialized form of endocytosis) and in maintaining membrane lipid composition, as well as acting as signaling platforms. New data also support a model in which the central function of caveolae could be related to the protection of cells from mechanical stress within the plasma membrane. We present evidence for these diverse roles and consider in vitro and in vivo experiments confirming a mechanoprotective role. We conclude by highlighting current gaps in our knowledge of how mechanical signals may be transduced by caveolae.

Caveolae protect endothelial cells from membrane rupture during increased cardiac output.

Caveolae are strikingly abundant in endothelial cells, yet the physiological functions of caveolae in endothelium and other tissues remain incompletely understood. Previous studies suggest a mechanoprotective role, but whether this is relevant under the mechanical forces experienced by endothelial cells in vivo is unclear. In this study we have sought to determine whether endothelial caveolae disassemble under increased hemodynamic forces, and whether caveolae help prevent acute rupture of the plasma membrane under these conditions. Experiments in cultured cells established biochemical assays for disassembly of caveolar protein complexes, and assays for acute loss of plasma membrane integrity. In vivo, we demonstrate that caveolae in endothelial cells of the lung and cardiac muscle disassemble in response to acute increases in cardiac output. Electron microscopy and two-photon imaging reveal that the plasma membrane of microvascular endothelial cells in caveolin 1(-/-) mice is much more susceptible to acute rupture when cardiac output is increased. These data imply that mechanoprotection through disassembly of caveolae is important for endothelial function in vivo.

Ancestral population genomics using coalescence hidden Markov models and heuristic optimisation algorithms.

With full genome data from several closely related species now readily available, we have the ultimate data for demographic inference. Exploiting these full genomes, however, requires models that can explicitly model recombination along alignments of full chromosomal length. Over the last decade a class of models, based on the sequential Markov coalescence model combined with hidden Markov models, has been developed and used to make inference in simple demographic scenarios. To move forward to more complex demographic modelling we need better and more automated ways of specifying these models and efficient optimisation algorithms for inferring the parameters in complex and often high-dimensional models. In this paper we present a framework for building such coalescence hidden Markov models for pairwise alignments and present results for using heuristic optimisation algorithms for parameter estimation. We show that we can build more complex demographic models than our previous frameworks and that we obtain more accurate parameter estimates using heuristic optimisation algorithms than when using our previous gradient based approaches. Our new framework provides a flexible way of constructing coalescence hidden Markov models almost automatically. While estimating parameters in more complex models is still challenging we show that using heuristic optimisation algorithms we still get a fairly good accuracy.

Opposing structural changes in two symmetrical polypeptides bring about opposing changes to the thermal stability of a complex integral membrane protein.

The relationship between membrane protein structure and thermal stability has been examined in the reaction centre from the bacterium Rhodobacter sphaeroides, a complex membrane protein comprising three polypeptide chains and 10 cofactors. The core of this protein exhibits an approximate twofold symmetry, the cofactors being held in two membrane-spanning branches by two polypeptides, termed L and M, that have very similar folds. In assays of the thermal stability of wild-type and mutant reaction centres embedded in the native bilayer membrane, replacement of a Phe at position 197 of the M polypeptide by His produced an increase in stability, whereas an opposing replacement of His by Phe at the symmetrical position 168 of the L-polypeptide produced a decrease in stability. In light of the known X-ray crystal structures of wild-type and mutant variants of this protein, and further mutagenesis, it is concluded that these stability changes result from the introduction or removal, respectively, of a hydrogen bond between the side-chain of the His and that of an Asn located two positions along the M or L polypeptide chain, in addition to a hydrogen bond between the His side-chain and an adjacent bacteriochlorophyll cofactor.

Organelle dynamics and membrane trafficking in apoptosis and autophagy.

The accurate control of cell death is a vital aspect of development in metazoans and plays crucial roles in the prevention of disease. Apoptosis is the main form of regulated cell death in multicellular organisms, although there are other contributory pathways. During apoptosis, mammalian cells undergo dramatic changes in organelle structure ad organisation that define the apoptotic execution phase. Although the roles of apoptotic protease machinery (the caspases) in these rearrangements are quite well understood, the purpose of organelle disruption during cell death is not yet entirely appreciated. Indeed, recent evidence implicates caspase targeting of organellar proteins and subsequent organelle disruption upstream of apoptotic execution proper, suggesting the existence of pathways linking organelle damage to cell death. In this review, we describe the changes to the endomembrane system that are inherent during the apoptotic execution phase, and examine the evidence for endomembrane-mediated pathways towards apoptotic execution. We also discuss aspects of the molecular control of autophagy - an important contributor to a cell's response to stress, and a membrane trafficking process whose regulation is linked to the apoptotic machinery at multiple levels.

Structural responses to cavity-creating mutations in an integral membrane protein.

X-ray crystallography has been used to investigate the extent of structural changes in mutants of the purple bacterial reaction center that assemble without a particular ubiquinone or bacteriopheophytin cofactor. In the case of the bacteriopheophytin-exclusion mutant, in which Ala M149 was replaced by Trp (AM149W), the quality of protein crystals was improved over that seen in previous work by minimizing illumination, time, and temperature during the purification protocol and carrying out crystal growth at 4 degrees C after overnight incubation at 18 degrees C. The X-ray crystal structure of the AM149W mutant, determined to a resolution of 2.2 A, showed very little change in protein structure despite the absence of the bacteriopheophytin cofactor. Changes in the electron density map in the region of the cofactor binding site could be accounted for by changes in the conformation of the phytol side chains of adjacent cofactors and the presence of a buried water molecule. Residues lining the vacated binding pocket did not show any significant changes in conformation or increases in disorder as assessed through crystallographic atomic displacement parameters (B-factors). The X-ray crystal structure of a reaction center lacking the primary acceptor ubiquinone through mutation of Ala M248 to Trp (AM248W) was also determined, to a resolution of 2.8 A. Again, despite the absence of an internal cofactor only very minor changes in protein structure were observed. This is in contrast to a previous report on a reaction center lacking this ubiquinone through mutation of Ala M260 to Trp (AM260W) where more extensive changes in structure were apparent. All three mutant reaction centers showed a decrease in thermal stability when housed in the native membrane, but this decrease was smaller for the AM260W mutant than the AM248W complex, possibly due to beneficial effects of the observed changes in protein structure. The lack of major changes in protein structure despite the absence of large internal cofactors is discussed in terms of protein rigidity, the protective influence of the adaptable membrane environment, and the role of small molecules and ions as packing material in the internal cavities created by this type of mutation.