Large socioeconomic gap in period life expectancy and life years spent with complications of diabetes in the Scottish population with type 1 diabetes, 2013-2018.

BACKGROUND: We report the first study to estimate the socioeconomic gap in period life expectancy (LE) and life years spent with and without complications in a national cohort of individuals with type 1 diabetes. METHODS: This retrospective cohort study used linked healthcare records from SCI-Diabetes, the population-based diabetes register of Scotland. We studied all individuals aged 50 and older with a diagnosis of type 1 diabetes who were alive and residing in Scotland on 1 January 2013 (N = 8591). We used the Scottish Index of Multiple Deprivation (SIMD) 2016 as an area-based measure of socioeconomic deprivation. For each individual, we constructed a history of transitions by capturing whether individuals developed retinopathy/maculopathy, cardiovascular disease, chronic kidney disease, and diabetic foot, or died throughout the study period, which lasted until 31 December 2018. Using parametric multistate survival models, we estimated total and state-specific LE at an attained age of 50. RESULTS: At age 50, remaining LE was 22.2 years (95% confidence interval (95% CI): 21.6 - 22.8) for males and 25.1 years (95% CI: 24.4 - 25.9) for females. Remaining LE at age 50 was around 8 years lower among the most deprived SIMD quintile when compared with the least deprived SIMD quintile: 18.7 years (95% CI: 17.5 - 19.9) vs. 26.3 years (95% CI: 24.5 - 28.1) among males, and 21.2 years (95% CI: 19.7 - 22.7) vs. 29.3 years (95% CI: 27.5 - 31.1) among females. The gap in life years spent without complications was around 5 years between the most and the least deprived SIMD quintile: 4.9 years (95% CI: 3.6 - 6.1) vs. 9.3 years (95% CI: 7.5 - 11.1) among males, and 5.3 years (95% CI: 3.7 - 6.9) vs. 10.3 years (95% CI: 8.3 - 12.3) among females. SIMD differences in transition rates decreased marginally when controlling for time-updated information on risk factors such as HbA1c, blood pressure, BMI, or smoking. CONCLUSIONS: In addition to societal interventions, tailored support to reduce the impact of diabetes is needed for individuals from low socioeconomic backgrounds, including access to innovations in management of diabetes and the prevention of complications.

Flash monitor initiation is associated with improvements in HbA1c levels and DKA rates among people with type 1 diabetes in Scotland: a retrospective nationwide observational study.

AIMS/HYPOTHESIS: We assessed the real-world effect of flash monitor (FM) usage on HbA1c levels and diabetic ketoacidosis (DKA) and severe hospitalised hypoglycaemia (SHH) rates among people with type 1 diabetes in Scotland and across sociodemographic strata within this population. METHODS: This study was retrospective, observational and registry based. Using the national diabetes registry, 14,682 individuals using an FM at any point between 2014 and mid-2020 were identified. Within-person change from baseline in HbA1c following FM initiation was modelled using linear mixed models accounting for within-person pre-exposure trajectory. DKA and SHH events were captured through linkage to hospital admission and mortality data. The difference in DKA and SHH rates between FM-exposed and -unexposed person-time was assessed among users, using generalised linear mixed models with a Poisson likelihood. In a sensitivity analysis, we tested whether changes in these outcomes were seen in an age-, sex- and baseline HbA1c-matched sample of non-users over the same time period. RESULTS: Prevalence of ever-FM use was 45.9% by mid-2020, with large variations by age and socioeconomic status: 64.3% among children aged <13 years vs 32.7% among those aged ≥65 years; and 54.4% vs 36.2% in the least-deprived vs most-deprived quintile. Overall, the median (IQR) within-person change in HbA1c in the year following FM initiation was -2.5 (-9.0, 2.5) mmol/mol (-0.2 [-0.8, 0.2]%). The change varied widely by pre-usage HbA1c: -15.5 (-31.0, -4.0) mmol/mol (-1.4 [-2.8, -0.4]%) in those with HbA1c > 84 mmol/mol [9.8%] and 1.0 (-2.0, 5.5) mmol/mol (0.1 [-0.2, 0.5]%) in those with HbA1c < 54 mmol/mol (7.1%); the corresponding estimated fold change (95% CI) was 0.77 (0.76, 0.78) and 1.08 (1.07, 1.09). Significant reductions in HbA1c were found in all age bands, sexes and socioeconomic strata, and regardless of prior/current pump use, completion of a diabetes education programme or early FM adoption. Variation between the strata of these factors beyond that driven by differing HbA1c at baseline was slight. No change in HbA1c in matched non-users was observed in the same time period (median [IQR] within-person change = 0.5 [-5.0, 5.5] mmol/mol [0.0 (-0.5, 0.5)%]). DKA rates decreased after FM initiation overall and in all strata apart from the adolescents. Estimated overall reduction in DKA event rates (rate ratio) was 0.59 [95% credible interval (CrI) 0.53, 0.64]) after FM vs before FM initiation, accounting for pre-exposure trend. Finally, among those at higher risk for SHH, estimated reduction in event rates was rate ratio 0.25 (95%CrI 0.20, 0.32) after FM vs before FM initiation. CONCLUSIONS/INTERPRETATION: FM initiation is associated with clinically important reductions in HbA1c and striking reduction in DKA rate. Increasing uptake among the socioeconomically disadvantaged offers considerable potential for tightening the current socioeconomic disparities in glycaemia-related outcomes.

