Search

Until very recently, complete characterization of the megagenomes of conifers has remained elusive. The diploid genome of sugar pine (Pinus lambertiana Dougl.) has a highly repetitive, 31 billion bp genome. It is the largest genome sequenced and assembled to date, and the first from the subgenus Strobus, or white pines, a group that is notable for having the largest genomes among the pines. The genome represents a unique opportunity to investigate genome "obesity" in conifers and white pines. Comparative analysis of P. lambertiana and P. taeda L. reveals new insights on the conservation, age, and diversity of the highly abundant transposable elements, the primary factor determining genome size. Like most North American white pines, the principal pathogen of P. lambertiana is white pine blister rust (Cronartium ribicola J.C. Fischer ex Raben.). Identification of candidate genes for resistance to this pathogen is of great ecological importance. The genome sequence afforded us the opportunity to make substantial progress on locating the major dominant gene for simple resistance hypersensitive response, Cr1 We describe new markers and gene annotation that are both tightly linked to Cr1 in a mapping population, and associated with Cr1 in unrelated sugar pine individuals sampled throughout the species' range, creating a solid foundation for future mapping. This genomic variation and annotated candidate genes characterized in our study of the Cr1 region are resources for future marker-assisted breeding efforts as well as for investigations of fundamental mechanisms of invasive disease and evolutionary response.

The 1966 GENETICS papers by John Hubby and Richard Lewontin were a landmark in the study of genome-wide levels of variability. They used the technique of gel electrophoresis of enzymes and proteins to study variation in natural populations of Drosophila pseudoobscura, at a set of loci that had been chosen purely for technical convenience, without prior knowledge of their levels of variability. Together with the independent study of human populations by Harry Harris, this seminal study provided the first relatively unbiased picture of the extent of genetic variability in protein sequences within populations, revealing that many genes had surprisingly high levels of diversity. These papers stimulated a large research program that found similarly high electrophoretic variability in many different species and led to statistical tools for interpreting the data in terms of population genetics processes such as genetic drift, balancing and purifying selection, and the effects of selection on linked variants. The current use of whole-genome sequences in studies of variation is the direct descendant of this pioneering work.

Hundreds of wild-derived Drosophila melanogaster genomes have been published, but rigorous comparisons across data sets are precluded by differences in alignment methodology. The most common approach to reference-based genome assembly is a single round of alignment followed by quality filtering and variant detection. We evaluated variations and extensions of this approach and settled on an assembly strategy that utilizes two alignment programs and incorporates both substitutions and short indels to construct an updated reference for a second round of mapping prior to final variant detection. Utilizing this approach, we reassembled published D. melanogaster population genomic data sets and added unpublished genomes from several sub-Saharan populations. Most notably, we present aligned data from phase 3 of the Drosophila Population Genomics Project (DPGP3), which provides 197 genomes from a single ancestral range population of D. melanogaster (from Zambia). The large sample size, high genetic diversity, and potentially simpler demographic history of the DPGP3 sample will make this a highly valuable resource for fundamental population genetic research. The complete set of assemblies described here, termed the Drosophila Genome Nexus, presently comprises 623 consistently aligned genomes and is publicly available in multiple formats with supporting documentation and bioinformatic tools. This resource will greatly facilitate population genomic analysis in this model species by reducing the methodological differences between data sets.

Conifers are the predominant gymnosperm. The size and complexity of their genomes has presented formidable technical challenges for whole-genome shotgun sequencing and assembly. We employed novel strategies that allowed us to determine the loblolly pine (Pinus taeda) reference genome sequence, the largest genome assembled to date. Most of the sequence data were derived from whole-genome shotgun sequencing of a single megagametophyte, the haploid tissue of a single pine seed. Although that constrained the quantity of available DNA, the resulting haploid sequence data were well-suited for assembly. The haploid sequence was augmented with multiple linking long-fragment mate pair libraries from the parental diploid DNA. For the longest fragments, we used novel fosmid DiTag libraries. Sequences from the linking libraries that did not match the megagametophyte were identified and removed. Assembly of the sequence data were aided by condensing the enormous number of paired-end reads into a much smaller set of longer "super-reads," rendering subsequent assembly with an overlap-based assembly algorithm computationally feasible. To further improve the contiguity and biological utility of the genome sequence, additional scaffolding methods utilizing independent genome and transcriptome assemblies were implemented. The combination of these strategies resulted in a draft genome sequence of 20.15 billion bases, with an N50 scaffold size of 66.9 kbp.

