Optimization of Mycobacterium tuberculosis DNA processing prior to whole genome sequencing.

The process of whole genome sequencing of the Mycobacterium tuberculosis complex is dependent on complete the inactivation of the strain and subsequent DNA extraction. The objective of this study was to optimise the two steps. Firstly, the efficacy of Triton X-100 as a solvent for the inactivation step was evaluated. This solvent has been demonstrated to be effective in killing bacteria and is less toxic than the previously employed chloroform. For the extraction step, two lysis methods were evaluated: enzymatic (B1 protocol) and mechanical (B2 protocol). For whole genome sequencing, the Nextera XT DNA library preparation protocol was performed for both the B1 and B2 protocols. Subsequently, each library was subjected to whole-genome sequencing. The results demonstrated that heat lysis inactivation with Triton was effective, with no bacteria remaining viable following this treatment. The enzymatic and mechanical extraction protocols yielded comparable results in terms of DNA quantity and quality. The sequencing results showed that there was no significant difference in read depths between the two protocols. In conclusion, for MTBC strains, we recommend the use of our Triton inactivation method, which meets biosafety expectations.

Implementation of the WHO Tricycle protocol for surveillance of extended-spectrum ß-lactamase producing Escherichia coli in humans, chickens, and the environment in Madagascar: a prospective genomic epidemiology study.

BACKGROUND: Antimicrobial resistance (AMR) is a major public health threat, affecting not only people but also animals and the environment. The One Health dimension of AMR is well known; however, data are lacking on the circulation of resistance-conferring genes, particularly in low-income countries. In 2017, WHO proposed a protocol called Tricycle, focusing on extended-spectrum ß-lactamase (ESBL)-Escherichia coli surveillance in the three sectors (humans, animals, and the environment). We implemented Tricycle in Madagascar to assess ESBL-E coli prevalence and describe intrasector and intersector circulation of ESBL-E coli and plasmids. METHODS: In this prospective study, we collected blood culture data from hospitalised patients with a suspected bloodstream infection processed from May 1, 2018, to April 30, 2019, and rectal swabs from healthy pregnant women from July 30, 2018, to April 27, 2019, both from three hospitals in Antananarivo, Madagascar; and caeca from farm chickens and surface waters from the Ikopa river, wastewater, and slaughterhouse effluents in the Antananarivo area, Madagascar, from April 9, 2018, to April 30, 2019. All samples were tested for ESBL-E coli. The genomes of all isolates were sequenced using a short-read method on NextSeq 500 and NovaSeq 6000 platforms (Illumina, San Diego, CA, USA) and those carrying plasmid replicons using an additional long-read method on a MinION platform (Oxford Nanopore Technologies, Oxford, UK). We characterised genomes of isolated strains (sequence type, resistance and virulence gene content, and plasmid replicons). We then compared isolates using the variant calling method (single-nucleotide polymorphism). FINDINGS: Data from 1056 blood cultures were collected and 289 pregnant women, 246 chickens, and 28 surface waters were sampled. Of the blood cultures, 18 contained E coli, of which seven (39%) were ESBL. ESBL-E coli was present in samples from 86 (30%) of 289 pregnant women, 140 (57%) of 246 chickens, and 28 (100%) of 28 surface water samples. The wet season (November to April) was associated with higher rates of carriage in humans (odds ratio 3·08 [1·81-5·27]) and chickens (2·79 [1·65-4·81]). Sequencing of 277 non-duplicated isolates (82 from pregnant women, 118 from chickens, and 77 from environmental samples) showed high genetic diversity (90 sequence types identified) with sector-specific genomic features. Single nucleotide polymorphism (SNP) analysis revealed that 169 (61%) of 277 isolates grouped into 44 clusters (two or more isolates) of closely related isolates (<40 SNPs), of which 24 clusters contained isolates from two sectors and five contained isolates from all three sectors. ESBL genes were all blaCTX-M variants (215 [78%] of 277 being blaCTX-M-15) and were located on a plasmid in 113 (41%) of 277 isolates. These ESBL-carrying plasmids were mainly IncF (63 [55%] of 114; one strain carried two plasmids) and IncY (42 [37%] of 114). The F31/36:A4:B1 (n=13) and F-:A-:B53 (n=8) pMLST subtypes, and the IncY plasmids, which were all highly conserved, were observed in isolates of differing genetic backgrounds from all sectors and were transferable in vitro by conjugation. INTERPRETATION: Despite sector-specific population structures, both ESBL-E coli strains and plasmids are circulating among humans, chickens, and the environment in Antananarivo, Madagascar. The Tricycle protocol can be implemented in a low-income country and represents a powerful tool for investigating dissemination of AMR from a One Health perspective. FUNDING: Fondation Mérieux and INSERM, Université Paris Cité.

