Gene expression is controlled by the action of transcription factors that bind to DNA and influence the rate at which a gene is transcribed. The quantitative mapping between the regulator concentrations and the output of the gene is known as the cisregulatory input function (CRIF). Here, we show how the CRIF shapes the form of the joint probability distribution of molecular copy numbers of the regulators and the product of a gene. Namely, we derive a class of fluctuation-based relations that relate the moments of the distribution to the derivatives of the CRIF. These relations are useful because they enable statistics of naturally arising cell-to-cell variations in molecular copy numbers to substitute for traditional manipulations for probing regulatory mechanisms. We demonstrate that these relations can distinguish super-and subadditive gene regulatory scenarios (molecular analogs of AND and OR logic operations) in simulations that faithfully represent bacterial gene expression. Applications and extensions to other regulatory scenarios are discussed.gene expression | mathematical modeling | noise analysis inference | flow cytometry T ranscription factors regulate the expression of genes by binding to specific sites on the DNA that are typically spatially close to, sometimes even in, sequences that code for proteins. A group of such binding sites is collectively known as a cis-regulatory region. When a gene processes the effects of multiple transcription factors, one can view it as performing a computation. The inputs are the occupancies of the binding sites in the cis-regulatory region and the output is the rate of transcription. The simplest method of describing the relationship between the transcription factors and the output is by using boolean logic (1, 2). For example, 2 activators can regulate a gene with AND logic in which both are required or with OR logic in which either transcription factor is sufficient for transcription. Knowing which mode of combinatorial regulation a gene employs can be important for determining its function in regulatory networks. For example, the coherent feed forward loop, one of the most common motifs in gene regulatory networks (3), filters noise in upstream signals differently depending on how one of the participating genes integrates its inputs (4, 5).In general, the output of a gene will not be binary, and will depend on the concentrations of the transcription factors in a complex manner. The notion of logic operations can be generalized by introducing a continuous function that encodes the dependence of the rate of transcription on the concentrations of the inputs. Such "cis-regulatory input functions (CRIFs)" (5) have been evaluated experimentally for the well-studied lac operon. The CRIF of the wild type is complex (6, 7), but those of certain mutant operons represent molecular analogs of binary logic gates (8). In higher organisms, many genes are regulated by a large number of transcription factors, often with multiple binding sites for each factor (9-11). It is increa...