The largest genus in the conifer family Pinaceae is Pinus, with over 100 species. The size and complexity of their genomes (∼20-40 Gb, 2n = 24) have delayed the arrival of a well-annotated reference sequence. In this study, we present the annotation of the first whole-genome shotgun assembly of loblolly pine (Pinus taeda L.), which comprises 20.1 Gb of sequence. The MAKER-P annotation pipeline combined evidence-based alignments and ab initio predictions to generate 50,172 gene models, of which 15,653 are classified as high confidence. Clustering these gene models with 13 other plant species resulted in 20,646 gene families, of which 1554 are predicted to be unique to conifers. Among the conifer gene families, 159 are composed exclusively of loblolly pine members. The gene models for loblolly pine have the highest median and mean intron lengths of 24 fully sequenced plant genomes. Conifer genomes are full of repetitive DNA, with the most significant contributions from long-terminal-repeat retrotransposons. In depth analysis of the tandem and interspersed repetitive content yielded a combined estimate of 82%.

A primary goal of evolutionary genetics is to discover and explain the genetic basis of fitness-related traits and how this genetic basis evolves within natural populations. Unprecedented technological advances have fueled the discovery of genetic variants associated with ecologically relevant phenotypes in many different life forms, as well as the ability to scan genomes for deviations from selectively neutral models of evolution. Theoretically, the degree of overlap between lists of genomic regions identified using each approach is related to the genetic architecture of fitness-related traits and the strength and type of natural selection molding variation at these traits within natural populations. Here we address for the first time in a plant the degree of overlap between these lists, using patterns of nucleotide diversity and divergence for >7000 unique amplicons described from the extensive expressed sequence tag libraries generated for loblolly pine (Pinus taeda L.) in combination with the >1000 published genetic associations. We show that loci associated with phenotypic traits are distinct with regard to neutral expectations. Phenotypes measured at the whole plant level (e.g., disease resistance) exhibit an approximately twofold increase in the proportion of adaptive nonsynonymous substitutions over the genome-wide average. As expected for polygenic traits, these signals were apparent only when loci were considered at the level of functional sets. The ramifications of this result are discussed in light of the continued efforts to dissect the genetic basis of quantitative traits.

Heterozygosity is a major challenge to efficient, high-quality genomic assembly and to the full genomic survey of polymorphism and divergence. In Drosophila melanogaster lines derived from equatorial populations are particularly resistant to inbreeding, thus imposing a major barrier to the determination and analyses of genomic variation in natural populations of this model organism. Here we present a simple genome sequencing protocol based on the whole-genome amplification of the gynogenetically derived haploid genome of a progeny of females mated to males homozygous for the recessive male sterile mutation, ms(3)K81. A single "lane" of paired-end sequences (2 × 76 bp) provides a good syntenic assembly with >95% high-quality coverage (more than five reads). The amplification of the genomic DNA moderately inflates the variation in coverage across the euchromatic portion of the genome. It also increases the frequency of chimeric clones. But the low frequency and random genomic distribution of the chimeric clones limits their impact on the final assemblies. This method provides a solid path forward for population genomic sequencing and offers applications to many other systems in which small amounts of genomic DNA have unique experimental relevance.

While many functional elements of the meiotic process are well characterized in model organisms, the genetic basis of most of the natural phenotypic variation observed in meiotic pathways has not been determined. To begin to address this issue, we characterized patterns of polymorphism and divergence in the protein-coding regions of 33 genes across 31 lines of Drosophila melanogaster and 6 lines of Drosophila simulans. We sequenced genes known to be involved in chromosome segregation, recombination, DNA repair, and related heterochromatin binding. As expected, we found several of the genes to be highly conserved, consistent with purifying selection. However, a subset of genes showed patterns of polymorphism and divergence typical of other types of natural selection. Moreover, several intriguing differences between the two Drosophila lineages were evident: along the D. simulans lineage we consistently found evidence of adaptive protein evolution, whereas along the D. melanogaster lineage several loci exhibited patterns consistent with the maintenance of protein variation.

