Search

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., (Copyright © 2019 by the Genetics Society of America.)

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.

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.

Genetic background effects contribute to the phenotypic consequences of mutations and are pervasive across all domains of life that have been examined, yet little is known about how they modify genetic systems. In part this is due to the lack of tractable model systems that have been explicitly developed to study the genetic and evolutionary consequences of background effects. In this study we demonstrate that phenotypic expressivity of the scalloped(E3) (sd(E3)) mutation of Drosophila melanogaster is background dependent and is the result of at least one major modifier segregating between two standard lab wild-type strains. We provide evidence that at least one of the modifiers is linked to the vestigial region and demonstrate that the background effects modify the spatial distribution of known sd target genes in a genotype-dependent manner. In addition, microarrays were used to examine the consequences of genetic background effects on the global transcriptome. Expression differences between wild-type strains were found to be as large as or larger than the effects of mutations with substantial phenotypic effects, and expression differences between wild type and mutant varied significantly between genetic backgrounds. Significantly, we demonstrate that the epistatic interaction between sd(E3) and an optomotor blind mutation is background dependent. The results are discussed within the context of developing a complex but more realistic view of the consequences of genetic background effects with respect to mutational analysis and studies of epistasis and cryptic genetic variation segregating in natural populations.

The putative regulatory relationships between Antennapedia (Antp), spalt major (salm) and homothorax (hth) are tested with regard to the sensitive period of antenna-to-leg transformations. Although Antp expression repressed hth as predicted, contrary to expectations, hth did not show increased repression at higher Antp doses, whereas salm, a gene downstream of hth, did show such a dose response. Loss of hth allowed antenna-to-leg transformations but the relative timing of proximal-distal transformations was reversed, relative to transformations induced by ectopic Antp. Finally, overexpression of Hth was only partially able to rescue transformations induced by ectopic Antp. These results indicate that there may be additional molecules involved in antenna/leg identity and that spatial, temporal and dosage relationships are more subtle than suspected and must be part of a robust understanding of molecular network behaviour involved in determining appendage identity in Drosophila melanogaster.

Wing development in Drosophila is a common model system for the dissection of genetic networks and their roles during development. In particular, the RTK and TGF-beta regulatory networks appear to be involved with numerous aspects of wing development, including patterning, cell determination, growth, proliferation, and survival in the developing imaginal wing disc. However, little is known as to how subtle changes in the function of these genes may contribute to quantitative variation for wing shape, per se. In this study 50 insertional mutations, representing 43 loci in the RTK, Hedgehog, TGF-beta pathways, and their genetically interacting factors were used to study the role of these networks on wing shape. To concurrently examine how genetic background modulates the effects of the mutation, each insertion was introgressed into two wild-type genetic backgrounds. Using geometric morphometric methods, it is shown that the majority of these mutations have profound effects on shape but not size of the wing when measured as heterozygotes. To examine the relationships between how each mutation affects wing shape hierarchical clustering was used. Unlike previous observations of environmental canalization, these mutations did not generally increase within-line variation relative to their wild-type counterparts. These results provide an entry point into the genetics of wing shape and are discussed within the framework of the dissection of complex phenotypes.

Background: Quantitative differences between individuals stem from a combination of genetic and environmental factors, with the heritable variation being shaped by evolutionary forces. Drosophila wing shape has emerged as an attractive system for genetic dissection of multi-dimensional traits. We utilize several experimental genetic methods to validation of the contribution of several polymorphisms in the Epidermal growth factor receptor (Egfr) gene to wing shape and size, that were previously mapped in populations of Drosophila melanogaster from North Carolina (NC) and California (CA). This re-evaluation utilized different genetic testcrosses to generate heterozygous individuals with a variety of genetic backgrounds as well as sampling of new alleles from Kenyan stocks., Results: Only one variant, in the Egfr promoter, had replicable effects in all new experiments. However, expanded genotyping of the initial sample of inbred lines rendered the association non-significant in the CA population, while it persisted in the NC sample, suggesting population specific modification of the quantitative trait nucleotide QTN effect., Conclusion: Dissection of quantitative trait variation to the nucleotide level can identify sites with replicable effects as small as one percent of the segregating genetic variation. However, the testcross approach to validate QTNs is both labor intensive and time-consuming, and is probably less useful than resampling of large independent sets of outbred individuals.

The Epidermal growth factor receptor is an essential gene with diverse pleiotropic roles in development throughout the animal kingdom. Analysis of sequence diversity in 10.9 kb covering the complete coding region and 6.4 kb of potential regulatory regions in a sample of 250 alleles from three populations of Drosophila melanogaster suggests that the intensity of different population genetic forces varies along the locus. A total of 238 independent common SNPs and 20 indel polymorphisms were detected, with just six common replacements affecting >1475 amino acids, four of which are in the short alternate first exon. Sequence diversity is lowest in a 2-kb portion of intron 2, which is also highly conserved in comparison with D. simulans and D. pseudoobscura. Linkage disequilibrium decays to background levels within 500 bp of most sites, so haplotypes are generally restricted to up to 5 polymorphisms. The two North American samples from North Carolina and California have diverged in allele frequency at a handful of individual SNPs, but a Kenyan sample is both more divergent and more polymorphic. The effect of sample size on inference of the roles of population structure, uneven recombination, and weak selection in patterning nucleotide variation in the locus is discussed.

The morphogenetic field, a fundamental concept of classical embryology, is once again being invoked to describe developmental processes. Because the evolution of adult structures requires the modification of development, the ways in which morphogenetic fields can change over time may yield insights into evolutionary possibilities. We considered how the duplication/multiplication of a morphogenetic field in fruit flies, caused by the previously described obake (obk) mutation, is regulated by genetic and environmental factors. Mutations of genes in the canonical antenna-producing imaginal disc pathway suppressed duplication as expected, although the results suggested that other pathways might also be involved. Overgrowth mutations, expected to increase duplication, actually suppressed it. Mutations in the heat-shock protein gene Hsp83 did not uniformly enhance obk expressivity as hypothesized. Using third chromosomes extracted from wild-derived lines, natural genetic variation for modifiers of obk function was found to be extensive. Larval crowding suppressed the obk phenotype, but there was no evidence of trade-offs between body or head size and arista number. Our results suggest that a complex interplay of genetic and environmental factors in the regulation of fields may be responsible for ample natural variation in the expressivity of adult phenotypes, affording multiple opportunities for selection and evolutionary modification.