Using oral bioaccessibility measurements to refine risk assessment of potentially toxic elements in topsoils across an urban area.

Elevated concentrations of As, Cr, Cu, Ni, Pb, V and Zn in topsoils in Belfast, Northern Ireland have been found to exceed assessment criteria in the city and therefore may pose a risk to human health. Most generic assessment criteria (GAC) for potentially toxic elements (PTEs) in soils assume PTEs are 100% bioavailable to humans. Here we use in-vitro oral bioaccessibility testing using the Unified BARGE method (UBM) to measure what proportion of soil contamination dissolves in the digestive tract and therefore is available for absorption by the body. This study considers how PTE bioaccessibility in soils varies spatially across urban areas and refines human health risk assessment for these PTEs using site specific oral bioaccessibility results to present the first regional assessment of risk that incorporates bioaccessibility testing. A total of 103 urban soil samples were selected for UBM testing. Results showed low bioaccessible fraction (BAF) for the PTEs from geogenic sources: Cr (0.45-5.9%), Ni (1.1-46.3%) and V (2.2-23.9%). Higher BAF values were registered for PTEs from anthropogenic sources: As (8.0-86.9%), Cu (3.4-67.8%), Pb (9.1-106.2%) and Zn (2.4-77.5%). Graphs of bioaccessibility adjusted assessment criteria (BAAC) were derived for each urban land use type and PTE. These provide a visual representation of the significance of oral bioaccessibility when deriving BAAC and how this is affected by 1) dominant exposure pathways for each land use and 2) relative harm posed from exposure to PTEs via each pathway, allowing oral bioaccessibility research to be targeted to contaminants and pathways that most significantly impact risk assessment. Pb was the most widespread contaminant with 16.5% of sites exceeding the Pb GAC. Applying BAAC did not significantly change risk evaluation for these samples as many had Pb BAF>50%. In contrast, all samples that exceeded the As GAC were found to no longer exceed a minimal level of risk when oral bioaccessibility was considered. Oral bioaccessibility testing resulted in a 45% reduction in the number of sites identified as posing a potential risk to human health.

Interest and Satisfaction of Telemedicine Use Among Ambulatory Neurology Patients in Western North Carolina During the COVID-19 Pandemic.

Introduction: During the COVID-19 pandemic, care shifted from exclusively telemedicine to hybrid models with in-person, video, and telephone visits. We explored how patient satisfaction and visit preferences have changed by comparing in-person versus virtual visits (telephone and video) in an ambulatory neurology practice across three time points. Methods: Patients who completed a virtual visit in March 2020 (early-pandemic), May 2020 (mid-pandemic), and March 2021 (later-pandemic) were contacted. Patients were assessed for visit satisfaction and desire for future telemedicine. Univariate and multivariable logistic regression analysis was conducted to determine factors independently associated with video visit completion. Results: Four thousand seven hundred seventy-eight the number of ambulatory visits (n = 4,778) were performed (1,004 early; 1,265 mid; and 2,509 later); 1,724 patients (36%) assented to postvisit feedback; mean age 45.8 ± 24.4 years, 58% female, 79% white, and 56% with Medicare/Medicaid insurance. Patient satisfaction significantly increased (73% early, 79% mid, 81% later-pandemic, p = 0.008). Interest in telemedicine also increased for patients completing telephone visits (40% early, 50% mid, 59% later, p = 0.027) and video visits (52% early, 59% mid, 62% later, p = 0.035). Patients satisfied with telemedicine visits were younger (p < 0.001). White patients were more interested in future telemedicine (p = 0.037). Multivariable analysis showed that older patients (for each 1 year older), Black patients, and patients with Medicare/Medicaid were 2%, 45%, and 54% less likely to complete a video visit than telephone, respectively. Discussion: Patients, especially younger ones, have become more satisfied and more interested in hybrid care models during the COVID-19 pandemic. Barriers to conducting video visits persist for older, Black patients with Medicare or Medicaid insurance.