DNA sequence surveys in yeast and humans suggest that the forces shaping telomeric polymorphism and divergence are distinctly more dynamic than those in the euchromatic, gene-rich regions of the chromosomes. However, the generality of this pattern across outbreeding, multicellular eukaryotes has not been determined. To characterize the structure and evolution of Drosophila telomeres, we collected and analyzed molecular population genetics data from the X chromosome subtelomere in 58 lines of North American Drosophila melanogaster and 29 lines of African D. melanogaster. We found that Drosophila subtelomeres exhibit high levels of both structural and substitutional polymorphism relative to linked euchromatic regions. We also observed strikingly different patterns of variation in the North American and African samples. Moreover, our analyses of the polymorphism data identify a localized hotspot of recombination in the most-distal portion of the X subtelomere. While the levels of polymorphism decline sharply and in parallel with rates of crossing over per physical length over the distal first euchromatic megabase pairs of the X chromosome, our data suggest that they rise again sharply in the subtelomeric region ([approximately equal to] 80 kbp). These patterns of historical recombination and geographic differentiation indicate that, similar to yeast and humans, Drosophila subtelomeric DNA is evolving very differently from euchromatic DNA.

Disease resistance to the bacterial pathogen Pseudomonas syringae pv. tomato (Pst) in the host species Lycopersicon esculentum, the cultivated tomato, and the closely related L. pimpinellifolium is triggered by the physical interaction between the protein products of the host resistance (R) gene Pto and the pathogen avirulence genes AvrPto and AvrPtoB. Sequence variation at the Pto locus was surveyed in natural populations of seven species of Lycopersicon to test hypotheses of host-parasite coevolution and functional adaptation of the Pto gene. Pto shows significantly higher nonsynonymous polymorphism than 14 other non-R-gene loci in the same samples of Lycopersicon species, while showing no difference in synonymous polymorphism, suggesting that the maintenance of amino acid polymorphism at this locus is mediated by pathogen selection. Also, a larger proportion of ancestral variation is maintained at Pto as compared to these non-R-gene loci. The frequency spectrum of amino acid polymorphisms known to negatively affect Pto function is skewed toward low frequency compared to amino acid polymorphisms that do not affect function or silent polymorphisms. Therefore, the evolution of Pto appears to be influenced by a mixture of both purifying and balancing selection.

We analyzed a three-locus model of genetic hitchhiking with one locus experiencing positive directional selection and two partially linked neutral loci. Following the original hitchhiking approach by Maynard Smith and Haigh, our analysis is purely deterministic, In the first half of the selected phase after a favored mutation has entered the population, hitchhiking may lead to a strong increase of linkage disequilibrium (LD) between the two neutral sites if both are 0.3s), the decay rate of LD approaches that of neutrality. Averaging over a uniform distribution of initial gamete frequencies shows that the expected LD at the end of the hitchhiking phase is driven toward zero, while the variance is increased when the selected site is well outside the two neutral sites. When the direction of LD is polarized with respect to the more common allele at each neutral site, hitchhiking creates more positive than negative linkage disequilibrium. Thus, hitchhiking may have a distinctively patterned LD-reducing effect, in particular near the target of selection.

Disease resistance to the bacterial pathogen Pseudomonas syringae pv. tomato (Pst) in the cultivated tomato, Lycopersicon esculentum, and the closely related L. pimpinellifolium is triggered by the physical interaction between plant disease resistance protein, Pto, and the pathogen avirulence protein, AvrPto. To investigate the extent to which variation in the Pto gene is responsible for naturally occurring variation in resistance to Pst, we determined the resistance phenotype of 51 accessions from seven species of Lycopersicon to isogenic strains of Pst differing in the presence of avrPto. One-third of the plants displayed resistance specifically when the pathogen expressed AvrPto, consistent with a gene-for-gene interaction. To test whether this resistance in these species was conferred specifically by the Pto gene, alleles of Pto were amplified and sequenced from 49 individuals and a subset (16) of these alleles was tested in planta using Agrobacterium-mediated transient assays. Eleven alleles conferred a hypersensitive resistance response (HR) in the presence of AvrPto, while 5 did not. Ten amino acid substitutions associated with the absence of AvrPto recognition and HR were identified, none of which had been identified in previous structure-function studies. Additionally, 3 alleles encoding putative pseudogenes of Pto were isolated from two species of Lycopersicon. Therefore, a large proportion, but not all, of the natural variation in the reaction to strains of Pst expressing AvrPto can be attributed to sequence variation in the Pto gene.

