Search

Patients with chromosome 13q deletion or normal cytogenetics represent the majority of chronic lymphocytic leukaemia (CLL) cases, yet have relatively few driver mutations. To better understand their genomic landscape, here we perform whole-genome sequencing on a cohort of patients enriched with these cytogenetic characteristics. Mutations in known CLL drivers are seen in only 33% of this cohort, and associated with normal cytogenetics and unmutated IGHV. The most commonly mutated gene in our cohort, IGLL5, shows a mutational pattern suggestive of activation-induced cytidine deaminase (AID) activity. Unsupervised analysis of mutational signatures demonstrates the activities of canonical AID (c-AID), leading to clustered mutations near active transcriptional start sites; non-canonical AID (nc-AID), leading to genome-wide non-clustered mutations, and an ageing signature responsible for most mutations. Using mutation clonality to infer time of onset, we find that while ageing and c-AID activities are ongoing, nc-AID-associated mutations likely occur earlier in tumour evolution.

Mosquitoes in the Anopheles gambiae complex show rapid ecological and behavioral diversification, traits that promote malaria transmission and complicate vector control efforts. A high-density, genome-wide mosquito SNP-genotyping array allowed mapping of genomic differentiation between populations and species that exhibit varying levels of reproductive isolation. Regions near centromeres or within polymorphic inversions exhibited the greatest genetic divergence, but divergence was also observed elsewhere in the genomes. Signals of natural selection within populations were overrepresented among genomic regions that are differentiated between populations, implying that differentiation is often driven by population-specific selective events. Complex genomic differentiation among speciating vector mosquito populations implies that tools for genome-wide monitoring of population structure will prove useful for the advancement of malaria eradication.

We report a high-quality draft sequence of the genome of the horse (Equus caballus). The genome is relatively repetitive but has little segmental duplication. Chromosomes appear to have undergone few historical rearrangements: 53% of equine chromosomes show conserved synteny to a single human chromosome. Equine chromosome 11 is shown to have an evolutionary new centromere devoid of centromeric satellite DNA, suggesting that centromeric function may arise before satellite repeat accumulation. Linkage disequilibrium, showing the influences of early domestication of large herds of female horses, is intermediate in length between dog and human, and there is long-range haplotype sharing among breeds.

We describe a collection of expressed sequence tags (ESTs) for Saccoglossus kowalevskii, a direct-developing hemichordate valuable for evolutionary comparisons with chordates. The 202,175 ESTs represent 163,633 arrayed clones carrying cDNAs prepared from embryonic libraries, and they assemble into 13,677 continuous sequences (contigs), leaving 10,896 singletons (excluding mitochondrial sequences). Of the contigs, 53% had significant matches when BLAST was used to query the NCBI databases (< or = 10(-10)), as did 51% of the singletons. Contigs most frequently matched sequences from amphioxus (29%), chordates (67%), and deuterostomes (87%). From the clone array, we isolated 400 full-length sequences for transcription factors and signaling proteins of use for evolutionary and developmental studies. The set includes sequences for fox, pax, tbx, hox, and other homeobox-containing factors, and for ligands and receptors of the TGFbeta, Wnt, Hh, Delta/Notch, and RTK pathways. At least 80% of key sequences have been obtained, when judged against gene lists of model organisms. The median length of these cDNAs is 2.3 kb, including 1.05 kb of 3' untranslated region (UTR). Only 30% are entirely matched by single contigs assembled from ESTs. We conclude that an EST collection based on 150,000 clones is a rich source of sequences for molecular developmental work, and that the EST approach is an efficient way to initiate comparative studies of a new organism.

Infection with the malaria parasite Plasmodium falciparum leads to widely different clinical conditions in children, ranging from mild flu-like symptoms to coma and death. Despite the immense medical implications, the genetic and molecular basis of this diversity remains largely unknown. Studies of in vitro gene expression have found few transcriptional differences between different parasite strains. Here we present a large study of in vivo expression profiles of parasites derived directly from blood samples from infected patients. The in vivo expression profiles define three distinct transcriptional states. The biological basis of these states can be interpreted by comparison with an extensive compendium of expression data in the yeast Saccharomyces cerevisiae. The three states in vivo closely resemble, first, active growth based on glycolytic metabolism, second, a starvation response accompanied by metabolism of alternative carbon sources, and third, an environmental stress response. The glycolytic state is highly similar to the known profile of the ring stage in vitro, but the other states have not been observed in vitro. The results reveal a previously unknown physiological diversity in the in vivo biology of the malaria parasite, in particular evidence for a functional mitochondrion in the asexual-stage parasite, and indicate in vivo and in vitro studies to determine how this variation may affect disease manifestations and treatment.

Positive natural selection is the force that drives the increase in prevalence of advantageous traits, and it has played a central role in our development as a species. Until recently, the study of natural selection in humans has largely been restricted to comparing individual candidate genes to theoretical expectations. The advent of genome-wide sequence and polymorphism data brings fundamental new tools to the study of natural selection. It is now possible to identify new candidates for selection and to reevaluate previous claims by comparison with empirical distributions of DNA sequence variation across the human genome and among populations. The flood of data and analytical methods, however, raises many new challenges. Here, we review approaches to detect positive natural selection, describe results from recent analyses of genome-wide data, and discuss the prospects and challenges ahead as we expand our understanding of the role of natural selection in shaping the human genome.

