Energetic constraints imposed on trophic interaction strengths enhance resilience in empirical and model food webs.

Food web stability and resilience are at the heart of understanding the structure and functioning of ecosystems. Previous studies show that models of empirical food webs are substantially more stable than random ones, due to a few strong interactions embedded in a majority of weak interactions. Analyses of trophic interaction loops show that in empirical food webs the patterns of the interaction strengths prevent the occurrence of destabilizing heavy loops and thereby enhances resilience. Yet, it is still unexplored which biological mechanisms cause these patterns that enhance food web resilience. We quantified food web resilience using the real part of the maximum eigenvalue of the Jacobian matrix of the food web from a seagrass bed in the Yellow River Delta (YRD) wetland, that could be parametrized by the empirical data of the food web. We found that the empirically based Jacobian matrix of the YRD food web indicated a much higher resilience than random matrices with the same element values but arranged in random ways. Investigating the trophic interaction loops revealed that the high resilience was due to a negative correlation between the negative and positive interaction strengths (per capita top-down and bottom-up effects, respectively) within positive feedback loops with three species. The negative correlation showed that when the negative interaction strengths were strong the positive was weak, and vice versa. Our invented reformulation of loop weight in terms of biomasses and specific production rates showed that energetic properties of the trophic groups in the loop and mass-balance constraints, for example, the food uptake has to balance all losses, created the negative correlation between the interaction strengths. This result could be generalized using a dynamic intraguild predation model, which delivered the same pattern for a wide range of model parameters. Our results shed light on how energetic constraints at the trophic group and food web level create a pattern of interaction strengths within trophic interaction loops that enhances food web resilience.

First system for fully-automated multi-criterial treatment planning for a high-magnetic field MR-Linac applied to rectal cancer.

Background and purpose: In this study we developed a workflow for fully-automated generation of deliverable IMRT plans for a 1.5 T MR-Linac (MRL) based on contoured CT scans, and we evaluated automated MRL planning for rectal cancer.Methods: The Monte Carlo dose calculation engine used in the clinical MRL TPS (Monaco, Elekta AB, Stockholm, Sweden), suited for high accuracy dose calculations in a 1.5 T magnetic field, was coupled to our in-house developed Erasmus-iCycle optimizer. Clinically deliverable plans for 23 rectal cancer patients were automatically generated in a two-step process, i.e., multi-criterial fluence map optimization with Erasmus-iCycle followed by a conversion into a deliverable IMRT plan in the clinical TPS. Automatically generated plans (AUTOplans) were compared to plans that were manually generated with the clinical TPS (MANplans).Results: With AUTOplanning large reductions in planning time and workload were obtained; 4-6 h mainly hands-on planning for MANplans vs ∼1 h of mainly computer computation time for AUTOplans. For equal target coverage, the bladder and bowel bag Dmean was reduced in the AUTOplans by 1.3 Gy (6.9%) on average with a maximum reduction of 4.5 Gy (23.8%). Dosimetric measurements at the MRL demonstrated clinically acceptable delivery accuracy for the AUTOplans.Conclusions: A system for fully automated multi-criterial planning for a 1.5 T MR-Linac was developed and tested for rectal cancer patients. Automated planning resulted in major reductions in planning workload and time, while plan quality improved. Negative impact of the high magnetic field on the dose distributions could be avoided.

Low investment in sexual reproduction threatens plants adapted to phosphorus limitation.

Plant species diversity in Eurasian wetlands and grasslands depends not only on productivity but also on the relative availability of nutrients, particularly of nitrogen and phosphorus. Here we show that the impacts of nitrogen:phosphorus stoichiometry on plant species richness can be explained by selected plant life-history traits, notably by plant investments in growth versus reproduction. In 599 Eurasian sites with herbaceous vegetation we examined the relationship between the local nutrient conditions and community-mean life-history traits. We found that compared with plants in nitrogen-limited communities, plants in phosphorus-limited communities invest little in sexual reproduction (for example, less investment in seed, shorter flowering period, longer lifespan) and have conservative leaf economy traits (that is, a low specific leaf area and a high leaf dry-matter content). Endangered species were more frequent in phosphorus-limited ecosystems and they too invested little in sexual reproduction. The results provide new insight into how plant adaptations to nutrient conditions can drive the distribution of plant species in natural ecosystems and can account for the vulnerability of endangered species.

