Leukaemia disease genes: large-scale cloning and pathway predictions

Retroviral insertional mutagenesis in BXH2 and AKXD recombinant inbred mice induces a high incidence of myeloid or B- and T-cell leukaemia and the proviral integration sites in the leukaemias provide powerful genetic tags for disease gene identification. Some of the disease genes identified by proviral tagging are also associated with human disease, validating this approach for human disease gene identification. Although many leukaemia disease genes have been identified over the years, many more remain to be cloned. Here we describe an inverse PCR (IPCR) method for proviral tagging that makes use of automated DNA sequencing and the genetic tools provided by the Mouse Genome Project, which increases the throughput for disease gene identification. We also use this IPCR method to clone and analyse more than 400 proviral integration sites from AKXD and BXH2 leukaemias and, in the process, identify more than 90 candidate disease genes. Some of these genes function in pathways already implicated in leukaemia, whereas others are likely to define new disease pathways. Our studies underscore the power of the mouse as a tool for gene discovery and functional genomics.     

Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease

Integration of genome-wide expression profiling with linkage analysis is a new approach to identifying genes underlying complex traits. We applied this approach to the regulation of gene expression in the BXH/HXB panel of rat recombinant inbred strains, one of the largest available rodent recombinant inbred panels and a leading resource for genetic analysis of the highly prevalent metabolic syndrome. In two tissues important to the pathogenesis of the metabolic syndrome, we mapped cis- and trans-regulatory control elements for expression of thousands of genes across the genome. Many of the most highly linked expression quantitative trait loci are regulated in cis, are inherited essentially as monogenic traits and are good candidate genes for previously mapped physiological quantitative trait loci in the rat. By comparative mapping we generated a data set of 73 candidate genes for hypertension that merit testing in human populations. Mining of this publicly available data set is expected to lead to new insights into the genes and regulatory pathways underlying the extensive range of metabolic and cardiovascular disease phenotypes that segregate in these recombinant inbred strains.     

Network modeling links breast cancer susceptibility and centrosome dysfunction

Many cancer-associated genes remain to be identified to clarify the underlying molecular mechanisms of cancer susceptibility and progression. Better understanding is also required of how mutations in cancer genes affect their products in the context of complex cellular networks. Here we have used a network modeling strategy to identify genes potentially associated with higher risk of breast cancer. Starting with four known genes encoding tumor suppressors of breast cancer, we combined gene expression profiling with functional genomic and proteomic (or 'omic') data from various species to generate a network containing 118 genes linked by 866 potential functional associations. This network shows higher connectivity than expected by chance, suggesting that its components function in biologically related pathways. One of the components of the network is HMMR, encoding a centrosome subunit, for which we demonstrate previously unknown functional associations with the breast cancer-associated gene BRCA1. Two case-control studies of incident breast cancer indicate that the HMMR locus is associated with higher risk of breast cancer in humans. Our network modeling strategy should be useful for the discovery of additional cancer-associated genes.     

Fusion genes and rearranged genes as a linear function of chromosome aberrations in cancer

Cytogenetic aberrations have been reported in 45,000 human neoplasms. Structural balanced rearrangements are associated with distinct tumor subtypes with remarkable specificity and have been essential for identifying genes involved in tumorigenesis. All balanced rearrangements that have been characterized molecularly act by deregulating a gene in one of the breakpoints or by creating a fusion gene. Because most recurrent aberrations and rearranged genes have been found in hematological disorders, whereas numerous genomic imbalances have been identified in solid tumors, it has become generally accepted that there are pathogenetic differences between these neoplasms. We here show that in every tumor type, the numbers of recurrent balanced chromosome abnormalities, fusion genes and genes rearranged as a consequence of balanced aberrations are simply a function of the number of cases with an abnormal karyotype. Hence, there may not be any fundamental tissue-specific differences in the genetic mechanisms by which neoplasia is initiated.     

Systematic mapping of genetic interactions in Caenorhabditis elegans identifies common modifiers of diverse signaling pathways

Most heritable traits, including disease susceptibility, are affected by interactions between multiple genes. However, we understand little about how genes interact because very few possible genetic interactions have been explored experimentally. We have used RNA interference in Caenorhabditis elegans to systematically test approximately 65,000 pairs of genes for their ability to interact genetically. We identify approximately 350 genetic interactions between genes functioning in signaling pathways that are mutated in human diseases, including components of the EGF/Ras, Notch and Wnt pathways. Most notably, we identify a class of highly connected 'hub' genes: inactivation of these genes can enhance the phenotypic consequences of mutation of many different genes. These hub genes all encode chromatin regulators, and their activity as genetic hubs seems to be conserved across animals. We propose that these genes function as general buffers of genetic variation and that these hub genes may act as modifier genes in multiple, mechanistically unrelated genetic diseases in humans.     

