Inferences on the evolutionary history of the Drosophila americana polymorphic X/4 fusion from patterns of polymorphism at the X-linked paralytic and elav genes

In Drosophila there is limited evidence on the nature of evolutionary forces affecting chromosomal arrangements other than inversions. The study of the X/4 fusion polymorphism of Drosophila americana is thus of interest. Polymorphism patterns at the paralytic (para) gene, located at the base of the X chromosome, suggest that there is suppressed crossing over in this region between fusion and nonfusion chromosomes but not within fusion and nonfusion chromosomes. These data are thus compatible with previous claims that within fusion chromosomes the amino acid clines found at fused1 (also located at the base of the X chromosome) are likely maintained by local selection. The para data set also suggests a young age of the X/4 fusion. Polymorphism data on para and elav (located at the middle region of the X chromosome) suggest that there is no population structure other than that caused by the X/4 fusion itself. These findings are therefore compatible with previous claims that selection maintains the strong association observed between the methionine/threonine variants at fused1 and the status of the X chromosome as fused or unfused to the fourth chromosome.     

Clinal variation for amino acid polymorphisms at the Pgm locus in Drosophila melanogaster

Clinal variation is common for enzymes in the glycolytic pathway for Drosophila melanogaster and is generally accepted as an adaptive response to different climates. Although the enzyme phosphoglucomutase (PGM) possesses several allozyme polymorphisms, it is unique in that it had been reported to show no clinal variation. Our recent DNA sequence investigation of Pgm found extensive cryptic amino acid polymorphism segregating with the allozyme alleles. In this study, we characterize the geographic variation of Pgm amino acid polymorphisms at the nucleotide level along a latitudinal cline in the eastern United States. A survey of 15 SNPs across the Pgm gene finds significant clinal differentiation for the allozyme polymorphisms as well as for many of the cryptic amino acid polymorphisms. A test of independence shows that pervasive linkage disequilibrium across this gene region can explain many of the amino acid clines. A single Pgm haplotype defined by two amino acid polymorphisms shows the strongest correlation with latitude and the steepest change in allele frequency across the cline. We propose that clinal selection at Pgm may in part explain the extensive amino acid polymorphism at this locus and is consistent with a multilocus response to selection in the glycolytic pathway.     

Evidence for selection at the fused locus of Drosophila virilis

The genomic DNA sequence of a 2.4-kb region of the X-linked developmental gene fused was determined in 15 Drosophila virilis strains. One common replacement polymorphism is observed, where a negatively charged aspartic amino acid is replaced by the noncharged amino acid alanine. This replacement variant is located within the serine/threonine kinase domain of the fused gene and is present in approximately 50% of the sequences in our sample. Significant linkage disequilibrium is detected around this replacement site, although the fused gene is located in a region of the D. virilis X chromosome that seems to experience normal levels of recombination. In a 600-bp region around the replacement site, all eight alanine sequences are identical; of the six aspartic acid sequences, three are also identical. The occurrence of little or no variation within the aspartic acid and alanine haplotypes, coupled with the presence of several differences between them, is very unlikely under the usual equilibrium neutral model. Our results suggest that the fused alanine haplotypes have recently increased in frequency in the D. virilis population.     

On the location of the gene(s) harbouring the advantageous variant that maintains the X/4 fusion of Drosophila americana

Weak selection is maintaining the Drosophila americana X/4 fusion chromosomal frequency cline. The gene(s) harbouring the advantageous variant(s) that is responsible for the establishment and maintenance of this chromosomal frequency gradient must be located in a region of the X and/or 4th chromosome that is genetically isolated between the X/4 fusion and non-fusion forms. The limits of these regions must thus be determined before an attempt is made to identify these genes. For this purpose, the correspondence between the D. virilis X and 4th chromosome genome scaffolds sequence and the D. americana gene order was established. Polymorphism levels and patterns at seven genes located at the base of the D. americana X chromosome, as well as three genes located at the base of the 4th chromosome, were analysed. The data suggest that the D. americana X/4 fusion is no more than 29,000 years old. At the base of the X chromosome, there is suppression of recombination within X/4 fusion and non-fusion chromosomes, and little recombination between the two chromosomal forms. Apparent fixed silent and replacement differences are found in three of seven genes analysed located at the base of the X chromosome. There is no evidence for suppression of recombination between fusion and non-fusion chromosomes at the base of the 4th chromosome. The advantageous variant responsible for the establishment in frequency and maintenance of the X/4 fusion is thus inferred to be in the D. americana X centromere-inversion Xc basal breakpoint region.     