Several evolutionary models of linked selection (e.g., genetic hitchhiking, background selection, and random environment) predict a reduction in polymorphism relative to divergence in genomic regions where the rate of crossing over per physical distance is restricted. We tested this prediction near the telomere of the Drosophila melanogaster and D. simulans X chromosome at two loci, erect wing (ewg) and suppressor of sable [su(s)]. Consistent with this prediction, polymorphism is reduced at both loci, while divergence is normal. The reduction is greater at ewg, the more distal of the two regions. Two models can be discriminated by comparing the observed site frequency spectra with those predicted by the models. The hitchhiking model predicts a skew toward rare variants in a sample, while the spectra under the background-selection model are similar to those of the neutral model of molecular evolution. Statistical tests of the fit to the predictions of these models require many sampled alleles and segregating sites. Thus we used SSCP and stratified DNA sequencing to cover a large number of randomly sampled alleles (~50) from each of three populations. The result is a clear trend toward negative values of Tajima's D, indicating an excess of rare variants at ewg, the more distal of the two loci. One fixed difference among the populations and high [F.sub.ST] values indicate strong population subdivision among the three populations at ewg. These results indicate genetic hitchhiking at ewg, in particular, geographically localized hitchhiking events within Africa. The reduction of polymorphism at su(s) combined with the excess of high-frequency variants in D. simulans is inconsistent with the hitchhiking and background-selection models.

We have constructed a genetic map of Arabidopsis lyrata, a self-incompatible relative of the plant model species A. thaliana. A. lyrata is a diploid (n = 8) species that diverged from A. thaliana (n = 5) ~ 5 MYA. Mapping was conducted in a full-sib progeny of two unrelated [F.sub.1] hybrids between two European populations of A. lyrata ssp. petraea. We used the least-squares method of the Joinmap program for map construction. The gross chromosomal differences between the two species were most parsimoniously explained with three fusions, two reciprocal translocations, and one inversion. The total map length was 515 cM, and the distances were 12% larger than those between corresponding markers in the linkage map of A. thaliana. The 72 markers, consisting of microsatellites and gene-based markers, were spaced on average every 8 cM. Transmission ratio distortion was extensive, and most distortions were specific to each reciprocal cross, suggesting cytoplasmic interactions. We estimate locations and most probable genotype frequencies of transmission ratio distorting loci (TRDL) with a Bayesian method and discuss the possible reasons for the observed distortions.

We have used the naturally occurring plant-parasite system of Arabidopsis thaliana and its common parasite Peronospora parasitica (downy mildew) to study the evolution of resistance specificity in the host population. DNA sequence of the resistance gene, RPP13, from 24 accessions, including 20 from the United Kingdom, revealed amino acid sequence diversity higher than that of any protein coding gene reported so far in A. thaliana. A significant excess of amino acid polymorphism segregating within this species is localized within the leucine-rich repeat (LRR) domain of RPP13. These results indicate that single alleles of the gene have not swept through the population, but instead, a diverse collection of alleles have been maintained. Transgenic complementation experiments demonstrate functional differences among alleles in their resistance to various pathogen isolates, suggesting that the extreme amino acid polymorphism in RPP13 is maintained through continual reciprocal selection between host and pathogen.

Advances in medicine, agriculture, and an understanding of evolution depend on resolving the genetic architecture of quantitative traits, which is challenging since variation for complex traits is caused by multiple interacting quantitative trait loci (QTL) with small and conditional effects. Here, we show that the key developmental gene, hairy (h), is a QTL for Drosophila sternopleural bristle number, a model quantitative trait. Near-isoallelic lines (NIL) for the h gene region exhibited significant variation in bristle number and failed to complement a hairy mutation. Sequencing 10 h alleles from a single population revealed 330 polymorphic sites in ~10 kb. Genotypes for 25-of these and 14 additional sites in the flanking regions were determined for the 57 NIL and associated with variation in bristle number in four genetic backgrounds. A highly significant association was found for a complicated insertion/deletion polymorphism upstream of the transcription start site. This polymorphism, present in 17.5% of the h alleles, was associated with an increase of 0.5 bristle and accounted for 31% of the genetic variance in bristle number in the NIL.

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.

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[prime] 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[prime] 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 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.

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 genetic hitchhiking model predicted skewness in the case of directional selection on rare variants involving DNA sequence variation. The Tajima's D statistic gave a quantitative measure of the frequency spectrum. The statistic's expected value under the given model became large and negative with hitchhiking, indicating the model's insufficiency in explaining DNA polymorphism in restricted crossings.

