The Strength of Selection against Neanderthal Introgression

Hybridization between humans and Neanderthals has resulted in a low level of Neanderthal ancestry scattered across the genomes of many modern-day humans. After hybridization, on average, selection appears to have removed Neanderthal alleles from the human population. Quantifying the strength and causes of this selection against Neanderthal ancestry is key to understanding our relationship to Neanderthals and, more broadly, how populations remain distinct after secondary contact. Here, we develop a novel method for estimating the genome-wide average strength of selection and the density of selected sites using estimates of Neanderthal allele frequency along the genomes of modern-day humans. We confirm that East Asians had somewhat higher initial levels of Neanderthal ancestry than Europeans even after accounting for selection. We find that the bulk of purifying selection against Neanderthal ancestry is best understood as acting on many weakly deleterious alleles. We propose that the majority of these alleles were effectively neutral—and segregating at high frequency—in Neanderthals, but became selected against after entering human populations of much larger effective size. While individually of small effect, these alleles potentially imposed a heavy genetic load on the early-generation human–Neanderthal hybrids. This work suggests that differences in effective population size may play a far more important role in shaping levels of introgression than previously thought.     

Editor's choice: AnABlast: a new in silico strategy for the genome-wide search of novel genes and fossil regions

Genome annotation, assisted by computer programs, is one of the great advances in modern biology. Nevertheless, the in silico identification of small and complex coding sequences is still challenging. We observed that amino acid sequences inferred from coding—but rarely from non-coding—DNA sequences accumulated alignments in low-stringency BLAST searches, suggesting that this alignments accumulation could be used to highlight coding regions in sequenced DNA. To investigate this possibility, we developed a computer program (AnABlast) that generates profiles of accumulated alignments in query amino acid sequences using a low-stringency BLAST strategy. To validate this approach, all six-frame translations of DNA sequences between every two annotated exons of the fission yeast genome were analysed with AnABlast. AnABlast-generated profiles identified three new copies of known genes, and four new genes supported by experimental evidence. New pseudogenes, ancestral carboxyl- and amino-terminal subtractions, complex gene rearrangements, and ancient fragments of mitDNA and of bacterial origin, were also inferred. Thus, this novel in silico approach provides a powerful tool to uncover new genes, as well as fossil-coding sequences, thus providing insight into the evolutionary history of annotated genomes.     

Mammalian development does not recapitulate suspected key transformations in the evolutionary detachment of the mammalian middle ear

The ectotympanic, malleus and incus of the developing mammalian middle ear (ME) are initially attached to the dentary via Meckel''s cartilage, betraying their origins from the primary jaw joint of land vertebrates. This recapitulation has prompted mostly unquantified suggestions that several suspected—but similarly unquantified—key evolutionary transformations leading to the mammalian ME are recapitulated in development, through negative allometry and posterior/medial displacement of ME bones relative to the jaw joint. Here we show, using µCT reconstructions, that neither allometric nor topological change is quantifiable in the pre-detachment ME development of six marsupials and two monotremes. Also, differential ME positioning in the two monotreme species is not recapitulated. This challenges the developmental prerequisites of widely cited evolutionary scenarios of definitive mammalian middle ear (DMME) evolution, highlighting the requirement for further fossil evidence to test these hypotheses. Possible association between rear molar eruption, full ME ossification and ME detachment in marsupials suggests functional divergence between dentary and ME as a trigger for developmental, and possibly also evolutionary, ME detachment. The stable positioning of the dentary and ME supports suggestions that a ‘partial mammalian middle ear’ as found in many mammaliaforms—probably with a cartilaginous Meckel''s cartilage—represents the only developmentally plausible evolutionary DMME precursor.     

The beak of the other finch: coevolution of genetic covariance structure and developmental modularity during adaptive evolution