Only a small proportion of the mouse genome is transcribed into mature messenger RNA transcripts. There is an international collaborative effort to identify all full-length mRNA transcripts from the mouse, and to ensure that each is represented in a physical collection of clones. Here we report the manual annotation of 60,770 full-length mouse complementary DNA sequences. These are clustered into 33,409 'transcriptional units', contributing 90.1% of a newly established mouse transcriptome database. Of these transcriptional units, 4,258 are new protein-coding and 11,665 are new non-coding messages, indicating that non-coding RNA is a major component of the transcriptome. 41% of all transcriptional units showed evidence of alternative splicing. In protein-coding transcripts, 79% of splice variations altered the protein product. Whole-transcriptome analyses resulted in the identification of 2,431 sense-antisense pairs. The present work, completely supported by physical clones, provides the most comprehensive survey of a mammalian transcriptome so far, and is a valuable resource for functional genomics.

Type 2 diabetes mellitus is a heterogeneous inherited disorder characterized by chronic hyperglycemia resulting from pancreatic beta-cell dysfunction and insulin resistance. Although the pathogenic mechanisms are not fully understood, manifestation of the disease most likely requires interaction between both environmental and genetic factors. In the search for such susceptibility genes, we have performed a genomewide scan in 58 multiplex families (comprising 440 individuals, 229 of whom were affected) from the Botnia region in Finland. Initially, linkage between chromosome 12q24 and impaired insulin secretion had been reported, by Mahtani et al., in a subsample of 26 families. In the present study, we extend the initial genomewide scan to include 32 additional families, update the affectation status, and fine map regions of interest, and we try to replicate the initial stratification analysis. In our analysis of all 58 families, we identified suggestive linkage to one region, chromosome 9p13-q21 (nonparametric linkage [NPL] score 3.9; P<.0002). Regions with nominal P values <.05 include chromosomes 2p11 (NPL score 2.0 [P<.03]), 3p24-p22 (NPL score 2.2 [P<.02]), 4q32-q33 (NPL score 2.5 [P<.01]), 12q24 (NPL score 2.1 [P<.03]), 16p12-11 (NPL score 1.7 [P<.05]), and 17p12-p11 (NPL score 1.9 [P<.03]). When chromosome 12q24 was analyzed in only the 32 additional families, a nominal P value <.04 was observed. Together with data from other published genomewide scans, these findings lend support to the hypothesis that regions on chromosome 9p13-q21 and 12q24 may harbor susceptibility genes for type 2 diabetes.

Identification of the genetic bases for bipolar disorder remains a challenge for the understanding of this disease. Association between 76 candidate genes and bipolar disorder was tested by genotyping 90 single-nucleotide polymorphisms (SNPs) in these genes in 136 parent-proband trios. In this preliminary analysis, SNPs in two genes, brain-derived neurotrophic factor (BDNF) and the alpha subunit of the voltage-dependent calcium channel were associated with bipolar disorder at the P<0.05 level. In view of the large number of hypotheses tested, the two nominally positive associations were then tested in independent populations of bipolar patients and only BDNF remains a potential risk gene. In the replication samples, excess transmission of the valine allele of amino acid 66 of BDNF was observed in the direction of the original result in an additional sample of 334 parent-proband trios (T/U=108/87, P=0.066). Resequencing of 29 kb surrounding the BDNF gene identified 44 additional SNPs. Genotyping eight common SNPs identified three additional markers transmitted to bipolar probands at the P < 0.05 level. Strong LD was observed across this region and all adjacent pairwise haplotypes showed excess transmission to the bipolar proband. Analysis of these haplotypes using TRANSMIT revealed a global P value of 0.03. A single haplotype was identified that is shared by both the original dataset and the replication sample that is uniquely marked by both the rare A allele of the original SNP and a novel allele 11.5 kb 3'. Therefore, this study of 76 candidate genes has identified BDNF as a potential risk allele that will require additional study to confirm.

The optimal treatment of patients with cancer depends on establishing accurate diagnoses by using a complex combination of clinical and histopathological data. In some instances, this task is difficult or impossible because of atypical clinical presentation or histopathology. To determine whether the diagnosis of multiple common adult malignancies could be achieved purely by molecular classification, we subjected 218 tumor samples, spanning 14 common tumor types, and 90 normal tissue samples to oligonucleotide microarray gene expression analysis. The expression levels of 16,063 genes and expressed sequence tags were used to evaluate the accuracy of a multiclass classifier based on a support vector machine algorithm. Overall classification accuracy was 78%, far exceeding the accuracy of random classification (9%). Poorly differentiated cancers resulted in low-confidence predictions and could not be accurately classified according to their tissue of origin, indicating that they are molecularly distinct entities with dramatically different gene expression patterns compared with their well differentiated counterparts. Taken together, these results demonstrate the feasibility of accurate, multiclass molecular cancer classification and suggest a strategy for future clinical implementation of molecular cancer diagnostics.

We have generated a molecular taxonomy of lung carcinoma, the leading cause of cancer death in the United States and worldwide. Using oligonucleotide microarrays, we analyzed mRNA expression levels corresponding to 12,600 transcript sequences in 186 lung tumor samples, including 139 adenocarcinomas resected from the lung. Hierarchical and probabilistic clustering of expression data defined distinct subclasses of lung adenocarcinoma. Among these were tumors with high relative expression of neuroendocrine genes and of type II pneumocyte genes, respectively. Retrospective analysis revealed a less favorable outcome for the adenocarcinomas with neuroendocrine gene expression. The diagnostic potential of expression profiling is emphasized by its ability to discriminate primary lung adenocarcinomas from metastases of extra-pulmonary origin. These results suggest that integration of expression profile data with clinical parameters could aid in diagnosis of lung cancer patients.