Towards the reconstruction of a global TB history using a new pipeline "TB-Annotator".

Mycobacterium tuberculosis complex (MTBC) has a population structure consisting of 9 human and animal lineages. The genomic diversity within these lineages is a pathogenesis factor that affects virulence, transmissibility, host response, and antibiotic resistance. Hence it is important to develop improved information systems for tracking and understanding the spreading and evolution of genomes. We present results obtained thanks to a new informatics platform for computational biology of MTBC, that uses a convenience sample from public/private SRAs, designated as TB-Annotator. Version 1 was a first interactive graphic-based web tool based on 15,901 representative genomes. Version 2, still interactive, is a more sophisticated database, developed using the Snakemake Workflow Management System (WMS) that allows an unsupervised global and scalable analysis of the content of the USA National Center for Biotechnology Information Short Read Archives database. This platform analyzes nucleotide variants, the presence/absence of genes, known regions of difference and detect new deletions, the insertion sites of mobile genetic elements, and allows phylogenetic trees to be built, imported in a graphical interface and interactively analyzed between the data and the tree. The objective of TB-Annotator is triple: detect recent epidemiological links, reconstruct distant phylogeographical histories as well as perform more complex phenotypic/genotypic Genome-Wide Association Studies (GWAS). In this paper, we compare the various taxonomic SNPs-based labels and hierarchies previously described in recent reference papers for L1, and present a comparative analysis that allows identification of alias and thus provides the basis of a future unifying naming scheme for L1 sublineages. We present a global phylogenetic tree built with RAxML-NG, and one on L2; at the time of writing, we characterized about 200 sublineages, with many new ones; a detail tree for Modern L2 and a hierarchical scheme allowing to facilitate L2 lineage assignment are also presented.

The bacterial genetic determinants of Escherichia coli capacity to cause bloodstream infections in humans.

Escherichia coli is both a highly prevalent commensal and a major opportunistic pathogen causing bloodstream infections (BSI). A systematic analysis characterizing the genomic determinants of extra-intestinal pathogenic vs. commensal isolates in human populations, which could inform mechanisms of pathogenesis, diagnostic, prevention and treatment is still lacking. We used a collection of 912 BSI and 370 commensal E. coli isolates collected in France over a 17-year period (2000-2017). We compared their pangenomes, genetic backgrounds (phylogroups, STs, O groups), presence of virulence-associated genes (VAGs) and antimicrobial resistance genes, finding significant differences in all comparisons between commensal and BSI isolates. A machine learning linear model trained on all the genetic variants derived from the pangenome and controlling for population structure reveals similar differences in VAGs, discovers new variants associated with pathogenicity (capacity to cause BSI), and accurately classifies BSI vs. commensal strains. Pathogenicity is a highly heritable trait, with up to 69% of the variance explained by bacterial genetic variants. Lastly, complementing our commensal collection with an older collection from 1980, we predict that pathogenicity continuously increased through 1980, 2000, to 2010. Together our findings imply that E. coli exhibit substantial genetic variation contributing to the transition between commensalism and pathogenicity and that this species evolved towards higher pathogenicity.

Epistatic interactions between the high pathogenicity island and other iron uptake systems shape Escherichia coli extra-intestinal virulence.

The intrinsic virulence of extra-intestinal pathogenic Escherichia coli is associated with numerous chromosomal and/or plasmid-borne genes, encoding diverse functions such as adhesins, toxins, and iron capture systems. However, the respective contribution to virulence of those genes seems to depend on the genetic background and is poorly understood. Here, we analyze genomes of 232 strains of sequence type complex STc58 and show that virulence (quantified in a mouse model of sepsis) emerged in a sub-group of STc58 due to the presence of the siderophore-encoding high-pathogenicity island (HPI). When extending our genome-wide association study to 370 Escherichia strains, we show that full virulence is associated with the presence of the aer or sit operons, in addition to the HPI. The prevalence of these operons, their co-occurrence and their genomic location depend on strain phylogeny. Thus, selection of lineage-dependent specific associations of virulence-associated genes argues for strong epistatic interactions shaping the emergence of virulence in E. coli.

Using gamification to improve engagement and learning outcomes in medical microbiology: the case study of 'BacteriaGame'.

