Automated determination of the ion-recombination correction factor (ksat) in ultra-high dose rate electron radiation therapy

BACKGROUND

Plane-parallel ionization chambers are the recommended secondary standard systems for clinical reference dosimetry of electrons. Dosimetry in high dose rate and dose-per-pulse (DPP) is challenging as ionization chambers are subject to ion recombination, especially when dose rate and/or DPP is increased beyond the range of conventional radiotherapy. The lack of universally accepted models for correction of ion recombination in UDHR is still an issue as it is, especially in FLASH-RT research, which is crucial in order to be able to accurately measure the dose for a wide range of dose rates and DPPs.

PURPOSE

The objective of this study was to show the feasibility of developing an Artificial Intelligence model to predict the ion-recombination factor-ksat for a plane-parallel Advanced Markus ionization chamber for conventional and ultra-high dose rate electron beams based on machine parameters. In addition, the predicted ksat of the AI model was compared with the current applied analytical models for this correction factor.

METHODS

A total number of 425 measurements was collected with a balanced variety in machine parameter settings. The specific ksat values were determined by dividing the output of the reference dosimeter (optically stimulated luminescence [OSL]) by the output of the AM chamber. Subsequently, a XGBoost regression model was trained, which used the different machine parameters as input features and the corresponding ksat value as output. The prediction accuracy of this regression model was characterized by R2-coefficient of determination, mean absolute error and root mean squared error. In addition, the model was compared with the Two-Voltage (TVA) method and empirical Petersson model for 19 different dose-per-pulse values ranging from conventional to UDHR regimes. The Akiake Information criterion (AIC) was calculated for the three different models.

RESULTS

The XGBoost regression model reached a R2-score of 0.94 on the independent test set with a MAE of 0.067 and RMSE of 0.106. For the additional 19 random data points, the ksat values predicted by the XGBoost model showed to be in agreement, within the uncertainties, with the ones determined by the Petersson model and better than the TVA method for doses per pulse >3.5 Gy with a maximum deviation from the ground truth of 14.2%, 16.7%, and -36.0%, respectively, for DPP >4 Gy.

CONCLUSION

The proposed method of using AI for ksat determination displays efficiency. For the investigated DPPs, the ksat values obtained with the XGBoost model were in concurrence with the ones obtained with the current available analytical models within the boundaries of uncertainty, certainly for the DPP characterizing UDHR. But the overall performance of the AI model, taking the number of free parameters into account, lacked efficiency. Future research should optimize the determination of the experimental ksat, and investigate the determination the ksat for DPPs higher than the ones investigated in this study, while also evaluating the prediction of the proposed XGBoost model for UDHR machines of different centers.

Longitudinal alterations in brain perfusion and vascular reactivity in the zQ175DN mouse model of Huntington's disease

BACKGROUND

Huntington's disease (HD) is marked by a CAG-repeat expansion in the huntingtin gene that causes neuronal dysfunction and loss, affecting mainly the striatum and the cortex. Alterations in the neurovascular coupling system have been shown to lead to dysregulated energy supply to brain regions in several neurological diseases, including HD, which could potentially trigger the process of neurodegeneration. In particular, it has been observed in cross-sectional human HD studies that vascular alterations are associated to impaired cerebral blood flow (CBF). To assess whether whole-brain changes in CBF are present and follow a pattern of progression, we investigated both resting-state brain perfusion and vascular reactivity longitudinally in the zQ175DN mouse model of HD.

METHODS

Using pseudo-continuous arterial spin labelling (pCASL) MRI in the zQ175DN model of HD and age-matched wild-type (WT) mice, we assessed whole-brain, resting-state perfusion at 3, 6 and 9 and 13 months of age, and assessed hypercapnia-induced cerebrovascular reactivity (CVR), at 4.5, 6, 9 and 15 months of age.

RESULTS