Haplotypic Divergence Coupled with Lack of Diversity at the Arabidopsis Thaliana Alcohol Dehydrogenase Locus: Roles for Both Balancing and Directional Selection?

We designate a region of the alcohol dehydrogenase locus (Adh) of the weedy crucifer, Arabidopsis thaliana, as ``hypervariable'''' on the basis of a comparison of sequences from ecotypes Columbia and Landsberg. We found eight synonymous and two replacement mutations in the first 262 nucleotides of exon 4, and an additional two mutations in the contiguous region of intron 3. The rest of the sequence (2611 bp) has just three mutations, all of them confined to noncoding regions. Our survey of the hypervariable region among 37 ecotypes of A. thaliana revealed two predominant haplotypes, corresponding to the Columbia and Landsberg sequences. We identified five additional haplotypes and 4 additional segregating sites. The lack of haplotype diversity is presumably in part a function of low rates of recombination between haplotypes conferred by A. thaliana''s tendency to self-fertilize. However, an analysis in 32 ecotypes of 12 genome-wide polymorphic markers distinguishing Columbia and Landsberg ecotypes indicated levels of outcrossing sufficient at least to erode linkage disequilibrium between dispersed markers. We discuss possible evolutionary explanations for the coupled observation of marked divergence within the hypervariable region and a lack of haplotype diversity among ecotypes. The sequence of the region for closely related species argues against the possibility that one allele is the product of introgression. We note (1) that several loss of function mutations (both naturally and chemically induced) map to the hypervariable region, and (2) the presence of two amino acid replacement polymorphisms, one of which causes the mobility difference between the two major classes of A. thaliana Adh electrophoretic alleles. We argue that protein polymorphism in such a functionally significant part of the molecule may be subject to balancing selection. The observed pattern of extensive divergence between the alleles is consistent with this explanation because balancing selection on a particular site maintains linked neutral polymorphisms at intermediate frequencies.     

Polymorphic DNA haplotypes at the phenylalanine hydroxylase (PAH) locus in Asian families with phenylketonuria (PKU)

DNA polymorphisms at the phenylalanine hydroxylase (PAH) locus have proved highly effective in linkage diagnosis of phenylketonuria (PKU) in Caucasian families. More than 10 RFLP sites have been reported within the PAH structural locus in Caucasians. With information from affected and unaffected offspring in PKU families it is often possible to reconstruct complete RFLP haplotypes in parents and to use these haplotypes to follow the segregation of PKU within families and to determine the distribution of PKU chromosomes within populations. To establish the utility of these RFLPs in characterizing Asian families with PKU, we typed eight DNA sites in 21 Chinese families and 12 Japanese families with classical PKU. The eight RFLPs were chosen for their informativeness in Caucasians. From these families we reconstructed a total of 91 complete PAH haplotypes, 44 from non-PKU chromosomes and 47 from PKU-bearing chromosomes. Although all eight marker sites are polymorphic in both Chinese and Japanese, there is much less haplotypic variation in Asians than in Caucasians. In particular, one haplotype alone, haplotype 4, accounts for more than 77% of non-PKU chromosomes and for more than 80% of PKU-bearing chromosomes. Haplotype 4 is also relatively common in Caucasians. The next most common Asian haplotype is 10 times less frequent than haplotype 4. By contrast, in many Caucasian populations the sum of the frequencies of the five most common haplotypes is still less than 80%, and several of the most common haplotypes are equally frequent. Even though the extent of haplotypic variation in Asians is severely limited, the few haplotypes that are found often differ at a number of RFLP sites.(ABSTRACT TRUNCATED AT 250 WORDS)     

Molecular basis and haplotyping of the alphaII domain polymorphisms of spectrin: application to the study of hereditary elliptocytosis and pyropoikilocytosis.

Hereditary elliptocytosis (HE) and hereditary pyropoikilocytosis (HPP) are inherited disorders of erythrocyte shape that are frequently associated with abnormalities in alpha-spectrin, one of the principal structural proteins of the erythrocyte membrane skeleton. Five polymorphisms of the alpha-spectrin gene, located in a 6-kb interval of genomic DNA, were identified and analyzed in normal and mutant alpha-spectrin alleles. Three of these polymorphisms are due to single nucleotide substitutions in the alpha-spectrin gene coding region that lead to changes in the amino acid sequence. In combination, these three polymorphisms are responsible for the different peptide phenotypes of the alphaII domain previously observed following limited tryptic digestion of spectrin protein. The most common haplotype, type 1, was found predominantly in Caucasians and was the only haplotype identified in Asians. Haplotypes 2, 3, and 4 were identified predominantly in individuals of African ancestry and were commonly found in patients with HE or HPP. Analysis of coinheritance of alphaII domain polymorphisms with alpha-spectrin gene mutations causing HE or HPP in African-American patients with HE and HPP suggests that, with one exception, a given HE/HPP mutation is present in an alpha-spectrin gene of only one haplotype, indicating a founder effect. The other two polymorphisms located in this region of the alpha-spectrin gene do not change the amino acid sequence of the encoded alpha-spectrin chain and are not in linkage disequilibrium with three of the four alphaII domain haplotypes. A model is proposed for the evolutionary origin of the different haplotypes.     

Molecular Analysis of an Allozyme Cline: Alcohol Dehydrogenase in Drosophila Melanogaster on the East Coast of North America

Clines may either be selectively maintained or be the by-product of nonadaptive processes related to population structure and history. Drosophila melanogaster populations on the east coast of North America show a latitudinal cline in the frequencies of two common electrophoretically distinguishable alleles at the alcohol dehydrogenase locus (Adh), designated Adh-S and Adh-F. This cline may either be adaptive or an artifact of a possible recent dual founding of North American D. melanogaster populations in which frequencies of Adh alleles differed between founder populations. By means of a high resolution restriction-mapping technique, we studied the distribution of 113 haplotypes derived from 44 polymorphic DNA markers within the Adh region in 1533 individuals from 25 populations throughout the cline. We found significant clinal differentiation at the polymorphisms determining the mobility-difference causing amino acid replacement between Adh-F and Adh-S alleles. Hitchhiking was limited, despite extensive linkage disequilibrium, and other sites did not vary clinally. Such a pattern of differentiation implies that selection is responsible for the cline. To investigate whether selection acts only on the Adh-F/S site, we performed a ``selective equivalence'''' test under the assumption that all variability within the specified allelic class is selectively neutral. This revealed selective equivalence among Adh-S-bearing haplotypes, whose frequencies showed no differentiation throughout the cline, implying high levels of frequency-homogenizing gene flow. Geographical heterogeneity among Adh-F-bearing haplotypes implied the action of selection on one or more additional variants in linkage disequilibrium with Adh-F. In a further study of a subset of the data (n = 1076 from 18 populations), we found a combined insertion/deletion polymorphism, designated &1, located in the 5'' adult intron and in linkage disequilibrium with Adh-F, to show more marked clinal variation than Adh-F/S. Although the unequivocal identification of the precise target(s) of selection requires further study, we suggest that clinal selection may be acting epistatically on the Adh-F/S and &1 polymorphisms.     

