BACKGROUND
recently considerable attention has been given to the characterization and understanding of intrinsic molecular noise in biological systems  <cit> .
nearly all of these studies were performed using in vivo fluorescent reporters in single cell systems.
in multicellular organisms, however, most quantitative gene expression data are obtained from fixed tissues.
examples of such data for the drosophila segmentation system are contained in the flyex database, which provides spatiotemporal data on the expression of developmental segmentation genes  <cit> .
these data on protein expression levels are at cellular resolution and were obtained by means of immunofluorescence histochemistry and confocal scanning microscopy.
at a large spatial scale, expression levels in these embryos form expression domains characteristic for each gene, but smaller fluctuations in expression levels between adjacent nuclei appear random.
in this paper, we investigate the question of whether these fluctuations are a consequence of intrinsic molecular noise or stem from some type of measurement uncertainty.
these alternatives are, of course, not mutually exclusive.
a complicating factor in separating the above alternatives is that each one involves an unknown chemical mechanism.
intrinsic noise will have a major contribution from fluctuations in the rate of initiation of transcription, but the chemical mechanisms underlying this process in multicellular organisms are very poorly understood.
measurement uncertainty can stem from chemical causes such as fluctuations in the number of primary and secondary antibody molecules which bind to proteins in the fixed embryo.
the chemistry of this process is also very poorly understood.
if observations could be made on a process whose fluctuation properties could be reliably predicted by a numerical model, comparison of the predicted fluctuations with those observed will provide critical information for distinguishing whether observed nucleus-to-nucleus variations are a consequence of intrinsic biological noise or merely fluctuations arising from the staining procedure.
an excellent candidate for a process with predictable fluctuation properties is one that involves only diffusion and first order decay.
there is good evidence that the formation of the protein gradient of the morphogen bicoid  takes place by means of these two processes.
bcd protein is distributed in an exponential profile along the anterior-posterior  axis with higher concentrations towards the anterior  <cit> .
this gradient forms by translation of maternally deposited mrna at the anterior pole of the embryo, and the synthesized protein spreads through the syncytial embryo by diffusion accompanied by decay  <cit> .
the observed exponential profile corresponds to a solution of fick's equation for a substance undergoing first order decay and diffusing from a point source in one dimension, and hence it is reasonable to suppose that the reaction-diffusion mechanism leading to the formation of the gradient is reasonably well understood.
quantitative observations of this gradient in fixed tissue exhibit small fluctuations between neighboring nuclei, while the overall exponential profile ensures that such fluctuations can be monitored over a wide range of concentrations which extend to the lower limits of detectability.
the intrinsic fluctuations of the bcd gradient will be well described by a stochastic reaction-diffusion master equation .
such equations typically do not have analytic solutions and are usually solved by running repeated stochastic simulations.
a well known simulation algorithm due to gillespie  <cit>  performs an exact simulation of the chemical master equation for a well mixed system.
this method has been extended to spatially distributed systems by elf and others  <cit> .
these authors divide the spatially extended system into a series of subvolumes that are small enough to be regarded as well mixed.
in each subvolume chemical reactions are simulated by gillespie's algorithm, while diffusion between subvolumes is treated as a first order reaction.
in this study, we compare the results of such simulations to data in order to discover whether or not the data is sufficiently accurate as to exhibit the signature of a simple stochastic process.
stochastic processes underlying biological regulation can in general form complex patterns as a result of the reaction network, and for this reason a full consideration of  <dig> dimensional geometry is often necessary  <cit> .
in the case considered here fluctuations occur passively in the course of diffusion and the statistical signature that we seek is independent of detailed geometry.
for these reasons, we chose to model a  <dig> dimensional system.
RESULTS
in the following, we first characterize the statistical properties of random variations in expression level between adjacent nuclei of individual embryos and compare them to the results of stochastic simulations.
at the level of single embryos, we do not find a clear signature of stochastic processes in the data, but the need to separate spatially changing mean expression values from their variation limits the amount of statistical information that can be obtained from individual embryos.
to address this problem, we consider ensembles of embryos, both over the whole dataset and over a restricted subset with low embryo-to-embryo variation.
in order to interpret these data, we introduce a stochastic model of the immunostaining procedure.
we denote all random variables in this article with an upper hat but write a specific value without the hat.thus, for example, n^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ denotes a random variable for the number of bcd molecules in subvolume j, while nj denotes a particular value of this variable.
single embryo analysis in the logarithmic domain
we find that the bcd profiles of individual embryos observed by immunostaining  are exponential, as previously reported  <cit> .
this profile strongly supports the model of an effective fickian diffusion which gives an exponential decay solution at steady state  <cit> .
each embryo contains a collection of observations of expression  ij, and the observed properties of the exponential gradient will depend on the data through a function f embodying the least squares fitting procedure such that
f  = a exp.
the residuals of the exponential vary from embryo to embryo in our data.
in logarithmic coordinates the size of the residuals is independent of j for the anterior portion of the embryo, and in certain embryos  the size of the residuals is completely independent of j .
for these embryos, the residuals are well described by
ln⁡=ln⁡−j/λ+w_,
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagagiibawmaeioba4maeiikagiaemysak0aasbaasqaaiabdqgaqbqabagccqggpaqkcqgh9aqpcyggsbabcqggubgbcqggoaakcqwghbqycqggpaqkcqghsislcqwgqbgacqggvawliigacqwf7oabcqghrawkdaqiaaqaaiabdefaxbgaaypadagaeiilawcaaa@42b6@
where w_
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaaecaaeaacqwgxbwvaiaawkwaaaaa@2dcc@ is an unknown random variable independent of j, so that the variance of ln is equal to the variance of w_
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaaecaaeaacqwgxbwvaiaawkwaaaaa@2dcc@.
intrinsic noise is insufficient
in order to understand whether the observed characteristic variation is a consequence of intrinsic fluctuations in bcd molecular number, we performed stochastic simulations of an rdme which describes the time evolution of the bcd gradient and compared them to data.
we imagine the observed intensity ij to be a particular value of the random variable i^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@.
we further suppose that the random variable i^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@ is determined by a direct linear rescaling of the bcd molecular number such that
i^j=mn^j,
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9iabd2gatjqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqggsaalaaa@3530@
where the factor m represents the proportionality between one bcd molecule and the corresponding fluorescence intensity, which includes the combined effects of tissue fixation, first and second antibody binding, fluorescence excitation and image processing.
as we do not know the exact in vivo system parameters and molecular number within each nucleus, we performed a complete inspection of the behavior of the variance of mn^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ in the four dimensional parameter space f, where j is the synthesis rate of bcd in the anterior compartment, ω is its decay rate, and d is the diffusion coefficient.
it was always true that the residuals in logarithmic coordinate increased towards the posterior of the embryo, or in other words at lower levels of bcd.
this behavior strongly contrasts with the position-independent residuals seen in figure  <dig> and described in equation .
measurement rescaling noise dominates
the simulation results suggest that intrinsic noise cannot explain the pattern of variance seen in figure  <dig>  another possibility is the measurement process itself.
in order to analyze this process, we consider a simple model of the measurement of fluorescence intensity, where
i^j=α^n^j+β^.
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9ggaciqb=f7ahzaajagafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgucariqb=j7aizaajaacbagae4nla4caaa@381a@
here α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ is a spatially uniform random variable which replaces m in equation , and β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ is a spatially uniform random variable which represents nonspecific background staining.
this picture allows us to consider noise that arises from both intrinsic and measurement related sources.
the simplest way to understand the consequences of  is to imagine the consequences if only one of α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@, n^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@, and β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ are allowed to have finite variance while the other two are constrained to deterministic  behavior.
we have already pointed out that allowing finite variance for n^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ with α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ and β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ deterministic leads to an increase of variance towards the posterior.
the same is true if all noise comes from β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@, that is if noise from background staining dominates.
this will also lead to more noise in the logarithmic domain towards the posterior as β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ provides a larger proportion of the total detected signal.
let us next consider the case where all noise comes from α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@.
in order to understand the role of α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@, note that the spatial pattern of variance observed in figure  <dig> can be captured by an exponential function multiplied by a normal random variable.
this suggests that the simplest picture for  is given by assuming that there is no background noise β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ and no intrinsic noise for bcd so that n^j=〈n^j〉
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaakiabg2da9iabgmyihlqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqghqms8aaa@366c@.
moreover, it suggests that the rescaling noise α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ from measurement uncertainty is uniform across the embryo  and is normally distributed with
α^=m),
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaacqgh9aqpcqwgtbqbcqggoaakcqaixaqmcqghrawkcqwfdpwcdawgaawcbagae8xsdegabeaakiqbd6eaozaajagaeiikagiaegimaajaeiilawiaegymaejaeiykakiaeiykakiaeiilawcaaa@3dd7@
where n^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmota4kbakaacqggoaakcqaiwaamcqggsaalcqaixaqmcqggpaqkaaa@3178@ is a normal independent random variable with mean  <dig> and variance  <dig>  then we can model i^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@ by
i^j=α^〈n^j〉.
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9ggaciqb=f7ahzaajagaeyykjeuafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgqyixlabc6cauaaa@390a@
in steady state, 〈n^j〉
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeyykjeuafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgqyixdaa@325e@ = a/m exp.
taking logarithms allows us to write
ln⁡=ln⁡−j/λ+ln⁡).
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagagiibawmaeioba4maeiikagiafmysakkbakaadawgaawcbagaemoaaogabeaakiabcmcapiabg2da9igbcygasjabc6gaujabcicaoiabdggahjabcmcapiabgkhitiabdqgaqjabc+cavggaciab=t7asjabgucarigbcygasjabc6gaujabcicaoiabigdaxiabgucariab=n8aznaabaaaleaacqwfxoqyaeqaaogafmota4kbakaacqggoaakcqaiwaamcqggsaalcqaixaqmcqggpaqkcqggpaqkcqgguaglaaa@504e@
comparison with equation  indicates that w_
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaaecaaeaacqwgxbwvaiaawkwaaaaa@2dcc@ = ln
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmota4kbakaacqggoaakcqaiwaamcqggsaalcqaixaqmcqggpaqkaaa@3178@).
the reason we do not see intrinsic noise in this individual embryo is most likely low measurement sensitivity, that is to say a low molecule-to-fluorescence mean rescaling ratio m. when m is small enough, it can mask the variance of n^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ because the variance of the rescaled gradient is given by
var⁡=m2var⁡.
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagagiodaynaeiyyaemaeiocainaeiikagiaemyba0mafmoba4mbakaadawgaawcbagaemoaaogabeaakiabcmcapiabg2da9iabd2gatnaacaaaleqabagaegomaidaaogagiodaynaeiyyaemaeiocainaeiikagiafmoba4mbakaadawgaawcbagaemoaaogabeaakiabcmcapiabc6cauaaa@43c6@
as m <dig> →  <dig> with var bounded in a reasonable way throughout the embryo, we will get var →  <dig>  hence the rescaled gradient mn^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ can be treated as deterministic by letting
mn^j=〈mn^j〉=m〈n^j〉.
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaemyba0mafmoba4mbakaadawgaawcbagaemoaaogabeaakiabg2da9iabgmyihlabd2gatjqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqghqms8cqgh9aqpcqwgtbqbcqghpms4cuwgubgbgaqcamaabaaaleaacqwgqbgaaeqaaogaeyokjevaeiola4caaa@4358@
in summary, we suspect the nucleus-to-nucleus variation observed in our data comes chiefly from the experimental rescaling noise α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@, which is normally distributed.
if bcd intrinsic noise is to be observed, then the fluorescence noise intensity should be a function of the mean intensity in logarithm, instead of a constant as observed in embryo ms <dig>  nevertheless, the necessity of considering data in spatially resolved bins limits the amount of information that can be obtained from a single embryo.
more information can be obtained by pooling data from many embryos, and we discuss this point in the next section.
statistical analysis of an ensemble of embryos
statistical analysis of many embryos is required in order to take our analysis further.
this analysis will show how physical constraints on the model can be inferred from the ensemble dataset, and independent random variables separated.
we consider a set of embryos indexed by i with expression levels iij, and then pool data from corresponding bins to obtain the ensemble dataset ∪iiij
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaambuaeaacqwgjbqsdawgaawcbagaemyaakmaemoaaogabeaaaeaacqwgpbqaaeqaniablqiivbaaaa@32d6@.
since this dataset includes embryo-to-embryo variability, i.e.
the variation of system parameters and experimental conditions over different embryos, the variance of the ensemble profile will be an upper bound for the average variance within each embryo.
this allows us to identify physical constraints for system parameters and to determine if the model behaves properly within the permitted range of parameters.
we define independent global random variables i^j=α^n^j+β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9ggaciqb=f7ahzaajagafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgucariqb=j7aizaajaaaaa@36e3@ as described in the last section with normal distributed measurement uncertainty α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ = m
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmota4kbakaacqggoaakcqaiwaamcqggsaalcqaixaqmcqggpaqkaaa@3178@).
we also now assume that background noise is normally distributed with β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ = σβ n^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmota4kbakaacqggoaakcqaiwaamcqggsaalcqaixaqmcqggpaqkaaa@3178@.
we assume such global random variable represent the average variability for each embryo.
the statistics of simulated global random variables were then collected from  <dig> stochastic simulation runs and the bcd molecular number random variable n^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ were sampled after reaching steady state.
in this section we explicitly consider the effects of different choices for the molecule-to-fluorescence rescaling ratio m by comparing simulations to data at differing values of this parameter.
because m is not an explicit input to the model, this comparison is effected by varying the synthesis rate j, which varies the molecular number, and comparing the behavior of the model to fluorescence data which is on a fixed but arbitrary scale.
we seek optimal values of m such that the simulated global random variable i^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@ is constrained by the variation observed in the immunostained ensemble data ∪iiij
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaambuaeaacqwgjbqsdawgaawcbagaemyaakmaemoaaogabeaaaeaacqwgpbqaaeqaniablqiivbaaaa@32d6@ by the condition
〈i^j〉=〈∪iiij〉
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeyykjeuafmysakkbakaadawgaawcbagaemoaaogabeaakiabgqyixlabg2da9iabgmyihpaatafabagaemysak0aasbaasqaaiabdmgapjabdqgaqbqabaaabagaemyaakgabeqdcqwiqisvaogaeyokjepaaa@3df8@
var⁡≤var⁡.
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagagiodaynaeiyyaemaeiocainaeiikagiafmysakkbakaadawgaawcbagaemoaaogabeaakiabcmcapiabgsmijkgbcaha2jabcggahjabckhayjabcicaomaatafabagaemysak0aasbaasqaaiabdmgapjabdqgaqbqabaaabagaemyaakgabeqdcqwiqisvaogaeiykakiaeiola4caaa@4441@
because comparison with mean and variance of the ensemble data is not straightforward, we approximate the above conditions by comparing their exponential fit, f, and the normalized variance η <dig>  where η <dig> = var/⟨i⟩ <dig>  and i = i^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@ or iij for simulation and data respectively.
then  and  become
〈i^j〉=f
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeyykjeuafmysakkbakaadawgaawcbagaemoaaogabeaakiabgqyixlabg2da9iabbaeagjabcufabnaatafabagaemysak0aasbaasqaaiabdmgapjabdqgaqbqabaaabagaemyaakgabeqdcqwiqisvaogaeiyxa0faaa@3e08@
η2≤η <dig>
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcuwgjbqsgaqcamaabaaaleaacqwgqbgaaeqaaogaeiykakiaeyizimqae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakdaweqbqaaiabdmeajnaabaaaleaacqwgpbqacqwgqbgaaeqaaaqaaiabdmgapbqab0gaesokiufakiabcmcapiabc6cauaaa@4195@
in , ⟨i^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@⟩ = ma' exp from simulation and f=aexp⁡
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeeoraykaei4was1aambuaeaacqwgjbqsdawgaawcbagaemyaakmaemoaaogabeaaaeaacqwgpbqaaeqaniablqiivbgccqggdbqxcqgh9aqpcqwghbqycygglbqzcqgg4baecqggwbaccqggoaakcqghsislcqwgqbgacqggvawliigacqwf7oabcqggpaqkaaa@43e3@ from the ensemble data.
because λ' is only determined by d and ω, we first select combinations of d and ω such that λ' = λ. selecting a synthesis rate j determines a', and also m, because m = a/a'.
finally, biologically reasonable values of j are determined by the constraint .
we seek to establish which terms of equation  dominate the observed variance in different parts of the embryo.
to do this, we will graphically compare the total observed variance with a set of simulations in which variance arises from different subsets of the random variables in .
we denote these restricted models by placing brackets around the subset of random variables which contribute to the variance.
thus the full model in  can be denoted by i^j=
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9iabcufabnaahaaabaaccigae8xsdemaemoba4mae8nsdigacaglcmaacqggdbqxaaa@3780@.
a model with no variance contributed by background is denoted by =α^n^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4was1aaecaaeaaiigacqwfxoqycqwgubgbaiaawkwaaiabc2fadjabg2da9iqb=f7ahzaajagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@37ce@, while models in which all variance is contributed only by rescaling, background, or molecular number are denoted respectively by =α^〈n^j〉
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4wasfccigaf8xsdembakaacqggdbqxcqgh9aqpcuwfxoqygaqcaiabgmyihlqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqghqms8aaa@3944@, =m〈n^j〉+β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4wasfccigaf8nsdimbakaacqggdbqxcqgh9aqpcqwgtbqbcqghpms4cuwgubgbgaqcamaabaaaleaacqwgqbgaaeqaaogaeyokjevaey4kasiaf8nsdimbakaaaaa@3b8d@, and =mn^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4waslafmoba4mbakaacqggdbqxcqgh9aqpcqwgtbqbcuwgubgbgaqcamaabaaaleaacqwgqbgaaeqaaaaa@352f@.
using this notation, a comparison of the definition of η <dig> with  shows that we can write
η2=η2+η2=η2+σα2var⁡)/〈n^j〉2+η <dig>
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaqbaeaabiwaaaqaaggaciab=d7aonaacaaaleqabagaegomaidaaogaeiikagiafmysakkbakaadawgaawcbagaemoaaogabeaakiabcmcapaqaaiabg2da9aqaaiab=d7aonaacaaaleqabagaegomaidaaogaeiikagiaei4was1aaecaaeaacqwfxoqycqwgubgbaiaawkwaaiabc2fadjabcmcapiabgucariab=d7aonaacaaaleqabagaegomaidaaogaeiikagiaei4waslaf8nsdimbakaacqggdbqxcqggpaqkaeaaaeaacqgh9aqpaeaacqwf3oaadaahaawcbeqaaiabikdayaaakiabcicaoiabcufabjqbd6gauzaajagaeiyxa0laeiykakiaey4kasiae83wdm3aa0baasqaaiab=f7ahbqaaiabikdayaaakigbcaha2jabcggahjabckhayjabcicaoiqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccuwgobgtgaqcaiabcicaoiabicdawiabcycasiabigdaxiabcmcapiabcmcapiabc+caviabgmyihlqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqghqms8daahaawcbeqaaiabikdayaaakiabgucariab=d7aonaacaaaleqabagaegomaidaaogaeiikagiaei4waslaf8nsdimbakaacqggdbqxcqggpaqkcqgguaglaaaaaa@7951@
the middle term of the above equation simplifies even further when the molecular number is large so that n^j=〈n^j〉
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaakiabg2da9iabgmyihlqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqghqms8aaa@366c@ and hence σα2var⁡)=σα2
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83wdm3aa0baasqaaiab=f7ahbqaaiabikdayaaakigbcaha2jabcggahjabckhayjabcicaoiqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccuwgobgtgaqcaiabcicaoiabicdawiabcycasiabigdaxiabcmcapiabcmcapiabg2da9iab=n8aznaadaaaleaacqwfxoqyaeaacqaiyagmaaaaaa@4468@.
if the contribution of background noise β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ is also small then it is the case that
η2=σα2=η <dig>
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcuwgjbqsgaqcamaabaaaleaacqwgqbgaaeqaaogaeiykakiaeyypa0jae83wdm3aa0baasqaaiab=f7ahbqaaiabikdayaaakiabg2da9iab=d7aonaacaaaleqabagaegomaidaaogaeiikagiaei4waslaf8xsdembakaacqggdbqxcqggpaqkcqgguaglaaa@438a@
in summary, when bcd molecular number is large in the anterior region of the embryo, we expect to see a constant level of normalized variance in our fluorescence data, contributed solely by rescaling noise α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@.
thus, α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ can be identified independently from m and other random variables in our data.
physical constraints from a high-variance ensemble of embryos
we first show results of the above statistical comparison using all  <dig> flyex cycle  <dig> embryos.
this ensemble contains  <dig> nuclei, with about  <dig> nuclei per bin.
in figure 2a we see that the normalized variance of the ensemble data η2
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakdaweqbqaaiabdmeajnaabaaaleaacqwgpbqacqwgqbgaaeqaaaqaaiabdmgapbqab0gaesokiufakiabcmcapaaa@376e@ asymptotically approaches the simulation curve η2=σα2
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwcuwfxoqygaqcaiabc2fadjabcmcapiabg2da9iab=n8aznaadaaaleaacqwfxoqyaeaacqaiyagmaaaaaa@3a08@ in the anterior region of the embryo.
σα values that are too high will place the black model curve above the data points on the left, and this constrains σα to be less than  <dig> .
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakdaweqbqaaiabdmeajnaabaaaleaacqwgpbqacqwgqbgaaeqaaaqaaiabdmgapbqab0gaesokiufakiabcmcapaaa@376e@ from the high-variance ensemble of  <dig> embryos.
lines denote simulation results as shown in the key.
denotes the full model η2=η2
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcuwgjbqsgaqcamaabaaaleaacqwgqbgaaeqaaogaeiykakiaeyypa0jae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwdaqiaaqaaiab=f7ahjabd6gaujab=j7aibgaaypadagaeiyxa0laeiykakcaaa@4084@,  denotes η2
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwcuwfxoqygaqcaiabc2fadjabcmcapaaa@348b@,  denotes η2
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwcuwgubgbgaqcaiabc2fadjabcmcapaaa@3456@, and  denotes η2
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwcuwfyogygaqcaiabc2fadjabcmcapaaa@348d@.
the parameters used in simulation were σα =  <dig> , σβ =  <dig> , j =  <dig> , m =  <dig> , d =  <dig>   and ω =  <dig>  .
is a scatterplot of bcd molecular number n^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@, while the inset shows simulated fluorescence intensity for models and parameters used in panel a.  shows the residuals  of panel b. ,  and  show the same information as panel a, b, and c respectively, but from the low-variance ensemble of  <dig> embryos, with parameters used in simulation σα =  <dig> , σβ, =  <dig> , j =  <dig> , m =  <dig> , d =  <dig>   and ω =  <dig>  .
note that the axes in panel a and d are scaled differently, and the absolute molecule number is shown in the inset to panel e.
with regard to the molecular parameters of diffusion, adopting the value of d =  <dig>   given in the literature  <cit>  together with our observed value of λ =  <dig>   for this ensemble implies that ω =  <dig>   in order to satisfy .
these values, for any j, will cause the mean molecular number gradient, 〈n^j〉
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeyykjeuafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgqyixdaa@325e@, to be in steady state after about  <dig> seconds .
a lower bound of j >  <dig>  and an upper bound of m < <dig>  is imposed by the experimentally observed magnitude of η <dig>  finally, we estimate the constraint for background noise β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ by satisfying  in the posterior end of the embryo.
violation of the inequality  would cause the black model curve to be above the data on the right hand side of figure 2a, and hence we have an upper bound σβ <  <dig> .
in summary, analysis of our high-variance ensemble of embryos dataset implies that the bcd synthesis rate is higher than  <dig> .
figure 2b shows a scatterplot of the molecular number associated with this synthesis rate.
the panel indicates that most nuclei in the anterior fifth of the embryo contain more than  <dig> molecules of bcd after reaching steady state, and that these molecular numbers can fluctuate over a range of more than  <dig> molecules in this region .
this panel shows that these molecular fluctuations, even at the largest level of normalized variance compatible with data, still do not account for the observed variance in experimental observations.
the additional variance comes chiefly from rescaling noise α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@.
background noise β^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ only has a significant effect at the posterior end of the embryo, and indeed dominates the normalized variance in that region .
towards the posterior, η2
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwcuwfyogygaqcaiabc2fadjabcmcapaaa@348d@ rises faster than η2
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcuwgubgbgaqcamaabaaaleaacqwgqbgaaeqaaogaeiykakcaaa@3369@, and its sharp rise may in certain cases serve as a marker to distinguish regimes dominated by molecular noise from those dominated by background noise.
physical constraints from a low-variance ensemble of embryos
in the ensemble of embryos discussed above, a portion of the variance observed is likely to stem from embryo to embryo variation in staining conditions and inherent biological parameters.
in order to find a better upper limit for observed molecular fluctuations, it is desirable to analyze a set of embryos whose properties are as uniform as possible.
in order to generate such a set, we considered the value of λi in the exponential fit f  = ai exp to each embryo.
we then grouped embryos according to common values of 1/λi, taken to two decimal places.
in the largest such group, which contained  <dig> members, we rescaled the data to a common amplitude by letting
i′ij=〈ai〉aiiij.
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaqcfaoafmysakkbauaadawgaaqaaiabdmgapjabdqgaqbqabagaeyypa0zaasaaaeaacqghpms4cqwghbqydawgaaqaaiabdmgapbqabagaeyokjepabagaemyyae2aasbaaeaacqwgpbqaaeqaaaaacqwgjbqsdawgaaqaaiabdmgapjabdqgaqbqabagaeiola4caaa@3fae@
these  <dig> processed embryos constitute an ensemble ∪ii′ij
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaambuaeaacuwgjbqsgaqbamaabaaaleaacqwgpbqacqwgqbgaaeqaaaqaaiabdmgapbqab0gaesokiufaaaa@32e2@ for statistical analysis as described in the previous section.
a trade-off of this treatment is the loss of statistical sample size, with only around  <dig> nuclei in each bin.
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ is too small to be separated from n^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@.
we compared this data to simulations performed using the same parameter values reported in the previous section, except that λ  had a value of  <dig>  , leading to ω taking on a value of  <dig>  .
these values cause the deterministic system to relax to steady state after  <dig> seconds, as was the case with the high variance ensemble.
in the high variance ensemble , the choice of σβ was dictated by the necessity of matching the observed variance along the entire a-p axis.
in the case of the low variance ensemble this is not required, but measurements from nonexpressing nuclei indicate that σβ is equal to about  <dig>  it is important to have at least a rough estimate for σβ so that fluctuations from this source are not spuriously assigned to fluctuations in molecular number.
finally, these constraints require that the lower bound of synthesis rate j be  <dig>  in order that that η2≃η2
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcuwgjbqsgaqcamaabaaaleaacqwgqbgaaeqaaogaeiykakiaes4qisjae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakdaweqbqaaiqbdmeajzaafawaasbaasqaaiabdmgapjabdqgaqbqabaaabagaemyaakgabeqdcqwiqisvaogaeiykakcaaa@3feb@.
the corresponding upper bound of molecule-to-fluorescence rescaling ratio is m =  <dig> .
the lower limit of j imposed by this ensemble of low variance embryos in turn implies that there must be more than  <dig> bcd molecules per subvolume in the middle of the embryo  after reaching steady state .
moreover, it implies that the bcd molecular gradient does not drop to  <dig> as the fluorescence intensity reaches the presumed background level, but remains at a level of at least  <dig> molecules per subvolume at j =  <dig>  and  <dig> molecules per subvolume at the posterior pole .
note that the variance of fluorescence measurements is similar between the two embryo ensembles , but that the portion of that variance assigned to fluctuations in molecular number is smaller for the ensemble of  <dig> embryos .
the tighter constraints from the smaller ensemble make the lower limits on molecular number higher , and the lower limit on the variance of molecular number higher , although the higher limit of the normalized variance becomes smaller .
the bcd molecular number will thus vary by more than  <dig> molecules in the middle of the embryo.
the 13% rescaling noise α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ is still the main source of the characteristic variation observed in the anterior region of our fluorescence intensity data.
noise strength
in most applications the most important measure of fluctuation is the normalized variance η <dig>  <cit> .
a different quantity, known as the fano factor or noise strength  <cit> , has been used by some authors as a marker to distinguish different stochastic mechanisms.
the fano factor ν is given by ν=var⁡/〈i^j〉
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae8xvd4maeiikagiafmysakkbakaadawgaawcbagaemoaaogabeaakiabcmcapiabg2da9igbcaha2jabcggahjabckhayjabcicaoiqbdmeajzaajawaasbaasqaaiabdqgaqbqabagccqggpaqkcqggvawlcqghpms4cuwgjbqsgaqcamaabaaaleaacqwgqbgaaeqaaogaeyokjepaaa@42cb@, where ν =  <dig> in a poisson process.
all stochastic simulations of bcd intrinsic noise n^j
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ give ν =  <dig>  as shown in figure 3a.
by contrast, the full statistical model for either of the two ensembles of embryos examined  is obviously non-poisson, not only because ν ≠  <dig> but also because in the data, ν is a function of the mean.
this happens because at larger values of molecular number, the variance of α^
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ has a dominant role, even when its value is small.
even if we model our data without rescaling noise using the
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4waslafmoba4mbakaacqggdbqxaaa@2fc8@ random variable alone, uncertainty in the value of the rescaling constant m itself leads to ambiguity in the observed value of the fano factor .
this is a natural consequence of the dimensions of the fano factor.
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagagiodaynaeiyyaemaeiocainaeiikagiafmoba4mbakaadawgaawcbagaemoaaogabeaakiabcmcapaaa@34b9@ divided by molecular mean 〈n^j〉
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeyykjeuafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgqyixdaa@325e@.
the key indicates that parameters were obtained from  the high-variance ensemble of  <dig> embryos,  the low-variance ensemble of  <dig> embryos, and  the extreme condition of high diffusion rate d =  <dig> , decay rate ω =  <dig>   and synthesis rate j =  <dig> .
as shown in the key, noise strength of simulated fluorescence intensity i^j=
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9iabcufabnaahaaabaaccigae8xsdemaemoba4mae8nsdigacaglcmaacqggdbqxaaa@3780@ and rescaled gradient
mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4waslafmoba4mbakaacqggdbqxaaa@2fc8@ were obtained using parameters from  the high-variance ensemble of  <dig> embryos and  the low-variance ensemble of  <dig> embryos.
CONCLUSIONS
we have compared the nucleus to nucleus variation in expression levels of the exponentially distributed bcd gradient observed in fixed tissue in a steady state with a stochastic model of the diffusion equation.
the model is well supported, in the sense that there is a well-supported physical model for the spatial dependence of mean concentrations of bcd  <cit>  on the scale of the embryo.
the first major result of our analysis is to note that in many individual embryos the nucleus to nucleus variation in the log of concentration is independent of spatial position.
this pattern of variation, which amounts to multiplicative noise in concentration space, is completely incompatible with the stochastic behavior of the diffusion equation.
simulations of the diffusion equation over an exhaustively large region of parameter space without exception give rise to solutions in which nucleus to nucleus variation of the bcd gradient is a function of position in the embryo, whether this variation is measured directly in bcd levels or in their logarithms.
the data which we compare the model to is in the form of fluorescence levels, not concentrations.
although there is now good evidence that the specific batch of serum used to obtain this data has a mean response to bcd  <cit>  which is linear, there is no quantitative information about the variance of this sensitivity.
previous work on intrinsic molecular noise in yeast and bacteria utilized gfp  <cit> in vivo, a situation where fluorescence is detected without molecular amplification.
in the data reported here, and in most studies with fixed tissue, the signal from bound primary antibodies is amplified by incubation with secondary antibodies conjugated to a fluorescent dye.
it is easily imaginable that this amplification process itself is subject to molecular fluctuations.
these fluctuations would then give rise to rescaling noise in the proportionality between fluorescence levels and primary antigen molecular number from nucleus to nucleus.
we have shown that such variation can explain the multiplicative noise observed.
when data from fixed embryos are pooled, better statistics are obtained.
we analyzed data pooled from the entire dataset  as well as a smaller pool of data from a set of embryos selected to have nearly identical mean bcd profiles .
in order to analyze these data, we constructed an explicit statistical model.
the model considers three sources of variance: intrinsic noise, rescaling noise, and background noise.
by means of this statistical model it is possible to separate, at least roughly, the contributions of different sources of fluctuation.
we have discussed mechanisms for the first two of these; the third, background noise, represents small fluctuations in the quantity of nonspecific molecules  from nucleus to nucleus.
because mean background has been previously removed from this data  <cit> , the background noise has a mean of zero.
we have confirmed that the background removal method does not affect our results, and the mean background intensity before removal is independent from the bcd molecule-to-fluorescence rescaling ratio .
we also found from non-expressing areas of our data that the background noise  has about 54% positive correlation to the mean background intensity.
the results of this analysis constrained the physical parameters of the stochastic model considerably, with sharper constraints provided by the smaller dataset.
the data require that the synthesis rate, j, of bcd from its pool of anteriorly deposited mrna be greater than  <dig> .
we chose the subvolumes of the model to have the same volume as a nucleus, and hence the constraint on j also implies that there are a mean of at least  <dig> molecules of bcd per nucleus at the midpoint of the embryo, and mean levels of at least  <dig> molecules of bcd per nucleus at the posterior pole.
in terms of concentration, our results show that bcd concentration at the midpoint of the embryo is greater than  <dig> nm.
recently, a direct in vivo measurement of bcd concentration was performed and yielded a mid embryo bcd concentration of  <dig> nm  <cit> , fully compatible with our results.
although we are able to extract a clear signature of intrinsic molecular noise from the data by means of the statistical model, we also showed that at least one quantity diagnostic for different stochastic mechanisms, the fano factor, cannot be read out from the data.
although the fano factor v =  <dig> in the simulations, the full statistical model gives rise to a fano factor which is a function of the mean, and even if all noise is restricted to be intrinsic, the observed fano factor depends on the scale for conversion from fluorescence to molecular number.
we believe that our results in general demonstrate that fixed material processed with secondary fluorophores is not well suited to studies of molecular fluctuations.
this arises from three issues, which may be separable.
fixation obviously prevents repeated observations on the same cell.
while that is clearly a limitation, it need not affect an investigation of a molecule whose mean values are in steady state.
the other two issues concern amplification.
gfp is intrinsically fluorescent, but antibodies must bind to antigen, a process that is in itself subject to intrinsic molecular fluctuations.
in the present study, the level of molecular fluctuation is doubtless increased by the need to bind secondary antibodies conjugated to fluorophore to the already bound primary antibody.
it is thus possible that better data from fixed tissue could be obtained by conjugating dye directly to the primary antibody.
this may prove to be a useful measure in situations where constructing a gfp fusion that can functionally substitute for the native gene is difficult or impossible.
more generally, we suggest that it is important to know how precision and robustness of developmental control is achieved at the molecular number level throughout development.
indeed, there is little doubt that stochastic processes are important later in development.
adult drosophila normally have four scutellar bristles.
in mutants of the scute gene, the number of bristles varies between one and three  <cit> , strongly suggesting a stochastic process.
on a more theoretical level, it has been suggested that fluctuations can augment the operating capabilities of biological regulatory networks  <cit> .
our results indicate that in vivo monitoring of gene expression will be required to obtain high quality data on stochastic gene expression phenomena in eucaryotes.
the central technical problem that must be solved to conduct such studies is the complete replacement of the endogenous gene with a fluorescently tagged functional version  <cit> .

 recently considerable attention has been given to the characterization and understanding of intrinsic molecular noise in biological systems. nearly all of these studies were performed using in vivo fluorescent reporters in single cell systems. in multicellular organisms, however, most quantitative gene expression data are obtained from fixed tissues. examples of such data for the drosophila segmentation system are contained in the flyex database, which provides spatiotemporal data on the expression of developmental segmentation genes. these data on protein expression levels are at cellular resolution and were obtained by means of immunofluorescence histochemistry and confocal scanning microscopy. at a large spatial scale, expression levels in these embryos form expression domains characteristic for each gene, but smaller fluctuations in expression levels between adjacent nuclei appear random. in this paper, we investigate the question of whether these fluctuations are a consequence of intrinsic molecular noise or stem from some type of measurement uncertainty. these alternatives are, of course, not mutually exclusive. a complicating factor in separating the above alternatives is that each one involves an unknown chemical mechanism. intrinsic noise will have a major contribution from fluctuations in the rate of initiation of transcription, but the chemical mechanisms underlying this process in multicellular organisms are very poorly understood. measurement uncertainty can stem from chemical causes such as fluctuations in the number of primary and secondary antibody molecules which bind to proteins in the fixed embryo. the chemistry of this process is also very poorly understood. if observations could be made on a process whose fluctuation properties could be reliably predicted by a numerical model, comparison of the predicted fluctuations with those observed will provide critical information for distinguishing whether observed nucleus-to-nucleus variations are a consequence of intrinsic biological noise or merely fluctuations arising from the staining procedure. an excellent candidate for a process with predictable fluctuation properties is one that involves only diffusion and first order decay. there is good evidence that the formation of the protein gradient of the morphogen bicoid takes place by means of these two processes. bcd protein is distributed in an exponential profile along the anterior-posterior axis with higher concentrations towards the anterior. this gradient forms by translation of maternally deposited mrna at the anterior pole of the embryo, and the synthesized protein spreads through the syncytial embryo by diffusion accompanied by decay. the observed exponential profile corresponds to a solution of fick's equation for a substance undergoing first order decay and diffusing from a point source in one dimension, and hence it is reasonable to suppose that the reaction-diffusion mechanism leading to the formation of the gradient is reasonably well understood. quantitative observations of this gradient in fixed tissue exhibit small fluctuations between neighboring nuclei, while the overall exponential profile ensures that such fluctuations can be monitored over a wide range of concentrations which extend to the lower limits of detectability. the intrinsic fluctuations of the bcd gradient will be well described by a stochastic reaction-diffusion master equation. such equations typically do not have analytic solutions and are usually solved by running repeated stochastic simulations. a well known simulation algorithm due to gillespie performs an exact simulation of the chemical master equation for a well mixed system. this method has been extended to spatially distributed systems by elf and others. these authors divide the spatially extended system into a series of subvolumes that are small enough to be regarded as well mixed. in each subvolume chemical reactions are simulated by gillespie's algorithm, while diffusion between subvolumes is treated as a first order reaction. in this study, we compare the results of such simulations to data in order to discover whether or not the data is sufficiently accurate as to exhibit the signature of a simple stochastic process. stochastic processes underlying biological regulation can in general form complex patterns as a result of the reaction network, and for this reason a full consideration of dimensional geometry is often necessary. in the case considered here fluctuations occur passively in the course of diffusion and the statistical signature that we seek is independent of detailed geometry. for these reasons, we chose to model a dimensional system.  in the following, we first characterize the statistical properties of random variations in expression level between adjacent nuclei of individual embryos and compare them to the results of stochastic simulations. at the level of single embryos, we do not find a clear signature of stochastic processes in the data, but the need to separate spatially changing mean expression values from their variation limits the amount of statistical information that can be obtained from individual embryos. to address this problem, we consider ensembles of embryos, both over the whole dataset and over a restricted subset with low embryo-to-embryo variation. in order to interpret these data, we introduce a stochastic model of the immunostaining procedure. we denote all random variables in this article with an upper hat but write a specific value without the hat.thus, for example, n^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ denotes a random variable for the number of bcd molecules in subvolume j, while nj denotes a particular value of this variable. single embryo analysis in the logarithmic domain we find that the bcd profiles of individual embryos observed by immunostaining are exponential, as previously reported. this profile strongly supports the model of an effective fickian diffusion which gives an exponential decay solution at steady state. each embryo contains a collection of observations of expression ij, and the observed properties of the exponential gradient will depend on the data through a function f embodying the least squares fitting procedure such that f = a exp. the residuals of the exponential vary from embryo to embryo in our data. in logarithmic coordinates the size of the residuals is independent of j for the anterior portion of the embryo, and in certain embryos the size of the residuals is completely independent of j. for these embryos, the residuals are well described by ln⁡=ln⁡−j/λ+w_, mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagagiibawmaeioba4maeiikagiaemysak0aasbaasqaaiabdqgaqbqabagccqggpaqkcqgh9aqpcyggsbabcqggubgbcqggoaakcqwghbqycqggpaqkcqghsislcqwgqbgacqggvawliigacqwf7oabcqghrawkdaqiaaqaaiabdefaxbgaaypadagaeiilawcaaa@42b6@ where w_ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaaecaaeaacqwgxbwvaiaawkwaaaaa@2dcc@ is an unknown random variable independent of j, so that the variance of ln is equal to the variance of w_ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaaecaaeaacqwgxbwvaiaawkwaaaaa@2dcc@. intrinsic noise is insufficient in order to understand whether the observed characteristic variation is a consequence of intrinsic fluctuations in bcd molecular number, we performed stochastic simulations of an rdme which describes the time evolution of the bcd gradient and compared them to data. we imagine the observed intensity ij to be a particular value of the random variable i^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@. we further suppose that the random variable i^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@ is determined by a direct linear rescaling of the bcd molecular number such that i^j=mn^j, mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9iabd2gatjqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqggsaalaaa@3530@ where the factor m represents the proportionality between one bcd molecule and the corresponding fluorescence intensity, which includes the combined effects of tissue fixation, first and second antibody binding, fluorescence excitation and image processing. as we do not know the exact in vivo system parameters and molecular number within each nucleus, we performed a complete inspection of the behavior of the variance of mn^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ in the four dimensional parameter space f, where j is the synthesis rate of bcd in the anterior compartment, ω is its decay rate, and d is the diffusion coefficient. it was always true that the residuals in logarithmic coordinate increased towards the posterior of the embryo, or in other words at lower levels of bcd. this behavior strongly contrasts with the position-independent residuals seen in figure and described in equation. measurement rescaling noise dominates the simulation results suggest that intrinsic noise cannot explain the pattern of variance seen in figure another possibility is the measurement process itself. in order to analyze this process, we consider a simple model of the measurement of fluorescence intensity, where i^j=α^n^j+β^. mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9ggaciqb=f7ahzaajagafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgucariqb=j7aizaajaacbagae4nla4caaa@381a@ here α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ is a spatially uniform random variable which replaces m in equation, and β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ is a spatially uniform random variable which represents nonspecific background staining. this picture allows us to consider noise that arises from both intrinsic and measurement related sources. the simplest way to understand the consequences of is to imagine the consequences if only one of α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@, n^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@, and β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ are allowed to have finite variance while the other two are constrained to deterministic behavior. we have already pointed out that allowing finite variance for n^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ with α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ and β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ deterministic leads to an increase of variance towards the posterior. the same is true if all noise comes from β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@, that is if noise from background staining dominates. this will also lead to more noise in the logarithmic domain towards the posterior as β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ provides a larger proportion of the total detected signal. let us next consider the case where all noise comes from α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@. in order to understand the role of α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@, note that the spatial pattern of variance observed in figure can be captured by an exponential function multiplied by a normal random variable. this suggests that the simplest picture for is given by assuming that there is no background noise β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ and no intrinsic noise for bcd so that n^j=〈n^j〉 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaakiabg2da9iabgmyihlqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqghqms8aaa@366c@. moreover, it suggests that the rescaling noise α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ from measurement uncertainty is uniform across the embryo and is normally distributed with α^=m), mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaacqgh9aqpcqwgtbqbcqggoaakcqaixaqmcqghrawkcqwfdpwcdawgaawcbagae8xsdegabeaakiqbd6eaozaajagaeiikagiaegimaajaeiilawiaegymaejaeiykakiaeiykakiaeiilawcaaa@3dd7@ where n^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmota4kbakaacqggoaakcqaiwaamcqggsaalcqaixaqmcqggpaqkaaa@3178@ is a normal independent random variable with mean and variance then we can model i^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@ by i^j=α^〈n^j〉. mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9ggaciqb=f7ahzaajagaeyykjeuafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgqyixlabc6cauaaa@390a@ in steady state, 〈n^j〉 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeyykjeuafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgqyixdaa@325e@ = a/m exp. taking logarithms allows us to write ln⁡=ln⁡−j/λ+ln⁡). mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagagiibawmaeioba4maeiikagiafmysakkbakaadawgaawcbagaemoaaogabeaakiabcmcapiabg2da9igbcygasjabc6gaujabcicaoiabdggahjabcmcapiabgkhitiabdqgaqjabc+cavggaciab=t7asjabgucarigbcygasjabc6gaujabcicaoiabigdaxiabgucariab=n8aznaabaaaleaacqwfxoqyaeqaaogafmota4kbakaacqggoaakcqaiwaamcqggsaalcqaixaqmcqggpaqkcqggpaqkcqgguaglaaa@504e@ comparison with equation indicates that w_ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaaecaaeaacqwgxbwvaiaawkwaaaaa@2dcc@ = ln mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmota4kbakaacqggoaakcqaiwaamcqggsaalcqaixaqmcqggpaqkaaa@3178@). the reason we do not see intrinsic noise in this individual embryo is most likely low measurement sensitivity, that is to say a low molecule-to-fluorescence mean rescaling ratio m. when m is small enough, it can mask the variance of n^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ because the variance of the rescaled gradient is given by var⁡=m2var⁡. mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagagiodaynaeiyyaemaeiocainaeiikagiaemyba0mafmoba4mbakaadawgaawcbagaemoaaogabeaakiabcmcapiabg2da9iabd2gatnaacaaaleqabagaegomaidaaogagiodaynaeiyyaemaeiocainaeiikagiafmoba4mbakaadawgaawcbagaemoaaogabeaakiabcmcapiabc6cauaaa@43c6@ as m → with var bounded in a reasonable way throughout the embryo, we will get var → hence the rescaled gradient mn^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ can be treated as deterministic by letting mn^j=〈mn^j〉=m〈n^j〉. mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaemyba0mafmoba4mbakaadawgaawcbagaemoaaogabeaakiabg2da9iabgmyihlabd2gatjqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqghqms8cqgh9aqpcqwgtbqbcqghpms4cuwgubgbgaqcamaabaaaleaacqwgqbgaaeqaaogaeyokjevaeiola4caaa@4358@ in summary, we suspect the nucleus-to-nucleus variation observed in our data comes chiefly from the experimental rescaling noise α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@, which is normally distributed. if bcd intrinsic noise is to be observed, then the fluorescence noise intensity should be a function of the mean intensity in logarithm, instead of a constant as observed in embryo ms nevertheless, the necessity of considering data in spatially resolved bins limits the amount of information that can be obtained from a single embryo. more information can be obtained by pooling data from many embryos, and we discuss this point in the next section. statistical analysis of an ensemble of embryos statistical analysis of many embryos is required in order to take our analysis further. this analysis will show how physical constraints on the model can be inferred from the ensemble dataset, and independent random variables separated. we consider a set of embryos indexed by i with expression levels iij, and then pool data from corresponding bins to obtain the ensemble dataset ∪iiij mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaambuaeaacqwgjbqsdawgaawcbagaemyaakmaemoaaogabeaaaeaacqwgpbqaaeqaniablqiivbaaaa@32d6@. since this dataset includes embryo-to-embryo variability, i.e. the variation of system parameters and experimental conditions over different embryos, the variance of the ensemble profile will be an upper bound for the average variance within each embryo. this allows us to identify physical constraints for system parameters and to determine if the model behaves properly within the permitted range of parameters. we define independent global random variables i^j=α^n^j+β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9ggaciqb=f7ahzaajagafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgucariqb=j7aizaajaaaaa@36e3@ as described in the last section with normal distributed measurement uncertainty α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ = m mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmota4kbakaacqggoaakcqaiwaamcqggsaalcqaixaqmcqggpaqkaaa@3178@). we also now assume that background noise is normally distributed with β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ = σβ n^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmota4kbakaacqggoaakcqaiwaamcqggsaalcqaixaqmcqggpaqkaaa@3178@. we assume such global random variable represent the average variability for each embryo. the statistics of simulated global random variables were then collected from stochastic simulation runs and the bcd molecular number random variable n^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ were sampled after reaching steady state. in this section we explicitly consider the effects of different choices for the molecule-to-fluorescence rescaling ratio m by comparing simulations to data at differing values of this parameter. because m is not an explicit input to the model, this comparison is effected by varying the synthesis rate j, which varies the molecular number, and comparing the behavior of the model to fluorescence data which is on a fixed but arbitrary scale. we seek optimal values of m such that the simulated global random variable i^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@ is constrained by the variation observed in the immunostained ensemble data ∪iiij mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaambuaeaacqwgjbqsdawgaawcbagaemyaakmaemoaaogabeaaaeaacqwgpbqaaeqaniablqiivbaaaa@32d6@ by the condition 〈i^j〉=〈∪iiij〉 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeyykjeuafmysakkbakaadawgaawcbagaemoaaogabeaakiabgqyixlabg2da9iabgmyihpaatafabagaemysak0aasbaasqaaiabdmgapjabdqgaqbqabaaabagaemyaakgabeqdcqwiqisvaogaeyokjepaaa@3df8@ var⁡≤var⁡. mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagagiodaynaeiyyaemaeiocainaeiikagiafmysakkbakaadawgaawcbagaemoaaogabeaakiabcmcapiabgsmijkgbcaha2jabcggahjabckhayjabcicaomaatafabagaemysak0aasbaasqaaiabdmgapjabdqgaqbqabaaabagaemyaakgabeqdcqwiqisvaogaeiykakiaeiola4caaa@4441@ because comparison with mean and variance of the ensemble data is not straightforward, we approximate the above conditions by comparing their exponential fit, f, and the normalized variance η where η = var/⟨i⟩ and i = i^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@ or iij for simulation and data respectively. then and become 〈i^j〉=f mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeyykjeuafmysakkbakaadawgaawcbagaemoaaogabeaakiabgqyixlabg2da9iabbaeagjabcufabnaatafabagaemysak0aasbaasqaaiabdmgapjabdqgaqbqabaaabagaemyaakgabeqdcqwiqisvaogaeiyxa0faaa@3e08@ η2≤η mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcuwgjbqsgaqcamaabaaaleaacqwgqbgaaeqaaogaeiykakiaeyizimqae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakdaweqbqaaiabdmeajnaabaaaleaacqwgpbqacqwgqbgaaeqaaaqaaiabdmgapbqab0gaesokiufakiabcmcapiabc6cauaaa@4195@ in, ⟨i^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaaaaa@2e87@⟩ = ma' exp from simulation and f=aexp⁡ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeeoraykaei4was1aambuaeaacqwgjbqsdawgaawcbagaemyaakmaemoaaogabeaaaeaacqwgpbqaaeqaniablqiivbgccqggdbqxcqgh9aqpcqwghbqycygglbqzcqgg4baecqggwbaccqggoaakcqghsislcqwgqbgacqggvawliigacqwf7oabcqggpaqkaaa@43e3@ from the ensemble data. because λ' is only determined by d and ω, we first select combinations of d and ω such that λ' = λ. selecting a synthesis rate j determines a', and also m, because m = a/a'. finally, biologically reasonable values of j are determined by the constraint. we seek to establish which terms of equation dominate the observed variance in different parts of the embryo. to do this, we will graphically compare the total observed variance with a set of simulations in which variance arises from different subsets of the random variables in. we denote these restricted models by placing brackets around the subset of random variables which contribute to the variance. thus the full model in can be denoted by i^j= mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9iabcufabnaahaaabaaccigae8xsdemaemoba4mae8nsdigacaglcmaacqggdbqxaaa@3780@. a model with no variance contributed by background is denoted by =α^n^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4was1aaecaaeaaiigacqwfxoqycqwgubgbaiaawkwaaiabc2fadjabg2da9iqb=f7ahzaajagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@37ce@, while models in which all variance is contributed only by rescaling, background, or molecular number are denoted respectively by =α^〈n^j〉 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4wasfccigaf8xsdembakaacqggdbqxcqgh9aqpcuwfxoqygaqcaiabgmyihlqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqghqms8aaa@3944@, =m〈n^j〉+β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4wasfccigaf8nsdimbakaacqggdbqxcqgh9aqpcqwgtbqbcqghpms4cuwgubgbgaqcamaabaaaleaacqwgqbgaaeqaaogaeyokjevaey4kasiaf8nsdimbakaaaaa@3b8d@, and =mn^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4waslafmoba4mbakaacqggdbqxcqgh9aqpcqwgtbqbcuwgubgbgaqcamaabaaaleaacqwgqbgaaeqaaaaa@352f@. using this notation, a comparison of the definition of η with shows that we can write η2=η2+η2=η2+σα2var⁡)/〈n^j〉2+η mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaqbaeaabiwaaaqaaggaciab=d7aonaacaaaleqabagaegomaidaaogaeiikagiafmysakkbakaadawgaawcbagaemoaaogabeaakiabcmcapaqaaiabg2da9aqaaiab=d7aonaacaaaleqabagaegomaidaaogaeiikagiaei4was1aaecaaeaacqwfxoqycqwgubgbaiaawkwaaiabc2fadjabcmcapiabgucariab=d7aonaacaaaleqabagaegomaidaaogaeiikagiaei4waslaf8nsdimbakaacqggdbqxcqggpaqkaeaaaeaacqgh9aqpaeaacqwf3oaadaahaawcbeqaaiabikdayaaakiabcicaoiabcufabjqbd6gauzaajagaeiyxa0laeiykakiaey4kasiae83wdm3aa0baasqaaiab=f7ahbqaaiabikdayaaakigbcaha2jabcggahjabckhayjabcicaoiqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccuwgobgtgaqcaiabcicaoiabicdawiabcycasiabigdaxiabcmcapiabcmcapiabc+caviabgmyihlqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqghqms8daahaawcbeqaaiabikdayaaakiabgucariab=d7aonaacaaaleqabagaegomaidaaogaeiikagiaei4waslaf8nsdimbakaacqggdbqxcqggpaqkcqgguaglaaaaaa@7951@ the middle term of the above equation simplifies even further when the molecular number is large so that n^j=〈n^j〉 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaakiabg2da9iabgmyihlqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccqghqms8aaa@366c@ and hence σα2var⁡)=σα2 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83wdm3aa0baasqaaiab=f7ahbqaaiabikdayaaakigbcaha2jabcggahjabckhayjabcicaoiqbd6gauzaajawaasbaasqaaiabdqgaqbqabagccuwgobgtgaqcaiabcicaoiabicdawiabcycasiabigdaxiabcmcapiabcmcapiabg2da9iab=n8aznaadaaaleaacqwfxoqyaeaacqaiyagmaaaaaa@4468@. if the contribution of background noise β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ is also small then it is the case that η2=σα2=η mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcuwgjbqsgaqcamaabaaaleaacqwgqbgaaeqaaogaeiykakiaeyypa0jae83wdm3aa0baasqaaiab=f7ahbqaaiabikdayaaakiabg2da9iab=d7aonaacaaaleqabagaegomaidaaogaeiikagiaei4waslaf8xsdembakaacqggdbqxcqggpaqkcqgguaglaaa@438a@ in summary, when bcd molecular number is large in the anterior region of the embryo, we expect to see a constant level of normalized variance in our fluorescence data, contributed solely by rescaling noise α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@. thus, α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ can be identified independently from m and other random variables in our data. physical constraints from a high-variance ensemble of embryos we first show results of the above statistical comparison using all flyex cycle embryos. this ensemble contains nuclei, with about nuclei per bin. in figure 2a we see that the normalized variance of the ensemble data η2 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakdaweqbqaaiabdmeajnaabaaaleaacqwgpbqacqwgqbgaaeqaaaqaaiabdmgapbqab0gaesokiufakiabcmcapaaa@376e@ asymptotically approaches the simulation curve η2=σα2 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwcuwfxoqygaqcaiabc2fadjabcmcapiabg2da9iab=n8aznaadaaaleaacqwfxoqyaeaacqaiyagmaaaaaa@3a08@ in the anterior region of the embryo. σα values that are too high will place the black model curve above the data points on the left, and this constrains σα to be less than. mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakdaweqbqaaiabdmeajnaabaaaleaacqwgpbqacqwgqbgaaeqaaaqaaiabdmgapbqab0gaesokiufakiabcmcapaaa@376e@ from the high-variance ensemble of embryos. lines denote simulation results as shown in the key. denotes the full model η2=η2 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcuwgjbqsgaqcamaabaaaleaacqwgqbgaaeqaaogaeiykakiaeyypa0jae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwdaqiaaqaaiab=f7ahjabd6gaujab=j7aibgaaypadagaeiyxa0laeiykakcaaa@4084@, denotes η2 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwcuwfxoqygaqcaiabc2fadjabcmcapaaa@348b@, denotes η2 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwcuwgubgbgaqcaiabc2fadjabcmcapaaa@3456@, and denotes η2 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwcuwfyogygaqcaiabc2fadjabcmcapaaa@348d@. the parameters used in simulation were σα =, σβ =, j =, m =, d = and ω =. is a scatterplot of bcd molecular number n^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@, while the inset shows simulated fluorescence intensity for models and parameters used in panel a. shows the residuals of panel b., and show the same information as panel a, b, and c respectively, but from the low-variance ensemble of embryos, with parameters used in simulation σα =, σβ, =, j =, m =, d = and ω =. note that the axes in panel a and d are scaled differently, and the absolute molecule number is shown in the inset to panel e. with regard to the molecular parameters of diffusion, adopting the value of d = given in the literature together with our observed value of λ = for this ensemble implies that ω = in order to satisfy. these values, for any j, will cause the mean molecular number gradient, 〈n^j〉 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeyykjeuafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgqyixdaa@325e@, to be in steady state after about seconds. a lower bound of j > and an upper bound of m < is imposed by the experimentally observed magnitude of η finally, we estimate the constraint for background noise β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ by satisfying in the posterior end of the embryo. violation of the inequality would cause the black model curve to be above the data on the right hand side of figure 2a, and hence we have an upper bound σβ <. in summary, analysis of our high-variance ensemble of embryos dataset implies that the bcd synthesis rate is higher than. figure 2b shows a scatterplot of the molecular number associated with this synthesis rate. the panel indicates that most nuclei in the anterior fifth of the embryo contain more than molecules of bcd after reaching steady state, and that these molecular numbers can fluctuate over a range of more than molecules in this region. this panel shows that these molecular fluctuations, even at the largest level of normalized variance compatible with data, still do not account for the observed variance in experimental observations. the additional variance comes chiefly from rescaling noise α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@. background noise β^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ only has a significant effect at the posterior end of the embryo, and indeed dominates the normalized variance in that region. towards the posterior, η2 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcqggbbwwcuwfyogygaqcaiabc2fadjabcmcapaaa@348d@ rises faster than η2 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcuwgubgbgaqcamaabaaaleaacqwgqbgaaeqaaogaeiykakcaaa@3369@, and its sharp rise may in certain cases serve as a marker to distinguish regimes dominated by molecular noise from those dominated by background noise. physical constraints from a low-variance ensemble of embryos in the ensemble of embryos discussed above, a portion of the variance observed is likely to stem from embryo to embryo variation in staining conditions and inherent biological parameters. in order to find a better upper limit for observed molecular fluctuations, it is desirable to analyze a set of embryos whose properties are as uniform as possible. in order to generate such a set, we considered the value of λi in the exponential fit f = ai exp to each embryo. we then grouped embryos according to common values of 1/λi, taken to two decimal places. in the largest such group, which contained members, we rescaled the data to a common amplitude by letting i′ij=〈ai〉aiiij. mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xi8qivkypfjydhavhbbf9v8qqaqfr0xc9vqfj0dxdbba91qpepei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaqcfaoafmysakkbauaadawgaaqaaiabdmgapjabdqgaqbqabagaeyypa0zaasaaaeaacqghpms4cqwghbqydawgaaqaaiabdmgapbqabagaeyokjepabagaemyyae2aasbaaeaacqwgpbqaaeqaaaaacqwgjbqsdawgaaqaaiabdmgapjabdqgaqbqabagaeiola4caaa@3fae@ these processed embryos constitute an ensemble ∪ii′ij mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbawaambuaeaacuwgjbqsgaqbamaabaaaleaacqwgpbqacqwgqbgaaeqaaaqaaiabdmgapbqab0gaesokiufaaaa@32e2@ for statistical analysis as described in the previous section. a trade-off of this treatment is the loss of statistical sample size, with only around nuclei in each bin. mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8nsdimbakaaaaa@2d8b@ is too small to be separated from n^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@. we compared this data to simulations performed using the same parameter values reported in the previous section, except that λ had a value of, leading to ω taking on a value of. these values cause the deterministic system to relax to steady state after seconds, as was the case with the high variance ensemble. in the high variance ensemble, the choice of σβ was dictated by the necessity of matching the observed variance along the entire a-p axis. in the case of the low variance ensemble this is not required, but measurements from nonexpressing nuclei indicate that σβ is equal to about it is important to have at least a rough estimate for σβ so that fluctuations from this source are not spuriously assigned to fluctuations in molecular number. finally, these constraints require that the lower bound of synthesis rate j be in order that that η2≃η2 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakcuwgjbqsgaqcamaabaaaleaacqwgqbgaaeqaaogaeiykakiaes4qisjae83tdg2aawbaasqabeaacqaiyagmaagccqggoaakdaweqbqaaiqbdmeajzaafawaasbaasqaaiabdmgapjabdqgaqbqabaaabagaemyaakgabeqdcqwiqisvaogaeiykakcaaa@3feb@. the corresponding upper bound of molecule-to-fluorescence rescaling ratio is m =. the lower limit of j imposed by this ensemble of low variance embryos in turn implies that there must be more than bcd molecules per subvolume in the middle of the embryo after reaching steady state. moreover, it implies that the bcd molecular gradient does not drop to as the fluorescence intensity reaches the presumed background level, but remains at a level of at least molecules per subvolume at j = and molecules per subvolume at the posterior pole. note that the variance of fluorescence measurements is similar between the two embryo ensembles, but that the portion of that variance assigned to fluctuations in molecular number is smaller for the ensemble of embryos. the tighter constraints from the smaller ensemble make the lower limits on molecular number higher, and the lower limit on the variance of molecular number higher, although the higher limit of the normalized variance becomes smaller. the bcd molecular number will thus vary by more than molecules in the middle of the embryo. the 13% rescaling noise α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ is still the main source of the characteristic variation observed in the anterior region of our fluorescence intensity data. noise strength in most applications the most important measure of fluctuation is the normalized variance η. a different quantity, known as the fano factor or noise strength, has been used by some authors as a marker to distinguish different stochastic mechanisms. the fano factor ν is given by ν=var⁡/〈i^j〉 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigae8xvd4maeiikagiafmysakkbakaadawgaawcbagaemoaaogabeaakiabcmcapiabg2da9igbcaha2jabcggahjabckhayjabcicaoiqbdmeajzaajawaasbaasqaaiabdqgaqbqabagccqggpaqkcqggvawlcqghpms4cuwgjbqsgaqcamaabaaaleaacqwgqbgaaeqaaogaeyokjepaaa@42cb@, where ν = in a poisson process. all stochastic simulations of bcd intrinsic noise n^j mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmoba4mbakaadawgaawcbagaemoaaogabeaaaaa@2ed1@ give ν = as shown in figure 3a. by contrast, the full statistical model for either of the two ensembles of embryos examined is obviously non-poisson, not only because ν ≠ but also because in the data, ν is a function of the mean. this happens because at larger values of molecular number, the variance of α^ mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbaaccigaf8xsdembakaaaaa@2d89@ has a dominant role, even when its value is small. even if we model our data without rescaling noise using the mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4waslafmoba4mbakaacqggdbqxaaa@2fc8@ random variable alone, uncertainty in the value of the rescaling constant m itself leads to ambiguity in the observed value of the fano factor. this is a natural consequence of the dimensions of the fano factor. mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagagiodaynaeiyyaemaeiocainaeiikagiafmoba4mbakaadawgaawcbagaemoaaogabeaakiabcmcapaaa@34b9@ divided by molecular mean 〈n^j〉 mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaeyykjeuafmoba4mbakaadawgaawcbagaemoaaogabeaakiabgqyixdaa@325e@. the key indicates that parameters were obtained from the high-variance ensemble of embryos, the low-variance ensemble of embryos, and the extreme condition of high diffusion rate d =, decay rate ω = and synthesis rate j =. as shown in the key, noise strength of simulated fluorescence intensity i^j= mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagafmysakkbakaadawgaawcbagaemoaaogabeaakiabg2da9iabcufabnaahaaabaaccigae8xsdemaemoba4mae8nsdigacaglcmaacqggdbqxaaa@3780@ and rescaled gradient mathtype@mtef@5@5@+=feaafiart1ev1aaatcvaufkttlearuwrp9mdh5mbpbiqv92aaexatlxbi9gbaebbnrfifhhdyfgasaacpc6xni=xh8vivgi8gi=heeeu0xxdbba9frfj0xb9qqpg0dxdb9aspei8k8fii+fsy=rqgqvepae9pg0db9vqaivgfr0xfr=xfr=xc9adbaqaaegacagaaiaabeqaaeqabiwaaagcbagaei4waslafmoba4mbakaacqggdbqxaaa@2fc8@ were obtained using parameters from the high-variance ensemble of embryos and the low-variance ensemble of embryos.  we have compared the nucleus to nucleus variation in expression levels of the exponentially distributed bcd gradient observed in fixed tissue in a steady state with a stochastic model of the diffusion equation. the model is well supported, in the sense that there is a well-supported physical model for the spatial dependence of mean concentrations of bcd on the scale of the embryo. the first major result of our analysis is to note that in many individual embryos the nucleus to nucleus variation in the log of concentration is independent of spatial position. this pattern of variation, which amounts to multiplicative noise in concentration space, is completely incompatible with the stochastic behavior of the diffusion equation. simulations of the diffusion equation over an exhaustively large region of parameter space without exception give rise to solutions in which nucleus to nucleus variation of the bcd gradient is a function of position in the embryo, whether this variation is measured directly in bcd levels or in their logarithms. the data which we compare the model to is in the form of fluorescence levels, not concentrations. although there is now good evidence that the specific batch of serum used to obtain this data has a mean response to bcd which is linear, there is no quantitative information about the variance of this sensitivity. previous work on intrinsic molecular noise in yeast and bacteria utilized gfp in vivo, a situation where fluorescence is detected without molecular amplification. in the data reported here, and in most studies with fixed tissue, the signal from bound primary antibodies is amplified by incubation with secondary antibodies conjugated to a fluorescent dye. it is easily imaginable that this amplification process itself is subject to molecular fluctuations. these fluctuations would then give rise to rescaling noise in the proportionality between fluorescence levels and primary antigen molecular number from nucleus to nucleus. we have shown that such variation can explain the multiplicative noise observed. when data from fixed embryos are pooled, better statistics are obtained. we analyzed data pooled from the entire dataset as well as a smaller pool of data from a set of embryos selected to have nearly identical mean bcd profiles. in order to analyze these data, we constructed an explicit statistical model. the model considers three sources of variance: intrinsic noise, rescaling noise, and background noise. by means of this statistical model it is possible to separate, at least roughly, the contributions of different sources of fluctuation. we have discussed mechanisms for the first two of these; the third, background noise, represents small fluctuations in the quantity of nonspecific molecules from nucleus to nucleus. because mean background has been previously removed from this data, the background noise has a mean of zero. we have confirmed that the background removal method does not affect our results, and the mean background intensity before removal is independent from the bcd molecule-to-fluorescence rescaling ratio. we also found from non-expressing areas of our data that the background noise has about 54% positive correlation to the mean background intensity. the results of this analysis constrained the physical parameters of the stochastic model considerably, with sharper constraints provided by the smaller dataset. the data require that the synthesis rate, j, of bcd from its pool of anteriorly deposited mrna be greater than. we chose the subvolumes of the model to have the same volume as a nucleus, and hence the constraint on j also implies that there are a mean of at least molecules of bcd per nucleus at the midpoint of the embryo, and mean levels of at least molecules of bcd per nucleus at the posterior pole. in terms of concentration, our results show that bcd concentration at the midpoint of the embryo is greater than nm. recently, a direct in vivo measurement of bcd concentration was performed and yielded a mid embryo bcd concentration of nm, fully compatible with our results. although we are able to extract a clear signature of intrinsic molecular noise from the data by means of the statistical model, we also showed that at least one quantity diagnostic for different stochastic mechanisms, the fano factor, cannot be read out from the data. although the fano factor v = in the simulations, the full statistical model gives rise to a fano factor which is a function of the mean, and even if all noise is restricted to be intrinsic, the observed fano factor depends on the scale for conversion from fluorescence to molecular number. we believe that our results in general demonstrate that fixed material processed with secondary fluorophores is not well suited to studies of molecular fluctuations. this arises from three issues, which may be separable. fixation obviously prevents repeated observations on the same cell. while that is clearly a limitation, it need not affect an investigation of a molecule whose mean values are in steady state. the other two issues concern amplification. gfp is intrinsically fluorescent, but antibodies must bind to antigen, a process that is in itself subject to intrinsic molecular fluctuations. in the present study, the level of molecular fluctuation is doubtless increased by the need to bind secondary antibodies conjugated to fluorophore to the already bound primary antibody. it is thus possible that better data from fixed tissue could be obtained by conjugating dye directly to the primary antibody. this may prove to be a useful measure in situations where constructing a gfp fusion that can functionally substitute for the native gene is difficult or impossible. more generally, we suggest that it is important to know how precision and robustness of developmental control is achieved at the molecular number level throughout development. indeed, there is little doubt that stochastic processes are important later in development. adult drosophila normally have four scutellar bristles. in mutants of the scute gene, the number of bristles varies between one and three, strongly suggesting a stochastic process. on a more theoretical level, it has been suggested that fluctuations can augment the operating capabilities of biological regulatory networks. our results indicate that in vivo monitoring of gene expression will be required to obtain high quality data on stochastic gene expression phenomena in eucaryotes. the central technical problem that must be solved to conduct such studies is the complete replacement of the endogenous gene with a fluorescently tagged functional version. 
