MU variability in CBCT-guided online adaptive radiation therapy.

PURPOSE: CBCT-guided online-adaptive radiotherapy (oART) systems have been made possible by using artificial intelligence and automation to substantially reduce treatment planning time during on-couch adaptive sessions. Evaluating plans generated during an adaptive session presents significant challenges to the clinical team as the planning process gets compressed into a shorter window than offline planning. We identified MU variations up to 30% difference between the adaptive plan and the reference plan in several oART sessions that caused the clinical team to question the accuracy of the oART dose calculation. We investigated the cause of MU variation and the overall accuracy of the dose delivered when MU variations appear unnecessarily large. METHODS: Dosimetric and adaptive plan data from 604 adaptive sessions of 19 patients undergoing CBCT-guided oART were collected. The analysis included total MU per fraction, planning target volume (PTV) and organs at risk (OAR) volumes, changes in PTV-OAR overlap, and DVH curves. Sessions with MU greater than two standard deviations from the mean were reoptimized offline, verified by an independent calculation system, and measured using a detector array. RESULTS: MU variations relative to the reference plan were normally distributed with a mean of -1.0% and a standard deviation of 11.0%. No significant correlation was found between MU variation and anatomic changes. Offline reoptimization did not reliably reproduce either reference or on-couch total MUs, suggesting that stochastic effects within the oART optimizer are likely causing the variations. Independent dose calculation and detector array measurements resulted in acceptable agreement with the planned dose. CONCLUSIONS: MU variations observed between oART plans were not caused by any errors within the oART workflow. Providers should refrain from using MU variability as a way to express their confidence in the treatment planning accuracy. Clinical decisions during on-couch adaptive sessions should rely on validated secondary dose calculations to ensure optimal plan selection.

Case study with dosimetric analysis: Total body irradiation to a patient with a left ventricular assist device.

Key Clinical Message: A patient presented with cardiogenic shock, requiring the implantation of a left ventricular assist device (LVAD), and acute myeloblastic leukemia. This necessitated total body irradiation (TBI) while balancing dose reduction to the LVAD components to avoid potential radiation damage. Here we outline our treatment approach and dose estimates to the LVAD. Abstract: This case report discusses the delivery of TBI to a patient with an LVAD. This treatment required radiation-dose determinations and consequential reductions for the heart, LVAD, and an external controller connected to the LVAD. The patient was treated using a traditional 16MV anterior posterior (AP)/posterior anterior (PA) technique at a source-to-surface-distance of 515 cm for 400 cGy in two fractions. A 3 cm thick Cerrobend block was placed on the beam spoiler to reduce dose to the heart and LVAD to 150 cGy. The external controller was placed in a 1 cm thick acrylic box to reduce neutron dose and positioned as far from the treatment fields as achievable. In vivo measurements were made using optically stimulated luminescence dosimeters (OSLDs) placed inside the box at distances of 2 cm, 8.5 cm, and 14 cm from the field edge, and on the patient along the central axis and centered behind the LVAD block. Further ion chamber measurements were made using a solid water phantom to more accurately estimate the dose delivered to the LVAD. Neutron dose measurements were also conducted. The total estimated dose to the controller ranged from 135.3 cGy to 91.5 cGy. The LVAD block reduced the surface dose to the patient to 271.6 cGy (68.1%). The block transmission factors of the 3 cm Cerrobend block measured in the phantom were 45% at 1 cm depth and decreased asymptotically to around 30% at 3 cm depth. Applying these transmission factors to the in vivo measurements yielded a dose of 120 cGy to the implanted device. The neutron dose the LVAD region is estimated around 0.46 cGy. Physical limitations of the controller made it impossible to completely avoid dose. Shielding is recommended. The block had limited dose reduction to the surface, due to secondary particles, but appropriately reduced the dose at 3 cm and beyond. More research on LVADs dose limits would be beneficial.

Discovery of Small Molecule Interleukin 17A Inhibitors with Novel Binding Mode and Stoichiometry: Optimization of DNA-Encoded Chemical Library Hits to In Vivo Active Compounds.