Rising Rates and Widening Socioeconomic Disparities in Diabetic Ketoacidosis in Type 1 Diabetes in Scotland: A Nationwide Retrospective Cohort Observational Study.

OBJECTIVE: Whether advances in the management of type 1 diabetes are reducing rates of diabetic ketoacidosis (DKA) is unclear. We investigated time trends in DKA rates in a national cohort of individuals with type 1 diabetes monitored for 14 years, overall and by socioeconomic characteristics. RESEARCH DESIGN AND METHODS: All individuals in Scotland with type 1 diabetes who were alive and at least 1 year old between 1 January 2004 and 31 December 2018 were identified using the national register (N = 37,939). DKA deaths and hospital admissions were obtained through linkage to Scottish national death and morbidity records. Bayesian regression was used to test for DKA time trends and association with risk markers, including socioeconomic deprivation. RESULTS: There were 30,427 DKA admissions and 472 DKA deaths observed over 393,223 person-years at risk. DKA event rates increased over the study period (incidence rate ratio [IRR] per year 1.058 [95% credibility interval 1.054-1.061]). Males had lower rates than females (IRR male-to-female 0.814 [0.776-0.855]). DKA incidence rose in all age-groups other than 10- to 19-year-olds, in whom rates were the highest, but fell over the study. There was a large socioeconomic differential (IRR least-to-most deprived quintile 0.446 [0.406-0.490]), which increased during follow-up. Insulin pump use or completion of structured education were associated with lower DKA rates, and antidepressant and methadone prescription were associated with higher DKA rates. CONCLUSIONS: DKA incidence has risen since 2004, except in 10- to 19-year-olds. Of particular concern are the strong and widening socioeconomic disparities in DKA outcomes. Efforts to prevent DKA, especially in vulnerable groups, require strengthening.

Divergence-time estimates for hominins provide insight into encephalization and body mass trends in human evolution.

Quantifying speciation times during human evolution is fundamental as it provides a timescale to test for the correlation between key evolutionary transitions and extrinsic factors such as climatic or environmental change. Here, we applied a total evidence dating approach to a hominin phylogeny to estimate divergence times under different topological hypotheses. The time-scaled phylogenies were subsequently used to perform ancestral state reconstructions of body mass and phylogenetic encephalization quotient (PEQ). Our divergence-time estimates are consistent with other recent studies that analysed extant species. We show that the origin of the genus Homo probably occurred between 4.30 and 2.56 million years ago. The ancestral state reconstructions show a general trend towards a smaller body mass before the emergence of Homo, followed by a trend towards a greater body mass. PEQ estimations display a general trend of gradual but accelerating encephalization evolution. The obtained results provide a rigorous temporal framework for human evolution.

