Genome-wide detection and characterization of positive selection in human populations

With the advent of dense maps of human genetic variation, it is now possible to detect positive natural selection across the human genome. Here we report an analysis of over 3 million polymorphisms from the International HapMap Project Phase 2 (HapMap2). We used 'long-range haplotype' methods, which were developed to identify alleles segregating in a population that have undergone recent selection, and we also developed new methods that are based on cross-population comparisons to discover alleles that have swept to near-fixation within a population. The analysis reveals more than 300 strong candidate regions. Focusing on the strongest 22 regions, we develop a heuristic for scrutinizing these regions to identify candidate targets of selection. In a complementary analysis, we identify 26 non-synonymous, coding, single nucleotide polymorphisms showing regional evidence of positive selection. Examination of these candidates highlights three cases in which two genes in a common biological process have apparently undergone positive selection in the same population:LARGE and DMD, both related to infection by the Lassa virus, in West Africa;SLC24A5 and SLC45A2, both involved in skin pigmentation, in Europe; and EDAR and EDA2R, both involved in development of hair follicles, in Asia.     

Resolution of quantitative traits into Mendelian factors by using a complete linkage map of restriction fragment length polymorphisms

The conflict between the Mendelian theory of particulate inheritance and the observation of continuous variation for most traits in nature was resolved in the early 1900s by the concept that quantitative traits can result from segregation of multiple genes, modified by environmental effects. Although pioneering experiments showed that linkage could occasionally be detected to such quantitative trait loci (QTLs), accurate and systematic mapping of QTLs has not been possible because the inheritance of an entire genome could not be studied with genetic markers. The use of restriction fragment length polymorphisms (RFLPs) has made such investigations possible, at least in principle. Here, we report the first use of a complete RFLP linkage map to resolve quantitative traits into discrete Mendelian factors, in an interspecific back-cross of tomato. Applying new analytical methods, we mapped at least six QTLs controlling fruit mass, four QTLs for the concentration of soluble solids and five QTLs for fruit pH. This approach is broadly applicable to the genetic dissection of quantitative inheritance of physiological, morphological and behavioural traits in any higher plant or animal.     

Proof and evolutionary analysis of ancient genome duplication in the yeast Saccharomyces cerevisiae

Whole-genome duplication followed by massive gene loss and specialization has long been postulated as a powerful mechanism of evolutionary innovation. Recently, it has become possible to test this notion by searching complete genome sequence for signs of ancient duplication. Here, we show that the yeast Saccharomyces cerevisiae arose from ancient whole-genome duplication, by sequencing and analysing Kluyveromyces waltii, a related yeast species that diverged before the duplication. The two genomes are related by a 1:2 mapping, with each region of K. waltii corresponding to two regions of S. cerevisiae, as expected for whole-genome duplication. This resolves the long-standing controversy on the ancestry of the yeast genome, and makes it possible to study the fate of duplicated genes directly. Strikingly, 95% of cases of accelerated evolution involve only one member of a gene pair, providing strong support for a specific model of evolution, and allowing us to distinguish ancestral and derived functions.     

Cornelia de Lange syndrome is caused by mutations in NIPBL, the human homolog of Drosophila melanogaster Nipped-B

Cornelia de Lange syndrome (CdLS; OMIM 122470) is a dominantly inherited multisystem developmental disorder characterized by growth and cognitive retardation; abnormalities of the upper limbs; gastroesophageal dysfunction; cardiac, ophthalmologic and genitourinary anomalies; hirsutism; and characteristic facial features. Genital anomalies, pyloric stenosis, congenital diaphragmatic hernias, cardiac septal defects, hearing loss and autistic and self-injurious tendencies also frequently occur. Prevalence is estimated to be as high as 1 in 10,000 (ref. 4). We carried out genome-wide linkage exclusion analysis in 12 families with CdLS and identified four candidate regions, of which chromosome 5p13.1 gave the highest multipoint lod score of 2.7. This information, together with the previous identification of a child with CdLS with a de novo t(5;13)(p13.1;q12.1) translocation, allowed delineation of a 1.1-Mb critical region on chromosome 5 for the gene mutated in CdLS. We identified mutations in one gene in this region, which we named NIPBL, in four sporadic and two familial cases of CdLS. We characterized the genomic structure of NIPBL and found that it is widely expressed in fetal and adult tissues. The fly homolog of NIPBL, Nipped-B, facilitates enhancer-promoter communication and regulates Notch signaling and other developmental pathways in Drosophila melanogaster.     

Assessing the impact of population stratification on genetic association studies

