A general approach to single-nucleotide polymorphism discovery

Single-nucleotide polymorphisms (SNPs) are the most abundant form of human genetic variation and a resource for mapping complex genetic traits. The large volume of data produced by high-throughput sequencing projects is a rich and largely untapped source of SNPs (refs 2, 3, 4, 5). We present here a unified approach to the discovery of variations in genetic sequence data of arbitrary DNA sources. We propose to use the rapidly emerging genomic sequence as a template on which to layer often unmapped, fragmentary sequence data and to use base quality values to discern true allelic variations from sequencing errors. By taking advantage of the genomic sequence we are able to use simpler yet more accurate methods for sequence organization: fragment clustering, paralogue identification and multiple alignment. We analyse these sequences with a novel, Bayesian inference engine, POLYBAYES, to calculate the probability that a given site is polymorphic. Rigorous treatment of base quality permits completely automated evaluation of the full length of all sequences, without limitations on alignment depth. We demonstrate this approach by accurate SNP predictions in human ESTs aligned to finished and working-draft quality genomic sequences, a data set representative of the typical challenges of sequence-based SNP discovery.     

Dichotomy of single-nucleotide polymorphism haplotypes in olfactory receptor genes and pseudogenes

Substantial efforts are focused on identifying single-nucleotide polymorphisms (SNPs) throughout the human genome, particularly in coding regions (cSNPs), for both linkage disequilibrium and association studies. Less attention, however, has been directed to the clarification of evolutionary processes that are responsible for the variability in nucleotide diversity among different regions of the genome. We report here the population sequence diversity of genomic segments within a 450-kb cluster of olfactory receptor (OR) genes on human chromosome 17. We found a dichotomy in the pattern of nucleotide diversity between OR pseudogenes and introns on the one hand and the closely interspersed intact genes on the other. We suggest that weak positive selection is responsible for the observed patterns of genetic variation. This is inferred from a lower ratio of polymorphism to divergence in genes compared with pseudogenes or introns, high non-synonymous substitution rates in OR genes, and a small but significant overall reduction in variability in the entire OR gene cluster compared with other genomic regions. The dichotomy among functionally different segments within a short genomic distance requires high recombination rates within this OR cluster. Our work demonstrates the impact of weak positive selection on human nucleotide diversity, and has implications for the evolution of the olfactory repertoire.     

Determination of ancestral alleles for human single-nucleotide polymorphisms using high-density oligonucleotide arrays.

Here we report the application of high-density oligonucleotide array (DNA chip)-based analysis to determine the distant history of single nucleotide polymorphisms (SNPs) in current human populations. We analysed orthologues for 397 human SNP sites (identified in CEPH pedigrees from Amish, Venezuelan and Utah populations) from 23 common chimpanzee, 19 pygmy chimpanzee and 11 gorilla genomic DNA samples. From this data we determined 214 proposed ancestral alleles (the sequence found in the last common ancestor of humans and chimpanzees). In a diverse human population set, we found that SNP alleles with higher frequencies were more likely to be ancestral than less frequently occurring alleles. There were, however, exceptions. We also found three shared human/pygmy chimpanzee polymorphisms, all involving CpG dinucleotides, and two shared human/gorilla polymorphisms, one involving a CpG dinucleotide. We demonstrate that microarray-based assays allow rapid comparative sequence analysis of intra- and interspecies genetic variation.     

A high-resolution survey of deletion polymorphism in the human genome

Recent work has shown that copy number polymorphism is an important class of genetic variation in human genomes. Here we report a new method that uses SNP genotype data from parent-offspring trios to identify polymorphic deletions. We applied this method to data from the International HapMap Project to produce the first high-resolution population surveys of deletion polymorphism. Approximately 100 of these deletions have been experimentally validated using comparative genome hybridization on tiling-resolution oligonucleotide microarrays. Our analysis identifies a total of 586 distinct regions that harbor deletion polymorphisms in one or more of the families. Notably, we estimate that typical individuals are hemizygous for roughly 30-50 deletions larger than 5 kb, totaling around 550-750 kb of euchromatic sequence across their genomes. The detected deletions span a total of 267 known and predicted genes. Overall, however, the deleted regions are relatively gene-poor, consistent with the action of purifying selection against deletions. Deletion polymorphisms may well have an important role in the genetics of complex traits; however, they are not directly observed in most current gene mapping studies. Our new method will permit the identification of deletion polymorphisms in high-density SNP surveys of trio or other family data.     

The human Rhesus-associated RhAG protein and a kidney homologue promote ammonium transport in yeast

