Deeptime: a Python library for machine learning dynamical models from time series data

Generation and analysis of time-series data is relevant to many quantitative fields ranging from economics to fluid mechanics. In the physical sciences, structures such as metastable and coherent sets, slow relaxation processes, collective variables, dominant transition pathways or manifolds and channels of probability flow can be of great importance for understanding and characterizing the kinetic, thermodynamic and mechanistic properties of the system. Deeptime is a general purpose Python library offering various tools to estimate dynamical models based on time-series data including conventional linear learning methods, such as Markov state models (MSMs), Hidden Markov Models and Koopman models, as well as kernel and deep learning approaches such as VAMPnets and deep MSMs. The library is largely compatible with scikit-learn, having a range of Estimator classes for these different models, but in contrast to scikit-learn also provides deep Model classes, e.g. in the case of an MSM, which provide a multitude of analysis methods to compute interesting thermodynamic, kinetic and dynamical quantities, such as free energies, relaxation times and transition paths. The library is designed for ease of use but also easily maintainable and extensible code. In this paper we introduce the main features and structure of the deeptime software. Deeptime can be found under https://deeptime-ml.github.io/

Deep Sequencing of Candidate Genes Identified 14 Variants Associated with Smoking Abstinence in an Ethnically Diverse Sample

Despite the large public health toll of smoking, genetic studies of smoking cessation have been limited with few discoveries of risk or protective loci. We investigated common and rare variant associations with success in quitting smoking using a cohort from 8 randomized controlled trials involving 2231 participants and a total of 10,020 common and 24,147 rare variants. We identified 14 novel markers including 6 mapping to genes previously related to psychiatric and substance use disorders, 4 of which were protective (CYP2B6 (rs1175607105), HTR3B (rs1413172952; rs1204720503), rs80210037 on chr15), and 2 of which were associated with reduced cessation (PARP15 (rs2173763), SCL18A2 (rs363222)). The others mapped to areas associated with cancer including FOXP1 (rs1288980) and ZEB1 (rs7349). Network analysis identified significant canonical pathways for the serotonin receptor signaling pathway, nicotine and bupropion metabolism, and several related to tumor suppression. Two novel markers (rs6749438; rs6718083) on chr2 are flanked by genes associated with regulation of bodyweight. The identification of novel loci in this study can provide new targets of pharmacotherapy and inform efforts to develop personalized treatments based on genetic profiles

Myopathy associated BAG3 mutations lead to protein aggregation by stalling Hsp70 networks

BAG3 is a multi-domain hub that connects two classes of chaperones, small heat shock proteins (sHSPs) via two isoleucine-proline-valine (IPV) motifs and Hsp70 via a BAG domain.\ua0Mutations in either the IPV or BAG domain of BAG3 cause a dominant form of myopathy, characterized by protein aggregation in both skeletal and cardiac muscle tissues. Surprisingly, for both disease mutants, impaired chaperone binding is not sufficient to explain disease phenotypes. Recombinant mutants are correctly folded, show unaffected Hsp70 binding but are impaired in stimulating Hsp70-dependent client processing. As a consequence, the mutant BAG3 proteins become the node for a dominant gain of function causing aggregation of itself, Hsp70, Hsp70 clients and tiered interactors within the BAG3 interactome. Importantly, genetic and pharmaceutical interference with Hsp70 binding completely reverses stress-induced protein aggregation for both BAG3 mutations. Thus, the gain of function effects of BAG3 mutants act as Achilles heel of the HSP70 machinery

Cross-sectional analysis of a large cohort with X-linked Charcot-Marie-Tooth disease (CMTX1).

OBJECTIVE: To extend the phenotypic description of Charcot-Marie-Tooth disease (CMTX1) and to draw new genotype-phenotype relationships. METHODS: Mutations in GJB1 cause the main X-linked form of CMTX (CMTX1). We report cross-sectional data from 160 patients (from 120 different families, with 89 different mutations) seen at the Inherited Neuropathies Consortium centers. RESULTS: We evaluated 87 males who had a mean age of 41 years (range 10-78 years) and 73 females who had a mean age of 46 years (range 15-84 years). Sensory-motor polyneuropathy affects both sexes, more severely in males than in females, and there was a strong correlation between age and disease burden in males but not in females. Compared with females, males had more severe reduction in motor and sensory neurophysiology parameters. In contrast to females, the radial nerve sensory response in older males tended to be more severely affected compared with younger males. Median and ulnar nerve motor amplitudes were also more severely affected in older males, whereas ulnar nerve motor potentials tended to be more affected in older females. Conversely, there were no statistical differences between the sexes in other features of the disease, such as problems with balance and hand dexterity. CONCLUSIONS: In the absence of a phenotypic correlation with specific GJB1 mutations, sex-specific distinctions and clinically relevant attributes need to be incorporated into the measurements for clinical trials in people with CMTX1. CLINICALTRIALSGOV IDENTIFIER: NCT01193075

