Epidermal growth factor receptor and transforming growth factor-beta signaling contributes to variation for wing shape in Drosophila melanogaster

Odd-Paired is Involved in Morphological Divergence of Snail-Feeding Beetles

Body shape and size diversity and their evolutionary rates correlate with species richness at the macroevolutionary scale. However, the molecular genetic mechanisms underlying the morphological diversification across related species are poorly understood. In beetles, which account for one-fourth of the known species, adaptation to different trophic niches through morphological diversification appears to have contributed to species radiation. Here, we explored the key genes for the morphological divergence of the slender to stout body shape related to divergent feeding methods on large to small snails within the genus Carabus. We show that the zinc-finger transcription factor encoded by odd-paired (opa) controls morphological variation in the snail-feeding ground beetle Carabus blaptoides. Specifically, opa was identified as the gene underlying the slender to stout morphological difference between subspecies through genetic mapping and functional analysis via gene knockdown. Further analyses revealed that changes in opa cis-regulatory sequences likely contributed to the differences in body shape and size between C. blaptoides subspecies. Among opa cis-regulatory sequences, single nucleotide polymorphisms on the transcription factor binding sites may be associated with the morphological differences between C. blaptoides subspecies. opa was highly conserved in a wide range of taxa, especially in beetles. Therefore, opa may play an important role in adaptive morphological divergence in beetles.

Recombinant inbred line panels inform the genetic architecture and interactions of adaptive traits in Drosophila melanogaster

The distribution of allelic effects on traits, along with their gene-by-gene and gene-by-environment interactions, contributes to the phenotypes available for selection and the trajectories of adaptive variants. Nonetheless, uncertainty persists regarding the effect sizes underlying adaptations and the importance of genetic interactions. Herein, we aimed to investigate the genetic architecture and the epistatic and environmental interactions involving loci that contribute to multiple adaptive traits using two new panels of Drosophila melanogaster recombinant inbred lines (RILs). To better fit our data, we re-implemented functions from R/qtl (Broman et al. 2003) using additive genetic models. We found 14 quantitative trait loci (QTL) underlying melanism, wing size, song pattern, and ethanol resistance. By combining our mapping results with population genetic statistics, we identified potential new genes related to these traits. None of the detected QTLs showed clear evidence of epistasis, and our power analysis indicated that we should have seen at least one significant interaction if sign epistasis or strong positive epistasis played a pervasive role in trait evolution. In contrast, we did find roles for gene-by-environment interactions involving pigmentation traits. Overall, our data suggest that the genetic architecture of adaptive traits often involves alleles of detectable effect, that strong epistasis does not always play a role in adaptation, and that environmental interactions can modulate the effect size of adaptive alleles.

The Extraordinary Diversity of Merodon avidus Complex (Diptera: Syrphidae)-Adding New Areas, New Species and a New Molecular Marker

In this paper, the Merodon avidus (Diptera, Syrphidae) species complex was revised, whereupon we discovered and described four new species for science: Merodon atroavidus Vujić, Radenković et Likov sp. nov., M. magnus Vujić, Kočiš Tubić et Ačanski sp. nov., M. nigroscutum Vujić, Radenković et Likov sp. nov. and M. pseudomoenium Vujić, Kočiš Tubić et Ačanski sp. nov. An integrative taxonomy approach was used to delimit species boundaries. Two molecular markers (the mitochondrial COI gene and nuclear 28S rRNA gene-newly analysed marker for the complex) and geometric morphometry of the wing shape, together with morphological data and distribution, successfully separated all species from the complex. The morphological variability of the analysed species is described and discussed and an illustrated diagnostic key for typical morpho-forms of species from the M. avidus complex is presented. A distribution map of all investigated species from the complex is provided. The level of endemicity of the M. avidus complex was discussed.

Complexities of recapitulating polygenic effects in natural populations: replication of genetic effects on wing shape in artificially selected and wild-caught populations of Drosophila melanogaster

Identifying the genetic architecture of complex traits is important to many geneticists, including those interested in human disease, plant and animal breeding, and evolutionary genetics. Advances in sequencing technology and statistical methods for genome-wide association studies have allowed for the identification of more variants with smaller effect sizes, however, many of these identified polymorphisms fail to be replicated in subsequent studies. In addition to sampling variation, this failure to replicate reflects the complexities introduced by factors including environmental variation, genetic background, and differences in allele frequencies among populations. Using Drosophila melanogaster wing shape, we ask if we can replicate allelic effects of polymorphisms first identified in a genome-wide association studies in three genes: dachsous, extra-macrochaete, and neuralized, using artificial selection in the lab, and bulk segregant mapping in natural populations. We demonstrate that multivariate wing shape changes associated with these genes are aligned with major axes of phenotypic and genetic variation in natural populations. Following seven generations of artificial selection along the dachsous shape change vector, we observe genetic differentiation of variants in dachsous and genomic regions containing other genes in the hippo signaling pathway. This suggests a shared direction of effects within a developmental network. We also performed artificial selection with the extra-macrochaete shape change vector, which is not a part of the hippo signaling network, but showed a largely shared direction of effects. The response to selection along the emc vector was similar to that of dachsous, suggesting that the available genetic diversity of a population, summarized by the genetic (co)variance matrix (G), influenced alleles captured by selection. Despite the success with artificial selection, bulk segregant analysis using natural populations did not detect these same variants, likely due to the contribution of environmental variation and low minor allele frequencies, coupled with small effect sizes of the contributing variants.

The evolution of larger size in high-altitude Drosophila melanogaster has a variable genetic architecture

Important uncertainties persist regarding the genetic architecture of adaptive trait evolution in natural populations, including the number of genetic variants involved, whether they are drawn from standing genetic variation, and whether directional selection drives them to complete fixation. Here, we take advantage of a unique natural population of Drosophila melanogaster from the Ethiopian highlands, which has evolved larger body size than any other known population of this species. We apply a bulk segregant quantitative trait locus mapping approach to 4 unique crosses between highland Ethiopian and lowland Zambian populations for both thorax length and wing length. Results indicated a persistently variable genetic basis for these evolved traits (with largely distinct sets of quantitative trait loci for each cross), and at least a moderately polygenic architecture with relatively strong effects present. We complemented these mapping experiments with population genetic analyses of quantitative trait locus regions and gene ontology enrichment analysis, generating strong hypotheses for specific genes and functional processes that may have contributed to these adaptive trait changes. Finally, we find that the genetic architectures indicated by our quantitative trait locus mapping results for size traits mirror those from similar experiments on other recently evolved traits in this species. Collectively, these studies suggest a recurring pattern of polygenic adaptation in this species, in which causative variants do not approach fixation and moderately strong effect loci are present.