A comparison of telemedicine and in-person neurology visits: what are the factors that patients consider when selecting future visit type?

OBJECTIVE: To identify factors that patients consider when choosing between future in-person, video, or telephone visits. BACKGROUND: Telemedicine has been rapidly integrated into ambulatory neurology in response to the COVID-19 pandemic. METHODS: Ambulatory neurology patients at a single center were contacted via telephone to complete: (1) a survey quantifying likelihood of scheduling a future telemedicine visit, and (2) a semi-structured qualitative interview following their visit in March 2021. Data were processed using the principles of thematic analysis. RESULTS: Of 2493 visits, 39% assented to post-visit feedback; 74% were in-person visits and 13% video and telephone. Patients with in-person visits were less likely than those with video and telephone visits to "definitely" consider a future telemedicine visit (36 vs. 59 and 62%, respectively; p < 0.001). Patients considered five key factors when scheduling future visits: "Pros of Visit Type," "Barriers to Telemedicine," "Situational Context," "Inherent Beliefs," and "Extrinsic Variables." Patients with telemedicine visits considered convenience as a pro, while those with in-person visits cited improved quality of care. Accessibility and user familiarity were considered barriers to telemedicine by patients with in-person and telephone visits, whereas system limitations were prevalent among patients with video visits. Patients agreed that stable conditions can be monitored via telemedicine, whereas physical examination warrants an in-person visit. Telemedicine was inherently considered equivalent to in-person care by patients with telephone visits. Awareness of telemedicine must be improved for patients with in-person visits. CONCLUSION: Across visit types, patients agree that telemedicine is convenient and effective in many circumstances. Future care delivery models should incorporate the patient perspective to implement hybrid models where telemedicine is an adjunct to in-person visits in ambulatory neurology.

Chronic kidney disease of unknown origin is associated with environmental urbanisation in Belfast, UK.

Chronic kidney disease (CKD), a collective term for many causes of progressive renal failure, is increasing worldwide due to ageing, obesity and diabetes. However, these factors cannot explain the many environmental clusters of renal disease that are known to occur globally. This study uses data from the UK Renal Registry (UKRR) including CKD of uncertain aetiology (CKDu) to investigate environmental factors in Belfast, UK. Urbanisation has been reported to have an increasing impact on soils. Using an urban soil geochemistry database of elemental concentrations of potentially toxic elements (PTEs), we investigated the association of the standardised incidence rates (SIRs) of both CKD and CKD of uncertain aetiology (CKDu) with environmental factors (PTEs), controlling for social deprivation. A compositional data analysis approach was used through balances (a special class of log contrasts) to identify elemental balances associated with CKDu. A statistically significant relationship was observed between CKD with the social deprivation measures of employment, income and education (significance levels of 0.001, 0.01 and 0.001, respectively), which have been used as a proxy for socio-economic factors such as smoking. Using three alternative regression methods (linear, generalised linear and Tweedie models), the elemental balances of Cr/Ni and As/Mo were found to produce the largest correlation with CKDu. Geogenic and atmospheric pollution deposition, traffic and brake wear emissions have been cited as sources for these PTEs which have been linked to kidney damage. This research, thus, sheds light on the increasing global burden of CKD and, in particular, the environmental and anthropogenic factors that may be linked to CKDu, particularly environmental PTEs linked to urbanisation.

Incorporating oral bioaccessibility into human health risk assessment due to potentially toxic elements in extractive waste and contaminated soils from an abandoned mine site.

