Abstract P1-05-13: The metastatic breast cancer project: Translational genomics through direct patient engagement

Background: The Metastatic Breast Cancer Project is a nationwide research initiative that directly engages patients through social media and advocacy groups and seeks to empower them to share their samples and clinical information to accelerate research. Because the vast majority of patients are treated in the community setting, we sought to determine the feasibility of remotely obtaining tumor and saliva samples as well as medical records from a large cohort of metastatic breast cancer (MBC) patients who receive their care in diverse settings around the country.

Methods: In collaboration with patients and advocacy groups, we developed a website to allow MBC patients across the U.S. to participate. Enrolled patients are sent a saliva kit and asked to mail back a saliva sample, which is used to extract germline DNA. We contact participants' medical providers and obtain medical records and part of their tumor biopsy. Whole exome and transcriptome sequencing is performed on tumor and germline samples. Clinically annotated genomic data are used to identify mechanisms of response and resistance to therapies. The database will be shared widely with researchers. Study updates and discoveries are shared with participants regularly.

Results: In the first 8 months, 2285 MBC patients from all 50 states enrolled. 2163 (95%) completed the 16-question survey about their cancer, treatments, and demographic information. 1232 completed the online consent form permitting acquisition and analysis of medical records, tumor tissue, and saliva samples. 556 saliva samples have been received. Initial medical record and tumor sample requests have been made for patients who have provided saliva samples. To date, we have obtained medical records from 102 patients (93% success rate) and tumor samples from 32 patients (77% success rate). Whole exome and transcriptome sequencing has been successfully completed on initial samples received and is ongoing for additional samples.

Conclusions: Partnering directly with patients through social media and advocacy groups enables rapid identification of thousands of patients willing to share tumors, saliva, and medical records to accelerate research. This approach allows for rapid identification of patients with rare phenotypes such as extraordinary responders, who have been challenging to identify with traditional approaches. Remote acquisition of medical records, saliva samples, and tumor tissue for patients located throughout the U.S. is feasible. Genomic analysis and medical record abstraction for these patients is underway. As data is generated, a clinically annotated database will be shared widely with the research community.

