Life span and reproductive cost explain interspecific variation in the optimal onset of reproduction.

Fitness can be profoundly influenced by the age at first reproduction (AFR), but to date the AFR-fitness relationship only has been investigated intraspecifically. Here, we investigated the relationship between AFR and average lifetime reproductive success (LRS) across 34 bird species. We assessed differences in the deviation of the Optimal AFR (i.e., the species-specific AFR associated with the highest LRS) from the age at sexual maturity, considering potential effects of life history as well as social and ecological factors. Most individuals adopted the species-specific Optimal AFR and both the mean and Optimal AFR of species correlated positively with life span. Interspecific deviations of the Optimal AFR were associated with indices reflecting a change in LRS or survival as a function of AFR: a delayed AFR was beneficial in species where early AFR was associated with a decrease in subsequent survival or reproductive output. Overall, our results suggest that a delayed onset of reproduction beyond maturity is an optimal strategy explained by a long life span and costs of early reproduction. By providing the first empirical confirmations of key predictions of life-history theory across species, this study contributes to a better understanding of life-history evolution.

Global alcohol exposure estimates by country, territory and region for 2005--a contribution to the Comparative Risk Assessment for the 2010 Global Burden of Disease Study.

AIMS: This study aimed to estimate the prevalence of life-time abstainers, former drinkers and current drinkers, adult per-capita consumption of alcohol and pattern of drinking scores, by country and Global Burden of Disease region for 2005, and to forecast these indicators for 2010. DESIGN: Statistical modelling based on survey data and routine statistics. SETTING AND PARTICIPANTS: A total of 241 countries and territories. MEASUREMENTS: Per-capita consumption data were obtained with the help of the World Health Organization's Global Information System on Alcohol and Health. Drinking status data were obtained from Gender, Alcohol and Culture: An International Study, the STEPwise approach to Surveillance study, the World Health Survey/Multi-Country Study and other surveys. Consumption and drinking status data were triangulated to estimate alcohol consumption across multiple categories. FINDINGS: In 2005 adult per-capita annual consumption of alcohol was 6.1 litres, with 1.7 litres stemming from unrecorded consumption; 17.1 litres of alcohol were consumed per drinker, 45.8% of all adults were life-time abstainers, 13.6% were former drinkers and 40.6% were current drinkers. Life-time abstention was most prevalent in North Africa/Middle East and South Asia. Eastern Europe and Southern sub-Saharan Africa had the most detrimental pattern of drinking scores, while drinkers in Europe (Eastern and Central) and sub-Saharan Africa (Southern and West) consumed the most alcohol. CONCLUSIONS: Just over 40% of the world's adult population consumes alcohol and the average consumption per drinker is 17.1 litres per year. However, the prevalence of abstention, level of alcohol consumption and patterns of drinking vary widely across regions of the world.

Probability of achieving optimal molecular response to imatinib in chronic myeloid leukemia (CML) patients: Pharmacokinetic/Pharmacodynamic (PK/PD) relationships observed under field-conditions

Objectives: Imatinib has been increasingly proposed for therapeutic drug monitoring (TDM), as trough concentrations (Cmin) correlate with response rates in CML patients. This analysis aimed to evaluate the impact of imatinib exposure on optimal molecular response rates in a large European cohort of patients followed by centralized TDM.¦Methods: Sequential PK/PD analysis was performed in NONMEM 7 on 2230 plasma (PK) samples obtained along with molecular response (PD) data from 1299 CML patients. Model-based individual Bayesian estimates of exposure, parameterized as to initial dose adjusted and log-normalized Cmin (log-Cmin) or clearance (CL), were investigated as potential predictors of optimal molecular response, while accounting for time under treatment (stratified at 3 years), gender, CML phase, age, potentially interacting comedication, and TDM frequency. PK/PD analysis used mixed-effect logistic regression (iterative two-stage method) to account for intra-patient correlation.¦Results: In univariate analyses, CL, log-Cmin, time under treatment, TDM frequency, gender (all p<0.01) and CML phase (p=0.02) were significant predictors of the outcome. In multivariate analyses, all but log-Cmin remained significant (p<0.05). Our model estimates a 54.1% probability of optimal molecular response in a female patient with a median CL of 14.4 L/h, increasing by 4.7% with a 35% decrease in CL (percentile 10 of CL distribution), and decreasing by 6% with a 45% increased CL (percentile 90), respectively. Male patients were less likely than female to be in optimal response (odds ratio: 0.62, p<0.001), with an estimated probability of 42.3%.¦Conclusions: Beyond CML phase and time on treatment, expectedly correlated to the outcome, an effect of initial imatinib exposure on the probability of achieving optimal molecular response was confirmed in field-conditions by this multivariate analysis. Interestingly, male patients had a higher risk of suboptimal response, which might not exclusively derive from their 18.5% higher CL, but also from reported lower adherence to the treatment. A prospective longitudinal study would be desirable to confirm the clinical importance of identified covariates and to exclude biases possibly affecting this observational survey.

