QTL mapping reveals a striking coincidence in the positions of genomic regions associated with adaptive variation in body size in parallel clines of Drosophila melanogaster on different continents

Latitudinal genetic clines in body size are common in many ectotherm species and are attributed to climatic adaptation. Here, we use Quantitative Trait Loci (QTL) mapping to identify genomic regions associated with adaptive variation in body size in natural populations of Drosophila melanogaster from extreme ends of a cline in South America. Our results show that there is a significant association between the positions of QTL with strong effects on wing area in South America and those previously reported in a QTL mapping study of Australian cline end populations (P < 0.05). In both continents, the right arm of the third chromosome is associated with QTL with the strongest effect on wing area. We also show that QTL peaks for wing area and thorax length are associated with the same genomic regions, indicating that the clinal variation in the body size traits may have a similar genetic basis. The consistency of the results found for the South American and Australian cline end populations indicate that the genetic basis of the two clines may be similar and future efforts to identify the genes producing the response to selection should be focused on the genomic regions highlighted by the present work.     

The effect of developmental temperature on the genetic architecture underlying size and thermal clines in Drosophila melanogaster and D. simulans from the east coast of Australia

Body size and thermal tolerance clines in Drosophila melanogaster occur along the east coast of Australia. However the extent to which temperature affects the genetic architecture underlying the observed clinal divergence remains unknown. Clinal variation in these traits is associated with cosmopolitan chromosome inversions that cline in D. melanogaster. Whether this association influences the genetic architecture for these traits in D. melanogaster is unclear. Drosophila simulans shows linear clines in body size, but nonlinear clines in cold resistance. Clinally varying inversions are absent in D. simulans. Line-cross and clinal analyses were performed between tropical and temperate populations of D. melanogaster and D. simulans from the east coast of Australia to investigate whether clinal patterns and genetic effects contributing to clinal divergence in wing centroid size, thorax length, wing-to-thorax ratio, cold and heat resistance differed under different developmental temperatures (18 °C, 25 °C, and 29 °C). Developmental temperature influenced the genetic architecture in both species. Similarities between D. melanogaster and D. simulans suggest clinally varying inversion polymorphisms have little influence on the genetic architecture underlying clinal divergence in size in D. melanogaster. Differing genetic architectures across different temperatures highlight the need to consider different environments in future evolutionary and molecular studies of phenotypic divergence.     

An interspecific QTL study of Drosophila wing size and shape variation to investigate the genetic basis of morphological differences

The Drosophila wing has been used as a model in studies of morphogenesis and evolution; the use of such models can contribute to our understanding of mechanisms that promote morphological divergence among populations and species. We mapped quantitative trait loci (QTL) affecting wing size and shape traits using highly inbred introgression lines between D. simulans and D. sechellia, two sibling species of the melanogaster subgroup. Eighteen QTL peaks that are associated with 12 wing traits were identified, including two principal components. The wings of D. simulans and D. sechellia significantly diverged in size; two of the QTL peaks could account for part of this interspecific divergence. Both of these putative QTLs were mapped at the same cytological regions as other QTLs for intraspecific wing size variation identified in D. melanogaster studies. In these regions, one or more loci could account for intra- and interspecific variation in the size of Drosophila wings. Three other QTL peaks were related to a pattern of interspecific variation in wing size and shape traits that is summarized by one principal component. In addition, we observed that female wings are significantly larger and longer than male wings and the second, fourth and fifth longitudinal veins are closer together at the distal wing area. This pattern was summarized by another principal component, for which one QTL was mapped.     

Different cell size and cell number contribution in two newly established and one ancient body size cline of Drosophila subobscura

