Senckenberg - Leibniz Institution for Biodiversity and Earth System Research

Six DNA regions were evaluated as potential DNA barcodes for Fungi, the second largest kingdom of eukaryotic life, by a multinational, multilaboratory consortium. The region of the mitochondrial cytochrome c oxidase subunit 1 used as the animal barcode was excluded as a potential marker, because it is difficult to amplify in fungi, often includes large introns, and can be insufficiently variable. Three subunits from the nuclear ribosomal RNA cistron were compared together with regions of three representative protein-coding genes (largest subunit of RNA polymerase II, second largest subunit of RNA polymerase II, and minichromosome maintenance protein). Although the protein-coding gene regions often had a higher percent of correct identification compared with ribosomal markers, low PCR amplification and sequencing success eliminated them as candidates for a universal fungal barcode. Among the regions of the ribosomal cistron, the internal transcribed spacer (ITS) region has the highest probability of successful identification for the broadest range of fungi, with the most clearly defined barcode gap between inter- and intraspecific variation. The nuclear ribosomal large subunit, a popular phylogenetic marker in certain groups, had superior species resolution in some taxonomic groups, such as the early diverging lineages and the ascomycete yeasts, but was otherwise slightly inferior to the ITS. The nuclear ribosomal small subunit has poor species-level resolution in fungi. ITS will be formally proposed for adoption as the primary fungal barcode marker to the Consortium for the Barcode of Life, with the possibility that supplementary barcodes may be developed for particular narrowly circumscribed taxonomic groups.

Modern attempts to produce biogeographic maps focus on the distribution of species, and the maps are typically drawn without phylogenetic considerations. Here, we generate a global map of zoogeographic regions by combining data on the distributions and phylogenetic relationships of 21,037 species of amphibians, birds, and mammals. We identify 20 distinct zoogeographic regions, which are grouped into 11 larger realms. We document the lack of support for several regions previously defined based on distributional data and show that spatial turnover in the phylogenetic composition of vertebrate assemblages is higher in the Southern than in the Northern Hemisphere. We further show that the integration of phylogenetic information provides valuable insight on historical relationships among regions, permitting the identification of evolutionarily unique regions of the world.

Animal movement is fundamental for ecosystem functioning and species survival, yet the effects of the anthropogenic footprint on animal movements have not been estimated across species. Using a unique GPS-tracking database of 803 individuals across 57 species, we found that movements of mammals in areas with a comparatively high human footprint were on average one-half to one-third the extent of their movements in areas with a low human footprint. We attribute this reduction to behavioral changes of individual animals and to the exclusion of species with long-range movements from areas with higher human impact. Global loss of vagility alters a key ecological trait of animals that affects not only population persistence but also ecosystem processes such as predator-prey interactions, nutrient cycling, and disease transmission.

The novel coronavirus disease 2019 (COVID-19) has become a pandemic affecting health and wellbeing globally. In addition to the physical health, economic, and social implications, the psychological impacts of this pandemic are increasingly being reported in the scientific literature. This narrative review reflected on scholarly articles on the epidemiology of mental health problems in COVID-19. The current literature suggests that people affected by COVID-19 may have a high burden of mental health problems, including depression, anxiety disorders, stress, panic attack, irrational anger, impulsivity, somatization disorder, sleep disorders, emotional disturbance, posttraumatic stress symptoms, and suicidal behavior. Moreover, several factors associated with mental health problems in COVID-19 are found, which include age, gender, marital status, education, occupation, income, place of living, close contact with people with COVID-19, comorbid physical and mental health problems, exposure to COVID-19 related news and social media, coping styles, stigma, psychosocial support, health communication, confidence in health services, personal protective measures, risk of contracting COVID-19, and perceived likelihood of survival. Furthermore, the epidemiological distribution of mental health problems and associated factors were heterogeneous among the general public, COVID-19 patients, and healthcare providers. The current evidence suggests that a psychiatric epidemic is cooccurring with the COVID-19 pandemic, which necessitates the attention of the global health community. Future epidemiological studies should emphasize on psychopathological variations and temporality of mental health problems in different populations. Nonetheless, multipronged interventions should be developed and adopted to address the existing psychosocial challenges and promote mental health amid the COVID-19 pandemic.