The gut-microbiota-brain axis in autism: what Drosophila models can offer?

The idea that alterations in gut-microbiome-brain axis (GUMBA)-mediated communication play a crucial role in human brain disorders like autism remains a topic of intensive research in various labs. Gastrointestinal issues are a common comorbidity in patients with autism spectrum disorder (ASD). Although gut microbiome and microbial metabolites have been implicated in the etiology of ASD, the underlying molecular mechanism remains largely unknown. In this review, we have summarized recent findings in human and animal models highlighting the role of the gut-brain axis in ASD. We have discussed genetic and neurobehavioral characteristics of Drosophila as an animal model to study the role of GUMBA in ASD. The utility of Drosophila fruit flies as an amenable genetic tool, combined with axenic and gnotobiotic approaches, and availability of transgenic flies may reveal mechanistic insight into gut-microbiota-brain interactions and the impact of its alteration on behaviors relevant to neurological disorders like ASD.

Selective Disruption of Synaptic BMP Signaling by a Smad Mutation Adjacent to the Highly Conserved H2 Helix

Bone morphogenetic proteins (BMPs) shape normal development and function via canonical and noncanonical signaling pathways. BMPs initiate canonical signaling by binding to transmembrane receptors that phosphorylate Smad proteins and induce their translocation into the nucleus and regulation of target genes. Phosphorylated Smads also accumulate at cellular junctions, but this noncanonical, local BMP signaling modality remains less defined. We have recently reported that phosphorylated Smad (pMad in Drosophila) accumulates at synaptic junctions in protein complexes with genetically distinct composition and regulation. Here, we examined a wide collection of DrosophilaMad alleles and searched for molecular features relevant to pMad accumulation at synaptic junctions. We found that strong Mad alleles generally disrupt both synaptic and nuclear pMad, whereas moderate Mad alleles have a wider range of phenotypes and can selectively impact different BMP signaling pathways. Interestingly, regulatory Mad mutations reveal that synaptic pMad appears to be more sensitive to a net reduction in Mad levels than nuclear pMad. Importantly, a previously uncharacterized allele, Mad8 , showed markedly reduced synaptic pMad but only moderately diminished nuclear pMad. The postsynaptic composition and electrophysiological properties of Mad8 neuromuscular junctions (NMJs) were also altered. Using biochemical approaches, we examined how a single point mutation in Mad8 could influence the Mad-receptor interface and identified a key motif, the H2 helix. Our study highlights the biological relevance of Smad-dependent, synaptic BMP signaling and uncovers a highly conserved structural feature of Smads, critical for normal development and function.

Drosophila models of pathogenic copy-number variant genes show global and non-neuronal defects during development

While rare pathogenic copy-number variants (CNVs) are associated with both neuronal and non-neuronal phenotypes, functional studies evaluating these regions have focused on the molecular basis of neuronal defects. We report a systematic functional analysis of non-neuronal defects for homologs of 59 genes within ten pathogenic CNVs and 20 neurodevelopmental genes in Drosophila melanogaster. Using wing-specific knockdown of 136 RNA interference lines, we identified qualitative and quantitative phenotypes in 72/79 homologs, including 21 lines with severe wing defects and six lines with lethality. In fact, we found that 10/31 homologs of CNV genes also showed complete or partial lethality at larval or pupal stages with ubiquitous knockdown. Comparisons between eye and wing-specific knockdown of 37/45 homologs showed both neuronal and non-neuronal defects, but with no correlation in the severity of defects. We further observed disruptions in cell proliferation and apoptosis in larval wing discs for 23/27 homologs, and altered Wnt, Hedgehog and Notch signaling for 9/14 homologs, including AATF/Aatf, PPP4C/Pp4-19C, and KIF11/Klp61F. These findings were further supported by tissue-specific differences in expression patterns of human CNV genes, as well as connectivity of CNV genes to signaling pathway genes in brain, heart and kidney-specific networks. Our findings suggest that multiple genes within each CNV differentially affect both global and tissue-specific developmental processes within conserved pathways, and that their roles are not restricted to neuronal functions.

Rare copy-number variants (CNVs), or large deletions and duplications in the genome, are associated with both neuronal and non-neuronal clinical features. Previous functional studies for these disorders have primarily focused on understanding the cellular mechanisms for neurological and behavioral phenotypes. To understand how genes within these CNVs contribute to developmental defects in non-neuronal tissues, we assessed 79 homologs of CNV and known neurodevelopmental genes in Drosophila models. We found that most homologs showed developmental defects when knocked down in the adult fly wing, ranging from mild size changes to severe wrinkled wings or lethality. Although a majority of tested homologs showed defects when knocked down specifically in wings or eyes, we found no correlation in the severity of the observed defects in these two tissues. A subset of the homologs showed disruptions in cellular processes in the developing fly wing, including alterations in cell proliferation, apoptosis, and cellular signaling pathways. Furthermore, human CNV genes also showed differences in gene expression patterns and interactions with signaling pathway genes across multiple human tissues. Our findings suggest that genes within CNV disorders affect global developmental processes in both neuronal and non-neuronal tissues.

Integrating GWAS and Transcriptomics to Identify the Molecular Underpinnings of Thermal Stress Responses in Drosophila melanogaster