The fight against antibiotic resistance has become a true global public health challenge of gargantuan proportions. Amongst the myriad of approaches being explored to tackle this predicament, one strategy involves enhancing prescriber knowledge and in particular their basic knowledge of medical bacteriology. Yet, as we well know in medical microbiology teachings, traditional lectures can be arduous, attempting to cram in a vast array of information in a limited time. An alternative solution to improve student engagement and enhance learning outcomes is to utilize educational games in complementary approach. Such games are an effective means of inspiring students to learn, encouraging self-assessment, and injecting diversity into the teaching process. To this end, we have developed and evaluated an educational card game, the 'BacteriaGame,' aimed at our medical students in medical bacteriology. Designed for students at the basic level, it serves as activity at the end of their apprenticeship to their bacteriology education. Additionally, it can also be used as a review tool by more advanced students, with teachers able to impart additional knowledge as the game progresses. We also use it in continuous training of medical laboratory staff. In this study, we evaluated the game at various stages of medical education, collecting feedback and analysing its impact on knowledge acquisition, comparing it to traditional lectures. Feedback from the majority of students revealed that the rules were clear, the game was enjoyable, and neither too lengthy nor too challenging. The integration of 'BacteriaGame' into their future training piqued their interest. In terms of learning outcomes, we discovered a significant increase in knowledge acquisition among those who used the game (P < .05). 'BacteriaGame' is now published by the French Society of Microbiology (SFM) and distributed in all medical and pharmacy schools thanks to a funding of the French Health Ministry. An English edition of the game is also available for international use as a physical copy to be purchased from the SFM. This will allow a large-scale distribution to colleagues who would like to use this game in their teaching.

Cross sectional survey of Varicella-Zoster virus and measles seropositivity in people living with HIV in a Parisian suburb and a review of current immunization guidelines.

According to evidence-based guidelines, vaccines against measles and varicella are generally recommended to susceptible HIV-positive patients, as long as they are not severely immunocompromised. However, routine screening to determine serologic status is not recommended. We conducted a seroprevalence study of anti-measles and anti-Varicella-Zoster virus (VZV) antibodies in adults living with HIV (PLWHA) consulting at Avicenne University Hospital in a Parisian suburb. Sera were collected in years 2018-2020 and tested by commercial immunoassays in 268 patients. Most of the patients were born in Sub-Saharan Africa (55 %) and only 23 % in Europe. Measles and varicella seropositivity were present respectively in 91.4 % and 96.2 % of patients. One patient in ten was seronegative to at least one of tested diseases. In the univariate analysis, only younger age (p = 0.027) was associated with a higher risk of measles seronegativity, while shorter time since arrival in France (p < 0.001) and shorter time since HIV discovery (p = 0.007) were associated with a higher risk of VZV seronegativity. In multivariate analysis no association was found. This study highlights the absence of specific risk factors for VZV and measles seronegativity in PLWHA and supports the importance of routine screening, in order to increase immunization rates and reduce risk of complications.

Performances of bioinformatics tools for the analysis of sequencing data of Mycobacterium tuberculosis complex strains.

Whole-genome sequencing of Mycobacterium tuberculosis complex (MTBC) strains is a rapidly growing tool to obtain results regarding the resistance and phylogeny of the strains. We evaluated the performances of two bioinformatics tools for the analysis of whole-genome sequences of MTBC strains. Two hundred and twenty-seven MTBC strains were isolated and whole-genome sequenced at the laboratory of Avicenne Hospital between 2015 and 2021. We investigated the resistance and susceptibility status of strains using two online tools, Mykrobe and PhyResSE. We compared the genotypic and phenotypic resistance results obtained by drug susceptibility testing. Unlike with the Mykrobe tool, sequencing quality data were obtained using PhyResSE: average coverage of 98% and average depth of 119X. We found a similar concordance between phenotypic and genotypic results when determining susceptibility to first-line anti-tuberculosis drugs (95%) with both tools. The sensitivity and specificity of each tool compared to the phenotypic method were respectively 72% [52-87] and 98% [96-99] for Mykrobe and 76% [57-90] and 97% [94-99] for PhyResSE. Mykrobe and PhyResSE were easy to use and efficient. These platforms are accessible to people not trained in bioinformatics and constitute a complementary approach to phenotypic methods for the study of MTBC strains.

Use of Whole-Genome Sequencing to Explore Mycobacterium tuberculosis Complex Circulating in a Hotspot Department in France.

