Identification of acquired somatic mutations in the gene encoding chromatin-remodeling factor ATRX in the alpha-thalassemia myelodysplasia syndrome (ATMDS)

Inherited mutations of specific genes have elucidated the normal roles of the proteins they encode by relating specific mutations to particular phenotypes. But many potentially informative mutations in such genes are lethal early in development. Consequently, inherited mutations may not reflect all the functional roles of such proteins. Acquired, somatic defects should reflect a wider spectrum of mutations because they are not prone to negative selection in development. It has been difficult to identify such mutations so far, but microarray analysis provides a new opportunity to do so. Using this approach, we have shown that in individuals with myelodysplasia associated with alpha-thalassemia (ATMDS), somatic mutations of the gene encoding the chromatin remodeling factor ATRX cause an unexpectedly severe hematological phenotype compared with the wide spectrum of inherited mutations affecting this gene. These findings cast new light on this pleiotropic cofactor, which appears to be an essential component rather than a mere facilitator of globin gene expression.     

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.     

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.     

Extensive genomic duplication during early chordate evolution

Opinions on the hypothesis that ancient genome duplications contributed to the vertebrate genome range from strong skepticism to strong credence. Previous studies concentrated on small numbers of gene families or chromosomal regions that might not have been representative of the whole genome, or used subjective methods to identify paralogous genes and regions. Here we report a systematic and objective analysis of the draft human genome sequence to identify paralogous chromosomal regions (paralogons) formed during chordate evolution and to estimate the ages of duplicate genes. We found that the human genome contains many more paralogons than would be expected by chance. Molecular clock analysis of all protein families in humans that have orthologs in the fly and nematode indicated that a burst of gene duplication activity took place in the period 350 650 Myr ago and that many of the duplicate genes formed at this time are located within paralogons. Our results support the contention that many of the gene families in vertebrates were formed or expanded by large-scale DNA duplications in an early chordate. Considering the incompleteness of the sequence data and the antiquity of the event, the results are compatible with at least one round of polyploidy.     

Genome-wide retroviral insertional tagging of genes involved in cancer in Cdkn2a-deficient mice

We have used large-scale insertional mutagenesis to identify functional landmarks relevant to cancer in the recently completed mouse genome sequence. We infected Cdkn2a(-/-) mice with Moloney murine leukemia virus (MoMuLV) to screen for loci that can participate in tumorigenesis in collaboration with loss of the Cdkn2a-encoded tumor suppressors p16INK4a and p19ARF. Insertional mutagenesis by the latent retrovirus was synergistic with loss of Cdkn2a expression, as indicated by a marked acceleration in the development of both myeloid and lymphoid tumors. We isolated 747 unique sequences flanking retroviral integration sites and mapped them against the mouse genome sequence databases from Celera and Ensembl. In addition to 17 insertions targeting gene loci known to be cancer-related, we identified a total of 37 new common insertion sites (CISs), of which 8 encode components of signaling pathways that are involved in cancer. The effectiveness of large-scale insertional mutagenesis in a sensitized genetic background is demonstrated by the preference for activation of MAP kinase signaling, collaborating with Cdkn2a loss in generating the lymphoid and myeloid tumors. Collectively, our results show that large-scale retroviral insertional mutagenesis in genetically predisposed mice is useful both as a system for identifying genes underlying cancer and as a genetic framework for the assignment of such genes to specific oncogenic pathways.     

Systematic screen for human disease genes in yeast

High similarity between yeast and human mitochondria allows functional genomic study of Saccharomyces cerevisiae to be used to identify human genes involved in disease. So far, 102 heritable disorders have been attributed to defects in a quarter of the known nuclear-encoded mitochondrial proteins in humans. Many mitochondrial diseases remain unexplained, however, in part because only 40-60% of the presumed 700-1,000 proteins involved in mitochondrial function and biogenesis have been identified. Here we apply a systematic functional screen using the pre-existing whole-genome pool of yeast deletion mutants to identify mitochondrial proteins. Three million measurements of strain fitness identified 466 genes whose deletions impaired mitochondrial respiration, of which 265 were new. Our approach gave higher selection than other systematic approaches, including fivefold greater selection than gene expression analysis. To apply these advantages to human disorders involving mitochondria, human orthologs were identified and linked to heritable diseases using genomic map positions.     

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.     

