DNA methylation reprogramming in teleosts.

Early embryonic development is crucially important but also remarkably diverse among animal taxa. Axis formation and cell lineage specification occur due to both spatial and temporal control of gene expression. This complex system involves various signaling pathways and developmental genes such as transcription factors as well as other molecular interactants that maintain cellular states, including several types of epigenetic marks. 5mC DNA methylation, the chemical modification of cytosines in eukaryotes, represents one such mark. By influencing the compaction of chromatin (a high-order DNA structure), DNA methylation can either repress or induce transcriptional activity. Mammals exhibit a reprogramming of DNA methylation from the parental genomes in the zygote following fertilization, and later in primordial germ cells (PGCs). Whether these periods of methylation reprogramming are evolutionarily conserved, or an innovation in mammals, is an emerging question. Looking into these processes in other vertebrate lineages is thus important, and teleost fish, with their extensive species richness, phenotypic diversity, and multiple rounds of whole genome duplication, provide the perfect research playground for answering such a question. This review aims to present a concise state of the art of DNA methylation reprogramming in early development in fish by summarizing findings from different research groups investigating methylation reprogramming patterns in teleosts, while keeping in mind the ramifications of the methodology used, then comparing those patterns to reprogramming patterns in mammals.

The prey availability and diet of juvenile Atlantic salmon (Salmo salar L.) in low-productivity rivers in northern Europe.

The availability of resources varies across a species distributional range, and a low-productivity area can make a species more vulnerable. We investigated the invertebrate composition and prey choice of juvenile Atlantic salmon (Salmo salar L.) in low-productivity rivers in northeast Iceland, which is one of the species' most northerly distributions. By sampling benthic and drift invertebrate populations, we found that prey availability was similar within and between rivers. Gut content samples showed that the main prey choice for juvenile S. salar was the Chironomidae. The type of food items consumed varied across different weight groups of S. salar, with smaller juveniles having more diverse diet. S. salar did not have a selection preference for chironomids, which indicates that they were eating the highly available prey in their environment, rather than hunting high biomass items such as terrestrial invertebrates and large Dipterans. Estimates of dietary niche showed that S. salar in these low-productivity rivers relied on consuming what was most readily available, the chironomids, and that they must share resources with other salmonid species. This may be due to the low diversity of freshwater invertebrates (fewer prey options), whereas S. salar in nutrient-rich rivers could rely more on terrestrial invertebrates as an additional subsidy in their diet. In conclusion, with limited prey choices, juvenile S. salar in nutrient-poor rivers, especially in a biogeographically isolated region with low species diversity, may increase in vulnerability and decrease in adaptability to environmental change. Management methods that increase benthic prey abundance and diversity are recommended for conserving the S. salar population in a nutrient-poor river.

Variation in personality shaped by evolutionary history, genotype and developmental plasticity in response to feeding modalities in the Arctic charr.

Animal personality has been shown to be influenced by both genetic and environmental factors and shaped by natural selection. Currently, little is known about mechanisms influencing the development of personality traits. This study examines the extent to which personality development is genetically influenced and/or environmentally responsive (plastic). We also investigated the role of evolutionary history, assessing whether personality traits could be canalized along a genetic and ecological divergence gradient. We tested the plastic potential of boldness in juveniles of five Icelandic Arctic charr morphs (Salvelinus alpinus), including two pairs of sympatric morphs, displaying various degrees of genetic and ecological divergence from the ancestral anadromous charr, split between treatments mimicking benthic versus pelagic feeding modalities. We show that differences in mean boldness are mostly affected by genetics. While the benthic treatment led to bolder individuals overall, the environmental effect was rather weak, suggesting that boldness lies under strong genetic influence with reduced plastic potential. Finally, we found hints of differences by morphs in boldness canalization through reduced variance and plasticity, and higher consistency in boldness within morphs. These findings provide new insights on how behavioural development may impact adaptive diversification.

Adaptive cellular evolution or cellular system drift in hares.

