Analysis of the human protein interactome and comparison with yeast, worm and fly interaction datasets

We present the first analysis of the human proteome with regard to interactions between proteins. We also compare the human interactome with the available interaction datasets from yeast (Saccharomyces cerevisiae), worm (Caenorhabditis elegans) and fly (Drosophila melanogaster). Of >70,000 binary interactions, only 42 were common to human, worm and fly, and only 16 were common to all four datasets. An additional 36 interactions were common to fly and worm but were not observed in humans, although a coimmunoprecipitation assay showed that 9 of the interactions do occur in humans. A re-examination of the connectivity of essential genes in yeast and humans indicated that the available data do not support the presumption that the number of interaction partners can accurately predict whether a gene is essential. Finally, we found that proteins encoded by genes mutated in inherited genetic disorders are likely to interact with proteins known to cause similar disorders, suggesting the existence of disease subnetworks. The human interaction map constructed from our analysis should facilitate an integrative systems biology approach to elucidating the cellular networks that contribute to health and disease states.     

Automating sequence-based detection and genotyping of SNPs from diploid samples

The detection of sequence variation, for which DNA sequencing has emerged as the most sensitive and automated approach, forms the basis of all genetic analysis. Here we describe and illustrate an algorithm that accurately detects and genotypes SNPs from fluorescence-based sequence data. Because the algorithm focuses particularly on detecting SNPs through the identification of heterozygous individuals, it is especially well suited to the detection of SNPs in diploid samples obtained after DNA amplification. It is substantially more accurate than existing approaches and, notably, provides a useful quantitative measure of its confidence in each potential SNP detected and in each genotype called. Calls assigned the highest confidence are sufficiently reliable to remove the need for manual review in several contexts. For example, for sequence data from 47-90 individuals sequenced on both the forward and reverse strands, the highest-confidence calls from our algorithm detected 93% of all SNPs and 100% of high-frequency SNPs, with no false positive SNPs identified and 99.9% genotyping accuracy. This algorithm is implemented in a software package, PolyPhred version 5.0, which is freely available for academic use.     

Genome-wide genetic association of complex traits in heterogeneous stock mice

Difficulties in fine-mapping quantitative trait loci (QTLs) are a major impediment to progress in the molecular dissection of complex traits in mice. Here we show that genome-wide high-resolution mapping of multiple phenotypes can be achieved using a stock of genetically heterogeneous mice. We developed a conservative and robust bootstrap analysis to map 843 QTLs with an average 95% confidence interval of 2.8 Mb. The QTLs contribute to variation in 97 traits, including models of human disease (asthma, type 2 diabetes mellitus, obesity and anxiety) as well as immunological, biochemical and hematological phenotypes. The genetic architecture of almost all phenotypes was complex, with many loci each contributing a small proportion to the total variance. Our data set, freely available at http://gscan.well.ox.ac.uk, provides an entry point to the functional characterization of genes involved in many complex traits.     

A genome-wide scan identifies mutations in the gene encoding phosphodiesterase 11A4 (PDE11A) in individuals with adrenocortical hyperplasia

Phosphodiesterases (PDEs) regulate cyclic nucleotide levels. Increased cyclic AMP (cAMP) signaling has been associated with PRKAR1A or GNAS mutations and leads to adrenocortical tumors and Cushing syndrome. We investigated the genetic source of Cushing syndrome in individuals with adrenocortical hyperplasia that was not caused by known defects. We performed genome-wide SNP genotyping, including the adrenocortical tumor DNA. The region with the highest probability to harbor a susceptibility gene by loss of heterozygosity (LOH) and other analyses was 2q31-2q35. We identified mutations disrupting the expression of the PDE11A isoform-4 gene (PDE11A) in three kindreds. Tumor tissues showed 2q31-2q35 LOH, decreased protein expression and high cyclic nucleotide levels and cAMP-responsive element binding protein (CREB) phosphorylation. PDE11A codes for a dual-specificity PDE that is expressed in adrenal cortex and is partially inhibited by tadalafil and other PDE inhibitors; its germline inactivation is associated with adrenocortical hyperplasia, suggesting another means by which dysregulation of cAMP signaling causes endocrine tumors.     