Latitudinal genetic clines in body size occur in many ectotherms including Drosophila species. In the wing of D. melanogaster, these clines are generally based on latitudinal variation in cell number. In contrast, differences in wing area that evolve by thermal selection in the laboratory are in general based on cell size. To investigate possible reasons for the different cellular bases of these two types of evolutionary response, we compared the newly established North and South American wing size clines of Drosophila subobscura. The new clines are based on latitudinal variation in cell area in North America and cell number in South America. The ancestral European cline is also based on latitudinal variation in cell number. The difference in the cellular basis of wing size variation in the American clines, which are roughly the same age, together with the similar cellular basis of the new South American cline and the ancient European one, suggest that the antiquity of a cline does not explain its cellular basis. Furthermore, the results indicate that wing size as a whole, rather than its cellular basis, is under selection. The different cellular bases of different size clines are most likely explained either entirely by chance or by different patterns of genetic variance--or its expression--in founding populations.     

Comparative Analysis of Morphological Traits Among Drosophila Melanogaster and D. Simulans: Genetic Variability, Clines and Phenotypic Plasticity

The two sibling cosmopolitan species, Drosophila melanogaster and D. simulans, are able to proliferate under very different climatic conditions. This has resulted in local adaptations, which are often arranged in latitudinal clines. Such clines are documented for body weight, wing and thorax length, sternopleural and abdominal bristle number, ovariole number and thoracic pigmentation. The overall magnitude of geographical differentiation is, however, much less in D. simulans than in D. melanogaster, and latitudinal clines are less pronounced. The fact that natural populations live under different climates raises the problem of interaction between temperature and phenotype. The reaction norms of morphometrical traits have been investigated as a function of growth temperature. The shapes of the response curves vary according to the investigated trait. They are generally curvilinear and can be described by calculating characteristic values after polynomial adjustments. For a given trait, the reaction norms of the two species are similar in their shape, although some significant differences may be observed. Within each species, significant differences are also observed between geographic populations: reaction norms are not parallel and the divergence is better marked when more distant populations (e.g., temperate and tropical) are compared. It thus appears that besides mean trait value, phenotypic plasticity is also a target of natural selection. A specific analysis of wing shape variation according to growth temperature was also undertaken. Reaction norms with different shapes may be observed in various parts of the wing: the major effect is found between the basis and the tip of the wing, but in a similar way in the two species. By contrast, some ratios, called wing indices by taxonomists, may exhibit completely different reaction norms in the two species. For a single developmental temperature (25 degrees C) the phenotypic variability of morphometrical traits is generally similar in the two species, and also the genetic variability, estimated by the intraclass correlation. A difference exists, however, for the ovariole number which is less variable in D. simulans. Variance parameters may vary according to growth temperature, and a detailed analysis was made on wing dimensions. An increase of environmental variability at extreme, heat or cold temperatures, has been found in both species. Opposite trends were, however, observed for the genetic variability: a maximum heritability in D. simulans at middle temperatures, corresponding to a minimum heritability in D. melanogaster. Whether such a difference exists for other traits and in other populations deserves further investigations. In conclusion, morphometrical analyses reveal a large amount of significant differences which may be related to speciation and to the divergence of ecological niches. Within each species, numerous geographic variations are also observed which, in most cases, reflect some kinds of climatic adaptation.     

Contrasting patterns of phenotypic variation linked to chromosomal inversions in native and colonizing populations of Drosophila subobscura

In fewer than two decades after invading the Americas, the fly Drosophila subobscura evolved latitudinal clines for chromosomal inversion frequencies and wing size that are parallel to the long‐standing ones in native Palearctic populations. By sharp contrast, wing shape clines also evolved in the New World, but the relationship with latitude was opposite to that in the Old World. Previous work has suggested that wing trait differences among individuals are partially due to the association between chromosomal inversions and particular alleles which influence the trait under consideration. Furthermore, it is well documented that a few number of effective individuals founded the New World populations, which might have modified the biometrical effect of inversions on quantitative traits. Here we evaluate the relative contribution of chromosomal inversion clines in shaping the parallel clines in wing size and contrasting clines in wing shape in native and colonizing populations of the species. Our results reveal that inversion‐size and inversion‐shape associations in native and colonizing (South America) populations are generally different, probably due to the bottleneck effect. Contingent, unpredictable evolution was suggested as an explanation for the different details involved in the otherwise parallel wing size clines between Old and New World populations of D. subobscura. We challenge this assertion and conclude that contrasting wing shape clines came out as a correlated response of inversion clines that might have been predicted considering the genetic background of colonizers.     