Citation Format: Painter CA, Anastasio E, Krevalin M, Kim D, Larken K, Lennon N, Frank E, Winer EP, Lander ES, Golub T. The metastatic breast cancer project: Translational genomics through direct patient engagement [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr P1-05-13.

Abstract OT2-05-03: The metastatic breast cancer project: A national direct-to-patient research initiative to accelerate genomics research

Over the past decade, genomic characterization of tumors has shed enormous light on the molecular underpinnings of cancer. These discoveries have led to the development of novel therapies and preventive measures that have already revolutionized cancer care. Despite this progress, the genomics of metastatic breast cancer (MBC), one of the leading causes of cancer death in the U.S., remains poorly understood.

The challenge in studying tumor samples from patients with MBC has been that the tumors from most patients are not available for research, largely because the vast majority of patients are cared for in community settings where genomics studies are not typically conducted. To address this, we have launched a nationwide study, The Metastatic Breast Cancer Project, which seeks to empower patients to accelerate cancer research through sharing their samples and clinical information. We have developed an outreach program in collaboration with MBC advocacy organizations to connect MBC patients around the country with genomics research performed at the Broad Institute, allowing them to participate regardless of where they live.

Working with MBC patients and advocates, we designed a website (www.mbcproject.org) with an online questionnaire that allows patients with MBC to provide information about themselves and their cancer. Based on their answers, patients are offered an electronic consent form that explains the risks and benefits of the study and asks for permission to obtain a portion of their stored tumor tissue, a saliva sample, and copies of their medical records. For patients who consent, our clinical research team contacts their physicians and obtains copies of their medical records, which are reviewed to confirm eligibility. Enrolled patients are sent a saliva kit and asked to mail back a saliva sample, which is used to extract germline DNA. The clinical research team also contacts the patient's pathology department and requests a portion of the tumor to be sent to the Broad Institute for genomic analysis. Whole exome and transcriptome sequencing is performed on tumor and germline DNA. Sequencing data are linked to de-identified clinical information, and the resulting data are used to identify drivers of tumorigenesis, mechanisms of response and resistance to therapies, and diagnostic, prognostic, and therapeutic biomarkers. The database of clinically annotated genomic information will be shared with the NIH and the cancer research community. Study updates and discoveries are shared at regular intervals with all patients who complete the initial questionnaire.

This direct-to-patient approach should be particularly enabling for the identification of patients with rare phenotypes or clinical behavior. For this reason, the first cohorts being studied are patients with extraordinary responses to therapies and patients who present with de novo MBC. Additional cohorts will be added in the future, including young women with MBC and patients with drug-resistant MBC. This project seeks to establish a patient-researcher partnership to accelerate genomic discoveries and improve outcomes in MBC, and may ultimately serve as a means to build a new clinical and translational research model for all patients with cancer.

Citation Format: Wagle N, Painter CA, Ilzarbe M, Van Allen EM, Frank E, Oh C, Krevalin M, Lloyd M, Anderka K, Kryukov G, Boehm JS, Winer E, Lander ES, Golub TR. The metastatic breast cancer project: A national direct-to-patient research initiative to accelerate genomics research. [abstract]. In: Proceedings of the Thirty-Eighth Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2015 Dec 8-12; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(4 Suppl):Abstract nr OT2-05-03.

Whole-genome sequencing reveals activation-induced cytidine deaminase signatures during indolent chronic lymphocytic leukaemia evolution

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.

The oncogenic events driving indolent chronic lymphocytic leukaemia are relatively unknown. Here, the authors perform whole genome sequencing on 30 such tumours and identify recurrent mutations in IGLL5 and two activation induced cytidine deaminase signatures that are operative at different stages of CLL evolution.

SNP genotyping defines complex gene-flow boundaries among African malaria vector mosquitoes

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.

Genome sequence, comparative analysis, and population genetics of the domestic horse

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.

Breast epithelial differentiation is altered in BRCA1mut/+ carriers prior to the onset of cancer and contributes to the basal tumor phenotype

Abstract #3083

Human breast tumors are broadly divided into either luminal-like or basal-like cancers. This distinction is significant since basal-like tumors are more aggressive and afford a poor patient prognosis relative to luminal-like tumors. For reasons that are unclear, germline mutations in BRCA1 strongly predispose for poor prognosis basal-like tumors. The predisposition for basal-like tumors in BRCA1mut/+ patients could be due to (1) differences in underlying target cell populations between BRCA1mut/+ and BRCA1+/+ women or (2) differences in the genetic mutations arising within a single shared target cell type. This basic question has remained unresolved due to a lack of experimental models in which it can be addressed. We describe here a novel in vivo breast cancer system that enables the generation of tumors by introducing oncogenes into normal breast epithelium derived directly from human breast tissue. This system is unique in that it enables human-derived epithelial cells to be sorted for cell surface markers and transformed without requiring in vitro culture prior to implantation in vivo. Using this experimental system, we show that epithelial cells from BRCA1mut/+ patients give rise to tumors that exhibit multiple features of basal differentiation, in contrast to epithelial cells transformed with identical oncogenes from BRCA1+/+ patients. We show further that non-cancerous epithelial cells from BRCA1mut/+ patients already exhibit atypical differentiation even prior to the onset of cancer, in contrast to cells from BRCA1+/+ women. Remarkably, some of these differences are observable in the context of unperturbed breast tissue obtained from disease-free BRCA1mut/+ and BRCA1+/+ patients. Collectively, these findings show that the increased incidence of basal-like tumors in BRCA1mut/+ patients is a reflection of the altered differentiation of breast epithelial cells in BRCA1mut/+ patients.

Citation Information: Cancer Res 2009;69(2 Suppl):Abstract nr 3083.

cDNA sequences for transcription factors and signaling proteins of the hemichordate Saccoglossus kowalevskii: efficacy of the expressed sequence tag (EST) approach for evolutionary and developmental studies of a new organism

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.

Positive natural selection in the human lineage

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.

Genome-wide scan in Portuguese Island families identifies 5q31-5q35 as a susceptibility locus for schizophrenia and psychosis

Schizophrenia is a common psychiatric disorder with a complex genetic etiology. To understand the genetic basis of this syndrome in Portuguese Island populations, we performed a genome-wide scan of 29 families with schizophrenia, which identified a single region on 5q31-5q35 with strong linkage (NPL=3.09, P=0.0012 at D5S820). Empirical simulations set a genome-wide threshold of NPL=3.10 for significant linkage. Additional support for this locus in schizophrenia comes from higher-density mapping and mapping of 11 additional families. The combined set of 40 families had a peak NPL=3.28 (P=0.00066) at markers D5S2112-D5S820. These data and previous linkage findings from other investigators provide strong and consistent evidence for this genomic region as a susceptibility locus for schizophrenia. Exploratory analyses of a novel phenotype, psychosis, in families with schizophrenia and bipolar disorder detected evidence for linkage to the same markers as found in schizophrenia (peak NPL=3.03, P=0.0012 at D5S820), suggesting that this locus may be responsible for the psychotic symptoms observed in both diseases. Molecular Psychiatry (2004) 9, 213-218. doi:10.1038/sj.mp.4001418 Published online 30 December 2003

Family-based association study of 76 candidate genes in bipolar disorder: BDNF is a potential risk locus. Brain-derived neutrophic factor

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.

Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs

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.

Molecular classification of multiple tumor types

Using gene expression data to classify tumor types is a very promising tool in cancer diagnosis. Previous works show several pairs of tumor types can be successfully distinguished by their gene expression patterns (Golub et al. 1999, Ben-Dor et al. 2000, Alizadeh et al. 2000). However, the simultaneous classification across a heterogeneous set of tumor types has not been well studied yet. We obtained 190 samples from 14 tumor classes and generated a combined expression dataset containing 16063 genes for each of those samples. We performed multi-class classification by combining the outputs of binary classifiers. Three binary classifiers (k-nearest neighbors, weighted voting, and support vector machines) were applied in conjunction with three combination scenarios (one-vs-all, all-pairs, hierarchical partitioning). We achieved the best cross validation error rate of 18.75% and the best test error rate of 21.74% by using the one-vs-all support vector machine algorithm. The results demonstrate the feasibility of performing clinically useful classification from samples of multiple tumor types.

Genomewide search for type 2 diabetes mellitus susceptibility loci in Finnish families: the Botnia study

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.

Multiclass cancer diagnosis using tumor gene expression signatures

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.

Classification of human lung carcinomas by mRNA expression profiling reveals distinct adenocarcinoma subclasses

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.

The plasticity of dendritic cell responses to pathogens and their components

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.

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.

Genetic variation in the 5q31 cytokine gene cluster confers susceptibility to Crohn disease

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.

Chemosensitivity prediction by transcriptional profiling

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.

On the allelic spectrum of human disease

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.

Ducky mouse phenotype of epilepsy and ataxia is associated with mutations in the Cacna2d2 gene and decreased calcium channel current in cerebellar Purkinje cells

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 of stature in multiple populations reveals several regions with evidence of linkage to adult height

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.

Growth factor-specific signaling pathway stimulation and gene expression mediated by ErbB receptors

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.

The 1.4-Mb CMT1A duplication/HNPP deletion genomic region reveals unique genome architectural features and provides insights into the recent evolution of new genes

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.

Association analysis of NOTCH4 loci in schizophrenia using family and population-based controls

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.

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.

A susceptibility locus for asthma-related traits on chromosome 7 revealed by genome-wide scan in a founder population

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.

Deletion of cytosolic phospholipase A(2) suppresses Apc(Min)-induced tumorigenesis

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.

Axonemal beta heavy chain dynein DNAH9: cDNA sequence, genomic structure, and investigation of its role in primary ciliary dyskinesia

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.

A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms

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.

Initial sequencing and analysis of the human genome

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.

Remodeling of yeast genome expression in response to environmental changes

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.

Loss of genomic methylation causes p53-dependent apoptosis and epigenetic deregulation

Cytosine methylation of mammalian DNA is essential for the proper epigenetic regulation of gene expression and maintenance of genomic integrity. To define the mechanism through which demethylated cells die, and to establish a paradigm for identifying genes regulated by DNA methylation, we have generated mice with a conditional allele for the maintenance DNA methyltransferase gene Dnmt1. Cre-mediated deletion of Dnmt1 causes demethylation of cultured fibroblasts and a uniform p53-dependent cell death. Mutational inactivation of Trp53 partially rescues the demethylated fibroblasts for up to five population doublings in culture. Oligonucleotide microarray analysis showed that up to 10% of genes are aberrantly expressed in demethylated fibroblasts. Our results demonstrate that loss of Dnmt1 causes cell-type-specific changes in gene expression that impinge on several pathways, including expression of imprinted genes, cell-cycle control, growth factor/receptor signal transduction and mobilization of retroelements.

SBE-TAGS: an array-based method for efficient single-nucleotide polymorphism genotyping

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.

An SNP map of the human genome generated by reduced representation shotgun sequencing

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.

The common PPARgamma Pro12Ala polymorphism is associated with decreased risk of type 2 diabetes

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.

Loss-of-heterozygosity analysis of small-cell lung carcinomas using single-nucleotide polymorphism arrays

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.

Androgen metabolism and prostate cancer: establishing a model of genetic susceptibility

The prostate is an androgen-regulated organ, which has led to longstanding interest in the role of androgens in prostate carcinogenesis. Although evidence of a hormonal etiology for prostate cancer is strong, it is almost entirely circumstantial. Much of the problem in proving a causal relationship relates to the continued difficulties in reliably measuring human tissue-specific exposure to endogenous steroid hormones. The international and racial-ethnic variations in prostate cancer incidence, combined with the effects of migration on risk patterns, have suggested that genetic factors play a central role in determining prostate cancer risk. We are developing a polygenic model of prostate carcinogenesis, focused around a series of genes involved in androgen biosynthesis, transport and metabolism. We have begun to develop this model by utilizing sequence variants to study how polymorphic markers in two genes (SRD5A2 and AR) are related to prostate cancer risk within and between racial-ethnic groups. We are now collaborating with the Whitehead Institute/MIT, Center for Genome Research, to screen for single nucleotide polymorphisms in additional genes relevant to the androgen pathway and prostate cell growth. The model when fully developed can potentially provide a basis for targeting populations for screening interventions and for implementing primary preventive strategies.

Genomic analysis of metastasis reveals an essential role for RhoC

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.

Human and mouse gene structure: comparative analysis and application to exon prediction

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 Mom1AKR intestinal tumor resistance region consists of Pla2g2a and a locus distal to D4Mit64

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.

Sib-pair collection strategies for complex diseases

When planning an affected sib pair collection for use in a genomewide search for complex trait loci, researchers must ask: (a) Which family structures will yield the most informative pairs? and (b) Should recruitment extend beyond the index sib pair? The optimal collection strategy will depend on the trait's genetic architecture, but this is rarely known for non-Mendelian diseases. In the present report, we study the consequences of collecting only those sib pairs arising from pedigrees with a precisely specified structure as opposed to a strategy that collects all affected sib pairs at random (i.e., blind to the affection status of first-degree relatives). The former approach turns out to be risky because the power of specific pedigree structures can vary dramatically even among models producing identical observable parameters (such as population prevalence and sibling recurrence rate). In contrast, the latter approach typically involves only a modest loss of power as compared with the optimal (but unknowable) design. Further, we compare the strategy of collecting all affected sib pairs at random with the alternative of imposing some modest limitations on family structure (e.g., presence of at least one unaffected sib or parent). The latter approach generally provides some increase in power but entails additional effort to contact and phenotype relatives: the overall merit of imposing such requirements needs to be evaluated in the context of the specific disease to be studied and of the clinical and analytical resources available. In addition, these findings suggest that a further explanation for failure to replicate positive complex trait linkages lies in differences in ascertainment strategy between data sets.

Genomewide search in Canadian families with inflammatory bowel disease reveals two novel susceptibility loci

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.

Large-scale discovery and genotyping of single-nucleotide polymorphisms in the mouse

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.

Mutations in the human delta homologue, DLL3, cause axial skeletal defects in spondylocostal dysostosis

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