REVEL: An Ensemble Method for Predicting the Pathogenicity of Rare Missense Variants
AMERICAN JOURNAL OF HUMAN GENETICS
2016; 99 (4): 877-885
The vast majority of coding variants are rare, and assessment of the contribution of rare variants to complex traits is hampered by low statistical power and limited functional data. Improved methods for predicting the pathogenicity of rare coding variants are needed to facilitate the discovery of disease variants from exome sequencing studies. We developed REVEL (rare exome variant ensemble learner), an ensemble method for predicting the pathogenicity of missense variants on the basis of individual tools: MutPred, FATHMM, VEST, PolyPhen, SIFT, PROVEAN, MutationAssessor, MutationTaster, LRT, GERP, SiPhy, phyloP, and phastCons. REVEL was trained with recently discovered pathogenic and rare neutral missense variants, excluding those previously used to train its constituent tools. When applied to two independent test sets, REVEL had the best overall performance (p < 10(-12)) as compared to any individual tool and seven ensemble methods: MetaSVM, MetaLR, KGGSeq, Condel, CADD, DANN, and Eigen. Importantly, REVEL also had the best performance for distinguishing pathogenic from rare neutral variants with allele frequencies <0.5%. The area under the receiver operating characteristic curve (AUC) for REVEL was 0.046-0.182 higher in an independent test set of 935 recent SwissVar disease variants and 123,935 putatively neutral exome sequencing variants and 0.027-0.143 higher in an independent test set of 1,953 pathogenic and 2,406 benign variants recently reported in ClinVar than the AUCs for other ensemble methods. We provide pre-computed REVEL scores for all possible human missense variants to facilitate the identification of pathogenic variants in the sea of rare variants discovered as sequencing studies expand in scale.
View details for DOI 10.1016/j.ajhg.2016.08.016
View details for Web of Science ID 000385333700007
View details for PubMedID 27666373
Inferring transient particle transport dynamics in live cells
2015; 12 (9): 838-?
Live-cell imaging and particle tracking provide rich information on mechanisms of intracellular transport. However, trajectory analysis procedures to infer complex transport dynamics involving stochastic switching between active transport and diffusive motion are lacking. We applied Bayesian model selection to hidden Markov modeling to infer transient transport states from trajectories of mRNA-protein complexes in live mouse hippocampal neurons and metaphase kinetochores in dividing human cells. The software is available at http://hmm-bayes.org/.
View details for DOI 10.1038/NMETH.3483
View details for Web of Science ID 000360586700028
View details for PubMedID 26192083
Bayesian Approach to MSD-Based Analysis of Particle Motion in Live Cells
2012; 103 (3): 616-626
Quantitative tracking of particle motion using live-cell imaging is a powerful approach to understanding the mechanism of transport of biological molecules, organelles, and cells. However, inferring complex stochastic motion models from single-particle trajectories in an objective manner is nontrivial due to noise from sampling limitations and biological heterogeneity. Here, we present a systematic Bayesian approach to multiple-hypothesis testing of a general set of competing motion models based on particle mean-square displacements that automatically classifies particle motion, properly accounting for sampling limitations and correlated noise while appropriately penalizing model complexity according to Occam's Razor to avoid over-fitting. We test the procedure rigorously using simulated trajectories for which the underlying physical process is known, demonstrating that it chooses the simplest physical model that explains the observed data. Further, we show that computed model probabilities provide a reliability test for the downstream biological interpretation of associated parameter values. We subsequently illustrate the broad utility of the approach by applying it to disparate biological systems including experimental particle trajectories from chromosomes, kinetochores, and membrane receptors undergoing a variety of complex motions. This automated and objective Bayesian framework easily scales to large numbers of particle trajectories, making it ideal for classifying the complex motion of large numbers of single molecules and cells from high-throughput screens, as well as single-cell-, tissue-, and organism-level studies.
View details for DOI 10.1016/j.bpj.2012.06.029
View details for Web of Science ID 000307427700028
View details for PubMedID 22947879
Intracellular Transport by an Anchored Homogeneously Contracting F-Actin Meshwork
2011; 21 (7): 606-611
Actin-based contractility orchestrates changes in cell shape underlying cellular functions ranging from division to migration and wound healing. Actin also functions in intracellular transport, with the prevailing view that filamentous actin (F-actin) cables serve as tracks for motor-driven transport of cargo. We recently discovered an alternate mode of intracellular transport in starfish oocytes involving a contractile F-actin meshwork that mediates chromosome congression. The mechanisms by which this meshwork contracts and translates its contractile activity into directional transport of chromosomes remained open questions. Here, we use live-cell imaging with quantitative analysis of chromosome trajectories and meshwork velocities to show that the 3D F-actin meshwork contracts homogeneously and isotropically throughout the nuclear space. Centrifugation experiments reveal that this homogeneous contraction is translated into asymmetric, directional transport by mechanical anchoring of the meshwork to the cell cortex. Finally, by injecting inert particles of different sizes, we show that this directional transport activity is size-selective and transduced to chromosomal cargo at least in part by steric trapping or "sieving." Taken together, these results reveal mechanistic design principles of a novel and potentially versatile mode of intracellular transport based on sieving by an anchored homogeneously contracting F-actin meshwork.
View details for DOI 10.1016/j.cub.2011.03.002
View details for Web of Science ID 000289662600025
View details for PubMedID 21439825
Identification of Susceptibility Loci for Cutaneous Squamous Cell Carcinoma
JOURNAL OF INVESTIGATIVE DERMATOLOGY
2016; 136 (5): 930-937
We report a genome-wide association study of cutaneous squamous cell carcinoma conducted among non-Hispanic white members of the Kaiser Permanente Northern California health care system. The study includes a genome-wide screen of 61,457 members (6,891 cases and 54,566 controls) genotyped on the Affymetrix Axiom European array and a replication phase involving an independent set of 6,410 additional members (810 cases and 5,600 controls). Combined analysis of screening and replication phases identified 10 loci containing single-nucleotide polymorphisms (SNPs) with P-values < 5 × 10(-8). Six loci contain genes in the pigmentation pathway; SNPs at these loci appear to modulate squamous cell carcinoma risk independently of the pigmentation phenotypes. Another locus contains HLA class II genes studied in relation to elevated squamous cell carcinoma risk following immunosuppression. SNPs at the remaining three loci include an intronic SNP in FOXP1 at locus 3p13, an intergenic SNP at 3q28 near TP63, and an intergenic SNP at 9p22 near BNC2. These findings provide insights into the genetic factors accounting for inherited squamous cell carcinoma susceptibility.
View details for DOI 10.1016/j.jid.2016.01.013
View details for Web of Science ID 000375980600013
View details for PubMedID 26829030
Mapping translation 'hot-spots' in live cells by tracking single molecules of mRNA and ribosomes
Messenger RNA localization is important for cell motility by local protein translation. However, while single mRNAs can be imaged and their movements tracked in single cells, it has not yet been possible to determine whether these mRNAs are actively translating. Therefore, we imaged single β-actin mRNAs tagged with MS2 stem loops colocalizing with labeled ribosomes to determine when polysomes formed. A dataset of tracking information consisting of thousands of trajectories per cell demonstrated that mRNAs co-moving with ribosomes have significantly different diffusion properties from non-translating mRNAs that were exposed to translation inhibitors. These data indicate that ribosome load changes mRNA movement and therefore highly translating mRNAs move slower. Importantly, β-actin mRNA near focal adhesions exhibited sub-diffusive corralled movement characteristic of increased translation. This method can identify where ribosomes become engaged for local protein production and how spatial regulation of mRNA-protein interactions mediates cell directionality.
View details for DOI 10.7554/eLife.10415
View details for Web of Science ID 000373789000001
View details for PubMedID 26760529
An Arp2/3 Nucleated F-Actin Shell Fragments Nuclear Membranes at Nuclear Envelope Breakdown in Starfish Oocytes
2014; 24 (12): 1421-1428
Animal cells disassemble and reassemble their nuclear envelopes (NEs) upon each division. Nuclear envelope breakdown (NEBD) serves as a major regulatory mechanism by which mixing of cytoplasmic and nuclear compartments drives the complete reorganization of cellular architecture, committing the cell for division. Breakdown is initiated by phosphorylation-driven partial disassembly of the nuclear pore complexes (NPCs), increasing their permeability but leaving the overall NE structure intact. Subsequently, the NE is rapidly broken into membrane fragments, defining the transition from prophase to prometaphase and resulting in complete mixing of cyto- and nucleoplasm. However, the mechanism underlying this rapid NE fragmentation remains largely unknown. Here, we show that NE fragmentation during NEBD in starfish oocytes is driven by an Arp2/3 complex-nucleated F-actin "shell" that transiently polymerizes on the inner surface of the NE. Blocking the formation of this F-actin shell prevents membrane fragmentation and delays entry of large cytoplasmic molecules into the nucleus. We observe spike-like protrusions extending from the F-actin shell that appear to "pierce" the NE during the fragmentation process. Finally, we show that NE fragmentation is essential for successful reproduction, because blocking this process in meiosis leads to formation of aneuploid eggs.
View details for DOI 10.1016/j.cub.2014.05.019
View details for Web of Science ID 000337648200034
View details for PubMedID 24909322
The Kinetochore-Bound Ska1 Complex Tracks Depolymerizing Microtubules and Binds to Curved Protofilaments
2012; 23 (5): 968-980
To ensure equal chromosome segregation during mitosis, the macromolecular kinetochore must remain attached to depolymerizing microtubules, which drive chromosome movements. How kinetochores associate with depolymerizing microtubules, which undergo dramatic structural changes forming curved protofilaments, has yet to be defined in vertebrates. Here, we demonstrate that the conserved kinetochore-localized Ska1 complex tracks with depolymerizing microtubule ends and associates with both the microtubule lattice and curved protofilaments. In contrast, the Ndc80 complex, a central player in the kinetochore-microtubule interface, binds only to the straight microtubule lattice and lacks tracking activity. We demonstrate that the Ska1 complex imparts its tracking capability to the Ndc80 complex. Finally, we present a structure of the Ska1 microtubule-binding domain that reveals its interaction with microtubules and its regulation by Aurora B. This work defines an integrated kinetochore-microtubule interface formed by the Ska1 and Ndc80 complexes that associates with depolymerizing microtubules, potentially by interacting with curved microtubule protofilaments.
View details for DOI 10.1016/j.devcel.2012.09.012
View details for Web of Science ID 000311134100013
View details for PubMedID 23085020
Bayesian Approach to the Analysis of Fluorescence Correlation Spectroscopy Data II: Application to Simulated and In Vitro Data
2012; 84 (9): 3880-3888
Fluorescence correlation spectroscopy (FCS) is a powerful approach to characterizing the binding and transport dynamics of macromolecules. The unbiased interpretation of FCS data relies on the evaluation of multiple competing hypotheses to describe an underlying physical process under study, which is typically unknown a priori. Bayesian inference provides a convenient framework for this evaluation based on the temporal autocorrelation function (TACF), as previously shown theoretically using model TACF curves (He, J., Guo, S., and Bathe, M. Anal. Chem. 2012, 84). Here, we apply this procedure to simulated and experimentally measured photon-count traces analyzed using a multitau correlator, which results in complex noise properties in TACF curves that cannot be modeled easily. As a critical component of our technique, we develop two means of estimating the noise in TACF curves based either on multiple independent TACF curves themselves or a single raw underlying intensity trace, including a general procedure to ensure that independent, uncorrelated samples are used in the latter approach. Using these noise definitions, we demonstrate that the Bayesian approach selects the simplest hypothesis that describes the FCS data based on sampling and signal limitations, naturally avoiding overfitting. Further, we show that model probabilities computed using the Bayesian approach provide a reliability test for the downstream interpretation of model parameter values estimated from FCS data. Our procedure is generally applicable to FCS and image correlation spectroscopy and therefore provides an important advance in the application of these methods to the quantitative biophysical investigation of complex analytical and biological systems.
View details for DOI 10.1021/ac2034375
View details for Web of Science ID 000303349200005
View details for PubMedID 22455375
UMD-DYSF, a novel locus specific database for the compilation and interactive analysis of mutations in the dysferlin gene
2012; 33 (3): E2317-E2331
Mutations in the dysferlin gene (DYSF) lead to a complete or partial absence of the dysferlin protein in skeletal muscles and are at the origin of dysferlinopathies, a heterogeneous group of rare autosomal recessive inherited neuromuscular disorders. As a step towards a better understanding of the DYSF mutational spectrum, and towards possible inclusion of patients in future therapeutic clinical trials, we set up the Universal Mutation Database for Dysferlin (UMD-DYSF), a Locus-Specific Database developed with the UMD® software. The main objective of UMD-DYSF is to provide an updated compilation of mutational data and relevant interactive tools for the analysis of DYSF sequence variants, for diagnostic and research purposes. In particular, specific algorithms can facilitate the interpretation of newly identified intronic, missense- or isosemantic-exonic sequence variants, a problem encountered recurrently during genetic diagnosis in dysferlinopathies. UMD-DYSF v1.0 is freely accessible at www.umd.be/DYSF/. It contains a total of 742 mutational entries corresponding to 266 different disease-causing mutations identified in 558 patients worldwide diagnosed with dysferlinopathy. This article presents for the first time a comprehensive analysis of the dysferlin mutational spectrum based on all compiled DYSF disease-causing mutations reported in the literature to date, and using the main bioinformatics tools offered in UMD-DYSF.
View details for DOI 10.1002/humu.22015
View details for Web of Science ID 000300706000021
View details for PubMedID 22213072
Two-color fluorescence analysis of individual virions determines the distribution of the copy number of proteins in herpes simplex virus particles
2007; 93 (4): 1329-1337
We present a single virion method to determine absolute distributions of copy number in the protein composition of viruses and apply it to herpes simplex virus type 1. Using two-color coincidence fluorescence spectroscopy, we determine the virion-to-virion variability in copy numbers of fluorescently labeled tegument and envelope proteins relative to a capsid protein by analyzing fluorescence intensity ratios for ensembles of individual dual-labeled virions and fitting the resulting histogram of ratios. Using EYFP-tagged capsid protein VP26 as a reference for fluorescence intensity, we are able to calculate the mean and also, for the first time to our knowledge, the variation in numbers of gD, VP16, and VP22 tegument. The measurement of the number of glycoprotein D molecules was in good agreement with independent measurements of average numbers of these glycoproteins in bulk virus preparations, validating the method. The accuracy, straightforward data processing, and high throughput of this technique make it widely applicable to the analysis of the molecular composition of large complexes in general, and it is particularly suited to providing insights into virus structure, assembly, and infectivity.
View details for DOI 10.1529/biophysj.107.106351
View details for Web of Science ID 000248208800026
View details for PubMedID 17513380
Unique resistance of I/LnJ mice to a retrovirus is due to sustained interferon gamma-dependent production of virus-neutralizing antibodies
JOURNAL OF EXPERIMENTAL MEDICINE
2003; 197 (2): 233-243
Selection of immune escape variants impairs the ability of the immune system to sustain an efficient antiviral response and to control retroviral infections. Like other retroviruses, mouse mammary tumor virus (MMTV) is not efficiently eliminated by the immune system of susceptible mice. In contrast, MMTV-infected I/LnJ mice are capable of producing IgG2a virus-neutralizing antibodies, sustain this response throughout their life, and secrete antibody-coated virions into the milk, thereby preventing infection of their progeny. Antibodies were produced in response to several MMTV variants and were cross-reactive to them. Resistance to MMTV infection was recessive and was dependent on interferon (IFN)-gamma production, because I/LnJ mice with targeted deletion of the INF-gamma gene failed to produce any virus-neutralizing antibodies. These findings reveal a novel mechanism of resistance to retroviral infection that is based on a robust and sustained IFN-gamma-dependent humoral immune response.
View details for DOI 10.1084/jem.20021499
View details for Web of Science ID 000180688900010
View details for PubMedID 12538662