Detection of regulatory variation in mouse genes

Functional polymorphism in genes can be classified as coding variation, altering the amino-acid sequence of the encoded protein, or regulatory variation, affecting the level or pattern of expression of the gene. Coding variation can be recognized directly from DNA sequence, and consequently its frequency and characteristics have been extensively described. By contrast, virtually nothing is known about the extent to which gene regulation varies in populations. Yet it is likely that regulatory variants are important in modulating gene function: alterations in gene regulation have been proposed to influence disease susceptibility and to have been the primary substrate for the evolution of species. Here, we report a systematic study to assess the extent of cis-acting regulatory variation in 69 genes across four inbred mouse strains. We find that at least four of these genes show allelic differences in expression level of 1.5-fold or greater, and that some of these differences are tissue specific. The results show that the impact of regulatory variants can be detected at a significant frequency in a genomic survey and suggest that such variation may have important consequences for organismal phenotype and evolution. The results indicate that larger-scale surveys in both mouse and human could identify a substantial number of genes with common regulatory variation.     

Karyotypic abnormalities create discordance of germline genotype and cancer cell phenotypes

The nature of mendelian inheritance assumes that all tissues in which a phenotype of interest is expressed have a uniform diploid karyotype, which is often not the case in cancer cells. Owing to nonrandom gains of chromosomes, trisomies are present in many cases of leukemia and other malignances. We used polymorphisms in the genes encoding thiopurine S-methyltransferase (TPMT), gamma-glutamyl hydrolase (GGH) and the reduced folate carrier (SLC19A1) to assess the nature of chromosomal acquisition and its influence on genotype-phenotype concordance in cancer cells. TPMT and GGH activities in somatic cells were concordant with germline genotypes, whereas activities in leukemia cells were determined by chromosomal number and whether the acquired chromosomes contained a wild-type or variant allele. Leukemia cells that had acquired an additional chromosome containing a wild-type TPMT or GGH allele had significantly lower accumulation of thioguanine nucleotides or methotrexate polyglutamates, respectively. Among these genes, there was a comparable number of acquired chromosomes with wild-type and variant alleles. Therefore, chromosomal gain can alter the concordance of germline genotype and cancer cell phenotypes, indicating that allele-specific quantitative genotyping may be required to define cancer pharmacogenomics unequivocally.     

Exome sequencing of extreme phenotypes identifies DCTN4 as a modifier of chronic Pseudomonas aeruginosa infection in cystic fibrosis

Exome sequencing has become a powerful and effective strategy for the discovery of genes underlying Mendelian disorders. However, use of exome sequencing to identify variants associated with complex traits has been more challenging, partly because the sample sizes needed for adequate power may be very large. One strategy to increase efficiency is to sequence individuals who are at both ends of a phenotype distribution (those with extreme phenotypes). Because the frequency of alleles that contribute to the trait are enriched in one or both phenotype extremes, a modest sample size can potentially be used to identify novel candidate genes and/or alleles. As part of the National Heart, Lung, and Blood Institute (NHLBI) Exome Sequencing Project (ESP), we used an extreme phenotype study design to discover that variants in DCTN4, encoding a dynactin protein, are associated with time to first P. aeruginosa airway infection, chronic P. aeruginosa infection and mucoid P. aeruginosa in individuals with cystic fibrosis.     

Genomic sequencing of colorectal adenocarcinomas identifies a recurrent VTI1A-TCF7L2 fusion

Prior studies have identified recurrent oncogenic mutations in colorectal adenocarcinoma and have surveyed exons of protein-coding genes for mutations in 11 affected individuals. Here we report whole-genome sequencing from nine individuals with colorectal cancer, including primary colorectal tumors and matched adjacent non-tumor tissues, at an average of 30.7× and 31.9× coverage, respectively. We identify an average of 75 somatic rearrangements per tumor, including complex networks of translocations between pairs of chromosomes. Eleven rearrangements encode predicted in-frame fusion proteins, including a fusion of VTI1A and TCF7L2 found in 3 out of 97 colorectal cancers. Although TCF7L2 encodes TCF4, which cooperates with β-catenin in colorectal carcinogenesis, the fusion lacks the TCF4 β-catenin-binding domain. We found a colorectal carcinoma cell line harboring the fusion gene to be dependent on VTI1A-TCF7L2 for anchorage-independent growth using RNA interference-mediated knockdown. This study shows previously unidentified levels of genomic rearrangements in colorectal carcinoma that can lead to essential gene fusions and other oncogenic events.     

TIN2, a new regulator of telomere length in human cells