Dysregulation of IL17A drives numerous inflammatory and autoimmune disorders with inhibition of IL17A using antibodies proven as an effective treatment. Oral anti-IL17 therapies are an attractive alternative option, and several preclinical small molecule IL17 inhibitors have previously been described. Herein, we report the discovery of a novel class of small molecule IL17A inhibitors, identified via a DNA-encoded chemical library screen, and their subsequent optimization to provide in vivo efficacious inhibitors. These new protein-protein interaction (PPI) inhibitors bind in a previously undescribed mode in the IL17A protein with two copies binding symmetrically to the central cavities of the IL17A homodimer.

Genomic variation in montane bumblebees in Scandinavia: High levels of intraspecific diversity despite population vulnerability.

Populations of many bumblebee species are declining, with distributions shifting northwards to track suitable climates. Climate change is considered a major contributing factor. Arctic species are particularly vulnerable as they cannot shift further north, making assessment of their population viability important. Analysis of levels of whole-genome variation is a powerful way to analyse population declines and fragmentation. Here, we use genome sequencing to analyse genetic variation in seven species of bumblebee from the Scandinavian mountains, including two classified as vulnerable. We sequenced 333 samples from across the ranges of these species in Sweden. Estimates of effective population size (NE ) vary from ~55,000 for species with restricted high alpine distributions to 220,000 for more widespread species. Population fragmentation is generally very low or undetectable over large distances in the mountains, suggesting an absence of barriers to gene flow. The relatively high NE and low population structure indicate that none of the species are at immediate risk of negative genetic effects caused by high levels of genetic drift. However, reconstruction of historical fluctuations in NE indicates that the arctic specialist species Bombus hyperboreus has experienced population declines since the last ice age and we detected one highly inbred diploid male of this species close to the southern limit of its range, potentially indicating elevated genetic load. Although the levels of genetic variation in montane bumblebee populations are currently relatively high, their ranges are predicted to shrink drastically due to the effects of climate change and monitoring is essential to detect future population declines.

Screening Incarcerated Women for Opioid Use Disorder.

The overall aim of the present study is to examine the utility of the DSM OUD Checklist and the NM-ASSIST screening tools to identify symptoms consistent with OUD among incarcerated women in county jails. This study contributes to the existing literature because research on screening and assessment approaches for incarcerated women has been limited. The focus of the current study is to describe the screening procedures and study recruitment for a larger parent study focused on increasing treatment linkages. Study findings indicate a positive correlation between indicators of OUD using the two screening tools, as well as a high degree of correlation between street opioid misuse and other high-risk drug indicators (overdose and injection practices). These findings underscore the importance of outreach, screening, and intervention in real-world settings, including jails, in order to increase access to OUD treatment among this vulnerable sample of women.

Targeting the Tyrosine Kinase 2 (TYK2) Pseudokinase Domain: Discovery of the Selective TYK2 Inhibitor ABBV-712.

Tyrosine kinase 2 (TYK2) is a nonreceptor tyrosine kinase that belongs to the JAK family also comprising JAK1, JAK2, and JAK3. TYK2 is an attractive target for various autoimmune diseases as it regulates signal transduction downstream of IL-23 and IL-12 receptors. Selective TYK2 inhibition offers a differentiated clinical profile compared to currently approved JAK inhibitors. However, selectivity for TYK2 versus other JAK family members has been difficult to achieve with small molecules that inhibit the catalytically active kinase domain. Successful targeting of the TYK2 pseudokinase domain as a strategy to achieve isoform selectivity was recently exemplified with deucravacitinib. Described herein is the optimization of selective TYK2 inhibitors targeting the pseudokinase domain, resulting in the discovery of the clinical candidate ABBV-712 (21).

The Genomic Basis of Adaptation to High Elevations in Africanized Honey Bees.