The Rhesus blood-group antigens are defined by a complex association of membrane polypeptides that includes the non-glycosylated Rh proteins (RhD and RhCE) and the RHag glycoprotein, which is strictly required for cell surface expression of these antigens. RhAG and the Rh polypeptides are erythroid-specific transmembrane proteins belonging to the same family (36% identity). Despite their importance in transfusion medicine, the function of RhAG and Rh proteins remains unknown, except that their absence in Rh(null) individuals leads to morphological and functional abnormalities of erythrocytes, known as the Rh-deficiency syndrome. We recently found significant sequence similarity between the Rh family proteins, especially RhAG, and Mep/Amt ammonium transporters. We show here that RhAG and also RhGK, a new human homologue expressed in kidney cells only, function as ammonium transport proteins when expressed in yeast. Both specifically complement the growth defect of a yeast mutant deficient in ammonium uptake. Moreover, ammonium efflux assays and growth tests in the presence of toxic concentrations of the analogue methylammonium indicate that RhAG and RhGK also promote ammonium export. Our results provide the first experimental evidence for a direct role of RhAG and RhGK in ammonium transport. These findings are of high interest, because no specific ammonium transport system has been characterized so far in human.     

A road map for efficient and reliable human genome epidemiology

Networks of investigators have begun sharing best practices, tools and methods for analysis of associations between genetic variation and common diseases. A Network of Investigator Networks has been set up to drive the process, sponsored by the Human Genome Epidemiology Network. A workshop is planned to develop consensus guidelines for reporting results of genetic association studies. Published literature databases will be integrated, and unpublished data, including 'negative' studies, will be captured by online journals and through investigator networks. Systematic reviews will be expanded to include more meta-analyses of individual-level data and prospective meta-analyses. Field synopses will offer regularly updated overviews.     

A humanized model for multiple sclerosis using HLA-DR2 and a human T-cell receptor

Multiple sclerosis (MS) is a complex chronic neurologic disease with a suspected autoimmune pathogenesis. Although there is evidence that the development of MS is determined by both environmental influences and genes, these factors are largely undefined, except for major histocompatibility (MHC) genes. Linkage analyses and association studies have shown that susceptibility to MS is associated with genes in the human histocompatibility leukocyte antigens (HLA) class II region, but the contribution of these genes to MS disease development less compared with their contribution to disorders such as insulin-dependent diabetes mellitus. Due to the strong linkage disequilibrium in the MHC class II region, it has not been possible to determine which gene(s) is responsible for the genetic predisposition. In transgenic mice, we have expressed three human components involved in T-cell recognition of an MS-relevant autoantigen presented by the HLA-DR2 molecule: DRA*0101/DRB1*1501 (HLA-DR2), an MHC class II candidate MS susceptibility genes found in individuals of European descent; a T-cell receptor (TCR) from an MS-patient-derived T-cell clone specific for the HLA-DR2 bound immunodominant myelin basic protein (MBP) 4102 peptide; and the human CD4 coreceptor. The amino acid sequence of the MBP 84-102 peptide is the same in both human and mouse MBP. Following administration of the MBP peptide, together with adjuvant and pertussis toxin, transgenic mice developed focal CNS inflammation and demyelination that led to clinical manifestations and disease courses resembling those seen in MS. Spontaneous disease was observed in 4% of mice. When DR2 and TCR double-transgenic mice were backcrossed twice to Rag2 (for recombination-activating gene 2)-deficient mice, the incidence of spontaneous disease increased, demonstrating that T cells specific for the HLA-DR2 bound MBP peptide are sufficient and necessary for development of disease. Our study provides evidence that HLA-DR2 can mediate both induced and spontaneous disease resembling MS by presenting an MBP self-peptide to T cells.     

A literature network of human genes for high-throughput analysis of gene expression

We have carried out automated extraction of explicit and implicit biomedical knowledge from publicly available gene and text databases to create a gene-to-gene co-citation network for 13,712 named human genes by automated analysis of titles and abstracts in over 10 million MEDLINE records. The associations between genes have been annotated by linking genes to terms from the medical subject heading (MeSH) index and terms from the gene ontology (GO) database. The extracted database and accompanying web tools for gene-expression analysis have collectively been named 'PubGene'. We validated the extracted networks by three large-scale experiments showing that co-occurrence reflects biologically meaningful relationships, thus providing an approach to extract and structure known biology. We validated the applicability of the tools by analyzing two publicly available microarray data sets.     

A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC

The proteins encoded by the classical HLA class I and class II genes in the major histocompatibility complex (MHC) are highly polymorphic and are essential in self versus non-self immune recognition. HLA variation is a crucial determinant of transplant rejection and susceptibility to a large number of infectious and autoimmune diseases. Yet identification of causal variants is problematic owing to linkage disequilibrium that extends across multiple HLA and non-HLA genes in the MHC. We therefore set out to characterize the linkage disequilibrium patterns between the highly polymorphic HLA genes and background variation by typing the classical HLA genes and >7,500 common SNPs and deletion-insertion polymorphisms across four population samples. The analysis provides informative tag SNPs that capture much of the common variation in the MHC region and that could be used in disease association studies, and it provides new insight into the evolutionary dynamics and ancestral origins of the HLA loci and their haplotypes.     

A gene harvest revealing the archeology and complexity of human disease