Point mutations with positive selection were a major force during the evolution of a receptor-kinase resistance gene family of rice

The rice (Oryza sativa) Xa26 gene, which confers resistance to bacterial blight disease and encodes a leucine-rich repeat (LRR) receptor kinase, resides at a locus clustered with tandem homologous genes. To investigate the evolution of this family, four haplotypes from the two subspecies of rice, indica and japonica, were analyzed. Comparative sequence analysis of 34 genes of 10 types of paralogs of the family revealed haplotype polymorphisms and pronounced paralog diversity. The orthologs in different haplotypes were more similar than the paralogs in the same haplotype. At least five types of paralogs were formed before the separation of indica and japonica subspecies. Only 7% of amino acid sites were detected to be under positive selection, which occurred in the extracytoplasmic domain. Approximately 74% of the positively selected sites were solvent-exposed amino acid residues of the LRR domain that have been proposed to be involved in pathogen recognition, and 73% of the hypervariable sites detected in the LRR domain were subject to positive selection. The family is formed by tandem duplication followed by diversification through recombination, deletion, and point mutation. Most variation among genes in the family is caused by point mutations and positive selection.     

Emergence of Complex Haplotypes from Microevolutionary Variation in Sequence and Structure of <Emphasis Type="Italic">Colias</Emphasis> Phosphoglucose Isomerase

