Abundant genetic variation in transcript level during early Drosophila development

Genetic, Environmental, and Stochastic Components of Lifespan Variability: The Drosophila Paradigm

Lifespan is a complex quantitative trait involving genetic and non-genetic factors as well as the peculiarities of ontogenesis. As with all quantitative traits, lifespan shows considerable variation within populations and between individuals. Drosophila, a favourite object of geneticists, has greatly advanced our understanding of how different forms of variability affect lifespan. This review considers the role of heritable genetic variability, phenotypic plasticity and stochastic variability in controlling lifespan in Drosophila melanogaster. We discuss the major historical milestones in the development of the genetic approach to study lifespan, the breeding of long-lived lines, advances in lifespan QTL mapping, the environmental factors that have the greatest influence on lifespan in laboratory maintained flies, and the mechanisms, by which individual development affects longevity. The interplay between approaches to study ageing and lifespan limitation will also be discussed. Particular attention will be paid to the interaction of different types of variability in the control of lifespan.

Evolution and Developmental System Drift in the Endoderm Gene Regulatory Network of Caenorhabditis and Other Nematodes

Developmental gene regulatory networks (GRNs) underpin metazoan embryogenesis and have undergone substantial modification to generate the tremendous variety of animal forms present on Earth today. The nematode Caenorhabditis elegans has been a central model for advancing many important discoveries in fundamental mechanistic biology and, more recently, has provided a strong base from which to explore the evolutionary diversification of GRN architecture and developmental processes in other species. In this short review, we will focus on evolutionary diversification of the GRN for the most ancient of the embryonic germ layers, the endoderm. Early embryogenesis diverges considerably across the phylum Nematoda. Notably, while some species deploy regulative development, more derived species, such as C. elegans, exhibit largely mosaic modes of embryogenesis. Despite the relatively similar morphology of the nematode gut across species, widespread variation has been observed in the signaling inputs that initiate the endoderm GRN, an exemplar of developmental system drift (DSD). We will explore how genetic variation in the endoderm GRN helps to drive DSD at both inter- and intraspecies levels, thereby resulting in a robust developmental system. Comparative studies using divergent nematodes promise to unveil the genetic mechanisms controlling developmental plasticity and provide a paradigm for the principles governing evolutionary modification of an embryonic GRN.

Inference of Transcription Factor Regulation Patterns Using Gene Expression Covariation in Natural Populations of Drosophila melanogaster

Gene regulatory networks control the complex programs that drive development. Deciphering the connections between transcription factors (TFs) and target genes is challenging, in part because TFs bind to thousands of places in the genome but control expression through a subset of these binding events. We hypothesize that we can combine natural variation of expression levels and predictions of TF binding sites to identify TF targets. We gather RNA-seq data from 71 genetically distinct F1 Drosophila melanogaster embryos and calculate the correlations between TF and potential target genes' expression levels, which we call "regulatory strength." To separate direct and indirect TF targets, we hypothesize that direct TF targets will have a preponderance of binding sites in their upstream regions. Using 14 TFs active during embryogenesis, we find that 12 TFs showed a significant correlation between their binding strength and regulatory strength on downstream targets, and 10 TFs showed a significant correlation between the number of binding sites and the regulatory effect on target genes. The general roles, e.g. bicoid's role as an activator, and the particular interactions we observed between our TFs, e.g. twist's role as a repressor of sloppy paired and odd paired, generally coincide with the literature.

Parallel embryonic transcriptional programs evolve under distinct constraints and may enable morphological conservation amidst adaptation