New genes involved in cancer identified by retroviral tagging

Retroviral insertional mutagenesis in BXH2 and AKXD mice induces a high incidence of myeloid leukemia and B- and T-cell lymphoma, respectively. The retroviral integration sites (RISs) in these tumors thus provide powerful genetic tags for the discovery of genes involved in cancer. Here we report the first large-scale use of retroviral tagging for cancer gene discovery in the post-genome era. Using high throughput inverse PCR, we cloned and analyzed the sequences of 884 RISs from a tumor panel composed primarily of B-cell lymphomas. We then compared these sequences, and another 415 RIS sequences previously cloned from BXH2 myeloid leukemias and from a few AKXD lymphomas, against the recently assembled mouse genome sequence. These studies identified 152 loci that are targets of retroviral integration in more than one tumor (common retroviral integration sites, CISs) and therefore likely to encode a cancer gene. Thirty-six CISs encode genes that are known or predicted to be genes involved in human cancer or their homologs, whereas others encode candidate genes that have not yet been examined for a role in human cancer. Our studies demonstrate the power of retroviral tagging for cancer gene discovery in the post-genome era and indicate a largely unrecognized complexity in mouse and presumably human cancer.     

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

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

Intérêt de la cytogénétique dans les syndromes lymphoprolifératifs chroniques

Chronic B-cell lymphoid proliferations are characterised by the presence of multiple recurrent chromosomal aberrations. A certain number of these chromosomal anomalies have been found to represent disease-specific tumour markers : the t(11;14) in mantle cell lymphoma, the t(14;18) in follicular lymphoma and the t(11;18) in lowgrade MALT lymphoma, for instance. These tumour-specific genetic markers are not only useful in patient diagnosis and follow-up, but can also have great prognostic significance. Recent prospective studies in chronic lymphocytic leukaemias have shown that in multivariate analysis, 17p and 11q deletions are the two strongest independant predictors of survival followed by age and RAI stage.The occurrence of common genetic changes in different categories of chronic lymphoid disorders is suggestive of common disease mechanisms in these apparently distinct disease entities. Therefore, it is to be expected that further multidisciplinary studies on these lymphomas would identify other potentially disease- or progression-specific chromosomal aberrations useful for accurate classification.Finally, cytogenetic studies coupled with FISH mapping have already proved instrumental in identifying hot spots of rearrangement in human cancer which in turn have served as chromosomal guides for the cloning of new genes important in oncogenesis or tumour progression.     

More than 1,000 putative new human signalling proteins revealed by EST data mining

Cloning procedures aided by homology searches of EST databases have accelerated the pace of discovery of new genes, but EST database searching remains an involved and onerous task. More than 1.6 million human EST sequences have been deposited in public databases, making it difficult to identify ESTs that represent new genes. Compounding the problems of scale are difficulties in detection associated with a high sequencing error rate and low sequence similarity between distant homologues. We have developed a new method, coupling BLAST-based searches with a domain identification protocol, that filters candidate homologues. Application of this method in a large-scale analysis of 100 signalling domain families has led to the identification of ESTs representing more than 1,000 novel human signalling genes. The 4,206 publicly available ESTs representing these genes are a valuable resource for rapid cloning of novel human signalling proteins. For example, we were able to identify ESTs of at least 106 new small GTPases, of which 6 are likely to belong to new subfamilies. In some cases, further analyses of genomic DNA led to the discovery of previously unidentified full-length protein sequences. This is exemplified by the in silico cloning (prediction of a gene product sequence using only genomic and EST sequence data) of a new type of GTPase with two catalytic domains.     

A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance

Recent genome-wide association studies have described many loci implicated in type 2 diabetes (T2D) pathophysiology and β-cell dysfunction but have contributed little to the understanding of the genetic basis of insulin resistance. We hypothesized that genes implicated in insulin resistance pathways might be uncovered by accounting for differences in body mass index (BMI) and potential interactions between BMI and genetic variants. We applied a joint meta-analysis approach to test associations with fasting insulin and glucose on a genome-wide scale. We present six previously unknown loci associated with fasting insulin at P < 5 × 10(-8) in combined discovery and follow-up analyses of 52 studies comprising up to 96,496 non-diabetic individuals. Risk variants were associated with higher triglyceride and lower high-density lipoprotein (HDL) cholesterol levels, suggesting a role for these loci in insulin resistance pathways. The discovery of these loci will aid further characterization of the role of insulin resistance in T2D pathophysiology.     