The waste rock, tailings and soil around an abandoned mine site in Gorno (northwest Italy) contain elevated concentrations of potentially toxic elements (PTE) exceeding the permissible limits for residential uses. Specifically, the maximum concentrations of As, Cd, Pb, and Zn were 107 mg/kg, 340 mg/kg, 1064 mg/kg, and 148 433 mg/kg, respectively. A site-specific human health risk assessment (HHRA) was conducted for residential and recreational exposure scenarios, using an approach based on Risk Based Corrective Action (RBCA) method, refined by incorporating oral bioaccessibility data. Oral bioaccessibility analyses were performed by simulating the human digestion process in vitro (Unified BARGE Method). Detailed analysis of oral bioaccessible fraction (BAF i.e. ratio of bioaccessible concentrations to total concentrations on <250 µm fraction) indicated BAF of As (5-33%), Cd (72-98%), Co (24-42%), Cr (3-11%), Cu (25-90%), Ni (17-60%), Pb (16-88%) and Zn (73-94%). The solid phase distribution and mineralogical analyses showed that the variation of BAF is attributed to presence of alkaline calcareous rocks and association of PTE with a variety of minerals. The HHRA for ingestion pathway, suggested that bioaccessibility-corrected cancer risk reached up to 2.7 × 10-5 and 0.55 × 10-5 for residential and recreational senarios respectively (acceptable level is 1 × 10-5). The hazard index (HI) recalculated after incorporation of oral bioaccessible concentrations for a residential scenario ranged from 0.02 to 17.9. This was above the acceptable level (>1) for 50% of the samples, indicating potential human health risks. This study provides information for site-specific risk assessments and planning future research.

A comprehensive molecular phylogeny of Afrotropical white-eyes (Aves: Zosteropidae) highlights prior underestimation of mainland diversity and complex colonisation history.

White-eyes (Zosterops) are a hyper-diverse genus of passerine birds that have rapidly radiated across the Afrotropics and Southeast Asia. Despite their broad range, a disproportionately large number of species are currently recognised from islands compared to the mainland. Described species-level diversity of this 'great speciator' from continental Africa-Arabia is strikingly low, despite the vast size and environmental complexity of this region. However, efforts to identify natural groups using traditional approaches have been hindered by the remarkably uniform morphology and plumage of these birds. Here, we investigated the phylogenetic relationships and systematics of Afrotropical Zosterops, including the Gulf of Guinea and western Indian Ocean islands. We included exceptional sampling (~160 individuals) from all except one subspecies of the 55 taxa (32 species, plus 23 additional named sub-species) currently recognized throughout the region, in addition to a subset of extra-Afrotropical taxa, by exploiting blood and archival samples. Employing a multi-locus phylogenetic approach and applying quantitative species delimitation we tested: (1) if there has been a single colonisation event of the Afrotropical realm; (2) if constituent mainland and island birds are monophyletic; and (3) if mainland diversity has been underestimated. Our comprehensive regional phylogeny revealed a single recent colonisation of the Afrotropical realm c.1.30 Ma from Asia, but a subsequent complex colonisation history between constituent island and mainland lineages during their radiation across this vast area. Our findings suggest a significant previous underestimation of continental species diversity and, based on this, we propose a revised taxonomy. Our study highlights the need to densely sample species diversity across ranges, providing key findings for future conservation assessments and establishing a robust framework for evolutionary studies.

Apparent thinning of human visual cortex during childhood is associated with myelination.

Human cortex appears to thin during childhood development. However, the underlying microstructural mechanisms are unknown. Using functional magnetic resonance imaging (fMRI), quantitative MRI (qMRI), and diffusion MRI (dMRI) in children and adults, we tested what quantitative changes occur to gray and white matter in ventral temporal cortex (VTC) from childhood to adulthood, and how these changes relate to cortical thinning. T1 relaxation time from qMRI and mean diffusivity (MD) from dMRI provide independent and complementary measurements of microstructural properties of gray and white matter tissue. In face- and character-selective regions in lateral VTC, T1 and MD decreased from age 5 to adulthood in mid and deep cortex, as well as in their adjacent white matter. T1 reduction also occurred longitudinally in children's brain regions. T1 and MD decreases 1) were consistent with tissue growth related to myelination, which we verified with adult histological myelin stains, and 2) were correlated with apparent cortical thinning. In contrast, in place-selective cortex in medial VTC, we found no development of T1 or MD after age 5, and thickness was related to cortical morphology. These findings suggest that lateral VTC likely becomes more myelinated from childhood to adulthood, affecting the contrast of MR images and, in turn, the apparent gray-white boundary. These findings are important because they suggest that VTC does not thin during childhood but instead gets more myelinated. Our data have broad ramifications for understanding both typical and atypical brain development using advanced in vivo quantitative measurements and clinical conditions implicating myelin.