Marked improvements in glycaemic outcomes following insulin pump therapy initiation in people with type 1 diabetes: a nationwide observational study in Scotland.

AIMS/HYPOTHESIS: Our aim was to assess the use of continuous subcutaneous insulin infusion (CSII) in people with type 1 diabetes in Scotland and its association with glycaemic control, as measured by HbA1c levels, frequency of diabetic ketoacidosis (DKA) and severe hospitalised hypoglycaemia (SHH), overall and stratified by baseline HbA1c. METHODS: We included 4684 individuals with type 1 diabetes from the national Scottish register, who commenced CSII between 2004 and 2019. We presented crude within-person differences from baseline HbA1c over time since initiation, crude DKA and SHH event-rates pre-/post-CSII exposure. We then used mixed models to assess the significance of CSII exposure, taking into account: (1) the diffuse nature of the intervention (i.e. structured education often precedes initiation); (2) repeated within-person measurements; and (3) background time-trends occurring pre-intervention. RESULTS: HbA1c decreased after CSII initiation, with a median within-person change of -5.5 mmol/mol (IQR -12.0, 0.0) (-0.5% [IQR -1.1, 0.0]). Within-person changes were most substantial in those with the highest baseline HbA1c, with median -21.0 mmol/mol (-30.0, -11.0) (-1.9% [-2.7, -1.0]) change in those with a baseline >84 mmol/mol (9.8%) within a year of exposure, that was sustained: -19.0 mmol/mol (-27.6, -6.5) (-1.7% [-2.5, -0.6]) at ≥5 years. Statistical significance and magnitude of change were supported by the mixed models results. The crude DKA event-rate was significantly lower in post-CSII person-time compared with pre-CSII person-time: 49.6 events (95% CI 46.3, 53.1) per 1000 person-years vs 67.9 (64.1, 71.9); rate ratio from Bayesian mixed models adjusting for pre-exposure trend: 0.61 (95% credible interval [CrI] 0.47, 0.77; posterior probability of reduction pp = 1.00). The crude overall SHH event-rate in post-CSII vs pre-CSII person-time was also lower: 17.8 events (95% CI 15.8, 19.9) per 1000 person-years post-exposure vs 25.8 (23.5, 28.3) pre-exposure; rate ratio from Bayesian mixed models adjusting for pre-exposure trend: 0.67 (95% CrI 0.45, 1.01; pp = 0.97). CONCLUSIONS/INTERPRETATION: CSII therapy was associated with marked falls in HbA1c especially in those with high baseline HbA1c. CSII was independently associated with reduced DKA and SHH rates. CSII appears to be an effective option for intensive insulin therapy in people with diabetes for improving suboptimal glycaemic control.

The association of polypharmacy and high-risk drug classes with adverse health outcomes in the Scottish population with type 1 diabetes.

AIMS/HYPOTHESIS: The aim of this work was to map the number of prescribed drugs over age, sex and area-based socioeconomic deprivation, and to examine the association between the number of drugs and particular high-risk drug classes with adverse health outcomes among a national cohort of individuals with type 1 diabetes. METHODS: Utilising linked healthcare records from the population-based diabetes register of Scotland, we identified 28,245 individuals with a diagnosis of type 1 diabetes on 1 January 2017. For this population, we obtained information on health status, predominantly reflecting diabetes-related complications, and information on the total number of drugs and particular high-risk drug classes prescribed. We then studied the association of these baseline-level features with hospital admissions for falls, diabetic ketoacidosis (DKA), and hypoglycaemia or death within the subsequent year using multivariate Cox proportional hazards models. RESULTS: Not considering insulin and treatment for hypoglycaemia, the mean number of prescribed drugs was 4.00 (SD 4.35). The proportion of individuals being prescribed five or more drugs at baseline consistently increased with age (proportion [95% CI]: 0-19 years 2.04% [1.60, 2.49]; 40-49 years 28.50% [27.08, 29.93]; 80+ years 76.04% [67.73, 84.84]). Controlling for age, sex, area-based socioeconomic deprivation and health status, each additional drug at baseline was associated with an increase in the hazard for hospitalisation for falls, hypoglycaemia and death but not for DKA admissions (HR [95% CI]: falls 1.03 [1.01, 1.06]; DKA 1.01 [1.00, 1.03]; hypoglycaemia 1.05 [1.02, 1.07]; death 1.04 [1.02, 1.06]). We found a number of drug classes to be associated with an increased hazard of one or more of these adverse health outcomes, including antithrombotic/anticoagulant agents, corticosteroids, opioids, antiepileptics, antipsychotics, hypnotics and sedatives, and antidepressants. CONCLUSIONS: Polypharmacy is common among the Scottish population with type 1 diabetes and is strongly patterned by sociodemographic factors. The number of prescribed drugs and the prescription of particular high-risk drug classes are strong markers of an increased risk of adverse health outcomes, including acute complications of diabetes.