The use of phylogenies in ecology is increasingly common and has broadened our understanding of biological diversity. Ecological sub-disciplines, particularly conservation, community ecology and macroecology, all recognize the value of evolutionary relationships but the resulting development of phylogenetic approaches has led to a proliferation of phylogenetic diversity metrics. The use of many metrics across the sub-disciplines hampers potential meta-analyses, syntheses, and generalizations of existing results. Further, there is no guide for selecting the appropriate metric for a given question, and different metrics are frequently used to address similar questions. To improve the choice, application, and interpretation of phylo-diversity metrics, we organize existing metrics by expanding on a unifying framework for phylogenetic information. Generally, questions about phylogenetic relationships within or between assemblages tend to ask three types of question: how much; how different; or how regular? We show that these questions reflect three dimensions of a phylogenetic tree: richness, divergence, and regularity. We classify 70 existing phylo-diversity metrics based on their mathematical form within these three dimensions and identify 'anchor' representatives: for α-diversity metrics these are PD (Faith's phylogenetic diversity), MPD (mean pairwise distance), and VPD (variation of pairwise distances). By analysing mathematical formulae and using simulations, we use this framework to identify metrics that mix dimensions, and we provide a guide to choosing and using the most appropriate metrics. We show that metric choice requires connecting the research question with the correct dimension of the framework and that there are logical approaches to selecting and interpreting metrics. The guide outlined herein will help researchers navigate the current jungle of indices.

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.

Global change, especially land-use intensification, affects human well-being by impacting the delivery of multiple ecosystem services (multifunctionality). However, whether biodiversity loss is a major component of global change effects on multifunctionality in real-world ecosystems, as in experimental ones, remains unclear. Therefore, we assessed biodiversity, functional composition and 14 ecosystem services on 150 agricultural grasslands differing in land-use intensity. We also introduce five multifunctionality measures in which ecosystem services were weighted according to realistic land-use objectives. We found that indirect land-use effects, i.e. those mediated by biodiversity loss and by changes to functional composition, were as strong as direct effects on average. Their strength varied with land-use objectives and regional context. Biodiversity loss explained indirect effects in a region of intermediate productivity and was most damaging when land-use objectives favoured supporting and cultural services. In contrast, functional composition shifts, towards fast-growing plant species, strongly increased provisioning services in more inherently unproductive grasslands.

Abstract. The LPJ-GUESS dynamic vegetation model uniquely combines an individual- and patch-based representation of vegetation dynamics with ecosystem biogeochemical cycling from regional to global scales. We present an updated version that includes plant and soil N dynamics, analysing the implications of accounting for C–N interactions on predictions and performance of the model. Stand structural dynamics and allometric scaling of tree growth suggested by global databases of forest stand structure and development were well reproduced by the model in comparison to an earlier multi-model study. Accounting for N cycle dynamics improved the goodness of fit for broadleaved forests. N limitation associated with low N-mineralisation rates reduces productivity of cold-climate and dry-climate ecosystems relative to mesic temperate and tropical ecosystems. In a model experiment emulating free-air CO2 enrichment (FACE) treatment for forests globally, N limitation associated with low N-mineralisation rates of colder soils reduces CO2 enhancement of net primary production (NPP) for boreal forests, while some temperate and tropical forests exhibit increased NPP enhancement. Under a business-as-usual future climate and emissions scenario, ecosystem C storage globally was projected to increase by ca. 10%; additional N requirements to match this increasing ecosystem C were within the high N supply limit estimated on stoichiometric grounds in an earlier study. Our results highlight the importance of accounting for C–N interactions in studies of global terrestrial N cycling, and as a basis for understanding mechanisms on local scales and in different regional contexts.

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.

Human activities are fundamentally altering biodiversity. Projections of declines at the global scale are contrasted by highly variable trends at local scales, suggesting that biodiversity change may be spatially structured. Here, we examined spatial variation in species richness and composition change using more than 50,000 biodiversity time series from 239 studies and found clear geographic variation in biodiversity change. Rapid compositional change is prevalent, with marine biomes exceeding and terrestrial biomes trailing the overall trend. Assemblage richness is not changing on average, although locations exhibiting increasing and decreasing trends of up to about 20% per year were found in some marine studies. At local scales, widespread compositional reorganization is most often decoupled from richness change, and biodiversity change is strongest and most variable in the oceans.

Memory is critical to understanding animal movement but has proven challenging to study. Advances in animal tracking technology, theoretical movement models and cognitive sciences have facilitated research in each of these fields, but also created a need for synthetic examination of the linkages between memory and animal movement. Here, we draw together research from several disciplines to understand the relationship between animal memory and movement processes. First, we frame the problem in terms of the characteristics, costs and benefits of memory as outlined in psychology and neuroscience. Next, we provide an overview of the theories and conceptual frameworks that have emerged from behavioural ecology and animal cognition. Third, we turn to movement ecology and summarise recent, rapid developments in the types and quantities of available movement data, and in the statistical measures applicable to such data. Fourth, we discuss the advantages and interrelationships of diverse modelling approaches that have been used to explore the memory-movement interface. Finally, we outline key research challenges for the memory and movement communities, focusing on data needs and mathematical and computational challenges. We conclude with a roadmap for future work in this area, outlining axes along which focused research should yield rapid progress.

Quantifying animals' home ranges is a key problem in ecology and has important conservation and wildlife management applications. Kernel density estimation (KDE) is a workhorse technique for range delineation problems that is both statistically efficient and nonparametric. KDE assumes that the data are independent and identically distributed (IID). However, animal tracking data, which are routinely used as inputs to KDEs, are inherently autocorrelated and violate this key assumption. As we demonstrate, using realistically autocorrelated data in conventional KDEs results in grossly underestimated home ranges. We further show that the performance of conventional KDEs actually degrades as data quality improves, because autocorrelation strength increases as movement paths become more finely resolved. To remedy these flaws with the traditional KDE method, we derive an autocorrelated KDE (AKDE) from first principles to use autocorrelated data, making it perfectly suited for movement data sets. We illustrate the vastly improved performance of AKDE using analytical arguments, relocation data from Mongolian gazelles, and simulations based upon the gazelle's observed movement process. By yielding better minimum area estimates for threatened wildlife populations, we believe that future widespread use of AKDE will have significant impact on ecology and conservation biology.

Microorganisms synthesize small bioactive compounds through their secondary or specialized metabolism. Those compounds play an important role in microbial interactions and soil health, but are also crucial for the development of pharmaceuticals or agrochemicals. Over the past decades, advancements in genome sequencing have enabled the identification of large numbers of biosynthetic gene clusters directly from microbial genomes. Since its inception in 2011, antiSMASH (https://antismash.secondarymetabolites.org/), has become the leading tool for detecting and characterizing these gene clusters in bacteria and fungi. This paper introduces version 8 of antiSMASH, which has increased the number of detectable cluster types from 81 to 101, and has improved analysis support for terpenoids and tailoring enzymes, as well as improvements in the analysis of modular enzymes like polyketide synthases and nonribosomal peptide synthetases. These modifications keep antiSMASH up-to-date with developments in the field and extend its overall predictive capabilities for natural product genome mining.

BACKGROUND: For more than three decades, the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) has provided a framework to quantify health loss due to diseases, injuries, and associated risk factors. This paper presents GBD 2023 findings on disease and injury burden and risk-attributable health loss, offering a global audit of the state of world health to inform public health priorities. This work captures the evolving landscape of health metrics across age groups, sexes, and locations, while reflecting on the remaining post-COVID-19 challenges to achieving our collective global health ambitions. METHODS: The GBD 2023 combined analysis estimated years lived with disability (YLDs), years of life lost (YLLs), and disability-adjusted life-years (DALYs) for 375 diseases and injuries, and risk-attributable burden associated with 88 modifiable risk factors. Of the more than 310 000 total data sources used for all GBD 2023 (about 30% of which were new to this estimation round), more than 120 000 sources were used for estimation of disease and injury burden and 59 000 for risk factor estimation, and included vital registration systems, surveys, disease registries, and published scientific literature. Data were analysed using previously established modelling approaches, such as disease modelling meta-regression version 2.1 (DisMod-MR 2.1) and comparative risk assessment methods. Diseases and injuries were categorised into four levels on the basis of the established GBD cause hierarchy, as were risk factors using the GBD risk hierarchy. Estimates stratified by age, sex, location, and year from 1990 to 2023 were focused on disease-specific time trends over the 2010-23 period and presented as counts (to three significant figures) and age-standardised rates per 100 000 person-years (to one decimal place). For each measure, 95% uncertainty intervals [UIs] were calculated with the 2·5th and 97·5th percentile ordered values from a 250-draw distribution. FINDINGS: Total numbers of global DALYs grew 6·1% (95% UI 4·0-8·1), from 2·64 billion (2·46-2·86) in 2010 to 2·80 billion (2·57-3·08) in 2023, but age-standardised DALY rates, which account for population growth and ageing, decreased by 12·6% (11·0-14·1), revealing large long-term health improvements. Non-communicable diseases (NCDs) contributed 1·45 billion (1·31-1·61) global DALYs in 2010, increasing to 1·80 billion (1·63-2·03) in 2023, alongside a concurrent 4·1% (1·9-6·3) reduction in age-standardised rates. Based on DALY counts, the leading level 3 NCDs in 2023 were ischaemic heart disease (193 million [176-209] DALYs), stroke (157 million [141-172]), and diabetes (90·2 million [75·2-107]), with the largest increases in age-standardised rates since 2010 occurring for anxiety disorders (62·8% [34·0-107·5]), depressive disorders (26·3% [11·6-42·9]), and diabetes (14·9% [7·5-25·6]). Remarkable health gains were made for communicable, maternal, neonatal, and nutritional (CMNN) diseases, with DALYs falling from 874 million (837-917) in 2010 to 681 million (642-736) in 2023, and a 25·8% (22·6-28·7) reduction in age-standardised DALY rates. During the COVID-19 pandemic, DALYs due to CMNN diseases rose but returned to pre-pandemic levels by 2023. From 2010 to 2023, decreases in age-standardised rates for CMNN diseases were led by rate decreases of 49·1% (32·7-61·0) for diarrhoeal diseases, 42·9% (38·0-48·0) for HIV/AIDS, and 42·2% (23·6-56·6) for tuberculosis. Neonatal disorders and lower respiratory infections remained the leading level 3 CMNN causes globally in 2023, although both showed notable rate decreases from 2010, declining by 16·5% (10·6-22·0) and 24·8% (7·4-36·7), respectively. Injury-related age-standardised DALY rates decreased by 15·6% (10·7-19·8) over the same period. Differences in burden due to NCDs, CMNN diseases, and injuries persisted across age, sex, time, and location. Based on our risk analysis, nearly 50% (1·27 billion [1·18-1·38]) of the roughly 2·80 billion total global DALYs in 2023 were attributable to the 88 risk factors analysed in GBD. Globally, the five level 3 risk factors contributing the highest proportion of risk-attributable DALYs were high systolic blood pressure (SBP), particulate matter pollution, high fasting plasma glucose (FPG), smoking, and low birthweight and short gestation-with high SBP accounting for 8·4% (6·9-10·0) of total DALYs. Of the three overarching level 1 GBD risk factor categories-behavioural, metabolic, and environmental and occupational-risk-attributable DALYs rose between 2010 and 2023 only for metabolic risks, increasing by 30·7% (24·8-37·3); however, age-standardised DALY rates attributable to metabolic risks decreased by 6·7% (2·0-11·0) over the same period. For all but three of the 25 leading level 3 risk factors, age-standardised rates dropped between 2010 and 2023-eg, declining by 54·4% (38·7-65·3) for unsafe sanitation, 50·5% (33·3-63·1) for unsafe water source, and 45·2% (25·6-72·0) for no access to handwashing facility, and by 44·9% (37·3-53·5) for child growth failure. The three leading level 3 risk factors for which age-standardised attributable DALY rates rose were high BMI (10·5% [0·1 to 20·9]), drug use (8·4% [2·6 to 15·3]), and high FPG (6·2% [-2·7 to 15·6]; non-significant). INTERPRETATION: Our findings underscore the complex and dynamic nature of global health challenges. Since 2010, there have been large decreases in burden due to CMNN diseases and many environmental and behavioural risk factors, juxtaposed with sizeable increases in DALYs attributable to metabolic risk factors and NCDs in growing and ageing populations. This long-observed consequence of the global epidemiological transition was only temporarily interrupted by the COVID-19 pandemic. The substantially decreasing CMNN disease burden, despite the 2008 global financial crisis and pandemic-related disruptions, is one of the greatest collective public health successes known. However, these achievements are at risk of being reversed due to major cuts to development assistance for health globally, the effects of which will hit low-income countries with high burden the hardest. Without sustained investment in evidence-based interventions and policies, progress could stall or reverse, leading to widespread human costs and geopolitical instability. Moreover, the rising NCD burden necessitates intensified efforts to mitigate exposure to leading risk factors-eg, air pollution, smoking, and metabolic risks, such as high SBP, BMI, and FPG-including policies that promote food security, healthier diets, physical activity, and equitable and expanded access to potential treatments, such as GLP-1 receptor agonists. Decisive, coordinated action is needed to address long-standing yet growing health challenges, including depressive and anxiety disorders. Yet this can be only part of the solution. Our response to the NCD syndemic-the complex interaction of multiple health risks, social determinants, and systemic challenges-will define the future landscape of global health. To ensure human wellbeing, economic stability, and social equity, global action to sustain and advance health gains must prioritise reducing disparities by addressing socioeconomic and demographic determinants, ensuring equitable health-care access, tackling malnutrition, strengthening health systems, and improving vaccination coverage. We live in times of great opportunity. FUNDING: Gates Foundation and Bloomberg Philanthropies.

Dynamic global vegetation models (DGVMs) are powerful tools to project past, current and future vegetation patterns and associated biogeochemical cycles. However, most models are limited by how they define vegetation and by their simplistic representation of competition. We discuss how concepts from community assembly theory and coexistence theory can help to improve vegetation models. We further present a trait- and individual-based vegetation model (aDGVM2) that allows individual plants to adopt a unique combination of trait values. These traits define how individual plants grow and compete. A genetic optimization algorithm is used to simulate trait inheritance and reproductive isolation between individuals. These model properties allow the assembly of plant communities that are adapted to a site's biotic and abiotic conditions. The aDGVM2 simulates how environmental conditions influence the trait spectra of plant communities; that fire selects for traits that enhance fire protection and reduces trait diversity; and the emergence of life-history strategies that are suggestive of colonization-competition trade-offs. The aDGVM2 deals with functional diversity and competition fundamentally differently from current DGVMs. This approach may yield novel insights as to how vegetation may respond to climate change and we believe it could foster collaborations between functional plant biologists and vegetation modellers.

Over 140 Mha of restoration commitments have been pledged across the global tropics, yet guidance is needed to identify those landscapes where implementation is likely to provide the greatest potential benefits and cost-effective outcomes. By overlaying seven recent, peer-reviewed spatial datasets as proxies for socioenvironmental benefits and feasibility of restoration, we identified restoration opportunities (areas with higher potential return of benefits and feasibility) in lowland tropical rainforest landscapes. We found restoration opportunities throughout the tropics. Areas scoring in the top 10% (i.e., restoration hotspots) are located largely within conservation hotspots (88%) and in countries committed to the Bonn Challenge (73%), a global effort to restore 350 Mha by 2030. However, restoration hotspots represented only a small portion (19.1%) of the Key Biodiversity Area network. Concentrating restoration investments in landscapes with high benefits and feasibility would maximize the potential to mitigate anthropogenic impacts and improve human well-being.

The high species richness of tropical forests has long been recognized, yet there remains substantial uncertainty regarding the actual number of tropical tree species. Using a pantropical tree inventory database from closed canopy forests, consisting of 657,630 trees belonging to 11,371 species, we use a fitted value of Fisher's alpha and an approximate pantropical stem total to estimate the minimum number of tropical forest tree species to fall between ∼ 40,000 and ∼ 53,000, i.e., at the high end of previous estimates. Contrary to common assumption, the Indo-Pacific region was found to be as species-rich as the Neotropics, with both regions having a minimum of ∼ 19,000-25,000 tree species. Continental Africa is relatively depauperate with a minimum of ∼ 4,500-6,000 tree species. Very few species are shared among the African, American, and the Indo-Pacific regions. We provide a methodological framework for estimating species richness in trees that may help refine species richness estimates of tree-dependent taxa.

Abstract Question Are plant traits more closely correlated with mean annual temperature, or with mean annual precipitation? Location Global. Methods We quantified the strength of the relationships between temperature and precipitation and 21 plant traits from 447,961 species‐site combinations worldwide. We used meta‐analysis to provide an overall answer to our question. Results Mean annual temperature was significantly more strongly correlated with plant traits than was mean annual precipitation. Conclusions Our study provides support for some of the assumptions of classical vegetation theory, and points to many interesting directions for future research. The relatively low R 2 values for precipitation might reflect the weak link between mean annual precipitation and the availability of water to plants.

We analysed the responses of 11 ecosystem models to elevated atmospheric [CO2 ] (eCO2 ) at two temperate forest ecosystems (Duke and Oak Ridge National Laboratory (ORNL) Free-Air CO2 Enrichment (FACE) experiments) to test alternative representations of carbon (C)-nitrogen (N) cycle processes. We decomposed the model responses into component processes affecting the response to eCO2 and confronted these with observations from the FACE experiments. Most of the models reproduced the observed initial enhancement of net primary production (NPP) at both sites, but none was able to simulate both the sustained 10-yr enhancement at Duke and the declining response at ORNL: models generally showed signs of progressive N limitation as a result of lower than observed plant N uptake. Nonetheless, many models showed qualitative agreement with observed component processes. The results suggest that improved representation of above-ground-below-ground interactions and better constraints on plant stoichiometry are important for a predictive understanding of eCO2 effects. Improved accuracy of soil organic matter inventories is pivotal to reduce uncertainty in the observed C-N budgets. The two FACE experiments are insufficient to fully constrain terrestrial responses to eCO2 , given the complexity of factors leading to the observed diverging trends, and the consequential inability of the models to explain these trends. Nevertheless, the ecosystem models were able to capture important features of the experiments, lending some support to their projections.

Phytophthora infestans, the cause of potato late blight, is infamous for having triggered the Irish Great Famine in the 1840s. Until the late 1970s, P. infestans diversity outside of its Mexican center of origin was low, and one scenario held that a single strain, US-1, had dominated the global population for 150 years; this was later challenged based on DNA analysis of historical herbarium specimens. We have compared the genomes of 11 herbarium and 15 modern strains. We conclude that the 19th century epidemic was caused by a unique genotype, HERB-1, that persisted for over 50 years. HERB-1 is distinct from all examined modern strains, but it is a close relative of US-1, which replaced it outside of Mexico in the 20th century. We propose that HERB-1 and US-1 emerged from a metapopulation that was established in the early 1800s outside of the species' center of diversity. DOI:http://dx.doi.org/10.7554/eLife.00731.001.