The Seine-Saint-Denis is the French metropolitan department with the highest incidence of tuberculosis (TB). Our aim was to explore epidemiological and phylogenetic characteristics of TB strains in this hotspot department. We performed WGS on 227 strains of Mycobacterium tuberculosis complex isolated from patients at the Avicenne Hospital from 2016 to 2021 and randomly selected to represent the clinical diversity of French TB localization. Clinical and demographic data were recorded for each TB patient. The mean age of patients was 36 years old. They came from Africa (44%), Asia (27%), Europe (26%) and America (3%). Strains isolated from extrapulmonary samples were associated with Asian patients, whereas strains isolated from pulmonary samples were associated with European patients. We observed a high level of lineage diversity in line with the known worldwide diversity. Interestingly, lineage 3 was associated with lymph node TB. Additionally, the sensitivity of WGS for predicting resistance was 100% for rifampicin, isoniazid and ethambutol and 66.7% for pyrazinamide. The global concordance with drug-susceptibility testing using the phenotypic approach was 97%. In microbiology laboratories, WGS turns out to be an essential tool for better understanding local TB epidemiology, with direct access to circulating lineage identification and to drug susceptibilities to first- and second-line anti-TB drugs.

Cytokine Landscape in Central Nervous System Metastases.

The central nervous system is the location of metastases in more than 40% of patients with lung cancer, breast cancer and melanoma. These metastases are associated with one of the poorest prognoses in advanced cancer patients, mainly due to the lack of effective treatments. In this review, we explore the involvement of cytokines, including interleukins and chemokines, during the development of brain and leptomeningeal metastases from the epithelial-to-mesenchymal cell transition and blood-brain barrier extravasation to the interaction between cancer cells and cells from the brain microenvironment, including astrocytes and microglia. Furthermore, the role of the gut-brain axis on cytokine release during this process will also be addressed.

The Population Genomics of Increased Virulence and Antibiotic Resistance in Human Commensal Escherichia coli over 30 Years in France.

Escherichia coli is a commensal species of the lower intestine but is also a major pathogen causing intestinal and extraintestinal infections that is increasingly prevalent and resistant to antibiotics. Most studies on genomic evolution of E. coli used isolates from infections. Here, instead, we whole-genome sequenced a collection of 403 commensal E. coli isolates from fecal samples of healthy adult volunteers in France (1980 to 2010). These isolates were distributed mainly in phylogroups A and B2 (30% each) and belonged to 152 sequence types (STs), the five most frequent being ST10 (phylogroup A; 16.3%), ST73 and ST95 (phylogroup B2; 6.3 and 5.0%, respectively), ST69 (phylogroup D; 4.2%), and ST59 (phylogroup F; 3.9%), and 224 O:H serotypes. ST and serotype diversity increased over time. The O1, O2, O6, and O25 groups used in bioconjugate O-antigen vaccine against extraintestinal infections were found in 23% of the strains of our collection. The increase in frequency of virulence-associated genes and antibiotic resistance was driven by two evolutionary mechanisms. Evolution of virulence gene frequency was driven by both clonal expansion of STs with more virulence genes ("ST-driven") and increases in gene frequency within STs independent of changes in ST frequencies ("gene-driven"). In contrast, the evolution of resistance was dominated by increases in frequency within STs ("gene-driven"). This study provides a unique picture of the phylogenomic evolution of E. coli in its human commensal habitat over 30 years and will have implications for the development of preventive strategies. IMPORTANCE Escherichia coli is an opportunistic pathogen with the greatest burden of antibiotic resistance, one of the main causes of bacterial infections and an increasing concern in an aging population. Deciphering the evolutionary dynamics of virulence and antibiotic resistance in commensal E. coli is important to understand adaptation and anticipate future changes. The gut of vertebrates is the primary habitat of E. coli and probably where selection for virulence and resistance takes place. Unfortunately, most whole-genome-sequenced strains are isolated from pathogenic conditions. Here, we whole-genome sequenced 403 E. coli commensals isolated from healthy French subjects over a 30-year period. Virulence genes increased in frequency by both clonal expansion of clones carrying them and increases in frequency within clones, whereas resistance genes increased by within-clone increased frequency. Prospective studies of E. coli commensals should be performed worldwide to have a broader picture of evolution and adaptation of this species.

A global reptile assessment highlights shared conservation needs of tetrapods.