Wing trait–inversion associations in Drosophila subobscura can be generalized within continents,but may change through time

Clinal variation is one of the most emblematic examples of the action of natural selection at a wide geographical range. In Drosophila subobscura, parallel clines in body size and inversions, but not in wing shape, were found in Europe and South and North America. Previous work has shown that a bottleneck effect might be largely responsible for differences in wing trait–inversion association between one European and one South American population. One question still unaddressed is whether the associations found before are present across other populations of the European and South American clines. Another open question is whether evolutionary dynamics in a new environment can lead to relevant changes in wing traits–inversion association. To analyse geographical variation in these associations, we characterized three recently laboratory founded D. subobscura populations from both the European and South American latitudinal clines. To address temporal variation, we also characterized the association at a later generation in the European populations. We found that wing size and shape associations can be generalized across populations of the same continent, but may change through time for wing size. The observed temporal changes are probably due to changes in the genetic content of inversions, derived from adaptation to the new, laboratory environment. Finally, we show that it is not possible to predict clinal variation from intrapopulation associations. All in all this suggests that, at least in the present, wing traits–inversion associations are not responsible for the maintenance of the latitudinal clines in wing shape and size.     

THE GENETIC ARCHITECTURE OF WING SIZE DIVERGENCE AT VARYING SPATIAL SCALES ALONG A BODY SIZE CLINE IN DROSOPHILA MELANOGASTER

Latitudinal clines in quantitative traits are common, but surprisingly little is known about the genetic bases of these divergences and how they vary within and between clines. Here, we use line‐cross analysis to investigate the genetic architecture of wing size divergences at varying spatial scales along a body size cline in Drosophila melanogaster. Our results revealed that divergences in wing size along the cline were due to strong additive effects. Significant nonadditive genetic effects, including epistasis and maternal effects, were also detected, but they were relatively minor in comparison to the additive effects and none were common to all crosses. There was no evidence of increased epistasis in crosses between more geographically distant populations and, unlike in previous studies, we found no significant dominance effects on wing size in any cross. Our results suggest there is little variation in the genetic control of wing size along the length of the Australian cline. They also highlight marked inconsistencies in the magnitude of dominance effects across studies, which may reflect different opportunities for mutation accumulation while lines are in laboratory culture.     

Do quantitative trait loci (QTL) for a courtship song difference between Drosophila simulans and D. sechellia coincide with candidate genes and intraspecific QTL?

The genetic architecture of traits influencing sexual isolation can give insight into the evolution of reproductive isolation and hence speciation. Here we report a quantitative trait loci (QTL) analysis of the difference in mean interpulse interval (IPI), an important component of the male courtship song, between Drosophila simulans and D. sechellia. Using a backcross analysis, we find six QTL that explain a total of 40.7% of the phenotypic variance. Three candidate genes are located in the intervals bounded by two of the QTL and there are no significant QTL on the X chromosome. The values of mean IPI for hybrid individuals imply the presence of dominant alleles or epistasis. Because unisexual hybrid sterility prevents an F(2) analysis, we cannot distinguish dominant from additive genetic effects at the scale of QTL. A comparison with a study of QTL for intraspecific variation in D. melanogaster shows that, for these strains, the QTL we have identified for interspecific variation cannot be those that contribute to intraspecific variation. We find that the QTL have bidirectional effects, which indicates that the genetic architecture is compatible with divergence due to genetic drift, although other possibilities are discussed.     

Climatic selection on genes and traits after a 100 year-old invasion: a critical look at the temperate-tropical clines in Drosophila melanogaster from eastern Australia