Population stratification refers to differences in allele frequencies between cases and controls due to systematic differences in ancestry rather than association of genes with disease. It has been proposed that false positive associations due to stratification can be controlled by genotyping a few dozen unlinked genetic markers. To assess stratification empirically, we analyzed data from 11 case-control and case-cohort association studies. We did not detect statistically significant evidence for stratification but did observe that assessments based on a few dozen markers lack power to rule out moderate levels of stratification that could cause false positive associations in studies designed to detect modest genetic risk factors. After increasing the number of markers and samples in a case-cohort study (the design most immune to stratification), we found that stratification was in fact present. Our results suggest that modest amounts of stratification can exist even in well designed studies.     

High-resolution haplotype structure in the human genome

Linkage disequilibrium (LD) analysis is traditionally based on individual genetic markers and often yields an erratic, non-monotonic picture, because the power to detect allelic associations depends on specific properties of each marker, such as frequency and population history. Ideally, LD analysis should be based directly on the underlying haplotype structure of the human genome, but this structure has remained poorly understood. Here we report a high-resolution analysis of the haplotype structure across 500 kilobases on chromosome 5q31 using 103 single-nucleotide polymorphisms (SNPs) in a European-derived population. The results show a picture of discrete haplotype blocks (of tens to hundreds of kilobases), each with limited diversity punctuated by apparent sites of recombination. In addition, we develop an analytical model for LD mapping based on such haplotype blocks. If our observed structure is general (and published data suggest that it may be), it offers a coherent framework for creating a haplotype map of the human genome.     

A radiation hybrid map of mouse genes

A comprehensive gene-based map of a genome is a powerful tool for genetic studies and is especially useful for the positional cloning and positional candidate approaches. The availability of gene maps for multiple organisms provides the foundation for detailed conserved-orthology maps showing the correspondence between conserved genomic segments. These maps make it possible to use cross-species information in gene hunts and shed light on the evolutionary forces that shape the genome. Here we report a radiation hybrid map of mouse genes, a combined project of the Whitehead Institute/Massachusetts Institute of Technology Center for Genome Research, the Medical Research Council UK Mouse Genome Centre, and the National Center for Biotechnology Information. The map contains 11,109 genes, screened against the T31 RH panel and positioned relative to a reference map containing 2,280 mouse genetic markers. It includes 3,658 genes homologous to the human genome sequence and provides a framework for overlaying the human genome sequence to the mouse and for sequencing the mouse genome.     

Linkage disequilibrium in the human genome

With the availability of a dense genome-wide map of single nucleotide polymorphisms (SNPs), a central issue in human genetics is whether it is now possible to use linkage disequilibrium (LD) to map genes that cause disease. LD refers to correlations among neighbouring alleles, reflecting 'haplotypes' descended from single, ancestral chromosomes. The size of LD blocks has been the subject of considerable debate. Computer simulations and empirical data have suggested that LD extends only a few kilobases (kb) around common SNPs, whereas other data have suggested that it can extend much further, in some cases greater than 100 kb. It has been difficult to obtain a systematic picture of LD because past studies have been based on only a few (1-3) loci and different populations. Here, we report a large-scale experiment using a uniform protocol to examine 19 randomly selected genomic regions. LD in a United States population of north-European descent typically extends 60 kb from common alleles, implying that LD mapping is likely to be practical in this population. By contrast, LD in a Nigerian population extends markedly less far. The results illuminate human history, suggesting that LD in northern Europeans is shaped by a marked demographic event about 27,000-53,000 years ago.     

ARSACS, a spastic ataxia common in northeastern Québec, is caused by mutations in a new gene encoding an 11.5-kb ORF

Autosomal recessive spastic ataxia of Charlevoix-Saguenay (ARSACS or SACS) is an early onset neurodegenerative disease with high prevalence (carrier frequency 1/22) in the Charlevoix-Saguenay-Lac-Saint-Jean (CSLSJ) region of Quebec. We previously mapped the gene responsible for ARSACS to chromosome 13q11 and identified two ancestral haplotypes. Here we report the cloning of this gene, SACS, which encodes the protein sacsin. The ORF of SACS is 11,487 bp and is encoded by a single gigantic exon spanning 12,794 bp. This exon is the largest to be identified in any vertebrate organism. The ORF is conserved in human and mouse. The putative protein contains three large segments with sequence similarity to each other and to the predicted protein of an Arabidopsis thaliana ORF. The presence of heat-shock domains suggests a function for sacsin in chaperone-mediated protein folding. SACS is expressed in a variety of tissues, including the central nervous system. We identified two SACSmutations in ARSACS families that lead to protein truncation, consistent with haplotype analysis.