Priorities for research in soil ecology.

The ecological interactions that occur in and with soil are of consequence in many ecosystems on the planet. These interactions provide numerous essential ecosystem services, and the sustainable management of soils has attracted increasing scientific and public attention. Although soil ecology emerged as an independent field of research many decades ago, and we have gained important insights into the functioning of soils, there still are fundamental aspects that need to be better understood to ensure that the ecosystem services that soils provide are not lost and that soils can be used in a sustainable way. In this perspectives paper, we highlight some of the major knowledge gaps that should be prioritized in soil ecological research. These research priorities were compiled based on an online survey of 32 editors of Pedobiologia - Journal of Soil Ecology. These editors work at universities and research centers in Europe, North America, Asia, and Australia.The questions were categorized into four themes: (1) soil biodiversity and biogeography, (2) interactions and the functioning of ecosystems, (3) global change and soil management, and (4) new directions. The respondents identified priorities that may be achievable in the near future, as well as several that are currently achievable but remain open. While some of the identified barriers to progress were technological in nature, many respondents cited a need for substantial leadership and goodwill among members of the soil ecology research community, including the need for multi-institutional partnerships, and had substantial concerns regarding the loss of taxonomic expertise.

How habitat-modifying organisms structure the food web of two coastal ecosystems.

The diversity and structure of ecosystems has been found to depend both on trophic interactions in food webs and on other species interactions such as habitat modification and mutualism that form non-trophic interaction networks. However, quantification of the dependencies between these two main interaction networks has remained elusive. In this study, we assessed how habitat-modifying organisms affect basic food web properties by conducting in-depth empirical investigations of two ecosystems: North American temperate fringing marshes and West African tropical seagrass meadows. Results reveal that habitat-modifying species, through non-trophic facilitation rather than their trophic role, enhance species richness across multiple trophic levels, increase the number of interactions per species (link density), but decrease the realized fraction of all possible links within the food web (connectance). Compared to the trophic role of the most highly connected species, we found this non-trophic effects to be more important for species richness and of more or similar importance for link density and connectance. Our findings demonstrate that food webs can be fundamentally shaped by interactions outside the trophic network, yet intrinsic to the species participating in it. Better integration of non-trophic interactions in food web analyses may therefore strongly contribute to their explanatory and predictive capacity.

Soil food web properties explain ecosystem services across European land use systems.

Intensive land use reduces the diversity and abundance of many soil biota, with consequences for the processes that they govern and the ecosystem services that these processes underpin. Relationships between soil biota and ecosystem processes have mostly been found in laboratory experiments and rarely are found in the field. Here, we quantified, across four countries of contrasting climatic and soil conditions in Europe, how differences in soil food web composition resulting from land use systems (intensive wheat rotation, extensive rotation, and permanent grassland) influence the functioning of soils and the ecosystem services that they deliver. Intensive wheat rotation consistently reduced the biomass of all components of the soil food web across all countries. Soil food web properties strongly and consistently predicted processes of C and N cycling across land use systems and geographic locations, and they were a better predictor of these processes than land use. Processes of carbon loss increased with soil food web properties that correlated with soil C content, such as earthworm biomass and fungal/bacterial energy channel ratio, and were greatest in permanent grassland. In contrast, processes of N cycling were explained by soil food web properties independent of land use, such as arbuscular mycorrhizal fungi and bacterial channel biomass. Our quantification of the contribution of soil organisms to processes of C and N cycling across land use systems and geographic locations shows that soil biota need to be included in C and N cycling models and highlights the need to map and conserve soil biodiversity across the world.