Comprehensive assessments of species' extinction risks have documented the extinction crisis1 and underpinned strategies for reducing those risks2. Global assessments reveal that, among tetrapods, 40.7% of amphibians, 25.4% of mammals and 13.6% of birds are threatened with extinction3. Because global assessments have been lacking, reptiles have been omitted from conservation-prioritization analyses that encompass other tetrapods4-7. Reptiles are unusually diverse in arid regions, suggesting that they may have different conservation needs6. Here we provide a comprehensive extinction-risk assessment of reptiles and show that at least 1,829 out of 10,196 species (21.1%) are threatened-confirming a previous extrapolation8 and representing 15.6 billion years of phylogenetic diversity. Reptiles are threatened by the same major factors that threaten other tetrapods-agriculture, logging, urban development and invasive species-although the threat posed by climate change remains uncertain. Reptiles inhabiting forests, where these threats are strongest, are more threatened than those in arid habitats, contrary to our prediction. Birds, mammals and amphibians are unexpectedly good surrogates for the conservation of reptiles, although threatened reptiles with the smallest ranges tend to be isolated from other threatened tetrapods. Although some reptiles-including most species of crocodiles and turtles-require urgent, targeted action to prevent extinctions, efforts to protect other tetrapods, such as habitat preservation and control of trade and invasive species, will probably also benefit many reptiles.

Overexpression of NOX2 Exacerbates AngII-Mediated Cardiac Dysfunction and Metabolic Remodelling.

The present study aimed to examine the effects of low doses of angiotensin II (AngII) on cardiac function, myocardial substrate utilization, energetics, and mitochondrial function in C57Bl/6J mice and in a transgenic mouse model with cardiomyocyte specific upregulation of NOX2 (csNOX2 TG). Mice were treated with saline (sham), 50 or 400 ng/kg/min of AngII (AngII50 and AngII400) for two weeks. In vivo blood pressure and cardiac function were measured using plethysmography and echocardiography, respectively. Ex vivo cardiac function, mechanical efficiency, and myocardial substrate utilization were assessed in isolated perfused working hearts, and mitochondrial function was measured in left ventricular homogenates. AngII50 caused reduced mechanical efficiency despite having no effect on cardiac hypertrophy, function, or substrate utilization. AngII400 slightly increased systemic blood pressure and induced cardiac hypertrophy with no effect on cardiac function, efficiency, or substrate utilization. In csNOX2 TG mice, AngII400 induced cardiac hypertrophy and in vivo cardiac dysfunction. This was associated with a switch towards increased myocardial glucose oxidation and impaired mitochondrial oxygen consumption rates. Low doses of AngII may transiently impair cardiac efficiency, preceding the development of hypertrophy induced at higher doses. NOX2 overexpression exacerbates the AngII -induced pathology, with cardiac dysfunction and myocardial metabolic remodelling.

The genomic view of diversification.

The process of species diversification is traditionally summarized by a single tree, the species tree, whose reconstruction from molecular data is hindered by frequent conflicts between gene genealogies. Here, we argue that instead of seeing these conflicts as nuisances, we can exploit them to inform the diversification process itself. We adopt a gene-based view of diversification to model the ubiquitous presence of gene flow between diverging lineages, one of the most important processes explaining disagreements among gene trees. We propose a new framework for modelling the joint evolution of gene and species lineages relaxing the hierarchy between the species tree and gene trees inherent to the standard view, as embodied in a popular model known as the multispecies coalescent (MSC). We implement this framework in two alternative models called the gene-based diversification models (GBD): (a) GBD-forward following all evolving genomes through time and (b) GBD-backward based on coalescent theory. They feature four parameters tuning colonization, gene flow, genetic drift and genetic differentiation. We propose an inference method based on differences between gene trees. Applied to two empirical data sets prone to gene flow, we find better support for the GBD-backward model than for the MSC model. Along with the increasing awareness of the extent of gene flow, this work shows the importance of considering the richer signal contained in genomic histories, rather than in the mere species tree, to better apprehend the complex evolutionary history of species.

Model-Based Inference of Punctuated Molecular Evolution.

In standard models of molecular evolution, DNA sequences evolve through asynchronous substitutions according to Poisson processes with a constant rate (called the molecular clock) or a rate that can vary (relaxed clock). However, DNA sequences can also undergo episodes of fast divergence that will appear as synchronous substitutions affecting several sites simultaneously at the macroevolutionary timescale. Here, we develop a model, which we call the Relaxed Clock with Spikes model, combining basal, clock-like molecular substitutions with episodes of fast divergence called spikes arising at speciation events. Given a multiple sequence alignment and its time-calibrated species phylogeny, our model is able to detect speciation events (including hidden ones) cooccurring with spike events and to estimate the probability and amplitude of these spikes on the phylogeny. We identify the conditions under which spikes can be distinguished from the natural variance of the clock-like component of molecular substitutions and from variations of the clock. We apply the method to genes underlying snake venom proteins and identify several spikes at gene-specific locations in the phylogeny. This work should pave the way for analyses relying on whole genomes to inform on modes of species diversification.

