Senckenberg Research Institute and Natural History Museum Frankfurt/M

Collinearity refers to the non independence of predictor variables, usually in a regression‐type analysis. It is a common feature of any descriptive ecological data set and can be a problem for parameter estimation because it inflates the variance of regression parameters and hence potentially leads to the wrong identification of relevant predictors in a statistical model. Collinearity is a severe problem when a model is trained on data from one region or time, and predicted to another with a different or unknown structure of collinearity. To demonstrate the reach of the problem of collinearity in ecology, we show how relationships among predictors differ between biomes, change over spatial scales and through time. Across disciplines, different approaches to addressing collinearity problems have been developed, ranging from clustering of predictors, threshold‐based pre‐selection, through latent variable methods, to shrinkage and regularisation. Using simulated data with five predictor‐response relationships of increasing complexity and eight levels of collinearity we compared ways to address collinearity with standard multiple regression and machine‐learning approaches. We assessed the performance of each approach by testing its impact on prediction to new data. In the extreme, we tested whether the methods were able to identify the true underlying relationship in a training dataset with strong collinearity by evaluating its performance on a test dataset without any collinearity. We found that methods specifically designed for collinearity, such as latent variable methods and tree based models, did not outperform the traditional GLM and threshold‐based pre‐selection. Our results highlight the value of GLM in combination with penalised methods (particularly ridge) and threshold‐based pre‐selection when omitted variables are considered in the final interpretation. However, all approaches tested yielded degraded predictions under change in collinearity structure and the ‘folk lore’‐thresholds of correlation coefficients between predictor variables of |r| >0.7 was an appropriate indicator for when collinearity begins to severely distort model estimation and subsequent prediction. The use of ecological understanding of the system in pre‐analysis variable selection and the choice of the least sensitive statistical approaches reduce the problems of collinearity, but cannot ultimately solve them.

Abstract High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species 1–4 . To address this issue, the international Genome 10K (G10K) consortium 5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

Hybridization has many and varied impacts on the process of speciation. Hybridization may slow or reverse differentiation by allowing gene flow and recombination. It may accelerate speciation via adaptive introgression or cause near-instantaneous speciation by allopolyploidization. It may have multiple effects at different stages and in different spatial contexts within a single speciation event. We offer a perspective on the context and evolutionary significance of hybridization during speciation, highlighting issues of current interest and debate. In secondary contact zones, it is uncertain if barriers to gene flow will be strengthened or broken down due to recombination and gene flow. Theory and empirical evidence suggest the latter is more likely, except within and around strongly selected genomic regions. Hybridization may contribute to speciation through the formation of new hybrid taxa, whereas introgression of a few loci may promote adaptive divergence and so facilitate speciation. Gene regulatory networks, epigenetic effects and the evolution of selfish genetic material in the genome suggest that the Dobzhansky-Muller model of hybrid incompatibilities requires a broader interpretation. Finally, although the incidence of reinforcement remains uncertain, this and other interactions in areas of sympatry may have knock-on effects on speciation both within and outside regions of hybridization.

Plant traits-the morphological, anatomical, physiological, biochemical and phenological characteristics of plants-determine how plants respond to environmental factors, affect other trophic levels, and influence ecosystem properties and their benefits and detriments to people. Plant trait data thus represent the basis for a vast area of research spanning from evolutionary biology, community and functional ecology, to biodiversity conservation, ecosystem and landscape management, restoration, biogeography and earth system modelling. Since its foundation in 2007, the TRY database of plant traits has grown continuously. It now provides unprecedented data coverage under an open access data policy and is the main plant trait database used by the research community worldwide. Increasingly, the TRY database also supports new frontiers of trait-based plant research, including the identification of data gaps and the subsequent mobilization or measurement of new data. To support this development, in this article we evaluate the extent of the trait data compiled in TRY and analyse emerging patterns of data coverage and representativeness. Best species coverage is achieved for categorical traits-almost complete coverage for 'plant growth form'. However, most traits relevant for ecology and vegetation modelling are characterized by continuous intraspecific variation and trait-environmental relationships. These traits have to be measured on individual plants in their respective environment. Despite unprecedented data coverage, we observe a humbling lack of completeness and representativeness of these continuous traits in many aspects. We, therefore, conclude that reducing data gaps and biases in the TRY database remains a key challenge and requires a coordinated approach to data mobilization and trait measurements. This can only be achieved in collaboration with other initiatives.