The link between adaptation and evolutionary change remains the most central and least understood evolutionary problem. Rapid evolution and diversification of avian beaks is a textbook example of such a link, yet the mechanisms that enable beak''s precise adaptation and extensive adaptability are poorly understood. Often observed rapid evolutionary change in beaks is particularly puzzling in light of the neo-Darwinian model that necessitates coordinated changes in developmentally distinct precursors and correspondence between functional and genetic modularity, which should preclude evolutionary diversification. I show that during first 19 generations after colonization of a novel environment, house finches (Carpodacus mexicanus) express an array of distinct, but adaptively equivalent beak morphologies—a result of compensatory developmental interactions between beak length and width in accommodating microevolutionary change in beak depth. Directional selection was largely confined to the elimination of extremes formed by these developmental interactions, while long-term stabilizing selection along a single axis—beak depth—was mirrored in the structure of beak''s additive genetic covariance. These results emphasize three principal points. First, additive genetic covariance structure may represent a historical record of the most recurrent developmental and functional interactions. Second, adaptive equivalence of beak configurations shields genetic and developmental variation in individual components from depletion by natural selection. Third, compensatory developmental interactions among beak components can generate rapid reorganization of beak morphology under novel conditions and thus greatly facilitate both the evolution of precise adaptation and extensive diversification, thereby linking adaptation and adaptability in this classic example of Darwinian evolution.     

A bottom-up characterization of transfer functions for synthetic biology designs: lessons from enzymology

Within the field of synthetic biology, a rational design of genetic parts should include a causal understanding of their input-output responses—the so-called transfer function—and how to tune them. However, a commonly adopted strategy is to fit data to Hill-shaped curves without considering the underlying molecular mechanisms. Here we provide a novel mathematical formalization that allows prediction of the global behavior of a synthetic device by considering the actual information from the involved biological parts. This is achieved by adopting an enzymology-like framework, where transfer functions are described in terms of their input affinity constant and maximal response. As a proof of concept, we characterize a set of Lux homoserine-lactone-inducible genetic devices with different levels of Lux receptor and signal molecule. Our model fits the experimental results and predicts the impact of the receptor''s ribosome-binding site strength, as a tunable parameter that affects gene expression. The evolutionary implications are outlined.     

Segregational Drift Constrains the Evolutionary Rate of Prokaryotic Plasmids

Plasmids are extrachromosomal genetic elements in prokaryotes that have been recognized as important drivers of microbial ecology and evolution. Plasmids are found in multiple copies inside their host cell where independent emergence of mutations may lead to intracellular genetic heterogeneity. The intracellular plasmid diversity is thus subject to changes upon cell division. However, the effect of plasmid segregation on plasmid evolution remains understudied. Here, we show that genetic drift during cell division—segregational drift—leads to the rapid extinction of novel plasmid alleles. We established a novel experimental approach to control plasmid allele frequency at the levels of a single cell and the whole population. Following the dynamics of plasmid alleles in an evolution experiment, we find that the mode of plasmid inheritance—random or clustered—is an important determinant of plasmid allele dynamics. Phylogenetic reconstruction of our model plasmid in clinical isolates furthermore reveals a slow evolutionary rate of plasmid-encoded genes in comparison to chromosomal genes. Our study provides empirical evidence that genetic drift in plasmid evolution occurs at multiple levels: the host cell and the population of hosts. Segregational drift has implications for the evolutionary rate heterogeneity of extrachromosomal genetic elements.     

Neanderthal infant and adult infracranial remains from Marillac (Charente,France)

At the site of Marillac, near the Ligonne River in Marillac‐le‐Franc (Charente, France), a remarkable stratigraphic sequence has yielded a wealth of archaeological information, palaeoenvironmental data, as well as faunal and human remains. Marillac must have been a sinkhole used by Neanderthal groups as a hunting camp during MIS 4 (TL date 57,600 ± 4,600BP), where Quina Mousterian lithics and fragmented bones of reindeer predominate. This article describes three infracranial skeleton fragments. Two of them are from adults and consist of the incomplete shafts of a right radius (Marillac 24) and a left fibula (Marillac 26). The third fragment is the diaphysis of the right femur of an immature individual (Marillac 25), the size and shape of which resembles those from Teshik‐Tash and could be assigned to a child of a similar age. The three fossils have been compared with the remains of other Neanderthals or anatomically Modern Humans (AMH). Furthermore, the comparison of the infantile femora, Marillac 25 and Teshik‐Tash, with the remains of several European children from the early Middle Ages clearly demonstrates the robustness and rounded shape of both Neanderthal diaphyses. Evidence of peri‐mortem manipulations have been identified on all three bones, with spiral fractures, percussion pits and, in the case of the radius and femur, unquestionable cutmarks made with flint implements, probably during defleshing. Traces of periostosis appear on the fibula fragment and on the immature femoral diaphysis, although their aetiology remains unknown. Am J Phys Anthropol 155:99–113, 2014. © 2014 Wiley Periodicals, Inc.     