Embryonic development evolves by balancing stringent morphological constraints with genetic and environmental variation. The design principle that allows developmental transcriptional programs to conserve embryonic morphology while adapting to environmental changes is still not fully understood. To address this fundamental challenge, we compare developmental transcriptomes of two sea urchin species, Paracentrotus lividus and Strongylocentrotus purpuratus, that shared a common ancestor about 40 million years ago and are geographically distant yet show similar morphology. We find that both developmental and housekeeping genes show highly dynamic and strongly conserved temporal expression patterns. The expression of other gene sets, including homeostasis and response genes, show divergent expression which could result from either evolutionary drift or adaptation to local environmental conditions. The interspecies correlations of developmental gene expressions are highest between morphologically similar developmental time points whereas the interspecies correlations of housekeeping gene expression are high between all the late zygotic time points. Relatedly, the position of the phylotypic stage varies between these two groups of genes: developmental gene expression shows highest conservation at mid-developmental stage, in agreement with the hourglass model while the conservation of housekeeping genes keeps increasing with developmental time. When all genes are combined, the relationship between conservation of gene expression and morphological similarity is partially masked by housekeeping genes and genes with diverged expression. Our study illustrates various transcriptional programs that coexist in the developing embryo and evolve under different constraints. Apparently, morphological constraints underlie the conservation of developmental gene expression while embryonic fitness requires the conservation of housekeeping gene expression and the species-specific adjustments of homeostasis gene expression. The distinct evolutionary forces acting on these transcriptional programs enable the conservation of similar body plans while allowing adaption.

Genomic Characterization of the Evolutionary Potential of the Sea Urchin Strongylocentrotus droebachiensis Facing Ocean Acidification

Ocean acidification (OA) is increasing due to anthropogenic CO₂ emissions and poses a threat to marine species and communities worldwide. To better project the effects of acidification on organisms’ health and persistence, an understanding is needed of the 1) mechanisms underlying developmental and physiological tolerance and 2) potential populations have for rapid evolutionary adaptation. This is especially challenging in nonmodel species where targeted assays of metabolism and stress physiology may not be available or economical for large-scale assessments of genetic constraints. We used mRNA sequencing and a quantitative genetics breeding design to study mechanisms underlying genetic variability and tolerance to decreased seawater pH (-0.4 pH units) in larvae of the sea urchin Strongylocentrotus droebachiensis. We used a gene ontology-based approach to integrate expression profiles into indirect measures of cellular and biochemical traits underlying variation in larval performance (i.e., growth rates). Molecular responses to OA were complex, involving changes to several functions such as growth rates, cell division, metabolism, and immune activities. Surprisingly, the magnitude of pH effects on molecular traits tended to be small relative to variation attributable to segregating functional genetic variation in this species. We discuss how the application of transcriptomics and quantitative genetics approaches across diverse species can enrich our understanding of the biological impacts of climate change.

The impact of gene expression variation on the robustness and evolvability of a developmental gene regulatory network

Regulatory interactions buffer development against genetic and environmental perturbations, but adaptation requires phenotypes to change. We investigated the relationship between robustness and evolvability within the gene regulatory network underlying development of the larval skeleton in the sea urchin Strongylocentrotus purpuratus. We find extensive variation in gene expression in this network throughout development in a natural population, some of which has a heritable genetic basis. Switch-like regulatory interactions predominate during early development, buffer expression variation, and may promote the accumulation of cryptic genetic variation affecting early stages. Regulatory interactions during later development are typically more sensitive (linear), allowing variation in expression to affect downstream target genes. Variation in skeletal morphology is associated primarily with expression variation of a few, primarily structural, genes at terminal positions within the network. These results indicate that the position and properties of gene interactions within a network can have important evolutionary consequences independent of their immediate regulatory role.

Dissecting sources of quantitative gene expression pattern divergence between Drosophila species

Gene expression patterns can diverge between species due to changes in a gene's regulatory DNA or changes in the proteins, e.g., transcription factors (TFs), that regulate the gene. We developed a modeling framework to uncover the sources of expression differences in blastoderm embryos of three Drosophila species, focusing on the regulatory circuit controlling expression of the hunchback (hb) posterior stripe. Using this framework and cellular-resolution expression measurements of hb and its regulating TFs, we found that changes in the expression patterns of hb's TFs account for much of the expression divergence. We confirmed our predictions using transgenic D. melanogaster lines, which demonstrate that this set of orthologous cis-regulatory elements (CREs) direct similar, but not identical, expression patterns. We related expression pattern differences to sequence changes in the CRE using a calculation of the CRE's TF binding site content. By applying this calculation in both the transgenic and endogenous contexts, we found that changes in binding site content affect sensitivity to regulating TFs and that compensatory evolution may occur in circuit components other than the CRE.