Assessing Biodiversity Declines Understanding human impact on biodiversity depends on sound quantitative projection. Pereira et al. (p. 1496 , published online 26 October) review quantitative scenarios that have been developed for four main areas of concern: species extinctions, species abundances and community structure, habitat loss and degradation, and shifts in the distribution of species and biomes. Declines in biodiversity are projected for the whole of the 21st century in all scenarios, but with a wide range of variation. Hoffmann et al. (p. 1503 , published online 26 October) draw on the results of five decades' worth of data collection, managed by the International Union for Conservation of Nature Species Survival Commission. A comprehensive synthesis of the conservation status of the world's vertebrates, based on an analysis of 25,780 species (approximately half of total vertebrate diversity), is presented: Approximately 20% of all vertebrate species are at risk of extinction in the wild, and 11% of threatened birds and 17% of threatened mammals have moved closer to extinction over time. Despite these trends, overall declines would have been significantly worse in the absence of conservation actions.

Biodiversity is unevenly distributed on Earth and hotspots of biodiversity are often associated with areas that have undergone orogenic activity during recent geological history (i.e. tens of millions of years). Understanding the underlying processes that have driven the accumulation of species in some areas and not in others may help guide prioritization in conservation and may facilitate forecasts on ecosystem services under future climate conditions. Consequently, the study of the origin and evolution of biodiversity in mountain systems has motivated growing scientific interest. Despite an increasing number of studies, the origin and evolution of diversity hotspots associated with the Qinghai-Tibetan Plateau (QTP) remains poorly understood. We review literature related to the diversification of organisms linked to the uplift of the QTP. To promote hypothesis-based research, we provide a geological and palaeoclimatic scenario for the region of the QTP and argue that further studies would benefit from providing a complete set of complementary analyses (molecular dating, biogeographic, and diversification rates analyses) to test for a link between organismic diversification and past geological and climatic changes in this region. In general, we found that the contribution of biological interchange between the QTP and other hotspots of biodiversity has not been sufficiently studied to date. Finally, we suggest that the biological consequences of the uplift of the QTP would be best understood using a meta-analysis approach, encompassing studies on a variety of organisms (plants and animals) from diverse habitats (forests, meadows, rivers), and thermal belts (montane, subalpine, alpine, nival). Since the species diversity in the QTP region is better documented for some organismic groups than for others, we suggest that baseline taxonomic work should be promoted.

Our ability to correlate biological evolution with climate change, geological evolution, and other historical patterns is essential to understanding the processes that shape biodiversity. Combining data from the fossil record with molecular phylogenetics represents an exciting synthetic approach to this challenge. The first molecular divergence dating analysis (Zuckerkandl and Pauling 1962) was based on a measure of the amino acid differences in the hemoglobin molecule, with replacement rates established (calibrated) using paleontological age estimates from textbooks (e.g., Dodson 1960). Since that time, the amount of molecular sequence data has increased dramatically, affording ever-greater opportunities to apply molecular divergence approaches to fundamental problems in evolutionary biology.

Patterns of diversification and timing of evolution within Neoaves, which includes almost 95% of all bird species, are virtually unknown. On the other hand, molecular data consistently indicate a Cretaceous origin of many neoavian lineages and the fossil record seems to support an Early Tertiary diversification. Here, we present the first well-resolved molecular phylogeny for Neoaves, together with divergence time estimates calibrated with a large number of stratigraphically and phylogenetically well-documented fossils. Our study defines several well-supported clades within Neoaves. The calibration results suggest that Neoaves, after an initial split from Galloanseres in Mid-Cretaceous, diversified around or soon after the K/T boundary. Our results thus do not contradict palaeontological data and show that there is no solid molecular evidence for an extensive pre-Tertiary radiation of Neoaves.

Genetic diversity provides the basic substrate for evolution, yet few studies assess the impacts of global climate change (GCC) on intraspecific genetic variation. In this review, we highlight the importance of incorporating neutral and non-neutral genetic diversity when assessing the impacts of GCC, for example, in studies that aim to predict the future distribution and fate of a species or ecological community. Specifically, we address the following questions: Why study the effects of GCC on intraspecific genetic diversity? How does GCC affect genetic diversity? How is the effect of GCC on genetic diversity currently studied? Where is potential for future research? For each of these questions, we provide a general background and highlight case studies across the animal, plant and microbial kingdoms. We further discuss how cryptic diversity can affect GCC assessments, how genetic diversity can be integrated into studies that aim to predict species' responses on GCC and how conservation efforts related to GCC can incorporate and profit from inclusion of genetic diversity assessments. We argue that studying the fate of intraspecifc genetic diversity is an indispensable and logical venture if we are to fully understand the consequences of GCC on biodiversity on all levels.