Use of artificial genomes in assessing methods for atypical gene detection

Parametric methods for identifying laterally transferred genes exploit the directional mutational biases unique to each genome. Yet the development of new, more robust methods—as well as the evaluation and proper implementation of existing methods—relies on an arbitrary assessment of performance using real genomes, where the evolutionary histories of genes are not known. We have used the framework of a generalized hidden Markov model to create artificial genomes modeled after genuine genomes. To model a genome, “core” genes—those displaying patterns of mutational biases shared among large numbers of genes—are identified by a novel gene clustering approach based on the Akaike information criterion. Gene models derived from multiple “core” gene clusters are used to generate an artificial genome that models the properties of a genuine genome. Chimeric artificial genomes—representing those having experienced lateral gene transfer—were created by combining genes from multiple artificial genomes, and the performance of the parametric methods for identifying “atypical” genes was assessed directly. We found that a hidden Markov model that included multiple gene models, each trained on sets of genes representing the range of genotypic variability within a genome, could produce artificial genomes that mimicked the properties of genuine genomes. Moreover, different methods for detecting foreign genes performed differently—i.e., they had different sets of strengths and weaknesses—when identifying atypical genes within chimeric artificial genomes.     

Sexual antagonism for resistance and tolerance to infection in Drosophila melanogaster

A critical task in evolutionary genetics is to explain the persistence of heritable variation in fitness-related traits such as immunity. Ecological factors can maintain genetic variation in immunity, but less is known about the role of other factors, such as antagonistic pleiotropy, on immunity. Sexually dimorphic immunity—with females often being more immune-competent—may maintain variation in immunity in dioecious populations. Most eco-immunological studies assess host resistance to parasites rather than the host''s ability to maintain fitness during infection (tolerance). Distinguishing between resistance and tolerance is important as they are thought to have markedly different evolutionary and epidemiological outcomes. Few studies have investigated tolerance in animals, and the extent of sexual dimorphism in tolerance is unknown. Using males and females from 50 Drosophila melanogaster genotypes, we investigated possible sources of genetic variation for immunity by assessing both resistance and tolerance to the common bacterial pathogen Pseudomonas aeruginosa. We found evidence of sexual dimorphism and sexual antagonism for resistance and tolerance, and a trade-off between the two traits. Our findings suggest that antagonistic pleiotropy may be a major contributor to variation in immunity, with implications for host–parasite coevolution.     

An Analysis of the Alleged Skeletal Remains of Carin G?ring

In 1991, treasure hunters found skeletal remains in an area close to the destroyed country residence of former Nazi leader Hermann Göring in northeastern Berlin. The remains, which were believed to belong to Carin Göring, who was buried at the site, were examined to determine whether it was possible to make a positive identification. The anthropological analysis showed that the remains come from an adult woman. The DNA analysis of several bone elements showed female sex, and a reference sample from Carin''s son revealed mtDNA sequences identical to the remains. The profile has one nucleotide difference from the Cambridge reference sequence (rCRS), the common variant 263G. A database search resulted in a frequency of this mtDNA sequence of about 10% out of more than 7,000 European haplotypes. The mtDNA sequence found in the ulna, the cranium and the reference sample is, thus, very common among Europeans. Therefore, nuclear DNA analysis was attempted. The remains as well as a sample from Carin''s son were successfully analysed for the three nuclear markers TH01, D7S820 and D8S1179. The nuclear DNA analysis of the two samples revealed one shared allele for each of the three markers, supporting a mother and son relationship. This genetic information together with anthropological and historical files provides an additional piece of circumstantial evidence in our efforts to identify the remains of Carin Göring.     

Humble origins for a successful strategy: complete enrolment in early Cambrian olenellid trilobites

