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

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.

Using inbreeding to test the contribution of non-additive genetic effects to additive genetic variance: a case study in Drosophila serrata

Additive genetic variance, V_A, is the key parameter for predicting adaptive and neutral phenotypic evolution. Changes in demography (e.g. increased close-relative inbreeding) can alter V_A, but how they do so depends on the (typically unknown) gene action and allele frequencies across many loci. For example, V_A increases proportionally with the inbreeding coefficient when allelic effects are additive, but smaller (or larger) increases can occur when allele frequencies are unequal at causal loci with dominance effects. Here, we describe an experimental approach to assess the potential for dominance effects to deflate V_A under inbreeding. Applying a powerful paired pedigree design in Drosophila serrata, we measured 11 wing traits on half-sibling families bred via either random or sibling mating, differing only in homozygosity (not allele frequency). Despite close inbreeding and substantial power to detect small V_A, we detected no deviation from the expected additive effect of inbreeding on genetic (co)variances. Our results suggest the average dominance coefficient is very small relative to the additive effect, or that allele frequencies are relatively equal at loci affecting wing traits. We outline the further opportunities for this paired pedigree approach to reveal the characteristics of V_A, providing insight into historical selection and future evolutionary potential.

Genomics of altitude-associated wing shape in two tropical butterflies

Understanding how organisms adapt to their local environment is central to evolution. With new whole-genome sequencing technologies and the explosion of data, deciphering the genomic basis of complex traits that are ecologically relevant is becoming increasingly feasible. Here, we studied the genomic basis of wing shape in two Neotropical butterflies that inhabit large geographical ranges. Heliconius butterflies at high elevations have been shown to generally have rounder wings than those in the lowlands. We reared over 1,100 butterflies from 71 broods of H. erato and H. melpomene in common-garden conditions and showed that wing aspect ratio, that is, elongatedness, is highly heritable in both species and that elevation-associated wing aspect ratio differences are maintained. Genome-wide associations with a published data set of 666 whole genomes from across a hybrid zone, uncovered a highly polygenic basis to wing aspect ratio variation in the wild. We identified several genes that have roles in wing morphogenesis or wing aspect ratio variation in Drosophila flies, making them promising candidates for future studies. There was little evidence for molecular parallelism in the two species, with only one shared candidate gene, nor for a role of the four known colour pattern loci, except for optix in H. erato. Thus, we present the first insights into the heritability and genomic basis of within-species wing aspect ratio in two Heliconius species, adding to a growing body of evidence that polygenic adaptation may underlie many ecologically relevant traits.

A single-nucleotide change underlies the genetic assimilation of a plastic trait

Genetic assimilation-the evolutionary process by which an environmentally induced phenotype is made constitutive-represents a fundamental concept in evolutionary biology. Thought to reflect adaptive phenotypic plasticity, matricidal hatching in nematodes is triggered by maternal nutrient deprivation to allow for protection or resource provisioning of offspring. Here, we report natural Caenorhabditis elegans populations harboring genetic variants expressing a derived state of near-constitutive matricidal hatching. These variants exhibit a single amino acid change (V530L) in KCNL-1, a small-conductance calcium-activated potassium channel subunit. This gain-of-function mutation causes matricidal hatching by strongly reducing the sensitivity to environmental stimuli triggering egg-laying. We show that reestablishing the canonical KCNL-1 protein in matricidal isolates is sufficient to restore canonical egg-laying. While highly deleterious in constant food environments, KCNL-1 V530L is maintained under fluctuating resource availability. A single point mutation can therefore underlie the genetic assimilation-by either genetic drift or selection-of an ancestrally plastic trait.

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.

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.

Natural Genetic Variation and Candidate Genes for Morphological Traits in Drosophila melanogaster

Body size is a complex character associated to several fitness related traits that vary within and between species as a consequence of environmental and genetic factors. Latitudinal and altitudinal clines for different morphological traits have been described in several species of Drosophila and previous work identified genomic regions associated with such variation in D. melanogaster. However, the genetic factors that orchestrate morphological variation have been barely studied. Here, our main objective was to investigate genetic variation for different morphological traits associated to the second chromosome in natural populations of D. melanogaster along latitudinal and altitudinal gradients in Argentina. Our results revealed weak clinal signals and a strong population effect on morphological variation. Moreover, most pairwise comparisons between populations were significant. Our study also showed important within-population genetic variation, which must be associated to the second chromosome, as the lines are otherwise genetically identical. Next, we examined the contribution of different candidate genes to natural variation for these traits. We performed quantitative complementation tests using a battery of lines bearing mutated alleles at candidate genes located in the second chromosome and six second chromosome substitution lines derived from natural populations which exhibited divergent phenotypes. Results of complementation tests revealed that natural variation at all candidate genes studied, invected, Fasciclin 3, toucan, Reticulon-like1, jing and CG14478, affects the studied characters, suggesting that they are Quantitative Trait Genes for morphological traits. Finally, the phenotypic patterns observed suggest that different alleles of each gene might contribute to natural variation for morphological traits. However, non-additive effects cannot be ruled out, as wild-derived strains differ at myriads of second chromosome loci that may interact epistatically with mutant alleles.