Vocalizations of anuran amphibians have received much attention in studies of behavioral ecology and physiology, but also provide informative characters for identifying and delimiting species. We here review the terminology and variation of frog calls from a perspective of integrative taxonomy, and provide hands-on protocols for recording, analyzing, comparing, interpreting and describing these sounds. Our focus is on advertisement calls, which serve as premating isolation mechanisms and, therefore, convey important taxonomic information. We provide recommendations for terminology of frog vocalizations, with call, note and pulse being the fundamental subunits to be used in descriptions and comparisons. However, due to the complexity and diversity of these signals, an unequivocal application of the terms call and note can be challenging. We therefore provide two coherent concepts that either follow a note-centered approach (defining uninterrupted units of sound as notes, and their entirety as call) or a call-centered approach (defining uninterrupted units as call whenever they are separated by long silent intervals) in terminology. Based on surveys of literature, we show that numerous call traits can be highly variable within and between individuals of one species. Despite idiosyncrasies of species and higher taxa, the duration of calls or notes, pulse rate within notes, and number of pulses per note appear to be more static within individuals and somewhat less affected by temperature. Therefore, these variables might often be preferable as taxonomic characters over call rate or note rate, which are heavily influenced by various factors. Dominant frequency is also comparatively static and only weakly affected by temperature, but depends strongly on body size. As with other taxonomic characters, strong call divergence is typically indicative of species-level differences, whereas call similarities of two populations are no evidence for them being conspecific. Taxonomic conclusions can especially be drawn when the general advertisement call structure of two candidate species is radically different and qualitative call differences are thus observed. On the other hand, quantitative differences in call traits might substantially vary within and among conspecific populations, and require careful evaluation and analysis. We provide guidelines for the taxonomic interpretation of advertisement call differences in sympatric and allopatric situations, and emphasize the need for an integrative use of multiple datasets (bio-acoustics, morphology, genetics), particularly for allopatric scenarios. We show that small-sized frogs often emit calls with frequency components in the ultrasound spectrum, although it is unlikely that these high frequencies are of biological relevance for the majority of them, and we illustrate that detection of upper harmonics depends also on recording distance because higher frequencies are attenuated more strongly. Bioacoustics remains a prime approach in integrative taxonomy of anurans if uncertainty due to possible intraspecific variation and technical artifacts is adequately considered and acknowledged.

The aim of this study was to assess potential candidate gene regions and corresponding universal primer pairs as secondary DNA barcodes for the fungal kingdom, additional to ITS rDNA as primary barcode. Amplification efficiencies of 14 (partially) universal primer pairs targeting eight genetic markers were tested across > 1 500 species (1 931 strains or specimens) and the outcomes of almost twenty thousand (19 577) polymerase chain reactions were evaluated. We tested several well-known primer pairs that amplify: i) sections of the nuclear ribosomal RNA gene large subunit (D1-D2 domains of 26/28S); ii) the complete internal transcribed spacer region (ITS1/2); iii) partial β -tubulin II (TUB2); iv) γ-actin (ACT); v) translation elongation factor 1-α (TEF1α); and vi) the second largest subunit of RNA-polymerase II (partial RPB2, section 5-6). Their PCR efficiencies were compared with novel candidate primers corresponding to: i) the fungal-specific translation elongation factor 3 (TEF3); ii) a small ribosomal protein necessary for t-RNA docking; iii) the 60S L10 (L1) RP; iv) DNA topoisomerase I (TOPI); v) phosphoglycerate kinase (PGK); vi) hypothetical protein LNS2; and vii) alternative sections of TEF1α. Results showed that several gene sections are accessible to universal primers (or primers universal for phyla) yielding a single PCR-product. Barcode gap and multi-dimensional scaling analyses revealed that some of the tested candidate markers have universal properties providing adequate infra- and inter-specific variation that make them attractive barcodes for species identification. Among these gene sections, a novel high fidelity primer pair for TEF1α, already widely used as a phylogenetic marker in mycology, has potential as a supplementary DNA barcode with superior resolution to ITS. Both TOPI and PGK show promise for the Ascomycota, while TOPI and LNS2 are attractive for the Pucciniomycotina, for which universal primers for ribosomal subunits often fail.

