Local and global error models to improve uncertainty quantification

In groundwater applications, Monte Carlo methods are employed to model the uncertainty on geological parameters. However, their brute-force application becomes computationally prohibitive for highly detailed geological descriptions, complex physical processes, and a large number of realizations. The Distance Kernel Method (DKM) overcomes this issue by clustering the realizations in a multidimensional space based on the flow responses obtained by means of an approximate (computationally cheaper) model; then, the uncertainty is estimated from the exact responses that are computed only for one representative realization per cluster (the medoid). Usually, DKM is employed to decrease the size of the sample of realizations that are considered to estimate the uncertainty. We propose to use the information from the approximate responses for uncertainty quantification. The subset of exact solutions provided by DKM is then employed to construct an error model and correct the potential bias of the approximate model. Two error models are devised that both employ the difference between approximate and exact medoid solutions, but differ in the way medoid errors are interpolated to correct the whole set of realizations. The Local Error Model rests upon the clustering defined by DKM and can be seen as a natural way to account for intra-cluster variability; the Global Error Model employs a linear interpolation of all medoid errors regardless of the cluster to which the single realization belongs. These error models are evaluated for an idealized pollution problem in which the uncertainty of the breakthrough curve needs to be estimated. For this numerical test case, we demonstrate that the error models improve the uncertainty quantification provided by the DKM algorithm and are effective in correcting the bias of the estimate computed solely from the MsFV results. The framework presented here is not specific to the methods considered and can be applied to other combinations of approximate models and techniques to select a subset of realizations

Threats from climate change to terrestrial vertebrate hotspots in europe.

We identified hotspots of terrestrial vertebrate species diversity in Europe and adjacent islands. Moreover, we assessed the extent to which by the end of the 21(st) century such hotspots will be exposed to average monthly temperature and precipitation patterns which can be regarded as extreme if compared to the climate experienced during 1950-2000. In particular, we considered the entire European sub-continent plus Turkey and a total of 1149 species of terrestrial vertebrates. For each species, we developed species-specific expert-based distribution models (validated against field data) which we used to calculate species richness maps for mammals, breeding birds, amphibians, and reptiles. Considering four global circulation model outputs and three emission scenarios, we generated an index of risk of exposure to extreme climates, and we used a bivariate local Moran's I to identify the areas with a significant association between hotspots of diversity and high risk of exposure to extreme climates. Our results outline that the Mediterranean basin represents both an important hotspot for biodiversity and especially for threatened species for all taxa. In particular, the Iberian and Italian peninsulas host particularly high species richness as measured over all groups, while the eastern Mediterranean basin is particularly rich in amphibians and reptiles; the islands (both Macaronesian and Mediterranean) host the highest richness of threatened species for all taxa occurs. Our results suggest that the main hotspots of biodiversity for terrestrial vertebrates may be extensively influenced by the climate change projected to occur over the coming decades, especially in the Mediterranean bioregion, posing serious concerns for biodiversity conservation.

Spatial relationship between the atmospheric circulation and the precipitation measured in the western Swiss Alps by means of the analogue method

An adaptation technique based on the synoptic atmospheric circulation to forecast local precipitation, namely the analogue method, has been implemented for the western Swiss Alps. During the calibration procedure, relevance maps were established for the geopotential height data. These maps highlight the locations were the synoptic circulation was found of interest for the precipitation forecasting at two rain gauge stations (Binn and Les Marécottes) that are located both in the alpine Rhône catchment, at a distance of about 100 km from each other. These two stations are sensitive to different atmospheric circulations. We have observed that the most relevant data for the analogue method can be found where specific atmospheric circulation patterns appear concomitantly with heavy precipitation events. Those skilled regions are coherent with the atmospheric flows illustrated, for example, by means of the back trajectories of air masses. Indeed, the circulation recurrently diverges from the climatology during days with strong precipitation on the southern part of the alpine Rhône catchment. We have found that for over 152 days with precipitation amount above 50 mm at the Binn station, only 3 did not show a trajectory of a southerly flow, meaning that such a circulation was present for 98% of the events. Time evolution of the relevance maps confirms that the atmospheric circulation variables have significantly better forecasting skills close to the precipitation period, and that it seems pointless for the analogue method to consider circulation information days before a precipitation event as a primary predictor. Even though the occurrence of some critical circulation patterns leading to heavy precipitation events can be detected by precursors at remote locations and 1 week ahead (Grazzini, 2007; Martius et al., 2008), time extrapolation by the analogue method seems to be rather poor. This would suggest, in accordance with previous studies (Obled et al., 2002; Bontron and Obled, 2005), that time extrapolation should be done by the Global Circulation Model, which can process atmospheric variables that can be used by the adaptation method.