Dendritic cells are involved in the initiation of both innate and adaptive immunity. To systematically explore how dendritic cells modulate the immune system in response to different pathogens, we used oligonucleotide microarrays to measure gene expression profiles of dendritic cells in response to Escherichia coli, Candida albicans, and influenza virus as well as to their molecular components. Both a shared core response and pathogen-specific programs of gene expression were observed upon exposure to each of these pathogens. These results reveal that dendritic cells sense diverse pathogens and elicit tailored pathogen-specific immune responses.

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 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 (LD) mapping provides a powerful method for fine-structure localization of rare disease genes, but has not yet been widely applied to common disease. We sought to design a systematic approach for LD mapping and apply it to the localization of a gene (IBD5) conferring susceptibility to Crohn disease. The key issues are: (i) to detect a significant LD signal (ii) to rigorously bound the critical region and (iii) to identify the causal genetic variant within this region. We previously mapped the IBD5 locus to a large region spanning 18 cM of chromosome 5q31 (P<10(-4)). Using dense genetic maps of microsatellite markers and single-nucleotide polymorphisms (SNPs) across the entire region, we found strong evidence of LD. We bound the region to a common haplotype spanning 250 kb that shows strong association with the disease (P< 2 x 10(-7)) and contains the cytokine gene cluster. This finding provides overwhelming evidence that a specific common haplotype of the cytokine region in 5q31 confers susceptibility to Crohn disease. However, genetic evidence alone is not sufficient to identify the causal mutation within this region, as strong LD across the region results in multiple SNPs having equivalent genetic evidence-each consistent with the expected properties of the IBD5 locus. These results have important implications for Crohn disease in particular and LD mapping in general.

In an effort to develop a genomics-based approach to the prediction of drug response, we have developed an algorithm for classification of cell line chemosensitivity based on gene expression profiles alone. Using oligonucleotide microarrays, the expression levels of 6,817 genes were measured in a panel of 60 human cancer cell lines (the NCI-60) for which the chemosensitivity profiles of thousands of chemical compounds have been determined. We sought to determine whether the gene expression signatures of untreated cells were sufficient for the prediction of chemosensitivity. Gene expression-based classifiers of sensitivity or resistance for 232 compounds were generated and then evaluated on independent sets of data. The classifiers were designed to be independent of the cells' tissue of origin. The accuracy of chemosensitivity prediction was considerably better than would be expected by chance. Eighty-eight of 232 expression-based classifiers performed accurately (with P < 0.05) on an independent test set, whereas only 12 of the 232 would be expected to do so by chance. These results suggest that at least for a subset of compounds genomic approaches to chemosensitivity prediction are feasible.

Understanding the pattern of linkage disequilibrium (LD) in the human genome is important both for successful implementation of disease-gene mapping approaches and for inferences about human demographic histories. Previous studies have examined LD between loci within single genes or confined genomic regions, which may not be representative of the genome; between loci separated by large distances, where little LD is seen; or in population groups that differ from one study to the next. We measured LD in a large set of locus pairs distributed throughout the genome, with loci within each pair separated by short distances (average 124 bp). Given current models of the history of the human population, nearly all pairs of loci at such short distances would be expected to show complete LD as a consequence of lack of recombination in the short interval. Contrary to this expectation, a significant fraction of pairs showed incomplete LD. A standard model of recombination applied to these data leads to an estimate of effective human population size of 110,000. This estimate is an order of magnitude higher than most estimates based on nucleotide diversity. The most likely explanation of this discrepancy is that gene conversion increases the apparent rate of recombination between nearby loci.

Human disease genes show enormous variation in their allelic spectra; that is, in the number and population frequency of the disease-predisposing alleles at the loci. For some genes, there are a few predominant disease alleles. For others, there is a wide range of disease alleles, each relatively rare. The allelic spectrum is important: disease genes with only a few deleterious alleles can be more readily identified and are more amenable to clinical testing. Here, we weave together strands from the human mutation and population genetics literature to provide a framework for understanding and predicting the allelic spectra of disease genes. The theory does a reasonable job for diseases where the genetic etiology is well understood. It also has bearing on the Common Disease/Common Variants (CD/CV) hypothesis, predicting that at loci where the total frequency of disease alleles is not too small, disease loci will have relatively simple spectra.

The mouse mutant ducky, a model for absence epilepsy, is characterized by spike-wave seizures and ataxia. The ducky gene was mapped previously to distal mouse chromosome 9. High-resolution genetic and physical mapping has resulted in the identification of the Cacna2d2 gene encoding the alpha2delta2 voltage-dependent calcium channel subunit. Mutations in Cacna2d2 were found to underlie the ducky phenotype in the original ducky (du) strain and in a newly identified strain (du(2J)). Both mutations are predicted to result in loss of the full-length alpha2delta2 protein. Functional analysis shows that the alpha2delta2 subunit increases the maximum conductance of the alpha1A/beta4 channel combination when coexpressed in vitro in Xenopus oocytes. The Ca(2+) channel current in acutely dissociated du/du cerebellar Purkinje cells was reduced, with no change in single-channel conductance. In contrast, no effect on Ca(2+) channel current was seen in cerebellar granule cells, results consistent with the high level of expression of the Cacna2d2 gene in Purkinje, but not granule, neurons. Our observations document the first mammalian alpha2delta mutation and complete the association of each of the major classes of voltage-dependent Ca(2+) channel subunits with a phenotype of ataxia and epilepsy in the mouse.

