Nucleotide variation in the Egfr locus of Drosophila melanogaster

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.

Two Drosophilids exhibit distinct EGF pathway patterns in oogenesis

Deciphering the evolution of morphological structures is a remaining challenge in the field of developmental biology. The respiratory structures of insect eggshells, called the dorsal appendages, provide an outstanding system for exploring these processes since considerable information is known about their patterning and morphogenesis in Drosophila melanogaster and dorsal appendage number and morphology vary widely across Drosophilid species. We investigated the patterning differences that might facilitate morphogenetic differences between D. melanogaster, which produces two oar-like structures first by wrapping and then elongating the tubes via cell intercalation and cell crawling, and Scaptodrosophila lebanonensis, which produces a variable number of appendages simply by cell intercalation and crawling. Analyses of BMP pathway components thickveins and P-Mad demonstrate that anterior patterning is conserved between these species. In contrast, EGF signaling exhibits significant differences. Transcripts for the ligand encoded by gurken localize similarly in the two species, but this morphogen creates a single dorsolateral primordium in S. lebanonensis as defined by activated MAP kinase and the downstream marker broad. Expression patterns of pointed, argos, and Capicua, early steps in the EGF pathway, exhibit a heterochronic shift in S. lebanonensis relative to those seen in D. melanogaster. We demonstrate that the S. lebanonensis Gurken homolog is active in D. melanogaster but is insufficient to alter downstream patterning responses, indicating that Gurken-EGF receptor interactions do not distinguish the two species' patterning. Altogether, these results differentiate EGF signaling patterns between species and shed light on how changes to the regulation of patterning genes may contribute to different tube-forming mechanisms.

Naturally occurring deletions of hunchback binding sites in the even-skipped stripe 3+7 enhancer

Changes in regulatory DNA contribute to phenotypic differences within and between taxa. Comparative studies show that many transcription factor binding sites (TFBS) are conserved between species whereas functional studies reveal that some mutations segregating within species alter TFBS function. Consistently, in this analysis of 13 regulatory elements in Drosophila melanogaster populations, single base and insertion/deletion polymorphism are rare in characterized regulatory elements. Experimentally defined TFBS are nearly devoid of segregating mutations and, as has been shown before, are quite conserved. For instance 8 of 11 Hunchback binding sites in the stripe 3+7 enhancer of even-skipped are conserved between D. melanogaster and Drosophila virilis. Oddly, we found a 72 bp deletion that removes one of these binding sites (Hb8), segregating within D. melanogaster. Furthermore, a 45 bp deletion polymorphism in the spacer between the stripe 3+7 and stripe 2 enhancers, removes another predicted Hunchback site. These two deletions are separated by ∼250 bp, sit on distinct haplotypes, and segregate at appreciable frequency. The Hb8Δ is at 5 to 35% frequency in the new world, but also shows cosmopolitan distribution. There is depletion of sequence variation on the Hb8Δ-carrying haplotype. Quantitative genetic tests indicate that Hb8Δ affects developmental time, but not viability of offspring. The Eve expression pattern differs between inbred lines, but the stripe 3 and 7 boundaries seem unaffected by Hb8Δ. The data reveal segregating variation in regulatory elements, which may reflect evolutionary turnover of characterized TFBS due to drift or co-evolution.

Recombination modulates how selection affects linked sites in Drosophila

Recombination rate in Drosophila species shapes the impact of selection in the genome and is positively correlated with nucleotide diversity.

One of the most influential observations in molecular evolution has been a strong association between local recombination rate and nucleotide polymorphisms across the genome. This is interpreted as evidence for ubiquitous natural selection. The alternative explanation, that recombination is mutagenic, has been rejected by the absence of a similar association between local recombination rate and nucleotide divergence between species. However, many recent studies show that recombination rates are often very different even in closely related species, questioning whether an association between recombination rate and divergence between species has been tested satisfactorily. To circumvent this problem, we directly surveyed recombination across approximately 43% of the D. pseudoobscura physical genome in two separate recombination maps and 31% of the D. miranda physical genome, and we identified both global and local differences in recombination rate between these two closely related species. Using only regions with conserved recombination rates between and within species and accounting for multiple covariates, our data support the conclusion that recombination is positively related to diversity because recombination modulates Hill–Robertson effects in the genome and not because recombination is predominately mutagenic. Finally, we find evidence for dips in diversity around nonsynonymous substitutions. We infer that at least some of this reduction in diversity resulted from selective sweeps and examine these dips in the context of recombination rate.

