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.     

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.     

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.     

Swift laboratory thermal evolution of wing shape (but not size) in Drosophila subobscura and its relationship with chromosomal inversion polymorphism

Latitudinal clinal variation in wing size and shape has evolved in North American populations of Drosophila subobscura within about 20 years since colonization. While the size cline is consistent to that found in original European populations (and globally in other Drosophila species), different parts of the wing have evolved on the two continents. This clearly suggests that 'chance and necessity' are simultaneously playing their roles in the process of adaptation. We report here rapid and consistent thermal evolution of wing shape (but not size) that apparently is at odds with that suggestion. Three replicated populations of D. subobscura derived from an outbred stock at Puerto Montt (Chile) were kept at each of three temperatures (13, 18 and 22 degrees C) for 1 year and have diverged for 27 generations at most. We used the methods of geometric morphometrics to study wing shape variation in both females and males from the thermal stocks, and rates of genetic divergence for wing shape were found to be as fast or even faster than those previously estimated for wing size on a continental scale. These shape changes did not follow a neat linear trend with temperature, and are associated with localized shifts of particular landmarks with some differences between sexes. Wing shape variables were found to differ in response to male genetic constitution for polymorphic chromosomal inversions, which strongly suggests that changes in gene arrangement frequencies as a response to temperature underlie the correlated changes in wing shape because of gene-inversion linkage disequilibria. In fact, we also suggest that the shape cline in North America likely predated the size cline and is consistent with the quite different evolutionary rates between inversion and size clines. These findings cast strong doubts on the supposed 'unpredictability' of the geographical cline for wing traits in D. subobscura North American colonizing populations.     

Population differentiation in the beetleTribolium castaneum. I. Genetic architecture

