Uncovering precision phenotype-biomarker associations in traumatic brain injury using topological data analysis.

BACKGROUND: Traumatic brain injury (TBI) is a complex disorder that is traditionally stratified based on clinical signs and symptoms. Recent imaging and molecular biomarker innovations provide unprecedented opportunities for improved TBI precision medicine, incorporating patho-anatomical and molecular mechanisms. Complete integration of these diverse data for TBI diagnosis and patient stratification remains an unmet challenge. METHODS AND FINDINGS: The Transforming Research and Clinical Knowledge in Traumatic Brain Injury (TRACK-TBI) Pilot multicenter study enrolled 586 acute TBI patients and collected diverse common data elements (TBI-CDEs) across the study population, including imaging, genetics, and clinical outcomes. We then applied topology-based data-driven discovery to identify natural subgroups of patients, based on the TBI-CDEs collected. Our hypothesis was two-fold: 1) A machine learning tool known as topological data analysis (TDA) would reveal data-driven patterns in patient outcomes to identify candidate biomarkers of recovery, and 2) TDA-identified biomarkers would significantly predict patient outcome recovery after TBI using more traditional methods of univariate statistical tests. TDA algorithms organized and mapped the data of TBI patients in multidimensional space, identifying a subset of mild TBI patients with a specific multivariate phenotype associated with unfavorable outcome at 3 and 6 months after injury. Further analyses revealed that this patient subset had high rates of post-traumatic stress disorder (PTSD), and enrichment in several distinct genetic polymorphisms associated with cellular responses to stress and DNA damage (PARP1), and in striatal dopamine processing (ANKK1, COMT, DRD2). CONCLUSIONS: TDA identified a unique diagnostic subgroup of patients with unfavorable outcome after mild TBI that were significantly predicted by the presence of specific genetic polymorphisms. Machine learning methods such as TDA may provide a robust method for patient stratification and treatment planning targeting identified biomarkers in future clinical trials in TBI patients. TRIAL REGISTRATION: ClinicalTrials.gov Identifier NCT01565551.

Targeted mutation screening panels expose systematic population bias in detection of cystic fibrosis risk.

PURPOSE: Carrier screening for mutations contributing to cystic fibrosis (CF) is typically accomplished with panels composed of variants that are clinically validated primarily in patients of European descent. This approach has created a static genetic and phenotypic profile for CF. An opportunity now exists to reevaluate the disease profile of CFTR at a global population level. METHODS: CFTR allele and genotype frequencies were obtained from a nonpatient cohort with more than 60,000 unrelated personal genomes collected by the Exome Aggregation Consortium. Likely disease-contributing mutations were identified with the use of public database annotations and computational tools. RESULTS: We identified 131 previously described and likely pathogenic variants and another 210 untested variants with a high probability of causing protein damage. None of the current genetic screening panels or existing CFTR mutation databases covered a majority of deleterious variants in any geographical population outside of Europe. CONCLUSIONS: Both clinical annotation and mutation coverage by commercially available targeted screening panels for CF are strongly biased toward detection of reproductive risk in persons of European descent. South and East Asian populations are severely underrepresented, in part because of a definition of disease that preferences the phenotype associated with European-typical CFTR alleles.

Topological data analysis for discovery in preclinical spinal cord injury and traumatic brain injury.

Data-driven discovery in complex neurological disorders has potential to extract meaningful syndromic knowledge from large, heterogeneous data sets to enhance potential for precision medicine. Here we describe the application of topological data analysis (TDA) for data-driven discovery in preclinical traumatic brain injury (TBI) and spinal cord injury (SCI) data sets mined from the Visualized Syndromic Information and Outcomes for Neurotrauma-SCI (VISION-SCI) repository. Through direct visualization of inter-related histopathological, functional and health outcomes, TDA detected novel patterns across the syndromic network, uncovering interactions between SCI and co-occurring TBI, as well as detrimental drug effects in unpublished multicentre preclinical drug trial data in SCI. TDA also revealed that perioperative hypertension predicted long-term recovery better than any tested drug after thoracic SCI in rats. TDA-based data-driven discovery has great potential application for decision-support for basic research and clinical problems such as outcome assessment, neurocritical care, treatment planning and rapid, precision-diagnosis.

Multidimensional endotypes of asthma: topological data analysis of cross-sectional clinical, pathological, and immunological data.

