Nucleotide variation in the tinman and bagpipe homeobox genes of Drosophila melanogaster

The tinman (tin) and bagpipe (bap) genes are members of the NK homeobox gene family of Drosophila, so that tin occupies a higher position than bap in the regulatory hierarchy. Little is known about the level and pattern of genetic polymorphism in homeobox genes. We have analyzed nucleotide polymorphism in 27 strains of Drosophila melanogaster and one each of D. simulans and D. sechellia, within two closely linked regions encompassing a partial sequence of tin and the complete sequence of bap. The two genes exhibit different levels and patterns of nucleotide diversity. Two sets of sharply divergent sequence types are detected for tin. The haplotype structure of bap is more complex: about half of the sequences are identical (or virtually so), while the rest are fairly heterogeneous. The level of silent nucleotide variability is 0.0063 for tin but significantly higher, 0.0141, for bap, a level of polymorphism comparable to the most polymorphic structural genes of D. melanogaster. Recombination rate and gene conversion are also higher for bap than for tin. There is strong linkage disequilibrium, with the highest values in the introns of both genes and exon II of bap. The patterns of polymorphism in tin and bap are not compatible with an equilibrium model of selective neutrality. We suggest that negative selection and demographic history are the major factors shaping the pattern of nucleotide polymorphism in the tin and bap genes; moreover, there are clear indications of positive selection in the bap gene.     

Unusual pattern of nucleotide sequence variation at the OS-E and OS-F genomic regions of Drosophila simulans

Nucleotide variation at the genomic region encompassing the odorant-binding protein genes OS-E and OS-F (OS region) was surveyed in two populations of Drosophila simulans, one from Europe and the other from Africa. We found that the European population shows an atypical and large haplotype structure, which extends throughout the approximately 5-kb surveyed genomic region. This structure is depicted by two major haplotype groups segregating at intermediate frequency in the sample, one haplogroup with nearly no variation, and the other at levels more typical for this species. This pattern of variation was incompatible with neutral predictions for a population at a stationary equilibrium. Nevertheless, neutrality tests contrasting polymorphism and divergence data fail to detect any departure from the standard neutral model in this species, whereas they confirm the non-neutral behavior previously observed at the OS-E gene in D. melanogaster. Although positive Darwinian selection may have been responsible for the observed unusual nucleotide variation structure, coalescent simulation results do not allow rejecting the hypothesis that the pattern was generated by a recent bottleneck in the history of European populations of D. simulans.     

Molecular variation at the In(2L)t proximal breakpoint site in natural populations of Drosophila melanogaster and D. simulans

A previous study of nucleotide polymorphism in a Costa Rican population of Drosophila melanogaster found evidence for a nonneutral deficiency in the number of haplotypes near the proximal breakpoint of In(2L)t, a common inversion polymorphism in this species. Another striking feature of the data was a window of unusually high nucleotide diversity spanning the breakpoint site. To distinguish between selective and neutral demographic explanations for the observed patterns in the data, we sample alleles from three additional populations of D. melanogaster and one population of D. simulans. We find that the strength of associations among sites found at the breakpoint varies between populations of D. melanogaster. In D. simulans, analysis of the homologous region reveals unusually elevated levels of nucleotide polymorphism spanning the breakpoint site. As with American populations of D. melanogaster, our D. simulans sample shows a marked reduction in the number of haplotypes but not in nucleotide diversity. Haplotype tests reveal a significant deficiency in the number of haplotypes relative to the neutral expectation in the D. simulans sample and some populations of D. melanogaster. At the breakpoint site, the level of divergence between haplotype classes is comparable to interspecific divergence. The observation of interspecific polymorphisms that differentiate major haplotype classes in both species suggests that haplotype classes at this locus are considerably old. When considered in the context of other studies on patterns of variation within and between populations of D. melanogaster and D. simulans, our data appear more consistent with the operation of selection than with simple demographic explanations.     

Actin gene family of Caenorhabditis elegans

Four actin genes have been isolated from Caenorhabditis elegans that account for all of the major actin hybridization to total genomic DNA. Actin genes I, II and III are clustered within a 12 X 10(3) base region; gene IV is unlinked to the others. All four genes have been sequenced from at least nucleotide -109 to +250. Genes I and III are identical for the first 307 coding nucleotides. Genes I and II differ in 14 positions within the first 250 coding nucleotides; one difference substitutes an aspartic acid for a glutamic acid at codon 5. Genes I and IV differ in 18 positions within the first 259 coding nucleotides without causing any amino acid differences. Genes I, II and III have introns after the first nucleotide of codon 64 and gene IV has an intron between codons 19 and 20. The four nucleotide sequences thus far define two different amino acid sequences. Both of the amino acid sequences resemble vertebrate cytoplasmic actin more than vertebrate muscle actin. A DNA polymorphism between the Bristol and Bergerac strains has been used as a phenotypic marker in genetic crosses to map the cluster of actin genes within a 2% recombination interval on linkage group V between unc-23 and sma-1 in order to begin a molecular genetic analysis of the actin loci.     