Homéostasie, marqueurs somatiques et intentionnalité : le statut de la représentation et ses implications en psychanalyse et en neuroscience

Résumé : Si la psychanalyse est régulièrement remise en cause pour son manque de fondement scientifique, de nombreux travaux issus des neurosciences en appellent paradoxalement à un retour aux thèses freudiennes. Cette étude se propose, dans ce contexte, de revisiter les bases neurobiologiques de la métapsychologie et de montrer comment celle-ci repose sur une conception homéostatique du vivant. On en tire plusieurs conséquences. Premièrement, contre Brentano, selon qui l'intentionnalité est la marque du psychologique, on propose que celle-ci s'avère être, plus originellement, la marque des régulations biologiques. Deuxièmement, on montre comment Freud et Damasio développent des conceptions similaires sur le fonctionnement homéostatique de la cognition. Sur le plan épistémologique, on plaide en faveur de la complémentarité de ces deux approches. Si Damasio met au jour les variables somatiques permettant de mesurer les effets de la vie psychique, la théorie pulsionnelle permet de dépasser les problèmes liés à l'ambiguïté de la notion de «représentation cérébrale ». C'est à partir de cette thèse que l'on peut, selon nous, appréhender la psychanalyse comme un double projet visant à articuler une naturalisation de l'appareil psychique à une herméneutique de ses productions. On considère cette double visée comme une réponse à la difficulté centrale que rencontrent les sciences de l'esprit quant aux conditions de possibilité d'une naturalisation de l'intentionnalité. A partir de cette problématique, il devient pertinent d'évaluer la solution freudienne à la lumière de la posture intentionnelle dennettienne et de la théorie gazzaniguienne de l'interpreter. II ne s'agit pas, selon ces dernières approches, de naturaliser directement l'intentionnalité, au risque de commettre une erreur de catégorie. Il est question, au contraire, d'utiliser les principes du fonctionnement homéostatique comme autant de «règles »interprétatives permettant, sur le plan psychologique, de construire du sens. On trouve enfin dans le modèle de l'espace global de travail un moyen de remettre au goût du jour la conception freudienne, dynamique et conflictuelle, du fonctionnement psychique. On avance enfin, à partir de ce modèle, des hypothèses sur les mécanismes cérébraux susceptibles de sous-tendre l'efficacité thérapeutique de la cure analytique. Abstract : While Psychoanalysis' scientific basis is open to question, some neuroscience works are paradoxically calling for a return to Freud. Therefore, the first aim of our study consists in revisiting metapychology's neurobiological roots and to show how it is grounded on a homeostatic theory of life. Several consequences can be drawn from this statement. Firstly, in opposition to Brentanian definition of intentionality, as the specific mark of psychology, we argue that it is more specifically the mark of biological balance. Secondly, this statement allows to show how Freud's and Damasio's theories share common views on the homeostatic functioning of cognition. From an epistemological point of view, they complete one another. If Damasio has highlighted the measurable somatic variables from which we can infer psychic life, the theory of drives allows conceptual difficulties linked to the use of the ambiguous notion of "cerebral representation" to be overcome. According to us, this thesis leads psychoanalysis to be approached as a twin project aiming to articulate psyche naturalization, to a hermeneutic of its productions. It maybe seen as a way to respond to the central issue of mind sciences, that is, to account for the naturalization of intentionality. Given this theoretical framework, it seems relevant to reconsider freudism in the light of the dennettian intentional stance and the gazzanigan theory of interpreter. Then, freudism can be seen as a way to avoid a category mistake. Its solution rejects direct intentionality naturalization in favor of a construction of sense. In this framework, interpretation is regulated by the rules abstracted from the homeostatic functioning of life. Furthermore, we show how the dynamic and conflictual Freudian psyche can be evaluated using the Global workspace model, which allows us to put forth hypotheses on the cerebral mechanisms that may underlie the efficiency of analytical cure.

Modeling resting-state functional networks when the cortex falls asleep: local and global changes.