We found increased perfusion in cortical regions of zQ175DN HET mice at 3 months of age, and a reduction of this anomaly at 6 and 9 months, ages at which behavioural deficits have been reported. On the other hand, under hypercapnia, CBF was reduced in zQ175DN HET mice as compared to the WT: for multiple brain regions at 6 months of age, for only somatosensory and retrosplenial cortices at 9 months of age, and brain-wide by 15 months. CVR impairments in cortical regions, the thalamus and globus pallidus were observed in zQ175DN HET mice at 9 months, with whole brain reactivity diminished at 15 months of age. Interestingly, blood vessel density was increased in the motor cortex at 3 months, while average vessel length was reduced in the lateral portion of the caudate putamen at 6 months of age.

CONCLUSION

Our findings reveal early cortical resting-state hyperperfusion and impaired CVR at ages that present motor anomalies in this HD model, suggesting that further characterization of brain perfusion alterations in animal models is warranted as a potential therapeutic target in HD.

Intraoperative electron beam intercomparison of 6 sites using mailed thermoluminescence dosimetry: Absolute dose and energy

PURPOSE

In 2018, the Netherlands Commission on Radiation Dosimetry subcommittee on IORT initiated a limited intercomparison of electron IORT (IOERT) in Belgium and The Netherlands. The participating institutions have enough variability in age, type of equipment, and in dose calibration protocols.

METHODS

In this study, three types of IOERT-dedicated mobile accelerators were represented: Mobetron 2000, LIAC HWL and LIAC. Mobetron produces electron beams with energies of 6, 9 and 12 MeV, while LIAC HWL and LIAC can deliver 6, 8, 10 and 12 MeV electron beams. For all energies, the reference beam (10 cm diameter, 0° incidence) and 5 cm diameter beams were measured, as these smaller beams are used more frequently in clinic. The mailed TLD service from the Radiation Dosimetry Services (RDS, Houston, USA) has been used. Following RDS' standard procedures, each beam was irradiated to 300 cGy at dmax with TLDs around dmax and around depth of 50 % dose (R50). Absolute dose at 100 % and beam energy, expressed as R50, could be verified in this way.

RESULTS

All absolute doses and energies under reference conditions were well within RDS-stated uncertainties: dose deviations were <5 % and deviations in R50 were <5 mm. For the small 5 cm beams, all results were also within acceptance levels except one absolute dose value. Deviations were not significantly dependent on manufacturer, energy, diameter and calibration protocol.

CONCLUSIONS

All absolute dose values, except one of a non-reference beam, and all energy values were well within the measurement accuracy of RDS TLDs.

Multi-institutional generalizability of a plan complexity machine learning model for predicting pre-treatment quality assurance results in radiotherapy

Background and purpose

Treatment plans in radiotherapy are subject to measurement-based pre-treatment verifications. In this study, plan complexity metrics (PCMs) were calculated per beam and used as input features to develop a predictive model. The aim of this study was to determine the robustness against differences in machine type and institutional-specific quality assurance (QA).

Material and methods

A number of 567 beams were collected, where 477 passed and 90 failed the pre-treatment QA. Treatment plans of different anatomical regions were included. One type of linear accelerator was represented. For all beams, 16 PCMs were calculated. A random forest classifier was trained to distinct between acceptable and non-acceptable beams. The model was validated on other datasets to investigate its robustness. Firstly, plans for another machine type from the same institution were evaluated. Secondly, an inter-institutional validation was conducted on three datasets from different centres with their associated QA.

Results

Intra-institutionally, the PCMs beam modulation, mean MLC gap, Q1 gap, and Modulation Complexity Score were the most informative to detect failing beams. Eighty-tree percent of the failed beams (15/18) were detected correctly. The model could not detect over-modulated beams of another machine type. Inter-institutionally, the model performance reached higher accuracy for centres with comparable equipment both for treatment and QA as the local institute.

Conclusions

The study demonstrates that the robustness decreases when major differences appear in the QA platform or in planning strategies, but that it is feasible to extrapolate institutional-specific trained models between centres with similar clinical practice. Predictive models should be developed for each machine type.

Optically stimulated luminescence system as an alternative for radiochromic film for 2D reference dosimetry in UHDR electron beams