Adaptations occur at many levels, for example, from DNA sequence of regulatory elements and cellular homeostatic systems to organismal physiology and behaviour (Mayr, 1997). Established adaptations are maintained by purifying and stabilizing selection. Students of animal diversity tend to focus on higher order traits, anatomy, physiology, organismal function and interactions. The core cellular and metabolic systems of metazoans evolved early in their history and are assumed to be rather similar between groups. The housekeeping functions and core metabolic functions of cells are generally considered relatively static, especially among closely related species. The extent to which evolution shapes core cellular metabolism and physiology in animals is largely unexplored. Ecological opportunities or strong positive selection can alter basal metabolic rate, activity levels and life-history traits (e.g., life span, age of maturity, offspring number) and potentially lead to divergence in core cellular and metabolic trait systems (Norin & Metcalfe, 2019; Speakman, 2005). Furthermore, systems under stabilizing selection can also change. Developmental systems of related species may produce the same phenotype or structure, but experience drift that can alter connections and even lead to turnover of cogs in the system (True & Haag, 2001). Are the cellular functions of animals highly constrained, subject to cellular system drift or affected by positive selection? This was tackled by a new study by Kateryna Gaertner and colleagues in a From the Cover manuscript in this issue of Molecular Ecology (Gaertner et al., 2022), using fibroblasts from the closely related but ecologically distinct brown and mountain hares.

DNA methylation differences during development distinguish sympatric morphs of Arctic charr (Salvelinus alpinus).

Changes in DNA methylation in specific coding or non-coding regions can influence development and potentially divergence in traits within species and groups. While the impact of epigenetic variation on developmental pathways associated with evolutionary divergence is the focus of intense investigation, few studies have looked at recently diverged systems. Phenotypic diversity between closely related populations of Arctic charr (Salvelinus alpinus), which diverged within the last 10,000 years, offers an interesting ecological model to address such effects. Using bisulphite sequencing, we studied general DNA methylation patterns during development in the four sympatric morphs of Arctic charr from Lake Thingvallavatn. The data revealed strong differences between developmental timepoints and between morphs (mainly along the benthic-limnetic axis), both at single CpG sites and in 1000 bp-regions. Genes located close to differentially methylated CpG sites were involved in nucleosome assembly, regulation of osteoclast differentiation, and cell-matrix adhesion. Differentially methylated regions were enriched in tRNA and rRNA sequences, and half of them were located close to transcription start sites. The expression of 14 genes showing methylation differences over time or between morphs was further investigated by qPCR and nine of these were found to be differentially expressed between morphs. Four genes (ARHGEF37-like, H3-like, MPP3 and MEGF9) showed a correlation between methylation and expression. Lastly, histone gene clusters displayed interesting methylation differences between timepoints and morphs, as well as intragenic methylation variation. The results presented here provide a motivation for further studies on the contribution of epigenetic traits, such as DNA methylation, to phenotypic diversity and developmental mechanisms.

Mitf is a master regulator of the v-ATPase, forming a control module for cellular homeostasis with v-ATPase and TORC1.

The v-ATPase is a fundamental eukaryotic enzyme that is central to cellular homeostasis. Although its impact on key metabolic regulators such as TORC1 is well documented, our knowledge of mechanisms that regulate v-ATPase activity is limited. Here, we report that the Drosophila transcription factor Mitf is a master regulator of this holoenzyme. Mitf directly controls transcription of all 15 v-ATPase components through M-box cis-sites and this coordinated regulation affects holoenzyme activity in vivo. In addition, through the v-ATPase, Mitf promotes the activity of TORC1, which in turn negatively regulates Mitf. We provide evidence that Mitf, v-ATPase and TORC1 form a negative regulatory loop that maintains each of these important metabolic regulators in relative balance. Interestingly, direct regulation of v-ATPase genes by human MITF also occurs in cells of the melanocytic lineage, showing mechanistic conservation in the regulation of the v-ATPase by MITF family proteins in fly and mammals. Collectively, this evidence points to an ancient module comprising Mitf, v-ATPase and TORC1 that serves as a dynamic modulator of metabolism for cellular homeostasis.