Genomewide linkage analysis has been extremely successful at identification of the genetic variation underlying single-gene disorders. However, linkage analysis has been less successful for common human diseases and other complex traits in which multiple genetic and environmental factors interact to influence disease risk. We hypothesized that a highly heritable complex trait, in which the contribution of environmental factors was relatively limited, might be more amenable to linkage analysis. We therefore chose to study stature (adult height), for which heritability is approximately 75%-90% (Phillips and Matheny 1990; Carmichael and McGue 1995; Preece 1996; Silventoinen et al. 2000). We reanalyzed genomewide scans from four populations for which genotype and height data were available, using a variance-components method implemented in GENEHUNTER 2.0 (Pratt et al. 2000). The populations consisted of 408 individuals in 58 families from the Botnia region of Finland, 753 individuals in 183 families from other parts of Finland, 746 individuals in 179 families from Southern Sweden, and 420 individuals in 63 families from the Saguenay-Lac-St.-Jean region of Quebec. Four regions showed evidence of linkage to stature: 6q24-25, multipoint LOD score 3.85 at marker D6S1007 in Botnia (genomewide P<.06), 7q31.3-36 (LOD 3.40 at marker D7S2195 in Sweden, P<.02), 12p11.2-q14 (LOD 3.35 at markers D12S10990-D12S398 in Finland, P<.05) and 13q32-33 (LOD 3.56 at markers D13S779-D13S797 in Finland, P<.05). In a companion article (Perola et al. 2001 [in this issue]), strong supporting evidence is obtained for linkage to the region on chromosome 7. These studies suggest that highly heritable complex traits such as stature may be genetically tractable and provide insight into the genetic architecture of complex traits.

The mechanisms by which receptor tyrosine kinases (RTKs) utilize intracellular signaling pathways to direct gene expression and cellular response remain unclear. A current question is whether different RTKs within a single cell target similar or different sets of genes. In this study we have used the ErbB receptor network to explore the relationship between RTK activation and gene expression. We profiled growth factor-stimulated signaling pathway usage and broad gene expression patterns in two human mammary tumor cell lines expressing different complements of ErbB receptors. Although the growth factors epidermal growth factor (EGF) and neuregulin (NRG) 1 similarly stimulated Erk1/2 in MDA-MB-361 cells, EGF acting through an EGF receptor/ErbB2 heterodimer preferentially stimulated protein kinase C, and NRG1beta acting through an ErbB2/ErbB3 heterodimer preferentially stimulated Akt. The two growth factors regulated partially overlapping yet distinct sets of genes in these cells. In MDA-MB-453 cells, NRG1beta acting through an ErbB2/ErbB3 heterodimer stimulated prolonged signaling of all pathways examined relative to NRG2beta acting through the same heterodimeric receptor species. Surprisingly, NRG1beta and NRG2beta also regulated partially overlapping but distinct sets of genes in these cells. These results demonstrate that the activation of different RTKs, or activation of the same RTKs with different ligands, can lead to distinct profiles of gene regulation within a single cell type. Our observations also suggest that the identity and kinetics of signaling pathway usage by RTKs may play a role in the selection of regulated genes.

A genetic association between NOTCH4 and schizophrenia has previously been proposed. Unsing all markers previously shown to be associated, we found no evidence for such in three independent family-based samples (n=519 parent-offspring trios), and a case-control sample derived from the same ethnic background as the original observation. These data strongly suggest that NOTCH4 is not a significant susceptibility allele for schizophrenia.

Duplication and deletion of the 1.4-Mb region in 17p12 that is delimited by two 24-kb low copy number repeats (CMT1A-REPs) represent frequent genomic rearrangements resulting in two common inherited peripheral neuropathies, Charcot-Marie-Tooth disease type 1A (CMT1A) and hereditary neuropathy with liability to pressure palsy (HNPP). CMT1A and HNPP exemplify a paradigm for genomic disorders wherein unique genome architectural features result in susceptibility to DNA rearrangements that cause disease. A gene within the 1.4-Mb region, PMP22, is responsible for these disorders through a gene-dosage effect in the heterozygous duplication or deletion. However, the genomic structure of the 1.4-Mb region, including other genes contained within the rearranged genomic segment, remains essentially uncharacterized. To delineate genomic structural features, investigate higher-order genomic architecture, and identify genes in this region, we constructed PAC and BAC contigs and determined the complete nucleotide sequence. This CMT1A/HNPP genomic segment contains 1,421,129 bp of DNA. A low copy number repeat (LCR) was identified, with one copy inside and two copies outside of the 1.4-Mb region. Comparison between physical and genetic maps revealed a striking difference in recombination rates between the sexes with a lower recombination frequency in males (0.67 cM/Mb) versus females (5.5 cM/Mb). Hypothetically, this low recombination frequency in males may enable a chromosomal misalignment at proximal and distal CMT1A-REPs and promote unequal crossing over, which occurs 10 times more frequently in male meiosis. In addition to three previously described genes, five new genes (TEKT3, HS3ST3B1, NPD008/CGI-148, CDRT1, and CDRT15) and 13 predicted genes were identified. Most of these predicted genes are expressed only in embryonic stages. Analyses of the genomic region adjacent to proximal CMT1A-REP indicated an evolutionary mechanism for the formation of proximal CMT1A-REP and the creation of novel genes by DNA rearrangement during primate speciation.

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.