Drosophila melanogaster invaded Australia around 100 years ago, most likely through a northern invasion. The wide range of climatic conditions in eastern Australia across which D. melanogaster is now found provides an opportunity for researchers to identify traits and genes that are associated with climatic adaptation. Allozyme studies indicate clinal patterns for at least four loci including a strong linear cline in Adh and a non-linear cline in alpha-Gpdh. Inversion clines were initially established from cytological studies but have now been validated with larger sample sizes using molecular markers for breakpoints. Recent collections indicate that some genetic markers (Adh and In(3R)Payne) have changed over the last 20 years reflecting continuing evolution. Heritable clines have been established for quantitative traits including wing length/area, thorax length and cold and heat resistance. A cline in egg size independent of body size and a weak cline in competitive ability have also been established. Postulated clinal patterns for resistance to desiccation and starvation have not been supported by extensive sampling. Experiments under laboratory and semi-natural conditions have suggested selective factors generating clinal patterns, particularly for reproductive patterns over winter. Attempts are being made to link clinal variation in traits to specific genes using QTL analysis and the candidate locus approach, and to identify the genetic architecture of trait variation along the cline. This is proving difficult because of inversion polymorphisms that generate disequilibrium among genes. Substantial gaps still remain in linking clines to field selection and understanding the genetic and physiological basis of the adaptive shifts. However D. melanogaster populations in eastern Australia remain an excellent resource for understanding past and future evolutionary responses to climate change.     

Naturally segregating quantitative trait loci affecting wing shape of Drosophila melanogaster

Variation in vein position and wing shape of Drosophila melanogaster depends on many genes. In the following, we report the results of a QTL analysis of wing shape in D. melanogaster. We identified QTL responsible for natural variation for wing shape and analyzed their interactions with developmental genetic signaling pathways important for vein positioning. The QTL analysis indicated that the total number of QTL segregating in this population is likely to be very large. The locations of putative QTL identified in this study were compared to those identified in previous studies and, while there is more correspondence across studies than expected by chance on the third chromosome, the studies appear to have identified different QTL. Using a complementation design, we tested for interactions among these QTL with the Hedgehog and Decapentaplegic signaling pathways, which are important for the development and position of vein pairs L3-L4 and L2-L5. Three QTL showed strong interactions with these two pathways, supporting the hypothesis that these QTL are involved in these pathways. Naturally segregating variation can therefore act through known signaling pathways to produce variation in vein position.     

A comparison of the genetic basis of wing size divergence in three parallel body size clines of Drosophila melanogaster

Body size clines in Drosophila melanogaster have been documented in both Australia and South America, and may exist in Southern Africa. We crossed flies from the northern and southern ends of each of these clines to produce F(1), F(2), and first backcross generations. Our analysis of generation means for wing area and wing length produced estimates of the additive, dominance, epistatic, and maternal effects underlying divergence within each cline. For both females and males of all three clines, the generation means were adequately described by these parameters, indicating that linkage and higher order interactions did not contribute significantly to wing size divergence. Marked differences were apparent between the clines in the occurrence and magnitude of the significant genetic parameters. No cline was adequately described by a simple additive-dominance model, and significant epistatic and maternal effects occurred in most, but not all, of the clines. Generation variances were also analyzed. Only one cline was described sufficiently by a simple additive variance model, indicating significant epistatic, maternal, or linkage effects in the remaining two clines. The diversity in genetic architecture of the clines suggests that natural selection has produced similar phenotypic divergence by different combinations of gene action and interaction.     

The genetic basis of the interspecific differences in wing size in Nasonia (Hymenoptera; Pteromalidae): major quantitative trait loci and epistasis

There is a 2.5-fold difference in male wing size between two haplodiploid insect species, Nasonia vitripennis and N. giraulti. The haploidy of males facilitated a full genomic screen for quantitative trait loci (QTL) affecting wing size and the detection of epistatic interactions. A QTL analysis of the interspecific wing-size difference revealed QTL with major effects and epistatic interactions among loci affecting the trait. We analyzed 178 hybrid males and initially found two major QTL for wing length, one for wing width, three for a normalized wing-size variable, and five for wing seta density. One QTL for wing width explains 38.1% of the phenotypic variance, and the same QTL explains 22% of the phenotypic variance in normalized wing size. This corresponds to a region previously introgressed from N. giraulti into N. vitripennis that accounts for 44% of the normalized wing-size difference between the species. Significant epistatic interactions were also found that affect wing size and density of setae on the wing. Screening for pairwise epistatic interactions between loci on different linkage groups revealed four additional loci for wing length and four loci for normalized wing size that were not detected in the original QTL analysis. We propose that the evolution of smaller wings in N. vitripennis males is primarily the result of major mutations at few genomic regions and involves epistatic interactions among some loci.     