Bones in motion: Ontogeny of craniofacial development in sympatric arctic charr morphs.

BACKGROUND: The impressive diversity in the feeding apparatus often seen among related fish species clearly reflects differences in feeding modes and habitat utilization. Such variation can also be found within species. One example of such intraspecific diversity is the Arctic charr in Lake Thingvallavatn, where four distinct morphs coexist: two limnetic, with evenly protruding jaws, and two benthic, with subterminal lower jaws. We used these recently evolved morphs to study the role of ontogenetic variation in shaping craniofacial diversity. RESULTS: The segmental development of the pharyngeal arches and the order of events in craniofacial development is the same as has been described for teleosts, emphasizing the conserved nature of this process. However, our morphometric analyses reveal differences between morphs. Hatching is accompanied by increase in size and allometric shape changes in Arctic charr. Ontogenetic trajectories of craniofacial shape also differ significantly between morphs. CONCLUSIONS: The results point to hatching as a significant developmental event in Arctic charr and possibly other fishes. Also, the developmental origins of limnetic and benthic specializations in the craniofacial elements of Arctic charr may stem from events around hatching. This calls for investigations of the mechanisms and consequences of hatching, in the context of development and evolution. Developmental Dynamics 244:1168-1178, 2015. © 2015 Wiley Periodicals, Inc.

Diversity in the internal functional feeding elements of sympatric morphs of Arctic charr (Salvelinus alpinus).

The diversity of functional feeding anatomy is particularly impressive in fishes and correlates with various interspecific ecological specializations. Intraspecific polymorphism can manifest in divergent feeding morphology and ecology, often along a benthic-pelagic axis. Arctic charr (Salvelinus alpinus) is a freshwater salmonid known for morphological variation and sympatric polymorphism and in Lake Þingvallavatn, Iceland, four morphs of charr coexist that differ in preferred prey, behaviour, habitat use, and external feeding morphology. We studied variation in six upper and lower jaw bones in adults of these four morphs using geometric morphometrics and univariate statistics. We tested for allometric differences in bone size and shape among morphs, morph effects on bone size and shape, and divergence along the benthic-pelagic axis. We also examined the degree of integration between bone pairs. We found differences in bone size between pelagic and benthic morphs for two bones (dentary and premaxilla). There was clear bone shape divergence along a benthic-pelagic axis in four bones (dentary, articular-angular, premaxilla and maxilla), as well as allometric shape differences between morphs in the dentary. Notably for the dentary, morph explained more shape variation than bone size. Comparatively, benthic morphs possess a compact and taller dentary, with shorter dentary palate, consistent with visible (but less prominent) differences in external morphology. As these morphs emerged in the last 10,000 years, these results indicate rapid functional evolution of specific feeding structures in arctic charr. This sets the stage for studies of the genetics and development of rapid and parallel craniofacial evolution.

Variants conferring risk of atrial fibrillation on chromosome 4q25.

Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia in humans and is characterized by chaotic electrical activity of the atria. It affects one in ten individuals over the age of 80 years, causes significant morbidity and is an independent predictor of mortality. Recent studies have provided evidence of a genetic contribution to AF. Mutations in potassium-channel genes have been associated with familial AF but account for only a small fraction of all cases of AF. We have performed a genome-wide association scan, followed by replication studies in three populations of European descent and a Chinese population from Hong Kong and find a strong association between two sequence variants on chromosome 4q25 and AF. Here we show that about 35% of individuals of European descent have at least one of the variants and that the risk of AF increases by 1.72 and 1.39 per copy. The association with the stronger variant is replicated in the Chinese population, where it is carried by 75% of individuals and the risk of AF is increased by 1.42 per copy. A stronger association was observed in individuals with typical atrial flutter. Both variants are adjacent to PITX2, which is known to have a critical function in left-right asymmetry of the heart.

Genetic structure and relatedness of brown trout (Salmo trutta) populations in the drainage basin of the Ölfusá river, South-Western Iceland.