In vivo selection using a cell-growth switch

A major obstacle to stem-cell gene therapy rests in the inability to deliver a gene into a therapeutically relevant fraction of stem cells. One way to circumvent this obstacle is to use selection. Vectors containing two linked genes serve as the basis for selection, with one gene encoding a selectable product and the other, a therapeutic protein. Applying selection in vivo has the potential to bring a minor population of genetically corrected cells into the therapeutic range. But strategies for achieving in vivo selection have traditionally relied on genes that confer resistance to cytotoxic drugs and are encumbered by toxicity. Here we describe a new system for in vivo selection that uses a 'cell-growth switch', allowing a minor population of genetically corrected cells into the therapeutic range. But strategies for achieving in vivo selection have traditionally relied on genes that confer resistance to cytotoxic drugs and are encumbered by toxicity. Here we describe a new system for in vivo selection that uses a 'cell-growth switch', allowing a minor population of genetically modified cells to be inducibly amplified, thereby averting the risks associated with cytotoxic drugs. This system provides a general platform for conditionally expanding genetically modified cell populations in vivo, and may have widespread applications in gene and cell therapy.     

Mutations in the skeletal muscle alpha-actin gene in patients with actin myopathy and nemaline myopathy.

Muscle contraction results from the force generated between the thin filament protein actin and the thick filament protein myosin, which causes the thick and thin muscle filaments to slide past each other. There are skeletal muscle, cardiac muscle, smooth muscle and non-muscle isoforms of both actin and myosin. Inherited diseases in humans have been associated with defects in cardiac actin (dilated cardiomyopathy and hypertrophic cardiomyopathy), cardiac myosin (hypertrophic cardiomyopathy) and non-muscle myosin (deafness). Here we report that mutations in the human skeletal muscle alpha-actin gene (ACTA1) are associated with two different muscle diseases, 'congenital myopathy with excess of thin myofilaments' (actin myopathy) and nemaline myopathy. Both diseases are characterized by structural abnormalities of the muscle fibres and variable degrees of muscle weakness. We have identified 15 different missense mutations resulting in 14 different amino acid changes. The missense mutations in ACTA1 are distributed throughout all six coding exons, and some involve known functional domains of actin. Approximately half of the patients died within their first year, but two female patients have survived into their thirties and have children. We identified dominant mutations in all but 1 of 14 families, with the missense mutations being single and heterozygous. The only family showing dominant inheritance comprised a 33-year-old affected mother and her two affected and two unaffected children. In another family, the clinically unaffected father is a somatic mosaic for the mutation seen in both of his affected children. We identified recessive mutations in one family in which the two affected siblings had heterozygous mutations in two different exons, one paternally and the other maternally inherited. We also identified de novo mutations in seven sporadic probands for which it was possible to analyse parental DNA.     

Sequence variation in the human angiotensin converting enzyme.

Angiotensin converting enzyme (encoded by the gene DCP1, also known as ACE) catalyses the conversion of angiotensin I to the physiologically active peptide angiotensin II, which controls fluid-electrolyte balance and systemic blood pressure. Because of its key function in the renin-angiotensin system, many association studies have been performed with DCP1. Nearly all studies have associated the presence (insertion, I) or absence (deletion, D) of a 287-bp Alu repeat element in intron 16 with the levels of circulating enzyme or cardiovascular pathophysiologies. Many epidemiological studies suggest that the DCP1*D allele confers increased susceptibility to cardiovascular disease; however, other reports have found no such association or even a beneficial effect. We present here the complete genomic sequence of DCP1 from 11 individuals, representing the longest contiguous scan (24 kb) for sequence variation in human DNA. We identified 78 varying sites in 22 chromosomes that resolved into 13 distinct haplotypes. Of the variant sites, 17 were in absolute linkage disequilibrium with the commonly typed Alu insertion/deletion polymorphism, producing two distinct and distantly related clades. We also identified a major subdivision in the Alu deletion clade that enables further analysis of the traits associated with this gene. The diversity uncovered in DCP1 is comparable to that described for other regions in the human genome. The highly correlated structure in DCP1 raises important issues for the determination of functional DNA variants within genes and genetic studies in humans based on marker association.     