Trilobites are typified by the behavioural and morphological ability to enrol their bodies, most probably as a defence mechanism against adverse environmental conditions or predators. Although most trilobites could enrol at least partially, there is uncertainty about whether olenellids—among the most phylogenetically and stratigraphically basal representatives—could perform this behaviour because of their poorly caudalized trunk and scarcity of coaptative devices. Here, we report complete—but not encapsulating—enrolment for the olenellid genus Mummaspis from the early Cambrian Mural Formation in Alberta, the earliest direct evidence of this strategy in the fossil record of polymerid trilobites. Complete enrolment in olenellids was achieved through a combination of ancestral morphological features, and thus provides new information on the character polarity associated with this key trilobite adaptation.     

The future of evolutionary diversity in reef corals

One-third of the world''s reef-building corals are facing heightened extinction risk from climate change and other anthropogenic impacts. Previous studies have shown that such threats are not distributed randomly across the coral tree of life, and future extinctions have the potential to disproportionately reduce the phylogenetic diversity of this group on a global scale. However, the impact of such losses on a regional scale remains poorly known. In this study, we use phylogenetic metrics in conjunction with geographical distributions of living reef coral species to model how extinctions are likely to affect evolutionary diversity across different ecoregions. Based on two measures—phylogenetic diversity and phylogenetic species variability—we highlight regions with the largest losses of evolutionary diversity and hence of potential conservation interest. Notably, the projected loss of evolutionary diversity is relatively low in the most species-rich areas such as the Coral Triangle, while many regions with fewer species stand to lose much larger shares of their diversity. We also suggest that for complex ecosystems like coral reefs it is important to consider changes in phylogenetic species variability; areas with disproportionate declines in this measure should be of concern even if phylogenetic diversity is not as impacted. These findings underscore the importance of integrating evolutionary history into conservation planning for safeguarding the future diversity of coral reefs.     

Genetic pest management and the background genetics of release strains

Genetic pest management (GPM) methods involve releasing modified versions of a pest species to mate with wild pests in the target area. Proposed for a wide range of applications in public health, agriculture and conservation, most progress has been made with pest insects. Offspring of the released modified insects and wild pests carry the modification—which might be transgenes, artificially introduced Wolbachia or genetic damage from radiation, for example—but they also carry a complete haploid genome from their laboratory-reared parent, as well as one from their wild parent. Unless these F₁ hybrids are completely unable to reproduce, further mating will lead to introgression of DNA sequences from the release strain into the wild population. We discuss issues around strain selection and the potential consequences of such introgression. We conclude that such introgression is probably harmless in almost all circumstances, and could, in theory, provide specific additional benefits to the release programme. We outline population monitoring approaches that could be used, going forward, to determine how background genetics may affect GPM.This article is part of the theme issue ‘Novel control strategies for mosquito-borne diseases’.     

The Synthesis Paradigm in Genetics