Background: Lake Þingvallavatn in Iceland, a part of the river Ölfusá drainage basin, was presumably populated by brown trout soon after it formed at the end of the last Ice Age. The genetic relatedness of the brown trout in Þingvallavatn to other populations in the Ölfusá drainage basin is unknown. After the building of a dam at the outlet of the lake in 1959 brown trout catches declined, though numbers have now increased. The aim of this study was to assess effects of geographic isolation and potential downstream gene flow on the genetic structure and diversity in brown trout sampled in several locations in the western side of the watershed of River Ölfusá. We hypothesized that brown trout in Lake Þingvallavatn constituted several local spawning populations connected by occasional gene flow before the damming of the lake. We also estimated the effective population size (NE) of some of these populations and tested for signs of a recent population bottleneck in Lake Þingvallavatn. Methods: We sampled brown trout inhabiting four lakes and 12 rivers within and near the watershed of River Ölfusá by means of electro- and net- fishing. After stringent data filtering, 2,597 polymorphic loci obtained from ddRADseq data from 317 individuals were ascertained as putative neutral markers. Results: Overall, the genetic relatedness of brown trout in the Ölfusá watershed reflected the connectivity and topography of the waterways. Ancestry proportion analyses and a phylogenetic tree revealed seven distinct clusters, some of which corresponded to small populations with reduced genetic diversity. There was no evidence of downstream gene flow from Lake Þingvallavatn, although gene flow was observed from much smaller mountain populations. Most locations showed low NE values (i.e., ~14.6 on average) while the putative anadromous trout from River Sog and the spawning population from River Öxará, that flows into Lake Þingvallavatn, showed notably higher NE values (i.e., 71.2 and 56.5, respectively). No signals of recent population bottlenecks were detected in the brown trout of Lake Þingvallavatn. Discussion: This is the first time that the genetic structure and diversity of brown trout in the watershed of River Ölfusá have been assessed. Our results point towards the presence of a metapopulation in the watershed of Lake Þingvallavatn, which has been influenced by restoration efforts and is now dominated by a genetic component originated in River Öxará. Many of the locations studied represent different populations. Those that are isolated in headwater streams and lakes are genetically distinct presenting low genetic diversity, yet they can be important in increasing the genetic variation in downstream populations. These populations should be considered for conservation and direct management.

DNA methylation programming and reprogramming in primate embryonic stem cells.

DNA methylation is an important epigenetic mechanism, affecting normal development and playing a key role in reprogramming epigenomes during stem cell derivation. Here we report on DNA methylation patterns in native monkey embryonic stem cells (ESCs), fibroblasts, and ESCs generated through somatic cell nuclear transfer (SCNT), identifying and comparing epigenome programming and reprogramming. We characterize hundreds of regions that are hyper- or hypomethylated in fibroblasts compared to native ESCs and show that these are conserved in human cells and tissues. Remarkably, the vast majority of these regions are reprogrammed in SCNT ESCs, leading to almost perfect correlation between the epigenomic profiles of the native and reprogrammed lines. At least 58% of these changes are correlated in cis to transcription changes, Polycomb Repressive Complex-2 occupancy, or binding by the CTCF insulator. We also show that while epigenomic reprogramming is extensive and globally accurate, the efficiency of adding and stripping DNA methylation during reprogramming is regionally variable. In several cases, this variability results in regions that remain methylated in a fibroblast-like pattern even after reprogramming.

Canalization of segmentation and its evolution in Drosophila.

Segmentation in Drosophila embryogenesis occurs through a hierarchical cascade of regulatory gene expression driven by the establishment of a diffusion-mediated morphogen gradient. Here, we investigate the response of this pattern formation process to genetic variation and evolution in egg size. Specifically, we ask whether spatial localization of gap genes Kruppel (Kr) and giant (gt) and the pair-rule gene even-skipped (eve) during cellularization is robust to genetic variation in embryo length in three Drosophila melanogaster isolines and two closely related species. We identified two wild-derived strains of D. melanogaster whose eggs differ by approximately 25% in length when reared under identical conditions. These two lines, a D. melanogaster laboratory stock (w1118), and offspring from crosses between the lines all exhibit precise scaling in the placement of gap and pair-rule gene expression along the anterior-posterior axis in relation to embryo length. Genetic analysis indicates that this scaling is maternally controlled. Maternal regulation of scaling must be required for consistent localization of segmentation gene expression because embryo size, a genetically variable and adaptive trait, is maternally inherited. We also investigated spatial scaling between these D. melanogaster lines and single lines of Drosophila sechellia and Drosophila simulans, the latter two differing by approximately 25% in egg length. In contrast to the robust scaling we observed within species, localization of gene expression relative to embryo length differs significantly between the three species. Thus, the developmental mechanism that assures robust scaling within a species does not prevent rapid evolution between species.