Intensive agriculture reduces soil biodiversity across Europe.

Soil biodiversity plays a key role in regulating the processes that underpin the delivery of ecosystem goods and services in terrestrial ecosystems. Agricultural intensification is known to change the diversity of individual groups of soil biota, but less is known about how intensification affects biodiversity of the soil food web as a whole, and whether or not these effects may be generalized across regions. We examined biodiversity in soil food webs from grasslands, extensive, and intensive rotations in four agricultural regions across Europe: in Sweden, the UK, the Czech Republic and Greece. Effects of land-use intensity were quantified based on structure and diversity among functional groups in the soil food web, as well as on community-weighted mean body mass of soil fauna. We also elucidate land-use intensity effects on diversity of taxonomic units within taxonomic groups of soil fauna. We found that between regions soil food web diversity measures were variable, but that increasing land-use intensity caused highly consistent responses. In particular, land-use intensification reduced the complexity in the soil food webs, as well as the community-weighted mean body mass of soil fauna. In all regions across Europe, species richness of earthworms, Collembolans, and oribatid mites was negatively affected by increased land-use intensity. The taxonomic distinctness, which is a measure of taxonomic relatedness of species in a community that is independent of species richness, was also reduced by land-use intensification. We conclude that intensive agriculture reduces soil biodiversity, making soil food webs less diverse and composed of smaller bodied organisms. Land-use intensification results in fewer functional groups of soil biota with fewer and taxonomically more closely related species. We discuss how these changes in soil biodiversity due to land-use intensification may threaten the functioning of soil in agricultural production systems.

Infectious disease agents mediate interaction in food webs and ecosystems.

Infectious agents are part of food webs and ecosystems via the relationship with their host species that, in turn, interact with both hosts and non-hosts. Through these interactions, infectious agents influence food webs in terms of structure, functioning and stability. The present literature shows a broad range of impacts of infectious agents on food webs, and by cataloguing that range, we worked towards defining the various mechanisms and their specific effects. To explore the impact, a direct approach is to study changes in food-web properties with infectious agents as separate species in the web, acting as additional nodes, with links to their host species. An indirect approach concentrates not on adding new nodes and links, but on the ways that infectious agents affect the existing links across host and non-host nodes, by influencing the 'quality' of consumer-resource interaction as it depends on the epidemiological state host involved. Both approaches are natural from an ecological point of view, but the indirect approach may connect more straightforwardly to commonly used tools in infectious disease dynamics.

Reconciling complexity with stability in naturally assembling food webs.

Understanding how complex food webs assemble through time is fundamental both for ecological theory and for the development of sustainable strategies of ecosystem conservation and restoration. The build-up of complexity in communities is theoretically difficult, because in random-pattern models complexity leads to instability. There is growing evidence, however, that nonrandom patterns in the strengths of the interactions between predators and prey strongly enhance system stability. Here we show how such patterns explain stability in naturally assembling communities. We present two series of below-ground food webs along natural productivity gradients in vegetation successions. The complexity of the food webs increased along the gradients. The stability of the food webs was captured by measuring the weight of feedback loops of three interacting 'species' locked in omnivory. Low predator-prey biomass ratios in these omnivorous loops were shown to have a crucial role in preserving stability as productivity and complexity increased during succession. Our results show the build-up of food-web complexity in natural productivity gradients and pin down the feedback loops that govern the stability of whole webs. They show that it is the heaviest three-link feedback loop in a network of predator-prey effects that limits its stability. Because the weight of these feedback loops is kept relatively low by the biomass build-up in the successional process, complexity does not lead to instability.

Spatial vegetation patterns and imminent desertification in Mediterranean arid ecosystems.

