In this report, we compare and contrast three previously published Bayesian methods for inferring haplotypes from genotype data in a population sample. We review the methods, emphasizing the differences between them in terms of both the models ("priors") they use and the computational strategies they employ. We introduce a new algorithm that combines the modeling strategy of one method with the computational strategies of another. In comparisons using real and simulated data, this new algorithm outperforms all three existing methods. The new algorithm is included in the software package PHASE, version 2.0, available online (http://www.stat.washington.edu/stephens/software.html).

Retinitis pigmentosa (RP), the main cause of adult blindness, is a genetically heterogeneous disorder characterized by progressive loss of photoreceptors through apoptosis. Up to now, 39 genes and loci have been implicated in nonsyndromic RP, yet the genetic bases of >50% of the cases, particularly of the recessive forms, remain unknown. Previous linkage analysis in a Spanish consanguineous family allowed us to define a novel autosomal recessive RP (arRP) locus, RP26, within an 11-cM interval (17.4 Mb) on 2q31.2-q32.3. In the present study, we further refine the RP26 locus down to 2.5 Mb, by microsatellite and single-nucleotide polymorphism (SNP) homozygosity mapping. After unsuccessful mutational analysis of the nine genes initially reported in this region, a detailed gene search based on expressed-sequence-tag data was undertaken. We finally identified a novel gene encoding a ceramide kinase (CERKL), which encompassed 13 exons. All of the patients from the RP26 family bear a homozygous mutation in exon 5, which generates a premature termination codon. The same mutation was also characterized in another, unrelated, Spanish pedigree with arRP. Human CERKL is expressed in the retina, among other adult and fetal tissues. A more detailed analysis by in situ hybridization on adult murine retina sections shows expression of Cerkl in the ganglion cell layer. Ceramide kinases convert the sphingolipid metabolite ceramide into ceramide-1-phosphate, both key mediators of cellular apoptosis and survival. Ceramide metabolism plays an essential role in the viability of neuronal cells, the membranes of which are particularly rich in sphingolipids. Therefore, CERKL deficiency could shift the relative levels of the signaling sphingolipid metabolites and increase sensitivity of photoreceptor and other retinal cells to apoptotic stimuli. This is the first genetic report suggesting a direct link between retinal neurodegeneration in RP and sphingolipid-mediated apoptosis.

It is a basic principle of genetics that each chromosome is transmitted from parent to offspring with a probability that is given by Mendel's laws. However, several known biological processes lead to skewed transmission probabilities among surviving offspring and, therefore, to excess genetic sharing among relatives. Examples include in utero selection against deleterious mutations, meiotic drive, and maternal-fetal incompatibility. Although these processes affect our basic understanding of inheritance, little is known about their overall impact in humans or other mammals. In this study, we examined genome screen data from 148 nuclear families, collected without reference to phenotype, to look for departures from Mendelian transmission proportions. Using single-point and multipoint linkage analysis, we detected a modest but significant genomewide shift towards excess genetic sharing among siblings (average sharing of 50.43% for the autosomes; P=.009). Our calculations indicate that many loci with skewed transmission are required to produce a genomewide shift of this magnitude. Since transmission distortion loci are subject to strong selection, this raises interesting questions about the evolutionary forces that keep them polymorphic. Finally, our results also have implications for mapping disease genes and for the genetics of fertility.

Peak bone mineral density (BMD) is a highly heritable trait and is a good predictor of the risk of osteoporosis and fracture in later life. Recent studies have sought to identify the genes underlying peak BMD. Linkage analysis in a sample of 464 premenopausal white sister pairs detected linkage of spine BMD to chromosome 1q (LOD 3.6). An independent sample of 254 white sister pairs has now been genotyped, and it also provides evidence of linkage to chromosome 1q (LOD 2.5) for spine BMD. Microsatellite markers were subsequently genotyped for a 4-cM map in the chromosome 1q region in all available white sister pairs (n=938), and a LOD score of 4.3 was obtained near the marker D1S445. Studies in the mouse have also detected evidence of linkage to BMD phenotypes in the region syntenic to our linkage finding on chromosome 1q. Thus, we have replicated a locus on 1q contributing to BMD at the spine and have found further support for the region in analyses employing an enlarged sample. Studies are now ongoing to identify the gene(s) contributing to peak spine BMD in women.

Between the 15th and 19th centuries ad, the Atlantic slave trade resulted in the forced movement of approximately 13 million people from Africa, mainly to the Americas. Only approximately 11 million survived the passage, and many more died in the early years of captivity. We have studied 481 mitochondrial DNAs (mtDNAs) of recent African ancestry in the Americas and in Eurasia, in an attempt to trace them back to particular regions of Africa. Our results show that mtDNAs in America and Eurasia can, in many cases, be traced to broad geographical regions within Africa, largely in accordance with historical evidence, and raise the possibility that a greater resolution may be possible in the future. However, they also indicate that, at least for the moment, considerable caution is warranted when assessing claims to be able to trace the ancestry of particular lineages to a particular locality within modern-day Africa.

Hereditary leiomyomatosis and renal cell cancer (HLRCC) is an autosomal dominant disorder characterized by smooth-muscle tumors of the skin and uterus and/or renal cancer. Although the identification of germline mutations in the fumarate hydratase (FH) gene in European families supports it as the susceptibility gene for HLRCC, its role in families in North America has not been studied. We screened for germline mutations in FH in 35 families with cutaneous leiomyomas. Sequence analysis revealed mutations in FH in 31 families (89%). Twenty different mutations in FH were identified, of which 18 were novel. Of these 20 mutations, 2 were insertions, 5 were small deletions that caused frameshifts leading to premature truncation of the protein, and 13 were missense mutations. Eleven unrelated families shared a common mutation: R190H. Eighty-one individuals (47 women and 34 men) had cutaneous leiomyomas. Ninety-eight percent (46/47) of women with cutaneous leiomyomas also had uterine leiomyomas. Eighty-nine percent (41/46) of women with cutaneous and uterine leiomyomas had a total hysterectomy, 44% at age < or =30 years. We identified 13 individuals in 5 families with unilateral and solitary renal tumors. Seven individuals from four families had papillary type II renal cell carcinoma, and another individual from one of these families had collecting duct carcinoma of the kidney. The present study shows that mutations in FH are associated with HLRCC in North America. HLRCC is associated with clinically significant uterine fibroids and aggressive renal tumors. The present study also expands the histologic spectrum of renal tumors and FH mutations associated with HLRCC.

Juvenile myoclonic epilepsy (JME) is a common form of generalized epilepsy that starts in adolescence. A major JME susceptibility locus (EJM1) was mapped to chromosomal region 6p21 in three independent linkage studies, and association was reported between JME and a microsatellite marker in the 6p21 region. The critical region for EJM1 is delimited by obligate recombinants at HLA-DQ and HLA-DP. In the present study, we found highly significant linkage disequilibrium (LD) between JME and a core haplotype of five single-nucleotide-polymorphism (SNP) and microsatellite markers in this critical region, with LD peaking in the BRD2 (RING3) gene (odds ratio 6.45; 95% confidence interval 2.36-17.58). DNA sequencing revealed two JME-associated SNP variants in the BRD2 (RING3) promoter region but no other potentially causative coding mutations in 20 probands from families with positive LOD scores. BRD2 (RING3) is a putative nuclear transcriptional regulator from a family of genes that are expressed during development. Our findings strongly suggest that BRD2 (RING3) is EJM1, the first gene identified for a common idiopathic epilepsy. These findings also suggest that abnormalities of neural development may be a cause of common idiopathic epilepsy, and the findings have implications for the generalizability of proposed pathogenetic mechanisms, derived from diseases that show Mendelian transmission, to their complex counterparts.

Numerous studies have clearly indicated a role for the major histocompatibility complex (MHC) in susceptibility to autoimmune diseases. Such studies have focused on the genetic variation of a small number of classical human-leukocyte-antigen (HLA) genes in the region. Although these genes represent good candidates, given their immunological roles, linkage disequilibrium (LD) surrounding these genes has made it difficult to rule out neighboring genes, many with immune function, as influencing disease susceptibility. It is likely that a comprehensive analysis of the patterns of LD and variation, by using a high-density map of single-nucleotide polymorphisms (SNPs), would enable a greater understanding of the nature of the observed associations, as well as lead to the identification of causal variation. We present herein an initial analysis of this region, using 201 SNPs, nine classical HLA loci, two TAP genes, and 18 microsatellites. This analysis suggests that LD and variation in the MHC, aside from the classical HLA loci, are essentially no different from those in the rest of the genome. Furthermore, these data show that multi-SNP haplotypes will likely be a valuable means for refining association signals in this region.

To avoid problems related to unknown population substructure, association studies may be conducted in founder populations. In such populations, however, the relatedness among individuals may be considerable. Neglecting such correlations among individuals can lead to seriously spurious associations. Here, we propose a method for case-control association studies of binary traits that is suitable for any set of related individuals, provided that their genealogy is known. Although we focus here on large inbred pedigrees, this method may also be used in outbred populations for case-control studies in which some individuals are relatives. We base inference on a quasi-likelihood score (QLS) function and construct a QLS test for allelic association. This approach can be used even when the pedigree structure is far too complex to use an exact-likelihood calculation. We also present an alternative approach to this test, in which we use the known genealogy to derive a correction factor for the case-control association chi2 test. We perform analytical power calculations for each of the two tests by deriving their respective noncentrality parameters. The QLS test is more powerful than the corrected chi2 test in every situation considered. Indeed, under certain regularity conditions, the QLS test is asymptotically the locally most powerful test in a general class of linear tests that includes the corrected chi2 test. The two methods are used to test for associations between three asthma-associated phenotypes and 48 SNPs in 35 candidate genes in the Hutterites. We report a highly significant novel association (P=2.10-6) between atopy and an amino acid polymorphism in the P-selectin gene, detected with the QLS test and also, but less significantly (P=.0014), with the transmission/disequilibrium test.

We report the results of a genomewide scan for age-related macular degeneration (AMD) in 158 multiplex families. AMD classification was based on fundus photography and was assigned a grade ranging from 1 (no disease) to 5 (exudative disease). Genotyping was performed by the National Heart, Lung, and Blood Institute Mammalian Genotyping Service at Marshfield (404 short tandem repeat markers). The sample included 158 families with two or more siblings with AMD, 490 affected individuals, 101 unaffected individuals, and 38 whose affection status was unknown. Relative pairs included 511 affected sibling, 28 avuncular, 53 cousin, 7 grandparent-grandchild, and 9 grand-avuncular pairs. Two-point parametric and multipoint parametric and nonparametric analyses were performed. Maximum two-point LOD scores of 1.0-2.0 were found for markers on chromosomes 1, 2, 8, 10, 14, 15, and 22. Multipoint analyses were consistent with the two-point results for chromosomes 1, 2, 8, 10, and 22 and provided evidence for additional linkage regions on chromosomes 3, 6, 8, 12, 16, and X. Our signals on chromosomes 1q, 6p, and 10q are consistent with some other previously published results. Significant linkage to AMD was found for one marker on chromosome 2, two adjacent markers on chromosome 3, two adjacent markers on chromosome 6, and seven contiguous markers on chromosome 8, with empirical P values of .00001. The consistency of many of the other signals across both two-point and multipoint, as well as parametric and nonparametric, analyses indicate several other regions worthy of follow-up.

Myotonic dystrophy (DM), the most common form of muscular dystrophy in adults, is a clinically and genetically heterogeneous neuromuscular disorder. DM is characterized by autosomal dominant inheritance, muscular dystrophy, myotonia, and multisystem involvement. Type 1 DM (DM1) is caused by a (CTG)(n) expansion in the 3' untranslated region of DMPK in 19q13.3. Multiple families, predominantly of German descent and with clinically variable presentation that included proximal myotonic myopathy (PROMM) and type 2 DM (DM2) but without the DM1 mutation, showed linkage to the 3q21 region and were recently shown to segregate a (CCTG)(n) expansion mutation in intron 1 of ZNF9. Here, we present linkage to 3q21 and mutational confirmation in 17 kindreds of European origin with PROMM and proximal myotonic dystrophy, from geographically distinct populations. All patients have the DM2 (CCTG)(n) expansion. To study the evolution of this mutation, we constructed a comprehensive physical map of the DM2 region around ZNF9. High-resolution haplotype analysis of disease chromosomes with five microsatellite and 22 single-nucleotide polymorphism markers around the DM2 mutation identified extensive linkage disequilibrium and a single shared haplotype of at least 132 kb among patients from the different populations. With the exception of the (CCTG)(n) expansion, the available markers indicate that the DM2 haplotype is identical to the most common haplotype in normal individuals. This situation is reminiscent of that seen in DM1. Taken together, these data suggest a single founding mutation in DM2 patients of European origin. We estimate the age of the founding haplotype and of the DM2 (CCTG) expansion mutation to be approximately 200-540 generations.

Prader-Willi and Angelman syndromes (PWS and AS) typically result from an approximately 4-Mb deletion of human chromosome 15q11-q13, with clustered breakpoints (BP) at either of two proximal sites (BP1 and BP2) and one distal site (BP3). HERC2 and other duplicons map to these BP regions, with the 2-Mb PWS/AS imprinted domain just distal of BP2. Previously, the presence of genes and their imprinted status have not been examined between BP1 and BP2. Here, we identify two known (CYFIP1 and GCP5) and two novel (NIPA1 and NIPA2) genes in this region in human and their orthologs in mouse chromosome 7C. These genes are expressed from a broad range of tissues and are nonimprinted, as they are expressed in cells derived from normal individuals, patients with PWS or AS, and the corresponding mouse models. However, replication-timing studies in the mouse reveal that they are located in a genomic domain showing asynchronous replication, a feature typically ascribed to monoallelically expressed loci. The novel genes NIPA1 and NIPA2 each encode putative polypeptides with nine transmembrane domains, suggesting function as receptors or as transporters. Phylogenetic analyses show that NIPA1 and NIPA2 are highly conserved in vertebrate species, with ancestral members in invertebrates and plants. Intriguingly, evolutionary studies show conservation of the four-gene cassette between BP1 and BP2 in human, including NIPA1/2, CYFIP1, and GCP5, and proximity to the Herc2 gene in both mouse and Fugu. These observations support a model in which duplications of the HERC2 gene at BP3 in primates first flanked the four-gene cassette, with subsequent transposition of these four unique genes by a HERC2 duplicon-mediated process to form the BP1-BP2 region. Duplicons therefore appear to mediate genomic fluidity in both disease and evolutionary processes.

Supernumerary marker chromosomes (SMCs) of chromosome 15, designated "SMC(15)s," are the most common SMC in humans, accounting for as much as 60% of all those observed. We report the characterization of 46 large SMC(15)s, using both fluorescence in situ hybridization and polymerase chain reaction analysis within and distal to the Prader-Willi/Angelman syndrome critical region (PWACR). Our aim was to establish detailed information on origin, content, and breakpoints, to address the formation of SMC(15)s, and to facilitate genotype-phenotype correlations. For all patients in whom we were able to establish the parental origin, the SMC(15)s were maternally derived. Two patients were observed who had familial SMC(15)s, both inherited from the mother; however, in all remaining patients for whom parental samples were available, the SMC(15)s were shown to have arisen de novo. With one exception, all the SMC(15)s were shown to include the entire PWACR. Detailed investigations of the distal breakpoints categorized the SMC(15)s into two groups. Group A, representing approximately two-thirds of the SMC(15)s, had a breakpoint beyond the standard distal PWS/AS deletion breakpoint BP3, at a position close to the microsatellite marker D15S1010 and the bacterial artificial chromosome 10I10. The group B SMC(15)s were shorter, with more variable breakpoints located around BP3. The majority of the SMC(15)s were shown to have asymmetrical breakpoints, with the two inverted arms of the SMC being unequal in length. Our study revealed an unexpected level of complexity and heterogeneity among SMC(15)s that is not seen in other chromosome 15 rearrangements, such as deletions and duplications. This suggests that multiple mechanisms are involved in the formation of large SMC(15)s.

Alzheimer disease (AD) is a complex disorder characterized by a wide range, within and between families, of ages at onset of symptoms. Consideration of age at onset as a covariate in genetic-linkage studies may reduce genetic heterogeneity and increase statistical power. Ordered-subsets analysis includes continuous covariates in linkage analysis by rank ordering families by a covariate and summing LOD scores to find a subset giving a significantly increased LOD score relative to the overall sample. We have analyzed data from 336 markers in 437 multiplex (>/=2 sampled individuals with AD) families included in a recent genomic screen for AD loci. To identify genetic heterogeneity by age at onset, families were ordered by increasing and decreasing mean and minimum ages at onset. Chromosomewide significance of increases in the LOD score in subsets relative to the overall sample was assessed by permutation. A statistically significant increase in the nonparametric multipoint LOD score was observed on chromosome 2q34, with a peak LOD score of 3.2 at D2S2944 (P=.008) in 31 families with a minimum age at onset between 50 and 60 years. The LOD score in the chromosome 9p region previously linked to AD increased to 4.6 at D9S741 (P=.01) in 334 families with minimum age at onset between 60 and 75 years. LOD scores were also significantly increased on chromosome 15q22: a peak LOD score of 2.8 (P=.0004) was detected at D15S1507 (60 cM) in 38 families with minimum age at onset >/=79 years, and a peak LOD score of 3.1 (P=.0006) was obtained at D15S153 (62 cM) in 43 families with mean age at onset >80 years. Thirty-one families were contained in both 15q22 subsets, indicating that these results are likely detecting the same locus. There is little overlap in these subsets, underscoring the utility of age at onset as a marker of genetic heterogeneity. These results indicate that linkage to chromosome 9p is strongest in late-onset AD and that regions on chromosome 2q34 and 15q22 are linked to early-onset AD and very-late-onset AD, respectively.

Carriers of premutation alleles (55-200 CGG repeats) of the fragile-X mental retardation 1 (FMR1) gene are often regarded as being clinically uninvolved. However, it is now apparent that such individuals can present with one (or more) of three distinct clinical disorders: mild cognitive and/or behavioral deficits on the fragile-X spectrum; premature ovarian failure; and a newly described, neurodegenerative disorder of older adult carriers, fragile-X-associated tremor/ataxia syndrome (FXTAS). Awareness of these clinical presentations is important for proper diagnosis and therapeutic intervention, not only among families with known cases of fragile-X syndrome but also more broadly for adults with tremor, gait ataxia, and parkinsonism who are seen in movement-disorders clinics.

We investigated a family with an autosomal recessive syndrome of cafe-au-lait patches and childhood malignancy, notably supratentorial primitive neuroectodermal tumor. There was no cancer predisposition in heterozygotes; nor was there bowel cancer in any individual. However, autozygosity mapping indicated linkage to a region of 7p22 surrounding the PMS2 mismatch-repair gene. Sequencing of genomic PCR products initially failed to identify a PMS2 mutation. Genome searches then revealed a previously unrecognized PMS2 pseudogene, corresponding to exons 9-15, within a 100-kb inverted duplication situated 600 kb centromeric from PMS2 itself. This information allowed a redesigned sequence analysis, identifying a homozygous mutation (R802X) in PMS2 exon 14. Furthermore, in the family with Turcot syndrome, in which the first inherited PMS2 mutation (R134X) was described, a further truncating mutation was identified on the other allele, in exon 13. Further whole-genome analysis shows that the complexity of PMS2 pseudogenes is greater than appreciated and may have hindered previous mutation studies. Several previously reported PMS2 polymorphisms are, in fact, pseudogene sequence variants. Although PMS2 mutations may be rare in colorectal cancer, they appear, for the most part, to behave as recessive traits. For technical reasons, their involvement in childhood cancer, particularly in primitive neuroectodermal tumor, may have been underestimated.

Admixture mapping (also known as "mapping by admixture linkage disequilibrium," or MALD) provides a way of localizing genes that cause disease, in admixed ethnic groups such as African Americans, with approximately 100 times fewer markers than are required for whole-genome haplotype scans. However, it has not been possible to perform powerful scans with admixture mapping because the method requires a dense map of validated markers known to have large frequency differences between Europeans and Africans. To create such a map, we screened through databases containing approximately 450000 single-nucleotide polymorphisms (SNPs) for which frequencies had been estimated in African and European population samples. We experimentally confirmed the frequencies of the most promising SNPs in a multiethnic panel of unrelated samples and identified 3011 as a MALD map (1.2 cM average spacing). We estimate that this map is approximately 70% informative in differentiating African versus European origins of chromosomal segments. This map provides a practical and powerful tool, which is freely available without restriction, for screening for disease genes in African American patient cohorts. The map is especially appropriate for those diseases that differ in incidence between the parental African and European populations.

Juvenile amyotrophic lateral sclerosis (ALS4) is a rare autosomal dominant form of juvenile amyotrophic lateral sclerosis (ALS) characterized by distal muscle weakness and atrophy, normal sensation, and pyramidal signs. Individuals affected with ALS4 usually have an onset of symptoms at age <25 years, a slow rate of progression, and a normal life span. The ALS4 locus maps to a 1.7-Mb interval on chromosome 9q34 flanked by D9S64 and D9S1198. To identify the molecular basis of ALS4, we tested 19 genes within the ALS4 interval and detected missense mutations (T3I, L389S, and R2136H) in the Senataxin gene (SETX). The SETX gene encodes a novel 302.8-kD protein. Although its function remains unknown, SETX contains a DNA/RNA helicase domain with strong homology to human RENT1 and IGHMBP2, two genes encoding proteins known to have roles in RNA processing. These observations of ALS4 suggest that mutations in SETX may cause neuronal degeneration through dysfunction of the helicase activity or other steps in RNA processing.

In most human populations, the ability to digest lactose contained in milk usually disappears in childhood, but in European-derived populations, lactase activity frequently persists into adulthood (Scrimshaw and Murray 1988). It has been suggested (Cavalli-Sforza 1973; Hollox et al. 2001; Enattah et al. 2002; Poulter et al. 2003) that a selective advantage based on additional nutrition from dairy explains these genetically determined population differences (Simoons 1970; Kretchmer 1971; Scrimshaw and Murray 1988; Enattah et al. 2002), but formal population-genetics-based evidence of selection has not yet been provided. To assess the population-genetics evidence for selection, we typed 101 single-nucleotide polymorphisms covering 3.2 Mb around the lactase gene. In northern European-derived populations, two alleles that are tightly associated with lactase persistence (Enattah et al. 2002) uniquely mark a common (~77%) haplotype that extends largely undisrupted for >1 Mb. We provide two new lines of genetic evidence that this long, common haplotype arose rapidly due to recent selection: (1) by use of the traditional F(ST) measure and a novel test based on p(excess), we demonstrate large frequency differences among populations for the persistence-associated markers and for flanking markers throughout the haplotype, and (2) we show that the haplotype is unusually long, given its high frequency--a hallmark of recent selection. We estimate that strong selection occurred within the past 5,000-10,000 years, consistent with an advantage to lactase persistence in the setting of dairy farming; the signals of selection we observe are among the strongest yet seen for any gene in the genome.

The hemoglobin E variant (HbE; ( beta )26Glu-->Lys) is concentrated in parts of Southeast Asia where malaria is endemic, and HbE carrier status has been shown to confer some protection against Plasmodium falciparum malaria. To examine the effect of natural selection on the pattern of linkage disequilibrium (LD) and to infer the evolutionary history of the HbE variant, we analyzed biallelic markers surrounding the HbE variant in a Thai population. Pairwise LD analysis of HbE and 43 surrounding biallelic markers revealed LD of HbE extending beyond 100 kb, whereas no LD was observed between non-HbE variants and the same markers. The inferred haplotype network suggests a single origin of the HbE variant in the Thai population. Forward-in-time computer simulations under a variety of selection models indicate that the HbE variant arose 1,240-4,440 years ago. These results support the conjecture that the HbE mutation occurred recently, and the allele frequency has increased rapidly. Our study provides another clear demonstration that a high-resolution LD map across the human genome can detect recent variants that have been subjected to positive selection.