Individuals within a species differ in the DNA sequences of their genes. This sequence variation affects how well individuals survive or reproduce and is transmitted to their offspring. Genes near each other on individual chromosomes tend to be passed to offspring together—neighboring genes are unlikely to be separated by exchanges of genetic material derived from different parents during meiotic recombination. When genes are inherited together, however, the evolutionary forces acting on one gene can interfere with variation at its neighbors. Thus, variation at multiple genes can be lost if natural selection acts on one gene in close proximity. Recombination can prevent or reduce this loss of variation, but previous tests of this phenomenon failed to account for recombination rate differences between species. In this study, we show that some parts of the genome differ in recombination rate between two species of fruit fly, Drosophila pseudoobscura and D. miranda. Avoiding an assumption made in previous studies, we then examine sequence variation within and between fly species in those parts of the genome that have conserved recombination rates. Based on the results, we conclude that recombination indeed preserves variation within species that would otherwise have been eliminated by natural selection.

Genomic differentiation between temperate and tropical Australian populations of Drosophila melanogaster

Determining the genetic basis of environmental adaptation is a central problem of evolutionary biology. This issue has been fruitfully addressed by examining genetic differentiation between populations that are recently separated and/or experience high rates of gene flow. A good example of this approach is the decades-long investigation of selection acting along latitudinal clines in Drosophila melanogaster. Here we use next-generation genome sequencing to reexamine the well-studied Australian D. melanogaster cline. We find evidence for extensive differentiation between temperate and tropical populations, with regulatory regions and unannotated regions showing particularly high levels of differentiation. Although the physical genomic scale of geographic differentiation is small--on the order of gene sized--we observed several larger highly differentiated regions. The region spanned by the cosmopolitan inversion polymorphism In(3R)P shows higher levels of differentiation, consistent with the major difference in allele frequencies of Standard and In(3R)P karyotypes in temperate vs. tropical Australian populations. Our analysis reveals evidence for spatially varying selection on a number of key biological processes, suggesting fundamental biological differences between flies from these two geographic regions.

Estimation of fine-scale recombination intensity variation in the white-echinus interval of D. melanogaster

Accurate assessment of local recombination rate variation is crucial for understanding the recombination process and for determining the impact of natural selection on linked sites. In Drosophila, local recombination intensity has been estimated primarily by statistical approaches, by estimating the local slope of the relationship between the physical and genetic maps. However, these estimates are limited in resolution and, as a result, the physical scale at which recombination intensity varies in Drosophila is largely unknown. Although there is some evidence suggesting as much as a 40-fold variation in crossover rate at a local scale in D. pseudoobscura, little is known about the fine-scale structure of recombination rate variation in D. melanogaster. Here we experimentally examine the fine-scale distribution of crossover events in a 1.2-Mb region on the D. melanogaster X chromosome using a classic genetic mapping approach. Our results show that crossover frequency is significantly heterogeneous within this region, varying approximately 3.5-fold. Simulations suggest that this degree of heterogeneity is sufficient to affect levels of standing nucleotide diversity, although the magnitude of this effect is small. We recover no statistical association between empirical estimates of nucleotide diversity and recombination intensity, which is likely due to the limited number of loci sampled in our population genetic data set. However, codon bias is significantly negatively correlated with fine-scale recombination intensity estimates, as expected. Our results shed light on the relevant physical scale to consider in evolutionary analyses relating to recombination rate and highlight the motivations to increase the resolution of the recombination map in Drosophila.

Recombinational landscape and population genomics of Caenorhabditis elegans