BACKGROUND: Incomplete understanding of mechanisms and clinicopathobiological heterogeneity in asthma hinders research progress. Pathogenic roles for T-helper-type 17 (Th17) cells and invariant T cells implied by murine data have yet to be assessed in man. We aimed to investigate the role of Th17 and mucosal associated invariant T (MAIT) cells in airway inflammation; to characterise associations between diverse clinical and immunological features of asthma; and to identify novel multidimensional asthma endotypes. METHODS: In this single-centre, cross-sectional observational study in the UK, we assessed volunteers with mild-to-severe asthma and healthy non-atopic controls using clinical and physiological assessment and immunological sampling of blood, induced sputum, endobronchial biopsy, and bronchoalveolar lavage for flow cytometry and multiplex-electrochemiluminescence assays. Primary outcomes were changes in frequencies of Th17 and MAIT cells between health and asthma using Mann-Whitney U tests and the Jonckheere-Terpstra test (linear trend across ranked groups). The study had 80% power to detect 60% differences in T-cell frequencies at p<0·05. Bayesian Network Analysis (BNA) was used to explore associations between parameters. Topological Data Analysis (TDA) was used to identify multidimensional endotypes. The study had local research ethics approval. All participants provided informed consent. FINDINGS: Participants were 84 male and female volunteers (60 with mild-to-severe asthma and 24 healthy, non-atopic controls) aged 18-70 years recruited from clinics and research cohorts. Th17 cells and γδ17 cells were not associated with asthma, even in severe neutrophilic forms. MAIT-cell frequencies were strikingly reduced in asthma compared with health (median frequency in blood 0·9% of CD3+ cells [IQR 0·3-1·8] in asthma vs 1·6 [1·2-2·6] in health, p=0·005; in sputum 1·1 [0·7-2·0] vs 1·8 [1·6-2·3], p=0·002; and in biopsy samples 1·3 [0·7-2·3] vs 3·9% [1·3-5·3%], p=0·02), especially in severe asthma where BAL regulatory T cells were also reduced compared with those in health (4·4, 3·1-6·1, vs 8·1, 5·6-10; p=0·02). BNA and TDA identified six novel clinicopathobiological clusters of underlying disease mechanisms, with elevated mast cell mediators tryptase (p<0·0001), chymase (p=0·02), and carboxypeptidase A3 (p=0·02) in severe asthma. INTERPRETATION: This study suggests that Th17 cells do not have a major pathogenic role in human asthma. We describe a novel deficiency of MAIT cells in severe asthma. We also provide proof of concept for application of TDA to identification of multidimensional clinicopathobiological endotypes. Endotypes will require validation in further cohorts. FUNDING: Wellcome Trust.

Innate and adaptive T cells in asthmatic patients: Relationship to severity and disease mechanisms.

BACKGROUND: Asthma is a chronic inflammatory disease involving diverse cells and mediators whose interconnectivity and relationships to asthma severity are unclear. OBJECTIVE: We performed a comprehensive assessment of TH17 cells, regulatory T cells, mucosal-associated invariant T (MAIT) cells, other T-cell subsets, and granulocyte mediators in asthmatic patients. METHODS: Sixty patients with mild-to-severe asthma and 24 control subjects underwent detailed clinical assessment and provided induced sputum, endobronchial biopsy, bronchoalveolar lavage, and blood samples. Adaptive and invariant T-cell subsets, cytokines, mast cells, and basophil mediators were analyzed. RESULTS: Significant heterogeneity of T-cell phenotypes was observed, with levels of IL-13-secreting T cells and type 2 cytokines increased at some, but not all, asthma severities. TH17 cells and γδ-17 cells, proposed drivers of neutrophilic inflammation, were not strongly associated with asthma, even in severe neutrophilic forms. MAIT cell frequencies were strikingly reduced in both blood and lung tissue in relation to corticosteroid therapy and vitamin D levels, especially in patients with severe asthma in whom bronchoalveolar lavage regulatory T-cell numbers were also reduced. Bayesian network analysis identified complex relationships between pathobiologic and clinical parameters. Topological data analysis identified 6 novel clusters that are associated with diverse underlying disease mechanisms, with increased mast cell mediator levels in patients with severe asthma both in its atopic (type 2 cytokine-high) and nonatopic forms. CONCLUSION: The evidence for a role for TH17 cells in patients with severe asthma is limited. Severe asthma is associated with a striking deficiency of MAIT cells and high mast cell mediator levels. This study provides proof of concept for disease mechanistic networks in asthmatic patients with clusters that could inform the development of new therapies.