SNP Haplotype Mapping in a Small ALS Family

The identification of genes for monogenic disorders has proven to be highly effective for understanding disease mechanisms, pathways and gene function in humans. Nevertheless, while thousands of Mendelian disorders have not yet been mapped there has been a trend away from studying single-gene disorders. In part, this is due to the fact that many of the remaining single-gene families are not large enough to map the disease locus to a single site in the genome. New tools and approaches are needed to allow researchers to effectively tap into this genetic gold-mine. Towards this goal, we have used haploid cell lines to experimentally validate the use of high-density single nucleotide polymorphism (SNP) arrays to define genome-wide haplotypes and candidate regions, using a small amyotrophic lateral sclerosis (ALS) family as a prototype. Specifically, we used haploid-cell lines to determine if high-density SNP arrays accurately predict haplotypes across entire chromosomes and show that haplotype information significantly enhances the genetic information in small families. Panels of haploid-cell lines were generated and a 5 centimorgan (cM) short tandem repeat polymorphism (STRP) genome scan was performed. Experimentally derived haplotypes for entire chromosomes were used to directly identify regions of the genome identical-by-descent in 5 affected individuals. Comparisons between experimentally determined and in silico haplotypes predicted from SNP arrays demonstrate that SNP analysis of diploid DNA accurately predicted chromosomal haplotypes. These methods precisely identified 12 candidate intervals, which are shared by all 5 affected individuals. Our study illustrates how genetic information can be maximized using readily available tools as a first step in mapping single-gene disorders in small families.     

Molecular characterization of major histocompatibility complex (B) haplotypes in broiler chickens

In Leghorn (laying) chickens, susceptibility to a number of infectious diseases is strongly associated with the major histocompatibility ( B ) complex. Nucleotide sequence data have been published for six class I ( B-F ) alleles and for class II ( B-Lβ ) alleles or isotypes from 17 Leghorn haplotypes. It is not known if classical B-L or B-F alleles in broilers are identical, at the sequence level, to any Leghorn alleles. This report describes molecular and immunogenetic characterization of two haplotypes from commercial broiler breeder chickens that were originally identified by serology as a single haplotype, but were differentiated serologically in the present work. The two haplotypes, designated B ^A4 and B ^A4variant, shared identical B-G restriction fragment length polymorphism patterns, but differed in one B-Lβ fragment that cosegregated with the serological B haplotype. Furthermore, the nucleotide sequences of the highly variable exons of an expressed B-LβII family gene and B-F gene from the two haplotypes were markedly different from each other. Both the B-LβII family and B-F gene sequences from the B ^A4 haplotype were identical to the sequences obtained from the reference B ²¹ haplotype in Leghorns; however, in the B ^A4 haplotype the B-Lβ ²¹ and B-F ²¹ alleles were in linkage with B-G alleles that were not G ²¹. The nucleotide sequences from B ^A4variant were unique among the reported chicken B-LβII family and B-F alleles.     

Identification by denaturing high-performance liquid chromatography of numerous polymorphisms in a candidate region for multiple sclerosis susceptibility

Genetic association analysis of candidate regions where evidence of linkage has accumulated is becoming a key issue in the study of complex diseases. A high density of markers, at least one per centimorgan, is required to improve the chances of observing linkage disequilibrium with disease alleles. A recently available single nucleotide polymorphism (SNP) map designed to cover the whole genome provides an average density of one marker per 2 cM. In the present study we show that the number of markers can be approximately doubled in a selected region, thus reaching a density suitable for association studies, by applying a completely automated technique for polymorphism detection, denaturing high-performance liquid chromatography (DHPLC). A systematic search for SNPs was performed in the region 5ptel-q13, where weak but convergent evidence for linkage with multiple sclerosis has accumulated. Screening for polymorphisms was performed on 124 sequence tagged sites (STSs) in the 3'UTR ends of expressed sequence tags totaling about 30,000 bp. Thirty SNPs in 28 STSs were found with less than 10% overlap with the markers already detected in the same region. The data confirm the validity of the approach using DHPLC on expressed gene sequences tagged by a set of standard commercially available primers.     

Major-histocompatibility-complex DRB genes of a New-World monkey, the cottontop tamarin (Saguinus oedipus).