Induced pseudoscalar coupling of the proton weak interaction

The induced pseudoscalar coupling $g_p$ is the least well known of the weak coupling constants of the proton's charged--current interaction. Its size is dictated by chiral symmetry arguments, and its measurement represents an important test of quantum chromodynamics at low energies. During the past decade a large body of new data relevant to the coupling $g_p$ has been accumulated. This data includes measurements of radiative and non radiative muon capture on targets ranging from hydrogen and few--nucleon systems to complex nuclei. Herein the authors review the theoretical underpinnings of $g_p$, the experimental studies of $g_p$, and the procedures and uncertainties in extracting the coupling from data. Current puzzles are highlighted and future opportunities are discussed.Comment: 58 pages, Latex, Revtex4, prepared for Reviews of Modern Physic

How a Diverse Research Ecosystem Has Generated New Rehabilitation Technologies: Review of NIDILRR’s Rehabilitation Engineering Research Centers

Over 50 million United States citizens (1 in 6 people in the US) have a developmental, acquired, or degenerative disability. The average US citizen can expect to live 20% of his or her life with a disability. Rehabilitation technologies play a major role in improving the quality of life for people with a disability, yet widespread and highly challenging needs remain. Within the US, a major effort aimed at the creation and evaluation of rehabilitation technology has been the Rehabilitation Engineering Research Centers (RERCs) sponsored by the National Institute on Disability, Independent Living, and Rehabilitation Research. As envisioned at their conception by a panel of the National Academy of Science in 1970, these centers were intended to take a “total approach to rehabilitation”, combining medicine, engineering, and related science, to improve the quality of life of individuals with a disability. Here, we review the scope, achievements, and ongoing projects of an unbiased sample of 19 currently active or recently terminated RERCs. Specifically, for each center, we briefly explain the needs it targets, summarize key historical advances, identify emerging innovations, and consider future directions. Our assessment from this review is that the RERC program indeed involves a multidisciplinary approach, with 36 professional fields involved, although 70% of research and development staff are in engineering fields, 23% in clinical fields, and only 7% in basic science fields; significantly, 11% of the professional staff have a disability related to their research. We observe that the RERC program has substantially diversified the scope of its work since the 1970’s, addressing more types of disabilities using more technologies, and, in particular, often now focusing on information technologies. RERC work also now often views users as integrated into an interdependent society through technologies that both people with and without disabilities co-use (such as the internet, wireless communication, and architecture). In addition, RERC research has evolved to view users as able at improving outcomes through learning, exercise, and plasticity (rather than being static), which can be optimally timed. We provide examples of rehabilitation technology innovation produced by the RERCs that illustrate this increasingly diversifying scope and evolving perspective. We conclude by discussing growth opportunities and possible future directions of the RERC program

Finishing a whole-genome shotgun: Release 3 of the Drosophila melanogaster euchromatic genome sequence

BACKGROUND: The Drosophila melanogaster genome was the first metazoan genome to have been sequenced by the whole-genome shotgun (WGS) method. Two issues relating to this achievement were widely debated in the genomics community: how correct is the sequence with respect to base-pair (bp) accuracy and frequency of assembly errors? And, how difficult is it to bring a WGS sequence to the accepted standard for finished sequence? We are now in a position to answer these questions. RESULTS: Our finishing process was designed to close gaps, improve sequence quality and validate the assembly. Sequence traces derived from the WGS and draft sequencing of individual bacterial artificial chromosomes (BACs) were assembled into BAC-sized segments. These segments were brought to high quality, and then joined to constitute the sequence of each chromosome arm. Overall assembly was verified by comparison to a physical map of fingerprinted BAC clones. In the current version of the 116.9 Mb euchromatic genome, called Release 3, the six euchromatic chromosome arms are represented by 13 scaffolds with a total of 37 sequence gaps. We compared Release 3 to Release 2; in autosomal regions of unique sequence, the error rate of Release 2 was one in 20,000 bp. CONCLUSIONS: The WGS strategy can efficiently produce a high-quality sequence of a metazoan genome while generating the reagents required for sequence finishing. However, the initial method of repeat assembly was flawed. The sequence we report here, Release 3, is a reliable resource for molecular genetic experimentation and computational analysis

Development and validation of a computerized expert system for evaluation of automated visual fields from the Ischemic Optic Neuropathy Decompression Trial