A molecular evolutionary explanation of natural genetic variation requires analysis of specific variants’ evolutionary dynamics. To pursue this for phosphoglucose isomerase (PGI) of Colias butterflies, whose polymorphism is maintained by strong natural selection, we assembled a large data set of wild haplotypes, highly variable at the amino acid and DNA levels. The most common electrophoretic, i.e., charge macrostate, allele class, 3, is conserved in its pattern of charged amino acid residues. The next most common macrostate, 4, has multiple patterns of charge, i.e., microstates, while less common (1, 2, 5, 6) macrostates are very diverse. Macrostate 4 shows significant linkage disequilibrium (LD) among its variants, especially for two groups of five haplotypes each. We find extensive intragenic recombination among all haplotypes except the two high-LD groups of macrostate 4, which display none. Phyletic relations among haplotypes are largely reticulate, again except for the high-LD groups of macrostate 4, which form clades with strong bootstrap support. Charge-changing and linked charge-neutral amino acid variants occur in diverse parts of PGI’s sequence. Homology-based modeling of PGI’s structure shows that these regions are related spatially in ways suggesting functional interaction. The high-LD groups of macrostate 4 display parallel amino acid variation in these regions. This pattern of haplotype clades with high LD among multiple varying sites, emerging from chaotically recombining variation, may be a “signature” of refinement of complex adaptive sequences by recombination and selection. It should be tested further in this study system and others as a possibly general feature of the evolution of living complexity.     

Clines and adaptive evolution in the methuselah gene region in Drosophila melanogaster

In an effort to characterize further the patterns of selection and adaptive evolution at the methuselah locus in Drosophila species, we extended an analysis of geographical variation to include single nucleotide polymorphisms (SNPs) in adjacent genes on either side of the mth locus, and examined the molecular variation in a neighbouring methuselah paralogue (mth2). An analysis of 13 SNPs spanning a region of nearly 19 kilobases surrounding the mth locus demonstrated that a clinal pattern associated with the most common mth haplotype does not extend to adjacent gene loci, providing compelling evidence that the clinal pattern results from selection on as yet unidentified sites associated with the functional mth locus. mth2 exhibited a significant pattern of adaptive divergence among D. melanogaster, D. simulans and D. yakuba similar to that seen at mth. However, Ka : Ks ratios indicate a difference in levels of functional constraint at the two methuselah, loci with mth2 exhibiting a five- to six-fold reduction in levels of amino acid divergence relative to mth.     

Polymorphism and divergence at the prune locus in Drosophila melanogaster and D. simulans

The prune locus of Drosophila melanogaster lies at the tip of the X chromosome, in a region of reduced recombination in which nearby loci show reduced variation relative to evolutionary divergence from D. simulans. DNA sequencing of prune alleles from D. melanogaster and D. simulans reveals extremely low variation in D. melanogaster but greater variation in D. simulans. Divergence between the two species is not reduced. This pattern may be explained by either positive selection leading to hitchhiking of neutral variation or background selection against deleterious mutations. The pattern of silent versus replacement polymorphism and divergence at prune is consistent with either a model of weakly deleterious selection against amino acid substitutions or balancing selection.      

DNA sequence polymorphisms in exonic and intronic regions of the human phenylalanine hydroxylase gene aid in the identification of alleles

Five sequence polymorphisms at the phenylalanine hydroxylase (PAH) gene locus were observed to be in tight association with specific alleles of this locus. Since these polymorphisms can be detected using polymerase chain reaction (PCR) methodology, application of a combination of these polymorphisms reduces the effort involved in PAH DNA haplotype analysis, which is needed for population genetic analysis or diagnosis of the disease status. In addition our results indicate the evolution of haplotype 3, 4 and 7 PAH alleles from a common ancestor, whereas PAH haplotypes 5, 6, and 11 arose form another common ancestor allele. These data reveal that two of the polymorphisms investigated originated before the separation of races.     

Balancing Selection on a Regulatory Region Exhibiting Ancient Variation That Predates Human–Neandertal Divergence