Thermal tolerance of an organism depends on both the ability to dynamically adjust to a thermal stress and preparatory developmental processes that enhance thermal resistance. However, the extent to which standing genetic variation in thermal tolerance alleles influence dynamic stress responses vs. preparatory processes is unknown. Here, using the model species Drosophila melanogaster, we used a combination of Genome Wide Association mapping (GWAS) and transcriptomic profiling to characterize whether genes associated with thermal tolerance are primarily involved in dynamic stress responses or preparatory processes that influence physiological condition at the time of thermal stress. To test our hypotheses, we measured the critical thermal minimum (CTmin) and critical thermal maximum (CTmax) of 100 lines of the Drosophila Genetic Reference Panel (DGRP) and used GWAS to identify loci that explain variation in thermal limits. We observed greater variation in lower thermal limits, with CTmin ranging from 1.81 to 8.60°C, while CTmax ranged from 38.74 to 40.64°C. We identified 151 and 99 distinct genes associated with CTmin and CTmax, respectively, and there was strong support that these genes are involved in both dynamic responses to thermal stress and preparatory processes that increase thermal resistance. Many of the genes identified by GWAS were involved in the direct transcriptional response to thermal stress (72/151 for cold; 59/99 for heat), and overall GWAS candidates were more likely to be differentially expressed than other genes. Further, several GWAS candidates were regulatory genes that may participate in the regulation of stress responses, and gene ontologies related to development and morphogenesis were enriched, suggesting many of these genes influence thermal tolerance through effects on development and physiological status. Overall, our results suggest that thermal tolerance alleles can influence both dynamic plastic responses to thermal stress and preparatory processes that improve thermal resistance. These results also have utility for directly comparing GWAS and transcriptomic approaches for identifying candidate genes associated with thermal tolerance.

G. Teixeira J, Ferreira H, J. Neves C, Francoy TM, Pinto MA. Wing Geometric Morphometrics of Workers and Drones and Single Nucleotide Polymorphisms Provide Similar Genetic Structure in the Iberian Honey Bee (Apis mellifera iberiensis)

Wing geometric morphometrics has been applied to honey bees (Apis mellifera) in identification of evolutionary lineages or subspecies and, to a lesser extent, in assessing genetic structure within subspecies. Due to bias in the production of sterile females (workers) in a colony, most studies have used workers leaving the males (drones) as a neglected group. However, considering their importance as reproductive individuals, the use of drones should be incorporated in these analyses in order to better understand diversity patterns and underlying evolutionary processes. Here, we assessed the usefulness of drone wings, as well as the power of wing geometric morphometrics, in capturing the signature of complex evolutionary processes by examining wing shape data, integrated with geographical information, from 711 colonies sampled across the entire distributional range of Apis mellifera iberiensis in Iberia. We compared the genetic patterns reconstructed from spatially-explicit shape variation extracted from wings of both sexes with that previously reported using 383 genome-wide SNPs (single nucleotide polymorphisms). Our results indicate that the spatial structure retrieved from wings of drones and workers was similar (r = 0.93) and congruent with that inferred from SNPs (r = 0.90 for drones; r = 0.87 for workers), corroborating the clinal pattern that has been described for A. m. iberiensis using other genetic markers. In addition to showing that drone wings carry valuable genetic information, this study highlights the capability of wing geometric morphometrics in capturing complex genetic patterns, offering a reliable and low-cost alternative for preliminary estimation of population structure.

A role for triglyceride lipase brummer in the regulation of sex differences in Drosophila fat storage and breakdown

Triglycerides are the major form of stored fat in all animals. One important determinant of whole-body fat storage is whether an animal is male or female. Here, we use Drosophila, an established model for studies on triglyceride metabolism, to gain insight into the genes and physiological mechanisms that contribute to sex differences in fat storage. Our analysis of triglyceride storage and breakdown in both sexes identified a role for triglyceride lipase brummer (bmm) in the regulation of sex differences in triglyceride homeostasis. Normally, male flies have higher levels of bmm mRNA both under normal culture conditions and in response to starvation, a lipolytic stimulus. We find that loss of bmm largely eliminates the sex difference in triglyceride storage and abolishes the sex difference in triglyceride breakdown via strongly male-biased effects. Although we show that bmm function in the fat body affects whole-body triglyceride levels in both sexes, in males, we identify an additional role for bmm function in the somatic cells of the gonad and in neurons in the regulation of whole-body triglyceride homeostasis. Furthermore, we demonstrate that lipid droplets are normally present in both the somatic cells of the male gonad and in neurons, revealing a previously unrecognized role for bmm function, and possibly lipid droplets, in these cell types in the regulation of whole-body triglyceride homeostasis. Taken together, our data reveal a role for bmm function in the somatic cells of the gonad and in neurons in the regulation of male–female differences in fat storage and breakdown and identify bmm as a link between the regulation of triglyceride homeostasis and biological sex.

An investigation of the genetic and physiological mechanisms underlying sex differences in fat storage and breakdown in the fruit fly Drosophila identifies previously unrecognized sex- and cell type-specific roles for the conserved triglyceride lipase brummer.

BmBlimp-1 gene encoding a C2H2 zinc finger protein is required for wing development in the silkworm Bombyx mori

Cys2-His2 zinc finger (C2H2-ZF) proteins represent the most common class of transcription factors. These factors have great potential for the management of developmental progression by regulating the specific spatiotemporal expression of genes. In this study, we cloned one C2H2-ZF protein gene of Bombyx mori, BGIBMGA000319, that is orthologous to B-lymphocyte-induced maturation protein-1 (Blimp-1); we thus named it as Bombyx mori Blimp-1 (BmBlimp-1). In the silkworm, the BmBlimp-1 gene is specifically upregulated during day 2 of the pupal to adult stage and is highly expressed in wing discs on day 3 of the pupa. Knockdown of its expression level in the pupal stage results in a crumpled-winged silkworm moth. Using the predicted DNA-binding sequences of BmBlimp-1 to search the silkworm genome to screen target genes of BmBlimp-1, 7049 genes were identified to have at least one binding site of BmBlimp-1 on their 1 kb upstream and downstream genome regions. Comparisons of those genes with a reported pupal wing disc transcriptome data resulted in 4065 overlapping genes being retrieved. GO enrichment analysis of the overlapping genes showed that most of the genes were enriched in the binding term. Combining functional annotation and real-time quantitative PCR, 15 genes were identified as the candidate target genes of BmBlimp-1, including several wing cuticular protein genes, chitin synthase A, and wing disc development genes, such as Wnt1, cubitus interruptus (ci) and engrailed (en). Moreover, the amino acid sequence of the zinc finger motif of Blimp-1 gene was highly conserved among the 15 insect species. We propose that BmBlimp-1 is an important regulatory factor in silkworm wing development.

A Multivariate Genome-Wide Association Study of Wing Shape in Drosophila melanogaster