The transition from wakefulness to sleep represents the most conspicuous change in behavior and the level of consciousness occurring in the healthy brain. It is accompanied by similarly conspicuous changes in neural dynamics, traditionally exemplified by the change from "desynchronized" electroencephalogram activity in wake to globally synchronized slow wave activity of early sleep. However, unit and local field recordings indicate that the transition is more gradual than it might appear: On one hand, local slow waves already appear during wake; on the other hand, slow sleep waves are only rarely global. Studies with functional magnetic resonance imaging also reveal changes in resting-state functional connectivity (FC) between wake and slow wave sleep. However, it remains unclear how resting-state networks may change during this transition period. Here, we employ large-scale modeling of the human cortico-cortical anatomical connectivity to evaluate changes in resting-state FC when the model "falls asleep" due to the progressive decrease in arousal-promoting neuromodulation. When cholinergic neuromodulation is parametrically decreased, local slow waves appear, while the overall organization of resting-state networks does not change. Furthermore, we show that these local slow waves are structured macroscopically in networks that resemble the resting-state networks. In contrast, when the neuromodulator decrease further to very low levels, slow waves become global and resting-state networks merge into a single undifferentiated, broadly synchronized network.

A general model to predict individual exposure to solar UV by using ambient irradiance data

Excessive exposure to solar ultraviolet (UV) is the main cause of skin cancer. Specific prevention should be further developed to target overexposed or highly vulnerable populations. A better characterisation of anatomical UV exposure patterns is however needed for specific prevention. To develop a regression model for predicting the UV exposure ratio (ER, ratio between the anatomical dose and the corresponding ground level dose) for each body site without requiring individual measurements. A 3D numeric model (SimUVEx) was used to compute ER for various body sites and postures. A multiple fractional polynomial regression analysis was performed to identify predictors of ER. The regression model used simulation data and its performance was tested on an independent data set. Two input variables were sufficient to explain ER: the cosine of the maximal daily solar zenith angle and the fraction of the sky visible from the body site. The regression model was in good agreement with the simulated data ER (R(2)=0.988). Relative errors up to +20% and -10% were found in daily doses predictions, whereas an average relative error of only 2.4% (-0.03% to 5.4%) was found in yearly dose predictions. The regression model predicts accurately ER and UV doses on the basis of readily available data such as global UV erythemal irradiance measured at ground surface stations or inferred from satellite information. It renders the development of exposure data on a wide temporal and geographical scale possible and opens broad perspectives for epidemiological studies and skin cancer prevention.

Development of glycopeptide-intermediate resistance by Staphylococcus aureus leads to attenuated infectivity in a rat model of endocarditis.

Glycopeptide-intermediate resistant Staphylococcus aureus (GISA) are characterized by multiple changes in the cell wall and an altered expression of global virulence regulators. We investigated whether GISA are affected in their infectivity in a rat model of experimental endocarditis. The glycopeptide-susceptible, methicillin-resistant S. aureus M1V2 and its laboratory-derived GISA M1V16 were examined for their ability to (i) adhere to fibrinogen and fibronectin in vitro, (ii) persist in the bloodstream after intravenous inoculation, (iii) colonize aortic vegetations in rats, and (iv) compete for valve colonization by co-inoculation. Both GISA M1V16 and M1V2 adhered similarly to fibrinogen and fibronectin in vitro. In rats, GISA M1V16 was cleared faster from the blood (P < 0.05) and required 100-times more bacteria than parent M1V2 (10(6) versus 10(4)CFU) to infect 90% of vegetations. GISA M1V16 also had 100 to 1000-times lower bacterial densities in vegetations. Moreover, after co-inoculation with GISA M1V16 and M1V2Rif, a rifampin-resistant variant of M1V2 to discriminate them in organ cultures, GISA M1V16 was out-competed by the glycopeptide-susceptible counterpart. Thus, in rats with experimental endocarditis, GISA showed an attenuated virulence, likely due to a faster clearance from the blood and a reduced fitness in cardiac vegetations. The GISA phenotype appeared globally detrimental to infectivity.

Identifying genetic signatures of selection in a non-model species, alpine gentian (Gentiana nivalis L.), using a landscape genetic approach