Recombination rate and linkage disequilibrium, the latter a function of population genomic processes, are the critical parameters for mapping by linkage and association, and their patterns in Caenorhabditis elegans are poorly understood. We performed high-density SNP genotyping on a large panel of recombinant inbred advanced intercross lines (RIAILs) of C. elegans to characterize the landscape of recombination and, on a panel of wild strains, to characterize population genomic patterns. We confirmed that C. elegans autosomes exhibit discrete domains of nearly constant recombination rate, and we show, for the first time, that the pattern holds for the X chromosome as well. The terminal domains of each chromosome, spanning about 7% of the genome, exhibit effectively no recombination. The RIAILs exhibit a 5.3-fold expansion of the genetic map. With median marker spacing of 61 kb, they are a powerful resource for mapping quantitative trait loci in C. elegans. Among 125 wild isolates, we identified only 41 distinct haplotypes. The patterns of genotypic similarity suggest that some presumed wild strains are laboratory contaminants. The Hawaiian strain, CB4856, exhibits genetic isolation from the remainder of the global population, whose members exhibit ample evidence of intercrossing and recombining. The population effective recombination rate, estimated from the pattern of linkage disequilibrium, is correlated with the estimated meiotic recombination rate, but its magnitude implies that the effective rate of outcrossing is extremely low, corroborating reports of selection against recombinant genotypes. Despite the low population, effective recombination rate and extensive linkage disequilibrium among chromosomes, which are techniques that account for background levels of genomic similarity, permit association mapping in wild C. elegans strains.

C. elegans is a model system for diverse fields of biology, but its ability to serve as a model for quantitative trait gene mapping depends on its recombination rate in the laboratory and in nature. The latter is a function of how worms mate and migrate in the wild. We examined the patterns of recombination in a population that we put through thousands of meioses in the laboratory and in a collection of strains isolated from nature. The data suggest that meiotic recombination rate is highly regular in worms, with discrete domains whose boundaries we identify. The pattern in natural strains suggests that population structure, population size, outcrossing rate, and selection combine to suppress the overall effects of recombination. Moreover, some “wild” strains appear to be laboratory contaminants. Nevertheless, the history of recombination in wild worms is sufficient to permit correlations between genotype and phenotype to pinpoint the loci responsible for phenotypic variation.

Mixture modeling of transcript abundance classes in natural populations

Expression profiling of Drosophila melanogaster adult female heads for 108 nearly isogenic lines from two different populations, and of CEPH lymphoblastoid lines, shows that differential expression of transcripts among individuals is due to a complex interplay of cis- and trans-acting factors.

Background

Populations diverge in genotype and phenotype under the influence of such evolutionary processes as genetic drift, mutation accumulation, and natural selection. Because genotype maps onto phenotype by way of transcription, it is of interest to evaluate how these evolutionary factors influence the structure of variation at the level of transcription. Here, we explore the distributions of cis-acting and trans-acting factors and their relative contributions to expression of transcripts that exhibit two or more classes of abundance among individuals within populations.

Results

Expression profiling using cDNA microarrays was conducted in Drosophila melanogaster adult female heads for 58 nearly isogenic lines from a North Carolina population and 50 from a California population. Using a mixture modeling approach, transcripts were identified that exhibit more than one mode of transcript abundance across the samples. Power studies indicate that sample sizes of 50 individuals will generally be sufficient to detect divergent transcript abundance classes. The distribution of transcript abundance classes is skewed toward low frequency minor classes, which is reminiscent of the typical skew in genotype frequencies. Similar results are observed in reported data on gene expression in human lymphoblast cell lines, in which analysis of association with linked polymorphisms implies that cis-acting single nucleotide polymorphisms make only a modest contribution to bimodal distributions of transcript abundance.

Conclusion

Population surveys of gene expression may complement genetical genomics as a general approach to quantifying sources of transcriptional variation. Differential expression of transcripts among individuals is due to a complex interplay of cis-acting and trans-acting factors.

Genetic variation for cardiac dysfunction in Drosophila

BACKGROUND

Common diseases may be attributed to combinations of variant alleles, but there are few model systems where the interactions among such variants can be studied in controlled genetic crosses. While association studies are designed to detect common polymorphisms of moderate effect, new approaches are required to characterize the impact on disease of interactions among rare alleles.

METHODOLOGY/PRINCIPAL FINDINGS