The DRB region of the human and great-ape major histocompatibility complex displays not only gene but also haplotype polymorphism. The number of genes in the human DRB region can vary from one to four, and even greater variability exists among the DRB haplotypes of chimpanzees, gorillas, and orangutans. Accumulating evidence indicates that, like gene polymorphism, part of the haplotype polymorphism predates speciation. In an effort to determine when the gene haplotype polymorphisms emerged in the primate lineage, we sequenced three cDNA clones of the New-World monkey, the cottontop tamarin (Saguinus oedipus). We could identify two DRB loci in this species, one (Saoe-DRB1) occupied by apparently functional alleles (*0101 and *0102) which differ by only two nucleotide substitutions and the other (Saoe-DRB2) occupied by an apparent pseudogene. The Saoe-DRB2 gene contains an extra sequence derived from the 3' portion of exon 2 and placed 5' to this exon. This sequence contains a stop codon which makes the translation of the bulk of the Saoe-DRB2 gene unlikely. Preliminary Southern blot hybridization analysis with probes derived from these two genes suggests that both the DRB gene polymorphism and the haplotype polymorphism in the cottontop tamarin may be low. In most individuals the DRB region of this species probably consists of three genes. Comparisons of the Saoe-DRB sequences with those of other primates suggest that probably all of the DRB genes found until now in the Catarrhini were derived from a common ancestor after the separation of the Catarrhini and Platyrrhini lineages. The extant DRB gene and haplotype polymorphism may therefore have been founded in the mid-Oligocene some 33 Mya.     

Recurrent deletion and gene presence/absence polymorphism: telomere dynamics dominate evolution at the tip of 3L in Drosophila melanogaster and D. simulans

Although Drosophila melanogaster has been the subject of intensive analysis of polymorphism and divergence, little is known about the distribution of variation at the most distal regions of chromosomes arms. Here we report a survey of genetic variation on the tip of 3L in D. melanogaster and D. simulans. Levels of single nucleotide polymorphism in the most distal euchromatic sequence are approximately one order of magnitude less than that typically observed in genomic regions of normal crossing over, consistent with what might be expected under models of linked selection in regions of low crossing over. However, despite this reduced level of nucleotide variation, we found abundant deletion polymorphism. These deletions create at least three gene presence/absence polymorphisms within D. melanogaster: the putative G-protein coupled receptor mthl-8 (which is the most distal known or predicted gene on 3L) and the unannotated mRNAs AY060886 and BT006009. Strikingly, D. simulans is also segregating deletions that cause mthl8 presence/absence polymorphism. Breakpoint sequencing and tests of correlations with segregating SNPs in D. melanogaster suggest that each deletion is unique. Cloned breakpoint sequences revealed the presence of Het-A elements just distal to unique, canonical euchromatic sequences. This pattern suggests a model in which repeated telomeric deficiencies cause deletions of euchromatic sequence followed by subsequent "healing" by retrotranposition of Het-A elements. These data reveal the dominance of telomeric dynamics on the evolution of closely linked sequences in Drosophila.     

Molecular evolution of two linked genes, Est-6 and Sod, in Drosophila melanogaster.

We have obtained 15 sequences of Est-6 from a natural population of Drosophila melanogaster to test whether linkage disequilibrium exists between Est-6 and the closely linked Sod, and whether natural selection may be involved. An early experiment with allozymes had shown linkage disequilibrium between these two loci, while none was detected between other gene pairs. The Sod sequences for the same 15 haplotypes were obtained previously. The two genes exhibit similar levels of nucleotide polymorphism, but the patterns are different. In Est-6, there are nine amino acid replacement polymorphisms, one of which accounts for the S-F allozyme polymorphism. In Sod, there is only one replacement polymorphism, which corresponds to the S-F allozyme polymorphism. The transversion/transition ratio is more than five times larger in Sod than in Est-6. At the nucleotide level, the S and F alleles of Est-6 make up two allele families that are quite different from each other, while there is relatively little variation within each of them. There are also two families of alleles in Sod, one consisting of a subset of F alleles, and the other consisting of another subset of F alleles, designed F(A), plus all the S alleles. The Sod F(A) and S alleles are completely or nearly identical in nucleotide sequence, except for the replacement mutation that accounts for the allozyme difference. The two allele families have independent evolutionary histories in the two genes. There are traces of statistically significant linkage disequilibrium between the two genes that, we suggest, may have arisen as a consequence of selection favoring one particular sequence at each locus.     

Constraints on Intron Evolution in the Gene Encoding the Myosin Alkali Light Chain in Drosophila

