National Institute for Mathematical and Biological Synthesis

Founder-event speciation, where a rare jump dispersal event founds a new genetically isolated lineage, has long been considered crucial by many historical biogeographers, but its importance is disputed within the vicariance school. Probabilistic modeling of geographic range evolution creates the potential to test different biogeographical models against data using standard statistical model choice procedures, as long as multiple models are available. I re-implement the Dispersal-Extinction-Cladogenesis (DEC) model of LAGRANGE in the R package BioGeoBEARS, and modify it to create a new model, DEC + J, which adds founder-event speciation, the importance of which is governed by a new free parameter, [Formula: see text]. The identifiability of DEC and DEC + J is tested on data sets simulated under a wide range of macroevolutionary models where geography evolves jointly with lineage birth/death events. The results confirm that DEC and DEC + J are identifiable even though these models ignore the fact that molecular phylogenies are missing many cladogenesis and extinction events. The simulations also indicate that DEC will have substantially increased errors in ancestral range estimation and parameter inference when the true model includes + J. DEC and DEC + J are compared on 13 empirical data sets drawn from studies of island clades. Likelihood-ratio tests indicate that all clades reject DEC, and AICc model weights show large to overwhelming support for DEC + J, for the first time verifying the importance of founder-event speciation in island clades via statistical model choice. Under DEC + J, ancestral nodes are usually estimated to have ranges occupying only one island, rather than the widespread ancestors often favored by DEC. These results indicate that the assumptions of historical biogeography models can have large impacts on inference and require testing and comparison with statistical methods.

Historical biogeography has been characterized by a large diversity of methods and unresolved debates about which processes, such as dispersal or vicariance, are most important for explaining distributions. A new R package, BioGeoBEARS, implements many models in a common likelihood framework, so that standard statistical model selection procedures can be applied to let the data choose the best model. Available models include a likelihood version of DIVA (“DIVALIKE”), LAGRANGE’s DEC model, and BAYAREA, as well as “+J” versions of these models which include founder-event speciation, an important process left out of most inference methods. I use BioGeoBEARS on a large sample of island and non-island clades (including two fossil clades) to show that founder-event speciation is a crucial process in almost every clade, and that most published datasets reject the non-J models currently in widespread use. BioGeoBEARS is open-source and freely available for installation at the Comprehensive R Archive Network at http://CRAN.R-project.org/package=BioGeoBEARS. A step-by-step tutorial is available at http://phylo.wikidot.com/biogeobears.

Uncovering the genetic and evolutionary basis of local adaptation is a major focus of evolutionary biology. The recent development of cost-effective methods for obtaining high-quality genome-scale data makes it possible to identify some of the loci responsible for adaptive differences among populations. Two basic approaches for identifying putatively locally adaptive loci have been developed and are broadly used: one that identifies loci with unusually high genetic differentiation among populations (differentiation outlier methods) and one that searches for correlations between local population allele frequencies and local environments (genetic-environment association methods). Here, we review the promises and challenges of these genome scan methods, including correcting for the confounding influence of a species' demographic history, biases caused by missing aspects of the genome, matching scales of environmental data with population structure, and other statistical considerations. In each case, we make suggestions for best practices for maximizing the accuracy and efficiency of genome scans to detect the underlying genetic basis of local adaptation. With attention to their current limitations, genome scan methods can be an important tool in finding the genetic basis of adaptive evolutionary change.

Biologists have long been fascinated by the exceptionally high diversity displayed by some evolutionary groups. Adaptive radiation in such clades is not only spectacular, but is also an extremely complex process influenced by a variety of ecological, genetic, and developmental factors and strongly dependent on historical contingencies. Using modeling approaches, we identify 10 general patterns concerning the temporal, spatial, and genetic/morphological properties of adaptive radiation. Some of these are strongly supported by empirical work, whereas for others, empirical support is more tentative. In almost all cases, more data are needed. Future progress in our understanding of adaptive radiation will be most successful if theoretical and empirical approaches are integrated, as has happened in other areas of evolutionary biology.

