Evolutionary insights into the emergence of virulent Leptospira spirochetes.

Pathogenic Leptospira are spirochete bacteria which cause leptospirosis, a re-emerging zoonotic disease of global importance. Here, we use a recently described lineage of environmental-adapted leptospires, which are evolutionarily the closest relatives of the highly virulent Leptospira species, to explore the key phenotypic traits and genetic determinants of Leptospira virulence. Through a comprehensive approach integrating phylogenomic comparisons with in vitro and in vivo phenotyping studies, we show that the evolution towards pathogenicity is associated with both a decrease of the ability to survive in the environment and the acquisition of strategies that enable successful host colonization. This includes the evasion of the mammalian complement system and the adaptations to avoid activation of the innate immune cells by the highly-virulent Leptospira species (also called P1+ species), unlike other species belonging to the phylogenetically related P1- and P2 groups, as well as saprophytes. Moreover, our analysis reveals specific genetic determinants that have undergone positive selection during the course of evolution in Leptospira, contributing directly to virulence and host adaptation as demonstrated by gain-of-function and knock-down studies. Taken together, our findings define a new vision on Leptospira pathogenicity, identifying virulence attributes associated with clinically relevant species, and provide insights into the evolution and emergence of these life-threatening pathogens.

Monkeypox Virus Infection in Humans across 16 Countries - April-June 2022.

BACKGROUND: Before April 2022, monkeypox virus infection in humans was seldom reported outside African regions where it is endemic. Currently, cases are occurring worldwide. Transmission, risk factors, clinical presentation, and outcomes of infection are poorly defined. METHODS: We formed an international collaborative group of clinicians who contributed to an international case series to describe the presentation, clinical course, and outcomes of polymerase-chain-reaction-confirmed monkeypox virus infections. RESULTS: We report 528 infections diagnosed between April 27 and June 24, 2022, at 43 sites in 16 countries. Overall, 98% of the persons with infection were gay or bisexual men, 75% were White, and 41% had human immunodeficiency virus infection; the median age was 38 years. Transmission was suspected to have occurred through sexual activity in 95% of the persons with infection. In this case series, 95% of the persons presented with a rash (with 64% having ≤10 lesions), 73% had anogenital lesions, and 41% had mucosal lesions (with 54 having a single genital lesion). Common systemic features preceding the rash included fever (62%), lethargy (41%), myalgia (31%), and headache (27%); lymphadenopathy was also common (reported in 56%). Concomitant sexually transmitted infections were reported in 109 of 377 persons (29%) who were tested. Among the 23 persons with a clear exposure history, the median incubation period was 7 days (range, 3 to 20). Monkeypox virus DNA was detected in 29 of the 32 persons in whom seminal fluid was analyzed. Antiviral treatment was given to 5% of the persons overall, and 70 (13%) were hospitalized; the reasons for hospitalization were pain management, mostly for severe anorectal pain (21 persons); soft-tissue superinfection (18); pharyngitis limiting oral intake (5); eye lesions (2); acute kidney injury (2); myocarditis (2); and infection-control purposes (13). No deaths were reported. CONCLUSIONS: In this case series, monkeypox manifested with a variety of dermatologic and systemic clinical findings. The simultaneous identification of cases outside areas where monkeypox has traditionally been endemic highlights the need for rapid identification and diagnosis of cases to contain further community spread.

The microbiological outcomes of culture-negative blood specimens using 16S rRNA broad-range PCR sequencing: a retrospective study in a Canadian province from 2018 to 2022.

Broad-range 16S rRNA PCR and sequencing of 1,183 blood specimens from 853 unique patients yielded an interpretable sequence and bacterial identification in 29%, 16S rRNA amplification with uninterpretable sequences in 53%, and no amplification in 18%. This study highlights the potential utility of this technique in identifying fastidious gram-negative and anaerobic bacteria but the frequent recovery of environmental and contaminant organisms argues for its judicious use. IMPORTANCE: The existing literature focuses on its performance compared to blood cultures in patients with sepsis, leaving a gap in the literature regarding other blood specimens in suspected infectious syndrome across the severity spectrum. We aimed to characterize its microbiological outcomes and provide insight into its potential clinical utility.

Monkeypox in Montréal: Epidemiology, Phylogenomics, and Public Health Response to a Large North American Outbreak.