The genome of the mesopolyploid crop species Brassica rapa

We report the annotation and analysis of the draft genome sequence of Brassica rapa accession Chiifu-401-42, a Chinese cabbage. We modeled 41,174 protein coding genes in the B. rapa genome, which has undergone genome triplication. We used Arabidopsis thaliana as an outgroup for investigating the consequences of genome triplication, such as structural and functional evolution. The extent of gene loss (fractionation) among triplicated genome segments varies, with one of the three copies consistently retaining a disproportionately large fraction of the genes expected to have been present in its ancestor. Variation in the number of members of gene families present in the genome may contribute to the remarkable morphological plasticity of Brassica species. The B. rapa genome sequence provides an important resource for studying the evolution of polyploid genomes and underpins the genetic improvement of Brassica oil and vegetable crops.     

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.     

The peripheral myelin protein gene PMP-22 is contained within the Charcot-Marie-Tooth disease type 1A duplication.

Charcot-Marie-Tooth disease (CMT1) is the most common form of inherited peripheral neuropathy. Although the disease is genetically heterogeneous, it has been demonstrated that the gene defect is the most frequent type (CMT1A) is the result of a partial duplication of band 17p11.2. Recent studies suggested that the peripheral hypomyelination syndrome in the trembler (Tr) mouse, a possible animal model for CMT1 disease, is associated with a point mutation in the peripheral myelin protein-22 gene (pmp-22). Expression of pmp-22 is particularly high in Schwann cells, and the protein is found in peripheral myelin. We now report that the human PMP-22 gene is contained within the CMT1A duplication. We therefore, suggest that increased dosage of the PMP-22 gene may be the cause of CMT1A neuropathy.     

Genome-wide in situ exon capture for selective resequencing

Increasingly powerful sequencing technologies are ushering in an era of personal genome sequences and raising the possibility of using such information to guide medical decisions. Genome resequencing also promises to accelerate the identification of disease-associated mutations. Roughly 98% of the human genome is composed of repeats and intergenic or non-protein-coding sequences. Thus, it is crucial to focus resequencing on high-value genomic regions. Protein-coding exons represent one such type of high-value target. We have developed a method of using flexible, high-density microarrays to capture any desired fraction of the human genome, in this case corresponding to more than 200,000 protein-coding exons. Depending on the precise protocol, up to 55-85% of the captured fragments are associated with targeted regions and up to 98% of intended exons can be recovered. This methodology provides an adaptable route toward rapid and efficient resequencing of any sizeable, non-repeat portion of the human genome.     

Screening a large reference sample to identify very low frequency sequence variants: comparisons between two genes

Most human sequence variation is in the form of single-nucleotide polymorphisms (SNPs). It has been proposed that coding-region SNPs (cSNPs) be used for direct association studies to determine the genetic basis of complex traits. The success of such studies depends on the frequency of disease-associated alleles, and their distribution in different ethnic populations. If disease-associated alleles are frequent in most populations, then direct genotyping of candidate variants could show robust associations in manageable study samples. This approach is less feasible if the genetic risk from a given candidate gene is due to many infrequent alleles. Previous studies of several genes demonstrated that most variants are relatively infrequent (<0.05). These surveys genotyped small samples (n<75) and thus had limited ability to identify rare alleles. Here we evaluate the prevalence and distribution of such rare alleles by genotyping an ethnically diverse reference sample that is more than six times larger than those used in previous studies (n=450). We screened for variants in the complete coding sequence and intron-exon junctions of two candidate genes for neuropsychiatric phenotypes: SLC6A4, encoding the serotonin transporter; and SLC18A2, encoding the vesicular monoamine transporter. Both genes have unique roles in neuronal transmission, and variants in either gene might be associated with neurobehavioral phenotypes.     

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.