Optimal robust estimates using the Hellinger distance

Optimal robust M-estimates of a multidimensional parameter are described using Hampel's infinitesimal approach. The optimal estimates are derived by minimizing a measure of efficiency under the model, subject to a bounded measure of infinitesimal robustness. To this purpose we define measures of efficiency and infinitesimal sensitivity based on the Hellinger distance.We show that these two measures coincide with similar ones defined by Yohai using the Kullback-Leibler divergence, and therefore the corresponding optimal estimates coincide too.We also give an example where we fit a negative binomial distribution to a real dataset of "days of stay in hospital" using the optimal robust estimates.

Hydrogeophysical data integration at larger scales: Application of Bayesian sequential simulation for the characterization of heterogeneous alluvial aquifers

Knowledge of the spatial distribution of hydraulic conductivity (K) within an aquifer is critical for reliable predictions of solute transport and the development of effective groundwater management and/or remediation strategies. While core analyses and hydraulic logging can provide highly detailed information, such information is inherently localized around boreholes that tend to be sparsely distributed throughout the aquifer volume. Conversely, larger-scale hydraulic experiments like pumping and tracer tests provide relatively low-resolution estimates of K in the investigated subsurface region. As a result, traditional hydrogeological measurement techniques contain a gap in terms of spatial resolution and coverage, and they are often alone inadequate for characterizing heterogeneous aquifers. Geophysical methods have the potential to bridge this gap. The recent increased interest in the application of geophysical methods to hydrogeological problems is clearly evidenced by the formation and rapid growth of the domain of hydrogeophysics over the past decade (e.g., Rubin and Hubbard, 2005).

Early broncoconstriction in anesthetized sheep: a simplified experimental model of asthma

Résumé Cette étude décrit un modèle expérimental de bronchoconstriction précoce induite par aérosolisation d'un extrait d'Ascaris suum chez des moutons anesthésiés par de l'isoflurane et ventilés mécaniquement. Dix moutons adultes ont été anesthésiés et ventilés mécaniquement puis ont été exposés à un stimulus bronchoconstrictif sous forme d'un aérosol d'extrait d'Ascaris suum durant 25 minutes. Tous les moutons ont été exposés deux fois à huit semaines d'intervalle à ce même stimulus. Les échanges gazeux ainsi que les paramètres respiratoires ont été mesurés régulièrement durant la période d'aérosolisation ainsi que durant les 60 minutes suivantes. A la fin de la période d'aérosolisation, une augmentation significative (p<0.05) des pressions de crête (+114%) et de plateau (+148%), de la résistance expiratoire (+93%) et de la pression partielle artérielle de gaz carbonique PaCO2 (+25%) a été constatée, de même qu'une diminution significative (p<0.05) de la compliance respiratoire (-41 %) et de la pression partielle artérielle d'oxygène PaO2 (-49%). Ces modifications sont restées stables durant toute la période d'observation. Ce modèle expérimental animal de bronchoconstriction offre de nombreux avantages : la stabilité hémodynamique et le confort de l'animal sont améliorés et la réaction de stress est inhibée. Il permet de plus une distribution optimale de l'antigène respiratoire et finalement évite l'utilisation d'un pléthysmographe corporel. Abstract This study describes a simplified experimental model of early bronchoconstriction induced by aerosolization of Ascaris suum extract in isoflurane-anesthetized and mechanically ventilated sheep. Ten adult sheep were anesthetized, mechanically ventilated and then challenged with an aerosol of Ascaris suum extract during 25 minutes. All of them were challenged twice at eight weeks intervals. During the bronchoconstrictive challenges and the following sixty minutes, gas exchange was measured and respiratory mechanics parameters computed from a lung mechanics calculator. At the end of the challenge, a significant increase (p<0.05) was observed in peak (+114%) and plateau (+148%) pressures, expiratory resistance (+93%) and PaCO2 (+25%) along with a significant decrease (p<0.05) in respiratory compliance (-41 %) and PaO2 (-49%). These changes remained stable throughout the 60 minutes study period. This model offers several advantages: hemodynamic stability and animal welfare are improved and the stress response is blunted. It allows an optimal distribution of the antigen and finally avoids the need of a body plethysmograph.

Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers

Recent technological advances in remote sensing have enabled investigation of the morphodynamics and hydrodynamics of large rivers. However, measuring topography and flow in these very large rivers is time consuming and thus often constrains the spatial resolution and reach-length scales that can be monitored. Similar constraints exist for computational fluid dynamics (CFD) studies of large rivers, requiring maximization of mesh-or grid-cell dimensions and implying a reduction in the representation of bedform-roughness elements that are of the order of a model grid cell or less, even if they are represented in available topographic data. These ``subgrid'' elements must be parameterized, and this paper applies and considers the impact of roughness-length treatments that include the effect of bed roughness due to ``unmeasured'' topography. CFD predictions were found to be sensitive to the roughness-length specification. Model optimization was based on acoustic Doppler current profiler measurements and estimates of the water surface slope for a variety of roughness lengths. This proved difficult as the metrics used to assess optimal model performance diverged due to the effects of large bedforms that are not well parameterized in roughness-length treatments. However, the general spatial flow patterns are effectively predicted by the model. Changes in roughness length were shown to have a major impact upon flow routing at the channel scale. The results also indicate an absence of secondary flow circulation cells in the reached studied, and suggest simpler two-dimensional models may have great utility in the investigation of flow within large rivers. Citation: Sandbach, S. D. et al. (2012), Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers, Water Resour. Res., 48, W12501, doi: 10.1029/2011WR011284.

Genetic structure of Gymnures (genus Hylomys; Erinaceidae) on continental islands of Southeast Asia: Historical effects of fragmentation

Eustatic sea level changes during Pleistocene climatic fluctuations produced several cycles of connection-isolation among continental islands of the Sunda shelf. To explore the potential effects of these fluctuations, we reconstructed a model of the vicariant events that separated these islands, based on bathymetric information. Among many possible scenarios, two opposite phylogenetic patterns of evolution were predicted for terrestrial organisms living in this region: one is based on the classical allopatric speciation mode of evolution, while the other is the outcome of a sequential dispersal colonization of the archipelago. We tested the applicability of these predictions with an analysis of sequence variation of the cytochrome b gene from several taxa of Hylomys. They were sampled throughout SE-Asia and the Sunda islands. High levels of haplotype differentiation characterize the different island taxa. Such levels of differentiation support the existence of several allopatric species, as was suggested by previous allozyme and morphological data. Also in accordance with previous results, the occurrence of two sympatric species from Sumatra is suggested by their strongly divergent haplotypes. One species, Hylomys suillus maxi, is found both on Sumatra and in Peninsular Malaysia, while the other, H. parvus, is endemic to Sumatra. Its closest relative is H. suillus dorsalis from Borneo. Phylogenetic reconstructions also demonstrate the existence of a Sundaic clade composed of all island taxa, as opposed to those from the continent. Although there is no statistical support for either proposed biogeographic model of evolution, we argue that the sequential dispersal scenario is more appropriate to describe the genetic variation found among the Hylomys taxa. However, despite strong differentiation among island haplotypes, the cladistic relationships between some island taxa could not be resolved. We argue that this is evidence of a rapid radiation, suggesting that the separation of the islands may have been perceived as a simultaneous event rather than as a succession of vicariant events. Furthermore, the estimates of divergence times between the haplotypes of these taxa suggest that this radiation may actually have predated the climatic fluctuations of the Pleistocene. Further refinement of the initial palaeogeographic models of evolution are therefore needed to account for these results.

Binary classification of dyslipidemia from the waist-to-hip ratio and body mass index: a comparison of linear, logistic, and CART models.