Topological data analysis of Escherichia coli O157:H7 and non-O157 survival in soils.

Shiga toxin-producing E. coli O157:H7 and non-O157 have been implicated in many foodborne illnesses caused by the consumption of contaminated fresh produce. However, data on their persistence in soils are limited due to the complexity in datasets generated from different environmental variables and bacterial taxa. There is a continuing need to distinguish the various environmental variables and different bacterial groups to understand the relationships among these factors and the pathogen survival. Using an approach called Topological Data Analysis (TDA); we reconstructed the relationship structure of E. coli O157 and non-O157 survival in 32 soils (16 organic and 16 conventionally managed soils) from California (CA) and Arizona (AZ) with a multi-resolution output. In our study, we took a community approach based on total soil microbiome to study community level survival and examining the network of the community as a whole and the relationship between its topology and biological processes. TDA produces a geometric representation of complex data sets. Network analysis showed that Shiga toxin negative strain E. coli O157:H7 4554 survived significantly longer in comparison to E. coli O157:H7 EDL 933, while the survival time of E. coli O157:NM was comparable to that of E. coli O157:H7 EDL 933 in all of the tested soils. Two non-O157 strains, E. coli O26:H11 and E. coli O103:H2 survived much longer than E. coli O91:H21 and the three strains of E. coli O157. We show that there are complex interactions between E. coli strain survival, microbial community structures, and soil parameters.

Uncovering the human methyltransferasome.

We present a comprehensive analysis of the human methyltransferasome. Primary sequences, predicted secondary structures, and solved crystal structures of known methyltransferases were analyzed by hidden Markov models, Fisher-based statistical matrices, and fold recognition prediction-based threading algorithms to create a model, or profile, of each methyltransferase superfamily. These profiles were used to scan the human proteome database and detect novel methyltransferases. 208 proteins in the human genome are now identified as known or putative methyltransferases, including 38 proteins that were not annotated previously. To date, 30% of these proteins have been linked to disease states. Possible substrates of methylation for all of the SET domain and SPOUT methyltransferases as well as 100 of the 131 seven-ß-strand methyltransferases were surmised from sequence similarity clusters based on alignments of the substrate-specific domains.

Multiple Motif Scanning to identify methyltransferases from the yeast proteome.

A new program (Multiple Motif Scanning) was developed to scan the Saccharomyces cerevisiae proteome for Class I S-adenosylmethionine-dependent methyltransferases. Conserved Motifs I, Post I, II, and III were identified and expanded in known methyltransferases by primary sequence and secondary structural analysis through hidden Markov model profiling of both a yeast reference database and a reference database of methyltransferases with solved three-dimensional structures. The roles of the conserved amino acids in the four motifs of the methyltransferase structure and function were then analyzed to expand the previously defined motifs. Fisher-based negative log statistical matrix sets were developed from the prevalence of amino acids in the motifs. Multiple Motif Scanning is able to scan the proteome and score different combinations of the top fitting sequences for each motif. In addition, the program takes into account the conserved number of amino acids between the motifs. The output of the program is a ranked list of proteins that can be used to identify new methyltransferases and to reevaluate the assignment of previously identified putative methyltransferases. The Multiple Motif Scanning program can be used to develop a putative list of enzymes for any type of protein that has one or more motifs conserved at variable spacings and is freely available (www.chem.ucla.edu/files/MotifSetup.Zip). Finally hidden Markov model profile clustering analysis was used to subgroup Class I methyltransferases into groups that reflect their methyl-accepting substrate specificity.

Bioinformatic Identification of Novel Methyltransferases.

Methylation of DNA, protein, and even RNA species are integral processes in epigenesis. Enzymes that catalyze these reactions using the donor S-adenosylmethionine fall into several structurally distinct classes. The members in each class share sequence similarity that can be used to identify additional methyltransferases. Here, we characterize these classes and in silico approaches to infer protein function. Computational methods such as hidden Markov model profiling and the Multiple Motif Scanning program can be used to analyze known methyltransferases and relay information into the prediction of new ones. In some cases, the substrate of methylation can be inferred from hidden Markov model sequence similarity networks. Functional identification of these candidate species is much more difficult; we discuss one biochemical approach.