Extensive genetic differentiation between recently evolved sympatric Arctic charr morphs.

The availability of diverse ecological niches can promote adaptation of trophic specializations and related traits, as has been repeatedly observed in evolutionary radiations of freshwater fish. The role of genetics, environment, and history in ecologically driven divergence and adaptation, can be studied on adaptive radiations or populations showing ecological polymorphism. Salmonids, especially the Salvelinus genus, are renowned for both phenotypic diversity and polymorphism. Arctic charr (Salvelinus alpinus) invaded Icelandic streams during the glacial retreat (about 10,000 years ago) and exhibits many instances of sympatric polymorphism. Particularly, well studied are the four morphs in Lake Þingvallavatn in Iceland. The small benthic (SB), large benthic (LB), planktivorous (PL), and piscivorous (PI) charr differ in many regards, including size, form, and life history traits. To investigate relatedness and genomic differentiation between morphs, we identified variable sites from RNA-sequencing data from three of those morphs and verified 22 variants in population samples. The data reveal genetic differences between the morphs, with the two benthic morphs being more similar and the PL-charr more genetically different. The markers with high differentiation map to all linkage groups, suggesting ancient and pervasive genetic separation of these three morphs. Furthermore, GO analyses suggest differences in collagen metabolism, odontogenesis, and sensory systems between PL-charr and the benthic morphs. Genotyping in population samples from all four morphs confirms the genetic separation and indicates that the PI-charr are less genetically distinct than the other three morphs. The genetic separation of the other three morphs indicates certain degree of reproductive isolation. The extent of gene flow between the morphs and the nature of reproductive barriers between them remain to be elucidated.

Functional evolution of a cis-regulatory module.

Lack of knowledge about how regulatory regions evolve in relation to their structure-function may limit the utility of comparative sequence analysis in deciphering cis-regulatory sequences. To address this we applied reverse genetics to carry out a functional genetic complementation analysis of a eukaryotic cis-regulatory module-the even-skipped stripe 2 enhancer-from four Drosophila species. The evolution of this enhancer is non-clock-like, with important functional differences between closely related species and functional convergence between distantly related species. Functional divergence is attributable to differences in activation levels rather than spatiotemporal control of gene expression. Our findings have implications for understanding enhancer structure-function, mechanisms of speciation and computational identification of regulatory modules.

Differential gene expression during early development in recently evolved and sympatric Arctic charr morphs.

Phenotypic differences between closely related taxa or populations can arise through genetic variation or be environmentally induced, leading to altered transcription of genes during development. Comparative developmental studies of closely related species or variable populations within species can help to elucidate the molecular mechanisms related to evolutionary divergence and speciation. Studies of Arctic charr (Salvelinus alpinus) and related salmonids have revealed considerable phenotypic variation among populations and in Arctic charr many cases of extensive variation within lakes (resource polymorphism) have been recorded. One example is the four Arctic charr morphs in the ∼10,000 year old Lake Thingvallavatn, which differ in numerous morphological and life history traits. We set out to investigate the molecular and developmental roots of this polymorphism by studying gene expression in embryos of three of the morphs reared in a common garden set-up. We performed RNA-sequencing, de-novo transcriptome assembly and compared gene expression among morphs during an important timeframe in early development, i.e., preceding the formation of key trophic structures. Expectedly, developmental time was the predominant explanatory variable. As the data were affected by some form of RNA-degradation even though all samples passed quality control testing, an estimate of 3'-bias was the second most common explanatory variable. Importantly, morph, both as an independent variable and as interaction with developmental time, affected the expression of numerous transcripts. Transcripts with morph effect, separated the three morphs at the expression level, with the two benthic morphs being more similar. However, Gene Ontology analyses did not reveal clear functional enrichment of transcripts between groups. Verification via qPCR confirmed differential expression of several genes between the morphs, including regulatory genes such as AT-Rich Interaction Domain 4A (arid4a) and translin (tsn). The data are consistent with a scenario where genetic divergence has contributed to differential expression of multiple genes and systems during early development of these sympatric Arctic charr morphs.