Risks of and risk factors for COVID-19 disease in people with diabetes: a cohort study of the total population of Scotland.

BACKGROUND: We aimed to ascertain the cumulative risk of fatal or critical care unit-treated COVID-19 in people with diabetes and compare it with that of people without diabetes, and to investigate risk factors for and build a cross-validated predictive model of fatal or critical care unit-treated COVID-19 among people with diabetes. METHODS: In this cohort study, we captured the data encompassing the first wave of the pandemic in Scotland, from March 1, 2020, when the first case was identified, to July 31, 2020, when infection rates had dropped sufficiently that shielding measures were officially terminated. The participants were the total population of Scotland, including all people with diabetes who were alive 3 weeks before the start of the pandemic in Scotland (estimated Feb 7, 2020). We ascertained how many people developed fatal or critical care unit-treated COVID-19 in this period from the Electronic Communication of Surveillance in Scotland database (on virology), the RAPID database of daily hospitalisations, the Scottish Morbidity Records-01 of hospital discharges, the National Records of Scotland death registrations data, and the Scottish Intensive Care Society and Audit Group database (on critical care). Among people with fatal or critical care unit-treated COVID-19, diabetes status was ascertained by linkage to the national diabetes register, Scottish Care Information Diabetes. We compared the cumulative incidence of fatal or critical care unit-treated COVID-19 in people with and without diabetes using logistic regression. For people with diabetes, we obtained data on potential risk factors for fatal or critical care unit-treated COVID-19 from the national diabetes register and other linked health administrative databases. We tested the association of these factors with fatal or critical care unit-treated COVID-19 in people with diabetes, and constructed a prediction model using stepwise regression and 20-fold cross-validation. FINDINGS: Of the total Scottish population on March 1, 2020 (n=5 463 300), the population with diabetes was 319 349 (5·8%), 1082 (0·3%) of whom developed fatal or critical care unit-treated COVID-19 by July 31, 2020, of whom 972 (89·8%) were aged 60 years or older. In the population without diabetes, 4081 (0·1%) of 5 143 951 people developed fatal or critical care unit-treated COVID-19. As of July 31, the overall odds ratio (OR) for diabetes, adjusted for age and sex, was 1·395 (95% CI 1·304-1·494; p<0·0001, compared with the risk in those without diabetes. The OR was 2·396 (1·815-3·163; p<0·0001) in type 1 diabetes and 1·369 (1·276-1·468; p<0·0001) in type 2 diabetes. Among people with diabetes, adjusted for age, sex, and diabetes duration and type, those who developed fatal or critical care unit-treated COVID-19 were more likely to be male, live in residential care or a more deprived area, have a COVID-19 risk condition, retinopathy, reduced renal function, or worse glycaemic control, have had a diabetic ketoacidosis or hypoglycaemia hospitalisation in the past 5 years, be on more anti-diabetic and other medication (all p<0·0001), and have been a smoker (p=0·0011). The cross-validated predictive model of fatal or critical care unit-treated COVID-19 in people with diabetes had a C-statistic of 0·85 (0·83-0·86). INTERPRETATION: Overall risks of fatal or critical care unit-treated COVID-19 were substantially elevated in those with type 1 and type 2 diabetes compared with the background population. The risk of fatal or critical care unit-treated COVID-19, and therefore the need for special protective measures, varies widely among those with diabetes but can be predicted reasonably well using previous clinical history. FUNDING: None.

Disparities in the analysis of morphological disparity.

Analyses of morphological disparity have been used to characterize and investigate the evolution of variation in the anatomy, function and ecology of organisms since the 1980s. While a diversity of methods have been employed, it is unclear whether they provide equivalent insights. Here, we review the most commonly used approaches for characterizing and analysing morphological disparity, all of which have associated limitations that, if ignored, can lead to misinterpretation. We propose best practice guidelines for disparity analyses, while noting that there can be no 'one-size-fits-all' approach. The available tools should always be used in the context of a specific biological question that will determine data and method selection at every stage of the analysis.

Time trends in deaths before age 50 years in people with type 1 diabetes: a nationwide analysis from Scotland 2004-2017.

AIMS/HYPOTHESIS: We aimed to examine whether crude mortality and mortality relative to the general population below 50 years of age have improved in recent years in those with type 1 diabetes. METHODS: Individuals with type 1 diabetes aged below 50 and at least 1 year old at any time between 2004 and 2017 in Scotland were identified using the national register. Death data were obtained by linkage to Scottish national death registrations. Indirect age standardisation was used to calculate sex-specific standardised mortality ratios (SMRs). Poisson regression was used to test for calendar-time effects as incidence rate ratios (IRRs). RESULTS: There were 1138 deaths in 251,143 person-years among 27,935 people with type 1 diabetes. There was a significant decline in mortality rate over time (IRR for calendar year 0.983 [95% CI 0.967, 0.998], p = 0.03), but the SMR remained approximately stable at 3.1 and 3.6 in men and 4.09 and 4.16 in women for 2004 and 2017, respectively. Diabetic ketoacidosis or coma (DKAoC) accounted for 22% of deaths and the rate did not decline significantly (IRR 0.975 [95% CI 0.94, 1.011], p = 0.168); 79.3% of DKAoC deaths occurred out of hospital. Circulatory diseases accounted for 27% of deaths and did decline significantly (IRR 0.946 [95% CI 0.914, 0.979], p = 0.002). CONCLUSIONS/INTERPRETATION: Absolute mortality has fallen, but the relative impact of type 1 diabetes on mortality below 50 years has not improved. There is scope to improve prevention of premature circulatory diseases and DKAoC and to develop more effective strategies for enabling people with type 1 diabetes to avoid clinically significant hyper- or hypoglycaemia. Graphical abstract.

The Effect of Fossil Sampling on the Estimation of Divergence Times with the Fossilized Birth-Death Process.

Timescales are of fundamental importance to evolutionary biology as they facilitate hypothesis tests of historical evolutionary processes. Through the incorporation of fossil occurrence data, the fossilized birth-death (FBD) process provides a framework for estimating divergence times using more paleontological data than traditional node calibration approaches have allowed. The inclusion of more data can refine evolutionary timescale estimates, but for many taxonomic groups it is computationally infeasible to include all available fossil occurrence data. Here, we utilize both empirical data and a simulation framework to identify approaches to subsampling fossil occurrence data that result in the most accurate estimates of divergence times. To achieve this we assess the performance of the FBD-Skyline model when implementing multiple approaches to incorporating subsampled fossil occurrence data. Our results demonstrate that it is necessary to account for all available fossil occurrence data to achieve the most accurate estimates of clade age. We show that this can be achieved if an empirical Bayes approach, accounting for fossil sampling through time, is applied to the FBD process. Random subsampling of occurrence data can lead to estimates of clade age that are incompatible with fossil evidence if no control over the affinities of fossil occurrences is enforced. Our results call into question the accuracy of previous divergence time studies incorporating the FBD process that have used only a subsample of all available fossil occurrence data.

Probabilistic methods surpass parsimony when assessing clade support in phylogenetic analyses of discrete morphological data.

Fossil taxa are critical to inferences of historical diversity and the origins of modern biodiversity, but realizing their evolutionary significance is contingent on restoring fossil species to their correct position within the tree of life. For most fossil species, morphology is the only source of data for phylogenetic inference; this has traditionally been analysed using parsimony, the predominance of which is currently challenged by the development of probabilistic models that achieve greater phylogenetic accuracy. Here, based on simulated and empirical datasets, we explore the relative efficacy of competing phylogenetic methods in terms of clade support. We characterize clade support using bootstrapping for parsimony and Maximum Likelihood, and intrinsic Bayesian posterior probabilities, collapsing branches that exhibit less than 50% support. Ignoring node support, Bayesian inference is the most accurate method in estimating the tree used to simulate the data. After assessing clade support, Bayesian and Maximum Likelihood exhibit comparable levels of accuracy, and parsimony remains the least accurate method. However, Maximum Likelihood is less precise than Bayesian phylogeny estimation, and Bayesian inference recaptures more correct nodes with higher support compared to all other methods, including Maximum Likelihood. We assess the effects of these findings on empirical phylogenies. Our results indicate probabilistic methods should be favoured over parsimony.

The Efficacy of Consensus Tree Methods for Summarizing Phylogenetic Relationships from a Posterior Sample of Trees Estimated from Morphological Data.

Consensus trees are required to summarize trees obtained through MCMC sampling of a posterior distribution, providing an overview of the distribution of estimated parameters such as topology, branch lengths, and divergence times. Numerous consensus tree construction methods are available, each presenting a different interpretation of the tree sample. The rise of morphological clock and sampled-ancestor methods of divergence time estimation, in which times and topology are coestimated, has increased the popularity of the maximum clade credibility (MCC) consensus tree method. The MCC method assumes that the sampled, fully resolved topology with the highest clade credibility is an adequate summary of the most probable clades, with parameter estimates from compatible sampled trees used to obtain the marginal distributions of parameters such as clade ages and branch lengths. Using both simulated and empirical data, we demonstrate that MCC trees, and trees constructed using the similar maximum a posteriori (MAP) method, often include poorly supported and incorrect clades when summarizing diffuse posterior samples of trees. We demonstrate that the paucity of information in morphological data sets contributes to the inability of MCC and MAP trees to accurately summarise of the posterior distribution. Conversely, majority-rule consensus (MRC) trees represent a lower proportion of incorrect nodes when summarizing the same posterior samples of trees. Thus, we advocate the use of MRC trees, in place of MCC or MAP trees, in attempts to summarize the results of Bayesian phylogenetic analyses of morphological data.

Uncertain-tree: discriminating among competing approaches to the phylogenetic analysis of phenotype data.

Morphological data provide the only means of classifying the majority of life's history, but the choice between competing phylogenetic methods for the analysis of morphology is unclear. Traditionally, parsimony methods have been favoured but recent studies have shown that these approaches are less accurate than the Bayesian implementation of the Mk model. Here we expand on these findings in several ways: we assess the impact of tree shape and maximum-likelihood estimation using the Mk model, as well as analysing data composed of both binary and multistate characters. We find that all methods struggle to correctly resolve deep clades within asymmetric trees, and when analysing small character matrices. The Bayesian Mk model is the most accurate method for estimating topology, but with lower resolution than other methods. Equal weights parsimony is more accurate than implied weights parsimony, and maximum-likelihood estimation using the Mk model is the least accurate method. We conclude that the Bayesian implementation of the Mk model should be the default method for phylogenetic estimation from phenotype datasets, and we explore the implications of our simulations in reanalysing several empirical morphological character matrices. A consequence of our finding is that high levels of resolution or the ability to classify species or groups with much confidence should not be expected when using small datasets. It is now necessary to depart from the traditional parsimony paradigms of constructing character matrices, towards datasets constructed explicitly for Bayesian methods.

Tips and nodes are complementary not competing approaches to the calibration of molecular clocks.

Molecular clock methodology provides the best means of establishing evolutionary timescales, the accuracy and precision of which remain reliant on calibration, traditionally based on fossil constraints on clade (node) ages. Tip calibration has been developed to obviate undesirable aspects of node calibration, including the need for maximum age constraints that are invariably very difficult to justify. Instead, tip calibration incorporates fossil species as dated tips alongside living relatives, potentially improving the accuracy and precision of divergence time estimates. We demonstrate that tip calibration yields node calibrations that violate fossil evidence, contributing to unjustifiably young and ancient age estimates, less precise and (presumably) accurate than conventional node calibration. However, we go on to show that node and tip calibrations are complementary, producing meaningful age estimates, with node minima enforcing realistic ages and fossil tips interacting with node calibrations to objectively define maximum age constraints on clade ages. Together, tip and node calibrations may yield evolutionary timescales that are better justified, more precise and accurate than either calibration strategy can achieve alone.

Bayesian methods outperform parsimony but at the expense of precision in the estimation of phylogeny from discrete morphological data.

Different analytical methods can yield competing interpretations of evolutionary history and, currently, there is no definitive method for phylogenetic reconstruction using morphological data. Parsimony has been the primary method for analysing morphological data, but there has been a resurgence of interest in the likelihood-based Mk-model. Here, we test the performance of the Bayesian implementation of the Mk-model relative to both equal and implied-weight implementations of parsimony. Using simulated morphological data, we demonstrate that the Mk-model outperforms equal-weights parsimony in terms of topological accuracy, and implied-weights performs the most poorly. However, the Mk-model produces phylogenies that have less resolution than parsimony methods. This difference in the accuracy and precision of parsimony and Bayesian approaches to topology estimation needs to be considered when selecting a method for phylogeny reconstruction.

The Interrelationships of Placental Mammals and the Limits of Phylogenetic Inference.

Placental mammals comprise three principal clades: Afrotheria (e.g., elephants and tenrecs), Xenarthra (e.g., armadillos and sloths), and Boreoeutheria (all other placental mammals), the relationships among which are the subject of controversy and a touchstone for debate on the limits of phylogenetic inference. Previous analyses have found support for all three hypotheses, leading some to conclude that this phylogenetic problem might be impossible to resolve due to the compounded effects of incomplete lineage sorting (ILS) and a rapid radiation. Here we show, using a genome scale nucleotide data set, microRNAs, and the reanalysis of the three largest previously published amino acid data sets, that the root of Placentalia lies between Atlantogenata and Boreoeutheria. Although we found evidence for ILS in early placental evolution, we are able to reject previous conclusions that the placental root is a hard polytomy that cannot be resolved. Reanalyses of previous data sets recover Atlantogenata + Boreoeutheria and show that contradictory results are a consequence of poorly fitting evolutionary models; instead, when the evolutionary process is better-modeled, all data sets converge on Atlantogenata. Our Bayesian molecular clock analysis estimates that marsupials diverged from placentals 157-170 Ma, crown Placentalia diverged 86-100 Ma, and crown Atlantogenata diverged 84-97 Ma. Our results are compatible with placental diversification being driven by dispersal rather than vicariance mechanisms, postdating early phases in the protracted opening of the Atlantic Ocean.

Dating Tips for Divergence-Time Estimation.

The molecular clock is the only viable means of establishing an accurate timescale for Life on Earth, but it remains reliant on a capricious fossil record for calibration. 'Tip-dating' promises a conceptual advance, integrating fossil species among their living relatives using molecular/morphological datasets and evolutionary models. Fossil species of known age establish calibration directly, and their phylogenetic uncertainty is accommodated through the co-estimation of time and topology. However, challenges remain, including a dearth of effective models of morphological evolution, rate correlation, the non-random nature of missing characters in fossil data, and, most importantly, accommodating uncertainty in fossil age. We show uncertainty in fossil-dating propagates to divergence-time estimates, yielding estimates that are older and less precise than those based on traditional node calibration. Ultimately, node and tip calibrations are not mutually incompatible and may be integrated to achieve more accurate and precise evolutionary timescales.