Ecological opportunity--through entry into a new environment, the origin of a key innovation or extinction of antagonists--is widely thought to link ecological population dynamics to evolutionary diversification. The population-level processes arising from ecological opportunity are well documented under the concept of ecological release. However, there is little consensus as to how these processes promote phenotypic diversification, rapid speciation and adaptive radiation. We propose that ecological opportunity could promote adaptive radiation by generating specific changes to the selective regimes acting on natural populations, both by relaxing effective stabilizing selection and by creating conditions that ultimately generate diversifying selection. We assess theoretical and empirical evidence for these effects of ecological opportunity and review emerging phylogenetic approaches that attempt to detect the signature of ecological opportunity across geological time. Finally, we evaluate the evidence for the evolutionary effects of ecological opportunity in the diversification of Caribbean Anolis lizards. Some of the processes that could link ecological opportunity to adaptive radiation are well documented, but others remain unsupported. We suggest that more study is required to characterize the form of natural selection acting on natural populations and to better describe the relationship between ecological opportunity and speciation rates.

Comparative methods used to study patterns of evolutionary change in a continuous trait on a phylogeny range from Brownian motion processes to models where the trait is assumed to evolve according to an Ornstein-Uhlenbeck (OU) process. Although these models have proved useful in a variety of contexts, they still do not cover all the scenarios biologists want to examine. For models based on the OU process, model complexity is restricted in current implementations by assuming that the rate of stochastic motion and the strength of selection do not vary among selective regimes. Here, we expand the OU model of adaptive evolution to include models that variously relax the assumption of a constant rate and strength of selection. In its most general form, the methods described here can assign each selective regime a separate trait optimum, a rate of stochastic motion parameter, and a parameter for the strength of selection. We use simulations to show that our models can detect meaningful differences in the evolutionary process, especially with larger sample sizes. We also illustrate our method using an empirical example of genome size evolution within a large flowering plant clade.

Human cooperation is highly unusual. We live in large groups composed mostly of non-relatives. Evolutionists have proposed a number of explanations for this pattern, including cultural group selection and extensions of more general processes such as reciprocity, kin selection, and multi-level selection acting on genes. Evolutionary processes are consilient; they affect several different empirical domains, such as patterns of behavior and the proximal drivers of that behavior. In this target article, we sketch the evidence from five domains that bear on the explanatory adequacy of cultural group selection and competing hypotheses to explain human cooperation. Does cultural transmission constitute an inheritance system that can evolve in a Darwinian fashion? Are the norms that underpin institutions among the cultural traits so transmitted? Do we observe sufficient variation at the level of groups of considerable size for group selection to be a plausible process? Do human groups compete, and do success and failure in competition depend upon cultural variation? Do we observe adaptations for cooperation in humans that most plausibly arose by cultural group selection? If the answer to one of these questions is "no," then we must look to other hypotheses. We present evidence, including quantitative evidence, that the answer to all of the questions is "yes" and argue that we must take the cultural group selection hypothesis seriously. If culturally transmitted systems of rules (institutions) that limit individual deviance organize cooperation in human societies, then it is not clear that any extant alternative to cultural group selection can be a complete explanation.

The distribution of diversity can vary considerably from clade to clade. Attempts to understand these patterns often employ state-dependent speciation and extinction models to determine whether the evolution of a particular novel trait has increased speciation rates and/or decreased extinction rates. It is still unclear, however, whether these models are uncovering important drivers of diversification, or whether they are simply pointing to more complex patterns involving many unmeasured and co-distributed factors. Here we describe an extension to the popular state-dependent speciation and extinction models that specifically accounts for the presence of unmeasured factors that could impact diversification rates estimated for the states of any observed trait, addressing at least one major criticism of BiSSE (Binary State Speciation and Extinction) methods. Specifically, our model, which we refer to as HiSSE (Hidden State Speciation and Extinction), assumes that related to each observed state in the model are "hidden" states that exhibit potentially distinct diversification dynamics and transition rates than the observed states in isolation. We also demonstrate how our model can be used as character-independent diversification models that allow for a complex diversification process that is independent of the evolution of a character. Under rigorous simulation tests and when applied to empirical data, we find that HiSSE performs reasonably well, and can at least detect net diversification rate differences between observed and hidden states and detect when diversification rate differences do not correlate with the observed states. We discuss the remaining issues with state-dependent speciation and extinction models in general, and the important ways in which HiSSE provides a more nuanced understanding of trait-dependent diversification.

Historical biogeography has been characterized by a large diversity of methods and unresolved debates about which processes, such as dispersal or vicariance, are most important for explaining distributions. A new R package, BioGeoBEARS, implements many models in a common likelihood framework, so that standard statistical model selection procedures can be applied to let the data choose the best model. Available models include a likelihood version of DIVA (“DIVALIKE”), LAGRANGE’s DEC model, and BAYAREA, as well as “+J” versions of these models which include founder-event speciation, an important process left out of most inference methods. I use BioGeoBEARS on a large sample of island and non-island clades (including two fossil clades) to show that founder-event speciation is a crucial process in almost every clade, and that most published datasets reject the non-J models currently in widespread use. BioGeoBEARS is open-source and freely available for installation at the Comprehensive R Archive Network at http://CRAN.R-project.org/package=BioGeoBEARS. A step-by-step tutorial is available at http://phylo.wikidot.com/biogeobears.