BACKGROUND: Monkeypox, a viral zoonotic disease, is causing a global outbreak outside of endemic areas. OBJECTIVE: To characterize the outbreak of monkeypox in Montréal, the first large outbreak in North America. DESIGN: Epidemiologic and laboratory surveillance data and a phylogenomic analysis were used to describe and place the outbreak in a global context. SETTING: Montréal, Canada. PATIENTS: Probable or confirmed cases of monkeypox. MEASUREMENTS: Epidemiologic, clinical, and demographic data were aggregated. Whole-genome sequencing and phylogenetic analysis were performed for a set of outbreak sequences. The public health response and its evolution are described. RESULTS: Up to 18 October 2022, a total of 402 cases of monkeypox were reported mostly among men who have sex with men (MSM), most of which were suspected to be acquired through sexual contact. All monkeypox genomes nested within the B.1 lineage. Montréal Public Health worked closely with the affected communities to control the outbreak, becoming the first jurisdiction to offer 1 dose of the Modified Vaccinia Ankara-Bavarian Nordic vaccine as preexposure prophylaxis (PrEP) to those at risk in early June 2022. Two peaks of cases were seen in early June and July (43 and 44 cases per week, respectively) followed by a decline toward near resolution of the outbreak in October. Reasons for the biphasic peak are not fully elucidated but may represent the tempo of vaccination and/or several factors related to transmission dynamics and case ascertainment. LIMITATIONS: Clinical data are self-reported. Limited divergence among sequences limited genomic epidemiologic conclusions. CONCLUSION: A large outbreak of monkeypox occurred in Montréal, primarily among MSM. Successful control of the outbreak rested on early and sustained engagement with the affected communities and rapid offer of PrEP vaccination to at-risk persons. PRIMARY FUNDING SOURCE: None.

Laboratory Features of Trichinellosis and Eosinophilia Threshold for Testing, Nunavik, Quebec, Canada, 2009-2019.

Prolonged eosinophilia is characteristic of trichinellosis. To determine the optimal eosinophil threshold for reflex Trichinella testing, we examined all 43 cases in Nunavik, Quebec, Canada, during 2009-2019. Using receiver operating characteristic analysis, we determined that eosinophil counts >0.8 × 109 cells/L should prompt consideration of trichinellosis and testing to rapidly identify potential outbreaks.

Specificity of SARS-CoV-2 Antibody Detection Assays against S and N Proteins among Pre-COVID-19 Sera from Patients with Protozoan and Helminth Parasitic Infections.

We aimed to assess the specificity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) antibody detection assays among people with tissue-borne parasitic infections. We tested three SARS-CoV-2 antibody-detection assays (cPass SARS-CoV-2 neutralization antibody detection kit [cPass], Abbott SARS-CoV-2 IgG assay [Abbott Architect], and Standard Q COVID-19 IgM/IgG combo rapid diagnostic test [SD RDT IgM/SD RDT IgG]) among 559 pre-COVID-19 seropositive sera for several parasitic infections. The specificity of assays was 95 to 98% overall. However, lower specificity was observed among sera from patients with protozoan infections of the reticuloendothelial system, such as human African trypanosomiasis (Abbott Architect; 88% [95% CI, 75 to 95]) and visceral leishmaniasis (SD RDT IgG; 80% [95% CI, 30 to 99]), and from patients with recent malaria in areas of Senegal where malaria is holoendemic (ranging from 91% for Abbott Architect and SD RDT IgM to 98 to 99% for cPass and SD RDT IgG). For specimens from patients with evidence of past or present helminth infection overall, test specificity estimates were all ≥96%. Sera collected from patients clinically suspected of parasitic infections that tested negative for these infections yielded a specificity of 98 to 100%. The majority (>85%) of false-positive results were positive by only one assay. The specificity of SARS-CoV-2 serological assays among sera from patients with tissue-borne parasitic infections was below the threshold required for decisions about individual patient care. Specificity is markedly increased by the use of confirmatory testing with a second assay. Finally, the SD RDT IgG proved similarly specific to laboratory-based assays and provides an option in low-resource settings when detection of anti-SARS-CoV-2 IgG is indicated.

Remdesivir for the treatment of patients in hospital with COVID-19 in Canada: a randomized controlled trial.