Genetic experiments enable the mapping of factors in Drosophila melanogaster that influence the response to 25 generations of divergent selection from a large base population for abdominal bristle number. Derivation of isogenic chromosome substitution lines and in situ hybridization of polymorphic neutral roo transposable elements are involved in the mapping. Mapping of the factors influencing bristle number and relative viability of the chromosome recombinant isogenic lines using the multiple regression technique helps identify two factors on the X chromosome and five on the third chromosome that affect abdominal bristle number.

Sperm displacement tests on 152 homozygous lines of Drosophila melanogaster reveal that they exhibit high heterogeneity in sperm displacement in both the defense and offense directions. Analysis of sperm displacement indicates that the sperm displacement phenomenon is mediated by accessory gland protein genes. Capacity to displace resident sperm exhibits no correlation with the efficiency to resist being displaced by another sperm, indicating the involvement of distinct mechanisms in the two components of sperm displacement.

Studies covering DNA sequence polymorphism and divergence in Drosophila male accessory gland proteins are analyzed. Drosophila males transfer protein substances to females in the process of mating. The proteins exhibit influence in female reproduction. Mst26Aa and Mst26Ab transcription units are fixed in 1.6kb segments. Results show that Mst26A polymorphisms and divergence do not follow expected patterns of mutation-drift equilibrium.

Homologous recombination is the primary process involved in repairing double strand breaks. Such a repair procedure will be affected by heterozygosity because of mismatch detection within the system. Molecular evolution can therefore be influenced by polymorphism-associated recombinational repair inhibition. Mismatch detection also poses a limitation on the average heterozygosity, leads to the generation of more rare variants and diminish the rate of neutral molecular evolution.

The DdC region encompassing 65 kilobases in the second chromosome of Drosophila melanogaster was surveyed in 46 individuals from five populations of fruit flies. Thirteen polymorphic restriction sites were discovered out of the 69 scored. Heterozygosity, insertions and deletions in the region were determined. The variation in the restriction sites observed were consistent with normal distribution in neutral populations, however, insertions and deletions were significantly less frequent. Variation in dopa decarboxylase activity and linkage disequilibrium were also reported.

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.sup.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.

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 tot 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.

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.sub.0] tetrads, is observed within both natural populations. Because most nondisjunction arises from [E.sub.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.sub.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.sub.0] tetrads.

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 ootid 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.

Heterozygosity is a major challenge to efficient, high-quality genomic assembly and to the full genomic survey of polymorphism and divergence. In Drosophila melanogaster lines derived from equatorial populations are particularly resistant to inbreeding, thus imposing a major barrier to the determination and analyses of genomic variation in natural populations of this model organism. Here we present a simple genome sequencing protocol based on the whole-genome amplification of the gynogenetically derived haploid genome of a progeny of females mated to males homozygous for the recessive male sterile mutation, ms(3)K81. A single "lane" of paired-end sequences (2x76 bp) provides a good syntenic assembly with >95% high-quality coverage (more than five reads). The amplification of the genomic DNA moderately inflates the variation in coverage across the euchromatic portion of the genome. It also increases the frequency of chimeric clones. But the low frequency and random genomic distribution of the chimeric clones limits their impact on the final assemblies. This method provides a solid path forward for population genomic sequencing and offers applications to many other systems in which small amounts of genomic DNA have unique experimental relevance. [ABSTRACT FROM AUTHOR]

Examines the linkage disequilibria and site frequency spectra in the Drosophila melanogaster X chromosome regions. Reduction of DNA sequence polymorphisms; High gene conversions in regions with low crossing over; Association of hitchhiking effect and background selection models with population size reduction; Mechanism of historical recombination.