BACKGROUND: We sought to improve upon previously published statistical modeling strategies for binary classification of dyslipidemia for general population screening purposes based on the waist-to-hip circumference ratio and body mass index anthropometric measurements. METHODS: Study subjects were participants in WHO-MONICA population-based surveys conducted in two Swiss regions. Outcome variables were based on the total serum cholesterol to high density lipoprotein cholesterol ratio. The other potential predictor variables were gender, age, current cigarette smoking, and hypertension. The models investigated were: (i) linear regression; (ii) logistic classification; (iii) regression trees; (iv) classification trees (iii and iv are collectively known as "CART"). Binary classification performance of the region-specific models was externally validated by classifying the subjects from the other region. RESULTS: Waist-to-hip circumference ratio and body mass index remained modest predictors of dyslipidemia. Correct classification rates for all models were 60-80%, with marked gender differences. Gender-specific models provided only small gains in classification. The external validations provided assurance about the stability of the models. CONCLUSIONS: There were no striking differences between either the algebraic (i, ii) vs. non-algebraic (iii, iv), or the regression (i, iii) vs. classification (ii, iv) modeling approaches. Anticipated advantages of the CART vs. simple additive linear and logistic models were less than expected in this particular application with a relatively small set of predictor variables. CART models may be more useful when considering main effects and interactions between larger sets of predictor variables.

Neuroanatomical and neuropsychological features of euthymic patients with bipolar disorder.

OBJECTIVE: Previous studies reported that the severity of cognitive deficits in euthymic patients with bipolar disorder (BD) increases with the duration of illness and postulated that progressive neuronal loss or shrinkage and white matter changes may be at the origin of this phenomenon. To explore this issue, the authors performed a case-control study including detailed neuropsychological and magnetic resonance imaging analyses in 17 euthymic elderly patients with BD and 17 healthy individuals. METHODS: Neuropsychological evaluation concerned working memory, episodic memory, processing speed, and executive functions. Volumetric estimates of the amygdala, hippocampus, entorhinal cortex, and anterior cingulate cortex were obtained using both voxel-based and region of interest morphometric methods. Periventricular and deep white matter were assessed semiquantitatively. Differences in cognitive performances and structural data between BD and comparison groups were analyzed using paired t-test or analysis of variance. Wilcoxon test was used in the absence of normal distribution. RESULTS: Compared with healthy individuals, patients with BD obtained significantly lower performances in processing speed, working memory, and episodic memory but not in executive functions. Morphometric analyses did not show significant volumetric or white matter differences between the two groups. CONCLUSIONS: Our results revealed impairment in verbal memory, working memory, and processing speed in euthymic older adults with BD. These cognitive deficits are comparable both in terms of affected functions and size effects to those previously reported in younger cohorts with BD. Both this observation and the absence of structural brain abnormalities in our cohort do not support a progressively evolving neurotoxic effect in BD.

Estimation of the lateral correlation structure of subsurface water content from surface-based ground-penetrating radar reflection images

Over the past decade, significant interest has been expressed in relating the spatial statistics of surface-based reflection ground-penetrating radar (GPR) data to those of the imaged subsurface volume. A primary motivation for this work is that changes in the radar wave velocity, which largely control the character of the observed data, are expected to be related to corresponding changes in subsurface water content. Although previous work has indeed indicated that the spatial statistics of GPR images are linked to those of the water content distribution of the probed region, a viable method for quantitatively analyzing the GPR data and solving the corresponding inverse problem has not yet been presented. Here we address this issue by first deriving a relationship between the 2-D autocorrelation of a water content distribution and that of the corresponding GPR reflection image. We then show how a Bayesian inversion strategy based on Markov chain Monte Carlo sampling can be used to estimate the posterior distribution of subsurface correlation model parameters that are consistent with the GPR data. Our results indicate that if the underlying assumptions are valid and we possess adequate prior knowledge regarding the water content distribution, in particular its vertical variability, this methodology allows not only for the reliable recovery of lateral correlation model parameters but also for estimates of parameter uncertainties. In the case where prior knowledge regarding the vertical variability of water content is not available, the results show that the methodology still reliably recovers the aspect ratio of the heterogeneity.

Tightening force and torque of nonlocking screws in a reverse shoulder prosthesis.