Due to the complexity of genotype-phenotype relationships, simultaneous analyses of genomic associations with multiple traits will be more powerful and informative than a series of univariate analyses. However, in most cases, studies of genotype-phenotype relationships have been analyzed only one trait at a time. Here, we report the results of a fully integrated multivariate genome-wide association analysis of the shape of the Drosophila melanogaster wing in the Drosophila Genetic Reference Panel. Genotypic effects on wing shape were highly correlated between two different laboratories. We found 2396 significant SNPs using a 5% false discovery rate cutoff in the multivariate analyses, but just four significant SNPs in univariate analyses of scores on the first 20 principal component axes. One quarter of these initially significant SNPs retain their effects in regularized models that take into account population structure and linkage disequilibrium. A key advantage of multivariate analysis is that the direction of the estimated phenotypic effect is much more informative than a univariate one. We exploit this fact to show that the effects of knockdowns of genes implicated in the initial screen were on average more similar than expected under a null model. A subset of SNP effects were replicable in an unrelated panel of inbred lines. Association studies that take a phenomic approach, considering many traits simultaneously, are an important complement to the power of genomics.

Comparative Transcriptome Analyses Uncover Key Candidate Genes Mediating Flight Capacity in Bactrocera dorsalis (Hendel) and Bactrocera correcta (Bezzi) (Diptera: Tephritidae)

Flight capacity is important for invasive pests during entry, establishment and spreading. Both Bactroceradorsalis Hendel and Bactroceracorrecta Bezzi are invasive fruit flies but their flight capacities differ. Here, a tethered flight mill test demonstrated that B. dorsalis exhibits a greater flight capacity than B. correcta. RNA-Seq was used to determine the transcriptomic differences associated with the flight capacity of two Bactrocera species. Transcriptome data showed that 6392 unigenes were differentially expressed between the two species in the larval stage, whereas in the adult stage, 4104 differentially expressed genes (DEGs) were identified in females, and 3445 DEGs were observed in males. The flight capacity appeared to be correlated with changes in the transcriptional levels of genes involved in wing formation, flight muscle structure, energy metabolism, and hormonal control. Using RNA interference (RNAi) to verify the function of one DEG, the epidermal growth factor receptor (EGFR), we confirmed the role of this gene in regulating wing development, and thereby flight capacity, in both species. This work reveals the flight mechanism of fruit flies and provides insight into fundamental transcriptomics for further studies on the flight performance of insects.

The Interaction of Genetic Background and Mutational Effects in Regulation of Mouse Craniofacial Shape

Inbred genetic background significantly influences the expression of phenotypes associated with known genetic perturbations and can underlie variation in disease severity between individuals with the same mutation. However, the effect of epistatic interactions on the development of complex traits, such as craniofacial morphology, is poorly understood. Here, we investigated the effect of three inbred backgrounds (129X1/SvJ, C57BL/6J, and FVB/NJ) on the expression of craniofacial dysmorphology in mice (Mus musculus) with loss of function in three members of the Sprouty family of growth factor negative regulators (Spry1, Spry2, or Spry4) in order to explore the impact of epistatic interactions on skull morphology. We found that the interaction of inbred background and the Sprouty genotype explains as much craniofacial shape variation as the Sprouty genotype alone. The most severely affected genotypes display a relatively short and wide skull, a rounded cranial vault, and a more highly angled inferior profile. Our results suggest that the FVB background is more resilient to Sprouty loss of function than either C57 or 129, and that Spry4 loss is generally less severe than loss of Spry1 or Spry2 While the specific modifier genes responsible for these significant background effects remain unknown, our results highlight the value of intercrossing mice of multiple inbred backgrounds to identify the genes and developmental interactions that modulate the severity of craniofacial dysmorphology. Our quantitative results represent an important first step toward elucidating genetic interactions underlying variation in robustness to known genetic perturbations in mice.

Leger RJ. The genetic basis for variation in resistance to infection in the Drosophila melanogaster genetic reference panel

Individuals vary extensively in the way they respond to disease but the genetic basis of this variation is not fully understood. We found substantial individual variation in resistance and tolerance to the fungal pathogen Metarhizium anisopliae Ma549 using the Drosophila melanogaster Genetic Reference Panel (DGRP). In addition, we found that host defense to Ma549 was correlated with defense to the bacterium Pseudomonas aeruginosa Pa14, and several previously published DGRP phenotypes including oxidative stress sensitivity, starvation stress resistance, hemolymph glucose levels, and sleep indices. We identified polymorphisms associated with differences between lines in both their mean survival times and microenvironmental plasticity, suggesting that lines differ in their ability to adapt to variable pathogen exposures. The majority of polymorphisms increasing resistance to Ma549 were sex biased, located in non-coding regions, had moderately large effect and were rare, suggesting that there is a general cost to defense. Nevertheless, host defense was not negatively correlated with overall longevity and fecundity. In contrast to Ma549, minor alleles were concentrated in the most Pa14-susceptible as well as the most Pa14-resistant lines. A pathway based analysis revealed a network of Pa14 and Ma549-resistance genes that are functionally connected through processes that encompass phagocytosis and engulfment, cell mobility, intermediary metabolism, protein phosphorylation, axon guidance, response to DNA damage, and drug metabolism. Functional testing with insertional mutagenesis lines indicates that 12/13 candidate genes tested influence susceptibility to Ma549. Many candidate genes have homologs identified in studies of human disease, suggesting that genes affecting variation in susceptibility are conserved across species.

Differing roles for sur-2/MED23 in C. elegans and C. briggsae vulval development

Normal vulval development in the nematode Caenorhabditis briggsae is identical to that in the related Caenorhabditis elegans. However, several experiments suggest that there are differences between the two species with respect to the contribution of EGF/Ras signaling. To investigate these differences genetically, we have characterized a C. briggsae mutant strain that phenocopies the effect observed when C. briggsae animals are treated with U0126, an inhibitor of the EGF pathway component MEK. We identify that the gene affected in the mutant strain is Cbr-sur-2, which encodes a MED23 mediator complex protein that acts downstream of EGF signaling in C. elegans and other organisms, such as mammals. When Cbr-sur-2 and Cel-sur-2 mutants are compared, we find that the production of additional vulval cells from P5.p and P7.p in C. elegans is dependent on proper development of P6.p, while C. briggsae does not have a similar requirement. Combined chemical and genetic interference with the EGF pathway completely eliminates vulval development in C. elegans but not in C. briggsae. Our results provide genetic evidence for the differing requirements for EGF signaling in the two species.

Disintegrating the fly: A mutational perspective on phenotypic integration and covariation

