Linkage disequilibria and the site frequency spectra in the su(s) and su(w(a)) regions of the Drosophila melanogaster X chromosome.

Over the last decade, surveys of DNA sequence variation in natural populations of several Drosophila species and other taxa have established that polymorphism is reduced in genomic regions characterized by low rates of crossing over per physical length. Parallel studies have also established that divergence between species is not reduced in these same genomic regions, thus eliminating explanations that rely on a correlation between the rates of mutation and crossing over. Several theoretical models (directional hitchhiking, background selection, and random environment) have been proposed as population genetic explanations. In this study samples from an African population (n = 50) and a European population (n = 51) were surveyed at the su(s) (1955 bp) and su(w(a)) (3213 bp) loci for DNA sequence polymorphism, utilizing a stratified SSCP/DNA sequencing protocol. These loci are located near the telomere of the X chromosome, in a region of reduced crossing over per physical length, and exhibit a significant reduction in DNA sequence polymorphism. Unlike most previously surveyed, these loci reveal substantial skews toward rare site frequencies, consistent with the predictions of directional hitchhiking and random environment models and inconsistent with the general predictions of the background selection model (or neutral theory). No evidence for excess geographic differentiation at these loci is observed. Although linkage disequilibrium is observed between closely linked sites within these loci, many recombination events in the genealogy of the sampled alleles can be inferred and the genomic scale of linkage disequilibrium, measured in base pairs between sites, is the same as that observed for loci in regions of normal crossing over. We conclude that gene conversion must be high in these regions of low crossing over.

Both naturally occurring insertions of transposable elements and intermediate frequency polymorphisms at the achaete-scute complex are associated with variation in bristle number in Drosophila melanogaster.

A restriction enzyme survey of a 110-kb region including the achaete scute complex (ASC) examined 14 polymorphic molecular markers in a sample of 56 naturally occurring chromosomes. Large insertions as a class were associated with a reduction in both sternopleural and abdominal bristle number, supporting deleterious mutation-selection equilibrium models for the maintenance of quantitative genetic variation. Two polymorphic sites were independently associated with variation in bristle number measured in two genetic backgrounds as assessed by a permutation test. A 6-bp deletion near sc alpha is associated with sternopleural bristle number variation in both sexes and a 3.4-kb insertion between sc beta and sc gamma is associated with abdominal bristle number variation in females. Under an additive genetic model, the small deletion polymorphism near sc alpha accounts for 25% of the total X chromosome genetic variation in sternopleural bristle number, and the 3.4 kb insertion accounts for 22% of the total X chromosome variation in female abdominal bristle number. The observation of common polymorphisms associated with variation in bristle number is more parsimoniously explained by models that incorporate balancing selection or assume variants affecting bristle number are neutral, than mutation-selection equilibrium models.

Contrasting patterns of nucleotide polymorphism at the alcohol dehydrogenase locus in the outcrossing Arabidopsis lyrata and the selfing Arabidopsis thaliana.