Interspecific comparisons of intron sequences reveal conserved blocks of invariant nucleotides and several other departures from the strictly neutral model of molecular evolution. To distinguish the past action of evolutionary forces in introns known to have regulatory information, we examined nucleotide sequence variation at 991 sites in a random sample of 16 Drosophila melanogaster alleles of the gene encoding the myosin alkali light chain (Mlc1). The Mlc1 gene of D. melanogaster encodes two MLC1 isoforms via developmentally regulated alternative pre-mRNA splicing. Analyses of these data reveal that introns 4 and 5, which flank the alternatively spliced exon 5, have reduced levels of both intraspecific polymorphism and interspecific divergence relative to intron 3. No polymorphism was observed in any of the exons examined in D. melanogaster. A genealogical analysis clearly demonstrates the occurrence of intragenic recombination in the ancestral history of Mlc1. Recombination events are estimated to be 13 times more likely than mutation events over the span of the sequenced region. Although there is little evidence for pairwise linkage disequilibrium in the Mlc1 region, higher order disequilibrium does seem to be present in the 5' half of the portion of the gene that was examined. Predictions of the folding free energy of the pre-mRNA reveal that sampled alleles have a significantly higher (less stable) free energy than do randomly permuted sequences. These results are consistent with the hypothesis that introns surrounding an alternatively spliced exon are subjected to additional constraints, perhaps due to specific aspects of secondary structure required for appropriate splicing of the pre-mRNA molecule.     

Complex signatures of selection and gene conversion in the duplicated globin genes of house mice

Results of electrophoretic surveys have suggested that hemoglobin polymorphism may be maintained by balancing selection in natural populations of house mice, Mus musculus. Here we report a survey of nucleotide variation in the adult globin genes of house mice from South America. We surveyed nucleotide polymorphism in two closely linked alpha-globin paralogs and two closely linked beta-globin paralogs to test whether patterns of variation are consistent with a model of long-term balancing selection. Surprisingly high levels of nucleotide polymorphism at the two beta-globin paralogs were attributable to the segregation of two highly divergent haplotypes, Hbbs (which carries two identical beta-globin paralogs) and Hbbd (which carries two functionally divergent beta-globin paralogs). Interparalog gene conversion on the Hbbs haplotype has produced a highly unusual situation in which the two paralogs are more similar to one another than either one is to its allelic counterpart on the Hbbd haplotype. Levels of nucleotide polymorphism and linkage disequilibrium at the two beta-globin paralogs suggest a complex history of diversity-enhancing selection that may be responsible for long-term maintenance of alternative protein alleles. The alternative two-locus beta-globin haplotypes are associated with pronounced differences in intraerythrocyte glutathione and nitric oxide metabolism, suggesting a possible mechanism for selection on hemoglobin function.     

An efficient haplotyping method with DNA pools

Determination of haplotype frequencies (the joint distribution of genetic markers) in large population samples is a powerful tool for association studies. This is due to their greater extent of polymorphism since any two bi-allelic single nucleotide polymorphisms (SNPs) generate a potential four-allele genetic marker. Therefore, a haplotype may capture a given functional polymorphism with higher statistical power than its SNP components. The statistical estimation of haplotype frequencies, usually employed in linkage disequilibrium studies, requires individual genotyping for each SNP in the haplotype, thus making it an expensive process. In this study, we describe a new method for direct measurement of haplotype frequencies in DNA pools by allele-specific, long-range haplotype amplification. The proposed method allows the efficient determination of haplotypes composed of two SNPs in close vicinity (up to 20 kb).     

Positive selection driving the evolution of a gene of male reproduction, Acp26Aa, of Drosophila: II. Divergence versus polymorphism

The evolution of the gene for a male ejaculatory protein, Acp26Aa, has been shown to be driven by positive selection when nonsibling species in the Drosophila melanogaster subgroup are compared. To know if selection has been operating in the recent past and to understand the details of its dynamics, we obtained DNA sequences of Acp26Aa and the nearby Acp26Ab gene from 39 D. melanogaster chromosomes. Together with the 10 published sequences, we analyzed 49 sequences from five populations in four continents. The southern African population is somewhat differentiated from all other populations, but its nucleotide diversity is lower at these two loci. We find the following results for Acp26Aa: (1) The R: S (replacement : silent changes) ratio is significantly higher in the between-species comparisons than in the within-species data by the McDonald and Kreitman test. Positive selection is probably responsible for the excess of amino acid replacements between species. (2) However, within-species nucleotide diversity is high. Neither the Tajima test nor the Fu and Li test indicates a reduction in nucleotide diversity due to positive selection in the recent past. (3) The newly derived nucleotides in D. melanogaster are at high frequency significantly more often than predicted by the neutral equilibrium. Since the nearby Acp26Ab gene does not show these patterns, these observations cannot be attributed to the characteristics of this chromosomal region. We suggest that positive selection is active, but may be weak, for each amino acid change in the Acp26Aa gene.      