A model-based method for gene dependency measurement

Many computational methods have been widely used to identify transcription regulatory interactions based on gene expression profiles. The selection of dependency measure is very important for successful regulatory network inference. In this paper, we develop a new method-DBoMM (Difference in BIC of Mixture Models)-for estimating dependency of gene by fitting the gene expression profiles into mixture Gaussian models. We show that DBoMM out-performs 4 other existing methods, including Kendall's tau correlation (TAU), Pearson Correlation (COR), Euclidean distance (EUC) and Mutual information (MI) using Escherichia coli, Saccharomyces cerevisiae, Drosophila melanogaster, Arabidopsis thaliana data and synthetic data. DBoMM can also identify condition-dependent regulatory interactions and is robust to noisy data. Of the 741 Escherichia coli regulatory interactions inferred by DBoMM at a 60% true positive rate, 65 are previously known interactions and 676 are novel predictions. To validate the new prediction, the promoter sequences of target genes regulated by the same transcription factors were analyzed and significant motifs were identified.

The effects of weak genetic perturbations on the transcriptome of the wing imaginal disc and its association with wing shape in Drosophila melanogaster

A major objective of genomics is to elucidate the mapping between genotypic and phenotypic space as a step toward understanding how small changes in gene function can lead to elaborate phenotypic changes. One approach that has been utilized is to examine overall patterns of covariation between phenotypic variables of interest, such as morphology, physiology, and behavior, and underlying aspects of gene activity, in particular transcript abundance on a genome-wide scale. Numerous studies have demonstrated that such patterns of covariation occur, although these are often between samples with large numbers of unknown genetic differences (different strains or even species) or perturbations of large effect (sexual dimorphism or strong loss-of-function mutations) that may represent physiological changes outside of the normal experiences of the organism. We used weak mutational perturbations in genes affecting wing development in Drosophila melanogaster that influence wing shape relative to a co-isogenic wild type. We profiled transcription of 1150 genes expressed during wing development in 27 heterozygous mutants, as well as their co-isogenic wild type and one additional wild-type strain. Despite finding clear evidence of expression differences between mutants and wild type, transcriptional profiles did not covary strongly with shape, suggesting that information from transcriptional profiling may not generally be predictive of final phenotype. We discuss these results in the light of possible attractor states of gene expression and how this would affect interpretation of covariation between transcriptional profiles and other phenotypes.

The strength of transcription-factor binding modulates co-variation in transcriptional networks

An appreciable fraction of the transcriptome differs in level of expression among individuals. Transcription factor (TF) expression and DNA binding causes cell-specific activation and repression of downstream targets, and TF expression levels vary across individuals. However, it is not clear how the strength of DNA binding for individual TFs translates into regulatory control, or whether a different set of binding motifs is used for strongly regulated modules. Here we integrate two publicly available data sets in Drosophila melanogaster, as well as conduct novel analyses, to address these questions.

Genotype-phenotype mapping developmental biology confronts the toolkit paradox

The quest to understand the relationship between an organism's DNA sequence and three-dimensional form is an interdisciplinary task, integrating diverse fields of the life sciences. The relevance of the metaphor of a genotype-phenotype map is explored from a developmental perspective, in light of the recent concept of a "molecular toolkit" of protein-coding genes, and the widespread view that analyzing the logic and mechanics of gene regulation at multiple levels is key to explaining how morphology is genetically encoded. We discuss the challenges of decoding genomes despite variable genetic backgrounds, the dynamically changing physical and molecular contexts of the internal environment during development, and the impact of external forces on morphogenesis.