The structure of environmentally induced phenotypic covariation can influence the effective strength and magnitude of natural selection. Yet our understanding of the factors that contribute to and influence the evolutionary lability of such covariation is poor. Most studies have either examined environmental variation without accounting for covariation, or examined phenotypic and genetic covariation without distinguishing the environmental component. In this study, we examined the effect of mutational perturbations on different properties of environmental covariation, as well as mean shape. We use strains of Drosophila melanogaster bearing well-characterized mutations known to influence wing shape, as well as naturally derived strains, all reared under carefully controlled conditions and with the same genetic background. We find that mean shape changes more freely than the covariance structure, and that different properties of the covariance matrix change independently from each other. The perturbations affect matrix orientation more than they affect matrix eccentricity or total variance. Yet, mutational effects on matrix orientation do not cluster according to the developmental pathway that they target. These results suggest that it might be useful to consider a more general concept of "decanalization," involving all aspects of variation and covariation.

Disruption of snRNP biogenesis factors Tgs1 and pICln induces phenotypes that mirror aspects of SMN-Gemins complex perturbation in Drosophila, providing new insights into spinal muscular atrophy

The neuromuscular disorder, spinal muscular atrophy (SMA), results from insufficient levels of the survival motor neuron (SMN) protein. Together with Gemins 2-8 and Unrip, SMN forms the large macromolecular SMN-Gemins complex, which is known to be indispensable for chaperoning the assembly of spliceosomal small nuclear ribonucleoproteins (snRNPs). It remains unclear whether disruption of this function is responsible for the selective neuromuscular degeneration in SMA. In the present study, we first show that loss of wmd, the Drosophila Unrip orthologue, has a negative impact on the motor system. However, due to lack of a functional relationship between wmd/Unrip and Gemin3, it is likely that Unrip joined the SMN-Gemins complex only recently in evolution. Second, we uncover that disruption of either Tgs1 or pICln, two cardinal players in snRNP biogenesis, results in viability and motor phenotypes that closely resemble those previously uncovered on loss of the constituent members of the SMN-Gemins complex. Interestingly, overexpression of both factors leads to motor dysfunction in Drosophila, a situation analogous to that of Gemin2. Toxicity is conserved in the yeast S. pombe where pICln overexpression induces a surplus of Sm proteins in the cytoplasm, indicating that a block in snRNP biogenesis is partly responsible for this phenotype. Importantly, we show a strong functional relationship and a physical interaction between Gemin3 and either Tgs1 or pICln. We propose that snRNP biogenesis is the pathway connecting the SMN-Gemins complex to a functional neuromuscular system, and its disturbance most likely leads to the motor dysfunction that is typical in SMA.

The sex-limited effects of mutations in the EGFR and TGF-β signaling pathways on shape and size sexual dimorphism and allometry in the Drosophila wing

Much of the morphological diversity in nature-including among sexes within a species-is a direct consequence of variation in size and shape. However, disentangling variation in sexual dimorphism for both shape (SShD), size (SSD), and their relationship with one another remains complex. Understanding how genetic variation influences both size and shape together, and how this in turn influences SSD and SShD, is challenging. In this study, we utilize Drosophila wing size and shape as a model system to investigate how mutations influence size and shape as modulated by sex. Previous work has demonstrated that mutations in epidermal growth factor receptor (EGFR) and transforming growth factor-β (TGF-β) signaling components can influence both wing size and shape. In this study, we re-analyze this data to specifically address how they impact the relationship between size and shape in a sex-specific manner, in turn altering the pattern of sexual dimorphism. While most mutations influence shape overall, only a subset have a genotypic specific effect that influences SShD. Furthermore, while we observe sex-specific patterns of allometric shape variation, the effects of most mutations on allometry tend to be small. We discuss this within the context of using mutational analysis to understand sexual size and shape dimorphism.

C-terminal Src Kinase Gates Homeostatic Synaptic Plasticity and Regulates Fasciclin II Expression at the Drosophila Neuromuscular Junction

Forms of homeostatic plasticity stabilize neuronal outputs and promote physiologically favorable synapse function. A well-studied homeostatic system operates at the Drosophila melanogaster larval neuromuscular junction (NMJ). At the NMJ, impairment of postsynaptic glutamate receptor activity is offset by a compensatory increase in presynaptic neurotransmitter release. We aim to elucidate how this process operates on a molecular level and is preserved throughout development. In this study, we identified a tyrosine kinase-driven signaling system that sustains homeostatic control of NMJ function. We identified C-terminal Src Kinase (Csk) as a potential regulator of synaptic homeostasis through an RNAi- and electrophysiology-based genetic screen. We found that Csk loss-of-function mutations impaired the sustained expression of homeostatic plasticity at the NMJ, without drastically altering synapse growth or baseline neurotransmission. Muscle-specific overexpression of Src Family Kinase (SFK) substrates that are negatively regulated by Csk also impaired NMJ homeostasis. Surprisingly, we found that transgenic Csk-YFP can support homeostatic plasticity at the NMJ when expressed either in the muscle or in the nerve. However, only muscle-expressed Csk-YFP was able to localize to NMJ structures. By immunostaining, we found that Csk mutant NMJs had dysregulated expression of the Neural Cell Adhesion Molecule homolog Fasciclin II (FasII). By immunoblotting, we found that levels of a specific isoform of FasII were decreased in homeostatically challenged GluRIIA mutant animals–but markedly increased in Csk mutant animals. Additionally, we found that postsynaptic overexpression of FasII from its endogenous locus was sufficient to impair synaptic homeostasis, and genetically reducing FasII levels in Csk mutants fully restored synaptic homeostasis. Based on these data, we propose that Csk and its SFK substrates impinge upon homeostatic control of NMJ function by regulating downstream expression or localization of FasII.

Homeostasis is a fundamental topic in biology. Individual cells and systems of cells constantly monitor their environments and adjust their outputs in order to maintain physiological properties within ranges that can support life. The nervous system is no exception. Synapses and circuits are endowed with a capacity to respond to environmental challenges in a homeostatic fashion. As a result, synaptic output stays within an appropriate physiological range. We know that homeostasis is a fundamental form of regulation in animal nervous systems, but we have very little information about how it works. In this study, we examine the fruit fly Drosophila melanogaster and its ability to maintain normal levels of synaptic output over long periods of developmental time. We identify new roles in this process for classical signaling molecules called C-terminal Src kinase, Src family kinases, as well as a neuronal cell adhesion molecule called Fasciclin II, which was previously shown to stabilize synaptic contacts between neurons and muscles. Our work contributes to a broader understanding of how neurons work to maintain stable outputs. Ultimately, this type of knowledge could have important implications for neurological disorders in which stability is lost, such as forms of epilepsy or ataxia.