Genome-wide variation and identification of vaccine targets in the Plasmodium falciparum genome

One goal in sequencing the Plasmodium falciparum genome, the agent of the most lethal form of malaria, is to discover vaccine and drug targets. However, identifying those targets in a genome in which approximately 60% of genes have unknown functions is an enormous challenge. Because the majority of known malaria antigens and drug-resistant genes are highly polymorphic and under various selective pressures, genome-wide analysis for signatures of selection may lead to discovery of new vaccine and drug candidates. Here we surveyed 3,539 P. falciparum genes ( approximately 65% of the predicted genes) for polymorphisms and identified various highly polymorphic loci and genes, some of which encode new antigens that we confirmed using human immune sera. Our collections of genome-wide SNPs ( approximately 65% nonsynonymous) and polymorphic microsatellites and indels provide a high-resolution map (one marker per approximately 4 kb) for mapping parasite traits and studying parasite populations. In addition, we report new antigens, providing urgently needed vaccine candidates for disease control.     

Distribution and functional impact of DNA copy number variation in the rat

The abundance and dynamics of copy number variants (CNVs) in mammalian genomes poses new challenges in the identification of their impact on natural and disease phenotypes. We used computational and experimental methods to catalog CNVs in rat and found that they share important functional characteristics with those in human. In addition, 113 one-to-one orthologous genes overlap CNVs in both human and rat, 80 of which are implicated in human disease. CNVs are nonrandomly distributed throughout the genome. Chromosome 18 is a cold spot for CNVs as well as evolutionary rearrangements and segmental duplications, suggesting stringent selective mechanisms underlying CNV genesis or maintenance. By exploiting gene expression data available for rat recombinant inbred lines, we established the functional relationship of CNVs underlying 22 expression quantitative trait loci. These characteristics make the rat an excellent model for studying phenotypic effects of structural variation in relation to human complex traits and disease.     

Bayesian inference of epistatic interactions in case-control studies

Epistatic interactions among multiple genetic variants in the human genome may be important in determining individual susceptibility to common diseases. Although some existing computational methods for identifying genetic interactions have been effective for small-scale studies, we here propose a method, denoted 'bayesian epistasis association mapping' (BEAM), for genome-wide case-control studies. BEAM treats the disease-associated markers and their interactions via a bayesian partitioning model and computes, via Markov chain Monte Carlo, the posterior probability that each marker set is associated with the disease. Testing this on an age-related macular degeneration genome-wide association data set, we demonstrate that the method is significantly more powerful than existing approaches and that genome-wide case-control epistasis mapping with many thousands of markers is both computationally and statistically feasible.     

Priorities and standards in pharmacogenetic research

The current enthusiasm for pharmacogenetics draws much of its inspiration from the relatively few examples of polymorphisms that have marked and seemingly clinically relevant effects on drug response. In this regard, pharmacogenetic research has paralleled the study of human disease, which has enjoyed success in identifying mutations underlying mendelian conditions. Progress in deciphering the genetics of complex diseases, involving the interaction of multiple genes with each other and with the environment has been considerably less successful. In most instances, drug responses will probably also prove to be complex, influenced by both the environment and multiple genetic factors. For pharmacogenetics to deliver on its potential, this complexity will need to be recognized and accommodated, both in basic research and in clinical application of pharmacogenetics. As the attention of researchers begins to shift toward more systematic pharmacogenetic investigations, we suggest some priorities and standards for pharmacogenetic research.     

Global assessment of promoter methylation in a mouse model of cancer identifies ID4 as a putative tumor-suppressor gene in human leukemia

DNA methylation is associated with malignant transformation, but limitations imposed by genetic variability, tumor heterogeneity, availability of paired normal tissues and methodologies for global assessment of DNA methylation have limited progress in understanding the extent of epigenetic events in the initiation and progression of human cancer and in identifying genes that undergo methylation during cancer. We developed a mouse model of T/natural killer acute lymphoblastic leukemia that is always preceded by polyclonal lymphocyte expansion to determine how aberrant promoter DNA methylation and consequent gene silencing might be contributing to leukemic transformation. We used restriction landmark genomic scanning with this mouse model of preleukemia reproducibly progressing to leukemia to show that specific genomic methylation is associated with only the leukemic phase and is not random. We also identified Idb4 as a putative tumor-suppressor gene that is methylated in most mouse and human leukemias but in only a minority of other human cancers.     