Ancient population structure shaping contemporary genetic variation has been recently appreciated and has important implications regarding our understanding of the structure of modern human genomes. We identified a ∼36-kb DNA segment in the human genome that displays an ancient substructure. The variation at this locus exists primarily as two highly divergent haplogroups. One of these haplogroups (the NE1 haplogroup) aligns with the Neandertal haplotype and contains a 4.6-kb deletion polymorphism in perfect linkage disequilibrium with 12 single nucleotide polymorphisms (SNPs) across diverse populations. The other haplogroup, which does not contain the 4.6-kb deletion, aligns with the chimpanzee haplotype and is likely ancestral. Africans have higher overall pairwise differences with the Neandertal haplotype than Eurasians do for this NE1 locus (p<10⁻¹⁵). Moreover, the nucleotide diversity at this locus is higher in Eurasians than in Africans. These results mimic signatures of recent Neandertal admixture contributing to this locus. However, an in-depth assessment of the variation in this region across multiple populations reveals that African NE1 haplotypes, albeit rare, harbor more sequence variation than NE1 haplotypes found in Europeans, indicating an ancient African origin of this haplogroup and refuting recent Neandertal admixture. Population genetic analyses of the SNPs within each of these haplogroups, along with genome-wide comparisons revealed significant F_ST (p = 0.00003) and positive Tajima''s D (p = 0.00285) statistics, pointing to non-neutral evolution of this locus. The NE1 locus harbors no protein-coding genes, but contains transcribed sequences as well as sequences with putative regulatory function based on bioinformatic predictions and in vitro experiments. We postulate that the variation observed at this locus predates Human–Neandertal divergence and is evolving under balancing selection, especially among European populations.     

Dimorphic DNA variation in the anionic peroxidase gene AtPrx53 of Arabidopsis thaliana

The DNA polymorphism in the AtPrx53 gene which encodes anionic peroxidase was analyzed in 20 Arabidopsis thaliana accessions. There are two divergent sequence types (Col and Dj-like haplotypes) in the AtPrx53 gene that differ by 2 indel and 16 non-singleton nucleotide polymorphisms including 5 nucleotide polymorphic sites responsible for 4 deduced amino acid replacements. Two of the amino acid substitutions (Phe/Ser180and Asp/Asn270) could be responsible for the difference in electrophoretic mobility of AtPrx53 allozymes. One of them (Phe/Ser180) lies within the hypervariable region, indicating that this amino acid polymorphism is subjected to balancing selection. The revealed difference between deduced allozymes is related to the dimorphism in mobility of three major anionic peroxidase isoforms which according to previously established data encoded by AtPrx53 gene. The haplotype Col which included 12 accessions from three different continents is characterized by faster mobility of three isoforms in comparison with the Dj haplotype represented by eight accessions. There is a significant association between the haplotype and several developmental traits: leaf number, flowering time, main stem height etc. Lines of the Dj haplotype have shorter duration of vegetative stages and flower earlier than most of Col haplotype accessions. The reasons of this association are discussed.     

The functional impact of Pgm amino acid polymorphism on glycogen content in Drosophila melanogaster

Earlier studies of the common PGM allozymes in Drosophila melanogaster reported no in vitro activity differences. However, our study of nucleotide variation observed that PGM allozymes are a heterogeneous mixture of amino acid polymorphisms. In this study, we analyze 10 PGM protein haplotypes with respect to PGM activity, thermostability, and adult glycogen content. We find a twofold difference in activity among PGM protein haplotypes that is associated with a threefold difference in glycogen content. The latitudinal clines for several Pgm amino acid polymorphisms show that high PGM activity, and apparently higher flux to glycogen synthesis, parallel the low activity clines at G6PD for reduced pentose shunt flux in northern latitudes. This suggests that amino acid polymorphism is under selection at this branch point and may be favored for increased metabolic storage associated with stress resistance and adaptation to temperate regions.     

Apolipoprotein E variation at the sequence haplotype level: implications for the origin and maintenance of a major human polymorphism