an estimate based on a low variance ensemble of embryos suggests that the steady-state bicoid molecular number in a nucleus should be larger than in the middle of the embryo, and hence the gradient should extend to the posterior end of the embryo, beyond the previously assumed background limit. we show that the intrinsic noise of an autonomous reaction-diffusion gradient is poisson distributed. we demonstrate how experimental noise can be identified in the logarithm domain from single embryo analysis, and then separated from intrinsic noise in the normalized variance domain of an ensemble statistical analysis. intrinsic noise is clearly detectable as well, although the staining process introduces certain limits in resolution. we compared quantitative observations of bicoid levels in immunostained drosophila embryos with a spatially extended master equation model which represents diffusion, decay, and anterior synthesis. stochastic simulations can be applied to further the understanding of the dynamics of bicoid morphogen gradient formation at the molecular number level, and determine the source of the nucleus-to-nucleus expression variation observed in the bicoid gradient. we exhibit the predicted molecular number gradient together with measurement effects, and make a comparison between conditions of higher and lower variance respectively. we show how measurement sensitivity affects our observations, and how small amounts of rescaling noise can perturb the noise strength observed. because different developmental fates can be elicited by different concentrations of bicoid, it is important to probe the limits of this specification by analyzing intrinsic fluctuations of the bicoid gradient arising from small molecular number.
 the concentration gradient of bicoid protein which determines the developmental pathways in early drosophila embryo is the best characterized morphogen gradient at the molecular level. because different developmental fates can be elicited by different concentrations of bicoid, it is important to probe the limits of this specification by analyzing intrinsic fluctuations of the bicoid gradient arising from small molecular number. stochastic simulations can be applied to further the understanding of the dynamics of bicoid morphogen gradient formation at the molecular number level, and determine the source of the nucleus-to-nucleus expression variation observed in the bicoid gradient.  we compared quantitative observations of bicoid levels in immunostained drosophila embryos with a spatially extended master equation model which represents diffusion, decay, and anterior synthesis. we show that the intrinsic noise of an autonomous reaction-diffusion gradient is poisson distributed. we demonstrate how experimental noise can be identified in the logarithm domain from single embryo analysis, and then separated from intrinsic noise in the normalized variance domain of an ensemble statistical analysis. we show how measurement sensitivity affects our observations, and how small amounts of rescaling noise can perturb the noise strength observed. we demonstrate that the biological noise level in data can serve as a physical constraint for restricting the model's parameter space, and for predicting the bicoid molecular number and variation range. an estimate based on a low variance ensemble of embryos suggests that the steady-state bicoid molecular number in a nucleus should be larger than in the middle of the embryo, and hence the gradient should extend to the posterior end of the embryo, beyond the previously assumed background limit. we exhibit the predicted molecular number gradient together with measurement effects, and make a comparison between conditions of higher and lower variance respectively.  quantitative comparison of master equation simulations with immunostained data enabled us to determine narrow ranges for key biophysical parameters, which for this system can be independently validated. intrinsic noise is clearly detectable as well, although the staining process introduces certain limits in resolution. 