Humans and climate affect ecosystems and their services, which may involve continuous and discontinuous transitions from one stable state to another. Discontinuous transitions are abrupt, irreversible and among the most catastrophic changes of ecosystems identified. For terrestrial ecosystems, it has been hypothesized that vegetation patchiness could be used as a signature of imminent transitions. Here, we analyse how vegetation patchiness changes in arid ecosystems with different grazing pressures, using both field data and a modelling approach. In the modelling approach, we extrapolated our analysis to even higher grazing pressures to investigate the vegetation patchiness when desertification is imminent. In three arid Mediterranean ecosystems in Spain, Greece and Morocco, we found that the patch-size distribution of the vegetation follows a power law. Using a stochastic cellular automaton model, we show that local positive interactions among plants can explain such power-law distributions. Furthermore, with increasing grazing pressure, the field data revealed consistent deviations from power laws. Increased grazing pressure leads to similar deviations in the model. When grazing was further increased in the model, we found that these deviations always and only occurred close to transition to desert, independent of the type of transition, and regardless of the vegetation cover. Therefore, we propose that patch-size distributions may be a warning signal for the onset of desertification.

Benchmarking daily adaptation using fully automated radiotherapy treatment plan optimization for rectal cancer.

Background/purpose: In daily plan adaptation the radiotherapy treatment plan is adjusted just prior to delivery. A simple approach is taking the planning objectives of the reference plan and directly applying these in re-optimization. Here we present a tested method to verify whether daily adaptation without tweaking of the objectives can maintain the plan quality throughout treatment. Materials/methods: For fifteen rectal cancer patients, automated treatment planning was used to generate plans mimicking manual reference plans on the planning scans. For 74 fraction scans (4-5 per patient) an automated plan and a daily adapted plan were generated, where the latter re-optimizes the reference plan objectives without any tweaking. To evaluate the robustness of the daily adaptation, the adapted plans were compared to the autoplanning plans. Results: Median differences between the autoplanning plans on the planning scans and the reference plans were between -1 and 0.2 Gy. The largest interquartile range (1 Gy) was seen for the Lumbar Skin D2%. For the daily scans the PTV D2% and D98% differences between autoplanning and adapted plans were within ± 0.7 Gy, with mean differences within ± 0.3 Gy. Positive differences indicate higher values were obtained using autoplanning. For the Bowelarea + Bladder and the Lumbar Skin the D2% and Dmean differences were all within ± 2.6 Gy, with mean differences between -0.9 and 0.1 Gy. Conclusion: Automated treatment planning can be used to benchmark daily adaptation techniques. The investigated adaptation workflow can robustly perform high quality adaptations without daily adjusting of the patient-specific planning objectives for rectal cancer radiotherapy.

Robust scaling in ecosystems and the meltdown of patch size distributions before extinction.

Robust critical systems are characterized by power laws which occur over a broad range of conditions. Their robust behaviour has been explained by local interactions. While such systems could be widespread in nature, their properties are not well understood. Here, we study three robust critical ecosystem models and a null model that lacks spatial interactions. In all these models, individuals aggregate in patches whose size distributions follow power laws which melt down under increasing external stress. We propose that this power-law decay associated with the connectivity of the system can be used to evaluate the level of stress exerted on the ecosystem. We identify several indicators along the transition to extinction. These indicators give us a relative measure of the distance to extinction, and have therefore potential application to conservation biology, especially for ecosystems with self-organization and critical transitions.

MR-Linac Radiotherapy - The Beam Angle Selection Problem.