The genetics of asthma and atopy have been difficult to determine because these diseases are genetically heterogeneous and modified by environment. The pedigrees in our study (n=86) originate in eastern central Finland (Kainuu province). According to census records, this region had only 200 households (2,000 inhabitants) in the mid sixteenth to mid seventeenth centuries. The current population of 100,000 represents the expansion of these founders within the past 400 years. Because this population is relatively homogeneous, we hypothesized that the molecular genetic mechanisms underlying asthma might also have reduced heterogeneity and therefore be easier to dissect than in mixed populations. A recent twin family study supported a strong genetic component for asthma in Finland. We carried out a genome-wide scan for susceptibility loci in asthma in the Kainuu subpopulation. We identified two regions of suggestive linkage and studied them further with higher-density mapping. We obtained evidence for linkage in a 20-cM region of chromosome 7p14-p15 for three phenotypes: asthma, a high level of immunoglobulin E (IgE; atopy) and the combination of the phenotypes. The strongest linkage was seen for high serum IgE (non-parametric linkage (NPL) score 3.9, P=0.0001), exceeding the threshold for genome-wide significance based on simulations. We also observed linkage between this locus and asthma or atopy in two independent data sets.

Although nonsteroidal antiinflammatory drugs (NSAIDs) show great promise as therapies for colon cancer, a dispute remains regarding their mechanism of action. NSAIDs are known to inhibit cyclooxygenase (COX) enzymes, which convert arachidonic acid (AA) to prostaglandins (PGs). Therefore, NSAIDs may suppress tumorigenesis by inhibiting PG synthesis. However, various experimental studies have suggested the possibility of PG-independent mechanisms. Notably, disruption of the mouse group IIA secretory phospholipase A(2) locus (Pla2g2a), a potential source of AA for COX-2, increases tumor number despite the fact that the mutation has been predicted to decrease PG production. Some authors have attempted to reconcile the results by suggesting that the level of the precursor (AA), not the products (PGs), is the critical factor. To clarify the role of AA in tumorigenesis, we have examined the effect of deleting the group IV cytosolic phospholipase A(2) (cPLA(2)) locus (Pla2g4). We report that Apc(Min/+), cPLA(2)(-/-) mice show an 83% reduction in tumor number in the small intestine compared with littermates with genotypes Apc(Min/+), cPLA(2)(+/-) and Apc(Min/+), cPLA(2)(+/+). This tumor phenotype parallels that of COX-2 knockout mice, suggesting that cPLA(2) is the predominant source of AA for COX-2 in the intestine. The protective effect of cPLA(2) deletion is thus most likely attributed to a decrease in the AA supply to COX-2 and a resultant decrease in PG synthesis. The tumorigenic effect of sPLA(2) mutations is likely to be through a completely different pathway.

The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.

We describe a map of 1.42 million single nucleotide polymorphisms (SNPs) distributed throughout the human genome, providing an average density on available sequence of one SNP every 1.9 kilobases. These SNPs were primarily discovered by two projects: The SNP Consortium and the analysis of clone overlaps by the International Human Genome Sequencing Consortium. The map integrates all publicly available SNPs with described genes and other genomic features. We estimate that 60,000 SNPs fall within exon (coding and untranslated regions), and 85% of exons are within 5 kb of the nearest SNP. Nucleotide diversity varies greatly across the genome, in a manner broadly consistent with a standard population genetic model of human history. This high-density SNP map provides a public resource for defining haplotype variation across the genome, and should help to identify biomedically important genes for diagnosis and therapy.

Dyneins are multisubunit protein complexes that couple ATPase activity with conformational changes. They are involved in the cytoplasmatic movement of organelles (cytoplasmic dyneins) and the bending of cilia and flagella (axonemal dyneins). Here we present the first complete cDNA and genomic sequences of a human axonemal dynein beta heavy chain gene, DNAH9, which maps to 17p12. The 14-kb-long cDNA is divided into 69 exons spread over 390 kb. The cDNA sequence of DNAH9 was determined using a combination of methods including 5' rapid amplification of cDNA ends, RT-PCR, and cDNA library screening. RT-PCR using nasal epithelium and testis RNA revealed several alternatively spliced transcripts. The genomic structure was determined using three overlapping BACs sequenced by the Whitehead Institute/MIT Center for Genome Research. The predicted protein, of 4486 amino acids, is highly homologous to sea urchin axonemal beta heavy chain dyneins (67% identity). It consists of an N-terminal stem and a globular C-terminus containing the four P-loops that constitute the motor domain. Lack of proper ciliary and flagellar movement characterizes primary ciliary dyskinesia (PCD), a genetically heterogeneous autosomal recessive disorder with respiratory tract infections, bronchiectasis, male subfertility, and, in 50% of cases, situs inversus (Kartagener syndrome, KS). Dyneins are excellent candidate genes for PCD and KS because in over 50% of cases the ultrastructural defects of cilia are related to the dynein complex. Genotype analysis was performed in 31 PCD families with two or more affected siblings using a highly informative dinucleotide polymorphism located in intron 26 of DNAH9. Two families with concordant inheritance of DNAH9 alleles in affected individuals were observed. A mutation search was performed in these two "candidate families," but only polymorphic variants were found. In the absence of pathogenic mutations, the DNAH9 gene has been excluded as being responsible for autosomal recessive PCD in these families., (Copyright 2001 Academic Press.)