We show that wild populations of Drosophila melanogaster harbor rare polymorphisms of major effect (RAME) that predispose flies to a specific disease phenotype, age-dependent cardiac dysfunction. A screen of fifty inbred wild-type lines revealed a continuous spectrum of pacing-induced heart failure that generally increases in frequency with age. High-speed video analysis of the inbred lines with high rates of inducible heart failure indicates specific defects in cardiac function, including arrhythmias and contractile disorders ('cardiomyopathies'). A combination of bulked segregant analysis and single feature polymorphism (SFP) detection localizes one of the cardiac susceptibility loci to the 97C interval on the fly genome.

CONCLUSIONS/SIGNIFICANCE

Wild-type Drosophila, like humans, are predisposed to cardiac dysfunction. Identification of factors associated with these naturally occurring cardiac traits promises to provide important insights into the epidemiology of cardiac disease.

How repeatable are associations between polymorphisms in achaete-scute and bristle number variation in Drosophila?

Currently, the relevance of common genetic variants--particularly those significantly associated with phenotypic variation in laboratory studies--to standing phenotypic variation in the wild is poorly understood. To address this, we quantified the relationship between achaete-scute complex (ASC) polymorphisms and Drosophila bristle number phenotypes in several new population samples. MC22 is a biallelic, nonrepetitive-length polymorphism 97 bp downstream of the scute transcript. It has been previously shown to be associated with sternopleural bristle number variation in both sexes in a set of isogenic lines. We replicated this association in a large cohort of wild-caught Drosophila melanogaster. We also detected a significant association at MC22 in an outbred population maintained under laboratory conditions for approximately 25 years, but the phenotypic effects in this sample were opposite from the direction estimated in the initial study. Finally, no significant associations were detected in a second large wild-caught cohort or in a set of 134 nearly isogenic lines. Our ability to repeat the initial association in wild samples suggests that it was not spurious. Nevertheless, inconsistent results from the other three panels suggest that the relationship between polymorphic genetic markers and loci contributing to continuous variation is not a simple one.

Quantitative trait transcripts for nicotine resistance in Drosophila melanogaster

Although most genetic association studies are performed with the intention of detecting nucleotide polymorphisms that are correlated with a complex trait, transcript abundance should also be expected to associate with diseases or phenotypes. We performed a scan for such quantitative trait transcripts in adult female heads of the fruit fly (Drosophila melanogaster) that might explain variation for nicotine resistance. The strongest association was seen for abundance of ornithine aminotransferase transcripts, implicating detoxification and neurotransmitter biosynthesis as mediators of the quantitative response to the drug. Subsequently, genetic analysis and metabolite profiling confirmed a complex role for ornithine and GABA levels in modification of survival time upon chronic nicotine exposure. Differences between populations from North Carolina and California suggest that the resistance mechanism may be an evolved response to environmental exposure.

Fine-scale crossover rate heterogeneity in Drosophila pseudoobscura

Broad-scale differences in crossover rate across the genome have been characterized in most genomes studied. Fine-scale differences, however, have only been examined in a few taxa, such as Arabidopsis, yeast, humans, and mice. No prior studies have directly looked for fine-scale recombination rate heterogeneity in Drosophila. We produced 370 Drosophila pseudoobscura containing a crossover event within the 2-megabase (MB) region between the genes yellow and white. We then examined 19 intervals within this region and determined where the crossovers occurred. We found that recombination events occur nonrandomly on a small scale and that mild "hotspots" of a few kilobases exist in Drosophila. Among the regions studied, recombination rates varied from 1.4 to 52 cM/MB. We also observed a trend toward high codon bias in regions of high recombination. Finally, we identified a significantly positive correlation between recombination rate and simple repeats, as well as the motif CACAC. These sequence features may contribute to broad-scale variation in crossover rate and, thus, shed light on features associated with crossover rate heterogeneity at a genome-wide scale.

Tests for the replication of an association between Egfr and natural variation in Drosophila melanogaster wing morphology

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.

Genetic variation for dorsal-ventral patterning of the Drosophila melanogaster eggshell