An image database of Drosophila melanogaster wings for phenomic and biometric analysis

BACKGROUND

Extracting important descriptors and features from images of biological specimens is an ongoing challenge. Features are often defined using landmarks and semi-landmarks that are determined a priori based on criteria such as homology or some other measure of biological significance. An alternative, widely used strategy uses computational pattern recognition, in which features are acquired from the image de novo. Subsets of these features are then selected based on objective criteria. Computational pattern recognition has been extensively developed primarily for the classification of samples into groups, whereas landmark methods have been broadly applied to biological inference.

RESULTS

To compare these approaches and to provide a general community resource, we have constructed an image database of Drosophila melanogaster wings - individually identifiable and organized by sex, genotype and replicate imaging system - for the development and testing of measurement and classification tools for biological images. We have used this database to evaluate the relative performance of current classification strategies. Several supervised parametric and nonparametric machine learning algorithms were used on principal components extracted from geometric morphometric shape data (landmarks and semi-landmarks). For comparison, we also classified phenotypes based on de novo features extracted from wing images using several computer vision and pattern recognition methods as implemented in the Bioimage Classification and Annotation Tool (BioCAT).

CONCLUSIONS

Because we were able to thoroughly evaluate these strategies using the publicly available Drosophila wing database, we believe that this resource will facilitate the development and testing of new tools for the measurement and classification of complex biological phenotypes.

Wing sexual dimorphism of pathogen-vector culicids

Background

Sexual dimorphism in animals has been studied from different perspectives for decades. In 1874 Darwin hypothesized that it was related to sexual selection, and even after nearly 140 years, when additional empirical data has become available and the subject has been investigated from a contemporary viewpoint, this idea is still supported. Although mosquito (Culicidae) wings are of great importance as they play a sex-specific role, little is known about wing sexual dimorphism in these pathogen-vector insects. Detection and characterization of wing sexual dimorphism in culicids may indirectly enhance our knowledge of their epidemiology or reveal sex-linked genes, aspects that have been discussed by vector control initiatives and developers of genetically modified mosquitoes.

Methods

Using geometric morphometrics, we carried out a comparative assessment of wing sexual dimorphism in ten culicid species of medical/veterinary importance from genera Culex, Aedes, Anopheles and Ochlerotatus collected in Brazil.

Results

Discriminant analysis revealed significant sexual dimorphism in all the species studied, indicating that phenotypic expression of wing shape in mosquitoes is indeed sex-specific. A cross-validated test performed to reclassify the sexes with and without allometry yielded very similar results. Mahalanobis distances among the ten species showed that the species had different patterns of shape sexual dimorphism and that females are larger than males in some species.

Conclusion

Wing morphology differed significantly between species. The finding of sexual dimorphism in all the species would suggest that the wing geometry of Culicidae is canalized. Although sexual dimorphism is prevalent, species-specific patterns occur. Allometry was not the main determinant of sexual dimorphism, which suggests that sexual selection or other evolutionary mechanisms underlie wing sexual dimorphism in these insects.

Causes and consequences of genetic background effects illuminated by integrative genomic analysis

The phenotypic consequences of individual mutations are modulated by the wild-type genetic background in which they occur. Although such background dependence is widely observed, we do not know whether general patterns across species and traits exist or about the mechanisms underlying it. We also lack knowledge on how mutations interact with genetic background to influence gene expression and how this in turn mediates mutant phenotypes. Furthermore, how genetic background influences patterns of epistasis remains unclear. To investigate the genetic basis and genomic consequences of genetic background dependence of the scalloped(E3) allele on the Drosophila melanogaster wing, we generated multiple novel genome-level datasets from a mapping-by-introgression experiment and a tagged RNA gene expression dataset. In addition we used whole genome resequencing of the parental lines-two commonly used laboratory strains-to predict polymorphic transcription factor binding sites for SD. We integrated these data with previously published genomic datasets from expression microarrays and a modifier mutation screen. By searching for genes showing a congruent signal across multiple datasets, we were able to identify a robust set of candidate loci contributing to the background-dependent effects of mutations in sd. We also show that the majority of background-dependent modifiers previously reported are caused by higher-order epistasis, not quantitative noncomplementation. These findings provide a useful foundation for more detailed investigations of genetic background dependence in this system, and this approach is likely to prove useful in exploring the genetic basis of other traits as well.

Genetic complexity in a Drosophila model of diabetes-associated misfolded human proinsulin

Drosophila melanogaster has been widely used as a model of human Mendelian disease, but its value in modeling complex disease has received little attention. Fly models of complex disease would enable high-resolution mapping of disease-modifying loci and the identification of novel targets for therapeutic intervention. Here, we describe a fly model of permanent neonatal diabetes mellitus and explore the complexity of this model. The approach involves the transgenic expression of a misfolded mutant of human preproinsulin, hINS^C96Y, which is a cause of permanent neonatal diabetes. When expressed in fly imaginal discs, hINS^C96Y causes a reduction of adult structures, including the eye, wing, and notum. Eye imaginal discs exhibit defects in both the structure and the arrangement of ommatidia. In the wing, expression of hINS^C96Y leads to ectopic expression of veins and mechano-sensory organs, indicating disruption of wild-type signaling processes regulating cell fates. These readily measurable “disease” phenotypes are sensitive to temperature, gene dose, and sex. Mutant (but not wild-type) proinsulin expression in the eye imaginal disc induces IRE1-mediated XBP1 alternative splicing, a signal for endoplasmic reticulum stress response activation, and produces global change in gene expression. Mutant hINS transgene tester strains, when crossed to stocks from the Drosophila Genetic Reference Panel, produce F₁ adults with a continuous range of disease phenotypes and large broad-sense heritability. Surprisingly, the severity of mutant hINS-induced disease in the eye is not correlated with that in the notum in these crosses, nor with eye reduction phenotypes caused by the expression of two dominant eye mutants acting in two different eye development pathways, Drop (Dr) or Lobe (L), when crossed into the same genetic backgrounds. The tissue specificity of genetic variability for mutant hINS-induced disease has, therefore, its own distinct signature. The genetic dominance of disease-specific phenotypic variability in our model of misfolded human proinsulin makes this approach amenable to genome-wide association study in a simple F₁ screen of natural variation.