Radiotherapy is part of the treatment of over 50% of cancer patients. Its efficacy is limited by the radiotoxicity to the healthy tissue. FLASH-RT is based on the biological effect that ultra-high dose rates (UHDR) and very short treatment times strongly reduce normal tissue toxicity, while preserving the anti-tumoral effect. Despite many positive preclinical results, the translation of FLASH-RT to the clinic is hampered by the lack of accurate dosimetry for UHDR beams. To date radiochromic film is commonly used for dose assessment but has the drawback of lengthy and cumbersome read out procedures. In this work, we investigate the equivalence of a 2D OSL system to radiochromic film dosimetry in terms of dose rate independency. The comparison of both systems was done using the ElectronFlash linac. We investigated the dose rate dependence by variation of the (1) modality, (2) pulse repetition frequency, (3) pulse length and (4) source to surface distance. Additionally, we compared the 2D characteristics by field size measurements. The OSL calibration showed transferable between conventional and UHDR modality. Both systems are equally independent of average dose rate, pulse length and instantaneous dose rate. The OSL system showed equivalent in field size determination within 3 sigma. We show the promising nature of the 2D OSL system to serve as alternative for radiochromic film in UHDR electron beams. However, more in depth characterization is needed to assess its full potential.

Altered dynamics of glymphatic flow in a mature-onset Tet-off APP mouse model of amyloidosis

BACKGROUND

Alzheimer's disease (AD) is an incurable neurodegenerative disorder characterised by the progressive buildup of toxic amyloid-beta (Aβ) and tau protein aggregates eventually leading to cognitive decline. Recent lines of evidence suggest that an impairment of the glymphatic system (GS), a brain waste clearance pathway, plays a key role in the pathology of AD. Moreover, a relationship between GS function and neuronal network integrity has been strongly implicated. Here, we sought to assess the efficacy of the GS in a transgenic Tet-Off APP mouse model of amyloidosis, in which the expression of mutant APP was delayed until maturity, mimicking features of late-onset AD-the most common form of dementia in humans.

METHODS

To evaluate GS function, we used dynamic contrast-enhanced MRI (DCE-MRI) in 14-month-old Tet-Off APP (AD) mice and aged-matched littermate controls. Brain-wide transport of the Gd-DOTA contrast agent was monitored over time after cisterna magna injection. Region-of-interest analysis and computational modelling were used to assess GS dynamics while characterisation of brain tissue abnormalities at the microscale was performed ex vivo by immunohistochemistry.

RESULTS

We observed reduced rostral glymphatic flow and higher accumulation of the contrast agent in areas proximal to the injection side in the AD group. Clustering and subsequent computational modelling of voxel time courses revealed significantly lower influx time constants in AD relative to the controls. Ex vivo evaluation showed abundant amyloid plaque burden in the AD group coinciding with extensive astrogliosis and microgliosis. The neuroinflammatory responses were also found in plaque-devoid regions, potentially impacting brain-fluid circulation.

CONCLUSIONS

In a context resembling late-onset AD in humans, we demonstrate the disruption of glymphatic function and particularly a reduction in brain-fluid influx in the AD group. We conjecture that the hindered circulation of cerebrospinal fluid is potentially caused by wide-spread astrogliosis and amyloid-related obstruction of the normal routes of glymphatic flow resulting in redirection towards caudal regions. In sum, our study highlights the translational potential of alternative approaches, such as targeting brain-fluid circulation as potential therapeutic strategies for AD.

A new solution for UHDP and UHDR (Flash) measurements: Theory and conceptual design of ALLS chamber

Ultra-High dose-per-pulse regimens (UHDP), necessary to trigger the "FLASH" effect, still pose serious challenges to dosimetry. Dosimetry plays a crucial role, both to significantly improve the accuracy of the radiobiological experiments necessary to fully understand the mechanisms underlying the effect and its dependencies on the beam parameters, and to be able to translate such effect into clinical practice. The standard ionization chamber in UHDP region is significantly affected by the effects of the electric field generated by the enormous density of charges produced by the dose pulse. This work describes the theory and the conceptual design of a gas chamber (the ALLS chamber) which overcomes the above-mentioned problems.

Application of a novel diamond detector for commissioning of FLASH radiotherapy electron beams

Purpose