Linking oral bioaccessibility and solid phase distribution of potentially toxic elements in extractive waste and soil from an abandoned mine site: Case study in Campello Monti, NW Italy.

Oral bioaccessibility and solid phase distribution of potentially toxic elements (PTE) from extractive waste streams were investigated to assess the potential human health risk posed by abandoned mines. The solid phase distribution along with micro-X-ray fluorescence (micro-XRF) and scanning electron microscopy (SEM) analysis were also performed. The results showed that the total concentrations of PTE were higher in <250 µm size fractions of waste rock and soil samples in comparison to the <2 mm size fractions. Mean value of total concentrations of chromium(Cr), copper (Cu), and nickel (Ni) in waste rocks (size fractions <250 µm) were found to be 1299, 1570, and 4010 mg/kg respectively due to the parent material. However, only 11% of Ni in this sample was orally bioaccessible. Detailed analysis of the oral bioaccessible fraction (BAF, reported as the ratio of highest bioaccessible concentration compared with the total concentration from the 250 µm fraction) across all samples showed that Cr, Cu, and Ni varied from 1 to 6%, 14 to 47%, and 5 to 21%, respectively. The variation can be attributed to the difference in pH, organic matter content and mineralogical composition of the samples. Non-specific sequential extraction showed that the non-mobile forms of PTE were associated with the clay and Fe oxide components of the environmental matrices. The present study demonstrates how oral bioaccessibility, solid phase distribution and mineralogical analysis can provide insights into the distribution, fate and behaviour of PTE in waste streams from abandoned mine sites and inform human health risk posed by such sites .

Evidence for co-selection of antibiotic resistance genes and mobile genetic elements in metal polluted urban soils.

Antibiotic resistance genes (ARGs) constitute emerging environmental pollutants and pose risks to public health. Toxic metals are known to select for metal-resistant bacteria in metal-contaminated soils, but there is growing concern that metal contaminants can also act as co-selective agents thereby causing environmental proliferation of antibiotic resistance. In this study, we quantified ARGs and selected mobile genetic elements (MGEs) known to constitute potential ARG hosts in 50 archived urban and suburban soils from the Belfast metropolitan area using a high-throughput qPCR ARG chip. ARG prevalence was linked to concentrations of individual metals and a soil metal toxicity index calculated based on the relative toxicity of different metals to soil microbial processes. A total of 164 ARGs were detected across the 50 soils analyzed with an average absolute abundance of 3.4 × 107 ARG gene copies per gram of soil. A significant correlation between abundance of ARGs and MGEs was observed, suggesting the importance of horizontal gene transfer for ARG dissemination. Network analysis revealed significant co-occurrence patterns between specific metals (As, Cd, Co, Cr, Cu. Hg, Ni and Zn) and associated ARGs. Path analysis further indicated that the soil metal toxicity index significantly affected the number of detected ARGs (λ = 0.32, P < 0.001) and the abundance of metal co-occurring ARGs (λ = 0.612, P < 0.001) via effects on MGEs. Collectively, our results indicate a role of soil metals in co-selection of ARGs and MGEs in urban and semi-urban soils and suggest a risk for environmental ARG dissemination via horizontal gene transfer.

The relationship between historical development and potentially toxic element concentrations in urban soils.