A range of different genetic architectures underpin local adaptation in nature. Honey bees (Apis mellifera) in the Eastern African Mountains harbor high frequencies of two chromosomal inversions that likely govern adaptation to this high-elevation habitat. In the Americas, honey bees are hybrids of European and African ancestries and adaptation to latitudinal variation in climate correlates with the proportion of these ancestries across the genome. It is unknown which, if either, of these forms of genetic variation governs adaptation in honey bees living at high elevations in the Americas. Here, we performed whole-genome sequencing of 29 honey bees from both high- and low-elevation populations in Colombia. Analysis of genetic ancestry indicated that both populations were predominantly of African ancestry, but the East African inversions were not detected. However, individuals in the higher elevation population had significantly higher proportions of European ancestry, likely reflecting local adaptation. Several genomic regions exhibited particularly high differentiation between highland and lowland bees, containing candidate loci for local adaptation. Genes that were highly differentiated between highland and lowland populations were enriched for functions related to reproduction and sperm competition. Furthermore, variation in levels of European ancestry across the genome was correlated between populations of honey bees in the highland population and populations at higher latitudes in South America. The results are consistent with the hypothesis that adaptation to both latitude and elevation in these hybrid honey bees are mediated by variation in ancestry at many loci across the genome.

A quick guide on implementing and quality assuring 3D printing in radiation oncology.

As three-dimensional (3D) printing becomes increasingly common in radiation oncology, proper implementation, usage, and ongoing quality assurance (QA) are essential. While there have been many reports on various clinical investigations and several review articles, there is a lack of literature on the general considerations of implementing 3D printing in radiation oncology departments, including comprehensive process establishment and proper ongoing QA. This review aims to guide radiation oncology departments in effectively using 3D printing technology for routine clinical applications and future developments. We attempt to provide recommendations on 3D printing equipment, software, workflow, and QA, based on existing literature and our experience. Specifically, we focus on three main applications: patient-specific bolus, high-dose-rate (HDR) surface brachytherapy applicators, and phantoms. Additionally, cost considerations are briefly discussed. This review focuses on point-of-care (POC) printing in house, and briefly touches on outsourcing printing via mail-order services.

Influence of air mapping errors on the dosimetric accuracy of prostate CBCT-guided online adaptive radiation therapy.

PURPOSE: CBCT-guided online adaptive radiotherapy (oART) plans presently utilize daily synthetic CTs (sCT) that are automatically generated using deformable registration algorithms. These algorithms may have poor performance at reproducing variable volumes of gas present during treatment. Therefore, we have analyzed the air mapping error between the daily CBCTs and the corresponding sCT and explored its dosimetric effect on oART plan calculation. METHODS: Abdominopelvic air volume was contoured on both the daily CBCT images and the corresponding synthetic images for 207 online adaptive pelvic treatments. Air mapping errors were tracked over all fractions. For two case studies representing worst case scenarios, dosimetric effects of air mapping errors were corrected in the sCT images using the daily CBCT air contours, then recalculating dose. Dose volume histogram statistics and 3D gamma passing rates were used to compare the original and air-corrected sCT-based dose calculations. RESULTS: All analyzed patients showed observable air pocket contour differences between the sCT and the CBCT images. The largest air volume difference observed in daily CBCT images for a given patient was 276.3 cc, a difference of more than 386% compared to the sCT. For the two case studies, the largest observed change in DVH metrics was a 2.6% reduction in minimum PTV dose, with all other metrics varying by less than 1.5%. 3D gamma passing rates using 1%/1 mm criteria were above 90% when comparing the uncorrected and corrected dose distributions. CONCLUSION: Current CBCT-based oART workflow can lead to inaccuracies in the mapping of abdominopelvic air pockets from daily CBCT to the sCT images used for the optimization and calculation of the adaptive plan. Despite the large observed mapping errors, the dosimetric effects of such differences on the accuracy of the adapted plan dose calculation are unlikely to cause differences greater than 3% for prostate treatments.

Interview-Based Patient- and Caregiver-Reported Experiences of Hunger and Improved Quality of Life with Setmelanotide Treatment in Bardet-Biedl Syndrome.