It is generally accepted that most plant populations are locally adapted. Yet, understanding how environmental forces give rise to adaptive genetic variation is a challenge in conservation genetics and crucial to the preservation of species under rapidly changing climatic conditions. Environmental variation, phylogeographic history, and population demographic processes all contribute to spatially structured genetic variation, however few current models attempt to separate these confounding effects. To illustrate the benefits of using a spatially-explicit model for identifying potentially adaptive loci, we compared outlier locus detection methods with a recently-developed landscape genetic approach. We analyzed 157 loci from samples of the alpine herb Gentiana nivalis collected across the European Alps. Principle coordinates of neighbor matrices (PCNM), eigenvectors that quantify multi-scale spatial variation present in a data set, were incorporated into a landscape genetic approach relating AFLP frequencies with 23 environmental variables. Four major findings emerged. 1) Fifteen loci were significantly correlated with at least one predictor variable (R (adj) (2) > 0.5). 2) Models including PCNM variables identified eight more potentially adaptive loci than models run without spatial variables. 3) When compared to outlier detection methods, the landscape genetic approach detected four of the same loci plus 11 additional loci. 4) Temperature, precipitation, and solar radiation were the three major environmental factors driving potentially adaptive genetic variation in G. nivalis. Techniques presented in this paper offer an efficient method for identifying potentially adaptive genetic variation and associated environmental forces of selection, providing an important step forward for the conservation of non-model species under global change.

Instruction of haematopoietic lineage choices, evolution of transcriptional landscapes and cancer stem cell hierarchies derived from an AML1-ETO mouse model.

The t(8;21) chromosomal translocation activates aberrant expression of the AML1-ETO (AE) fusion protein and is commonly associated with core binding factor acute myeloid leukaemia (CBF AML). Combining a conditional mouse model that closely resembles the slow evolution and the mosaic AE expression pattern of human t(8;21) CBF AML with global transcriptome sequencing, we find that disease progression was characterized by two principal pathogenic mechanisms. Initially, AE expression modified the lineage potential of haematopoietic stem cells (HSCs), resulting in the selective expansion of the myeloid compartment at the expense of normal erythro- and lymphopoiesis. This lineage skewing was followed by a second substantial rewiring of transcriptional networks occurring in the trajectory to manifest leukaemia. We also find that both HSC and lineage-restricted granulocyte macrophage progenitors (GMPs) acquired leukaemic stem cell (LSC) potential being capable of initiating and maintaining the disease. Finally, our data demonstrate that long-term expression of AE induces an indolent myeloproliferative disease (MPD)-like myeloid leukaemia phenotype with complete penetrance and that acute inactivation of AE function is a potential novel therapeutic option.

National, regional, and global trends in fasting plasma glucose and diabetes prevalence since 1980: systematic analysis of health examination surveys and epidemiological studies with 370 country-years and 2·7 million participants.

BACKGROUND: Data for trends in glycaemia and diabetes prevalence are needed to understand the effects of diet and lifestyle within populations, assess the performance of interventions, and plan health services. No consistent and comparable global analysis of trends has been done. We estimated trends and their uncertainties in mean fasting plasma glucose (FPG) and diabetes prevalence for adults aged 25 years and older in 199 countries and territories. METHODS: We obtained data from health examination surveys and epidemiological studies (370 country-years and 2·7 million participants). We converted systematically between different glycaemic metrics. For each sex, we used a Bayesian hierarchical model to estimate mean FPG and its uncertainty by age, country, and year, accounting for whether a study was nationally, subnationally, or community representative. FINDINGS: In 2008, global age-standardised mean FPG was 5·50 mmol/L (95% uncertainty interval 5·37-5·63) for men and 5·42 mmol/L (5·29-5·54) for women, having risen by 0·07 mmol/L and 0·09 mmol/L per decade, respectively. Age-standardised adult diabetes prevalence was 9·8% (8·6-11·2) in men and 9·2% (8·0-10·5) in women in 2008, up from 8·3% (6·5-10·4) and 7·5% (5·8-9·6) in 1980. The number of people with diabetes increased from 153 (127-182) million in 1980, to 347 (314-382) million in 2008. We recorded almost no change in mean FPG in east and southeast Asia and central and eastern Europe. Oceania had the largest rise, and the highest mean FPG (6·09 mmol/L, 5·73-6·49 for men; 6·08 mmol/L, 5·72-6·46 for women) and diabetes prevalence (15·5%, 11·6-20·1 for men; and 15·9%, 12·1-20·5 for women) in 2008. Mean FPG and diabetes prevalence in 2008 were also high in south Asia, Latin America and the Caribbean, and central Asia, north Africa, and the Middle East. Mean FPG in 2008 was lowest in sub-Saharan Africa, east and southeast Asia, and high-income Asia-Pacific. In high-income subregions, western Europe had the smallest rise, 0·07 mmol/L per decade for men and 0·03 mmol/L per decade for women; North America had the largest rise, 0·18 mmol/L per decade for men and 0·14 mmol/L per decade for women. INTERPRETATION: Glycaemia and diabetes are rising globally, driven both by population growth and ageing and by increasing age-specific prevalences. Effective preventive interventions are needed, and health systems should prepare to detect and manage diabetes and its sequelae. FUNDING: Bill & Melinda Gates Foundation and WHO.