Leigh syndrome (LS) affects 1/40,000 newborn infants in the worldwide population and is characterized by the presence of developmental delay and lactic acidosis and by a mean life expectancy variously estimated at 3-5 years. Saguenay-Lac-Saint-Jean (SLSJ) cytochrome oxidase (COX) deficiency (LS French-Canadian type [LSFC] [MIM 220111]), an autosomal recessive form of congenital lactic acidosis, presents with developmental delay and hypotonia. It is an LS variant that is found in a geographically isolated region of Quebec and that occurs in 1/2,178 live births. Patients with LSFC show a phenotype similar to that of patients with LS, but the two groups differ in clinical presentation. We studied DNA samples from 14 patients with LSFC and from their parents, representing a total of 13 families. Because of founder effects in the SLSJ region, considerable linkage disequilibrium (LD) was expected to surround the LSFC mutation. We therefore performed a genomewide screen for LD, using 290 autosomal microsatellite markers. A single marker, D2S1356, located on 2p16, showed significant (P < 10(-5)) genomewide LD. Using high-resolution genetic mapping with additional markers and four additional families with LSFC, we were able to identify a common ancestral haplotype and to limit the critical region to approximately 2 cM between D2S119 and D2S2174. COX7AR, a gene encoding a COX7a-related protein, had previously been mapped to this region. We determined the genomic structure and resequenced this gene in patients with LSFC and in controls but found no functional mutations. Although the LSFC gene remains to be elucidated, the present study demonstrates the feasibility of using a genomewide LD strategy to localize the critical region for a rare genetic disease in a founder population.

We used genome-wide expression analysis to explore how gene expression in Saccharomyces cerevisiae is remodeled in response to various changes in extracellular environment, including changes in temperature, oxidation, nutrients, pH, and osmolarity. The results demonstrate that more than half of the genome is involved in various responses to environmental change and identify the global set of genes induced and repressed by each condition. These data implicate a substantial number of previously uncharacterized genes in these responses and reveal a signature common to environmental responses that involves approximately 10% of yeast genes. The results of expression analysis with MSN2/MSN4 mutants support the model that the Msn2/Msn4 activators induce the common response to environmental change. These results provide a global description of the transcriptional response to environmental change and extend our understanding of the role of activators in effecting this response.

Coronary heart disease (CHD) is a complex disorder constituting a major health problem in Western societies. To assess the genetic background of CHD, we performed a genomewide linkage scan in two study samples from the genetically isolated population of Finland. An initial study sample consisted of family material from the northeastern part of Finland, settled by a small number of founders approximately 300 years ago. A second study sample originated from the southwestern region of Finland, settled approximately 2,000 years ago. Families were ascertained through probands exhibiting premature CHD, defined as >50% stenosis of at least two coronary arteries at a young age, as verified by coronary angiography. Both study samples and the pooled data set provided evidence for linkage in two chromosomal regions. A region on chromosome 2q21.1-22 yielded two-point LOD scores of 3.2, 1.9, and 3.7, in the affected sib-pair (ASP) analyses of the northeastern, southwestern, and pooled study samples. The corresponding multipoint maximum-likelihood scores (MLSs) for these three study samples were 2.4, 1.3, and 3.0. In addition, a region on chromosome Xq23-26 resulted in two-point LOD scores of 1.9, 3.5, and 2.9 and in multipoint MLSs of 3.4, 3.1, and 2.5, respectively. In conclusion, this study identifies two loci likely to contribute to premature CHD: one on chromosome 2q21.1-22 and another on chromosome Xq23-26.

Generating human single-nucleotide polymorphisms (SNPs) is no longer a rate-limiting step for genetic studies of disease. The number of SNPs in public databases already exceeds 200,000, and the total is expected to exceed 1,000,000 within a year. Rather, progress is limited by the inability to genotype large numbers of SNPs. Current genotyping methods are suitable for studying individual loci or at most a handful at a time. Here, we describe a method for parallel genotyping of SNPs, called single base extension-tag array on glass slides, SBE-TAGS. The principle is as follows. SNPs are genotyped by single base extension (SBE), using bifunctional primers carrying a unique sequence tag in addition to a locus-specific sequence. Because each locus has a distinct tag, the genotyping reactions can be performed in a highly multiplexed fashion, and the resulting product can then be "demultiplexed" by hybridization to the reverse complements of the sequence tags arrayed on a glass slide. SBE-TAGS is simple and inexpensive because of the high degree of multiplexing and the use of an easily generated, generic tag array. The method is also highly accurate: we genotyped over 100 SNPs, obtaining over 5, 000 genotypes, with approximately 99% accuracy.

Most genomic variation is attributable to single nucleotide polymorphisms (SNPs), which therefore offer the highest resolution for tracking disease genes and population history. It has been proposed that a dense map of 30,000-500,000 SNPs can be used to scan the human genome for haplotypes associated with common diseases. Here we describe a simple but powerful method, called reduced representation shotgun (RRS) sequencing, for creating SNP maps. RRS re-samples specific subsets of the genome from several individuals, and compares the resulting sequences using a highly accurate SNP detection algorithm. The method can be extended by alignment to available genome sequence, increasing the yield of SNPs and providing map positions. These methods are being used by The SNP Consortium, an international collaboration of academic centres, pharmaceutical companies and a private foundation, to discover and release at least 300,000 human SNPs. We have discovered 47,172 human SNPs by RRS, and in total the Consortium has identified 148,459 SNPs. More broadly, RRS facilitates the rapid, inexpensive construction of SNP maps in biomedically and agriculturally important species. SNPs discovered by RRS also offer unique advantages for large-scale genotyping.