INTRODUCTION: Bardet-Biedl syndrome (BBS) is a rare genetic disease associated with hyperphagia, a pathologic insatiable hunger, due to impaired signaling in the melanocortin-4 receptor (MC4R) pathway. The impact of hyperphagia on the lives of patients with BBS and their families has not been fully characterized. METHODS: Patients with BBS or their caregivers who participated in clinical trials of the MC4R agonist setmelanotide (NCT03013543 and NCT03746522) were included in this qualitative study. Telephone interviews were conducted using a semistructured interview guide to explore patient experience and caregiver observations of hyperphagia before and during setmelanotide treatment. RESULTS: Nineteen interviews (8 patients, 11 caregivers) were conducted. The term "hunger" (rather than "hyperphagia") was used in interviews to ensure common terminology. Before setmelanotide treatment, all participants described their (or their child's) hunger as all-consuming, leading to an obsessive focus on food. Nine participants recalled intense, continuous hunger, and most participants (5 patients, 10 caregivers) reported lack of control with eating. Negative impacts on patients' lives included difficulties with concentration, emotional and physical manifestations, and impaired relationships. All participants experienced or observed improvements in hunger and health outcomes during treatment, the most meaningful of which included weight loss and decrease in obsessive focus on food and food-seeking behaviors. All participants reported improvements in either physical and/or emotional well-being and being satisfied with setmelanotide. CONCLUSIONS: Hyperphagia and resulting food-seeking behaviors have notable negative impacts on quality of life in patients with BBS and caregivers. Setmelanotide improved hyperphagia, reduced body weight and obsessive focus on food, and facilitated improvements in physical and emotional well-being for both patients and caregivers. TRIAL REGISTRATION: NCT03013543 and NCT03746522.

Population Genomics for Insect Conservation.

Insects constitute vital components of ecosystems. There is alarming evidence for global declines in insect species diversity, abundance, and biomass caused by anthropogenic drivers such as habitat degradation or loss, agricultural practices, climate change, and environmental pollution. This raises important concerns about human food security and ecosystem functionality and calls for more research to assess insect population trends and identify threatened species and the causes of declines to inform conservation strategies. Analysis of genetic diversity is a powerful tool to address these goals, but so far animal conservation genetics research has focused strongly on endangered vertebrates, devoting less attention to invertebrates, such as insects, that constitute most biodiversity. Insects' shorter generation times and larger population sizes likely necessitate different analytical methods and management strategies. The availability of high-quality reference genome assemblies enables population genomics to address several key issues. These include precise inference of past demographic fluctuations and recent declines, measurement of genetic load levels, delineation of evolutionarily significant units and cryptic species, and analysis of genetic adaptation to stressors. This enables identification of populations that are particularly vulnerable to future threats, considering their potential to adapt and evolve. We review the application of population genomics to insect conservation and the outlook for averting insect declines.

Genome Evolution of a Symbiont Population for Pathogen Defense in Honeybees.

The honeybee gut microbiome is thought to be important for bee health, but the role of the individual members is poorly understood. Here, we present closed genomes and associated mobilomes of 102 Apilactobacillus kunkeei isolates obtained from the honey crop (foregut) of honeybees sampled from beehives in Helsingborg in the south of Sweden and from the islands Gotland and Åland in the Baltic Sea. Each beehive contained a unique composition of isolates and repeated sampling of similar isolates from two beehives in Helsingborg suggests that the bacterial community is stably maintained across bee generations during the summer months. The sampled bacterial population contained an open pan-genome structure with a high genomic density of transposons. A subset of strains affiliated with phylogroup A inhibited growth of the bee pathogen Melissococcus plutonius, all of which contained a 19.5 kb plasmid for the synthesis of the antimicrobial compound kunkecin A, while a subset of phylogroups B and C strains contained a 32.9 kb plasmid for the synthesis of a putative polyketide antibiotic. This study suggests that the mobile gene pool of A. kunkeei plays a key role in pathogen defense in honeybees, providing new insights into the evolutionary dynamics of defensive symbiont populations.

A macroscopic link between interhemispheric tract myelination and cortico-cortical interactions during action reprogramming.

Myelination has been increasingly implicated in the function and dysfunction of the adult human brain. Although it is known that axon myelination shapes axon physiology in animal models, it is unclear whether a similar principle applies in the living human brain, and at the level of whole axon bundles in white matter tracts. Here, we hypothesised that in humans, cortico-cortical interactions between two brain areas may be shaped by the amount of myelin in the white matter tract connecting them. As a test bed for this hypothesis, we use a well-defined interhemispheric premotor-to-motor circuit. We combined TMS-derived physiological measures of cortico-cortical interactions during action reprogramming with multimodal myelin markers (MT, R1, R2* and FA), in a large cohort of healthy subjects. We found that physiological metrics of premotor-to-motor interaction are broadly associated with multiple myelin markers, suggesting interindividual differences in tract myelination may play a role in motor network physiology. Moreover, we also demonstrate that myelination metrics link indirectly to action switching by influencing local primary motor cortex dynamics. These findings suggest that myelination levels in white matter tracts may influence millisecond-level cortico-cortical interactions during tasks. They also unveil a link between the physiology of the motor network and the myelination of tracts connecting its components, and provide a putative mechanism mediating the relationship between brain myelination and human behaviour.

Massive genome inversion drives coexistence of divergent morphs in common quails.

The presence of population-specific phenotypes often reflects local adaptation or barriers to gene flow. The co-occurrence of phenotypic polymorphisms that are restricted within the range of a highly mobile species is more difficult to explain. An example of such polymorphisms is in the common quail Coturnix coturnix, a small migratory bird that moves widely during the breeding season in search of new mating opportunities, following ephemeral habitats,1,2 and whose females may lay successive clutches at different locations while migrating.3 In spite of this vagility, previous studies reported a higher frequency of heavier males with darker throat coloration in the southwest of the distribution (I. Jiménez-Blasco et al., 2015, Int. Union Game Biol., conference). We used population genomics and cytogenetics to explore the basis of this polymorphism and discovered a large inversion in the genome of the common quail. This inversion extends 115 Mbp in length and encompasses more than 7,000 genes (about 12% of the genome), producing two very different forms. Birds with the inversion are larger, have darker throat coloration and rounder wings, are inferred to have poorer flight efficiency, and are geographically restricted despite the high mobility of the species. Stable isotope analyses confirmed that birds carrying the inversion have shorter migratory distances or do not migrate. However, we found no evidence of pre- or post-zygotic isolation, indicating the two forms commonly interbreed and that the polymorphism remains locally restricted because of the effect on behavior. This illustrates a genomic mechanism underlying maintenance of geographically structured polymorphisms despite interbreeding with a lineage with high mobility.

A genomic and morphometric analysis of alpine bumblebees: Ongoing reductions in tongue length but no clear genetic component.

Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyse population structure in these populations. We generated a genome assembly of B. balteatus. We then analysed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species ("incognitus"). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait probably has a complex genetic and environmental basis.

Single-Crystal Bismuth Thiophosphate, BiPS4, as a Nontoxic and Mechanically Robust X-ray Detector.

Bismuth thiophosphate, BiPS4, is a promising nontoxic, high-density ternary chalcogenide semiconductor. Polycrystalline BiPS4 was synthesized from the stoichiometric melt of Bi, P, and S. Phosphorus was purified via an in-situ sublimation method. Single crystals of BiPS4 were grown using a modified vertical Bridgman method with a thermal gradient of 18 °C/cm. The material exhibits an electrical resistivity of 2 × 109 Ω·cm. The Knoop hardness of the single crystals is 128 ± 0.8 kg mm-2. A blocking contact behavior was observed with asymmetric contacts of Ga/BiPS4/Ag. A clear photocurrent response was observed from a BiPS4 crystal under an electric field as low as 1.14 V mm-1. Using a tungsten X-ray source, an X-ray sensitivity of 52 µ Gy-1 cm-2 was measured under an electric field of 80 V mm-1. When a single-crystal BiPS4 radiation detector device was used in a pulse-height radiation detection system, a clear charge collection response was observed under a 241Am α-particle radiation source.

New dosimetric guidelines for linear Boltzmann transport equations through comparative evaluation of stereotactic body radiation therapy for lung treatment planning.

PURPOSE: To propose guidelines for lung stereotactic body radiation therapy (SBRT) when using Acuros XB (AXB) equivalent to the existing ones developed for convolution algorithms such as analytic anisotropic algorithm (AAA), considering the difference between the algorithms. METHODS: A retrospective analysis was performed on 30 lung patients previously treated with SBRT. The original AAA plans, which were developed using dynamic conformal arcs, were recalculated and then renormalized for planning target volume (PTV) coverage using AXB. The recalculated and renormalized plans were compared to the original plans based on V100% and V90% PTV coverage, as well as V105%, conformality index, D2cm , Rx/Dmax , R50, and Dmin . These metrics were analyzed nominally and on variations according to RTOG and NRG guidelines. Based on the relative difference between each metric in the AAA and AXB plans, new guidelines were developed. The relative differences in our cohort were compared to previously documented AAA to AXB comparisons found in the literature. RESULTS: AAA plans recalculated in AXB had a significant reduction in most dosimetric metrics. The most notable changes were in V100% (4%) and the conformality index (7.5%). To achieve equal PTV coverage, AXB required an average of 1.8% more monitor units (MU). This fits well with previously published data. Applying the new guidelines to the AXB plans significantly increased the number of minor violations with no change in major violations, making them comparable to those of the original AAA plans. CONCLUSION: The relative difference found between AAA and AXB for SBRT lung plans has been shown to be consistent with previous works. Based on these findings, new guidelines for lung SBRT are recommended when planning with AXB.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients.

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

Informed breeding dispersal following stochastic changes to patch quality in a pond-breeding amphibian.

The unidirectional movement of animals between breeding patches (i.e. breeding dispersal) has profound implications for the ecological and evolutionary dynamics of spatially structured populations. In spatiotemporally variable environments, individuals are expected to adjust their dispersal decisions according to information gathered on the environmental and/or social cues that reflect the fitness prospects in a given breeding patch (i.e. informed dispersal). A paucity of empirical work limited our understanding of the ability of animals to depart from low-quality breeding patches and settle in high-quality breeding patches. We examined the capacity of individuals to respond to stochastic changes in habitat quality via informed breeding dispersal in a pond-breeding amphibian. We conducted a 5-year (2015-2019) capture-recapture study of boreal toads Anaxyrus boreas boreas (n = 1,100) that breed in beaver ponds in western Wyoming, USA. During early spring of 2017, an extreme flooding event destroyed several beaver dams and resulted in the loss of breeding habitat. We used multi-state models to investigate how temporal changes in pond characteristics influenced breeding dispersal, and determine whether movement decisions were in accordance with prospects for reproductive fitness. Boreal toads more often departed from low-quality breeding ponds (without successful metamorphosis) and settled in high-quality breeding ponds (with successful metamorphosis). Movement decisions were context-dependent and associated with pond characteristics altered by beaver dam destruction. Individuals were more likely to depart from shallow ponds with high vegetation cover and settle in deep ponds with low vegetation cover. The probability of metamorphosis was related to the same environmental cues, suggesting that boreal toads assess the fitness prospects of a breeding patch and adjust movement decisions accordingly (i.e. informed breeding dispersal). We demonstrated that stochastic variability in environmental conditions and habitat quality can underpin dispersal behaviour in amphibians. Our study highlighted the mechanistic linkages between habitat change, movement behaviour and prospects for reproductive performance, which is critical for understanding how wild animals respond to rapid environmental change.

Genetic Barriers to Historical Gene Flow between Cryptic Species of Alpine Bumblebees Revealed by Comparative Population Genomics.

Evidence is accumulating that gene flow commonly occurs between recently diverged species, despite the existence of barriers to gene flow in their genomes. However, we still know little about what regions of the genome become barriers to gene flow and how such barriers form. Here, we compare genetic differentiation across the genomes of bumblebee species living in sympatry and allopatry to reveal the potential impact of gene flow during species divergence and uncover genetic barrier loci. We first compared the genomes of the alpine bumblebee Bombus sylvicola and a previously unidentified sister species living in sympatry in the Rocky Mountains, revealing prominent islands of elevated genetic divergence in the genome that colocalize with centromeres and regions of low recombination. This same pattern is observed between the genomes of another pair of closely related species living in allopatry (B. bifarius and B. vancouverensis). Strikingly however, the genomic islands exhibit significantly elevated absolute divergence (dXY) in the sympatric, but not the allopatric, comparison indicating that they contain loci that have acted as barriers to historical gene flow in sympatry. Our results suggest that intrinsic barriers to gene flow between species may often accumulate in regions of low recombination and near centromeres through processes such as genetic hitchhiking, and that divergence in these regions is accentuated in the presence of gene flow.