Novel approach to quantitative spatial gene expression uncovers genetic stochasticity in the developing Drosophila eye.

Robustness in development allows for the accumulation of genetically based variation in expression. However, this variation is usually examined in response to large perturbations, and examination of this variation has been limited to being spatial, or quantitative, but because of technical restrictions not both. Here we bridge these gaps by investigating replicated quantitative spatial gene expression using rigorous statistical models, in different genotypes, sexes, and species (Drosophila melanogaster and D. simulans). Using this type of quantitative approach with molecular developmental data allows for comparison among conditions, such as different genetic backgrounds. We apply this approach to the morphogenetic furrow, a wave of differentiation that patterns the developing eye disc. Within the morphogenetic furrow, we focus on four genes, hairy, atonal, hedgehog, and Delta. Hybridization chain reaction quantitatively measures spatial gene expression, co-staining for all four genes simultaneously. We find considerable variation in the spatial expression pattern of these genes in the eye between species, genotypes, and sexes. We also find that there has been evolution of the regulatory relationship between these genes, and that their spatial interrelationships have evolved between species. This variation has no phenotypic effect, and could be buffered by network thresholds or compensation from other genes. Both of these mechanisms could potentially be contributing to long term developmental systems drift.

Threshold response to stochasticity in morphogenesis.

During development of biological organisms, multiple complex structures are formed. In many instances, these structures need to exhibit a high degree of order to be functional, although many of their constituents are intrinsically stochastic. Hence, it has been suggested that biological robustness ultimately must rely on complex gene regulatory networks and clean-up mechanisms. Here we explore developmental processes that have evolved inherent robustness against stochasticity. In the context of the Drosophila eye disc, multiple optical units, ommatidia, develop into crystal-like patterns. During the larva-to-pupa stage of metamorphosis, the centers of the ommatidia are specified initially through the diffusion of morphogens, followed by the specification of R8 cells. Establishing the R8 cell is crucial in setting up the geometric, and functional, relationships of cells within an ommatidium and among neighboring ommatidia. Here we study an PDE mathematical model of these spatio-temporal processes in the presence of parametric stochasticity, defining and applying measures that quantify order within the resulting spatial patterns. We observe a universal sigmoidal response to increasing transcriptional noise. Ordered patterns persist up to a threshold noise level in the model parameters. In accordance with prior qualitative observations, as the noise is further increased past a threshold point of no return, these ordered patterns rapidly become disordered. Such robustness in development allows for the accumulation of genetic variation without any observable changes in phenotype. We argue that the observed sigmoidal dependence introduces robustness allowing for sizable amounts of genetic variation and transcriptional noise to be tolerated in natural populations without resulting in phenotype variation.

Genetic Determinants of RNA Editing Levels of ADAR Targets in Drosophila melanogaster.

RNA editing usually affects only a fraction of expressed transcripts and there is a vast amount of variation in editing levels of ADAR (adenosine deaminase, RNA-specific) targets. Here we explore natural genetic variation affecting editing levels of particular sites in 81 natural strains of Drosophila melanogaster. The analysis of associations between editing levels and single-nucleotide polymorphisms allows us to map putative cis-regulatory regions affecting editing of 16 A-to-I editing sites (cis-RNA editing quantitative trait loci or cis-edQTLs, P < 10(-8)). The observed changes in editing levels are validated by independent molecular technique. All identified regulatory variants are located in close proximity of modulated editing sites. Moreover, colocalized editing sites are often regulated by same loci. Similar to expression and splicing QTL studies, the characterization of edQTLs will greatly expand our understanding of cis-regulatory evolution of gene expression.