Human cancers arise by a combination of discrete mutations and chromosomal alterations. Loss of heterozygosity (LOH) of chromosomal regions bearing mutated tumor suppressor genes is a key event in the evolution of epithelial and mesenchymal tumors. Global patterns of LOH can be understood through allelotyping of tumors with polymorphic genetic markers. Simple sequence length polymorphisms (SSLPs, or microsatellites) are reliable genetic markers for studying LOH, but only a modest number of SSLPs are used in LOH studies because the genotyping procedure is rather tedious. Here, we report the use of a highly parallel approach to genotype large numbers of single-nucleotide polymorphisms (SNPs) for LOH, in which samples are genotyped for nearly 1,500 loci by performing 24 polymerase chain reactions (PCR), pooling the resulting amplification products and hybridizing the mixture to a high-density oligonucleotide array. We characterize the results of LOH analyses on human small-cell lung cancer (SCLC) and control DNA samples by hybridization. We show that the patterns of LOH are consistent with those obtained by analysis with both SSLPs and comparative genomic hybridization (CGH), whereas amplifications rarely are detected by the SNP array. The results validate the use of SNP array hybridization for tumor studies.

Genetic association studies are viewed as problematic and plagued by irreproducibility. Many associations have been reported for type 2 diabetes, but none have been confirmed in multiple samples and with comprehensive controls. We evaluated 16 published genetic associations to type 2 diabetes and related sub-phenotypes using a family-based design to control for population stratification, and replication samples to increase power. We were able to confirm only one association, that of the common Pro12Ala polymorphism in peroxisome proliferator-activated receptor-gamma(PPARgamma) with type 2 diabetes. By analysing over 3,000 individuals, we found a modest (1.25-fold) but significant (P=0.002) increase in diabetes risk associated with the more common proline allele (85% frequency). Moreover, our results resolve a controversy about common variation in PPARgamma. An initial study found a threefold effect, but four of five subsequent publications failed to confirm the association. All six studies are consistent with the odds ratio we describe. The data implicate inherited variation in PPARgamma in the pathogenesis of type 2 diabetes. Because the risk allele occurs at such high frequency, its modest effect translates into a large population attributable risk-influencing as much as 25% of type 2 diabetes in the general population.

The most damaging change during cancer progression is the switch from a locally growing tumour to a metastatic killer. This switch is believed to involve numerous alterations that allow tumour cells to complete the complex series of events needed for metastasis. Relatively few genes have been implicated in these events. Here we use an in vivo selection scheme to select highly metastatic melanoma cells. By analysing these cells on DNA arrays, we define a pattern of gene expression that correlates with progression to a metastatic phenotype. In particular, we show enhanced expression of several genes involved in extracellular matrix assembly and of a second set of genes that regulate, either directly or indirectly, the actin-based cytoskeleton. One of these, the small GTPase RhoC, enhances metastasis when overexpressed, whereas a dominant-negative Rho inhibits metastasis. Analysis of the phenotype of cells expressing dominant-negative Rho or RhoC indicates that RhoC is important in tumour cell invasion. The genomic approach allows us to identify families of genes involved in a process, not just single genes, and can indicate which molecular and cellular events might be important in complex biological processes such as metastasis.

We describe a novel analytical approach to gene recognition based on cross-species comparison. We first undertook a comparison of orthologous genomic loci from human and mouse, studying the extent of similarity in the number, size and sequence of exons and introns. We then developed an approach for recognizing genes within such orthologous regions by first aligning the regions using an iterative global alignment system and then identifying genes based on conservation of exonic features at aligned positions in both species. The alignment and gene recognition are performed by new programs called and, respectively. performed well at exact identification of coding exons in 117 orthologous pairs tested.

The Mom1 (Modifier of Min-1) region of distal chromosome 4 was identified during a screen for polymorphic modifiers of intestinal tumorigenesis in ApcMin/+ mice. Here, we demonstrate that the Mom1AKR allele consists of two genetic components. These include the secretory phospholipase Pla2g2a, whose candidacy as a Mom1 resistance modifier has now been tested with several transgenic lines. A second region, distal to Pla2g2a, has also been identified using fine structure recombinants. Pla2g2aAKR transgenic mice demonstrate a modest resistance to tumorigenesis in the small intestine and a very robust resistance in the large intestine. Moreover, the tumor resistance in the colon of Pla2g2aAKR animals is dosage-dependent, a finding that is consistent with our observation that Pla2g2a is expressed in goblet cells. By contrast, mice carrying the distal Mom1 modifier demonstrate a modest tumor resistance that is confined to the small intestine. Thus, the phenotypes of these two modifier loci are complementary, both in their quantitative and regional effects. The additive effects and tight linkage of these modifiers may have been necessary for the initial identification of the Mom1 region.

The chronic inflammatory bowel diseases (IBDs)-Crohn disease (CD) and ulcerative colitis (UC)-are idiopathic, inflammatory disorders of the gastrointestinal tract. These conditions have a peak incidence in early adulthood and a combined prevalence of approximately 100-200/100,000. Although the etiology of IBD is multifactorial, a significant genetic contribution to disease susceptibility is implied by epidemiological data revealing a sibling risk of approximately 35-fold for CD and approximately 15-fold for UC. To elucidate the genetic basis for these disorders, we undertook a genomewide scan in 158 Canadian sib-pair families and identified three regions of suggestive linkage (3p, 5q31-33, and 6p) and one region of significant linkage to 19p13 (LOD score 4.6). Higher-density mapping in the 5q31-q33 region revealed a locus of genomewide significance (LOD score 3.9) that contributes to CD susceptibility in families with early-onset disease. Both of these genomic regions contain numerous genes that are important to the immune and inflammatory systems and that provide good targets for future candidate-gene studies.