Telomeres are DNA-protein structures that cap linear chromosomes and are essential for maintaining genomic stability and cell phenotype. We identified a novel human telomere-associated protein, TIN2, by interaction cloning using the telomeric DNA-binding-protein TRF1 as a bait. TIN2 interacted with TRF1 in vitro and in cells, and co-localized with TRF1 in nuclei and metaphase chromosomes. A mutant TIN2 that lacks amino-terminal sequences effects elongated human telomeres in a telomerase-dependent manner. Our findings suggest that TRF1 is insufficient for control of telomere length in human cells, and that TIN2 is an essential mediator of TRF1 function.     

Biochemical networking contributes more to genetic buffering in human and mouse metabolic pathways than does gene duplication

During evolution different genes evolve at unequal rates, reflecting the varying functional constraints on phenotype. An important contributor to this variation is genetic buffering, which reduces the potential detrimental effects of mutations. We studied whether gene duplication and redundant metabolic networks affect genetic buffering by comparing the evolutionary rate of 242 human and mouse orthologous genes involved in metabolic pathways. A gene with a redundant network is defined as one for which the structural layout of metabolic pathways provides an alternative metabolic route that can, in principle, compensate for the loss of a protein function encoded by the gene. We found that genes with redundant networks evolve at similar rates as did genes without redundant networks, [corrected] but no significant difference was detected between single-copy genes and gene families. This implies that redundancy in metabolic networks provides significantly more genetic buffering than do gene families. We also found that genes encoding proteins involved in glycolysis and gluconeogenesis showed as a group a distinct pattern of variation, in contrast to genes involved in other pathways. These results suggest that redundant networks provide genetic buffering and contribute to the functional diversification of metabolic pathways.     

Structural variants deconstruct the genome

Common genomic structural variants predispose to deleterious de novo genomic rearrangements. Understanding how they do so will require population studies across the continuum of genomic variation and ethical discussion of the nature and uses of human variation.     

Caenorhabditis elegans ABL-1 antagonizes p53-mediated germline apoptosis after ionizing irradiation

c-Abl, a conserved nonreceptor tyrosine kinase, integrates genotoxic stress responses, acting as a transducer of both pro- and antiapoptotic effector pathways. Nuclear c-Abl seems to interact with the p53 homolog p73 to elicit apoptosis. Although several observations suggest that cytoplasmic localization of c-Abl is required for antiapoptotic function, the signals that mediate its antiapoptotic effect are largely unknown. Here we show that worms carrying an abl-1 deletion allele, abl-1(ok171), are specifically hypersensitive to radiation-induced apoptosis in the Caenorhabditis elegans germ line. Our findings delineate an apoptotic pathway antagonized by ABL-1, which requires sequentially the cell cycle checkpoint genes clk-2, hus-1 and mrt-2; the C. elegans p53 homolog, cep-1; and the genes encoding the components of the conserved apoptotic machinery, ced-3, ced-9 and egl-1. ABL-1 does not antagonize germline apoptosis induced by the DNA-alkylating agent ethylnitrosourea. Furthermore, worms treated with the c-Abl inhibitor STI-571 (Gleevec; used in human cancer therapy), two newly synthesized STI-571 variants or PD166326 had a phenotype similar to that generated by abl-1(ok171). These studies indicate that ABL-1 distinguishes proapoptotic signals triggered by two different DNA-damaging agents and suggest that C. elegans might provide tissue models for development of anticancer drugs.     

Progress and prospects in rat genetics: a community view

The rat is an important system for modeling human disease. Four years ago, the rich 150-year history of rat research was transformed by the sequencing of the rat genome, ushering in an era of exceptional opportunity for identifying genes and pathways underlying disease phenotypes. Genome-wide association studies in human populations have recently provided a direct approach for finding robust genetic associations in common diseases, but identifying the precise genes and their mechanisms of action remains problematic. In the context of significant progress in rat genomic resources over the past decade, we outline achievements in rat gene discovery to date, show how these findings have been translated to human disease, and document an increasing pace of discovery of new disease genes, pathways and mechanisms. Finally, we present a set of principles that justify continuing and strengthening genetic studies in the rat model, and further development of genomic infrastructure for rat research.     

Conserved noncoding sequences are selectively constrained and not mutation cold spots

Noncoding genetic variants are likely to influence human biology and disease, but recognizing functional noncoding variants is difficult. Approximately 3% of noncoding sequence is conserved among distantly related mammals, suggesting that these evolutionarily conserved noncoding regions (CNCs) are selectively constrained and contain functional variation. However, CNCs could also merely represent regions with lower local mutation rates. Here we address this issue and show that CNCs are selectively constrained in humans by analyzing HapMap genotype data. Specifically, new (derived) alleles of SNPs within CNCs are rarer than new alleles in nonconserved regions (P = 3 x 10(-18)), indicating that evolutionary pressure has suppressed CNC-derived allele frequencies. Intronic CNCs and CNCs near genes show greater allele frequency shifts, with magnitudes comparable to those for missense variants. Thus, conserved noncoding variants are more likely to be functional. Allele frequency distributions highlight selectively constrained genomic regions that should be intensively surveyed for functionally important variation.     

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.     