Abstract Understanding how and why populations evolve is of fundamental importance to molecular ecology. Restriction site‐associated DNA sequencing ( RAD seq), a popular reduced representation method, has ushered in a new era of genome‐scale research for assessing population structure, hybridization, demographic history, phylogeography and migration. RAD seq has also been widely used to conduct genome scans to detect loci involved in adaptive divergence among natural populations. Here, we examine the capacity of those RAD seq‐based genome scan studies to detect loci involved in local adaptation. To understand what proportion of the genome is missed by RAD seq studies, we developed a simple model using different numbers of RAD ‐tags, genome sizes and extents of linkage disequilibrium (length of haplotype blocks). Under the best‐case modelling scenario, we found that RAD seq using six‐ or eight‐base pair cutting restriction enzymes would fail to sample many regions of the genome, especially for species with short linkage disequilibrium. We then surveyed recent studies that have used RAD seq for genome scans and found that the median density of markers across these studies was 4.08 RAD ‐tag markers per megabase (one marker per 245 kb). The length of linkage disequilibrium for many species is one to three orders of magnitude less than density of the typical recent RAD seq study. Thus, we conclude that genome scans based on RAD seq data alone, while useful for studies of neutral genetic variation and genetic population structure, will likely miss many loci under selection in studies of local adaptation.

The science and management of infectious disease are entering a new stage. Increasingly public policy to manage epidemics focuses on motivating people, through social distancing policies, to alter their behavior to reduce contacts and reduce public disease risk. Person-to-person contacts drive human disease dynamics. People value such contacts and are willing to accept some disease risk to gain contact-related benefits. The cost-benefit trade-offs that shape contact behavior, and hence the course of epidemics, are often only implicitly incorporated in epidemiological models. This approach creates difficulty in parsing out the effects of adaptive behavior. We use an epidemiological-economic model of disease dynamics to explicitly model the trade-offs that drive person-to-person contact decisions. Results indicate that including adaptive human behavior significantly changes the predicted course of epidemics and that this inclusion has implications for parameter estimation and interpretation and for the development of social distancing policies. Acknowledging adaptive behavior requires a shift in thinking about epidemiological processes and parameters.

Summary Climate refugia, locations where taxa survive periods of regionally adverse climate, are thought to be critical for maintaining biodiversity through the glacial–interglacial climate changes of the Q uaternary. A critical research need is to better integrate and reconcile the three major lines of evidence used to infer the existence of past refugia – fossil records, species distribution models and phylogeographic surveys – in order to characterize the complex spatiotemporal trajectories of species and populations in and out of refugia. Here we review the complementary strengths, limitations and new advances for these three approaches. We provide case studies to illustrate their combined application, and point the way towards new opportunities for synthesizing these disparate lines of evidence. Case studies with E uropean beech, Q inghai spruce and D ouglas‐fir illustrate how the combination of these three approaches successfully resolves complex species histories not attainable from any one approach. Promising new statistical techniques can capitalize on the strengths of each method and provide a robust quantitative reconstruction of species history. Studying past refugia can help identify contemporary refugia and clarify their conservation significance, in particular by elucidating the fine‐scale processes and the particular geographic locations that buffer species against rapidly changing climate. Contents Summary 38 I. Climate refugia: biogeographical and conservation significance 38 II. Approaches for reconstructing refugia: strengths, limitations and recent advances 39 III. Climate refugia of the past: three case studies 46 IV. New integrative approaches to reconstructing refugia 47 V. How can historical refugia inform us about future refugia? 48 VI. Concluding thoughts 49 Acknowledgements 49 References 49