BACKGROUND: With the large-scale introduction of volumetric modulated arc therapy (VMAT), selection of optimal beam angles for coplanar static-beam IMRT has increasingly become obsolete. Due to unavailability of VMAT in current MR-linacs, the problem has re-gained importance. An application for automated IMRT treatment planning with integrated, patient-specific computer-optimization of beam angles (BAO) was used to systematically investigate computer-aided generation of beam angle class solutions (CS) for replacement of computationally expensive patient-specific BAO. Rectal cancer was used as a model case. MATERIALS AND METHODS: 23 patients treated at a Unity MR-linac were included. BAOx plans (x=7-12 beams) were generated for all patients. Analyses of BAO12 plans resulted in CSx class solutions. BAOx plans, CSx plans, and plans with equi-angular setups (EQUIx, x=9-56) were mutually compared. RESULTS: For x>7, plan quality for CSx and BAOx was highly similar, while both were superior to EQUIx. E.g. with CS9, bowel/bladder Dmean reduced by 22% [11%, 38%] compared to EQUI9 (p<0.001). For equal plan quality, the number of EQUI beams had to be doubled compared to BAO and CS. CONCLUSIONS: Computer-generated beam angle CS could replace individualized BAO without loss in plan quality, while reducing planning complexity and calculation times, and resulting in a simpler clinical workflow. CS and BAO largely outperformed equi-angular treatment. With the developed CS, time consuming beam angle re-optimization in daily adaptive MR-linac treatment could be avoided. Further systematic research on computerized development of beam angle class solutions for MR-linac treatment planning is warranted.

Resource contrast in patterned peatlands increases along a climatic gradient.

Spatial patterning of ecosystems can be explained by several mechanisms. One approach to disentangling the influence of these mechanisms is to study a patterned ecosystem along a gradient of environmental conditions. This study focused on hummock-hollow patterning of peatlands. Previous models predicted that patterning in drainage-dominated peatlands is driven by a peat-accumulation mechanism, reflected by higher nutrient availability in hollows relative to hummocks. Alternatively, patterning in evapotranspiration (ET)-dominated peatlands may be driven by a nutrient-accumulation mechanism, reflected by reversed nutrient distribution, namely, higher nutrient availability in hummocks relative to hollows. Here, we tested these predictions by comparing nutrient distributions among patterned peatlands in maritime (Scotland), humid temperate (Sweden), and humid continental (Siberia) climates. The areas comprise a climatic gradient from very wet and drainage-dominated (Scotland) to less wet and ET-dominated (Siberia) peatlands. Nutrient distribution was quantified as resource contrast, a measure for hummock-hollow difference in nutrient availability. We tested the hypothesis that the climatic gradient shows a trend in the resource contrast; from negative (highest nutrient availability in hollows) in Scotland to positive (highest nutrient availability in hummocks) in Siberia. The resource contrasts as measured in vegetation indeed showed a trend along the climatic gradient: contrasts were negative to slightly positive in Scotland, positive in Sweden, and strongly positive in Siberia. This finding corroborates the main prediction of previous models. Our results, however, also provided indications for further model development. The low concentrations of nutrients in the water suggest that existing models could be improved by considering both the dissolved and adsorbed phase and explicit inclusion of both nutrient-uptake and nutrient-storage processes. Our study suggests that future climate change may affect the ecosystem functioning of patterned peatlands by altering the contribution of pattern-forming mechanisms to redistribution of water and nutrients within these systems.

Knowledge-Based Assessment of Focal Dose Escalation Treatment Plans in Prostate Cancer.

PURPOSE: In a randomized focal dose escalation radiation therapy trial for prostate cancer (FLAME), up to 95 Gy was prescribed to the tumor in the dose-escalated arm, with 77 Gy to the entire prostate in both arms. As dose constraints to organs at risk had priority over dose escalation and suboptimal planning could occur, we investigated how well the dose to the tumor was boosted. We developed an anatomy-based prediction model to identify plans with suboptimal tumor dose and performed replanning to validate our model. METHODS AND MATERIALS: We derived dose-volume parameters from planned dose distributions of 539 FLAME trial patients in 4 institutions and compared them between both arms. In the dose-escalated arm, we determined overlap volume histograms and derived features representing patient anatomy. We predicted tumor D98% with a linear regression on anatomic features and performed replanning on 21 plans. RESULTS: In the dose-escalated arm, the median tumor D50% and D98% were 93.0 and 84.7 Gy, and 99% of the tumors had a dose escalation greater than 82.4 Gy (107% of 77 Gy). In both arms organs at risk constraints were met. Five out of 73 anatomic features were found to be predictive for tumor D98%. Median predicted tumor D98% was 4.4 Gy higher than planned D98%. Upon replanning, median tumor D98% increased by 3.0 Gy. A strong correlation between predicted increase in D98% and realized increase upon replanning was found (ρ = 0.86). CONCLUSIONS: Focal dose escalation in prostate cancer was feasible with a dose escalation to 99% of the tumors. Replanning resulted in an increased tumor dose that correlated well with the prediction model. The model was able to identify tumors on which a higher boost dose could be planned. The model has potential as a quality assessment tool in focal dose escalated treatment plans.