Looking Beyond the Genes: The Interplay Between Signaling Pathways and Mechanics in the Shaping and Diversification of Epithelial Tissues

The core of Evo-Devo lies in the intuition that the way tissues grow during embryonic development, the way they sustain their structure and function throughout lifetime, and the way they evolve are closely linked. Epithelial tissues are ubiquitous in metazoans, covering the gut and internal branched organs, as well as the skin and its derivatives (ie, teeth). Here, we discuss in vitro, in vivo, and in silico studies on epithelial tissues to illustrate the conserved, dynamical, and complex aspects of their development. We then explore the implications of the dynamical and nonlinear nature of development on the evolution of their size and shape at the phenotypic and genetic levels. In rare cases, when the interplay between signaling and mechanics is well understood at the cell level, it is becoming clear that the structure of development leads to covariation of characters, an integration which in turn provides some predictable structure to evolutionary changes. We suggest that such nonlinear systems are prone to genetic drift, cryptic genetic variation, and context-dependent mutational effects. We argue that experimental and theoretical studies at the cell level are critical to our understanding of the phenotypic and genetic evolution of epithelial tissues, including carcinomas.

Abundance and distribution of transposable elements in two Drosophila QTL mapping resources

Here we present computational machinery to efficiently and accurately identify transposable element (TE) insertions in 146 next-generation sequenced inbred strains of Drosophila melanogaster. The panel of lines we use in our study is composed of strains from a pair of genetic mapping resources: the Drosophila Genetic Reference Panel (DGRP) and the Drosophila Synthetic Population Resource (DSPR). We identified 23,087 TE insertions in these lines, of which 83.3% are found in only one line. There are marked differences in the distribution of elements over the genome, with TEs found at higher densities on the X chromosome, and in regions of low recombination. We also identified many more TEs per base pair of intronic sequence and fewer TEs per base pair of exonic sequence than expected if TEs are located at random locations in the euchromatic genome. There was substantial variation in TE load across genes. For example, the paralogs derailed and derailed-2 show a significant difference in the number of TE insertions, potentially reflecting differences in the selection acting on these loci. When considering TE families, we find a very weak effect of gene family size on TE insertions per gene, indicating that as gene family size increases the number of TE insertions in a given gene within that family also increases. TEs are known to be associated with certain phenotypes, and our data will allow investigators using the DGRP and DSPR to assess the functional role of TE insertions in complex trait variation more generally. Notably, because most TEs are very rare and often private to a single line, causative TEs resulting in phenotypic differences among individuals may typically fail to replicate across mapping panels since individual elements are unlikely to segregate in both panels. Our data suggest that “burden tests” that test for the effect of TEs as a class may be more fruitful.

Identifying permethrin resistance loci in malaria vectors by genetic mapping

Identification of the major loci responsible for insecticide resistance in malaria vectors would aid the development and implementation of effective resistance management strategies, which are urgently needed to tackle the growing threat posed by resistance to the limited insecticides available for malaria control. Genome-wide association studies in the major malaria vector, Anopheles gambiae, have been hindered by the high degree of within-population structuring and very low levels of linkage disequilibrium hence we revisited the use of quantitative trait loci (QTL) mapping to study resistance phenotypes in this vector species. Earlier work, identified two major QTL associated with pyrethroid resistance in A. gambiae s.s. from East Africa using genetic crossing of laboratory-colonized resistant and susceptible strains. In this study, we report the results from genetic mapping of pyrethroid resistance in three isofemale pedigrees established from wild-caught female A. gambiae s.s. mosquitoes from Benin. We identified two QTL on chromosomes 2L and 3R in these field populations, in similar genomic locations to the QTL identified in laboratory strains. The relative merits of two alternative study designs are discussed and suggestions made for future genetic mapping studies of insecticide resistance in mosquitoes.