A perspective on micro-evo-devo: progress and potential

The term "micro-evo-devo" refers to the combined study of the genetic and developmental bases of natural variation in populations and the evolutionary forces that have shaped this variation. It thus represents a synthesis of the fields of evolutionary developmental biology and population genetics. As has been pointed out by several others, this synthesis can provide insights into the evolution of organismal form and function that have not been possible within these individual disciplines separately. Despite a number of important successes in micro-evo-devo, however, it appears that evo devo and population genetics remain largely separate spheres of research, limiting their ability to address evolutionary questions. This also risks pushing contemporary evo devo to the fringes of evolutionary biology because it does not describe the causative molecular changes underlying evolution or the evolutionary forces involved. Here we reemphasize the theoretical and practical importance of micro-evo-devo as a strategy for understanding phenotypic evolution, review the key recent insights that it has provided, and present a perspective on both the potential and the remaining challenges of this exciting interdisciplinary field.

The conditional nature of genetic interactions: the consequences of wild-type backgrounds on mutational interactions in a genome-wide modifier screen

The phenotypic outcome of a mutation cannot be simply mapped onto the underlying DNA variant. Instead, the phenotype is a function of the allele, the genetic background in which it occurs and the environment where the mutational effects are expressed. While the influence of genetic background on the expressivity of individual mutations is recognized, its consequences on the interactions between genes, or the genetic network they form, is largely unknown. The description of genetic networks is essential for much of biology; yet if, and how, the topologies of such networks are influenced by background is unknown. Furthermore, a comprehensive examination of the background dependent nature of genetic interactions may lead to identification of novel modifiers of biological processes. Previous work in Drosophila melanogaster demonstrated that wild-type genetic background influences the effects of an allele of scalloped (sd), with respect to both its principal consequence on wing development and its interactions with a mutation in optomotor blind. In this study we address whether the background dependence of mutational interactions is a general property of genetic systems by performing a genome wide dominant modifier screen of the sd(E3) allele in two wild-type genetic backgrounds using molecularly defined deletions. We demonstrate that ~74% of all modifiers of the sd(E3) phenotype are background-dependent due in part to differential sensitivity to genetic perturbation. These background dependent interactions include some with qualitative differences in the phenotypic outcome, as well as instances of sign epistasis. This suggests that genetic interactions are often contingent on genetic background, with flexibility in genetic networks due to segregating variation in populations. Such background dependent effects can substantially alter conclusions about how genes influence biological processes, the potential for genetic screens in alternative wild-type backgrounds identifying new loci that contribute to trait expression, and the inferences of the topology of genetic networks.

Gene-by-temperature interactions and candidate plasticity genes for morphological traits in Drosophila melanogaster

Understanding the genetic architecture of any quantitative trait requires identifying the genes involved in its expression in different environmental conditions. This goal can be achieved by mutagenesis screens in genetically tractable model organisms such as Drosophila melanogaster. Temperature during ontogenesis is an important environmental factor affecting development and phenotypic variation in holometabolous insects. In spite of the importance of phenotypic plasticity and genotype by environment interaction (GEI) for fitness related traits, its genetic basis has remained elusive. In this context, we analyzed five different adult morphological traits (face width, head width, thorax length, wing size and wing shape) in 42 co-isogenic single P-element insertional lines of Drosophila melanogaster raised at 17°C and 25°C. Our analyses showed that all lines differed from the control for at least one trait in males or females at either temperature. However, no line showed those differences for all traits in both sexes and temperatures simultaneously. In this sense, the most pleiotropic candidate genes were CG34460, Lsd-2 and Spn. Our analyses also revealed extensive genetic variation for all the characters mostly indicated by strong GEIs. Further, our results indicate that GEIs were predominantly explained by changes in ranking order in all cases suggesting that a moderate number of genes are involved in the expression of each character at both temperatures. Most lines displayed a plastic response for at least one trait in either sex. In this regard, P-element insertions affecting plasticity of a large number of traits were associated to the candidate genes Btk29A, CG43340, Drak and jim. Further studies will help to elucidate the relevance of these genes on the morphogenesis of different body structures in natural populations of D. melanogaster.

Specific-gene studies of evolutionary mechanisms in an age of genome-wide surveying

The molecular tools of genomics have great power to reveal patterns of genetic difference within or among species, but must be complemented by the mechanistic study of the genetic variants found if these variants' evolutionary meaning is to be well understood. Central to this purpose is knowledge of the organisms' genotype-phenotype-environment interactions, which embody biological adaptation and constraint and thus drive natural selection. The history of this approach is briefly reviewed. Strategies embracing the complementarity of genomics and specific-gene studies in evolution are considered. Implementation of these strategies, and examples showing their feasibility and power, are discussed. Initial generalizations emphasize: (1) reproducibility of adaptive mechanisms; (2) evolutionary co-importance of variation in protein sequences and expression; (3) refinement of rudimentary molecular functions as an origin of evolutionary innovations; (4) identification of specific-gene mechanisms as underpinnings of genomic or quantitative genetic variation; and (5) multiple forms of adaptive or constraining epistasis among genes. Progress along these lines will advance understanding of evolution and support its use in addressing urgent medical and environmental applications.

Does your gene need a background check? How genetic background impacts the analysis of mutations, genes, and evolution

The premise of genetic analysis is that a causal link exists between phenotypic and allelic variation. However, it has long been documented that mutant phenotypes are not a simple result of a single DNA lesion, but are instead due to interactions of the focal allele with other genes and the environment. Although an experimentally rigorous approach focused on individual mutations and isogenic control strains has facilitated amazing progress within genetics and related fields, a glimpse back suggests that a vast complexity has been omitted from our current understanding of allelic effects. Armed with traditional genetic analyses and the foundational knowledge they have provided, we argue that the time and tools are ripe to return to the underexplored aspects of gene function and embrace the context-dependent nature of genetic effects. We assert that a broad understanding of genetic effects and the evolutionary dynamics of alleles requires identifying how mutational outcomes depend upon the 'wild type' genetic background. Furthermore, we discuss how best to exploit genetic background effects to broaden genetic research programs.

Shape variation in outline shapes