Nutrients and hydrology indicate the driving mechanisms of peatland surface patterning.

Peatland surface patterning motivates studies that identify underlying structuring mechanisms. Theoretical studies so far suggest that different mechanisms may drive similar types of patterning. The long time span associated with peatland surface pattern formation, however, limits possibilities for empirically testing model predictions by field manipulations. Here, we present a model that describes spatial interactions between vegetation, nutrients, hydrology, and peat. We used this model to study pattern formation as driven by three different mechanisms: peat accumulation, water ponding, and nutrient accumulation. By on-and-off switching of each mechanism, we created a full-factorial design to see how these mechanisms affected surface patterning (pattern of vegetation and peat height) and underlying patterns in nutrients and hydrology. Results revealed that different combinations of structuring mechanisms lead to similar types of peatland surface patterning but contrasting underlying patterns in nutrients and hydrology. These contrasting underlying patterns suggest that the presence or absence of the structuring mechanisms can be identified by relatively simple short-term field measurements of nutrients and hydrology, meaning that longer-term field manipulations can be circumvented. Therefore, this study provides promising avenues for future empirical studies on peatland patterning.

Complexity in quantitative food webs.

Food webs depict who eats whom in communities. Ecologists have examined statistical metrics and other properties of food webs, but mainly due to the uneven quality of the data, the results have proved controversial. The qualitative data on which those efforts rested treat trophic interactions as present or absent and disregard potentially huge variation in their magnitude, an approach similar to analyzing traffic without differentiating between highways and side roads. More appropriate data are now available and were used here to analyze the relationship between trophic complexity and diversity in 59 quantitative food webs from seven studies (14-202 species) based on recently developed quantitative descriptors. Our results shed new light on food-web structure. First, webs are much simpler when considered quantitatively, and link density exhibits scale invariance or weak dependence on food-web size. Second, the "constant connectance" hypothesis is not supported: connectance decreases with web size in both qualitative and quantitative data. Complexity has occupied a central role in the discussion of food-web stability, and we explore the implications for this debate. Our findings indicate that larger webs are more richly endowed with the weak trophic interactions that recent theories show to be responsible for food-web stability.

Investigating the impact of patient arm position in an MR-linac on liver SBRT treatment plans.