Experimental genetics with model organisms and mathematically explicit genetic theory are generally considered to be the major paradigms by which progress in genetics is achieved. Here I argue that this view is incomplete and that pivotal advances in genetics—and other fields of biology—are also made by synthesizing disparate threads of extant information rather than generating new information from experiments or formal theory. Because of the explosive expansion of information in numerous “-omics” data banks, and the fragmentation of genetics into numerous subdisciplines, the importance of the synthesis paradigm will likely expand with time.MAJOR advances in the field of genetics have been developed on a foundation supported by three major pillars (i.e., paradigms, by which I mean a framework of basic assumptions, logical approaches, and methodologies), two of which are widely known and appreciated while the third is rarely even acknowledged. The first major paradigm is experimental genetics, especially in the context of model organisms. The work of Thomas H. Morgan and his colleagues at Cal Tech during the early 20th century is a classic example of this approach. A succession of elegant experimental studies by this research team led to the development of the Drosophila melanogaster model system, which Morgan et al. (1915) used to construct the first genomic map that included genes assigned to precise locations on all of an organism’s chromosomes. Their accumulated experimental results also contributed importantly to their book, The Mechanism of Mendelian Heredity (1915), which many consider to be the catalyst that launched the modern era of genetics.The second paradigm is mathematically explicit genetic theory. The succession of genetical theory papers published throughout the first half of the 20th century by Ronald A. Fisher is a classic example of this approach. Fisher’s work reconciled a fundamental rift in the early history of modern genetics—i.e., the genetic approaches of the Mendelians (advocated by William Bateson and Hugo de Vries) vs. the Galtonians (also known as the biometricians, represented in particular by Karl Pearson and Walter F. Weldon)—by showing that Mendelian particulate inheritance could be unified with the quantitative genetics used to analyze continuously varying traits such as height and weight (Fisher 1918). Although Darwin developed the basic framework of evolution, it was Fisher—and contemporary theoreticians Sewall Wright and J. B. S. Haldane—who integrated this qualitative idea into a quantitatively explicit genetic theory that led to the modern synthesis of evolution and launched the field of evolutionary genetics (also known as population genetics and summarized in Fisher’s now classic book, The Genetical Theory of Natural Selection, first published in 1930). Of course, some theory in genetics is not mathematically explicit, such as the “chromosomal theory of inheritance” or the “central dogma.” But this form of theory usually represents the culmination of studies using the experimental genetics paradigm rather than a unique approach to genetics.Most major advances in genetics have been achieved via one, the other, or a combination of these experimental and theoretical paradigms. But there is a well-known exception: Watson and Crick’s discovery of the structure of DNA (Watson and Crick 1953a,b, ∼11,000 combined citations—throughout, numbers of citations are taken from Google Scholar—, and arguably the pivotal publications that launched the modern field of molecular genetics). Watson and Crick used no mathematical genetic theory, nor did they do any critical experiments; instead, they integrated many threads of established information (some unpublished) to deduce the chemical structure of the hereditary material, i.e., the DNA double helix and how this structure could explain gene replication. Although later experiments, such as those of Meselson and Stahl (1958) on DNA replication, would ultimately confirm the deduced structure and replication of DNA that was proposed by Watson and Crick, the pivotal publications of these researchers used neither the experimental nor the theory paradigms of genetics. Their approach exemplifies what I will call the “synthesis paradigm.” Watson and Crick’s work demonstrates that there is actually a trichotomy of approaches—the experimental, theoretical-mathematical, and theoretical-synthetic approaches—that combine like interwoven, reinforcing strands in a cord of historical advances in genetics.In the next few sections I describe other instances in which the synthesis paradigm has been of critical importance in the field of genetics. This set of examples is meant to be illustrative and by no means exhaustive. Next I illustrate how the synthesis paradigm has been of critical importance in other fields of biology. Finally, I describe how a fuller appreciation of the synthesis paradigm can influence the training of the next cohort of geneticists and the career trajectory of current geneticists.     

Evaluation of Two Approaches for Assessing the Genetic Similarity of Virioplankton Populations as Defined by Genome Size

Viral production estimates show that virioplankton communities turn over rapidly in aquatic ecosystems. Thus, it is likely that the genetic identity of viral populations comprising the virioplankton also change over temporal and spatial scales, reflecting shifts in viral-host interactions. However, there are few approaches that can provide data on the genotypic identity of viral populations at low cost and with the sample throughput necessary to assess dynamic changes in the virioplankton. This study examined two of these approaches—T4-like major capsid protein (g23) gene polymorphism and randomly amplified polymorphic DNA-PCR (RAPD-PCR) fingerprinting—to ask how well each technique could track differences in virioplankton populations over time and geographic location. Seasonal changes in overall virioplankton composition were apparent from pulsed-field gel electrophoresis (PFGE) analysis. T4-like phages containing similar g23 proteins were found within both small- and large-genome populations, including populations from different geographic locations and times. The surprising occurrence of T4-like g23 within small genomic groups (23 to 64 kb) indicated that the genome size range of T4-like phages may be broader than previously believed. In contrast, RAPD-PCR fingerprinting detected high genotypic similarity within PFGE bands from the same location, time, and genome size class without the requirement for DNA sequencing. Unlike g23 polymorphism, RAPD-PCR fingerprints showed a greater temporal than geographic variation. Thus, while polymorphism in a viral signature gene, such as g23, can be a powerful tool for inferring evolutionary relationships, the degree to which this approach can capture fine-scale variability within virioplankton populations is less clear.     

Hypervariable Bkm DNA Loci in a Moth, Ephestia kuehniella : Does Transposition Cause Restriction Fragment Length Polymorphism?