Increasing urbanisation has a direct impact on soil quality, resulting in elevated concentrations of potentially toxic elements (PTEs) in soils. This research aims to assess if soil PTE concentrations can be used as an 'urbanisation tracer' by investigating geogenic and anthropogenic source contributions and controls, and considering PTE enrichment across historical urban development zones. The UK cities of Belfast and Sheffield are chosen as study areas, where available shallow and deep concentrations of PTEs in soil are compared to identify geogenic and anthropogenic contributions to PTEs. Cluster analysis and principal component analysis are used to elucidate the main controls over PTE concentrations. Pollution indices indicate that different periods of historical development are linked to enrichment of different PTEs. Urban subdomains are identified and background values calculated using various methodologies and compared to generic site assessment criteria. Exceedances for a number of the PTEs considered suggest a potential human health risk could be posed across subdomains of both Belfast and Sheffield. This research suggests that airborne diffuse contamination from often historical sources such as traffic, domestic combustion and industrial processes contribute greatly to soil contamination within urban environments. The relationship between historical development and differing PTEs is a novel finding, suggesting that PTEs have the potential for use as 'urbanisation tracers'. The investigative methodology employed has potential applications for decision makers, urban planners, regulators and developers of urban areas.

When are total concentrations not total? Factors affecting geochemical analytical techniques for measuring element concentrations in soil.

Inductively coupled plasma (ICP) following aqua regia digestion and X-ray fluorescence (XRF) are both geochemical techniques used to determine 'total' concentrations of elements in soil. The aim of this study is to compare these techniques, identify elements for which inconsistencies occur and investigate why they arise. A study area (∼14,000 km(2)) with a variety of total concentration controls and a large geochemical dataset (n = 7950) was selected. Principal component analysis determined underlying variance in a dataset composed of both geogenic and anthropogenic elements. Where inconsistencies between the techniques were identified, further numerical and spatial analysis was completed. The techniques are more consistent for elements of geogenic sources and lead, whereas other elements of anthropogenic sources show less consistency within rural samples. XRF is affected by sample matrix, while the form of element affects ICP concentrations. Depending on their use in environmental studies, different outcomes would be expected from the techniques employed, suggesting the choice of analytical technique for geochemical analyses may be more critical than realised.

The effects of lead sources on oral bioaccessibility in soil and implications for contaminated land risk management.

Lead (Pb) is a non-threshold toxin capable of inducing toxic effects at any blood level but availability of soil screening criteria for assessing potential health risks is limited. The oral bioaccessibility of Pb in 163 soil samples was attributed to sources through solubility estimation and domain identification. Samples were extracted following the Unified BARGE Method. Urban, mineralisation, peat and granite domains accounted for elevated Pb concentrations compared to rural samples. High Pb solubility explained moderate-high gastric (G) bioaccessible fractions throughout the study area. Higher maximum G concentrations were measured in urban (97.6 mg kg(-1)) and mineralisation (199.8 mg kg(-1)) domains. Higher average G concentrations occurred in mineralisation (36.4 mg kg(-1)) and granite (36.0 mg kg(-1)) domains. Findings suggest diffuse anthropogenic and widespread geogenic contamination could be capable of presenting health risks, having implications for land management decisions in jurisdictions where guidance advises these forms of pollution should not be regarded as contaminated land.

Niche divergence promotes rapid diversification of East African sky island white-eyes (Aves: Zosteropidae).