An inflection point in genetics has been reached and was witnessed at the Genomics of Common Diseases meeting (http://www.nature.com/ng/meetings/genomics) co-sponsored by the Wellcome Trust, held at the Wellcome Trust Hinxton Conference Centre July 7-10, 2007 and co-organized by Tim Aitman, David Altshuler and Eddy Rubin.     

Loss of DMP1 causes rickets and osteomalacia and identifies a role for osteocytes in mineral metabolism

The osteocyte, a terminally differentiated cell comprising 90%-95% of all bone cells, may have multiple functions, including acting as a mechanosensor in bone (re)modeling. Dentin matrix protein 1 (encoded by DMP1) is highly expressed in osteocytes and, when deleted in mice, results in a hypomineralized bone phenotype. We investigated the potential for this gene not only to direct skeletal mineralization but also to regulate phosphate (P(i)) homeostasis. Both Dmp1-null mice and individuals with a newly identified disorder, autosomal recessive hypophosphatemic rickets, manifest rickets and osteomalacia with isolated renal phosphate-wasting associated with elevated fibroblast growth factor 23 (FGF23) levels and normocalciuria. Mutational analyses showed that autosomal recessive hypophosphatemic rickets family carried a mutation affecting the DMP1 start codon, and a second family carried a 7-bp deletion disrupting the highly conserved DMP1 C terminus. Mechanistic studies using Dmp1-null mice demonstrated that absence of DMP1 results in defective osteocyte maturation and increased FGF23 expression, leading to pathological changes in bone mineralization. Our findings suggest a bone-renal axis that is central to guiding proper mineral metabolism.     

A worldwide survey of haplotype variation and linkage disequilibrium in the human genome

Recent genomic surveys have produced high-resolution haplotype information, but only in a small number of human populations. We report haplotype structure across 12 Mb of DNA sequence in 927 individuals representing 52 populations. The geographic distribution of haplotypes reflects human history, with a loss of haplotype diversity as distance increases from Africa. Although the extent of linkage disequilibrium (LD) varies markedly across populations, considerable sharing of haplotype structure exists, and inferred recombination hotspot locations generally match across groups. The four samples in the International HapMap Project contain the majority of common haplotypes found in most populations: averaging across populations, 83% of common 20-kb haplotypes in a population are also common in the most similar HapMap sample. Consequently, although the portability of tag SNPs based on the HapMap is reduced in low-LD Africans, the HapMap will be helpful for the design of genome-wide association mapping studies in nearly all human populations.     

A genome-wide association scan identifies the hepatic cholesterol transporter ABCG8 as a susceptibility factor for human gallstone disease

With an overall prevalence of 10-20%, gallstone disease (cholelithiasis) represents one of the most frequent and economically relevant health problems of industrialized countries. We performed an association scan of >500,000 SNPs in 280 individuals with gallstones and 360 controls. A follow-up study of the 235 most significant SNPs in 1,105 affected individuals and 873 controls replicated the disease association of SNP A-1791411 in ABCG8 (allelic P value P(CCA) = 4.1 x 10(-9)), which was subsequently attributed to coding variant rs11887534 (D19H). Additional replication was achieved in 728 German (P = 2.8 x 10(-7)) and 167 Chilean subjects (P = 0.02). The overall odds ratio for D19H carriership was 2.2 (95% confidence interval: 1.8-2.6, P = 1.4 x 10(-14)) in the full German sample. Association was stronger in subjects with cholesterol gallstones (odds ratio = 3.3), suggesting that His19 might be associated with a more efficient transport of cholesterol into the bile.     

A genetic interaction network of five genes for human polycystic kidney and liver diseases defines polycystin-1 as the central determinant of cyst formation

Autosomal dominant polycystic liver disease results from mutations in PRKCSH or SEC63. The respective gene products, glucosidase IIβ and SEC63p, function in protein translocation and quality control pathways in the endoplasmic reticulum. Here we show that glucosidase IIβ and Sec63p are required in mice for adequate expression of a functional complex of the polycystic kidney disease gene products, polycystin-1 and polycystin-2. We find that polycystin-1 is the rate-limiting component of this complex and that there is a dose-response relationship between cystic dilation and levels of functional polycystin-1 following mutation of Prkcsh or Sec63. Reduced expression of polycystin-1 also serves to sensitize the kidney to cyst formation resulting from mutations in Pkhd1, the recessive polycystic kidney disease gene. Finally, we show that proteasome inhibition increases steady-state levels of polycystin-1 in cells lacking glucosidase IIβ and that treatment with a proteasome inhibitor reduces cystic disease in orthologous gene models of human autosomal dominant polycystic liver disease.     

A novel GC-rich human macrosatellite VNTR in Xq24 is differentially methylated on active and inactive X chromosomes.

A new X chromosome-specific repetitive sequence, a 3 kilobase HindIII clone with a base composition of 63% C+G, has been isolated. The sequence is organized as a hypervariable tandem repeat cluster ranging in size from 150-350 kilobases, with outlying single copies. This locus, designated DXZ4 and mapped to chromosome band Xq24, may consist of as many as 50 variable-length alleles. It represents a class of variable number of tandem repeat polymorphism which may be termed 'macrosatellite'. The cluster is highly methylated on the active X chromosome and hypomethylated on the inactive X.