PURPOSE: The superior soft-tissue contrast offered by the integrated magnetic resonance imaging (MRI) of the Unity MR-linac compared to the x-ray-based image guidance on conventional linacs potentially allows for liver stereotactic body radiation therapy (SBRT) without the need for implanted markers or other surrogates. On conventional linacs, liver SBRT patients are typically positioned with their arms above their heads (arms-up) to minimize exposure to healthy tissue. However, the spatial confinement of the MRI-bore and increased treatment delivery times can make the arms-up position straining. Therefore, we assessed the plan quality for MR-linac treatment plans with the patient in the arms-up and in the arms-down position. Additionally, we compared the MR-linac plans with clinically used arms-up treatment plans made for a conventional linac. METHODS: Fifteen consecutively treated patients with oligometastatic liver disease were included in this retrospective study. For each patient, a planning computed tomography (CT) with delineations, a diagnostic MRI, and a 3 × 20 Gy dual-arc volumetric modulated arc therapy (VMAT) plan, which was used to treat the patient in an arms-up position on the conventional linac, were available. For the MR-linac, 15-beam step-and-shoot intensity-modulated radiation therapy (IMRT) plans were created for four patient positioning scenarios: arms-up, mimicking current clinical practice; arms-down, with treatment beams avoiding the arms on the entrance side; arms-through, arms are down but not avoided, and right-arm-up; only the right arm is up and the left arm is avoided on the entrance side. Resulting treatment plans were compared. Bonferroni-corrected two-sided Wilcoxon signed-ranks tests were used to assess statistical significance (P < 0.05). RESULTS: No significant differences were found in gross tumour volume (GTV) coverage (D 2 % , D 50 % , and D 98 % ) or liver sparing (liver-GTV V < 15 Gy ) between the clinical plans and any of the MR-linac plans. The median target conformity [exterior V 40 % /planning target volume (PTV)] was significantly better in the clinical plans (5.8) than in the MR-linac scenarios (arms-down: 6.6, arms-up/right-arm-up: 6.2, arms-through: 6.3). No MR-linac plan violated any additional organ-at-risk (OAR) constraint that was not already violated in the clinical plans. In the arms-down scenario a significantly increased median spinal cord D 1 % (14.5 Gy) was detected compared to the clinical setup (7.2 Gy). For the arms-down (arms-through) scenario, the median left arm D 1 % was 1.5 (2.7) Gy, the median right arm D 1 % was 5.8 (22.7) Gy, and the median right arm V 20 Gy was 0.0 (14.7) cc. These differences were statistically significant. For the right-arm-up scenario, the median left arm D 1 % (2.3 Gy) and V 5 Gy (0.0) were not significantly different compared to the arms-down scenario. CONCLUSIONS: Mimicking the current clinical practice by treating patients in the arms-up/right-arm-up position on the MR-linac leads to plans which are dosimetrically very similar to the conventional linac plans. Treating in the arms-down position is expected to increase patient comfort at the cost of a small reduction in OAR sparing for individual patients. Treating through the arms is not encouraged due to substantial dose deposition in the arms.

Evaluation of plan quality in radiotherapy planning with an MR-linac.

BACKGROUND & PURPOSE: Clinical introduction of magnetic resonance (MR)-guided radiotherapy involves treatment planning while taking into account machine-specific characteristics. Our aim was to investigate the feasibility of high-quality MR-linac treatment planning for an MR-linac and to benchmark MR-linac plan quality (IMRT) against current clinical practice (VMAT). MATERIALS & METHODS: Data of eight rectal and eight prostate cancer patients, who received radiotherapy on a conventional CBCT-integrated linac, were selected. Clinically acquired CTs and associated delineations of target volumes and organs-at-risk (OARs) were used for MR-linac treatment planning in Monaco. To investigate treatment planning software bias 'quasi MR-linac plans' were generated in Pinnacle3 by mimicking MR-linac specific beam characteristics. MR-linac, quasi MR-linac, and clinical plans were compared and differences in target and OAR doses assessed. Differences in plan complexity were determined by the number of segments and monitor units. RESULTS: Compared to clinical plans, MR-linac plans showed a statistically significant decrease in plan homogeneity, an increase in PTV Dmean (prostate: 0.6 Gy; rectum: 0.8 Gy) and D1% (prostate: 1.9 Gy; rectum: 2.0 Gy), and increases in OAR dose. Quasi MR-linac plans were comparable to MR-linac plans with respect to OAR dose and plan homogeneity. For rectal cancer an increase was seen in PTV Dmean (0.12 Gy) and D1% (0.5 Gy) compared to regular MR-linac plans. All created plans were clinically equivalent to current clinical practice. CONCLUSIONS: This study demonstrates the feasibility of creating high-quality MR-linac treatment plans. The results supported the clinical introduction of an MR-linac.