A general morphometric method for describing shape variation in a sample consisting of landmarks and multiple outline shapes is developed in this article. A distance metric is developed for such data and is used to embed the data in a low-dimensional Euclidean space. The Euclidean space is used to generate summary statistics such as mean and principal shape variation which are implicitly represented in the original space using elements of the sample. A new distance metric for outline shapes is proposed based on Procrustes distance that does not require the extraction of discrete points along the curve. The outline distance metric can be naturally combined with distances between landmarks. A method for aligning outlines and multiple outlines is developed that minimizes the distance metric. The method is compared with semilandmarks on synthetic data and 2 real data sets. Outline methods produce useful and valid results when suitably constrained by landmarks and are useful visualization aids, but questions remain about their suitability for answering biological questions until appropriate distance metrics can be biologically validated.

Multiple capacitors for natural genetic variation in Drosophila melanogaster

Cryptic genetic variation (CGV) or a standing genetic variation that is not ordinarily expressed as a phenotype is released when the robustness of organisms is impaired under environmental or genetic perturbations. Evolutionary capacitors modulate the amount of genetic variation exposed to natural selection and hidden cryptically; they have a fundamental effect on the evolvability of traits on evolutionary timescales. In this study, I have demonstrated the effects of multiple genomic regions of Drosophila melanogaster on CGV in wing shape. I examined the effects of 61 genomic deficiencies on quantitative and qualitative natural genetic variation in the wing shape of D. melanogaster. I have identified 10 genomic deficiencies that do not encompass a known candidate evolutionary capacitor, Hsp90, exposing natural CGV differently depending on the location of the deficiencies in the genome. Furthermore, five genomic deficiencies uncovered qualitative CGV in wing morphology. These findings suggest that CGV in wing shape of wild-type D. melanogaster is regulated by multiple capacitors with divergent functions. Future analysis of genes encompassed by these genomic regions would help elucidate novel capacitor genes and better understand the general features of capacitors regarding natural genetic variation.

The Measurement of Local Variation in Shape

Geometric morphometrics comprises tools for measuring and analyzing shape as captured by an entire set of landmark configurations. Many interesting questions in evolutionary, genetic, and developmental research, however, are only meaningful at a local level, where a focus on "parts" or "traits" takes priority over properties of wholes. To study variational properties of such traits, current approaches partition configurations into subsets of landmarks which are then studied separately. This approach is unable to fully capture both variational and spatial characteristics of these subsets because interpretability of shape differences is context-dependent. Landmarks omitted from a partition usually contain information about that partition's shape. We present an interpolation-based approach that can be used to model shape differences at a local, infinitesimal level as a function of information available globally. This approach belongs in a large family of methods that see shape differences as continuous "fields" spanning an entire structure, for which landmarks serve as reference parameters rather than as data. We show, via analyses of simulated and real data, how interpolation models provide a more accurate representation of regional shapes than partitioned data. A key difference of this interpolation approach from current morphometric practice is that one must assume an explicit interpolation model, which in turn implies a particular kind of behavior of the regions between landmarks. This choice presents novel methodological challenges, but also an opportunity to incorporate and test biomechanical models that have sought to explain tissue-level processes underlying the generation of morphological shape.

Wg signaling via Zw3 and mad restricts self-renewal of sensory organ precursor cells in Drosophila

It is well known that the Dpp signal transducer Mad is activated by phosphorylation at its carboxy-terminus. The role of phosphorylation on other regions of Mad is not as well understood. Here we report that the phosphorylation of Mad in the linker region by the Wg antagonist Zw3 (homolog of vertebrate Gsk3-β) regulates the development of sensory organs in the anterior-dorsal quadrant of the wing. Proneural expression of Mad-RNA interference (RNAi) or a Mad transgene with its Zw3/Gsk3-β phosphorylation sites mutated (MGM) generated wings with ectopic sensilla and chemosensory bristle duplications. Studies with pMad-Gsk (an antibody specific to Zw3/Gsk3-β-phosphorylated Mad) in larval wing disks revealed that this phosphorylation event is Wg dependent (via an unconventional mechanism), is restricted to anterior-dorsal sensory organ precursors (SOP) expressing Senseless (Sens), and is always co-expressed with the mitotic marker phospho-histone3. Quantitative analysis in both Mad-RNAi and MGM larval wing disks revealed a significant increase in the number of Sens SOP. We conclude that the phosphorylation of Mad by Zw3 functions to prevent the self-renewal of Sens SOP, perhaps facilitating their differentiation via asymmetric division. The conservation of Zw3/Gsk3-β phosphorylation sites in vertebrate homologs of Mad (Smads) suggests that this pathway, the first transforming growth factor β-independent role for any Smad protein, may be widely utilized for regulating mitosis during development.

Genetic basis of wing morphogenesis in Drosophila: sexual dimorphism and non-allometric effects of shape variation

BACKGROUND

The Drosophila wing represents a particularly appropriate model to investigate the developmental control of phenotypic variation. Previous studies which aimed to identify candidate genes for wing morphology demonstrated that the genetic basis of wing shape variation in D. melanogaster is composed of numerous genetic factors causing small, additive effects. In this study, we analyzed wing shape in males and females from 191 lines of D. melanogaster, homozygous for a single P-element insertion, using geometric morphometrics techniques. The analysis allowed us to identify known and novel candidate genes that may contribute to the expression of wing shape in each sex separately and to compare them to candidate genes affecting wing size which have been identified previously using the same lines.

RESULTS

Our results indicate that more than 63% of induced mutations affected wing shape in one or both sexes, although only 33% showed significant differences in both males and females. The joint analysis of wing size and shape revealed that only 19% of the P-element insertions caused coincident effects on both components of wing form in one or both sexes. Further morphometrical analyses revealed that the intersection between veins showed the smallest displacements in the proximal region of the wing. Finally, we observed that mutations causing general deformations were more common than expected in both sexes whereas the opposite occurred with those generating local changes. For most of the 94 candidate genes identified, this seems to be the first record relating them with wing shape variation.

CONCLUSIONS

Our results support the idea that the genetic architecture of wing shape is complex with many different genes contributing to the trait in a sexually dimorphic manner. This polygenic basis, which is relatively independent from that of wing size, is composed of genes generally involved in development and/or metabolic functions, especially related to the regulation of different cellular processes such as motility, adhesion, communication and signal transduction. This study suggests that understanding the genetic basis of wing shape requires merging the regulation of vein patterning by signalling pathways with processes that occur during wing development at the cellular level.

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.