Mapping genes that control hormone-induced ovulation rate in mice.

The present study mapped quantitative trait loci (QTL) that control 6-fold genetic differences in hormone-induced ovulation rate (HIOR) between C57BL/6J (B6) (HIOR = 54) and A/J strain mice (HIOR = 9). (The gene name is Ovulation Rate Induced [ORI] QTL and the gene symbol is Oriq.) QTL linkage analysis was conducted on 167 (B6xA)xA backcross mice at 165 loci. Suggestive B6 ORI QTL that control the number of eggs in cumulus mapped, as follows, near: Cyp19 and D9Mit4 on chromosome (Chr) 9 (Oriq1); D2Mit433 on Chr2 (Oriq2); D6Mit316 on Chr6 (Oriq3); DXMit22 on ChrX (Oriq4) and were associated with a 2.7, 2.7, 2.6, and 4.2 egg increases in HIOR, respectively. Oriq3 was significant (LOD = 3.45) based on composite interval mapping. QTL linkage analysis of the number of eggs matured by endogenous gonadotropins and ovulated by eCG mapped a significant Oriq5 to Chr 10 and suggestive Oriq to Chr 6, 7, and X. These data provide the first molecular genetic markers for reproductive QTL that control major differences in ovarian responsiveness to gonadotropins. These and closely linked syntenic molecular markers will enable a more accurate prediction of ovarian responsiveness to gonadotropins and provide selection criteria for improving reproductive performance in diverse mammalian species.     

QTL linkage mapping of wing length in zebra finch using genome-wide single nucleotide polymorphisms markers

Avian wing length is an important trait that covaries with the ecology and migratory behaviour of a species and tends to change rapidly when the conditions are altered. Long-distance migrants typically have longer wings than short-distance migrants and sedentary species, and long-winged species also tend to be more dispersive. Although the substantial heritability of avian wing length is well established, the identification of causal genes has remained elusive. Based on large-scale genotyping of 1404 informative single nucleotide polymorphisms (SNP) in a captive population of 1067 zebra finches, we here show that the within-population variation of relative wing length (h(2) = 0.74 ± 0.05) is associated with standing genetic variation in at least six genomic regions (one genome-wide significant and five suggestive). The variance explained by these six quantitative trait loci (QTL) sums to 36.8% of the phenotypic variance (half of the additive genetic variance), although this likely is an overestimate attributable to the Beavis effect. As avian wing length is primarily determined by the length of the primary feathers, we then searched for candidate genes that are related to feather growth. Interestingly, all of the QTL signals co-locate with Wnt growth factors and closely interacting genes (Wnt3a, Wnt5a, Wnt6, Wnt7a, Wnt9a, RhoU and RhoV). Our findings therefore suggest that standing genetic variation in the Wnt genes might be linked to avian wing morphology, although there are many other genes that also fall within the confidence regions.     

Quantitative trait loci for sexual isolation between Drosophila simulans and D. mauritiana

Sexual isolating mechanisms that act before fertilization are often considered the most important genetic barriers leading to speciation in animals. While recent progress has been made toward understanding the genetic basis of the postzygotic isolating mechanisms of hybrid sterility and inviability, little is known about the genetic basis of prezygotic sexual isolation. Here, we map quantitative trait loci (QTL) contributing to prezygotic reproductive isolation between the sibling species Drosophila simulans and D. mauritiana. We mapped at least seven QTL affecting discrimination of D. mauritiana females against D. simulans males, three QTL affecting D. simulans male traits against which D. mauritiana females discriminate, and six QTL affecting D. mauritiana male traits against which D. simulans females discriminate. QTL affecting sexual isolation act additively, are largely different in males and females, and are not disproportionately concentrated on the X chromosome: The QTL of greatest effect are located on chromosome 3. Unlike the genetic components of postzygotic isolation, the loci for prezygotic isolation do not interact epistatically. The observation of a few QTL with moderate to large effects will facilitate positional cloning of genes underlying sexual isolation.     