BACKGROUND: Reversed shoulder arthroplasty is an accepted treatment for glenohumeral arthritis associated to rotator cuff deficiency. For most reversed shoulder prostheses, the baseplate of the glenoid component is uncemented and its primary stability is provided by a central peg and peripheral screws. Because of the importance of the primary stability for a good osteo-integration of the baseplate, the optimal fixation of the screws is crucial. In particular, the amplitude of the tightening force of the nonlocking screws is clearly associated to this stability. Since this force is unknown, it is currently not accounted for in experimental or numerical analyses. Thus, the primary goal of this work is to measure this tightening force experimentally. In addition, the tightening torque was also measured, to estimate an optimal surgical value. METHODS: An experimental setup with an instrumented baseplate was developed to measure simultaneously the tightening force, tightening torque and screwing angle, of the nonlocking screws of the Aquealis reversed prosthesis. In addition, the amount of bone volume around each screw was measured with a micro-CT. Measurements were performed on 6 human cadaveric scapulae. FINDINGS: A statistically correlated relationship (p<0.05, R=0.83) was obtained between the maximal tightening force and the bone volume. The relationship between the tightening torque and the bone volume was not statistically significant. INTERPRETATION: The experimental relationship presented in this paper can be used in numerical analyses to improve the baseplate fixation in the glenoid bone.

Robust response transformations based on optimal prediction

Nonlinear regression problems can often be reduced to linearity by transforming the response variable (e.g., using the Box-Cox family of transformations). The classic estimates of the parameter defining the transformation as well as of the regression coefficients are based on the maximum likelihood criterion, assuming homoscedastic normal errors for the transformed response. These estimates are nonrobust in the presence of outliers and can be inconsistent when the errors are nonnormal or heteroscedastic. This article proposes new robust estimates that are consistent and asymptotically normal for any unimodal and homoscedastic error distribution. For this purpose, a robust version of conditional expectation is introduced for which the prediction mean squared error is replaced with an M scale. This concept is then used to develop a nonparametric criterion to estimate the transformation parameter as well as the regression coefficients. A finite sample estimate of this criterion based on a robust version of smearing is also proposed. Monte Carlo experiments show that the new estimates compare favorably with respect to the available competitors.

Looking for orders of magnitude : five essays in trade and the environment and economic geography