Bkm sequences, originally isolated from snake satellite DNA, are a component of eukaryote genomes with a preferential location on sex chromosomes. In the Ephestia genome, owing to the presence of only a few Bkm-positive BamHI restriction fragments and to extensive restriction fragment length polymorphisms between and within inbred strains, a genetic crossbreeding analysis was feasible. No sex linkage of Bkm was detected. Instead—depending on the strain—two or three autosomal Bkm DNA loci were identified. All three loci were located on different chromosomes. Fragment length and transmission of fragments was stable in some crosses. In others, changes in fragment length or loss of the Bkm component were observed, probably depending on the source strain of the fragment. The anomalous genetic behaviour is best accounted for by the assumption that Bkm sequences are included in mobile genetic elements.     

Adaptive basis of geographic variation: genetic,phenotypic and environmental differences among beach mouse populations

A major goal in evolutionary biology is to understand how and why populations differentiate, both genetically and phenotypically, as they invade a novel habitat. A classical example of adaptation is the pale colour of beach mice, relative to their dark mainland ancestors, which colonized the isolated sandy dunes and barrier islands on Florida''s Gulf Coast. However, much less is known about differentiation among the Gulf Coast beach mice, which comprise five subspecies linearly arrayed on Florida''s shoreline. Here, we test the role of selection in maintaining variation among these beach mouse subspecies at multiple levels—phenotype, genotype and the environments they inhabit. While all beach subspecies have light pelage, they differ significantly in colour pattern. These subspecies are also genetically distinct: pair-wise F_st-values range from 0.23 to 0.63 and levels of gene flow are low. However, we did not find a correlation between phenotypic and genetic distance. Instead, we find a significant association between the average ‘lightness’ of each subspecies and the brightness of the substrate it inhabits: the two most genetically divergent subspecies occupy the most similar habitats and have converged on phenotype, whereas the most genetically similar subspecies occupy the most different environments and have divergent phenotypes. Moreover, allelic variation at the pigmentation gene, Mc1r, is statistically correlated with these colour differences but not with variation at other genetic loci. Together, these results suggest that natural selection for camouflage—via changes in Mc1r allele frequency—contributes to pigment differentiation among beach mouse subspecies.     

Cooperation,clumping and the evolution of multicellularity

The evolution of multicellular organisms represents one of the major evolutionary transitions in the history of life. A potential advantage of forming multicellular clumps is that it provides an efficiency benefit to pre-existing cooperation, such as the production of extracellular ‘public goods’. However, this is complicated by the fact that cooperation could jointly evolve with clumping, and clumping could have multiple consequences for the evolution of cooperation. We model the evolution of clumping and a cooperative public good, showing that (i) when considered separately, both clumping and public goods production gradually increase with increasing genetic relatedness; (ii) in contrast, when the traits evolve jointly, a small increase in relatedness can lead to a major shift in evolutionary outcome—from a non-clumping state with low public goods production to a cooperative clumping state with high values of both traits; (iii) high relatedness makes it easier to get to the cooperative clumping state and (iv) clumping can be inhibited when it increases the number of cells that the benefits of cooperation must be shared with, but promoted when it increases relatedness between those cells. Overall, our results suggest that public goods sharing can facilitate the formation of well-integrated cooperative clumps as a first step in the evolution of multicellularity.     

Joint Estimation of Contamination,Error and Demography for Nuclear DNA from Ancient Humans

When sequencing an ancient DNA sample from a hominin fossil, DNA from present-day humans involved in excavation and extraction will be sequenced along with the endogenous material. This type of contamination is problematic for downstream analyses as it will introduce a bias towards the population of the contaminating individual(s). Quantifying the extent of contamination is a crucial step as it allows researchers to account for possible biases that may arise in downstream genetic analyses. Here, we present an MCMC algorithm to co-estimate the contamination rate, sequencing error rate and demographic parameters—including drift times and admixture rates—for an ancient nuclear genome obtained from human remains, when the putative contaminating DNA comes from present-day humans. We assume we have a large panel representing the putative contaminant population (e.g. European, East Asian or African). The method is implemented in a C++ program called ‘Demographic Inference with Contamination and Error’ (DICE). We applied it to simulations and genome data from ancient Neanderthals and modern humans. With reasonable levels of genome sequence coverage (>3X), we find we can recover accurate estimates of all these parameters, even when the contamination rate is as high as 50%.