Spondylocostal dysostosis (SD, MIM 277300) is a group of vertebral malsegmentation syndromes with reduced stature resulting from axial skeletal defects. SD is characterized by multiple hemivertebrae, rib fusions and deletions with a non-progressive kyphoscoliosis. Cases may be sporadic or familial, with both autosomal dominant and autosomal recessive modes of inheritance reported. Autosomal recessive SD maps to a 7.8-cM interval on chromosome 19q13.1-q13.3 that is homologous with a mouse region containing a gene encoding the Notch ligand delta-like 3 (Dll3). Dll3 is mutated in the X-ray-induced mouse mutant pudgy (pu), causing a variety of vertebrocostal defects similar to SD phenotypes. Here we have cloned and sequenced human DLL3 to evaluate it as a candidate gene for SD and identified mutations in three autosomal recessive SD families. Two of the mutations predict truncations within conserved extracellular domains. The third is a missense mutation in a highly conserved glycine residue of the fifth epidermal growth factor (EGF) repeat, which has revealed an important functional role for this domain. These represent the first mutations in a human Delta homologue, thus highlighting the critical role of the Notch signalling pathway and its components in patterning the mammalian axial

Single-nucleotide polymorphisms (SNPs) have been the focus of much attention in human genetics because they are extremely abundant and well-suited for automated large-scale genotyping. Human SNPs, however, are less informative than other types of genetic markers (such as simple-sequence length polymorphisms or microsatellites) and thus more loci are required for mapping traits. SNPs offer similar advantages for experimental genetic organisms such as the mouse, but they entail no loss of informativeness because bi-allelic markers are fully informative in analysing crosses between inbred strains. Here we report a large-scale analysis of SNPs in the mouse genome. We characterized the rate of nucleotide polymorphism in eight mouse strains and identified a collection of 2,848 SNPs located in 1,755 sequence-tagged sites (STSs) using high-density oligonucleotide arrays. Three-quarters of these SNPs have been mapped on the mouse genome, providing a first-generation SNP map of the mouse. We have also developed a multiplex genotyping procedure by which a genome scan can be performed with only six genotyping reactions per animal.

MYC affects normal and neoplastic cell proliferation by altering gene expression, but the precise pathways remain unclear. We used oligonucleotide microarray analysis of 6,416 genes and expressed sequence tags to determine changes in gene expression caused by activation of c-MYC in primary human fibroblasts. In these experiments, 27 genes were consistently induced, and 9 genes were repressed. The identity of the genes revealed that MYC may affect many aspects of cell physiology altered in transformed cells: cell growth, cell cycle, adhesion, and cytoskeletal organization. Identified targets possibly linked to MYC's effects on cell growth include the nucleolar proteins nucleolin and fibrillarin, as well as the eukaryotic initiation factor 5A. Among the cell cycle genes identified as targets, the G1 cyclin D2 and the cyclin-dependent kinase binding protein CksHs2 were induced whereas the cyclin-dependent kinase inhibitor p21(Cip1) was repressed. A role for MYC in regulating cell adhesion and structure is suggested by repression of genes encoding the extracellular matrix proteins fibronectin and collagen, and the cytoskeletal protein tropomyosin. A possible mechanism for MYC-mediated apoptosis was revealed by identification of the tumor necrosis factor receptor associated protein TRAP1 as a MYC target. Finally, two immunophilins, peptidyl-prolyl cis-trans isomerase F and FKBP52, the latter of which plays a role in cell division in Arabidopsis, were up-regulated by MYC. We also explored pattern-matching methods as an alternative approach for identifying MYC target genes. The genes that displayed an expression profile most similar to endogenous Myc in microarray-based expression profiling of myeloid differentiation models were highly enriched for MYC target genes.

Many valuable animal models of human disease are known and new models are continually being generated in existing inbred strains,. Some disease models are simple mendelian traits, but most have a polygenic basis. The current approach to identifying quantitative trait loci (QTLs) that underlie such traits is to localize them in crosses, construct congenic strains carrying individual QTLs, and finally map and clone the genes. This process is time-consuming and expensive, requiring the genotyping of large crosses and many generations of breeding. Here we describe a different approach in which a panel of chromosome substitution strains (CSSs) is used for QTL mapping. Each of these strains has a single chromosome from the donor strain substituting for the corresponding chromosome in the host strain. We discuss the construction, applications and advantages of CSSs compared with conventional crosses for detecting and analysing QTLs, including those that have weak phenotypic effects.

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.

One approach to sequencing a large genome is (1) to sequence a collection of nonoverlapping "seeds" chosen from a genomic library of large-insert clones [such as bacterial artificial chromosomes (BACs)] and then (2) to take successive "walking" steps by selecting and sequencing minimally overlapping clones, using information such as clone-end sequences to identify the overlaps. In this paper we analyze the strategic issues involved in using this approach. We derive formulas showing how two key factors, the initial density of seed clones and the depth of the genomic library used for walking, affect the cost and time of a sequencing project-that is, the amount of redundant sequencing and the number of steps to cover the vast majority of the genome. We also discuss a variant strategy in which a second genomic library with clones having a somewhat smaller insert size is used to close gaps. This approach can dramatically decrease the amount of redundant sequencing, without affecting the rate at which the genome is covered.