Different instability of the CAG microsatellite in two haplotype groups of human mitochondrial DNA polymerase gamma

Two single nucleotide polymorphisms of the mitochondrial DNA polymerase gamma gene (POLG1), rs2238296 (T/C) and rs758130 (T/C), were analyzed in individuals of different ethnicity (Russians and Buryats) with known genotypes of the CAG microsatellite located in the same gene. It was shown that microsatellite alleles with repeat numbers other than 10 were significantly more frequent within the TT haplotype. A phylogenetic analysis of human and chimpanzee POLG1 intron 2 sequences suggested that the haplotype TT, which is more heterogeneous regarding the CAG repeat polymorphism, is evolutionally younger than the haplotype CC. These data may be useful in the further research of the association between the CAG microsatellite polymorphism of POLG1 and male infertility.     

Molecular characterization and evolutionary pattern of the 9- cis -epoxycarotenoid dioxygenase NCED1 gene in grapevine

This study provides the first analysis of the level and patterns of nucleotide polymorphism of the NCED1 gene in grapevine (Vitis vinifera L.). A total of 123 sequences of the gene were analyzed to give a sample of 50 wild accessions and 73 cultivars. A high single nucleotide polymorphism and haplotype diversity was revealed in the cultivars studied, especially Tunisian germplasms which present an important and diverse reservoir of genetic diversity for grape breeding and conservation. The haplotype distribution highlights two origins of the cultivars studied: one may be related to primary grapevine gene pool domestication while the second seems to be more recent. Thus, besides domestication, gene introgression has also played a role in shaping the current varietal landscape of grape cultivars. Higher nucleotide and haplotype polymorphism was recorded for cultivars. This was accompanied by a higher recombination rate in cultivated grapevines for this gene, a recent selective sweep in wild samples and a balancing selection in cultivars. The conservation of genetic diversity of the endangered wild germplasm is important to ensure that the wild population can be used in future breeding programs of the domesticated cultivars. The high number of alleles discovered can be used as a valuable source for association studies between allele frequencies and phenotypic variations in this gene. In addition to natural selection, molecular evidence shows that genetic variation in this locus appears to be shaped by a combination of mutation and recombination events.     

Inferences on the evolutionary history of the S-element family of Drosophila melanogaster

The S-element family of transposable elements has been characterized in D. melanogaster. Attempts to find it in other Drosophila-related species have failed, suggesting that this element family may have recently invaded the D. melanogaster genome by horizontal transfer. In order to investigate its evolutionary history, we analyzed the patterns of DNA polymorphism among the S-element copies present in a sample genome (Drosophila Genome Project). The observed levels of nucleotide diversity are significantly lower than theoretical expectations based on the neutral model. This is consistent with evidence for ongoing gene conversion among copies and for purifying selection on the elements' sequences, particularly on the terminal inverted repeats. A phylogenetic analysis revealed that the members of the S-element family can be grouped into at least two genetically differentiated clusters. The level of divergence between these clusters suggests that the S elements invaded the genome of the ancestor of D. melanogaster before the speciation of the D. melanogaster complex. However, other relevant scenarios are also discussed.     

Extensive amino acid polymorphism at the pgm locus is consistent with adaptive protein evolution in Drosophila melanogaster

PGM plays a central role in the glycolytic pathway at the branch point leading to glycogen metabolism and is highly polymorphic in allozyme studies of many species. We have characterized the nucleotide diversity across the Pgm gene in Drosophila melanogaster and D. simulans to investigate the role that protein polymorphism plays at this crucial metabolic branch point shared with several other enzymes. Although D. melanogaster and D. simulans share common allozyme mobility alleles, we find these allozymes are the result of many different amino acid changes at the nucleotide level. In addition, specific allozyme classes within species contain several amino acid changes, which may explain the absence of latitudinal clines for PGM allozyme alleles, the lack of association of PGM allozymes with the cosmopolitan In(3L)P inversion, and the failure to detect differences between PGM allozymes in functional studies. We find a significant excess of amino acid polymorphisms within D. melanogaster when compared to the complete absence of fixed replacements with D. simulans. There is also strong linkage disequilibrium across the 2354 bp of the Pgm locus, which may be explained by a specific amino acid haplotype that is high in frequency yet contains an excess of singleton polymorphisms. Like G6pd, Pgm shows strong evidence for a branch point enzyme that exhibits adaptive protein evolution.