Patterning of the insect eggshell is an excellent system for exploring the molecular basis of phenotypic variation. In Drosophila melanogaster, two dorsal-anterior respiratory appendages are produced in response to signaling through the Epidermal growth factor receptor (Egfr). Previous work implicates Egfr pathway function in both intraspecific variation for dorsal appendage spacing (DAS) on the eggshell, as well as interspecific differences in dorsal appendage number and location. To test the hypothesis that genetic variation in Egfr contributes to variation in eggshell patterning, we have made use of naturally occurring intraspecific variation for DAS as a model quantitative trait. We found that there is substantial segregating genetic variation for DAS in D. melanogaster, and have tested for associations with 289 common polymorphisms in the Egfr locus. A marginal association was seen with two polymorphic sites in Egfr; however, we failed to replicate these findings in a second population, or in a modified quantitative complementation test designed to specifically test the effects of the putative polymorphisms. Therefore, we conclude that the polymorphisms we have identified in Egfr do not contribute to variation in DAS, and further work is required to understand the genetic architecture of this trait.

Replication of an Egfr-wing shape association in a wild-caught cohort of Drosophila melanogaster

Linkage disequilibrium mapping has been used extensively in medical and evolutionary genetics to map causal polymorphisms within genes associated with disease status or phenotypic variation for a trait. However, the initial findings of most nonhuman studies have not been replicated in subsequent studies, due in part to false positives, as well as additional factors that can render true positives unreplicable. These factors may be more severe when the initial study is performed using an experimental population of organisms reared under controlled lab conditions. We demonstrate that despite considerable phenotypic differences for wing shape between a lab-reared experimental population and a wild-caught cohort of Drosophila melanogaster, an association between a putative regulatory polymorphism in Egfr and wing shape can be replicated. These results are discussed both within the framework of future association-mapping studies and within the context of the evolutionary dynamics of alleles in populations.

Association between nucleotide variation in Egfr and wing shape in Drosophila melanogaster

As part of an effort to dissect quantitative trait locus effects to the nucleotide level, association was assessed between 238 single-nucleotide and 20 indel polymorphisms spread over 11 kb of the Drosophila melanogaster Egfr locus and nine relative warp measures of wing shape. One SNP in a conserved potential regulatory site for a GAGA factor in the promoter of alternate first exon 2 approaches conservative experiment-wise significance (P < 0.00003) in the sample of 207 lines for association with the location of the crossveins in the central region of the wing. Several other sites indicate marginal association with one or more other aspects of shape. No strong effects of sex or population of origin were detected with measures of shape, but two different sites were strongly associated with overall wing size in interaction with these fixed factors. Whole-gene sequencing in very large samples, rather than selective genotyping, would appear to be the only strategy likely to be successful for detecting subtle associations in species with high polymorphism and little haplotype structure. However, these features severely limit the ability of linkage disequilibrium mapping in Drosophila to resolve quantitative effects to single nucleotides.

Effects of population structure and sex on association between serotonin receptors and Drosophila heart rate

As a first step toward population and quantitative genetic analysis of neurotransmitter receptors in Drosophila melanogaster, we describe the parameters of nucleotide variation in three serotonin receptors and their association with pupal heart rate. Thirteen kilobases of DNA including the complete coding regions of 5-HT1A, 5-HT1B, and 5-HT2 were sequenced in 216 highly inbred lines extracted from two North American populations in California and North Carolina. Nucleotide and amino acid polymorphism is in the normal range for Drosophila genes and proteins, and linkage disequilibrium decays rapidly such that haplotype blocks are typically only a few SNPs long. However, intron 1 of 5-HT1A consists of two haplotypes that are at significantly different frequencies in the two populations. Neither this region of the gene nor any of the common amino acid polymorphisms in the three loci associate with either heart rate or heart rate variability. A cluster of SNPs in intron 2 of 5-HT1A, including a triallelic site, do show a highly significant interaction between genotype, sex, and population. While it is likely that a combination of weak, complex selection pressures and population structure has helped shape variation in the serotonin receptors of Drosophila, much larger sampling strategies than are currently adopted in evolutionary genetics will be required to disentangle these effects.