Joint allelic effects on fitness and metric traits

Theoretical explanations of empirically observed standing genetic variation, mutation, and selection suggest that many alleles must jointly affect fitness and metric traits. However, there are few direct demonstrations of the nature and extent of these pleiotropic associations. We implemented a mutation accumulation (MA) divergence experimental design in Drosophila serrata to segregate genetic variants for fitness and metric traits. By exploiting naturally occurring MA line extinctions as a measure of line-level total fitness, manipulating sexual selection, and measuring productivity we were able to demonstrate genetic covariance between fitness and standard metric traits, wing size, and shape. Larger size was associated with lower total fitness and male sexual fitness, but higher productivity. Multivariate wing shape traits, capturing major axes of wing shape variation among MA lines, evolved only in the absence of sexual selection, and to the greatest extent in lines that went extinct, indicating that mutations contributing wing shape variation also typically had deleterious effects on both total fitness and male sexual fitness. This pleiotropic covariance of metric traits with fitness will drive their evolution, and generate the appearance of selection on the metric traits even in the absence of a direct contribution to fitness.

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.

ZIP8 zinc transporter: indispensable role for both multiple-organ organogenesis and hematopoiesis in utero

Previously this laboratory characterized Slc39a8-encoded ZIP8 as a Zn(2+)/(HCO(3)(-))(2) symporter; yet, the overall physiological importance of ZIP8 at the whole-organism level remains unclear. Herein we describe the phenotype of the hypomorphic Slc39a8(neo/neo) mouse which has retained the neomycin-resistance gene in intron 3, hence causing significantly decreased ZIP8 mRNA and protein levels in embryo, fetus, placenta, yolk sac, and several tissues of neonates. The Slc39a8(neo) allele is associated with diminished zinc and iron uptake in mouse fetal fibroblast and liver-derived cultures; consequently, Slc39a8(neo/neo) newborns exhibit diminished zinc and iron levels in several tissues. Slc39a8(neo/neo) homozygotes from gestational day(GD)-11.5 onward are pale, growth-stunted, and die between GD18.5 and 48 h postnatally. Defects include: severely hypoplastic spleen; hypoplasia of liver, kidney, lung, and lower limbs. Histologically, Slc39a8(neo/neo) neonates show decreased numbers of hematopoietic islands in yolk sac and liver. Low hemoglobin, hematocrit, red cell count, serum iron, and total iron-binding capacity confirmed severe anemia. Flow cytometry of fetal liver cells revealed the erythroid series strikingly affected in the hypomorph. Zinc-dependent 5-aminolevulinic acid dehydratase, required for heme synthesis, was not different between Slc39a8(+/+) and Slc39a8(neo/neo) offspring. To demonstrate further that the mouse phenotype is due to ZIP8 deficiency, we bred Slc39a8(+/neo) with BAC-transgenic BTZIP8-3 line (carrying three extra copies of the Slc39a8 allele); this cross generated viable Slc39a8(neo/neo)_BTZIP8-3(+/+) pups showing none of the above-mentioned congenital defects-proving Slc39a8(neo/neo) causes the described phenotype. Our study demonstrates that ZIP8-mediated zinc transport plays an unappreciated critical role during in utero and neonatal growth, organ morphogenesis, and hematopoiesis.

Developmental stability: a major role for cyclin G in drosophila melanogaster

Morphological consistency in metazoans is remarkable given the pervasive occurrence of genetic variation, environmental effects, and developmental noise. Developmental stability, the ability to reduce developmental noise, is a fundamental property of multicellular organisms, yet its genetic bases remains elusive. Imperfect bilateral symmetry, or fluctuating asymmetry, is commonly used to estimate developmental stability. We observed that Drosophila melanogaster overexpressing Cyclin G (CycG) exhibit wing asymmetry clearly detectable by sight. Quantification of wing size and shape using geometric morphometrics reveals that this asymmetry is a genuine-but extreme-fluctuating asymmetry. Overexpression of CycG indeed leads to a 40-fold increase of wing fluctuating asymmetry, which is an unprecedented effect, for any organ and in any animal model, either in wild populations or mutants. This asymmetry effect is not restricted to wings, since femur length is affected as well. Inactivating CycG by RNAi also induces fluctuating asymmetry but to a lesser extent. Investigating the cellular bases of the phenotypic effects of CycG deregulation, we found that misregulation of cell size is predominant in asymmetric flies. In particular, the tight negative correlation between cell size and cell number observed in wild-type flies is impaired when CycG is upregulated. Our results highlight the role of CycG in the control of developmental stability in D. melanogaster. Furthermore, they show that wing developmental stability is normally ensured via compensatory processes between cell growth and cell proliferation. We discuss the possible role of CycG as a hub in a genetic network that controls developmental stability.