BACKGROUND: The objective of this report is to describe the methods used to develop and validate a computerized system to analyze Humphrey visual fields obtained from patients with non-arteritic anterior ischemic optic neuropathy (NAION) and enrolled in the Ischemic Optic Neuropathy Decompression Trial (IONDT). The IONDT was a multicenter study that included randomized and non-randomized patients with newly diagnosed NAION in the study eye. At baseline, randomized eyes had visual acuity of 20/64 or worse and non-randomized eyes had visual acuity of better than 20/64 or were associated with patients refusing randomization. Visual fields were measured before treatment using the Humphrey Field Analyzer with the 24-2 program, foveal threshold, and size III stimulus. METHODS: We used visual fields from 189 non-IONDT eyes with NAION to develop the computerized classification system. Six neuro-ophthalmologists ("expert panel") described definitions for visual field patterns defects using 19 visual fields representing a range of pattern defect types. The expert panel then used 120 visual fields, classified using these definitions, to refine the rules, generating revised definitions for 13 visual field pattern defects and 3 levels of severity. These definitions were incorporated into a rule-based computerized classification system run on Excel(® )software. The computerized classification system was used to categorize visual field defects for an additional 95 NAION visual fields, and the expert panel was asked to independently classify the new fields and subsequently whether they agreed with the computer classification. To account for test variability over time, we derived an adjustment factor from the pooled short term fluctuation. We examined change in defects with and without adjustment in visual fields of study participants who demonstrated a visual acuity decrease within 30 days of NAION onset (progressive NAION). RESULTS: Despite an agreed upon set of rules, there was not good agreement among the expert panel when their independent visual classifications were compared. A majority did concur with the computer classification for 91 of 95 visual fields. Remaining classification discrepancies could not be resolved without modifying existing definitions. Without using the adjustment factor, visual fields of 63.6% (14/22) patients with progressive NAION and no central defect, and all (7/7) patients with a paracentral defect, worsened within 30 days of NAION onset. After applying the adjustment factor, the visual fields of the same patients with no initial central defect and 5/7 of the patients with a paracentral defect were seen to worsen. CONCLUSION: The IONDT developed a rule-based computerized system that consistently defines pattern and severity of visual fields of NAION patients for use in a research setting

The genomes of two key bumblebee species with primitive eusocial organization

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats. Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits. Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation

Genetic analysis and natural history of Charcot-Marie-Tooth disease CMTX1 due to GJB1 variants

Charcot-Marie-Tooth disease (CMT) due to GJB1 variants (CMTX1) is the second most common form of CMT. It is an X-linked disorder characterised by progressive sensory and motor neuropathy with males affected more severely than females. Many reported GJB1 variants remain classified as variants of uncertain significance (VUS). In this large, international, multicentre study we prospectively collected demographic, clinical and genetic data on patients with CMT associated with GJB1 variants. Pathogenicity for each variant was defined using adapted American College of Medical Genetics criteria. Baseline and longitudinal analyses were conducted to study genotype-phenotype correlations, to calculate longitudinal change using the CMT Examination Score (CMTES), to compare males versus females, and pathogenic/likely pathogenic (P/LP) variants versus VUS. We present 387 patients from 295 families harbouring 154 variants in GJB1. Of these, 319 patients (82.4%) were deemed to have P/LP variants, 65 had VUS (16.8%) and 3 benign variants (0.8%; excluded from analysis); an increased proportion of patients with P/LP variants compared with using ClinVar's classification (74.6%). Male patients (166/319, 52.0%, P/LP only) were more severely affected at baseline. Baseline measures in patients with P/LP variants and VUS showed no significant differences, and regression analysis suggested the disease groups were near identical at baseline. Genotype-phenotype analysis suggested c.-17G>A produces the most severe phenotype of the five most common variants, and missense variants in the intracellular domain are less severe than other domains. Progression of disease was seen with increasing CMTES over time up to 8 years follow-up. Standard response mean (SRM), a measure of outcome responsiveness, peaked at 3 years with moderate responsiveness (change in CMTES (ΔCMTES) = 1.3 ± 2.6, p = 0.00016, SRM = 0.50). Males and females progressed similarly up to 8 years, but baseline regression analysis suggested that over a longer period, females progress more slowly. Progression was most pronounced for mild phenotypes (CMTES = 0-7; 3-year ΔCMTES = 2.3 ± 2.5, p = 0.001, SRM = 0.90). Enhanced variant interpretation has yielded an increased proportion of GJB1 variants classified as P/LP and will aid future variant interpretation in this gene. Baseline and longitudinal analysis of this large cohort of CMTX1 patients describes the natural history of the disease including the rate of progression; CMTES showed moderate responsiveness for the whole group at 3 years and higher responsiveness for the mild group at 3, 4 and 5 years. These results have implications for patient selection for upcoming clinical trials