Weighted Shape-based Averaging with Neighborhood Prior Model for Multiple Atlas Fusion-based Medical Image Segmentation

In medical imaging, merging automated segmentations obtained from multiple atlases has become a standard practice for improving the accuracy. In this letter, we propose two new fusion methods: "Global Weighted Shape-Based Averaging" (GWSBA) and "Local Weighted Shape-Based Averaging" (LWSBA). These methods extend the well known Shape-Based Averaging (SBA) by additionally incorporating the similarity information between the reference (i.e., atlas) images and the target image to be segmented. We also propose a new spatially-varying similarity-weighted neighborhood prior model, and an edge-preserving smoothness term that can be used with many of the existing fusion methods. We first present our new Markov Random Field (MRF) based fusion framework that models the above mentioned information. The proposed methods are evaluated in the context of segmentation of lymph nodes in the head and neck 3D CT images, and they resulted in more accurate segmentations compared to the existing SBA.

Assessment of distribution and evolution of mechanical dyssynchrony in a porcine model of myocardial infarction by cardiovascular magnetic resonance.

BACKGROUND: We sought to investigate the relationship between infarct and dyssynchrony post- myocardial infarct (MI), in a porcine model. Mechanical dyssynchrony post-MI is associated with left ventricular (LV) remodeling and increased mortality. METHODS: Cine, gadolinium-contrast, and tagged cardiovascular magnetic resonance (CMR) were performed pre-MI, 9 ± 2 days (early post-MI), and 33 ± 10 days (late post-MI) post-MI in 6 pigs to characterize cardiac morphology, location and extent of MI, and regional mechanics. LV mechanics were assessed by circumferential strain (eC). Electro-anatomic mapping (EAM) was performed within 24 hrs of CMR and prior to sacrifice. RESULTS: Mean infarct size was 21 ± 4% of LV volume with evidence of post-MI remodeling. Global eC significantly decreased post MI (-27 ± 1.6% vs. -18 ± 2.5% (early) and -17 ± 2.7% (late), p < 0.0001) with no significant change in peri-MI and MI segments between early and late time-points. Time to peak strain (TTP) was significantly longer in MI, compared to normal and peri-MI segments, both early (440 ± 40 ms vs. 329 ± 40 ms and 332 ± 36 ms, respectively; p = 0.0002) and late post-MI (442 ± 63 ms vs. 321 ± 40 ms and 355 ± 61 ms, respectively; p = 0.012). The standard deviation of TTP in 16 segments (SD16) significantly increased post-MI: 28 ± 7 ms to 50 ± 10 ms (early, p = 0.012) to 54 ± 19 ms (late, p = 0.004), with no change between early and late post-MI time-points (p = 0.56). TTP was not related to reduction of segmental contractility. EAM revealed late electrical activation and greatly diminished conduction velocity in the infarct (5.7 ± 2.4 cm/s), when compared to peri-infarct (18.7 ± 10.3 cm/s) and remote myocardium (39 ± 20.5 cm/s). CONCLUSIONS: Mechanical dyssynchrony occurs early after MI and is the result of delayed electrical and mechanical activation in the infarct.

Determining the best population-level alcohol consumption model and its impact on estimates of alcohol-attributable harms.