Large-scale analysis of the regulatory architecture of the mouse genome with a transposon-associated sensor

We present here a Sleeping Beauty-based transposition system that offers a simple and efficient way to investigate the regulatory architecture of mammalian chromosomes in vivo. With this system, we generated several hundred mice and embryos, each with a regulatory sensor inserted at a random genomic position. This large sampling of the genome revealed the widespread presence of long-range regulatory activities along chromosomes, forming overlapping blocks with distinct tissue-specific expression potentials. The presence of tissue-restricted regulatory activities around genes with widespread expression patterns challenges the gene-centric view of genome regulation and suggests that most genes are modulated in a tissue-specific manner. The local hopping property of Sleeping Beauty provides a dynamic approach to map these regulatory domains at high resolution and, combined with Cre-mediated recombination, allows for the determination of their functions by engineering mice with specific chromosomal rearrangements.     

Widespread aneuploidy revealed by DNA microarray expression profiling

Expression profiling using DNA microarrays holds great promise for a variety of research applications, including the systematic characterization of genes discovered by sequencing projects. To demonstrate the general usefulness of this approach, we recently obtained expression profiles for nearly 300 Saccharomyces cerevisiae deletion mutants. Approximately 8% of the mutants profiled exhibited chromosome-wide expression biases, leading to spurious correlations among profiles. Competitive hybridization of genomic DNA from the mutant strains and their isogenic parental wild-type strains showed they were aneuploid for whole chromosomes or chromosomal segments. Expression profile data published by several other laboratories also suggest the use of aneuploid strains. In five separate cases, the extra chromosome harboured a close homologue of the deleted gene; in two cases, a clear growth advantage for cells acquiring the extra chromosome was demonstrated. Our results have implications for interpreting whole-genome expression data, particularly from cells known to suffer genomic instability, such as malignant or immortalized cells.     

Functional screening of an asthma QTL in YAC transgenic mice.

Many quantitative trait loci (QTLs) contributing to genetically complex conditions have been discovered, but few causative genes have been identified. This is mainly due to the large size of QTLs and the subtle connection between genotype and quantitative phenotype associated with these conditions. Transgenic mice have been successfully used to analyse well-characterized genes suspected of contributing to quantitative traits. Although this approach is powerful for examining one gene at a time, it can be impractical for surveying the large genomic intervals containing many genes that are typically associated with QTLs. To screen for genes contributing to an asthma QTL mapped to human chromosome 5q3 (refs 6,7), we characterized a panel of large-insert 5q31 transgenics based on studies demonstrating that altering gene dosage frequently affects quantitative phenotypes normally influenced by that gene. This panel of human YAC transgenics, propagating a 1-Mb interval of chromosome 5q31 containing 6 cytokine genes and 17 partially characterized genes, was screened for quantitative changes in several asthma-associated phenotypes. Multiple independent transgenic lines with altered IgE response to antigen treatment shared a 180-kb region containing 5 genes, including those encoding human interleukin 4 (IL4) and interleukin 13 (IL13 ), which induce IgE class switching in B cells. Further analysis of these mice and mice transgenic for mouse Il4 and Il13 demonstrated that moderate changes in Il4 and Il13 expression affect asthma-associated phenotypes in vivo. This functional screen of large-insert transgenics enabled us to identify genes that influence the QTL phenotype in vivo.     

Recommendations of the 2006 Human Variome Project meeting

Lists of variations in genomic DNA and their effects have been kept for some time and have been used in diagnostics and research. Although these lists have been carefully gathered and curated, there has been little standardization and coordination, complicating their use. Given the myriad possible variations in the estimated 24,000 genes in the human genome, it would be useful to have standard criteria for databases of variation. Incomplete collection and ascertainment of variants demonstrates a need for a universally accessible system. These and other problems led to the World Heath Organization-cosponsored meeting on June 20-23, 2006 in Melbourne, Australia, which launched the Human Variome Project. This meeting addressed all areas of human genetics relevant to collection of information on variation and its effects. Members of each of eight sessions (the clinic and phenotype, the diagnostic laboratory, the research laboratory, curation and collection, informatics, relevance to the emerging world, integration and federation and funding and sustainability) developed a number of recommendations that were then organized into a total of 96 recommendations to act as a foundation for future work worldwide. Here we summarize the background of the project, the meeting and its recommendations.     

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.