The bioassessment of aquatic ecosystems is currently based on various biotic indices that use the occurrence and/or abundance of selected taxonomic groups to define ecological status. These conventional indices have some limitations, often related to difficulties in morphological identification of bioindicator taxa. Recent development of DNA barcoding and metabarcoding could potentially alleviate some of these limitations, by using DNA sequences instead of morphology to identify organisms and to characterize a given ecosystem. In this paper, we review the structure of conventional biotic indices, and we present the results of pilot metabarcoding studies using environmental DNA to infer biotic indices. We discuss the main advantages and pitfalls of metabarcoding approaches to assess parameters such as richness, abundance, taxonomic composition and species ecological values, to be used for calculation of biotic indices. We present some future developments to fully exploit the potential of metabarcoding data and improve the accuracy and precision of their analysis. We also propose some recommendations for the future integration of DNA metabarcoding to routine biomonitoring programs.

Freshwater mussels of the Order Unionida provide important ecosystem functions and services, yet many of their populations are in decline. We comprehensively review the status of the 16 currently recognized species in Europe, collating for the first time their life-history traits, distribution, conservation status, habitat preferences, and main threats in order to suggest future management actions. In northern, central, and eastern Europe, a relatively homogeneous species composition is found in most basins. In southern Europe, despite the lower species richness, spatially restricted species make these basins a high conservation priority. Information on freshwater mussels in Europe is unevenly distributed with considerable differences in data quality and quantity among countries and species. To make conservation more effective in the future, we suggest greater international cooperation using standardized protocols and methods to monitor and manage European freshwater mussel diversity. Such an approach will not only help conserve this vulnerable group but also, through the protection of these important organisms, will offer wider benefits to freshwater ecosystems.

PREMISE: Bromeliaceae form a large, ecologically diverse family of angiosperms native to the New World. We use a bromeliad phylogeny based on eight plastid regions to analyze relationships within the family, test a new, eight-subfamily classification, infer the chronology of bromeliad evolution and invasion of different regions, and provide the basis for future analyses of trait evolution and rates of diversification. METHODS: We employed maximum-parsimony, maximum-likelihood, and Bayesian approaches to analyze 9341 aligned bases for four outgroups and 90 bromeliad species representing 46 of 58 described genera. We calibrate the resulting phylogeny against time using penalized likelihood applied to a monocot-wide tree based on plastid ndhF sequences and use it to analyze patterns of geographic spread using parsimony, Bayesian inference, and the program S-DIVA. RESULTS: Bromeliad subfamilies are related to each other as follows: (Brocchinioideae, (Lindmanioideae, (Tillandsioideae, (Hechtioideae, (Navioideae, (Pitcairnioideae, (Puyoideae, Bromelioideae))))))). Bromeliads arose in the Guayana Shield ca. 100 million years ago (Ma), spread centrifugally in the New World beginning ca. 16-13 Ma, and dispersed to West Africa ca. 9.3 Ma. Modern lineages began to diverge from each other roughly 19 Ma. CONCLUSIONS: Nearly two-thirds of extant bromeliads belong to two large radiations: the core tillandsioids, originating in the Andes ca. 14.2 Ma, and the Brazilian Shield bromelioids, originating in the Serro do Mar and adjacent regions ca. 9.1 Ma.

Local biodiversity trends over time are likely to be decoupled from global trends, as local processes may compensate or counteract global change. We analyze 161 long-term biological time series (15-91 years) collected across Europe, using a comprehensive dataset comprising ~6,200 marine, freshwater and terrestrial taxa. We test whether (i) local long-term biodiversity trends are consistent among biogeoregions, realms and taxonomic groups, and (ii) changes in biodiversity correlate with regional climate and local conditions. Our results reveal that local trends of abundance, richness and diversity differ among biogeoregions, realms and taxonomic groups, demonstrating that biodiversity changes at local scale are often complex and cannot be easily generalized. However, we find increases in richness and abundance with increasing temperature and naturalness as well as a clear spatial pattern in changes in community composition (i.e. temporal taxonomic turnover) in most biogeoregions of Northern and Eastern Europe.

Ning Li and colleagues report the whole-genome sequence of the duck, Anas platyrhynchos, a natural host of avian influenza viruses. They examine host response to infection by comparing the lung transcriptomes of ducks that were infected with influenza A viruses. The duck (Anas platyrhynchos) is one of the principal natural hosts of influenza A viruses. We present the duck genome sequence and perform deep transcriptome analyses to investigate immune-related genes. Our data indicate that the duck possesses a contractive immune gene repertoire, as in chicken and zebra finch, and this repertoire has been shaped through lineage-specific duplications. We identify genes that are responsive to influenza A viruses using the lung transcriptomes of control ducks and ones that were infected with either a highly pathogenic (A/duck/Hubei/49/05) or a weakly pathogenic (A/goose/Hubei/65/05) H5N1 virus. Further, we show how the duck's defense mechanisms against influenza infection have been optimized through the diversification of its β-defensin and butyrophilin-like repertoires. These analyses, in combination with the genomic and transcriptomic data, provide a resource for characterizing the interaction between host and influenza viruses.