Genome-wide association study identifies loci influencing concentrations of liver enzymes in plasma

Concentrations of liver enzymes in plasma are widely used as indicators of liver disease. We carried out a genome-wide association study in 61,089 individuals, identifying 42 loci associated with concentrations of liver enzymes in plasma, of which 32 are new associations (P = 10(-8) to P = 10(-190)). We used functional genomic approaches including metabonomic profiling and gene expression analyses to identify probable candidate genes at these regions. We identified 69 candidate genes, including genes involved in biliary transport (ATP8B1 and ABCB11), glucose, carbohydrate and lipid metabolism (FADS1, FADS2, GCKR, JMJD1C, HNF1A, MLXIPL, PNPLA3, PPP1R3B, SLC2A2 and TRIB1), glycoprotein biosynthesis and cell surface glycobiology (ABO, ASGR1, FUT2, GPLD1 and ST3GAL4), inflammation and immunity (CD276, CDH6, GCKR, HNF1A, HPR, ITGA1, RORA and STAT4) and glutathione metabolism (GSTT1, GSTT2 and GGT), as well as several genes of uncertain or unknown function (including ABHD12, EFHD1, EFNA1, EPHA2, MICAL3 and ZNF827). Our results provide new insight into genetic mechanisms and pathways influencing markers of liver function.     

Systematic meta-analyses of Alzheimer disease genetic association studies: the AlzGene database

The past decade has witnessed hundreds of reports declaring or refuting genetic association with putative Alzheimer disease susceptibility genes. This wealth of information has become increasingly difficult to follow, much less interpret. We have created a publicly available, continuously updated database that comprehensively catalogs all genetic association studies in the field of Alzheimer disease (http://www.alzgene.org). We performed systematic meta-analyses for each polymorphism with available genotype data in at least three case-control samples. In addition to identifying the epsilon4 allele of APOE and related effects, we pinpointed over a dozen potential Alzheimer disease susceptibility genes (ACE, CHRNB2, CST3, ESR1, GAPDHS, IDE, MTHFR, NCSTN, PRNP, PSEN1, TF, TFAM and TNF) with statistically significant allelic summary odds ratios (ranging from 1.11-1.38 for risk alleles and 0.92-0.67 for protective alleles). Our database provides a powerful tool for deciphering the genetics of Alzheimer disease, and it serves as a potential model for tracking the most viable gene candidates in other genetically complex diseases.     

A tiling resolution DNA microarray with complete coverage of the human genome

We constructed a tiling resolution array consisting of 32,433 overlapping BAC clones covering the entire human genome. This increases our ability to identify genetic alterations and their boundaries throughout the genome in a single comparative genomic hybridization (CGH) experiment. At this tiling resolution, we identified minute DNA alterations not previously reported. These alterations include microamplifications and deletions containing oncogenes, tumor-suppressor genes and new genes that may be associated with multiple tumor types. Our findings show the need to move beyond conventional marker-based genome comparison approaches, that rely on inference of continuity between interval markers. Our submegabase resolution tiling set for array CGH (SMRT array) allows comprehensive assessment of genomic integrity and thereby the identification of new genes associated with disease.     

Genome-wide association analysis identifies 20 loci that influence adult height

Adult height is a model polygenic trait, but there has been limited success in identifying the genes underlying its normal variation. To identify genetic variants influencing adult human height, we used genome-wide association data from 13,665 individuals and genotyped 39 variants in an additional 16,482 samples. We identified 20 variants associated with adult height (P < 5 x 10(-7), with 10 reaching P < 1 x 10(-10)). Combined, the 20 SNPs explain approximately 3% of height variation, with a approximately 5 cm difference between the 6.2% of people with 17 or fewer 'tall' alleles compared to the 5.5% with 27 or more 'tall' alleles. The loci we identified implicate genes in Hedgehog signaling (IHH, HHIP, PTCH1), extracellular matrix (EFEMP1, ADAMTSL3, ACAN) and cancer (CDK6, HMGA2, DLEU7) pathways, and provide new insights into human growth and developmental processes. Finally, our results provide insights into the genetic architecture of a classic quantitative trait.