Species diversity as a surrogate for conservation of phylogenetic and functional diversity in terrestrial vertebrates across the Americas.

Preserving the evolutionary history and ecological functions that different species embody, in addition to species themselves, is a growing concern for conservation. Recent studies warn that conservation priority regions identified using species diversity differ from those based on phylogenetic or functional diversity. However, spatial mismatches in conservation priority regions need not indicate low surrogacy among these dimensions in conservation planning. Here, we use data for 10,213 terrestrial vertebrate species across the Americas to evaluate surrogacy; that is, the proportion of phylogenetic or functional diversity represented in conservation plans targeting species. We find that most conservation plans targeting species diversity also represent phylogenetic and functional diversity well, despite spatial mismatches in the priority regions identified by each plan. However, not all phylogenetic and functional diversity is represented within species-based plans, with the highest-surrogacy conservation strategy depending on the proportion of land area included in plans. Our results indicate that targeting species diversity could be sufficient to preserve much of the phylogenetic and functional dimensions of biodiversity in terrestrial vertebrates of the Americas. Incorporating phylogenetic and functional data in broad-scale conservation planning may not always be necessary, especially when the cost of doing so is high.

Molecular evidence for the paraphyly of Scolecophidia and its evolutionary implications.

The phylogenetic relationships between the three main clades of worm snakes remain controversial. This question is, however, crucial to elucidate the origin of the successful snake radiation, as these burrowing and miniaturized wormlike organisms represent the earliest branching clades within the snake tree. The present molecular phylogenetic study, intended to minimize the amount of missing data, provides fully resolved inter-subfamilial relationships among Typhlopidae. It also brings robust evidence that worm snakes (Scolecophidia) are paraphyletic, with the scolecophidian family Anomalepididae recovered with strong support as sister clade to the 'typical snakes' (Alethinophidia). Ancestral state reconstructions applied to three different traits strongly associated to a burrowing life-style (scolecoidy, absence of retinal cones and microstomy) provide results in favour of a burrowing origin of snakes, and suggest that worm snakes might be the only extant fossorial representatives of the primordial snake incursion towards an underground environment.

Undersampling Genomes has Biased Time and Rate Estimates Throughout the Tree of Life.

Genomic data drive evolutionary research on the relationships and timescale of life but the genomes of most species remain poorly sampled. Phylogenetic trees can be reconstructed reliably using small data sets and the same has been assumed for the estimation of divergence time with molecular clocks. However, we show here that undersampling of molecular data results in a bias expressed as disproportionately shorter branch lengths and underestimated divergence times in the youngest nodes and branches, termed the small sample artifact. In turn, this leads to increasing speciation and diversification rates towards the present. Any evolutionary analyses derived from these biased branch lengths and speciation rates will be similarly biased. The widely used timetrees of the major species-rich studies of amphibians, birds, mammals, and squamate reptiles are all data-poor and show upswings in diversification rate, suggesting that their results were biased by undersampling. Our results show that greater sampling of genomes is needed for accurate time and rate estimation, which are basic data used in ecological and evolutionary research.

Evolutionary time drives global tetrapod diversity.

Global variation in species richness is widely recognized, but the explanation for what drives it continues to be debated. Previous efforts have focused on a subset of potential drivers, including evolutionary rate, evolutionary time (maximum clade age of species restricted to a region), dispersal (migration from one region to another), ecological factors and climatic stability. However, no study has evaluated these competing hypotheses simultaneously at a broad spatial scale. Here, we examine their relative contribution in determining the richness of the most comprehensive dataset of tetrapods to our knowledge (84% of the described species), distinguishing between the direct influences of evolutionary rate, evolutionary time and dispersal, and the indirect influences of ecological factors and climatic stability through their effect on direct factors. We found that evolutionary time exerted a primary influence on species richness, with evolutionary rate being of secondary importance. By contrast, dispersal did not significantly affect richness patterns. Ecological and climatic stability factors influenced species richness indirectly by modifying evolutionary time (i.e. persistence time) and rate. Overall, our findings suggest that global heterogeneity in tetrapod richness is explained primarily by the length of time species have had to diversify.