Leonid Peshkin

A method and server for predicting damaging missense mutations

To the Editor: Applications of rapidly advancing sequencing technologies exacerbate the need to interpret individual sequence variants. Sequencing of phenotyped clinical subjects will soon become a method of choice in studies of the genetic causes of Mendelian and complex diseases. New exon capture techniques will direct sequencing efforts towards the most informative and easily interpretable protein-coding fraction of the genome. Thus, the demand for computational predictions of the impact of protein sequence variants will continue to grow. Here we present a new method and the corresponding software tool, PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2/), which is different from the early tool PolyPhen1 in the set of predictive features, alignment pipeline, and the method of classification (Fig. 1a). PolyPhen-2 uses eight sequence-based and three structure-based predictive features (Supplementary Table 1) which were selected automatically by an iterative greedy algorithm (Supplementary Methods). Majority of these features involve comparison of a property of the wild-type (ancestral, normal) allele and the corresponding property of the mutant (derived, disease-causing) allele, which together define an amino acid replacement. Most informative features characterize how well the two human alleles fit into the pattern of amino acid replacements within the multiple sequence alignment of homologous proteins, how distant the protein harboring the first deviation from the human wild-type allele is from the human protein, and whether the mutant allele originated at a hypermutable site2. The alignment pipeline selects the set of homologous sequences for the analysis using a clustering algorithm and then constructs and refines their multiple alignment (Supplementary Fig. 1). The functional significance of an allele replacement is predicted from its individual features (Supplementary Figs. 2–4) by Naive Bayes classifier (Supplementary Methods). Figure 1 PolyPhen-2 pipeline and prediction accuracy. (a) Overview of the algorithm. (b) Receiver operating characteristic (ROC) curves for predictions made by PolyPhen-2 using five-fold cross-validation on HumDiv (red) and HumVar3 (light green). UniRef100 (solid ... We used two pairs of datasets to train and test PolyPhen-2. We compiled the first pair, HumDiv, from all 3,155 damaging alleles with known effects on the molecular function causing human Mendelian diseases, present in the UniProt database, together with 6,321 differences between human proteins and their closely related mammalian homologs, assumed to be non-damaging (Supplementary Methods). The second pair, HumVar3, consists of all the 13,032 human disease-causing mutations from UniProt, together with 8,946 human nsSNPs without annotated involvement in disease, which were treated as non-damaging. We found that PolyPhen-2 performance, as presented by its receiver operating characteristic curves, was consistently superior compared to PolyPhen (Fig. 1b) and it also compared favorably with the three other popular prediction tools4–6 (Fig. 1c). For a false positive rate of 20%, PolyPhen-2 achieves the rate of true positive predictions of 92% and 73% on HumDiv and HumVar, respectively (Supplementary Table 2). One reason for a lower accuracy of predictions on HumVar is that nsSNPs assumed to be non-damaging in HumVar contain a sizable fraction of mildly deleterious alleles. In contrast, most of amino acid replacements assumed non-damaging in HumDiv must be close to selective neutrality. Because alleles that are even mildly but unconditionally deleterious cannot be fixed in the evolving lineage, no method based on comparative sequence analysis is ideal for discriminating between drastically and mildly deleterious mutations, which are assigned to the opposite categories in HumVar. Another reason is that HumDiv uses an extra criterion to avoid possible erroneous annotations of damaging mutations. For a mutation, PolyPhen-2 calculates Naive Bayes posterior probability that this mutation is damaging and reports estimates of false positive (the chance that the mutation is classified as damaging when it is in fact non-damaging) and true positive (the chance that the mutation is classified as damaging when it is indeed damaging) rates. A mutation is also appraised qualitatively, as benign, possibly damaging, or probably damaging (Supplementary Methods). The user can choose between HumDiv- and HumVar-trained PolyPhen-2. Diagnostics of Mendelian diseases requires distinguishing mutations with drastic effects from all the remaining human variation, including abundant mildly deleterious alleles. Thus, HumVar-trained PolyPhen-2 should be used for this task. In contrast, HumDiv-trained PolyPhen-2 should be used for evaluating rare alleles at loci potentially involved in complex phenotypes, dense mapping of regions identified by genome-wide association studies, and analysis of natural selection from sequence data, where even mildly deleterious alleles must be treated as damaging.

Resveratrol Delays Age-Related Deterioration and Mimics Transcriptional Aspects of Dietary Restriction without Extending Life Span

A small molecule that safely mimics the ability of dietary restriction (DR) to delay age-related diseases in laboratory animals is greatly sought after. We and others have shown that resveratrol mimics effects of DR in lower organisms. In mice, we find that resveratrol induces gene expression patterns in multiple tissues that parallel those induced by DR and every-other-day feeding. Moreover, resveratrol-fed elderly mice show a marked reduction in signs of aging, including reduced albuminuria, decreased inflammation, and apoptosis in the vascular endothelium, increased aortic elasticity, greater motor coordination, reduced cataract formation, and preserved bone mineral density. However, mice fed a standard diet did not live longer when treated with resveratrol beginning at 12 months of age. Our findings indicate that resveratrol treatment has a range of beneficial effects in mice but does not increase the longevity of ad libitum-fed animals when started midlife.

Genome sequencing reveals insights into physiology and longevity of the naked mole rat

The naked mole rat (Heterocephalus glaber) is a strictly subterranean, extraordinarily long-lived eusocial mammal. Although it is the size of a mouse, its maximum lifespan exceeds 30 years, making this animal the longest-living rodent. Naked mole rats show negligible senescence, no age-related increase in mortality, and high fecundity until death. In addition to delayed ageing, they are resistant to both spontaneous cancer and experimentally induced tumorigenesis. Naked mole rats pose a challenge to the theories that link ageing, cancer and redox homeostasis. Although characterized by significant oxidative stress, the naked mole rat proteome does not show age-related susceptibility to oxidative damage or increased ubiquitination. Naked mole rats naturally reside in large colonies with a single breeding female, the ‘queen’, who suppresses the sexual maturity of her subordinates. They also live in full darkness, at low oxygen and high carbon dioxide concentrations, and are unable to sustain thermogenesis nor feel certain types of pain. Here we report the sequencing and analysis of the naked mole rat genome, which reveals unique genome features and molecular adaptations consistent with cancer resistance, poikilothermy, hairlessness and insensitivity to low oxygen, and altered visual function, circadian rythms and taste sensing. This information provides insights into the naked mole rat’s exceptional longevity and ability to live in hostile conditions, in the dark and at low oxygen. The extreme traits of the naked mole rat, together with the reported genome and transcriptome information, offer opportunities for understanding ageing and advancing other areas of biological and biomedical research.

Mapping Gene Expression in Two Xenopus Species: Evolutionary Constraints and Developmental Flexibility

Changes in gene expression are thought to be important for morphological evolution, though little is known about the nature or magnitude of the differences. Here, we examine Xenopus laevis and Xenopus tropicalis, two amphibians with very similar development, and ask how their transcriptomes compare. Despite separation for ~30-90 million years, there is strong conservation in gene expression in the vast majority of the expressed orthologs. Significant changes occur in the level of gene expression but changes in the timing of expression (heterochrony) were much less common. Differences in level were concentrated in the earliest embryonic stages. Changes in timing were prominently found in pathways that respond to selective features of the environment. We propose that different evolutionary rates across developmental stages may be explained by the stabilization of cell fate determination in the later stages.

A public resource facilitating clinical use of genomes

Rapid advances in DNA sequencing promise to enable new diagnostics and individualized therapies. Achieving personalized medicine, however, will require extensive research on highly reidentifiable, integrated datasets of genomic and health information. To assist with this, participants in the Personal Genome Project choose to forgo privacy via our institutional review board- approved “open consent” process. The contribution of public data and samples facilitates both scientific discovery and standardization of methods. We present our findings after enrollment of more than 1,800 participants, including whole-genome sequencing of 10 pilot participant genomes (the PGP-10). We introduce the Genome-Environment-Trait Evidence (GET-Evidence) system. This tool automatically processes genomes and prioritizes both published and novel variants for interpretation. In the process of reviewing the presumed healthy PGP-10 genomes, we find numerous literature references implying serious disease. Although it is sometimes impossible to rule out a late-onset effect, stringent evidence requirements can address the high rate of incidental findings. To that end we develop a peer production system for recording and organizing variant evaluations according to standard evidence guidelines, creating a public forum for reaching consensus on interpretation of clinically relevant variants. Genome analysis becomes a two-step process: using a prioritized list to record variant evaluations, then automatically sorting reviewed variants using these annotations. Genome data, health and trait information, participant samples, and variant interpretations are all shared in the public domain—we invite others to review our results using our participant samples and contribute to our interpretations. We offer our public resource and methods to further personalized medical research.

Roundup: a multi-genome repository of orthologs and evolutionary distances

SUMMARY We have created a tool for ortholog and phylogenetic profile retrieval called Roundup. Roundup is backed by a massive repository of orthologs and associated evolutionary distances that was built using the reciprocal smallest distance algorithm, an approach that has been shown to improve upon alternative approaches of ortholog detection, such as reciprocal blast. Presently, the Roundup repository contains all possible pair-wise comparisons for over 250 genomes, including 32 Eukaryotes, more than doubling the coverage of any similar resource. The orthologs are accessible through an intuitive web interface that allows searches by genome or gene identifier, presenting results as phylogenetic profiles together with gene and molecular function annotations. Results may be downloaded as phylogenetic matrices for subsequent analysis, including the construction of whole-genome phylogenies based on gene-content data. AVAILABILITY http://rodeo.med.harvard.edu/tools/roundup.

A Noncanonical Frizzled2 Pathway Regulates Epithelial-Mesenchymal Transition and Metastasis

Wnt signaling plays a critical role in embryonic development, and genetic aberrations in this network have been broadly implicated in colorectal cancer. We find that the Wnt receptor Frizzled2 (Fzd2) and its ligands Wnt5a/b are elevated in metastatic liver, lung, colon, and breast cancer cell lines and in high-grade tumors and that their expression correlates with markers of epithelial-mesenchymal transition (EMT). Pharmacologic and genetic perturbations reveal that Fzd2 drives EMT and cell migration through a previously unrecognized, noncanonical pathway that includes Fyn and Stat3. A gene signature regulated by this pathway predicts metastasis and overall survival in patients. We have developed an antibody to Fzd2 that reduces cell migration and invasion and inhibits tumor growth and metastasis in xenografts. We propose that targeting this pathway could provide benefit for patients with tumors expressing high levels of Fzd2 and Wnt5a/b.

Effects of temperature on gene expression in embryos of the coral Montastraea faveolata.

BackgroundCoral reefs are expected to be severely impacted by rising seawater temperatures associated with climate change. This study used cDNA microarrays to investigate transcriptional effects of thermal stress in embryos of the coral Montastraea faveolata. Embryos were exposed to 27.5°C, 29.0°C, and 31.5°C directly after fertilization. Differences in gene expression were measured after 12 and 48 hours.ResultsAnalysis of differentially expressed genes indicated that increased temperatures may lead to oxidative stress, apoptosis, and a structural reconfiguration of the cytoskeletal network. Metabolic processes were downregulated, and the action of histones and zinc finger-containing proteins may have played a role in the long-term regulation upon heat stress.ConclusionsEmbryos responded differently depending on exposure time and temperature level. Embryos showed expression of stress-related genes already at a temperature of 29.0°C, but seemed to be able to counteract the initial response over time. By contrast, embryos at 31.5°C displayed continuous expression of stress genes. The genes that played a role in the response to elevated temperatures consisted of both highly conserved and coral-specific genes. These genes might serve as a basis for research into coral-specific adaptations to stress responses and global climate change.

Resveratrol Prevents High Fat/Sucrose Diet-Induced Central Arterial Wall Inflammation and Stiffening in Nonhuman Primates

Central arterial wall stiffening, driven by a chronic inflammatory milieu, accompanies arterial diseases, the leading cause of cardiovascular (CV) morbidity and mortality in Western society. An increase in central arterial wall stiffening, measured as an increase in aortic pulse wave velocity (PWV), is a major risk factor for clinical CV disease events. However, no specific therapies to reduce PWV are presently available. In rhesus monkeys, a 2 year diet high in fat and sucrose (HFS) increases not only body weight and cholesterol, but also induces prominent central arterial wall stiffening and increases PWV and inflammation. The observed loss of endothelial cell integrity, lipid and macrophage infiltration, and calcification of the arterial wall were driven by genomic and proteomic signatures of oxidative stress and inflammation. Resveratrol prevented the HFS-induced arterial wall inflammation and the accompanying increase in PWV. Dietary resveratrol may hold promise as a therapy to ameliorate increases in PWV.

On the Relationship of Protein and mRNA Dynamics in Vertebrate Embryonic Development

A biochemical explanation of development from the fertilized egg to the adult requires an understanding of the proteins and RNAs expressed over time during embryogenesis. We present a comprehensive characterization of protein and mRNA dynamics across early development in Xenopus. Surprisingly, we find that most protein levels change little and duplicated genes are expressed similarly. While the correlation between protein and mRNA levels is poor, a mass action kinetics model parameterized using protein synthesis and degradation rates regresses protein dynamics to RNA dynamics, corrected for initial protein concentration. This study provides detailed data for absolute levels of ∼10,000 proteins and ∼28,000 transcripts via a convenient web portal, a rich resource for developmental biologists. It underscores the lasting impact of maternal dowry, finds surprisingly few cases where degradation alone drives a change in protein level, and highlights the importance of transcription in shaping the dynamics of the embryonic proteome.