Population genomics of wild and laboratory zebrafish (Danio rerio)

Understanding a wider range of genotype–phenotype associations can be achieved through ecological and evolutionary studies of traditional laboratory models. Here, we conducted the first large-scale geographic analysis of genetic variation within and among wild zebrafish (Danio rerio) populations occurring in Nepal, India, and Bangladesh, and we genetically compared wild populations to several commonly used lab strains. We examined genetic variation at 1832 polymorphic EST-based single nucleotide polymorphisms (SNPs) and the cytb mitochondrial gene in 13 wild populations and three lab strains. Natural populations were subdivided into three major mitochondrial DNA clades with an average among-clade sequence divergence of 5.8%. SNPs revealed five major evolutionarily and genetically distinct groups with an overall F_ST of 0.170 (95% CI 0.105–0.254). These genetic groups corresponded to discrete geographic regions and appear to reflect isolation in refugia during past climate cycles. We detected 71 significantly divergent outlier loci (3.4%) and nine loci (0.5%) with significantly low F_ST values. Valleys of reduced heterozygosity, consistent with selective sweeps, surrounded six of the 71 outliers (8.5%). The lab strains formed two additional groups that were genetically distinct from all wild populations. An additional subset of outlier loci was consistent with domestication selection within lab strains. Substantial genetic variation that exists in zebrafish as a whole is missing from lab strains that we analysed. A combination of laboratory and field studies that incorporates genetic variation from divergent wild populations along with the wealth of molecular information available for this model organism provides an opportunity to advance our understanding of genetic influences on phenotypic variation for a vertebrate species.

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.

Generation of a Slc39a8 hypomorph mouse: markedly decreased ZIP8 Zn²⁺/(HCO₃⁻)₂ transporter expression

Previously this laboratory has identified the mouse Slc39a8 gene encoding the ZIP8 transporter, important in cadmium uptake. ZIP8 functions endogenously as a electroneutral Zn(2+)/(HCO(3)(-))(2) symporter, moving both ions into the cell. The overall physiological importance of ZIP8 remains unclear. Herein we describe generation of a mouse line carrying the Slc39a8(neo) allele, containing the Frt-flanked neomycin-resistance (neo) mini-cassette in intron 3 and loxP sites in introns 3 and 6. Cre recombinase functions correctly in Escherichia coli and in adeno-Cre-infected mouse fetal fibroblasts, but does not function in the intact mouse for reasons not clear. Slc39a8(neo) is a hypomorphic allele, because Slc39a8(neo/neo) homozygotes exhibit dramatically decreased ZIP8 expression in embryo, fetus, and visceral yolk sac - in comparison to their littermate wild-type controls. This ZIP8 hypomorph will be instrumental in studying developmental and in utero physiological functions of the ZIP8 transporter.

Genome-wide deficiency mapping of the regions responsible for temporal canalization of the developmental processes of Drosophila melanogaster

Developmental processes of organisms are programmed to proceed in a finely regulated manner and finish within a certain period of time depending on the ambient environmental conditions. Therefore, variation in the developmental period under controlled genetic and environmental conditions indicates innate instability of the developmental process. In this study, we aimed to determine whether a molecular machinery exists that regulates the canalization of the developmental period and, if so, to test whether the same mechanism also stabilizes a morphological trait. To search for regions that influence the instability of the developmental period, we conducted genome-wide deficiency mapping with 441 isogenic deficiency strains covering 65.5% of the Drosophila melanogaster genome. We found that 11 independent deficiencies significantly increased the instability of the developmental period and 5 of these also significantly increased the fluctuating asymmetry of wing shape although there was no significant correlation between the instabilities of developmental period and wing shape in general. These results suggest that canalization processes of the developmental period and morphological traits are at least partially independent. Our findings emphasize the potential importance of temporal variation in development as an indicator of developmental stability and canalization and provide a novel perspective for understanding the regulation of phenotypic variability.

Getting real with real-time qPCR: a case study of reference gene selection for morphological variation in Drosophila melanogaster wings