How did human societies evolve from small groups, integrated by face-to-face cooperation, to huge anonymous societies of today, typically organized as states? Why is there so much variation in the ability of different human populations to construct viable states? Existing theories are usually formulated as verbal models and, as a result, do not yield sharply defined, quantitative predictions that could be unambiguously tested with data. Here we develop a cultural evolutionary model that predicts where and when the largest-scale complex societies arose in human history. The central premise of the model, which we test, is that costly institutions that enabled large human groups to function without splitting up evolved as a result of intense competition between societies-primarily warfare. Warfare intensity, in turn, depended on the spread of historically attested military technologies (e.g., chariots and cavalry) and on geographic factors (e.g., rugged landscape). The model was simulated within a realistic landscape of the Afroeurasian landmass and its predictions were tested against a large dataset documenting the spatiotemporal distribution of historical large-scale societies in Afroeurasia between 1,500 BCE and 1,500 CE. The model-predicted pattern of spread of large-scale societies was very similar to the observed one. Overall, the model explained 65% of variance in the data. An alternative model, omitting the effect of diffusing military technologies, explained only 16% of variance. Our results support theories that emphasize the role of institutions in state-building and suggest a possible explanation why a long history of statehood is positively correlated with political stability, institutional quality, and income per capita.

Animal acoustic communication often takes the form of complex sequences, made up of multiple distinct acoustic units. Apart from the well-known example of birdsong, other animals such as insects, amphibians, and mammals (including bats, rodents, primates, and cetaceans) also generate complex acoustic sequences. Occasionally, such as with birdsong, the adaptive role of these sequences seems clear (e.g. mate attraction and territorial defence). More often however, researchers have only begun to characterise - let alone understand - the significance and meaning of acoustic sequences. Hypotheses abound, but there is little agreement as to how sequences should be defined and analysed. Our review aims to outline suitable methods for testing these hypotheses, and to describe the major limitations to our current and near-future knowledge on questions of acoustic sequences. This review and prospectus is the result of a collaborative effort between 43 scientists from the fields of animal behaviour, ecology and evolution, signal processing, machine learning, quantitative linguistics, and information theory, who gathered for a 2013 workshop entitled, 'Analysing vocal sequences in animals'. Our goal is to present not just a review of the state of the art, but to propose a methodological framework that summarises what we suggest are the best practices for research in this field, across taxa and across disciplines. We also provide a tutorial-style introduction to some of the most promising algorithmic approaches for analysing sequences. We divide our review into three sections: identifying the distinct units of an acoustic sequence, describing the different ways that information can be contained within a sequence, and analysing the structure of that sequence. Each of these sections is further subdivided to address the key questions and approaches in that area. We propose a uniform, systematic, and comprehensive approach to studying sequences, with the goal of clarifying research terms used in different fields, and facilitating collaboration and comparative studies. Allowing greater interdisciplinary collaboration will facilitate the investigation of many important questions in the evolution of communication and sociality.

Arbuscular mycorrhizal fungi (AMF) can increase plant fitness under certain environmental conditions. Among the mechanisms that may drive this mutualism, the most studied is provisioning of nutrients by AMF in exchange for carbon from plant hosts. However, AMF may also provide a suite of non-nutritional benefits to plants including improved water uptake, disease resistance, plant chemical defense, soil aggregation, and allelochemical transport and protection. Here, we use a meta-analysis of 93 studies to assess the relative effect of AMF on nutritional and non-nutritional factors that may influence plant fitness. We find that the positive effects of AMF on soil aggregation, water flow and disease resistance are equal to the effect of AMF on plant nitrogen and phosphorus uptake. However, AMF had no effect on the uptake of other nutrients, plant water content, allelopathic transport or production of chemical defense compounds. We suggest future research directions, including experimentally assessing the relative contribution on plant fitness of AMF interactions by untangling the independence of alternative benefits of AMF from an increase in nutrient uptake. This will lead to a more holistic view of the mycorrhizal-plant association and a more accurate picture of the net impact on the plant or plant community in question.

Identifying and monitoring locally adaptive genetic variation can have direct utility for conserving species at risk, especially when management may include actions such as translocations for restoration, genetic rescue, or assisted gene flow. However, genomic studies of local adaptation require careful planning to be successful, and in some cases may not be a worthwhile use of resources. Here, we offer an adaptive management framework to help conservation biologists and managers decide when genomics is likely to be effective in detecting local adaptation, and how to plan assessment and monitoring of adaptive variation to address conservation objectives. Studies of adaptive variation using genomic tools will inform conservation actions in many cases, including applications such as assisted gene flow and identifying conservation units. In others, assessing genetic diversity, inbreeding, and demographics using selectively neutral genetic markers may be most useful. And in some cases, local adaptation may be assessed more efficiently using alternative approaches such as common garden experiments. Here, we identify key considerations of genomics studies of locally adaptive variation, provide a road map for successful collaborations with genomics experts including key issues for study design and data analysis, and offer guidelines for interpreting and using results from genomic assessments to inform monitoring programs and conservation actions.