General Summary Although the chapters of this thesis address a variety of issues, the principal aim is common: test economic ideas in an international economic context. The intention has been to supply empirical findings using the largest suitable data sets and making use of the most appropriate empirical techniques. This thesis can roughly be divided into two parts: the first one, corresponding to the first two chapters, investigates the link between trade and the environment, the second one, the last three chapters, is related to economic geography issues. Environmental problems are omnipresent in the daily press nowadays and one of the arguments put forward is that globalisation causes severe environmental problems through the reallocation of investments and production to countries with less stringent environmental regulations. A measure of the amplitude of this undesirable effect is provided in the first part. The third and the fourth chapters explore the productivity effects of agglomeration. The computed spillover effects between different sectors indicate how cluster-formation might be productivity enhancing. The last chapter is not about how to better understand the world but how to measure it and it was just a great pleasure to work on it. "The Economist" writes every week about the impressive population and economic growth observed in China and India, and everybody agrees that the world's center of gravity has shifted. But by how much and how fast did it shift? An answer is given in the last part, which proposes a global measure for the location of world production and allows to visualize our results in Google Earth. A short summary of each of the five chapters is provided below. The first chapter, entitled "Unraveling the World-Wide Pollution-Haven Effect" investigates the relative strength of the pollution haven effect (PH, comparative advantage in dirty products due to differences in environmental regulation) and the factor endowment effect (FE, comparative advantage in dirty, capital intensive products due to differences in endowments). We compute the pollution content of imports using the IPPS coefficients (for three pollutants, namely biological oxygen demand, sulphur dioxide and toxic pollution intensity for all manufacturing sectors) provided by the World Bank and use a gravity-type framework to isolate the two above mentioned effects. Our study covers 48 countries that can be classified into 29 Southern and 19 Northern countries and uses the lead content of gasoline as proxy for environmental stringency. For North-South trade we find significant PH and FE effects going in the expected, opposite directions and being of similar magnitude. However, when looking at world trade, the effects become very small because of the high North-North trade share, where we have no a priori expectations about the signs of these effects. Therefore popular fears about the trade effects of differences in environmental regulations might by exaggerated. The second chapter is entitled "Is trade bad for the Environment? Decomposing worldwide SO2 emissions, 1990-2000". First we construct a novel and large database containing reasonable estimates of SO2 emission intensities per unit labor that vary across countries, periods and manufacturing sectors. Then we use these original data (covering 31 developed and 31 developing countries) to decompose the worldwide SO2 emissions into the three well known dynamic effects (scale, technique and composition effect). We find that the positive scale (+9,5%) and the negative technique (-12.5%) effect are the main driving forces of emission changes. Composition effects between countries and sectors are smaller, both negative and of similar magnitude (-3.5% each). Given that trade matters via the composition effects this means that trade reduces total emissions. We next construct, in a first experiment, a hypothetical world where no trade happens, i.e. each country produces its imports at home and does no longer produce its exports. The difference between the actual and this no-trade world allows us (under the omission of price effects) to compute a static first-order trade effect. The latter now increases total world emissions because it allows, on average, dirty countries to specialize in dirty products. However, this effect is smaller (3.5%) in 2000 than in 1990 (10%), in line with the negative dynamic composition effect identified in the previous exercise. We then propose a second experiment, comparing effective emissions with the maximum or minimum possible level of SO2 emissions. These hypothetical levels of emissions are obtained by reallocating labour accordingly across sectors within each country (under the country-employment and the world industry-production constraints). Using linear programming techniques, we show that emissions are reduced by 90% with respect to the worst case, but that they could still be reduced further by another 80% if emissions were to be minimized. The findings from this chapter go together with those from chapter one in the sense that trade-induced composition effect do not seem to be the main source of pollution, at least in the recent past. Going now to the economic geography part of this thesis, the third chapter, entitled "A Dynamic Model with Sectoral Agglomeration Effects" consists of a short note that derives the theoretical model estimated in the fourth chapter. The derivation is directly based on the multi-regional framework by Ciccone (2002) but extends it in order to include sectoral disaggregation and a temporal dimension. This allows us formally to write present productivity as a function of past productivity and other contemporaneous and past control variables. The fourth chapter entitled "Sectoral Agglomeration Effects in a Panel of European Regions" takes the final equation derived in chapter three to the data. We investigate the empirical link between density and labour productivity based on regional data (245 NUTS-2 regions over the period 1980-2003). Using dynamic panel techniques allows us to control for the possible endogeneity of density and for region specific effects. We find a positive long run elasticity of density with respect to labour productivity of about 13%. When using data at the sectoral level it seems that positive cross-sector and negative own-sector externalities are present in manufacturing while financial services display strong positive own-sector effects. The fifth and last chapter entitled "Is the World's Economic Center of Gravity Already in Asia?" computes the world economic, demographic and geographic center of gravity for 1975-2004 and compares them. Based on data for the largest cities in the world and using the physical concept of center of mass, we find that the world's economic center of gravity is still located in Europe, even though there is a clear shift towards Asia. To sum up, this thesis makes three main contributions. First, it provides new estimates of orders of magnitudes for the role of trade in the globalisation and environment debate. Second, it computes reliable and disaggregated elasticities for the effect of density on labour productivity in European regions. Third, it allows us, in a geometrically rigorous way, to track the path of the world's economic center of gravity.

Optimization of preservation conditions of As (III) bioreporter bacteria.

Long-term preservation of bioreporter bacteria is essential for the functioning of cell-based detection devices, particularly when field application, e.g., in developing countries, is intended. We varied the culture conditions (i.e., the NaCl content of the medium), storage protection media, and preservation methods (vacuum drying vs. encapsulation gels remaining hydrated) in order to achieve optimal preservation of the activity of As (III) bioreporter bacteria during up to 12 weeks of storage at 4 degrees C. The presence of 2% sodium chloride during the cultivation improved the response intensity of some bioreporters upon reconstitution, particularly of those that had been dried and stored in the presence of sucrose or trehalose and 10% gelatin. The most satisfying, stable response to arsenite after 12 weeks storage was obtained with cells that had been dried in the presence of 34% trehalose and 1.5% polyvinylpyrrolidone. Amendments of peptone, meat extract, sodium ascorbate, and sodium glutamate preserved the bioreporter activity only for the first 2 weeks, but not during long-term storage. Only short-term stability was also achieved when bioreporter bacteria were encapsulated in gels remaining hydrated during storage.