Nucleotide variation at the alcohol dehydrogenase locus (Adh) was studied in the outcrossing Arabidopsis lyrata, a close relative of the selfing Arabidopsis thaliana. Overall, estimated nucleotide diversity in the North American ssp. lyrata and two European ssp. petraea populations was 0.0038, lower than the corresponding specieswide estimate for A. thaliana at the same set of nucleotide sites. The distribution of segregating sites across the gene differed between the two species. Estimated sequence diversity within an A. lyrata population with a large sample size (0.0023) was much higher than has previously been observed for A. thaliana. This North American population has an excess of sites at intermediate frequencies compared with neutral expectation (Tajima's D = 2.3, P < 0.005), suggestive of linked balancing selection or a recent population bottleneck. In contrast, an excess of rare polymorphisms has been found in A. thaliana. Polymorphism within A. lyrata and divergence from A. thaliana appear to be correlated across the Adh gene sequence. The geographic distribution of polymorphism was quite different from that of A. thaliana, for which earlier studies of several genes found low within-population nucleotide site polymorphism and no overall continental differentiation of variation despite large differences in site frequencies between local populations. Differences between the outcrossing A. lyrata and the selfing A. thaliana reflect the impact of differences in mating system and the influence of bottlenecks in A. thaliana during rapid colonization on DNA sequence polymorphism. The influence of additional variability-reducing mechanisms, such as background selection or hitchhiking, may not be discernible.

Genetic variation in rates of nondisjunction: association of two naturally occurring polymorphisms in the chromokinesin nod with increased rates of nondisjunction in Drosophila melanogaster.

Genetic variation in nondisjunction frequency among X chromosomes from two Drosophila melanogaster natural populations is examined in a sensitized assay. A high level of genetic variation is observed (a range of 0.006-0.241). Two naturally occurring variants at the nod locus, a chromokinesin required for proper achiasmate chromosome segregation, are significantly associated with an increased frequency of nondisjunction. Both of these polymorphisms are found at intermediate frequency in widely distributed natural populations. To account for these observations, we propose a general model incorporating unique opportunities for meiotic drive during female meiosis. The oötid competition model can account for both high mean rates of female-specific nondisjunction in Drosophila and humans as well as the standing genetic variation in this critical fitness character in natural populations.

Classic Weinstein: tetrad analysis, genetic variation and achiasmate segregation in Drosophila and humans.

A maximum-likelihood method for the estimation of tetrad frequencies from single-spore data is presented. The multilocus exchange with interference and viability (MEIV) model incorporates a clearly defined model of exchange, interference, and viability whose parameters define a multinomial distribution for single-spore data. Maximum-likelihood analysis of the MEIV model (MEIVLA) allows point estimation of tetrad frequencies and determination of confidence intervals. We employ MEIVLA to determine tetrad frequencies among 15 X chromosomes sampled at random from Drosophila melanogaster natural populations in Africa and North America. Significant variation in the frequency of nonexchange, or E(0) tetrads, is observed within both natural populations. Because most nondisjunction arises from E(0) tetrads, this observation is quite unexpected given both the prevalence and the deleterious consequences of nondisjunction in D. melanogaster. Use of MEIVLA is also demonstrated by reanalyzing a recently published human chromosome 21 dataset. Analysis of simulated datasets demonstrates that MEIVLA is superior to previous methods of tetrad frequency estimation and is particularly well suited to analyze samples where the E(0) tetrad frequency is low and sample sizes are small, conditions likely to be met in most samples from human populations. We discuss the implications of our analysis for determining whether an achiasmate system exists in humans to ensure the proper segregation of E(0) tetrads.

The power of association studies to detect the contribution of candidate genetic loci to variation in complex traits.

The statistical power of five association study test statistics (two haplotype-based tests, two marker-based tests, and the Transmission Disequilibrium Test-Q5) to detect single nucleotide polymorphism (SNP)/phenotype associations in a linkage-disequilibrium-based candidate gene scan employing a number of SNPs is examined. Power is estimated as a function of realistic parameters expected to affect the likelihood of detecting a significant association: the number of SNPs examined, the scaled recombination size of the region examined, the proportion of variance in the trait attributable to a hidden causative polymorphism within the region, and the number of individuals or families examined. For the different combinations of parameter values, power is estimated from a large number of realizations of a simulated coalescent describing a single random mating population with mutation, random genetic drift, and recombination. This explicit population genetics model results in a distribution of DNA marker heterozygosities and linkage disequilibria that are likely to resemble those expected in actual population samples. The study concludes that (1) marker-based permutation tests are more powerful than simple haplotype-based tests, (2) there is sufficient power to detect the presence of causative polymorphisms of small effect if on the order of 500 individuals are sampled, (3) greater power is achieved by increasing the sample size than by increasing the number of polymorphisms, (4) association studies are generally more powerful than transmission disequilibrium-based tests, and (5) for the range of parameters considered association studies have a low repeatability unless sample sizes are on the order of 500 individuals. Estimates of 4Nc for a number of gene regions and human populations will be of use in determining the density of SNPs that are likely to be required for successful association studies.

DNA polymorphism in lycopersicon and crossing-over per physical length.

Surveys in Drosophila have consistently found reduced levels of DNA sequence polymorphism in genomic regions experiencing low crossing-over per physical length, while these same regions exhibit normal amounts of interspecific divergence. Here we show that for 36 loci across the genomes of eight Lycopersicon species, naturally occurring DNA polymorphism (scaled by locus-specific divergence between species) is positively correlated with the density of crossing-over per physical length. Large between-species differences in the amount of DNA sequence polymorphism reflect breeding systems: selfing species show much less within-species polymorphism than outcrossing species. The strongest association of expected heterozygosity with crossing-over is found in species with intermediate levels of average nucleotide diversity. All of these observations appear to be in qualitative agreement with the hitchhiking effects caused by the fixation of advantageous mutations and/or "background selection" against deleterious mutations.

Two sites in the Delta gene region contribute to naturally occurring variation in bristle number in Drosophila melanogaster.

A restriction enzyme survey of a 57-kb region including the gene Delta uncovered 53 polymorphic molecular markers in a sample of 55 naturally occurring chromosomes. A permutation test, which assesses the significance of the molecular marker with the largest effect on bristle variation in four genetic backgrounds relative to permuted data-sets, found two sites that were independently associated with variation in bristle number. A common site in the second intron of Delta affected only sternopleural bristle number, and another common site in the fifth intron affected only abdominal bristle number in females. Under an additive genetic model, the polymorphism in the second intron may account for 12% of the total genetic variation in sternopleural bristle number due to third chromosomes, and the site in the fifth intron may account for 6% of the total variation in female abdominal bristle number due to the third chromosomes. These results suggest the following: (1) models that incorporate balancing selection are more consistent with observations than deleterious mutation-selection equilibrium models, (2) mapped quantitative trait loci of large effect may not represent a single variable site at a genetic locus, and (3) linkage disequilibrium can be used as a tool for understanding the molecular basis of quantitative variation.

quemao, a Drosophila bristle locus, encodes geranylgeranyl pyrophosphate synthase.

The quemao (qm) locus of Drosophila melanogaster is characterized by a P-element-associated mutant lacking most of the large bristles on the thorax and by several EMS-induced recessive lethals. quemao was cloned using a transposon tagging strategy. P-element-mediated transformation demonstrated that the cloned qm DNA sequence (from the 65F cytological region) rescues the mutant phenotype. A 2.3-kb qm transcript was identified by Northern blot analysis by sequencing of the isolated qm cDNA clones and by 5' rapid amplification cDNA end (RACE). The predicted amino acid sequence (338 residues) of the coding region of the qm transcript shares 42, 31, 13, 20, and 12% identical amino acid sequences with the geranylgeranyl pyrophosphate synthase (GGPPS) of fungi, yeast, plants, archaebacteria, and eubacteria, respectively. It also contains five highly conserved domains common among all known isoprenyl pyrophosphate synthases. The P element associated with the original qm mutant is inserted in the 5' untranslated region of the transcript. An EMS-induced qm nonsense mutation at the 12th codon leads to recessive lethality at the first larval instar, indicating the essential role of qm in the isoprenoid biosynthesis of insects.

A homologue of the 19 kDa signal recognition particle protein locus in Drosophila melanogaster.

A homologue of 19 kDa signal recognition particle locus (SRP19) was cloned and molecularly characterized in Drosophila melanogaster. It is located in the 65F region of the left arm on the third chromosome, approx. 500 bp 5' of the quemao locus. The SRP19 transcript was determined from cDNA clones, Northern blot analysis, and the 5' rapid amplification of cDNA end method. SRP19 was expressed in all the developmental stages of Drosophila. The predicted amino acid sequence (163 aa) shows that SRP19 of Drosophila shares 44%, 29%, 17% and 19% identity with the homologues from human, rice and two yeast species (Saccharomyces and Yarrowia), respectively. The most conserved amino acid residues across these species are located at those sites required for in vitro association with the 7S RNA component of the SRP.

Genetic interactions between naturally occurring alleles at quantitative trait loci and mutant alleles at candidate loci affecting bristle number in Drosophila melanogaster.

Previously, we mapped quantitative trait loci (QTL) affecting response to short-term selection for abdominal bristle number to seven suggestive regions that contain loci involved in bristle development and/or that have adult bristle number mutant phenotypes, and are thus candidates for bristle number QTL in natural populations. To test the hypothesis that the factors contributing to selection response genetically interact with these candidate loci, high and low chromosomes from selection lines were crossed to chromosomes containing wild-type or mutant alleles at the candidate loci, and the numbers of bristles were recorded in trans heterozygotes. Quantitative failure to complement, detected as a significant selection line*cross effect by analysis of variance, can be interpreted as evidence for allelism or epistasis between the factors on selected chromosomes and the candidate loci. Mutations at some candidate loci (bb, emc, h, Dl, Hairless) showed strong interactions with selected chromosomes, whereas others interacted weakly (ASC, abd, Scr) or not at all (N, mab, E(spl)). These results support the hypothesis that some candidate loci, initially identified through mutations of large effect on bristle number, either harbor or are close members in the same genetic pathway as variants that contribute to standing variation in bristle number.

The hitchhiking effect on the site frequency spectrum of DNA polymorphisms.

The level of DNA sequence variation is reduced in regions of the Drosophila melanogaster genome where the rate of crossing over per physical distance is also reduced. This observation has been interpreted as support for the simple model of genetic hitchhiking, in which directional selection on rare variants, e.g., newly arising advantageous mutants, sweeps linked neutral alleles to fixation, thus eliminating polymorphisms near the selected site. However, the frequency spectra of segregating sites of several loci from some populations exhibiting reduced levels of nucleotide diversity and reduced numbers of segregating sites did not appear different from what would be expected under a neutral equilibrium model. Specifically, a skew toward an excess of rare sites was not observed in these samples, as measured by Tajima's D. Because this skew was predicted by a simple hitchhiking model, yet it had never been expressed quantitatively and compared directly to DNA polymorphism data, this paper investigates the hitchhiking effect on the site frequency spectrum, as measured by Tajima's D and several other statistics, using a computer simulation model based on the coalescent process and recurrent hitchhiking events. The results presented here demonstrate that under the simple hitchhiking model (1) the expected value of Tajima's D is large and negative (indicating a skew toward rare variants), (2) that Tajima's test has reasonable power to detect a skew in the frequency spectrum for parameters comparable to those from actual data sets, and (3) that the Tajima's Ds observed in several data sets are very unlikely to have been the result of simple hitchhiking. Consequently, the simple hitchhiking model is not a sufficient explanation for the DNA polymorphism at those loci exhibiting a decreased number of segregating sites yet not exhibiting a skew in the frequency spectrum.

High resolution mapping of genetic factors affecting abdominal bristle number in Drosophila melanogaster.

Factors responsible for selection response for abdominal bristle number and correlated responses in sternopleural bristle number were mapped to the X and third chromosome of Drosophila melanogaster. Lines divergent for high and low abdominal bristle number were created by 25 generations of artificial selection from a large base population, with an intensity of 25 individuals of each sex selected from 100 individuals of each sex scored per generation. Isogenic chromosome substitution lines in which the high (H) X or third chromosome were placed in an isogenic low (L) background were derived from the selection lines and from the 93 recombinant isogenic (RI) HL X and 67 RI chromosome 3 lines constructed from them. Highly polymorphic neutral roo transposable elements were hybridized in situ to the polytene chromosomes of the RI lines to create a set of cytogenetic markers. These techniques yielded a dense map with an average spacing of 4 cM between informative markers. Factors affecting bristle number, and relative viability of the chromosome 3 RI lines, were mapped using a multiple regression interval mapping approach, conditioning on all markers > or = 10 cM from the tested interval. Two factors with large effects on abdominal bristle number were mapped on the X chromosome and five factors on the third chromosome. One factor with a large effect on sternopleural bristle number was mapped to the X and two were mapped to the third chromosome; all factors with sternopleural effects corresponded to those with effects on abdominal bristle number. Two of the chromosome 3 factors with large effects on abdominal bristle number were also associated with reduced viability. Significant sex-specific effects and epistatic interactions between mapped factors of the same order of magnitude as the additive effects were observed. All factors mapped to the approximate positions of likely candidate loci (ASC, bb, emc, h, mab, Dl and E(spl), previously characterized by mutations with large effects on bristle number.

Variation in sperm displacement and its association with accessory gland protein loci in Drosophila melanogaster.

Genes that influence mating and/or fertilization success may be targets for strong natural selection. If females remate frequently relative to the duration of sperm storage and rate of sperm use, sperm displacement may be an important component of male reproductive success. Although it has long been known that mutant laboratory stocks of Drosophila differ in sperm displacement, the magnitude of the naturally occurring genetic variation in this character has not been systematically quantified. Here we report the results of a screen for variation in sperm displacement among 152 lines of Drosophilia melanogaster that were made homozygous for second and/or third chromosomes recovered from natural populations. Sperm displacement was assayed by scoring the progeny of cn;bw females that had been mated sequentially to cn;bw and tested males in either order. Highly significant differences were seen in both the ability to displace sperm that is resident in the female's reproductive tract and in the ability to resist displacement by subsequent sperm. Most lines exhibited nearly complete displacement, having nearly all progeny sired by the second male, but several lines had as few as half the progeny fathered by the second male. Lines that were identified in the screen for naturally occurring variation in sperm displacement were also characterized for single-strand conformation polymorphisms (SSCP) at seven accessory gland protein (Acp) genes, Glucose dehydrogenase (Gld), and Esterase-6 (Est-6). Acp genes encode proteins that are in some cases known to be transmitted to the female in the seminal fluid and are likely candidates for genes that might mediate the phenomenon of sperm displacement. Significant associations were found between particular Acp alleles at four different loci (Acp26Aa/Ab, Acp29B, Acp36DE and Acp53E) and the ability of males to resist displacement by subsequent sperm. There was no correlation between the ability to displace resident sperm and the ability to resist being displaced by subsequent sperm. This lack of correlation, and the association of Acp alleles with resisting subsequent sperm only, suggests that different mechanisms mediate the two components of sperm displacement.

Naturally occurring variation in bristle number and DNA polymorphisms at the scabrous locus of Drosophila melanogaster.

The association between quantitative genetic variation in bristle number and molecular variation at a candidate neurogenic locus, scabrous, was examined in Drosophila melanogaster. Approximately 32 percent of the genetic variation in abdominal bristle number (21 percent for sternopleural bristle number) among 47 second chromosomes from a natural population was correlated with DNA sequence polymorphisms at this locus. Several polymorphic sites associated with large phenotypic effects occurred at intermediate frequency. Quantitative genetic variation in natural populations caused by alleles that have large effects at a few loci and that segregate at intermediate frequencies conflicts with the classical infinitesimal model of the genetic basis of quantitative variation.

Restriction map polymorphism in the forked and vermilion regions of Drosophila melanogaster.

Restriction map polymorphism at two X linked foci, forked and vermilion of Drosophila melanogaster was studied in three natural populations. The estimates of nucleotide variation were theta = 0.003 and pi = 0.002 for the forked region and theta = 0.004 and pi = 0.002 for the vermilion region. Three insertions (> 500 bp) were observed at each locus. Typical of other regions of this species each of these large insertions was unique in the sample. Non-random association among polymorphisms was common at the vermilion locus, while the forked locus was not polymorphic enough to test linkage disequilibrium. The amounts of restriction site and size variation in the vermilion and forked were within the range observed for other loci of D. melanogaster.