Human behavior is strongly affected by culturally transmitted norms and values. Certain norms are internalized (i.e., acting according to a norm becomes an end in itself rather than merely a tool in achieving certain goals or avoiding social sanctions). Humans' capacity to internalize norms likely evolved in our ancestors to simplify solving certain challenges-including social ones. Here we study theoretically the evolutionary origins of the capacity to internalize norms. In our models, individuals can choose to participate in collective actions as well as punish free riders. In making their decisions, individuals attempt to maximize a utility function in which normative values are initially irrelevant but play an increasingly important role if the ability to internalize norms emerges. Using agent-based simulations, we show that norm internalization evolves under a wide range of conditions so that cooperation becomes "instinctive." Norm internalization evolves much more easily and has much larger effects on behavior if groups promote peer punishment of free riders. Promoting only participation in collective actions is not effective. Typically, intermediate levels of norm internalization are most frequent but there are also cases with relatively small frequencies of "oversocialized" individuals willing to make extreme sacrifices for their groups no matter material costs, as well as "undersocialized" individuals completely immune to social norms. Evolving the ability to internalize norms was likely a crucial step on the path to large-scale human cooperation.

A fundamental assumption in invasion biology is that most invasive species exhibit enhanced performance in their introduced range relative to their home ranges. This idea has given rise to numerous hypotheses explaining "invasion success" by virtue of altered ecological and evolutionary pressures. There are surprisingly few data, however, testing the underlying assumption that the performance of introduced populations, including organism size, reproductive output, and abundance, is enhanced in their introduced compared to their native range. Here, we combined data from published studies to test this hypothesis for 26 plant and 27 animal species that are considered to be invasive. On average, individuals of these 53 species were indeed larger, more fecund, and more abundant in their introduced ranges. The overall mean, however, belied significant variability among species, as roughly half of the investigated species (N=27) performed similarly when compared to conspecific populations in their native range. Thus, although some invasive species are performing better in their new ranges, the pattern is not universal, and just as many are performing largely the same across ranges.

Gene tree topologies have proven a powerful data source for various tasks, including species tree inference and species delimitation. Consequently, methods for computing probabilities of gene trees within species trees have been developed and widely used in probabilistic inference frameworks. All these methods assume an underlying multispecies coalescent model. However, when reticulate evolutionary events such as hybridization occur, these methods are inadequate, as they do not account for such events. Methods that account for both hybridization and deep coalescence in computing the probability of a gene tree topology currently exist for very limited cases. However, no such methods exist for general cases, owing primarily to the fact that it is currently unknown how to compute the probability of a gene tree topology within the branches of a phylogenetic network. Here we present a novel method for computing the probability of gene tree topologies on phylogenetic networks and demonstrate its application to the inference of hybridization in the presence of incomplete lineage sorting. We reanalyze a Saccharomyces species data set for which multiple analyses had converged on a species tree candidate. Using our method, though, we show that an evolutionary hypothesis involving hybridization in this group has better support than one of strict divergence. A similar reanalysis on a group of three Drosophila species shows that the data is consistent with hybridization. Further, using extensive simulation studies, we demonstrate the power of gene tree topologies at obtaining accurate estimates of branch lengths and hybridization probabilities of a given phylogenetic network. Finally, we discuss identifiability issues with detecting hybridization, particularly in cases that involve extinction or incomplete sampling of taxa.

Accurate crop yield estimates are important for governments, farmers, scientists and agribusiness. This paper provides a novel demonstration of the use of freely available Sentinel-2 data to estimate within-field wheat yield variability in a single year. The impact of data resolution and availability on yield estimation is explored using different combinations of input data. This was achieved by combining Sentinel-2 with environmental data (e.g. meteorological, topographical, soil moisture) for different periods throughout the growing season. Yield was estimated using Random Forest (RF) regression models. They were trained and validated using a dataset containing over 8000 points collected by combine harvester yield monitors from 39 wheat fields in the UK. The results demonstrate that it is possible to produce accurate maps of within-field yield variation at 10 m resolution using Sentinel-2 data (RMSE 0.66 t/ha). When combined with environmental data further improvements in accuracy can be obtained (RMSE 0.61 t/ha). We demonstrate that with knowledge of crop-type distribution it is possible to use these models, trained with data from a few fields, to estimate within-field yield variability on a landscape scale. Applying this method gives us a range of crop yield across the landscape of 4.09 to 12.22 t/ha, with a total crop production of approx. 289,000 t.