Evidence that Egfr contributes to cryptic genetic variation for photoreceptor determination in natural populations of Drosophila melanogaster.

One objective of quantitative genetics is to identify the nucleotide variants within genes that contribute to phenotypic variation and susceptibility [1]. In an evolutionary context, this means characterizing the molecular polymorphisms that modify the penetrance and expressivity of perturbed traits. A survey of association between 267 SNPs in almost 11 kb of the D. melanogaster Egfr and the degree of eye roughening due to a gain-of-function Egfr(E1) allele crossed into 210 isogenic wild-type lines provides evidence that a handful of synonymous substitutions supply cryptic variation for photoreceptor determination. Ten sites exceed Bonferroni threshold for association in two sets of crosses to different Egfr(E1) backgrounds including a particularly significant cluster of sites in tight linkage disequilibrium toward the 3' end of the coding region. Epistatic interaction of this cluster with one other site enhances the expressivity of this haplotype. Replication of the strongest associations with an independent sample of 302 phenotypically extreme individuals derived from 1000 crosses of Egfr(E1) to freshly trapped males was achieved using modified case-control and transmission-disequilibrium tests. A tendency for the rarer alleles to have more disrupted eye development suggests that mutation-selection balance is a possible mechanism contributing to maintaining cryptic variation for Egfr.

Maedi-visna virus persistence: Antigenic variation and latency.

Maedi-visna virus (MVV), a lentivirus of sheep, shares with other lentiviruses the ability to establish a lifelong infection. In this study five sheep were infected intravenously with MVV and housed together with a number of uninfected sheep for natural transmission. All virus isolates from ten sheep that had been infected naturally had multiple mutations in the principal neutralization domain in Env and were antigenic variants, while three of four isolates from the carrier sheep had identical sequences to the infecting strain and were not antigenic variants. There was evidence of positive selection in the gene, particularly in amino acids comprising the neutralization epitope and some adjacent glycosylation sites. Together these results suggest that virus persistence is acquired by a reservoir of latent viruses, and that there is selection for antigenic variants of virus that is transmitted naturally.

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.

Fungal and cyanobacterial gene expression in a lichen symbiosis: Effect of temperature and location.

Organisms have evolved different cellular mechanisms to deal with environmental stress, primarily through complex molecular mechanisms including protein refolding and DNA repair. As mutualistic symbioses, lichens offer the possibility of analyzing molecular stress responses in a particularly tight interspecific relationship. We study the widespread cyanolichen Peltigera membranacea, a key player in carbon and nitrogen cycling in terrestrial ecosystems at northern latitudes. We ask whether increased temperature is reflected in mRNA levels of selected damage control genes, and do the response patterns show geographical associations? Using real-time PCR quantification of 38 transcripts, differential expression was demonstrated for nine cyanobacterial and nine fungal stress response genes (plus the fungal symbiosis-related lec2 gene) when the temperature was increased from 5 °C to 15 °C and 25 °C. Principle component analysis (PCA) revealed two gene groups with different response patterns. Whereas a set of cyanobacterial DNA repair genes and the fungal lec2 (PC1 group) showed an expression drop at 15 °C vs. 5 °C, most fungal candidates (PC2 group) showed increased expression at 25 °C vs. 5 °C. PC1 responses also correlated with elevation. The correlated downregulation of lec2 and cyanobacterial DNA repair genes suggests a possible interplay between the symbionts warranting further studies.

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.