Accurate estimation of gene expression differences during development requires sensitive techniques combined with gold-standard normalization procedures. This is particularly true in the case of quantitative traits, where expression changes might be small. Nevertheless, systematic selection and validation of reference genes has been overlooked, even in Drosophila studies. Here, we tested the stability of six traditional reference genes across samples of imaginal wing disks from morphologically divergent strains of Drosophila melanogaster, in a two-class comparison: quantitative or qualitative variation in wing morphology. Overall, we identified and validated a pair of genes (RpL32 and Tbp) as being stably expressed in both experimental comparisons. These genes might be considered as a bona fide pair of reference genes for gene expression analyses of morphological divergence in D. melanogaster wings. They might also be taken as good candidates for experimental identification of stable reference genes in other morphological comparisons using Drosophila or other insect species. Besides, we found that some genes traditionally used as reference in qPCR experiments were not stably expressed in wing disks from the different fly strains. In fact, a significant bias was observed when the expression of three genes of interest, which are involved in the regulation of growth and patterning during imaginal wing development, was normalized with such putative reference genes. Our results demonstrate how inaccurate findings and opposite conclusions might be drawn if traditional reference genes are arbitrarily used for internal normalization without proper validation in the given experimental condition, a practice still common in qPCR experiments.

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.

Association mapping in outbred populations: power and efficiency when genotyping parents and phenotyping progeny

We develop expressions for the power to detect associations between parental genotypes and offspring phenotypes for quantitative traits. Three different "indirect" experimental designs are considered: full-sib, half-sib, and full-sib-half-sib families. We compare the power of these designs to detect genotype-phenotype associations relative to the common, "direct," approach of genotyping and phenotyping the same individuals. When heritability is low, the indirect designs can outperform the direct method. However, the extra power comes at a cost due to an increased phenotyping effort. By developing expressions for optimal experimental designs given the cost of phenotyping relative to genotyping, we show how the extra costs associated with phenotyping a large number of individuals will influence experimental design decisions. Our results suggest that indirect association studies can be a powerful means of detecting allelic associations in outbred populations of species for which genotyping and phenotyping the same individuals is impractical and for life history and behavioral traits that are heavily influenced by environmental variance and therefore best measured on groups of individuals. Indirect association studies are likely to be favored only on purely economical grounds, however, when phenotyping is substantially less expensive than genotyping. A web-based application implementing our expressions has been developed to aid in the design of indirect association studies.

Intraspecific evolution of the intercellular signaling network underlying a robust developmental system

Many biological systems produce an invariant output when faced with stochastic or environmental variation. This robustness of system output to variation affecting the underlying process may allow for "cryptic" genetic evolution within the system without change in output. We studied variation of cell fate patterning of Caenorhabditis elegans vulva precursors, a developmental system that relies on a simple intercellular signaling network and yields an invariant output of cell fates and lineages among C. elegans wild isolates. We first investigated the system's genetic variation in C. elegans by means of genetic tools and cell ablation to break down its buffering mechanisms. We uncovered distinct architectures of quantitative variation along the Ras signaling cascade, including compensatory variation, and differences in cell sensitivity to induction along the anteroposterior axis. In the unperturbed system, we further found variation between isolates in spatio-temporal dynamics of Ras pathway activity, which can explain the phenotypic differences revealed upon perturbation. Finally, the variation mostly affects the signaling pathways in a tissue-specific manner. We thus demonstrate and characterize microevolution of a developmental signaling network. In addition, our results suggest that the vulva genetic screens would have yielded a different mutation spectrum, especially for Wnt pathway mutations, had they been performed in another C. elegans genetic background.

The genetic architecture of complex traits in teosinte (Zea mays ssp. parviglumis): new evidence from association mapping

Previous association analyses showed that variation at major regulatory genes contributes to standing variation for complex traits in Balsas teosinte, the progenitor of maize. This study expands our previous association mapping effort in teosinte by testing 123 markers in 52 candidate genes for association with 31 traits in a population of 817 individuals. Thirty-three significant associations for markers from 15 candidate genes and 10 traits survive correction for multiple testing. Our analyses suggest several new putative causative relationships between specific genes and trait variation in teosinte. For example, two ramosa genes (ra1 and ra2) associate with ear structure, and the MADS-box gene, zagl1, associates with ear shattering. Since zagl1 was previously shown to be a target of selection during maize domestication, we suggest that this gene was under selection for its effect on the loss of ear shattering, a key domestication trait. All observed effects were relatively small in terms of the percentage of phenotypic variation explained (<10%). We also detected several epistatic interactions between markers in the same gene that associate with the same trait. Candidate-gene-based association mapping appears to be a promising method for investigating the inheritance of complex traits in teosinte.