The Eastern Afromontane biodiversity hotspot composed of highly fragmented forested highlands (sky islands) harbours exceptional diversity and endemicity, particularly within birds. To explain their elevated diversity within this region, models founded on niche conservatism have been offered, although detailed phylogeographic studies are limited to a few avian lineages. Here, we focus on the recent songbird genus Zosterops, represented by montane and lowland members, to test the roles of niche conservatism versus niche divergence in the diversification and colonization of East Africa's sky islands. The species-rich white-eyes are a typically homogeneous family with an exceptional colonizing ability, but in contrast to their diversity on oceanic islands, continental diversity is considered depauperate and has been largely neglected. Molecular phylogenetic analysis of ~140 taxa reveals extensive polyphyly among different montane populations of Z. poliogastrus. These larger endemic birds are shown to be more closely related to taxa with divergent habitat types, altitudinal distributions and dispersal abilities than they are to populations of restricted endemics that occur in neighbouring montane forest fragments. This repeated transition between lowland and highland habitats over time demonstrate that diversification of the focal group is explained by niche divergence. Our results also highlight an underestimation of diversity compared to morphological studies that has implications for their taxonomy and conservation. Molecular dating suggests that the spatially extensive African radiation arose exceptionally rapidly (1-2.5 Ma) during the fluctuating Plio-Pleistocene climate, which may have provided the primary driver for lineage diversification.

Comparison of methods used to calculate typical threshold values for potentially toxic elements in soil.

The environmental quality of land can be assessed by calculating relevant threshold values, which differentiate between concentrations of elements resulting from geogenic and diffuse anthropogenic sources and concentrations generated by point sources of elements. A simple process allowing the calculation of these typical threshold values (TTVs) was applied across a region of highly complex geology (Northern Ireland) to six elements of interest; arsenic, chromium, copper, lead, nickel and vanadium. Three methods for identifying domains (areas where a readily identifiable factor can be shown to control the concentration of an element) were used: k-means cluster analysis, boxplots and empirical cumulative distribution functions (ECDF). The ECDF method was most efficient at determining areas of both elevated and reduced concentrations and was used to identify domains in this investigation. Two statistical methods for calculating normal background concentrations (NBCs) and upper limits of geochemical baseline variation (ULBLs), currently used in conjunction with legislative regimes in the UK and Finland respectively, were applied within each domain. The NBC methodology was constructed to run within a specific legislative framework, and its use on this soil geochemical data set was influenced by the presence of skewed distributions and outliers. In contrast, the ULBL methodology was found to calculate more appropriate TTVs that were generally more conservative than the NBCs. TTVs indicate what a "typical" concentration of an element would be within a defined geographical area and should be considered alongside the risk that each of the elements pose in these areas to determine potential risk to receptors.

The importance of solid-phase distribution on the oral bioaccessibility of Ni and Cr in soils overlying Palaeogene basalt lavas, Northern Ireland.

Potentially toxic elements (PTEs) including nickel and chromium are often present in soils overlying basalt at concentrations above regulatory guidance values due to the presence of these elements in underlying geology. Oral bioaccessibility testing allows the risk posed by PTEs to human health to be assessed; however, bioaccessibility is controlled by factors including mineralogy, particle size, solid-phase speciation and encapsulation. X-ray diffraction was used to characterise the mineralogy of 12 soil samples overlying Palaeogene basalt lavas in Northern Ireland, and non-specific sequential extraction coupled with chemometric analysis was used to determine the distribution of elements amongst soil components in 3 of these samples. The data obtained were related to total concentration and oral bioaccessible concentration to determine whether a relationship exists between the overall concentrations of PTEs, their bioaccessibility and the soils mineralogy and geochemistry. Gastric phase bioaccessible fraction (BAF %) ranged from 0.4 to 5.4 % for chromium in soils overlying basalt and bioaccessible and total chromium concentrations are positively correlated. In contrast, the range of gastric phase BAF for nickel was greater (1.4-43.8 %), while no significant correlation was observed between bioaccessible and total nickel concentrations. However, nickel BAF was inversely correlated with total concentration. Solid-phase fractionation information showed that bioaccessible nickel was associated with calcium carbonate, aluminium oxide, iron oxide and clay-related components, while bioaccessible chromium was associated with clay-related components. This suggests that weathering significantly affects nickel bioaccessibility, but does not have the same effect on the bioaccessibility of chromium.