Explaining the heritability of an ecologically significant trait in terms of individual quantitative trait loci

Most natural populations display substantial genetic variation in behaviour, morphology, physiology, life history and the susceptibility to disease. A major challenge is to determine the contributions of individual loci to variation in complex traits. Quantitative trait locus (QTL) mapping has identified genomic regions affecting ecologically significant traits of many species. In nearly all cases, however, the importance of these QTLs to population variation remains unclear. In this paper, we apply a novel experimental method to parse the genetic variance of floral traits of the annual plant Mimulus guttatus into contributions of individual QTLs. We first use QTL-mapping to identify nine loci and then conduct a population-based breeding experiment to estimate V(Q), the genetic variance attributable to each QTL. We find that three QTLs with moderate effects explain up to one-third of the genetic variance in the natural population. Variation at these loci is probably maintained by some form of balancing selection. Notably, the largest effect QTLs were relatively minor in their contribution to heritability.     

ALTITUDINAL CLINAL VARIATION IN WING SIZE AND SHAPE IN AFRICAN DROSOPHILA MELANOGASTER: ONE CLINE OR MANY?

Geographical patterns of morphological variation have been useful in addressing hypotheses about environmental adaptation. In particular, latitudinal clines in phenotypes have been studied in a number of Drosophila species. Some environmental conditions along latitudinal clines—for example, temperature—also vary along altitudinal clines, but these have been studied infrequently and it remains unclear whether these environmental factors are similar enough for convergence or parallel evolution. Most clinal studies in Drosophila have dealt exclusively with univariate phenotypes, allowing for the detection of clinal relationships, but not for estimating the directions of covariation between them. We measured variation in wing shape and size in D. melanogaster derived from populations at varying altitudes and latitudes across sub‐Saharan Africa. Geometric morphometrics allows us to compare shape changes associated with latitude and altitude, and manipulating rearing temperature allows us to quantify the extent to which thermal plasticity recapitulates clinal effects. Comparing effect vectors demonstrates that altitude, latitude, and temperature are only partly associated, and that the altitudinal shape effect may differ between Eastern and Western Africa. Our results suggest that selection responsible for these phenotypic clines may be more complex than just thermal adaptation.     

Molecular mapping of a new source of Fusarium wilt resistance in tetraploid cotton (Gossypium hirsutum L.)

Fusarium wilt (FW) disease is an economically important disease of cotton worldwide and a major cause of crop losses in Australia and many other cotton-producing countries. Symptoms include wilting, vascular browning and death. Australian races of the causal agent Fusarium oxysporum f. sp. vasinfectum (Fov) are genetically distinct from those in other countries and are thought to have evolved from indigenous races. New sources of resistance for breeding are rare, as cotton cultivars with significant FW resistance against Fov isolates from other cotton-producing regions are usually susceptible to Australian Fov races. MCU-5, an Upland Indian cotton cultivar, has been identified as having improved resistance to Australian Fov and is being used to breed new commercial cultivars with higher resistance to FW. To investigate the genetic basis of the FW resistance in MCU-5, QTL analysis was performed on 244 F3 and 244 F4 families derived from an intraspecific cross between MCU-5 and Siokra 1-4, a cultivar highly sensitive to Australian Fov races. Resistance, as measured by leaf symptoms, vascular browning and survival, showed low to moderate heritability between generations. MCU-5 resistance to FW was found to be complex with three quantitative trait loci (QTL) identified in the F3, and eight in the F4, that explained between 9 and 41% of the phenotypic variation. The QTL were located on four linkage groups including chromosomes A6 (Chr 6), D4 (Chr 22) and D6 (Chr 25), with two QTL located in similar regions to previously identified FW resistance from the Sea Island cultivar Pima 3-79. The QTL identified in this study represent the first targets for marker-assisted selection of FW resistance in Australia.     

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype−environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.