Associations between sperm competition and natural variation in male reproductive genes on the third chromosome of Drosophila melanogaster

We applied association analysis to elucidate the genetic basis for variation in phenotypes affecting postcopulatory sexual selection in a natural population of Drosophila melanogaster. We scored 96 third chromosome substitution lines for nine phenotypes affecting sperm competitive ability and genotyped them at 72 polymorphisms in 13 male reproductive genes. Significant heterogeneity among lines (P < 0.01) was detected for all phenotypes except male-induced refractoriness (P = 0.053). We identified 24 associations (8 single-marker associations, 12 three-marker haplotype associations, and 4 cases of epistasis revealed by single-marker interactions). Fewer than 9 of these associations are likely to be false positives. Several associations were consistent with previous findings [Acp70A with the male's influence on the female's refractoriness to remating (refractory), Esterase-6 with a male's remating probability (remating) and a measure of female offspring production (fecundity)], but many are novel associations with uncharacterized seminal fluid proteins. Four genes showed evidence for pleiotropic effects [CG6168 with a measure of sperm competition (P2') and refractory, CG14560 with a defensive measure of sperm competition (P1') and a measure of female fecundity, Acp62F with P2' and a measure of female fecundity, and Esterase-6 with remating and a measure of female fecundity]. Our findings provide evidence that pleiotropy and epistasis are important factors in the genetic architecture of male reproductive success and show that haplotype analyses can identify associations missed in the single-marker approach.

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.

Dissecting the complex genetic basis of mate choice

The genetic analysis of mate choice is fraught with difficulties. Males produce complex signals and displays that can consist of a combination of acoustic, visual, chemical and behavioural phenotypes. Furthermore, female preferences for these male traits are notoriously difficult to quantify. During mate choice, genes not only affect the phenotypes of the individual they are in, but can influence the expression of traits in other individuals. How can genetic analyses be conducted to encompass this complexity? Tighter integration of classical quantitative genetic approaches with modern genomic technologies promises to advance our understanding of the complex genetic basis of mate choice.

Epidermal growth factor receptor and transforming growth factor-beta signaling contributes to variation for wing shape 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.

Natural Malaria Infection in Anopheles gambiae Is Regulated by a Single Genomic Control Region

We surveyed an Anopheles gambiae population in a West African malaria transmission zone for naturally occurring genetic loci that control mosquito infection with the human malaria parasite, Plasmodium falciparum. The strongest Plasmodium resistance loci cluster in a small region of chromosome 2L and each locus explains at least 89% of parasite-free mosquitoes in independent pedigrees. Together, the clustered loci form a genomic Plasmodium-resistance island that explains most of the genetic variation for malaria parasite infection of mosquitoes in nature. Among the candidate genes in this chromosome region, RNA interference knockdown assays confirm a role in Plasmodium resistance for Anopheles Plasmodium-responsive leucine-rich repeat 1 (APL1), encoding a leucine-rich repeat protein that is similar to molecules involved in natural pathogen resistance mechanisms in plants and mammals.

The effect of polymorphisms in the enhancer of split gene complex on bristle number variation in a large wild-caught cohort of Drosophila melanogaster

The Enhancer of split complex [E(spl)-C] in Drosophila encompasses a variety of functional elements controlling bristle patterning and on the basis of prior work is a strong candidate for harboring alleles having subtle effects on bristle number variation. Here we extend earlier studies identifying associations between complex phenotypes and polymorphisms segregating among inbred laboratory lines of Drosophila and test the influence of E(spl)-C on bristle number variation in a natural cohort. We describe results from an association mapping study using 203 polymorphisms spread throughout the E(spl)-C genotyped in 2000 wild-caught Drosophila melanogaster. Despite power to detect associations accounting for as little as 2% of segregating variation for bristle number, and saturating the region with single-nucleotide polymorphisms (SNPs), we identified no single SNP marker showing a significant (additive over loci) effect after correcting for multiple tests. Using a newly developed test we conservatively identify six regions of the E(spl)-C in which the insertion of transposable elements as a class contributes to variation in bristle number, apparently in a sex- or trait-limited fashion. Finally, we carry out all possible 20,503 two-way tests for epistasis and identify a slight excess of marginally significant interactions, although none survive multiple-testing correction. It may not be straightforward to extend the results of laboratory-based association studies to natural populations.

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.