BACKGROUND: The role of remdesivir in the treatment of patients in hospital with COVID-19 remains ill defined in a global context. The World Health Organization Solidarity randomized controlled trial (RCT) evaluated remdesivir in patients across many countries, with Canada enrolling patients using an expanded data collection format in the Canadian Treatments for COVID-19 (CATCO) trial. We report on the Canadian findings, with additional demographics, characteristics and clinical outcomes, to explore the potential for differential effects across different health care systems. METHODS: We performed an open-label, pragmatic RCT in Canadian hospitals, in conjunction with the Solidarity trial. We randomized patients to 10 days of remdesivir (200 mg intravenously [IV] on day 0, followed by 100 mg IV daily), plus standard care, or standard care alone. The primary outcome was in-hospital mortality. Secondary outcomes included changes in clinical severity, oxygen- and ventilator-free days (at 28 d), incidence of new oxygen or mechanical ventilation use, duration of hospital stay, and adverse event rates. We performed a priori subgroup analyses according to duration of symptoms before enrolment, age, sex and severity of symptoms on presentation. RESULTS: Across 52 Canadian hospitals, we randomized 1282 patients between Aug. 14, 2020, and Apr. 1, 2021, to remdesivir (n = 634) or standard of care (n = 648). Of these, 15 withdrew consent or were still in hospital, for a total sample of 1267 patients. Among patients assigned to receive remdesivir, in-hospital mortality was 18.7%, compared with 22.6% in the standard-of-care arm (relative risk [RR] 0.83 (95% confidence interval [CI] 0.67 to 1.03), and 60-day mortality was 24.8% and 28.2%, respectively (95% CI 0.72 to 1.07). For patients not mechanically ventilated at baseline, the need for mechanical ventilation was 8.0% in those assigned remdesivir, and 15.0% in those receiving standard of care (RR 0.53, 95% CI 0.38 to 0.75). Mean oxygen-free and ventilator-free days at day 28 were 15.9 (± standard deviation [SD] 10.5) and 21.4 (± SD 11.3) in those receiving remdesivir and 14.2 (± SD 11) and 19.5 (± SD 12.3) in those receiving standard of care (p = 0.006 and 0.007, respectively). There was no difference in safety events of new dialysis, change in creatinine, or new hepatic dysfunction between the 2 groups. INTERPRETATION: Remdesivir, when compared with standard of care, has a modest but significant effect on outcomes important to patients and health systems, such as the need for mechanical ventilation. Trial registration: ClinicalTrials.gov, no. NCT04330690.

Among-character rate variation distributions in phylogenetic analysis of discrete morphological characters.

Likelihood-based methods are commonplace in phylogenetic systematics. Although much effort has been directed toward likelihood-based models for molecular data, comparatively less work has addressed models for discrete morphological character (DMC) data. Among-character rate variation (ACRV) may confound phylogenetic analysis, but there have been few analyses of the magnitude and distribution of rate heterogeneity among DMCs. Using 76 data sets covering a range of plants, invertebrate, and vertebrate animals, we used a modified version of MrBayes to test equal, gamma-distributed and lognormally distributed models of ACRV, integrating across phylogenetic uncertainty using Bayesian model selection. We found that in approximately 80% of data sets, unequal-rates models outperformed equal-rates models, especially among larger data sets. Moreover, although most data sets were equivocal, more data sets favored the lognormal rate distribution relative to the gamma rate distribution, lending some support for more complex character correlations than in molecular data. Parsimony estimation of the underlying rate distributions in several data sets suggests that the lognormal distribution is preferred when there are many slowly evolving characters and fewer quickly evolving characters. The commonly adopted four rate category discrete approximation used for molecular data was found to be sufficient to approximate a gamma rate distribution with discrete characters. However, among the two data sets tested that favored a lognormal rate distribution, the continuous distribution was better approximated with at least eight discrete rate categories. Although the effect of rate model on the estimation of topology was difficult to assess across all data sets, it appeared relatively minor between the unequal-rates models for the one data set examined carefully. As in molecular analyses, we argue that researchers should test and adopt the most appropriate model of rate variation for the data set in question. As discrete characters are increasingly used in more sophisticated likelihood-based phylogenetic analyses, it is important that these studies be built on the most appropriate and carefully selected underlying models of evolution.

Copper management strategies in obligate bacterial symbionts: balancing cost and benefit.

Bacteria employ diverse mechanisms to manage toxic copper in their environments, and these evolutionary strategies can be divided into two main categories: accumulation and rationalization of metabolic pathways. The strategies employed depend on the bacteria's lifestyle and environmental context, optimizing the metabolic cost-benefit ratio. Environmental and opportunistically pathogenic bacteria often possess an extensive range of copper regulation systems in order to respond to variations in copper concentrations and environmental conditions, investing in diversity and/or redundancy as a safeguard against uncertainty. In contrast, obligate symbiotic bacteria, such as Neisseria gonorrhoeae and Bordetella pertussis, tend to have specialized and more parsimonious copper regulation systems designed to function in the relatively stable host environment. These evolutionary strategies maintain copper homeostasis even in challenging conditions like encounters within phagocytic cells. These examples highlight the adaptability of bacterial copper management systems, tailored to their specific lifestyles and environmental requirements, in the context of an evolutionary the trade-off between benefits and energy costs.

An imputed ancestral reference genome for the Mycobacterium tuberculosis complex better captures structural genomic diversity for reference-based alignment workflows.

Reference-based alignment of short-reads is a widely used technique in genomic analysis of the Mycobacterium tuberculosis complex (MTBC) and the choice of reference sequence impacts the interpretation of analyses. The most widely used reference genomes include the ATCC type strain (H37Rv) and the putative MTBC ancestral sequence of Comas et al. both of which are based on a lineage 4 sequence. As such, these reference sequences do not capture all of the structural variation known to be present in the ancestor of the MTBC. To better represent the base of the MTBC, we generated an imputed ancestral genomic sequence, termed MTBC0 from reference-free alignments of closed MTBC genomes. When used as a reference sequence in alignment workflows, MTBC0 mapped more short sequencing reads and called more pairwise SNPs relative to the Comas et al. sequence while exhibiting minimal impact on the overall phylogeny of MTBC. The results also show that MTBC0 provides greater fidelity in capturing genomic variation and allows for the inclusion of regions absent from H37Rv in standard MTBC workflows without additional steps. The use of MTBC0 as an ancestral reference sequence in standard workflows modestly improved read mapping, SNP calling and intuitively facilitates the study of structural variation and evolution in MTBC.

Evolutionary insights into the emergence of virulent Leptospira spirochetes.

Pathogenic Leptospira are spirochete bacteria which cause leptospirosis, a re-emerging zoonotic disease of global importance. Here, we use a recently described lineage of environmental-adapted leptospires, which are evolutionarily the closest relatives of the highly virulent Leptospira species, to explore the key phenotypic traits and genetic determinants of Leptospira virulence. Through a comprehensive approach integrating phylogenomic comparisons with in vitro and in vivo phenotyping studies, we show that the evolution towards pathogenicity is associated with both a decrease of the ability to survive in the environment and the acquisition of strategies that enable successful host colonization. This includes the evasion of the human complement system and the adaptations to avoid activation of the innate immune cells. Moreover, our analysis reveals specific genetic determinants that have undergone positive selection during the course of evolution in Leptospira, contributing directly to virulence and host adaptation as demonstrated by gain-of-function and knock-down studies. Taken together, our findings define a new vision on Leptospira pathogenicity, identifying virulence attributes associated with clinically relevant species, and provide insights into the evolution and emergence of these life-threatening pathogens.

Ossification sequence heterochrony among amphibians.

Heterochrony is an important mechanism in the evolution of amphibians. Although studies have centered on the relationship between size and shape and the rates of development, ossification sequence heterochrony also may have been important. Rigorous, phylogenetic methods for assessing sequence heterochrony are relatively new, and a comprehensive study of the relative timing of ossification of skeletal elements has not been used to identify instances of sequence heterochrony across Amphibia. In this study, a new version of the program Parsimov-based genetic inference (PGi) was used to identify shifts in ossification sequences across all extant orders of amphibians, for all major structural units of the skeleton. PGi identified a number of heterochronic sequence shifts in all analyses, the most interesting of which seem to be tied to differences in metamorphic patterns among major clades. Early ossification of the vomer, premaxilla, and dentary is retained by Apateon caducus and members of Gymnophiona and Urodela, which lack the strongly biphasic development seen in anurans. In contrast, bones associated with the jaws and face were identified as shifting late in the ancestor of Anura. The bones that do not shift late, and thereby occupy the earliest positions in the anuran cranial sequence, are those in regions of the skull that undergo the least restructuring throughout anuran metamorphosis. Additionally, within Anura, bones of the hind limb and pelvic girdle were also identified as shifting early in the sequence of ossification, which may be a result of functional constraints imposed by the drastic metamorphosis of most anurans.

In-hospital outcomes of SARS-CoV-2-infected health care workers in the COVID-19 pandemic first wave, Quebec, Canada.

BACKGROUND: Health care workers (HCW), particularly immigrants and ethnic minorities are at increased risk for SARS-CoV-2 infection. Outcomes during a COVID-19 associated hospitalization are not well described among HCW. We aimed to describe the characteristics of HCW admitted with COVID-19 including immigrant status and ethnicity and the associated risk factors for Intensive Care unit (ICU) admission and death. METHODS: Adults with laboratory-confirmed community-acquired COVID-19 hospitalized from March 1 to June 30, 2020, at four tertiary-care hospitals in Montréal, Canada were included. Demographics, comorbidities, occupation, immigration status, country of birth, ethnicity, workplace exposures, and hospital outcomes (ICU admission and death) were obtained through a chart review and phone survey. A Fine and Gray competing risk proportional hazards model was used to estimate the risk of ICU admission among HCW stratified by immigrant status and region of birth. RESULTS: Among 1104 included persons, 150 (14%) were HCW, with a phone survey participation rate of 68%. HCWs were younger (50 vs 64 years; p<0.001), more likely to be female (61% vs 41%; p<0.001), migrants (68% vs 55%; p<0.01), non-White (65% vs 41%; p<0.001) and healthier (mean Charlson Comorbidity Index of 0.3 vs 1.2; p<0.001) compared to non-HCW. They were as likely to be admitted to the ICU (28% vs 31%; p = 0.40) but were less likely to die (4% vs. 17%; p<0.001). Immigrant HCW accounted for 68% of all HCW cases and, compared to Canadian HCW, were more likely to be personal support workers (PSW) (54% vs. 33%, p<0.01), to be Black (58% vs 4%) and to work in a Residential Care Facility (RCF) (59% vs 33%; p = 0.05). Most HCW believed that they were exposed at work, 55% did not always have access to personal protective equipment (PPE) and 40% did not receive COVID-19-specific Infection Control (IPAC) training. CONCLUSION: Immigrant HCW were particularly exposed to COVID-19 infection in the first wave of the pandemic in Quebec. Despite being young and healthy, one third of all HCW required ICU admission, highlighting the importance of preventing workplace transmission through strong infection prevention and control measures, including high COVID-19 vaccination coverage.

Remdesivir and systemic corticosteroids for the treatment of COVID-19: A Bayesian re-analysis.

BACKGROUND: The global death toll from coronavirus disease 2019 (COVID-19) has exceeded 2 million, and treatments to decrease mortality are needed urgently. OBJECTIVES: To examine the probabilities of a clinically meaningful reduction in mortality for remdesivir and systemic corticosteroids. DESIGN, SETTING AND PARTICIPANTS: This was a probabilistic re-analysis of clinical trial data for corticosteroids and remdesivir in the treatment of hospitalized patients with COVID-19 using a Bayesian random effects meta-analytic approach. Studies were identified from existing meta-analyses performed by the World Health Organization. MAIN OUTCOMES AND MEASURES: Posterior probabilities of an absolute decrease in mortality compared with control patients, by subgroups based on oxygen requirements, were calculated for corticosteroids and remdesivir. Probabilities of ≥1%, ≥2% and ≥5% absolute decrease in mortality were quantified. RESULTS: For patients needing mechanical ventilation, the probability of ≥1% absolute decrease in mortality was 4% for remdesivir and 93% for corticosteroids. For patients needing supplemental oxygen without mechanical ventilation, the probability of ≥1% absolute decrease in mortality was 81% for remdesivir and 93% for dexamethasone. Finally, for patients who did not need oxygen support, the probability of ≥1% absolute decrease in mortality was 29% for remdesivir and 4% for dexamethasone. CONCLUSIONS AND RELEVANCE: Using a Bayesian analytic approach, remdesivir had low probability of achieving a clinically meaningful reduction in mortality, except for patients needing supplemental oxygen without mechanical ventilation. Corticosteroids were more promising for patients needing oxygen support, especially mechanical ventilation. While awaiting more definitive studies, this probabilistic interpretation of the evidence will help to guide treatment decisions for clinicians, as well as guideline and policy makers.

Methods for the analysis of developmental sequence data.

Heterochrony, evolutionary changes in developmental timing, can be studied either by examining changes in growth or changes in the sequence of developmental events. Developmental sequence data has the potential to address many important questions in the field of developmental evolution, but methodological challenges remain due to the biological and logical dependence of events in a ranked sequence. In the past 10 years, the study of sequence heterochrony has undergone a rebirth, with the creation of several new methods for the analysis of this type of data. These methods can be divided into two broad categories: phenetic comparisons between terminal taxa that strive to uncover integrations within the developmental sequences and putative shared sequence heterochronies, and phylogeny-based methods that derive ancestor-descendent sequence heterochronies and establish statements of sequence evolution. In this review, we will discuss the strengths and weaknesses of the methodologies that have been proposed to quantitatively examine developmental sequence data, and studies that have attempted to implement them in an evolutionary context.

Evolution of budding yeast prion-determinant sequences across diverse fungi.

Prions are transmissible self-replicating alternative states of proteins. Four prions ([PSI+], [URE3], [RNQ+] and [NU+]) can be inherited cytoplasmically in Saccharomyces cerevisiae laboratory strains. In the case of [PSI+], there is increasing evidence that prion formation may engender mechanisms to uncover hidden genetic variation. Here, we have analysed the evolution of the prion-determinant (PD) domains across 21 fungi, focusing on compositional biases, repeats and substitution rates. We find evidence for constraint on all four PD domains, but each domain has its own evolutionary dynamics. For [PSI+], the Q/N bias is maintained in fungal clades that diverged one billion years ago, with purifying selection observed within the Saccharomyces species. The degree of Q/N bias is correlated with the degree of local homology to prion-associated repeats, which occur rarely in other proteins (<1% of sequences for the proteomes studied). The evolutionary conservation of Q/N bias in Sup35p is unusual, with only eight other S. cerevisiae proteins showing similar, phylogenetically deep patterns of bias conservation. The [URE3] PD domain is unique to Hemiascomycota; part of the PD domain shows purifying selection, whereas another part engenders bias changes between clades. Also, like for Sup35p, the [RNQ+] and [NU+] PD domains show purifying selection in Saccharomyces species. Additionally, in each proteome, we observe on average several hundred yeast-prion-like domains, with fewest in fission yeast. Our findings on yeast prion evolution provide further support for the functional significance of these molecules.

Wildlife population trends in protected areas predicted by national socio-economic metrics and body size.

Ensuring that protected areas (PAs) maintain the biodiversity within their boundaries is fundamental in achieving global conservation goals. Despite this objective, wildlife abundance changes in PAs are patchily documented and poorly understood. Here, we use linear mixed effect models to explore correlates of population change in 1,902 populations of birds and mammals from 447 PAs globally. On an average, we find PAs are maintaining populations of monitored birds and mammals within their boundaries. Wildlife population trends are more positive in PAs located in countries with higher development scores, and for larger-bodied species. These results suggest that active management can consistently overcome disadvantages of lower reproductive rates and more severe threats experienced by larger species of birds and mammals. The link between wildlife trends and national development shows that the social and economic conditions supporting PAs are critical for the successful maintenance of their wildlife populations.

Assessing the phylogenetic utility of sequence heterochrony: evolution of avian ossification sequences as a case study.

The evolution of developmental sequences, or sequence heterochrony, is an emerging field of study that addresses the temporal interplay between evolution and development. Some phylogenetic signal has been found in developmental sequence data, but sampling has generally been limited to small numbers of taxa and few developmental events. Here we present the largest ossification sequence dataset to date. The sequences are composed of ossification events throughout the avian skeleton, and are used to address the evolutionary signal of ossification sequence data within this clade. The results indicate that ossification sequences are conserved in birds, and show a stronger phylogenetic signal than previous studies, perhaps due to the volume of data. Phylogenetic signal is not strong enough, however, to consider ossification sequence data to be any better at resolving phylogenetic hypotheses than other morphological data and just as prone to evolutionary convergence. There is no one-to-one correlation between ossification sequence and developmental stage. We discuss some methodological implications of our findings, as well as commonalities in avian ossification sequences such as early ossification of the long bones relative to the dermatocranium, and of the hindlimb over the forelimb.