A diamond detector prototype was recently proposed by Marinelli et al. (Medical Physics 2022, https://doi.org/10.1002/mp.15473) for applications in ultrahigh‐dose‐per‐pulse (UH‐DPP) and ultrahigh‐dose‐rate (UH‐DR) beams, as used in FLASH radiotherapy (FLASH‐RT). In the present study, such so‐called flashDiamond (fD) was investigated from the dosimetric point of view, under pulsed electron beam irradiation. It was then used for the commissioning of an ElectronFlash linac (SIT S.p.A., Italy) both in conventional and UH‐DPP modalities.

Methods

Detector calibration was performed in reference conditions, under ⁶⁰Co and electron beam irradiation. Its response linearity was investigated in UH‐DPP conditions. For this purpose, the DPP was varied in the 1.2–11.9 Gy range, by changing either the beam applicator or the pulse duration from 1 to 4 μs. Dosimetric validation of the fD detector prototype was then performed in conventional modality, by measuring percentage depth dose (PDD) curves, beam profiles, and output factors (OFs). All such measurements were carried out in a motorized water phantom. The obtained results were compared with the ones from commercially available dosimeters, namely, a microDiamond, an Advanced Markus ionization chamber, a silicon diode detector, and EBT‐XD GAFchromic films. Finally, the fD detector was used to fully characterize the 7 and 9 MeV UH‐DPP electron beams delivered by the ElectronFlash linac. In particular, PDDs, beam profiles, and OFs were measured, for both energies and all the applicators, and compared with the ones from EBT‐XD films irradiated in the same experimental conditions.

Results

The fD calibration coefficient resulted to be independent from the investigated beam qualities. The detector response was found to be linear in the whole investigated DPP range. A very good agreement was observed among PDDs, beam profiles, and OFs measured by the fD prototype and reference detectors, both in conventional and UH‐DPP irradiation modalities.

Conclusions

The fD detector prototype was validated from the dosimetric point of view against several commercial dosimeters in conventional beams. It was proved to be suitable in UH‐DPP and UH‐DR conditions, for which no other commercial real‐time active detector is available to date. It was shown to be a very useful tool to perform fast and reproducible beam characterizations in standard clinical motorized water phantom setups. All of the previously mentioned demonstrate the suitability of the proposed detector for the commissioning of UH‐DR linac beams for preclinical FLASH‐RT applications.

Machine learning-based detection of aberrant deep learning segmentations of target and organs at risk for prostate radiotherapy using a secondary segmentation algorithm

Objective.The output of a deep learning (DL) auto-segmentation application should be reviewed, corrected if needed and approved before being used clinically. This verification procedure is labour-intensive, time-consuming and user-dependent, which potentially leads to significant errors with impact on the overall treatment quality. Additionally, when the time needed to correct auto-segmentations approaches the time to delineate target and organs at risk from scratch, the usability of the DL model can be questioned. Therefore, an automated quality assurance framework was developed with the aim to detect in advance aberrant auto-segmentations.Approach. Five organs (prostate, bladder, anorectum, femoral head left and right) were auto-delineated on CT acquisitions for 48 prostate patients by an in-house trained primary DL model. An experienced radiation oncologist assessed the correctness of the model output and categorised the auto-segmentations into two classes whether minor or major adaptations were needed. Subsequently, an independent, secondary DL model was implemented to delineate the same structures as the primary model. Quantitative comparison metrics were calculated using both models' segmentations and used as input features for a machine learning classification model to predict the output quality of the primary model.Main results. For every organ, the approach of independent validation by the secondary model was able to detect primary auto-segmentations that needed major adaptation with high sensitivity (recall = 1) based on the calculated quantitative metrics. The surface DSC and APL were found to be the most indicated parameters in comparison to standard quantitative metrics for the time needed to adapt auto-segmentations.Significance. This proposed method includes a proof of concept for the use of an independent DL segmentation model in combination with a ML classifier to improve time saving during QA of auto-segmentations. The integration of such system into current automatic segmentation pipelines can increase the efficiency of the radiotherapy contouring workflow.

Development of an ultra-thin parallel plate ionization chamber for dosimetry in flash radiotherapy

Background

Conventional air ionization chambers (ICs) exhibit ion recombination correction factors that deviate substantially from unity when irradiated with dose per pulse magnitudes higher than those used in conventional radiotherapy. This fact makes these devices unsuitable for the dosimetric characterization of beams in ultra‐high dose per pulse as used for FLASH radiotherapy.

Purpose

We present the design, development, and characterization of an ultra‐thin parallel plate IC that can be used in ultra‐high dose rate (UHDR) deliveries with minimal recombination.

Methods

The charge collection efficiency (CCE) of parallel plate ICs was modeled through a numerical solution of the coupled differential equations governing the transport of charged carriers produced by ionizing radiation. It was used to find out the optimal parameters for the purpose of designing an IC capable of exhibiting a linear response with dose (deviation less than 1%) up to 10 Gy per pulse at 4 μs pulse duration. As a proof of concept, two vented parallel plate IC prototypes have been built and tested in different ultra‐high pulse dose rate electron beams.

Results

It has been found that by reducing the distance between electrodes to a value of 0.25 mm it is possible to extend the dose rate operating range of parallel plate ICs to ultra‐high dose per pulse range, at standard voltage of clinical grade electrometers, well into several Gy per pulse. The two IC prototypes exhibit behavior as predicted by the numerical simulation. One of the so‐called ultra‐thin parallel plate ionization chamber (UTIC) prototypes was able to measure up to 10 Gy per pulse, 4 μs pulse duration, operated at 300 V with no significant deviation from linearity within the uncertainties (ElectronFlash Linac, SIT). The other prototype was tested up to 5.4 Gy per pulse, 2.5 μs pulse duration, operated at 250 V with CCE higher than 98.6% (Metrological Electron Accelerator Facility, MELAF at Physikalisch‐Technische Bundesanstalt, PTB).

Conclusions

This work demonstrates the ability to extend the dose rate operating range of ICs to ultra‐high dose per pulse range by reducing the spacing between electrodes. The results show that UTICs are suitable for measurement in UHDR electron beams.

Design, realization and characterization of a novel diamond detector prototype for flash radiotherapy dosimetry

Purpose

FLASH radiotherapy (RT) is an emerging technique in which beams with ultra‐high dose rates (UH‐DR) and dose per pulse (UH‐DPP) are used. Commercially available active real‐time dosimeters have been shown to be unsuitable in such conditions, due to severe response nonlinearities. In the present study, a novel diamond‐based Schottky diode detector was specifically designed and realized to match the stringent requirements of FLASH‐RT.

Methods

A systematic investigation of the main features affecting the diamond response in UH‐DPP conditions was carried out. Several diamond Schottky diode detector prototypes with different layouts were produced at Rome Tor Vergata University in cooperation with PTW‐Freiburg. Such devices were tested under electron UH‐DPP beams. The linearity of the prototypes was investigated up to DPPs of about 26 Gy/pulse and dose rates of approximately 1 kGy/s. In addition, percentage depth dose (PDD) measurements were performed in different irradiation conditions. Radiochromic films were used for reference dosimetry.

Results

The response linearity of the diamond prototypes was shown to be strongly affected by the size of their active volume as well as by their series resistance. By properly tuning the design layout, the detector response was found to be linear up to at least 20 Gy/pulse, well into the UH‐DPP range conditions. PDD measurements were performed by three different linac applicators, characterized by DPP values at the point of maximum dose of 3.5, 17.2, and 20.6 Gy/pulse, respectively. The very good superimposition of three curves confirmed the diamond response linearity. It is worth mentioning that UH‐DPP irradiation conditions may lead to instantaneous detector currents as high as several mA, thus possibly exceeding the electrometer specifications. This issue was properly addressed in the case of the PTW UNIDOS electrometers.

Conclusions

The results of the present study clearly demonstrate the feasibility of a diamond detector for FLASH‐RT applications.

Resting Brain Fluctuations Are Intrinsically Coupled to Visual Response Dynamics

How do intrinsic brain dynamics interact with processing of external sensory stimuli? We sought new insights using functional magnetic resonance imaging to track spatiotemporal activity patterns at the whole brain level in lightly anesthetized mice, during both resting conditions and visual stimulation trials. Our results provide evidence that quasiperiodic patterns (QPPs) are the most prominent component of mouse resting brain dynamics. These QPPs captured the temporal alignment of anticorrelation between the default mode (DMN)- and task-positive (TPN)-like networks, with global brain fluctuations, and activity in neuromodulatory nuclei of the reticular formation. Specifically, the phase of QPPs prior to stimulation could significantly stratify subsequent visual response magnitude, suggesting QPPs relate to brain state fluctuations. This is the first observation in mice that dynamics of the DMN- and TPN-like networks, and particularly their anticorrelation, capture a brain state dynamic that affects sensory processing. Interestingly, QPPs also displayed transient onset response properties during visual stimulation, which covaried with deactivations in the reticular formation. We conclude that QPPs appear to capture a brain state fluctuation that may be orchestrated through neuromodulation. Our findings provide new frontiers to understand the neural processes that shape functional brain states and modulate sensory input processing.

RealDRR - Rendering of realistic digitally reconstructed radiographs using locally trained image-to-image translation

INTRODUCTION

Digitally reconstructed radiographs (DRRs) represent valuable patient-specific pre-treatment training data for tumor tracking algorithms. However, using current rendering methods, the similarity of the DRRs to real X-ray images is limited, requires time-consuming measurements and/or are computationally expensive. In this study we present RealDRR, a novel framework for highly realistic and computationally efficient DRR rendering.

MATERIALS AND METHODS

RealDRR consists of two components applied sequentially to render a DRR. First, a raytracer is applied for forward projection from 3D CT data to a 2D image. Second, a conditional Generative Adverserial Network (cGAN) is applied to translate the 2D forward projection to a realistic 2D DRR. The planning CT and CBCT projections from a CIRS thorax phantom and 6 radiotherapy patients (3 prostate, 3 brain) were split in training and test sets for evaluating the intra-patient, inter-patient and inter-anatomical region generalization performance of the trained framework. Several image similarity metrics, as well as a verification based on template matching, were used between the rendered DRRs and respective CBCT projections in the test sets, and results were compared to those of a current state-of-the-art DRR rendering method.

RESULTS

When trained on 800 CBCT projection images from two patients and tested on a third unseen patient from either anatomical region, RealDRR outperformed the current state-of-the-art with statistical significance on all metrics (two-sample t-test, p < 0.05). Once trained, the framework is able to render 100 highly realistic DRRs in under two minutes.

CONCLUSION

A novel framework for realistic and efficient DRR rendering was proposed. As the framework requires a reasonable amount of computational resources, the internal parameters can be tailored to imaging systems and protocols through on-site training on retrospective imaging data.

Clinical and immunological control of experimental autoimmune encephalomyelitis by tolerogenic dendritic cells loaded with MOG-encoding mRNA

Background

Although effective in reducing relapse rate and delaying progression, current therapies for multiple sclerosis (MS) do not completely halt disease progression. T cell autoimmunity to myelin antigens is considered one of the main mechanisms driving MS. It is characterized by autoreactivity to disease-initiating myelin antigen epitope(s), followed by a cascade of epitope spreading, which are both strongly patient-dependent. Targeting a variety of MS-associated antigens by myelin antigen-presenting tolerogenic dendritic cells (tolDC) is a promising treatment strategy to re-establish tolerance in MS. Electroporation with mRNA encoding myelin proteins is an innovative technique to load tolDC with the full spectrum of naturally processed myelin-derived epitopes.

Methods

In this study, we generated murine tolDC presenting myelin oligodendrocyte glycoprotein (MOG) using mRNA electroporation and we assessed the efficacy of MOG mRNA-electroporated tolDC to dampen pathogenic T cell responses in experimental autoimmune encephalomyelitis (EAE). For this, MOG_35–55-immunized C57BL/6 mice were injected intravenously at days 13, 17, and 21 post-disease induction with 1α,25-dihydroxyvitamin D₃-treated tolDC electroporated with MOG-encoding mRNA. Mice were scored daily for signs of paralysis. At day 25, myelin reactivity was evaluated following restimulation of splenocytes with myelin-derived epitopes. Ex vivo magnetic resonance imaging (MRI) was performed to assess spinal cord inflammatory lesion load.

Results

Treatment of MOG_35–55-immunized C57BL/6 mice with MOG mRNA-electroporated or MOG_35–55-pulsed tolDC led to a stabilization of the EAE clinical score from the first administration onwards, whereas it worsened in mice treated with non-antigen-loaded tolDC or with vehicle only. In addition, MOG_35–55-specific pro-inflammatory pathogenic T cell responses and myelin antigen epitope spreading were inhibited in the peripheral immune system of tolDC-treated mice. Finally, magnetic resonance imaging analysis of hyperintense spots along the spinal cord was in line with the clinical score.

Conclusions

Electroporation with mRNA is an efficient and versatile tool to generate myelin-presenting tolDC that are capable to stabilize the clinical score in EAE. These results pave the way for further research into mRNA-electroporated tolDC treatment as a patient-tailored therapy for MS.

Electronic supplementary material

The online version of this article (10.1186/s12974-019-1541-1) contains supplementary material, which is available to authorized users.

Bottom-up sensory processing can induce negative BOLD responses and reduce functional connectivity in nodes of the default mode-like network in rats

The default mode network is a large-scale brain network that is active during rest and internally focused states and deactivates as well as desynchronizes during externally oriented (top-down) attention demanding cognitive tasks. However, it is not sufficiently understood if salient stimuli, able to trigger bottom-up attentional processes, could also result in similar reduction of activity and functional connectivity in the DMN. In this study, we investigated whether bottom-up sensory processing could influence the default mode-like network (DMLN) in rats. DMLN activity was examined using block-design visual functional magnetic resonance imaging (fMRI) while its synchronization was investigated by comparing functional connectivity during a resting versus a continuously stimulated brain state by unpredicted light flashes. We demonstrated that the BOLD response in DMLN regions was decreased during visual stimulus blocks and increased during blanks. Furthermore, decreased inter-network functional connectivity between the DMLN and visual networks as well as decreased intra-network functional connectivity within the DMLN was observed during the continuous visual stimulation. These results suggest that triggering of bottom-up attention mechanisms in sedated rats can lead to a cascade similar to top-down orienting of attention in humans and is able to deactivate and desynchronize the DMLN.

Dynamic resting state fMRI analysis in mice reveals a set of Quasi-Periodic Patterns and illustrates their relationship with the global signal

Time-resolved 'dynamic' over whole-period 'static' analysis of low frequency (LF) blood-oxygen level dependent (BOLD) fluctuations provides many additional insights into the macroscale organization and dynamics of neural activity. Although there has been considerable advancement in the development of mouse resting state fMRI (rsfMRI), very little remains known about its dynamic repertoire. Here, we report for the first time the detection of a set of recurring spatiotemporal Quasi-Periodic Patterns (QPPs) in mice, which show spatial similarity with known resting state networks. Furthermore, we establish a close relationship between several of these patterns and the global signal. We acquired high temporal rsfMRI scans under conditions of low (LA) and high (HA) medetomidine-isoflurane anesthesia. We then employed the algorithm developed by Majeed et al. (2011), previously applied in rats and humans, which detects and averages recurring spatiotemporal patterns in the LF BOLD signal. One type of observed patterns in mice was highly similar to those originally observed in rats, displaying propagation from lateral to medial cortical regions, which suggestively pertain to a mouse Task-Positive like network (TPN) and Default Mode like network (DMN). Other QPPs showed more widespread or striatal involvement and were no longer detected after global signal regression (GSR). This was further supported by diminished detection of subcortical dynamics after GSR, with cortical dynamics predominating. Observed QPPs were both qualitatively and quantitatively determined to be consistent across both anesthesia conditions, with GSR producing the same outcome. Under LA, QPPs were consistently detected at both group and single subject level. Under HA, consistency and pattern occurrence rate decreased, whilst cortical contribution to the patterns diminished. These findings confirm the robustness of QPPs across species and demonstrate a new approach to study mouse LF BOLD spatiotemporal dynamics and mechanisms underlying functional connectivity. The observed impact of GSR on QPPs might help better comprehend its controversial role in conventional resting state studies. Finally, consistent detection of QPPs at single subject level under LA promises a step forward towards more reliable mouse rsfMRI and further confirms the importance of selecting an optimal anesthesia regime.