Correlation analysis as a tool to investigate the bioaccessibility of nickel, vanadium and zinc in Northern Ireland soils.

Correlation analyses were conducted on nickel (Ni), vanadium (V) and zinc (Zn) oral bioaccessible fractions (BAFs) and selected geochemistry parameters to identify specific controls exerted over trace element bioaccessibility. BAFs were determined by previous research using the unified BARGE method. Total trace element concentrations and soil geochemical parameters were analysed as part of the Geological Survey of Northern Ireland Tellus Project. Correlation analysis included Ni, V and Zn BAFs against their total concentrations, pH, estimated soil organic carbon (SOC) and a further eight element oxides. BAF data were divided into three separate generic bedrock classifications of basalt, lithic arenite and mudstone prior to analysis, resulting in an increase in average correlation coefficients between BAFs and geochemical parameters. Sulphur trioxide and SOC, spatially correlated with upland peat soils, exhibited significant positive correlations with all BAFs in gastric and gastro-intestinal digestion phases, with such effects being strongest in the lithic arenite bedrock group. Significant negative relationships with bioaccessible Ni, V and Zn and their associated total concentrations were observed for the basalt group. Major element oxides were associated with reduced oral trace element bioaccessibility, with Al2O3 resulting in the highest number of significant negative correlations followed by Fe2O3. spatial mapping showed that metal oxides were present at reduced levels in peat soils. The findings illustrate how specific geology and soil geochemistry exert controls over trace element bioaccessibility, with soil chemical factors having a stronger influence on BAF results than relative geogenic abundance. In general, higher Ni, V and Zn bioaccessibility is expected in peat soil types.

Divergence times and colonization of the Canary Islands by Gallotia lizards.

The Canary Islands have become a model region for evolutionary studies. We obtained 1.8 Kbp of mtDNA sequence from all known island forms of the endemic lizard genus Gallotia and from its sister taxon Psammodromus in order to reanalyze phylogenetic relationships within the archipelago, estimate lineage divergence times, and reconstruct the colonization history of this group. Well-supported phylogenies were obtained using maximum parsimony and Bayesian inference. Previous studies have been unable to establish the branching pattern at the base of the tree. We found evidence that G. stehlini (Gran Canaria) originated from the most basal Gallotia node and G. atlantica from the subsequent node. Divergence times were estimated under a global clock using Bayesian Markov Chain Monte Carlo methods implemented by three different programs: BEAST, MCMCTREE, MULTIDIVTIME. Node constraints were derived from subaerial island appearance data and were incorporated into the analyses as soft or hard maximal bounds. Posterior node ages differed slightly between programs, possibly due to different priors on divergence times. The most eastern Canary Islands first emerged just over 20 mya and their colonization appears to have taken place relatively quickly, around 17-20 mya. The subsequent node is consistent with cladogenesis due to colonization of Gran Canaria from the eastern islands about 11-13 mya. The western islands appear to have been colonized by a dispersal event from Lanzarote/Fuerteventura in the east to either La Gomera or one of the ancient edifices that subsequently formed Tenerife in the west, about 9-10 mya. Within the western islands, the most recent node that is ancestral to both the G. intermedia/G. gomerana/G. simonyi and the G.galloti/G. caesaris clades is dated at about 5-6 mya. Subsequent dispersal events between ancient Tenerife islands and La Gomera are dated at around 3 mya in both clades, although the direction of dispersal cannot be determined. Finally, we show that G. galloti is likely to have colonized La Palma more than 0.5 Ma after emergence of the island 1.77 mya, while G. caesaris from the same clade may have colonized El Hierro very soon after it emerged 1.12 mya. There are tentative indications that the large-bodied endangered G. simonyi colonized El Hierro around the same time or even later than the smaller-bodied G. caesaris. This study demonstrates the effectiveness of Bayesian dating of a phylogeny in helping reconstruct the historical pattern of dispersal across an oceanic archipelago.