Abstract River networks are hierarchical dendritic habitats embedded within the terrestrial landscape, with varying connectivity between sites depending on their positions along the network. This physical organisation influences the dispersal of organisms, which ultimately affects metacommunity dynamics and biodiversity patterns. We provide a conceptual synthesis of the role of river networks in structuring metacommunities in relation to dispersal processes in riverine ecosystems. We explore where the river network best explains observed metacommunity structure compared to other measurements of physical connectivity. We mostly focus on invertebrates, but also consider other taxonomic groups, including microbes, fishes, plants, and amphibians. Synthesising studies that compared multiple spatial distance metrics, we found that the importance of the river network itself in explaining metacommunity patterns depended on a variety of factors, including dispersal mode (aquatic versus aerial versus terrestrial) and landscape type (arid versus mesic), as well as location‐specific factors, such as network connectivity, land use, topographic heterogeneity, and biotic interactions. The river network appears to be less important for strong aerial dispersers and insects in arid systems than for other groups and biomes, but there is considerable variability. Borrowing from other literature, particularly landscape genetics, we developed a conceptual model that predicts that the explanatory power of the river network peaks in mesic systems for obligate aquatic dispersers. We propose directions of future avenues of research, including the use of manipulative field and laboratory experiments that test metacommunity theory in river networks. While field and laboratory experiments have their own benefits and drawbacks (e.g. reality, control, cost), both are powerful approaches for understanding the mechanisms structuring metacommunities, by teasing apart dispersal and niche‐related factors. Finally, improving our knowledge of dispersal in river networks will benefit from expanding the breadth of cost‐distance modelling to better infer dispersal from observational data; an improved understanding of life‐history strategies rather than relying on independent traits; exploring individual‐level variation in dispersal through detailed genetic studies; detailed studies on fine‐scale environmental (e.g. daily hydrology) and organismal spatiotemporal variability; and synthesising comparative, experimental, and theoretical work. Expanding in these areas will help to push the current state of the science from a largely pattern‐detection mode into a new phase of more mechanistically driven research.

Abstract. Global-scale assessments of freshwater fluxes and storages by hydrological models under historic climate conditions are subject to a variety of uncertainties. Using the global hydrological model WaterGAP (Water – Global Assessment and Prognosis) 2.2, we investigated the sensitivity of simulated freshwater fluxes and water storage variations to five major sources of uncertainty: climate forcing, land cover input, model structure/refinements, consideration of human water use and calibration (or no calibration) against observed mean river discharge. In a modeling experiment, five variants of the standard version of WaterGAP 2.2 were generated that differed from the standard version only regarding the investigated source of uncertainty. The basin-specific calibration approach for WaterGAP was found to have the largest effect on grid cell fluxes as well as on global AET (actual evapotranspiration) and discharge into oceans for the period 1971–2000. Regarding grid cell fluxes, climate forcing ranks second before land cover input. Global water storage trends are most sensitive to model refinements (mainly modeling of groundwater depletion) and consideration of human water use. The best fit to observed time series of monthly river discharge or discharge seasonality is obtained with the standard WaterGAP 2.2 model version which is calibrated and driven by daily reanalysis-based WFD/WFDEI (combination of Watch Forcing Data based on ERA40 and Watch Forcing Data based on ERA-Interim) climate data. Discharge computed by a calibrated model version using monthly CRU TS (Climate Research Unit time-series) 3.2 and GPCC (Global Precipitation Climatology Center) v6 climate input reduced the fit to observed discharge for most stations. Taking into account uncertainties of climate and land cover data, global 1971–2000 discharge into oceans and inland sinks ranges between 40 000 and 42 000 km3 yr−1. Global actual evapotranspiration, with 70 000 km3 yr−1, is rather unaffected by climate and land cover uncertainties. Human water use reduced river discharge by 1000 km3 yr−1, such that global renewable water resources are estimated to range between 41 000 and 43 000 km3 yr−1. The climate data sets WFD (available until 2001) and WFDEI (starting in 1979) were found to be inconsistent with respect to shortwave radiation data, resulting in strongly different actual evapotranspiration. Global assessments of freshwater fluxes and storages would therefore benefit from the development of a global data set of consistent daily climate forcing from 1900 to present.