Three common protein isoforms of apolipoprotein E (apoE), encoded by the epsilon2, epsilon3, and epsilon4 alleles of the APOE gene, differ in their association with cardiovascular and Alzheimer's disease risk. To gain a better understanding of the genetic variation underlying this important polymorphism, we identified sequence haplotype variation in 5.5 kb of genomic DNA encompassing the whole of the APOE locus and adjoining flanking regions in 96 individuals from four populations: blacks from Jackson, MS (n=48 chromosomes), Mayans from Campeche, Mexico (n=48), Finns from North Karelia, Finland (n=48), and non-Hispanic whites from Rochester, MN (n=48). In the region sequenced, 23 sites varied (21 single nucleotide polymorphisms, or SNPs, 1 diallelic indel, and 1 multiallelic indel). The 22 diallelic sites defined 31 distinct haplotypes in the sample. The estimate of nucleotide diversity (site-specific heterozygosity) for the locus was 0.0005+/-0.0003. Sequence analysis of the chimpanzee APOE gene showed that it was most closely related to human epsilon4-type haplotypes, differing from the human consensus sequence at 67 synonymous (54 substitutions and 13 indels) and 9 nonsynonymous fixed positions. The evolutionary history of allelic divergence within humans was inferred from the pattern of haplotype relationships. This analysis suggests that haplotypes defining the epsilon3 and epsilon2 alleles are derived from the ancestral epsilon4s and that the epsilon3 group of haplotypes have increased in frequency, relative to epsilon4s, in the past 200,000 years. Substantial heterogeneity exists within all three classes of sequence haplotypes, and there are important interpopulation differences in the sequence variation underlying the protein isoforms that may be relevant to interpreting conflicting reports of phenotypic associations with variation in the common protein isoforms.     

Selective sweep at the Drosophila melanogaster Suppressor of Hairless locus and its association with the In(2L)t inversion polymorphism.

The hitchhiking model of population genetics predicts that an allele favored by Darwinian selection can replace haplotypes from the same locus previously established at a neutral mutation-drift equilibrium. This process, known as "selective sweep," was studied by comparing molecular variation between the polymorphic In(2L)t inversion and the standard chromosome. Sequence variation was recorded at the Suppressor of Hairless (Su[H]) gene in an African population of Drosophila melanogaster. We found 47 nucleotide polymorphisms among 20 sequences of 1.2 kb. Neutrality tests were nonsignificant at the nucleotide level. However, these sites were strongly associated, because 290 out of 741 observed pairwise combinations between them were in significant linkage disequilibrium. We found only seven haplotypes, two occurring in the 9 In(2L)t chromosomes, and five in the 11 standard chromosomes, with no shared haplotype. Two haplotypes, one in each chromosome arrangement, made up two-thirds of the sample. This low haplotype diversity departed from neutrality in a haplotype test. This pattern supports a selective sweep hypothesis for the Su(H) chromosome region.     

An anomalous haplotype distribution of the arrestin domain-containing 4 gene (ARRDC4) haplotypes in Caucasians

Little was known about the sequence variability of the human Arrestin domain-containing 4 gene (ARRDC4). We sequenced its DNA from exon 2 to exon 8 in a sample of 92 Russians. Seven variants were identified; one of them has not been described yet. It causes an amino acid change from Thr to Met. Identified variants were genotyped in the complete sample of 253 unrelated men and women to analyze haplotype distribution. Fifteen haplotypes were inferred. Nine haplotypes had estimated frequencies > 1%. Ninety-five percent of all haplotypes were determined by five haplotype-tagging single nucleotide polymorphisms. Haplotypes form two clades. The two most common haplotypes cover 76% of all haplotypes. The certainty of the haplotype reconstruction does not depend on the haplotype-inferring algorithms, but is a result of the anomalous haplotype distribution of ARRDC4, which makes this gene a suitable candidate gene for haplotype association studies. Interestingly, there is a great evolutionary distance between the two most common haplotypes, which could suggest a more complicated coalescent process with either past gene flow, selections, or bottlenecks.     

Novel PKU mutation on haplotype 2 in French-Canadians.

We analyzed DNA from nine French-Canadian probands from eastern Quebec province; all had hyperphenylalaninemia (phenylketonuria [PKU] or non-PKU forms) caused by mutations at the phenylalanine hydroxylase locus. Analysis of RFLP haplotypes and mutations revealed a novel mutation, an A-to-G transition (met----val) in codon 1 (the translation-initiation codon). It occurred on 5 of the 18 mutant chromosomes and was associated each time with haplotype 2. A proband homozygous for this mutation had the PKU phenotype. In other probands, the codon 1 mutation was inherited once with the splice junction mutation in exon 12 (on haplotype 3), conferring PKU, and was inherited twice with a mutation on haplotype 1, conferring PKU in one proband and non-PKU hyperphenylalaninemia in the other. The other five probands carried mutations, conferring PKU, on the following haplotype combinations: 1/3 (twice), 1/9, 3/4, and 1/1. The mutations on haplotypes 1, 4, and 9 are not yet characterized. This preliminary study reveals a novel PKU mutation and considerable genetic heterogeneity at the phenylalanine hydroxylase locus in French-Canadians.