Seborrhea-like dermatitis with psoriasiform elements caused by a mutation in ZNF750, encoding a putative C2H2 zinc finger protein

We describe an Israeli Jewish Moroccan family presenting with autosomal dominant seborrhea-like dermatosis with psoriasiform elements, including enhanced keratinocyte proliferation, parakeratosis, follicular plugging, Pityrosporum ovale overgrowth and dermal CD4 lymphocyte infiltrate. We mapped the disease gene to a 0.5-cM region overlapping the PSORS2 locus (17q25) and identified a frameshift mutation in ZNF750, which encodes a putative C2H2 zinc finger protein. ZNF750 is normally expressed in keratinocytes but not in fibroblasts and is barely detectable in CD4 lymphocytes.     

Heritable germline epimutation of MSH2 in a family with hereditary nonpolyposis colorectal cancer

Epimutations in the germline, such as methylation of the MLH1 gene, may contribute to hereditary cancer syndrome in human, but their transmission to offspring has never been documented. Here we report a family with inheritance, in three successive generations, of germline allele-specific and mosaic hypermethylation of the MSH2 gene, without evidence of DNA mismatch repair gene mutation. Three siblings carrying the germline methylation developed early-onset colorectal or endometrial cancers, all with microsatellite instability and MSH2 protein loss. Clonal bisulfite sequencing and pyrosequencing showed different methylation levels in different somatic tissues, with the highest level recorded in rectal mucosa and colon cancer tissue, and the lowest in blood leukocytes. This mosaic state of germline methylation with different tissue distribution could act as the first hit and provide a mechanism for genetic disease inheritance that may deviate from the mendelian pattern and be overlooked in conventional leukocyte-based genetic diagnosis strategy.     

Mutations in embryonic myosin heavy chain (MYH3) cause Freeman-Sheldon syndrome and Sheldon-Hall syndrome

The genetic basis of most conditions characterized by congenital contractures is largely unknown. Here we show that mutations in the embryonic myosin heavy chain (MYH3) gene cause Freeman-Sheldon syndrome (FSS), one of the most severe multiple congenital contracture (that is, arthrogryposis) syndromes, and nearly one-third of all cases of Sheldon-Hall syndrome (SHS), the most common distal arthrogryposis. FSS and SHS mutations affect different myosin residues, demonstrating that MYH3 genotype is predictive of phenotype. A structure-function analysis shows that nearly all of the MYH3 mutations are predicted to interfere with myosin's catalytic activity. These results add to the growing body of evidence showing that congenital contractures are a shared outcome of prenatal defects in myofiber force production. Elucidation of the genetic basis of these syndromes redefines congenital contractures as unique defects of the sarcomere and provides insights about what has heretofore been a poorly understood group of disorders.     

Evolution of new protein topologies through multistep gene rearrangements

New protein folds have emerged throughout evolution, but it remains unclear how a protein fold can evolve while maintaining its function, particularly when fold changes require several sequential gene rearrangements. Here, we explored hypothetical evolutionary pathways linking different topological families of the DNA-methyltransferase superfamily. These pathways entail successive gene rearrangements through a series of intermediates, all of which should be sufficiently active to maintain the organism's fitness. By means of directed evolution, and starting from HaeIII methyltransferase (M.HaeIII), we selected all the required intermediates along these paths (a duplicated fused gene and duplicates partially truncated at their 5' or 3' coding regions) that maintained function in vivo. These intermediates led to new functional genes that resembled natural methyltransferases from three known classes or that belonged to a new class first seen in our evolution experiments and subsequently identified in natural genomes. Our findings show that new protein topologies can evolve gradually through multistep gene rearrangements and provide new insights regarding these processes.