Abstract Mutated isocitrate dehydrogenase 1 ( IDH1 ) defines a molecularly distinct subtype of diffuse glioma 1–3 . The most common IDH1 mutation in gliomas affects codon 132 and encodes IDH1(R132H), which harbours a shared clonal neoepitope that is presented on major histocompatibility complex (MHC) class II 4,5 . An IDH1(R132H)-specific peptide vaccine (IDH1-vac) induces specific therapeutic T helper cell responses that are effective against IDH1(R132H) + tumours in syngeneic MHC-humanized mice 4,6–8 . Here we describe a multicentre, single-arm, open-label, first-in-humans phase I trial that we carried out in 33 patients with newly diagnosed World Health Organization grade 3 and 4 IDH1(R132H) + astrocytomas (Neurooncology Working Group of the German Cancer Society trial 16 (NOA16), ClinicalTrials.gov identifier NCT02454634). The trial met its primary safety endpoint, with vaccine-related adverse events restricted to grade 1. Vaccine-induced immune responses were observed in 93.3% of patients across multiple MHC alleles. Three-year progression-free and death-free rates were 0.63 and 0.84, respectively. Patients with immune responses showed a two-year progression-free rate of 0.82. Two patients without an immune response showed tumour progression within two years of first diagnosis. A mutation-specificity score that incorporates the duration and level of vaccine-induced IDH1(R132H)-specific T cell responses was associated with intratumoral presentation of the IDH1(R132H) neoantigen in pre-treatment tumour tissue. There was a high frequency of pseudoprogression, which indicates intratumoral inflammatory reactions. Pseudoprogression was associated with increased vaccine-induced peripheral T cell responses. Combined single-cell RNA and T cell receptor sequencing showed that tumour-infiltrating CD40LG + and CXCL13 + T helper cell clusters in a patient with pseudoprogression were dominated by a single IDH1(R132H)-reactive T cell receptor.

In order to test accretion simulations as well as planetary differentiation scenarios, we have integrated a multistage core–mantle differentiation model with N-body accretion simulations. Impacts between embryos and planetesimals are considered to result in magma ocean formation and episodes of core formation. The core formation model combines rigorous chemical mass balance with metal–silicate element partitioning data and requires that the bulk compositions of all starting embryos and planetesimals are defined as a function of their heliocentric distances of origin. To do this, we assume that non-volatile elements are present in Solar System (CI) relative abundances in all bodies and that oxygen and H2O contents are the main compositional variables. The primary constraint on the combined model is the composition of the Earth’s primitive mantle. In addition, we aim to reproduce the composition of the martian mantle and the mass fractions of the metallic cores of Earth and Mars. The model is refined by least squares minimization with up to five fitting parameters that consist of the metal–silicate equilibration pressure and 1–4 parameters that define the starting compositions of primitive bodies. This integrated model has been applied to six Grand Tack N-body accretion simulations. Investigations of a broad parameter space indicate that: (1) accretion of Earth was heterogeneous, (2) metal–silicate equilibration pressures increase as accretion progresses and are, on average, 60–70% of core–mantle boundary pressures at the time of each impact, and (3) a large fraction (70–100%) of the metal of impactor cores equilibrates with a small fraction of the silicate mantles of proto-planets during each core formation event. Results are highly sensitive to the compositional model for the primitive starting bodies and several accretion/core-formation models can thus be excluded. Acceptable fits to the Earth’s mantle composition are obtained only when bodies that originated close to the Sun, at <0.9–1.2 AU, are highly reduced and those from beyond this distance are increasingly oxidized. Reasonable concentrations of H2O in Earth’s mantle are obtained when bodies originating from beyond 6–7 AU contain 20 wt% water ice (icy bodies that originated between the snow line and this distance did not contribute to Earth’s accretion because they were swept up by Jupiter and Saturn). In the six models examined, water is added to the Earth mainly after 60–80% of its final mass has accreted. The compositional evolution of the mantles of Venus and Mars are also constrained by the model. The FeO content of the martian mantle depends critically on the heliocentric distance at which the Mars-forming embryo originated. Finally, the Earth’s core is predicted to contain 8–9 wt% silicon, 2–4 wt% oxygen and 10–60 ppm hydrogen, whereas the martian core is predicted to contain low concentrations (<1 wt%) of Si and O.