BACKGROUND: The goals of our study are to determine the most appropriate model for alcohol consumption as an exposure for burden of disease, to analyze the effect of the chosen alcohol consumption distribution on the estimation of the alcohol Population- Attributable Fractions (PAFs), and to characterize the chosen alcohol consumption distribution by exploring if there is a global relationship within the distribution. METHODS: To identify the best model, the Log-Normal, Gamma, and Weibull prevalence distributions were examined using data from 41 surveys from Gender, Alcohol and Culture: An International Study (GENACIS) and from the European Comparative Alcohol Study. To assess the effect of these distributions on the estimated alcohol PAFs, we calculated the alcohol PAF for diabetes, breast cancer, and pancreatitis using the three above-named distributions and using the more traditional approach based on categories. The relationship between the mean and the standard deviation from the Gamma distribution was estimated using data from 851 datasets for 66 countries from GENACIS and from the STEPwise approach to Surveillance from the World Health Organization. RESULTS: The Log-Normal distribution provided a poor fit for the survey data, with Gamma and Weibull distributions providing better fits. Additionally, our analyses showed that there were no marked differences for the alcohol PAF estimates based on the Gamma or Weibull distributions compared to PAFs based on categorical alcohol consumption estimates. The standard deviation of the alcohol distribution was highly dependent on the mean, with a unit increase in alcohol consumption associated with a unit increase in the mean of 1.258 (95% CI: 1.223 to 1.293) (R2 = 0.9207) for women and 1.171 (95% CI: 1.144 to 1.197) (R2 = 0. 9474) for men. CONCLUSIONS: Although the Gamma distribution and the Weibull distribution provided similar results, the Gamma distribution is recommended to model alcohol consumption from population surveys due to its fit, flexibility, and the ease with which it can be modified. The results showed that a large degree of variance of the standard deviation of the alcohol consumption Gamma distribution was explained by the mean alcohol consumption, allowing for alcohol consumption to be modeled through a Gamma distribution using only average consumption.

Predicting global change impacts on plant species' distributions: Future challenges

Given the rate of projected environmental change for the 21st century, urgent adaptation and mitigation measures are required to slow down the on-going erosion of biodiversity. Even though increasing evidence shows that recent human-induced environmental changes have already triggered species' range shifts, changes in phenology and species' extinctions, accurate projections of species' responses to future environmental changes are more difficult to ascertain. This is problematic, since there is a growing awareness of the need to adopt proactive conservation planning measures using forecasts of species' responses to future environmental changes. There is a substantial body of literature describing and assessing the impacts of various scenarios of climate and land-use change on species' distributions. Model predictions include a wide range of assumptions and limitations that are widely acknowledged but compromise their use for developing reliable adaptation and mitigation strategies for biodiversity. Indeed, amongst the most used models, few, if any, explicitly deal with migration processes, the dynamics of population at the "trailing edge" of shifting populations, species' interactions and the interaction between the effects of climate and land-use. In this review, we propose two main avenues to progress the understanding and prediction of the different processes A occurring on the leading and trailing edge of the species' distribution in response to any global change phenomena. Deliberately focusing on plant species, we first explore the different ways to incorporate species' migration in the existing modelling approaches, given data and knowledge limitations and the dual effects of climate and land-use factors. Secondly, we explore the mechanisms and processes happening at the trailing edge of a shifting species' distribution and how to implement them into a modelling approach. We finally conclude this review with clear guidelines on how such modelling improvements will benefit conservation strategies in a changing world. (c) 2007 Rubel Foundation, ETH Zurich. Published by Elsevier GrnbH. All rights reserved.

The In Vitro Mass-Produced Model Mycorrhizal Fungus, Rhizophagus irregularis, Significantly Increases Yields of the Globally Important Food Security Crop Cassava.

The arbuscular mycorrhizal symbiosis is formed between arbuscular mycorrhizal fungi (AMF) and plant roots. The fungi provide the plant with inorganic phosphate (P). The symbiosis can result in increased plant growth. Although most global food crops naturally form this symbiosis, very few studies have shown that their practical application can lead to large-scale increases in food production. Application of AMF to crops in the tropics is potentially effective for improving yields. However, a main problem of using AMF on a large-scale is producing cheap inoculum in a clean sterile carrier and sufficiently concentrated to cheaply transport. Recently, mass-produced in vitro inoculum of the model mycorrhizal fungus Rhizophagus irregularis became available, potentially making its use viable in tropical agriculture. One of the most globally important food plants in the tropics is cassava. We evaluated the effect of in vitro mass-produced R. irregularis inoculum on the yield of cassava crops at two locations in Colombia. A significant effect of R. irregularis inoculation on yield occurred at both sites. At one site, yield increases were observed irrespective of P fertilization. At the other site, inoculation with AMF and 50% of the normally applied P gave the highest yield. Despite that AMF inoculation resulted in greater food production, economic analyses revealed that AMF inoculation did not give greater return on investment than with conventional cultivation. However, the amount of AMF inoculum used was double the recommended dose and was calculated with European, not Colombian, inoculum prices. R. irregularis can also be manipulated genetically in vitro, leading to improved plant growth. We conclude that application of in vitro R. irregularis is currently a way of increasing cassava yields, that there is a strong potential for it to be economically profitable and that there is enormous potential to improve this efficiency further in the future.