We used joint-scaling analyses in conjunction with rearing temperature variation to investigate the contributions of additive, non-additive, and environmental effects to genetic divergence and incipient speciation among 12 populations of the red flour beetle, Tribolium castaneum, with small levels of pairwise nuclear genetic divergence (0.033 < Nei's D < 0.125). For 15 population pairs we created a full spectrum of line crosses (two parental, two reciprocal F1's, four F2's, and eight backcrosses), reared them at multiple temperatures, and analyzed the numbers and developmental defects of offspring. We assayed a total of 219,388 offspring from 5147 families. Failed crosses occurred predominately in F2's, giving evidence of F2 breakdown within this species. In all cases where a significant model could be fit to the data on offspring number, we observed at least one type of digenic epistasis. We also found maternal and cytoplasmic effects to be common components of divergence among T. castaneum populations. In some cases, the most complex model tested (additive, dominance, epistatic, maternal, and cytoplasmic effects) did not provide a significant fit to the data, suggesting that linkage or higher order epistasis is involved in differentiation between some populations. For the limb deformity data, we observed significant genotype-by-environment interaction in most crosses and pure parent crosses tended to have fewer deformities than hybrid crosses. Complexity of genetic architecture was not correlated with either geographic distance or genetic distance. Our results support the view that genetic incompatibilities responsible for postzygotic isolation, an important component of speciation, may be a natural but serendipitous consequence of nonadditive genetic effects and structured populations.     

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.     

Rapid evolution of wing size clines in Drosophila subobscura

Parallel latitudinal clines across species and continents provide dramatic evidence of the efficacy of natural selection, however little is known about the dynamics involved in cline formation. For example, several drosophilids and other ectotherms increase in body and wing size at higher latitudes. Here we compare evolution in an ancestral European and a recently introduced (North America) cline in wing size and shape in Drosophila subobscura. We show that clinal variation in wing size, spanning more than 15 degrees of latitude, has evolved in less than two decades. In females from Europe and North America, the clines are statistically indistinguishable however the cline for North American males is significantly shallower than that for European males. We document that while overall patterns of wing size are similar on two continents, the European cline is obtained largely through changing the proximal portion of the wing, whereas the North American cline is largely in the distal portion. We use data from sites collected in 1986/1988 (Pegueroles et al. 1995) and our 1997 collections to compare synchronic (divergence between contemporary populations that share a common ancestor) and allochronic (changes over time within a population) estimates of the rates of evolution. We find that, for these populations, allochronically estimated evolutionary rates within a single population are over 0.02 haldanes (2800 darwins), a value similar in magnitude to the synchronic estimates from the extremes of the cline. This paper represents an expanded analysis of data partially presented in Huey et al. (2000).     

B chromosome polymorphism in maize landraces: adaptive vs. demographic hypothesis of clinal variation

Cytogenetic analysis of maize landraces from northwestern Argentina has revealed an altitudinal cline in the mean number of B chromosomes (B's) per plant, with cultivars growing at higher altitudes exhibiting a higher number of B's. Altitudinal and longitudinal clines are frequently interpreted as evidence of selection, however, they can also be produced by the interplay between drift and spatially restricted gene flow or by admixture between previously isolated populations that have come into secondary contact. Here, we test the adaptive significance of the observed altitudinal gradient by comparing the levels of differentiation in the mean number of B's to those obtained from 18 selectively neutral loci [simple sequence repeats (SSRs)] among seven populations of the cline. The adequacy of alternative genetic-differentiation measures was determined, and associations between cytogenetic, genetic, and altitudinal distances were assessed by means of matrix- correspondence tests. No evidence for association between pairwise F(ST) and altitudinal distance or B-chromosome differentiation was found. The contrasting pattern of altitudinal divergence between the mean number of B's per plant and the genetic differentiation at SSR loci indicates that demographic processes cannot account for the observed levels of divergence in the mean number of B's.     

Epistasis and the Genetic Divergence of Photoperiodism between Populations of the Pitcher-Plant Mosquito, Wyeomyia Smithii

Parallel crosses between each of two southern (ancestral) and one northern (derived) population of the pitcher-plant mosquito, Wyeomyia smithii, were made to determine the genetic components of population divergence in critical photoperiod, a phenological trait that measures adaptation to seasonality along a climatic gradient. Joint scaling tests were used to analyze means and variances of first- and second-generation hybrids in order to determine whether nonadditive genetic variance, especially epistatic variance, contributed to divergence in critical photoperiod. In both crosses, digenic epistatic effects were highly significant, indicating that genetic divergence cannot have resulted solely from differences in additively acting loci. For one cross that could be tested directly for such effects, higher order epistasis and/or linkage did not contribute to the divergence of critical photoperiod between the constituent populations.     

Clines in polygenic traits

This article outlines theoretical models of clines in additive polygenic traits, which are maintained by stabilizing selection towards a spatially varying optimum. Clines in the trait mean can be accurately predicted, given knowledge of the genetic variance. However, predicting the variance is difficult, because it depends on genetic details. Changes in genetic variance arise from changes in allele frequency, and in linkage disequilibria. Allele frequency changes dominate when selection is weak relative to recombination, and when there are a moderate number of loci. With a continuum of alleles, gene flow inflates the genetic variance in the same way as a source of mutations of small effect. The variance can be approximated by assuming a Gaussian distribution of allelic effects; with a sufficiently steep cline, this is accurate even when mutation and selection alone are better described by the 'House of Cards' approximation. With just two alleles at each locus, the phenotype changes in a similar way: the mean remains close to the optimum, while the variance changes more slowly, and over a wider region. However, there may be substantial cryptic divergence at the underlying loci. With strong selection and many loci, linkage disequilibria are the main cause of changes in genetic variance. Even for strong selection, the infinitesimal model can be closely approximated by assuming a Gaussian distribution of breeding values. Linkage disequilibria can generate a substantial increase in genetic variance, which is concentrated at sharp gradients in trait means.     

The contrasting genetic architecture of wing size, viability, and development time in a rainforest species and its more widely distributed relative

Divergence among populations can occur via additive genetic effects and/or because of epistatic interactions among genes. Here we use line-cross analysis to compare the importance of epistasis in divergence among two sympatric Drosophila species from eastern Australia, one (D. serrata) distributed continuously and the other (D. birchii) confined to rainforest habitats that are often disjunct. For D. serrata, crosses indicated that development time and wing size differences were due to additive genetic effects, while for viability there were digenic epistatic effects. Crosses comparing geographically close populations as well as those involving the most geographically distant populations (including the southern species border) revealed epistatic interactions, whereas crosses at an intermediate distance showed no epistasis. In D. birchii, there was no evidence of epistasis for viability, although for development time and wing size there was epistasis in the cross between the most geographically diverged populations. Strong epistasis has not developed among the D. birchii populations, and this habitat specialist does not show stronger epistasis than D. serrata. Given that epistasis has been detected in crosses with other species from eastern Australia, including the recently introduced D. melanogaster, the results point to epistasis not being directly linked to divergence times among populations.     

Inference from Clines Stabilized by Frequency-Dependent Selection

Frequency-dependent selection against rare forms can maintain clines. For weak selection, s, in simple linear models of frequency-dependence, single locus clines are stabilized with a maximum slope of between square root of s/square root of 8 sigma and square root of s/square root of 12 delta, where sigma is the dispersal distance. These clines are similar to those maintained by heterozygote disadvantage. Using computer simulations, the weak-selection analytical results are extended to higher selection pressures with up to three unlinked genes. Graphs are used to display the effect of selection, migration, dominance, and number of loci on cline widths, speeds of cline movements, two-way gametic correlations ("linkage disequilibria"), and heterozygote deficits. The effects of changing the order of reproduction, migration, and selection, are also briefly explored. Epistasis can also maintain tension zones. We show that epistatic selection is similar in its effects to frequency-dependent selection, except that the disequilibria produced in the zone will be higher for a given level of selection. If selection consists of a mixture of frequency-dependence and epistasis, as is likely in nature, the error made in estimating selection is usually less than twofold. From the graphs, selection and migration can be estimated using knowledge of the dominance and number of genes, of gene frequencies and of gametic correlations from a hybrid zone.     

Adaptation to a steep environmental gradient and an associated barrier to gene exchange in Littorina saxatilis

Steep environmental gradients offer important opportunities to study the interaction between natural selection and gene flow. Allele frequency clines are expected to form at loci under selection, but unlinked neutral alleles may pass easily across these clines unless a generalized barrier evolves. Here we consider the distribution of forms of the intertidal gastropod Littorina saxatilis, analyzing shell shape and amplified fragment length polymorphism (AFLP) loci on two rocky shores in Britain. On the basis of previous work, the AFLP loci were divided into differentiated and undifferentiated groups. On both shores, we have shown a sharp cline in allele frequencies between the two morphs for differentiated AFLP loci. This is coincident with a habitat transition on the shore where the two habitats (cliff and boulder field) are immediately contiguous. The allele frequency clines coincide with a cline in shell morphology. In the middle of the cline, linkage disequilibrium for the differentiated loci rises in accordance with expectation. The clines are extremely narrow relative to dispersal, probably as a result of both strong selection and habitat choice. An increase in F(ST) for undifferentiated AFLPs between morphs, relative to within-morph comparisons, is consistent with there being a general barrier to gene flow across the contact zone. These features are consistent either with an episode of allopatric divergence followed by secondary contact or with primary, nonallopatric divergence. Further data will be needed to distinguish between these alternatives.     

Estimates of genetic variance in an F2 maize population

Maize (Zea mays L.) breeders have used several genetic-statistical models to study the inheritance of quantitative traits. These models provide information on the importance of additive, dominance, and epistatic genetic variance for a quantitative trait. Estimates of genetic variances are useful in understanding heterosis and determining the response to selection. The objectives of this study were to estimate additive and dominance genetic variances and the average level of dominance for an F2 population derived from the B73 x Mo17 hybrid and use weighted least squares to determine the importance of digenic epistatic variances relative to additive and dominance variances. Genetic variances were estimated using Design III and weighted least squares analyses. Both analyses determined that dominance variance was more important than additive variance for grain yield. For other traits, additive genetic variance was more important than dominance variance. The average level of dominance suggests either overdominant gene effects were present for grain yield or pseudo-overdominance because of linkage disequilibrium in the F2 population. Epistatic variances generally were not significantly different from zero and therefore were relatively less important than additive and dominance variances. For several traits estimates of additive by additive epistatic variance decreased estimates of additive genetic variance, but generally the decrease in additive genetic variance was not significant.     

家蚕茧质性状的QTL定位研究

采用QTLMapper 2．0 QTL作图软件,对F2群体的家蚕全茧量、茧层量、茧层率和蛹体重等性状进行了QTL定位分析,分别检测出7个、6个、2个、8个有显著效应分量的QTLs,分布于7个、5个、2个、7个不同的连锁群。控制全茧量、茧层量的QTLs一般存在复杂的上位性效应。对全茧量性状,有3对QTLs存在显著的加加上位性效应,其中1对还存在加显、显显互作;共有3个QTLs存在显著的显性效应,1个存在显著的加性效应。对茧层量QTLs,发现1对QTLs存在极显著的各项遗传效应,包括上位性效应;1对QTLs被检测到显著的显显互作,1个QTL具有显著的显性效应,并与另一个QTL存在显著的加加互作。茧层率、蛹体重主要受加性或显性的QTLs作用,没有发现茧层率QTLs的上位性效应,蛹体重的有效QTL大都呈现显著的负向显性效应,只有一对QTLs存在显著的加加上位性效应。第2、3、4、11、13、24、34、37、40连锁群是两个或多个性状QTLs分布的共同连锁群。全茧量和茧层量存在共同的QTL或染色体区域,育种上可通过适当选配,利用基因的互作效应,同步改良这两个性状。     

Genetic divergence in adaptive characters between sympatric species of stickleback

This study explored the genetic basis of phenotypic differences between two sympatric species of ecologically and morphologically divergent sticklebacks (Gasterosteus aculeatus complex). The aim was to understand how many loci determine the differences and to what extent the differences are due to additive or nonadditive gene action. I reared the two parental species, F1 and F2 hybrids, and both backcrosses in the laboratory and measured the following quantitative characters: gill raker number and length (both involved in feeding), lateral plate number and pelvic spine length (both involved in predator defense), and growth (a fitness component). I then applied joint-scaling regression models to estimate composite additive, dominance and epistatic effects, and their contribution to divergence of parental lines. A simple additive model was sufficient for gill raker number and growth; additive and dominance effects contributed significantly to divergence in plate number and pelvic spine length; and additive, dominance, and epistatic effects contributed significantly to divergence in gill raker length. Wright's estimator for the number of loci for the four morphological characters ranged from 1 to 50. My results suggest that adaptive divergence between limnetic and benthic sticklebacks has taken place through a variety of genetic mechanisms specific to different traits. Though interspecific hybrids are completely fertile and viable in the laboratory, they are selected against in the wild. The pattern of inheritance for the traits examined here directly influences how well hybrids can exploit the two major resource environments in the wild.     

Genetic covariances promote climatic adaptation in Australian Drosophila*

Evolutionary potential for adaptation hinges upon the orientation of genetic variation for traits under selection, captured by the additive genetic variance-covariance matrix (G), as well as the evolutionary stability of G. Yet studies that assess both the stability of G and its alignment with selection are extraordinarily rare. We evaluated the stability of G in three Drosophila melanogaster populations that have adapted to local climatic conditions along a latitudinal cline. We estimated population- and sex-specific G matrices for wing size and three climatic stress-resistance traits that diverge adaptively along the cline. To determine how G affects evolutionary potential within these populations, we used simulations to quantify how well G aligns with the direction of trait divergence along the cline (as a proxy for the direction of local selection) and how genetic covariances between traits and sexes influence this alignment. We found that G was stable across the cline, showing no significant divergence overall, or in sex-specific subcomponents, among populations. G also aligned well with the direction of clinal divergence, with genetic covariances strongly elevating evolutionary potential for adaptation to climatic extremes. These results suggest that genetic covariances between both traits and sexes should significantly boost evolutionary responses to environmental change.     

Parapatric diversification after post-glacial range expansion in the gall fly Urophora cardui (Tephritidae)

Aim Primary and secondary genetic clines in post-glacial colonized regions have different implications for biogeographic distributions and the origin of species. Primary clines arise in situ after colonization as adaptive responses to environmental gradients, while secondary clines are caused by contact between vicariant lineages. Here we analyse primary versus secondary origin of a genetic cline in the tephritid fly Urophora cardui in Jutland, Denmark, in a post-glacial landscape. Location Western Palaearctic. Methods Phylogeographic and demographic analyses of U. cardui based on mitochondrial DNA (mtDNA) genealogies, hierarchical genetic variance tests based on allozymes and distribution analysis of a rare allele from the Jutland cline. Results There was no phylogeographic divergence between the Jutland population of U. cardui north of the cline and neighbouring western European regional populations, which all shared the common western European mtDNA haplotype H1. At nuclear loci, by contrast, the North Jutland population was diverged above the mean level of divergence among regional populations and had no loss of genetic variation. A rare allozyme allele that was frequent in the cline area (up to 45%) and was missing north of the cline also occurred at low frequency (0–14%) elsewhere in the sampling range. Shallow phylogeographic divergence was observed between Russian and western European populations and between English and continental populations. Main conclusions The genetic variation patterns support primary cline evolution and parapatric divergence in Jutland following a demographic expansion of a western European ancestral source population of U. cardui, and suggest cryptic refugia and/or selection in other European population assemblages. The patterns of intra-specific regional divergence are discussed with respect to the interpretation of cryptic refugia in Europe after the most recent ice age.     

The quantitative genetics of fluctuating asymmetry

Fluctuating asymmetry (subtle departures from identical expression of a trait across an axis of symmetry) in many taxa is under stabilizing selection for reduced asymmetry. However, lack of reliable estimates of genetic parameters for asymmetry variation hampers our ability to predict the evolutionary outcome of this selection. Here we report on a study, based on analysis of variation within and between isofemale lines and of generation means (line-cross analysis), designed to dissect in detail the quantitative genetics of positional fluctuating asymmetry (PFA) in bristle number in natural populations of Drosophila falleni. PFA is defined as the difference between the two sides of the body in the placement or position of components of a meristic trait. Heritability (measured at 25 degrees C) of two related measures of PFA were 13% and 21%, both of which differed significantly from zero. In contrast, heritability estimates for fluctuating asymmetry in the total number of anterior (0.7%) and transverse (2.4%) sternopleural bristles were smaller, not significant, and in quantitative agreement with previously published estimates. Heritabilities for bristle number (trait size) were considerably greater than that for any asymmetry measure. The experimental design controlled for the potentially confounding effects of common familial environment, and repeated testing revealed that PFA differences between lines were genetically stable for up to 16 generations in the laboratory at 25 degrees C. We performed line cross analysis between strains at the extremes of the PFA distribution (highest and lowest values); parental strains, F1, F1r (reciprocal), F2, backcross, and backcross reciprocal generations were represented. The inheritance of PFA was described best by additive and dominance effects localized to the X-chromosomes, whereas autosomal dominance effects were also detected. Epistatic, maternal, and cytoplasmic effects were not detected. The inheritance of trait size was notably more complex and involved significant autosomal additive, dominance, and epistatic effects; maternal dominance effects; and additive and dominance effects localized to the X-chromosomes. The additive genetic correlation between PFA and its associated measure of trait size was negative (-0.049), but not statistically significant, indicating that the loci contributing additive genetic effects to these traits are probably different. It is suggested that PFA may be a sensitive measure of developmental instability because PFA taps the ability of an organism to integrate interconnected developmental pathways.     

The mapping of epistatic effects onto a genealogical tree in haploid populations

In this paper we present a model that maps epistatic effects onto a genealogical tree for a haploid population. Prior work has demonstrated that genealogical structure causes the genotypic values of individuals to covary. Our results indicate that epistasis can reduce genotypic covariance that is caused by genealogical structure. Genotypic effects (both additive and epistatic) occur along the branches of a genealogical tree, from the base of the tree to its tips. Epistasis reduces genotypic covariance because there is a reweighting of the contribution of branches to the states of genotypes compared to the additive case. Branches near the tips of a genealogical tree contribute proportionally more genetic effects with epistasis than without epistasis. Epistatic effects are most numerous at basal positions in a genealogical tree when a population is constant in size and experiencing no selection, optimizing selection, diversifying selection or directional selection, indicating that epistatic effects are typically old. For a population that is growing in size, epistatic effects are most numerous at midpoints in a genealogical tree, indicating epistatic effects are of moderate age. Our results are important in that they suggest epistatic effects may typically explain deep (old) divergences and broad patterns of divergence that exist in populations, except in growing populations. In a growing population, epistatic effects may cause more within group divergence higher up in a tree and less between group divergence that is deep in a tree. The distribution of the number of epistatic effects and the expected variance and covariance in the number of epistatic effects is also provided assuming neutrality.