1
|
Kharaghani A, Tio ES, Milic M, Bennett DA, De Jager PL, Schneider JA, Sun L, Felsky D. Association of whole-person eigen-polygenic risk scores with Alzheimer's disease. Hum Mol Genet 2024; 33:1315-1327. [PMID: 38679805 PMCID: PMC11262744 DOI: 10.1093/hmg/ddae067] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Revised: 03/06/2024] [Accepted: 04/05/2024] [Indexed: 05/01/2024] Open
Abstract
Late-Onset Alzheimer's Disease (LOAD) is a heterogeneous neurodegenerative disorder with complex etiology and high heritability. Its multifactorial risk profile and large portions of unexplained heritability suggest the involvement of yet unidentified genetic risk factors. Here we describe the "whole person" genetic risk landscape of polygenic risk scores for 2218 traits in 2044 elderly individuals and test if novel eigen-PRSs derived from clustered subnetworks of single-trait PRSs can improve the prediction of LOAD diagnosis, rates of cognitive decline, and canonical LOAD neuropathology. Network analyses revealed distinct clusters of PRSs with clinical and biological interpretability. Novel eigen-PRSs (ePRS) from these clusters significantly improved LOAD-related phenotypes prediction over current state-of-the-art LOAD PRS models. Notably, an ePRS representing clusters of traits related to cholesterol levels was able to improve variance explained in a model of the brain-wide beta-amyloid burden by 1.7% (likelihood ratio test P = 9.02 × 10-7). All associations of ePRS with LOAD phenotypes were eliminated by the removal of APOE-proximal loci. However, our association analysis identified modules characterized by PRSs of high cholesterol and LOAD. We believe this is due to the influence of the APOE region from both PRSs. We found significantly higher mean SNP effects for LOAD in the intersecting APOE region SNPs. Combining genetic risk factors for vascular traits and dementia could improve current single-trait PRS models of LOAD, enhancing the use of PRS in risk stratification. Our results are catalogued for the scientific community, to aid in generating new hypotheses based on our maps of clustered PRSs and associations with LOAD-related phenotypes.
Collapse
Affiliation(s)
- Amin Kharaghani
- Krembil Centre for Neuroinformatics, Centre for Addiction and Mental Health, 250 College Street, Toronto, ON M5T 1R8, Canada
- Dalla Lana School of Public Health, University of Toronto, 155 College Street, Toronto, ON M5T 3M7, Canada
| | - Earvin S Tio
- Krembil Centre for Neuroinformatics, Centre for Addiction and Mental Health, 250 College Street, Toronto, ON M5T 1R8, Canada
- Institute of Medical Science, Department of Psychiatry, University of Toronto, 1 King's College Circle, Toronto, ON M5S 1A8, Canada
| | - Milos Milic
- Krembil Centre for Neuroinformatics, Centre for Addiction and Mental Health, 250 College Street, Toronto, ON M5T 1R8, Canada
| | - David A Bennett
- Rush Alzheimer’s Disease Center, Rush University Medical Center, 1750 West Harrison Street, Chicago, IL 60612, United States
| | - Philip L De Jager
- Centre for Translational and Computational Neuroimmunology, Columbia University Medical Center, 622 West 168th Street, New York, NY 10032, United States
| | - Julie A Schneider
- Rush Alzheimer’s Disease Center, Rush University Medical Center, 1750 West Harrison Street, Chicago, IL 60612, United States
| | - Lei Sun
- Dalla Lana School of Public Health, University of Toronto, 155 College Street, Toronto, ON M5T 3M7, Canada
- Department of Statistical Sciences, University of Toronto, 700 University Avenue, Toronto, ON M5G 1X6, Canada
| | - Daniel Felsky
- Krembil Centre for Neuroinformatics, Centre for Addiction and Mental Health, 250 College Street, Toronto, ON M5T 1R8, Canada
- Dalla Lana School of Public Health, University of Toronto, 155 College Street, Toronto, ON M5T 3M7, Canada
- Institute of Medical Science, Department of Psychiatry, University of Toronto, 1 King's College Circle, Toronto, ON M5S 1A8, Canada
- Department of Psychiatry, University of Toronto, 250 College Street, Toronto, ON M5T 1R8, Canada
| |
Collapse
|
2
|
Bleker C, Grady SK, Langston MA. A Comparative Study of Gene Co-Expression Thresholding Algorithms. J Comput Biol 2024; 31:539-548. [PMID: 38781420 DOI: 10.1089/cmb.2024.0509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024] Open
Abstract
The thresholding problem is studied in the context of graph theoretical analysis of gene co-expression data. A number of thresholding methodologies are described, implemented, and tested over a large collection of graphs derived from real high-throughput biological data. Comparative results are presented and discussed.
Collapse
Affiliation(s)
- Carissa Bleker
- Department of Biotechnology and Systems Biology, National Institute of Biology, Ljubljana, Slovenia
| | - Stephen K Grady
- Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, Tennessee, USA
| | - Michael A Langston
- Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, Tennessee, USA
| |
Collapse
|
3
|
Grady SK, Dojcsak L, Harville EW, Wallace ME, Vilda D, Donneyong MM, Hood DB, Valdez RB, Ramesh A, Im W, Matthews-Juarez P, Juarez PD, Langston MA. Seminar: Scalable Preprocessing Tools for Exposomic Data Analysis. ENVIRONMENTAL HEALTH PERSPECTIVES 2023; 131:124201. [PMID: 38109119 PMCID: PMC10727037 DOI: 10.1289/ehp12901] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 11/22/2023] [Accepted: 11/28/2023] [Indexed: 12/19/2023]
Abstract
BACKGROUND The exposome serves as a popular framework in which to study exposures from chemical and nonchemical stressors across the life course and the differing roles that these exposures can play in human health. As a result, data relevant to the exposome have been used as a resource in the quest to untangle complicated health trajectories and help connect the dots from exposures to adverse outcome pathways. OBJECTIVES The primary aim of this methods seminar is to clarify and review preprocessing techniques critical for accurate and effective external exposomic data analysis. Scalability is emphasized through an application of highly innovative combinatorial techniques coupled with more traditional statistical strategies. The Public Health Exposome is used as an archetypical model. The novelty and innovation of this seminar's focus stem from its methodical, comprehensive treatment of preprocessing and its demonstration of the positive effects preprocessing can have on downstream analytics. DISCUSSION State-of-the-art technologies are described for data harmonization and to mitigate noise, which can stymie downstream interpretation, and to select key exposomic features, without which analytics may lose focus. A main task is the reduction of multicollinearity, a particularly formidable problem that frequently arises from repeated measurements of similar events taken at various times and from multiple sources. Empirical results highlight the effectiveness of a carefully planned preprocessing workflow as demonstrated in the context of more highly concentrated variable lists, improved correlational distributions, and enhanced downstream analytics for latent relationship discovery. The nascent field of exposome science can be characterized by the need to analyze and interpret a complex confluence of highly inhomogeneous spatial and temporal data, which may present formidable challenges to even the most powerful analytical tools. A systematic approach to preprocessing can therefore provide an essential first step in the application of modern computer and data science methods. https://doi.org/10.1289/EHP12901.
Collapse
Affiliation(s)
- Stephen K. Grady
- Graduate School of Genome Science and Technology, University of Tennessee, Knoxville, Tennessee, USA
| | - Levente Dojcsak
- Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, Tennessee, USA
| | - Emily W. Harville
- Department Epidemiology, Tulane University School of Public Health and Tropical Medicine, New Orleans, Louisiana, USA
| | - Maeve E. Wallace
- Department of Social, Behavioral, and Population Sciences, Tulane University School of Public Health and Tropical Medicine, New Orleans, Louisiana, USA
| | - Dovile Vilda
- Department of Social, Behavioral, and Population Sciences, Tulane University School of Public Health and Tropical Medicine, New Orleans, Louisiana, USA
| | | | - Darryl B. Hood
- Division of Environmental Health Sciences, College of Public Health, Ohio State University, Columbus, Ohio, USA
| | - R. Burciaga Valdez
- Department of Economics, University of New Mexico, Albuquerque, New Mexico, USA
| | - Aramandla Ramesh
- Department of Biochemistry, Cancer Biology, Neuroscience & Pharmacology, Meharry Medical College, Nashville, Tennessee, USA
| | - Wansoo Im
- Department of Family and Community Medicine, Meharry Medical College, Nashville, Tennessee, USA
| | | | - Paul D. Juarez
- Department of Family and Community Medicine, Meharry Medical College, Nashville, Tennessee, USA
- Institute on Health Disparities, Equity, and the Exposome, Meharry Medical College, Nashville, Tennessee, USA
| | - Michael A. Langston
- Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, Tennessee, USA
| |
Collapse
|
4
|
Nakamichi K, Akileswaran L, Meirick T, Lee MD, Chodosh J, Rajaiya J, Stroman D, Wolf-Yadlin A, Jackson Q, Holtz WB, Lee AY, Lee CS, Van Gelder RN. Machine Learning Prediction of Adenovirus D8 Conjunctivitis Complications from Viral Whole-Genome Sequence. OPHTHALMOLOGY SCIENCE 2022; 2:100166. [PMID: 36531578 PMCID: PMC9754964 DOI: 10.1016/j.xops.2022.100166] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/23/2021] [Revised: 04/24/2022] [Accepted: 04/25/2022] [Indexed: 11/16/2022]
Abstract
Objective To obtain complete DNA sequences of adenoviral (AdV) D8 genome from patients with conjunctivitis and determine the relation of sequence variation to clinical outcomes. Design This study is a post hoc analysis of banked conjunctival swab samples from the BAYnovation Study, a previously conducted, randomized controlled clinical trial for AdV conjunctivitis. Participants Ninety-six patients with AdV D8-positive conjunctivitis who received placebo treatment in the BAYnovation Study were included in the study. Methods DNA from conjunctival swabs was purified and subjected to whole-genome viral DNA sequencing. Adenovirus D8 variants were identified and correlated with clinical outcomes, including 2 machine learning methods. Main Outcome Measures Viral DNA sequence and development of subepithelial infiltrates (SEIs) were the main outcome measures. Results From initial sequencing of 80 AdV D8-positive samples, full adenoviral genome reconstructions were obtained for 71. A total of 630 single-nucleotide variants were identified, including 156 missense mutations. Sequence clustering revealed 3 previously unappreciated viral clades within the AdV D8 type. The likelihood of SEI development differed significantly between clades, ranging from 83% for Clade 1 to 46% for Clade 3. Genome-wide analysis of viral single-nucleotide polymorphisms failed to identify single-gene determinants of outcome. Two machine learning models were independently trained to predict clinical outcome using polymorphic sequences. Both machine learning models correctly predicted development of SEI outcomes in a newly sequenced validation set of 16 cases (P = 1.5 × 10-5). Prediction was dependent on ensemble groups of polymorphisms across multiple genes. Conclusions Adenovirus D8 has ≥ 3 prevalent molecular substrains, which differ in propensity to result in SEIs. Development of SEIs can be accurately predicted from knowledge of full viral sequence. These results suggest that development of SEIs in AdV D8 conjunctivitis is largely attributable to pathologic viral sequence variants within the D8 type and establishes machine learning paradigms as a powerful technique for understanding viral pathogenicity.
Collapse
Affiliation(s)
- Kenji Nakamichi
- Department of Ophthalmology, University of Washington School of Medicine, Seattle, Washington
- Roger and Angie Karalis Johnson Retina Center, University of Washington School of Medicine, Seattle, Washington
| | - Lakshmi Akileswaran
- Department of Ophthalmology, University of Washington School of Medicine, Seattle, Washington
- Roger and Angie Karalis Johnson Retina Center, University of Washington School of Medicine, Seattle, Washington
| | - Thomas Meirick
- Department of Ophthalmology, University of Washington School of Medicine, Seattle, Washington
| | - Michele D. Lee
- Department of Ophthalmology, University of Washington School of Medicine, Seattle, Washington
| | - James Chodosh
- Department of Ophthalmology, Massachusetts Eye and Ear Infirmary, Harvard Medical School, Boston, Massachusetts
| | - Jaya Rajaiya
- Department of Ophthalmology, Massachusetts Eye and Ear Infirmary, Harvard Medical School, Boston, Massachusetts
| | | | | | | | | | - Aaron Y. Lee
- Department of Ophthalmology, University of Washington School of Medicine, Seattle, Washington
- Roger and Angie Karalis Johnson Retina Center, University of Washington School of Medicine, Seattle, Washington
| | - Cecilia S. Lee
- Department of Ophthalmology, University of Washington School of Medicine, Seattle, Washington
- Roger and Angie Karalis Johnson Retina Center, University of Washington School of Medicine, Seattle, Washington
| | - Russell N. Van Gelder
- Department of Ophthalmology, University of Washington School of Medicine, Seattle, Washington
- Roger and Angie Karalis Johnson Retina Center, University of Washington School of Medicine, Seattle, Washington
- Department of Biological Structure, University of Washington School of Medicine, Seattle, Washington
- Department of Laboratory Medicine and Pathology, University of Washington School of Medicine, Seattle, Washington
| |
Collapse
|
5
|
A Clinical Investigation on the Theragnostic Effect of MicroRNA Biomarkers for Survival Outcome in Cervical Cancer: A PRISMA-P Compliant Protocol for Systematic Review and Comprehensive Meta-Analysis. Genes (Basel) 2022; 13:genes13030463. [PMID: 35328017 PMCID: PMC8948737 DOI: 10.3390/genes13030463] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Revised: 02/25/2022] [Accepted: 02/27/2022] [Indexed: 12/09/2022] Open
Abstract
Background: The most often diagnosed malignancy in women worldwide is cancer of the cervix. It is also the most prevalent kind of gynecological cancer in women. This cancer originates in the opening of the cervix and spreads through sexual contact. Even though human papillomavirus (HPV) may not cause cancer immediately, it does develop over time as a result of the virus’s lengthy persistence to cause dysplastic changes overtime, particularly in high-risk kinds. The primary objective of this research is to see if miRNAs are dysregulated as a result of treatment resistance in cervical cancer (CC). The aim is to see if these microRNAs may be utilized as biomarkers for detecting chemoresistance in CC, particularly for clinical applications. Methods: The recommended protocol for comprehensive study and meta-analysis (PRISMA-P) standards will be utilized for the analysis and data interpretation. The bibliographic databases will be methodically searched using a combination of search keywords. Based on established inclusion and exclusion criteria, the acquired findings will be reviewed, and data retrieved from the selected scientific papers for systematic review. We will then construct a forest from the pooled Hazard ratio (HR) and 95% C.I. values, data obtained using the random-effects model. Discussion: The focus of this study is to identify the function of miRNAs as a chemoresistance regulator and determine if they have the potential scope to be considered as biomarkers for cervical cancer. Through this systematic review and meta-analysis, the goal is to collect, compare, and analyze the data pertaining to the role of miRNAs in cervical cancer, thereby, enabling us to understand the role they play in chemosensitivity.
Collapse
|
6
|
Ray A. Machine learning in postgenomic biology and personalized medicine. WILEY INTERDISCIPLINARY REVIEWS. DATA MINING AND KNOWLEDGE DISCOVERY 2022; 12:e1451. [PMID: 35966173 PMCID: PMC9371441 DOI: 10.1002/widm.1451] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/23/2020] [Accepted: 12/22/2021] [Indexed: 06/15/2023]
Abstract
In recent years Artificial Intelligence in the form of machine learning has been revolutionizing biology, biomedical sciences, and gene-based agricultural technology capabilities. Massive data generated in biological sciences by rapid and deep gene sequencing and protein or other molecular structure determination, on the one hand, requires data analysis capabilities using machine learning that are distinctly different from classical statistical methods; on the other, these large datasets are enabling the adoption of novel data-intensive machine learning algorithms for the solution of biological problems that until recently had relied on mechanistic model-based approaches that are computationally expensive. This review provides a bird's eye view of the applications of machine learning in post-genomic biology. Attempt is also made to indicate as far as possible the areas of research that are poised to make further impacts in these areas, including the importance of explainable artificial intelligence (XAI) in human health. Further contributions of machine learning are expected to transform medicine, public health, agricultural technology, as well as to provide invaluable gene-based guidance for the management of complex environments in this age of global warming.
Collapse
Affiliation(s)
- Animesh Ray
- Riggs School of Applied Life Sciences, Keck Graduate Institute, 535 Watson Drive, Claremont, CA91711, USA
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, California, USA
| |
Collapse
|
7
|
Molecular Investigation of miRNA Biomarkers as Chemoresistance Regulators in Melanoma: A Protocol for Systematic Review and Meta-Analysis. Genes (Basel) 2022; 13:genes13010115. [PMID: 35052456 PMCID: PMC8775297 DOI: 10.3390/genes13010115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2021] [Revised: 12/18/2021] [Accepted: 01/05/2022] [Indexed: 01/25/2023] Open
Abstract
Introduction: Melanoma is a global disease that is predominant in Western countries. However, reliable data resources and comprehensive studies on the theragnostic efficiency of miRNAs in melanoma are scarce. Hence, a decisive study or comprehensive review is required to collate the evidence for profiling miRNAs as a theragnostic marker. This protocol details a comprehensive systematic review and meta-analysis on the impact of miRNAs on chemoresistance and their association with theragnosis in melanoma. Methods and analysis: The articles will be retrieved from online bibliographic databases, including Cochrane Review, EMBASE, MEDLINE, PubMed, Scopus, Science Direct, and Web of Science, with different permutations of ‘keywords’. To obtain full-text papers of relevant research, a stated search method will be used, along with selection criteria. The Preferred Reporting Items for Systematic Reviews and Meta-Analysis for Protocols 2015 (PRISMA-P) standards were used to create this study protocol. The hazard ratio (HR) with a 95% confidence interval will be analyzed using Comprehensive Meta-Analysis (CMA) software 3.0. (CI). The pooled effect size will be calculated using a random or fixed-effects meta-analysis model. Cochran’s Q test and the I2 statistic will be used to determine heterogeneity. Egger’s bias indicator test, Orwin’s and the classic fail-safe N tests, the Begg and Mazumdar rank collection test, and Duval and Tweedie’s trim and fill calculation will all be used to determine publication bias. The overall standard deviation will be evaluated using Z-statistics. Subgroup analyses will be performed according to the melanoma participants’ clinicopathological and biological characteristics and methodological factors if sufficient studies and retrieved data are identified and available. The source of heterogeneity will be assessed using a meta-regression analysis. A pairwise matrix could be developed using either a pairwise correlation or expression associations of miRNA with patients’ survival for the same studies.
Collapse
|
8
|
Barr JR, Shaw P, Abu-Khzam FN, Thatcher T, Yu S. Vulnerability Rating of Source Code with Token Embedding and Combinatorial Algorithms. INTERNATIONAL JOURNAL OF SEMANTIC COMPUTING 2021. [DOI: 10.1142/s1793351x20500087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
We present an empirical analysis of the source code of the Fluoride Bluetooth module, which is a part of standard Android OS distribution, by exhibiting a novel approach for classifying and scoring source code and vulnerability rating. Our workflow combines deep learning, combinatorial optimization, heuristics and machine learning. A combination of heuristics and deep learning is used to embed function (and method) labels into a low-dimensional Euclidean space. Because the corpus of the Fluoride source code is rather limited (containing approximately 12,000 functions), a straightforward embedding (using, e.g. code2vec) is untenable. To overcome the challenge of dearth of data, it is necessary to go through an intermediate step of Byte-Pair Encoding. Subsequently, we embed the tokens from which we assemble an embedding of function/method labels. Long short-term memory network (LSTM) is used to embed tokens. The next step is to form a distance matrix consisting of the cosines between every pairs of vectors (function embedding) which in turn is interpreted as a (combinatorial) graph whose vertices represent functions, and edges correspond to entries whose value exceed some given threshold. Cluster-Editing is then applied to partition the vertex set of the graph into subsets representing “dense graphs,” that are nearly complete subgraphs. Finally, the vectors representing the components, plus additional heuristic-based features are used as features to model the components for vulnerability risk.
Collapse
Affiliation(s)
| | - Peter Shaw
- Nanjing University of Information Science & Technology, Jiangsu, P. R. China
- School of Child Health, Menzies Darwin, Australia
| | | | | | - Sheng Yu
- University of California, Riverside, USA
| |
Collapse
|
9
|
Tabatabai M, Bailey S, Bursac Z, Tabatabai H, Wilus D, Singh KP. An introduction to new robust linear and monotonic correlation coefficients. BMC Bioinformatics 2021; 22:170. [PMID: 33789571 PMCID: PMC8011137 DOI: 10.1186/s12859-021-04098-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2020] [Accepted: 03/22/2021] [Indexed: 01/10/2023] Open
Abstract
BACKGROUND The most common measure of association between two continuous variables is the Pearson correlation (Maronna et al. in Safari an OMC. Robust statistics, 2019. https://login.proxy.bib.uottawa.ca/login?url=https://learning.oreilly.com/library/view/-/9781119214687/?ar&orpq&email=^u). When outliers are present, Pearson does not accurately measure association and robust measures are needed. This article introduces three new robust measures of correlation: Taba (T), TabWil (TW), and TabWil rank (TWR). The correlation estimators T and TW measure a linear association between two continuous or ordinal variables; whereas TWR measures a monotonic association. The robustness of these proposed measures in comparison with Pearson (P), Spearman (S), Quadrant (Q), Median (M), and Minimum Covariance Determinant (MCD) are examined through simulation. Taba distance is used to analyze genes, and statistical tests were used to identify those genes most significantly associated with Williams Syndrome (WS). RESULTS Based on the root mean square error (RMSE) and bias, the three proposed correlation measures are highly competitive when compared to classical measures such as P and S as well as robust measures such as Q, M, and MCD. Our findings indicate TBL2 was the most significant gene among patients diagnosed with WS and had the most significant reduction in gene expression level when compared with control (P value = 6.37E-05). CONCLUSIONS Overall, when the distribution is bivariate Log-Normal or bivariate Weibull, TWR performs best in terms of bias and T performs best with respect to RMSE. Under the Normal distribution, MCD performs well with respect to bias and RMSE; but TW, TWR, T, S, and P correlations were in close proximity. The identification of TBL2 may serve as a diagnostic tool for WS patients. A Taba R package has been developed and is available for use to perform all necessary computations for the proposed methods.
Collapse
Affiliation(s)
| | | | - Zoran Bursac
- Department of Biostatistics, Florida International University, Miami, FL 33199 USA
| | - Habib Tabatabai
- Department of Civil and Environmental Engineering, University of Wisconsin Milwaukee, Milwaukee, WI 53211 USA
| | - Derek Wilus
- Meharry Medical College, Nashville, TN 37208 USA
| | - Karan P. Singh
- Department of Epidemiology and Biostatistics, University of Texas Health Sciences Center at Tyler, Tyler, TX 75708 USA
| |
Collapse
|
10
|
Volpiano CG, Sant’Anna FH, Ambrosini A, de São José JFB, Beneduzi A, Whitman WB, de Souza EM, Lisboa BB, Vargas LK, Passaglia LMP. Genomic Metrics Applied to Rhizobiales ( Hyphomicrobiales): Species Reclassification, Identification of Unauthentic Genomes and False Type Strains. Front Microbiol 2021; 12:614957. [PMID: 33841347 PMCID: PMC8026895 DOI: 10.3389/fmicb.2021.614957] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Accepted: 03/04/2021] [Indexed: 12/13/2022] Open
Abstract
Taxonomic decisions within the order Rhizobiales have relied heavily on the interpretations of highly conserved 16S rRNA sequences and DNA-DNA hybridizations (DDH). Currently, bacterial species are defined as including strains that present 95-96% of average nucleotide identity (ANI) and 70% of digital DDH (dDDH). Thus, ANI values from 520 genome sequences of type strains from species of Rhizobiales order were computed. From the resulting 270,400 comparisons, a ≥95% cut-off was used to extract high identity genome clusters through enumerating maximal cliques. Coupling this graph-based approach with dDDH from clusters of interest, it was found that: (i) there are synonymy between Aminobacter lissarensis and Aminobacter carboxidus, Aurantimonas manganoxydans and Aurantimonas coralicida, "Bartonella mastomydis," and Bartonella elizabethae, Chelativorans oligotrophicus, and Chelativorans multitrophicus, Rhizobium azibense, and Rhizobium gallicum, Rhizobium fabae, and Rhizobium pisi, and Rhodoplanes piscinae and Rhodoplanes serenus; (ii) Chelatobacter heintzii is not a synonym of Aminobacter aminovorans; (iii) "Bartonella vinsonii" subsp. arupensis and "B. vinsonii" subsp. berkhoffii represent members of different species; (iv) the genome accessions GCF_003024615.1 ("Mesorhizobium loti LMG 6,125T"), GCF_003024595.1 ("Mesorhizobium plurifarium LMG 11,892T"), GCF_003096615.1 ("Methylobacterium organophilum DSM 760T"), and GCF_000373025.1 ("R. gallicum R-602 spT") are not from the genuine type strains used for the respective species descriptions; and v) "Xanthobacter autotrophicus" Py2 and "Aminobacter aminovorans" KCTC 2,477T represent cases of misuse of the term "type strain". Aminobacter heintzii comb. nov. and the reclassification of Aminobacter ciceronei as A. heintzii is also proposed. To facilitate the downstream analysis of large ANI matrices, we introduce here ProKlust ("Prokaryotic Clusters"), an R package that uses a graph-based approach to obtain, filter, and visualize clusters on identity/similarity matrices, with settable cut-off points and the possibility of multiple matrices entries.
Collapse
Affiliation(s)
- Camila Gazolla Volpiano
- Departamento de Genética, Instituto de Biociências, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | - Fernando Hayashi Sant’Anna
- Departamento de Genética, Instituto de Biociências, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | - Adriana Ambrosini
- Departamento de Genética, Instituto de Biociências, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | | | - Anelise Beneduzi
- Departamento de Diagnóstico e Pesquisa Agropecuária, Secretaria Estadual da Agricultura, Pecuária e Desenvolvimento Rural, Porto Alegre, Brazil
| | - William B. Whitman
- Department of Microbiology, University of Georgia, Athens, GA, United States
| | - Emanuel Maltempi de Souza
- Departamento de Bioquímica e Biologia Molecular, Centro Politécnico, Setor de Ciências Biológicas, Universidade Federal do Paraná, Curitiba, Brazil
| | - Bruno Brito Lisboa
- Departamento de Diagnóstico e Pesquisa Agropecuária, Secretaria Estadual da Agricultura, Pecuária e Desenvolvimento Rural, Porto Alegre, Brazil
| | - Luciano Kayser Vargas
- Departamento de Diagnóstico e Pesquisa Agropecuária, Secretaria Estadual da Agricultura, Pecuária e Desenvolvimento Rural, Porto Alegre, Brazil
| | | |
Collapse
|
11
|
Hagan RD, Langston MA. Molecular Subtyping and Outlier Detection in Human Disease Using the Paraclique Algorithm. ALGORITHMS 2021; 14:63. [PMID: 36092474 PMCID: PMC9455766 DOI: 10.3390/a14020063] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Recent discoveries of distinct molecular subtypes have led to remarkable advances in treatment for a variety of diseases. While subtyping via unsupervised clustering has received a great deal of interest, most methods rely on basic statistical or machine learning methods. At the same time, techniques based on graph clustering, particularly clique-based strategies, have been successfully used to identify disease biomarkers and gene networks. A graph theoretical approach based on the paraclique algorithm is described that can easily be employed to identify putative disease subtypes and serve as an aid in outlier detection as well. The feasibility and potential effectiveness of this method is demonstrated on publicly available gene co-expression data derived from patient samples covering twelve different disease families.
Collapse
|
12
|
Bannantine JP, Conde C, Bayles DO, Branger M, Biet F. Genetic Diversity Among Mycobacterium avium Subspecies Revealed by Analysis of Complete Genome Sequences. Front Microbiol 2020; 11:1701. [PMID: 32849358 PMCID: PMC7426613 DOI: 10.3389/fmicb.2020.01701] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2020] [Accepted: 06/29/2020] [Indexed: 11/13/2022] Open
Abstract
Mycobacterium avium comprises four subspecies that contain both human and veterinary pathogens. At the inception of this study, twenty-eight M. avium genomes had been annotated as RefSeq genomes, facilitating direct comparisons. These genomes represent strains from around the world and provided a unique opportunity to examine genome dynamics in this species. Each genome was confirmed to be classified correctly based on SNP genotyping, nucleotide identity and presence/absence of repetitive elements or other typing methods. The Mycobacterium avium subspecies paratuberculosis (Map) genome size and organization was remarkably consistent, averaging 4.8 Mb with a variance of only 29.6 kb among the 13 strains. Comparing recombination events along with the larger genome size and variance observed among Mycobacterium avium subspecies avium (Maa) and Mycobacterium avium subspecies hominissuis (Mah) strains (collectively termed non-Map) suggests horizontal gene transfer occurs in non-Map, but not in Map strains. Overall, M. avium subspecies could be divided into two major sub-divisions, with the Map type II (bovine strains) clustering tightly on one end of a phylogenetic spectrum and Mah strains clustering more loosely together on the other end. The most evolutionarily distinct Map strain was an ovine strain, designated Telford, which had >1,000 SNPs and showed large rearrangements compared to the bovine type II strains. The Telford strain clustered with Maa strains as an intermediate between Map type II and Mah. SNP analysis and genome organization analyses repeatedly demonstrated the conserved nature of Map versus the mosaic nature of non-Map M. avium strains. Finally, core and pangenomes were developed for Map and non-Map strains. A total of 80% Map genes belonged to the Map core genome, while only 40% of non-Map genes belonged to the non-Map core genome. These genomes provide a more complete and detailed comparison of these subspecies strains as well as a blueprint for how genetic diversity originated.
Collapse
Affiliation(s)
- John P Bannantine
- USDA-Agricultural Research Service, National Animal Disease Center, Ames, IA, United States
| | - Cyril Conde
- INRAE, Université de Tours, ISP, Nouzilly, France
| | - Darrell O Bayles
- USDA-Agricultural Research Service, National Animal Disease Center, Ames, IA, United States
| | | | - Franck Biet
- INRAE, Université de Tours, ISP, Nouzilly, France
| |
Collapse
|
13
|
Randhawa V, Pathania S. Advancing from protein interactomes and gene co-expression networks towards multi-omics-based composite networks: approaches for predicting and extracting biological knowledge. Brief Funct Genomics 2020; 19:364-376. [PMID: 32678894 DOI: 10.1093/bfgp/elaa015] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Revised: 05/31/2020] [Accepted: 06/15/2020] [Indexed: 01/17/2023] Open
Abstract
Prediction of biological interaction networks from single-omics data has been extensively implemented to understand various aspects of biological systems. However, more recently, there is a growing interest in integrating multi-omics datasets for the prediction of interactomes that provide a global view of biological systems with higher descriptive capability, as compared to single omics. In this review, we have discussed various computational approaches implemented to infer and analyze two of the most important and well studied interactomes: protein-protein interaction networks and gene co-expression networks. We have explicitly focused on recent methods and pipelines implemented to infer and extract biologically important information from these interactomes, starting from utilizing single-omics data and then progressing towards multi-omics data. Accordingly, recent examples and case studies are also briefly discussed. Overall, this review will provide a proper understanding of the latest developments in protein and gene network modelling and will also help in extracting practical knowledge from them.
Collapse
Affiliation(s)
- Vinay Randhawa
- Department of Biochemistry, Panjab University, Chandigarh, 160014, India
| | - Shivalika Pathania
- Department of Biotechnology, Panjab University, Chandigarh, 160014, India
| |
Collapse
|
14
|
Lu Y, Phillips CA, Langston MA. A robustness metric for biological data clustering algorithms. BMC Bioinformatics 2019; 20:503. [PMID: 31874625 PMCID: PMC6929270 DOI: 10.1186/s12859-019-3089-6] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2019] [Accepted: 09/10/2019] [Indexed: 02/05/2023] Open
Abstract
BACKGROUND Cluster analysis is a core task in modern data-centric computation. Algorithmic choice is driven by factors such as data size and heterogeneity, the similarity measures employed, and the type of clusters sought. Familiarity and mere preference often play a significant role as well. Comparisons between clustering algorithms tend to focus on cluster quality. Such comparisons are complicated by the fact that algorithms often have multiple settings that can affect the clusters produced. Such a setting may represent, for example, a preset variable, a parameter of interest, or various sorts of initial assignments. A question of interest then is this: to what degree do the clusters produced vary as setting values change? RESULTS This work introduces a new metric, termed simply "robustness", designed to answer that question. Robustness is an easily-interpretable measure of the propensity of a clustering algorithm to maintain output coherence over a range of settings. The robustness of eleven popular clustering algorithms is evaluated over some two dozen publicly available mRNA expression microarray datasets. Given their straightforwardness and predictability, hierarchical methods generally exhibited the highest robustness on most datasets. Of the more complex strategies, the paraclique algorithm yielded consistently higher robustness than other algorithms tested, approaching and even surpassing hierarchical methods on several datasets. Other techniques exhibited mixed robustness, with no clear distinction between them. CONCLUSIONS Robustness provides a simple and intuitive measure of the stability and predictability of a clustering algorithm. It can be a useful tool to aid both in algorithm selection and in deciding how much effort to devote to parameter tuning.
Collapse
Affiliation(s)
- Yuping Lu
- Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, 37996 TN USA
| | - Charles A. Phillips
- Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, 37996 TN USA
| | - Michael A. Langston
- Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, 37996 TN USA
| |
Collapse
|
15
|
Sun S, Miao Z, Ratcliffe B, Campbell P, Pasch B, El-Kassaby YA, Balasundaram B, Chen C. SNP variable selection by generalized graph domination. PLoS One 2019; 14:e0203242. [PMID: 30677030 PMCID: PMC6345469 DOI: 10.1371/journal.pone.0203242] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2018] [Accepted: 01/08/2019] [Indexed: 11/19/2022] Open
Abstract
BACKGROUND High-throughput sequencing technology has revolutionized both medical and biological research by generating exceedingly large numbers of genetic variants. The resulting datasets share a number of common characteristics that might lead to poor generalization capacity. Concerns include noise accumulated due to the large number of predictors, sparse information regarding the p≫n problem, and overfitting and model mis-identification resulting from spurious collinearity. Additionally, complex correlation patterns are present among variables. As a consequence, reliable variable selection techniques play a pivotal role in predictive analysis, generalization capability, and robustness in clustering, as well as interpretability of the derived models. METHODS AND FINDINGS K-dominating set, a parameterized graph-theoretic generalization model, was used to model SNP (single nucleotide polymorphism) data as a similarity network and searched for representative SNP variables. In particular, each SNP was represented as a vertex in the graph, (dis)similarity measures such as correlation coefficients or pairwise linkage disequilibrium were estimated to describe the relationship between each pair of SNPs; a pair of vertices are adjacent, i.e. joined by an edge, if the pairwise similarity measure exceeds a user-specified threshold. A minimum k-dominating set in the SNP graph was then made as the smallest subset such that every SNP that is excluded from the subset has at least k neighbors in the selected ones. The strength of k-dominating set selection in identifying independent variables, and in culling representative variables that are highly correlated with others, was demonstrated by a simulated dataset. The advantages of k-dominating set variable selection were also illustrated in two applications: pedigree reconstruction using SNP profiles of 1,372 Douglas-fir trees, and species delineation for 226 grasshopper mouse samples. A C++ source code that implements SNP-SELECT and uses Gurobi optimization solver for the k-dominating set variable selection is available (https://github.com/transgenomicsosu/SNP-SELECT).
Collapse
Affiliation(s)
- Shuzhen Sun
- Department of Biochemistry and Molecular Biology, Oklahoma State University, Stillwater, United States of America
- Department of Forest and Conservation Sciences, Faculty of Forestry, The University of British Columbia, Vancouver, B.C. Canada
| | - Zhuqi Miao
- Center for Health Systems Innovation, Oklahoma State University, Stillwater, United States of America
| | - Blaise Ratcliffe
- Department of Forest and Conservation Sciences, Faculty of Forestry, The University of British Columbia, Vancouver, B.C. Canada
| | - Polly Campbell
- Department of Integrative Biology, Oklahoma State University, Stillwater, United States of America
- Department of Evolution, Ecology and Organismal Biology, University of California, Riverside, Riverside, United States of America
| | - Bret Pasch
- Department of Biological Sciences, Northern Arizona University, Flagstaff, United States of America
| | - Yousry A. El-Kassaby
- Department of Forest and Conservation Sciences, Faculty of Forestry, The University of British Columbia, Vancouver, B.C. Canada
| | - Balabhaskar Balasundaram
- School of Industrial Engineering and Management, Oklahoma State University, Stillwater, United States of America
| | - Charles Chen
- Department of Biochemistry and Molecular Biology, Oklahoma State University, Stillwater, United States of America
- * E-mail:
| |
Collapse
|
16
|
Payá-Milans M, Olmstead JW, Nunez G, Rinehart TA, Staton M. Comprehensive evaluation of RNA-seq analysis pipelines in diploid and polyploid species. Gigascience 2018; 7:5168871. [PMID: 30418578 PMCID: PMC6275443 DOI: 10.1093/gigascience/giy132] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2018] [Accepted: 10/21/2018] [Indexed: 11/12/2022] Open
Abstract
Background The usual analysis of RNA sequencing (RNA-seq) reads is based on an existing reference genome and annotated gene models. However, when a reference for the sequenced species is not available, alternatives include using a reference genome from a related species or reconstructing transcript sequences with de novo assembly. In addition, researchers are faced with many options for RNA-seq data processing and limited information on how their decisions will impact the final outcome. Using both a diploid and polyploid species with a distant reference genome, we have tested the influence of different tools at various steps of a typical RNA-seq analysis workflow on the recovery of useful processed data available for downstream analysis. Findings At the preprocessing step, we found error correction has a strong influence on de novo assembly but not on mapping results. After trimming, a greater percentage of reads could be used in downstream analysis by selecting gentle quality trimming performed with Skewer instead of strict quality trimming with Trimmomatic. This availability of reads correlated with size, quality, and completeness of de novo assemblies and with number of mapped reads. When selecting a reference genome from a related species to map reads, outcome was significantly improved when using mapping software tolerant of greater sequence divergence, such as Stampy or GSNAP. Conclusions The selection of bioinformatic software tools for RNA-seq data analysis can maximize quality parameters on de novo assemblies and availability of reads in downstream analysis.
Collapse
Affiliation(s)
- Miriam Payá-Milans
- Department of Entomology and Plant Pathology, University of Tennessee, 370 PBB, 2505 EJ Chapman Blvd, Knoxville, TN, 37996, United States
| | - James W Olmstead
- Horticultural Sciences Department, University of Florida, 2550 Hull Rd, PO Box 110690, Gainesville, FL, 32611, United States
| | - Gerardo Nunez
- Horticultural Sciences Department, University of Florida, 2550 Hull Rd, PO Box 110690, Gainesville, FL, 32611, United States
| | - Timothy A Rinehart
- Thad Cochran Southern Horticultural Laboratory, USDA-Agricultural Research Service, PO Box 287, Poplarville, MS, 39470, United States.,Crop Production and Protection, USDA-Agricultural Research Service, 5601 Sunnyside Ave, Beltsville, MD, 20705, United States
| | - Margaret Staton
- Department of Entomology and Plant Pathology, University of Tennessee, 370 PBB, 2505 EJ Chapman Blvd, Knoxville, TN, 37996, United States
| |
Collapse
|
17
|
Li X, LeBlanc J, Elashoff D, McHardy I, Tong M, Roth B, Ippoliti A, Barron G, McGovern D, McDonald K, Newberry R, Graeber T, Horvath S, Goodglick L, Braun J. Microgeographic Proteomic Networks of the Human Colonic Mucosa and Their Association With Inflammatory Bowel Disease. Cell Mol Gastroenterol Hepatol 2016; 2:567-583. [PMID: 28174738 PMCID: PMC5042708 DOI: 10.1016/j.jcmgh.2016.05.003] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/20/2016] [Accepted: 05/06/2016] [Indexed: 12/28/2022]
Abstract
BACKGROUND & AIMS Interactions between mucosal cell types, environmental stressors, and intestinal microbiota contribute to pathogenesis in inflammatory bowel disease (IBD). Here, we applied metaproteomics of the mucosal-luminal interface to study the disease-related biology of the human colonic mucosa. METHODS We recruited a discovery cohort of 51 IBD and non-IBD subjects endoscopically sampled by mucosal lavage at 6 colonic regions, and a validation cohort of 38 no-IBD subjects. Metaproteome data sets were produced for each sample and analyzed for association with colonic site and disease state using a suite of bioinformatic approaches. Localization of select proteins was determined by immunoblot analysis and immunohistochemistry of human endoscopic biopsy samples. RESULTS Co-occurrence analysis of the discovery cohort metaproteome showed that proteins at the mucosal surface clustered into modules with evidence of differential functional specialization (eg, iron regulation, microbial defense) and cellular origin (eg, epithelial or hemopoietic). These modules, validated in an independent cohort, were differentially associated spatially along the gastrointestinal tract, and 7 modules were associated selectively with non-IBD, ulcerative colitis, and/or Crohn's disease states. In addition, the detailed composition of certain modules was altered in disease vs healthy states. We confirmed the predicted spatial and disease-associated localization of 28 proteins representing 4 different disease-related modules by immunoblot and immunohistochemistry visualization, with evidence for their distribution as millimeter-scale microgeographic mosaic. CONCLUSIONS These findings suggest that the mucosal surface is a microgeographic mosaic of functional networks reflecting the local mucosal ecology, whose compositional differences in disease and healthy samples may provide a unique readout of physiologic and pathologic mucosal states.
Collapse
Key Words
- ANOVA, analysis of variance
- CD, Crohn’s disease
- Ecology
- HBD, human β-defensin
- HD5, human alpha defensin 5
- HNP, human neutrophil peptide
- HPLC, high-performance liquid chromatography
- IBD, inflammatory bowel disease
- IHC, immunohistochemistry
- Inflammatory Bowel Disease
- MALDI, matrix-assisted laser desorption/ionization
- MFN, mucosal functional network
- MLI, mucosal–luminal interface
- MS/MS, tandem mass spectrometry
- Metaproteomics
- Mucosal
- NLME, nonlinear mixed-effect model
- Networks
- PVCA, principal variance component analysis
- TOF, time of flight
- UC, ulcerative colitis
- WGCNA, weighted correlation network analysis
Collapse
Affiliation(s)
- Xiaoxiao Li
- Department of Molecular and Medical Pharmacology, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California,Department of Pathology and Laboratory Medicine, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California,Inflammatory Bowel and Immunobiology Research Institute, Cedars-Sinai Medical Center, Los Angeles, California
| | - James LeBlanc
- Department of Pathology and Laboratory Medicine, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California
| | - David Elashoff
- Department of Medicine, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California
| | - Ian McHardy
- Department of Pathology and Laboratory Medicine, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California
| | - Maomeng Tong
- Department of Molecular and Medical Pharmacology, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California
| | - Bennett Roth
- Department of Medicine, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California
| | - Andrew Ippoliti
- Inflammatory Bowel and Immunobiology Research Institute, Cedars-Sinai Medical Center, Los Angeles, California
| | - Gildardo Barron
- Inflammatory Bowel and Immunobiology Research Institute, Cedars-Sinai Medical Center, Los Angeles, California
| | - Dermot McGovern
- Inflammatory Bowel and Immunobiology Research Institute, Cedars-Sinai Medical Center, Los Angeles, California
| | - Keely McDonald
- Department of Internal Medicine, Washington University School of Medicine, St. Louis, Missouri
| | - Rodney Newberry
- Department of Internal Medicine, Washington University School of Medicine, St. Louis, Missouri
| | - Thomas Graeber
- Department of Molecular and Medical Pharmacology, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California
| | - Steve Horvath
- Department of Human Genetics and Biostatistics, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California
| | - Lee Goodglick
- Department of Pathology and Laboratory Medicine, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California
| | - Jonathan Braun
- Department of Molecular and Medical Pharmacology, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California,Department of Pathology and Laboratory Medicine, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California,Correspondence Address correspondence to: Jonathan Braun, MD, PhD, Department of Pathology and Laboratory Medicine, University of California Los Angeles David Geffen School of Medicine, Los Angeles, California 90095. fax: (310) 267-4486.Department of Pathology and Laboratory Medicine, University of California Los Angeles David Geffen School of MedicineLos AngelesCalifornia 90095
| |
Collapse
|
18
|
Hagan RD, Langston MA, Wang K. Lower Bounds on Paraclique Density. DISCRETE APPLIED MATHEMATICS (AMSTERDAM, NETHERLANDS : 1988) 2016; 204:208-212. [PMID: 27057077 PMCID: PMC4820293 DOI: 10.1016/j.dam.2015.11.010] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
The scientific literature teems with clique-centric clustering strategies. In this paper we analyze one such method, the paraclique algorithm. Paraclique has found practical utility in a variety of application domains, and has been successfully employed to reduce the effects of noise. Nevertheless, its formal analysis and worst-case guarantees have remained elusive. We address this issue by deriving a series of lower bounds on paraclique densities.
Collapse
|
19
|
Jothi R, Mohanty SK, Ojha A. Functional grouping of similar genes using eigenanalysis on minimum spanning tree based neighborhood graph. Comput Biol Med 2016; 71:135-48. [PMID: 26945461 DOI: 10.1016/j.compbiomed.2016.02.007] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2015] [Revised: 01/16/2016] [Accepted: 02/12/2016] [Indexed: 10/22/2022]
Abstract
Gene expression data clustering is an important biological process in DNA microarray analysis. Although there have been many clustering algorithms for gene expression analysis, finding a suitable and effective clustering algorithm is always a challenging problem due to the heterogeneous nature of gene profiles. Minimum Spanning Tree (MST) based clustering algorithms have been successfully employed to detect clusters of varying shapes and sizes. This paper proposes a novel clustering algorithm using Eigenanalysis on Minimum Spanning Tree based neighborhood graph (E-MST). As MST of a set of points reflects the similarity of the points with their neighborhood, the proposed algorithm employs a similarity graph obtained from k(') rounds of MST (k(')-MST neighborhood graph). By studying the spectral properties of the similarity matrix obtained from k(')-MST graph, the proposed algorithm achieves improved clustering results. We demonstrate the efficacy of the proposed algorithm on 12 gene expression datasets. Experimental results show that the proposed algorithm performs better than the standard clustering algorithms.
Collapse
Affiliation(s)
- R Jothi
- Indian Institute of Information Technology, Design and Manufacturing Jabalpur, Madhya Pradesh, India.
| | - Sraban Kumar Mohanty
- Indian Institute of Information Technology, Design and Manufacturing Jabalpur, Madhya Pradesh, India.
| | - Aparajita Ojha
- Indian Institute of Information Technology, Design and Manufacturing Jabalpur, Madhya Pradesh, India.
| |
Collapse
|
20
|
Villalobos-Cid M, Chacón M, Zitko P, Instroza-Ponta M. A New Strategy to Evaluate Technical Efficiency in Hospitals Using Homogeneous Groups of Casemix : How to Evaluate When There is Not DRGs? J Med Syst 2016; 40:103. [PMID: 26880102 DOI: 10.1007/s10916-016-0458-9] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2015] [Accepted: 01/29/2016] [Indexed: 01/16/2023]
Abstract
The public health system has restricted economic resources. Because of that, it is necessary to know how the resources are being used and if they are properly distributed. Several works have applied classical approaches based in Data Envelopment Analysis (DEA) and Stochastic Frontier Analysis (SFA) for this purpose. However, if we have hospitals with different casemix, this is not the best approach. In order to avoid biases in the comparisons, other works have recommended the use of hospital production data corrected by the weights from Diagnosis Related Groups (DRGs), to adjust the casemix of hospitals. However, not all countries have this tool fully implemented, which limits the efficiency evaluation. This paper proposes a new approach for evaluating the efficiency of hospitals. It uses a graph-based clustering algorithm to find groups of hospitals that have similar production profiles. Then, DEA is used to evaluate the technical efficiency of each group. The proposed approach is tested using the production data from 2014 of 193 Chilean public hospitals. The results allowed to identify different performance profiles of each group, that differs from other studies that employs data from partially implemented DRGs. Our results are able to deliver a better description of the resource management of the different groups of hospitals. We have created a website with the results ( bioinformatic.diinf.usach.cl/publichealth ). Data can be requested to the authors.
Collapse
Affiliation(s)
- Manuel Villalobos-Cid
- Departamento de Ingeniería Informática, Facultad de Ingeniería, Universidad de Santiago de Chile, Santiago, Chile
| | - Max Chacón
- Departamento de Ingeniería Informática, Facultad de Ingeniería, Universidad de Santiago de Chile, Santiago, Chile
| | - Pedro Zitko
- Unidad de Estudios Asistenciales, Hospital Barros Luco Trudeau, Facultad de Medicina, Universidad Diego Portales, Santiago, Chile
| | - Mario Instroza-Ponta
- Departamento de Ingeniería Informática, Facultad de Ingeniería, Universidad de Santiago de Chile, Santiago, Chile.
| |
Collapse
|
21
|
Rodenas-Cuadrado P, Chen XS, Wiegrebe L, Firzlaff U, Vernes SC. A novel approach identifies the first transcriptome networks in bats: a new genetic model for vocal communication. BMC Genomics 2015; 16:836. [PMID: 26490347 PMCID: PMC4618519 DOI: 10.1186/s12864-015-2068-1] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2015] [Accepted: 10/13/2015] [Indexed: 12/15/2022] Open
Abstract
BACKGROUND Bats are able to employ an astonishingly complex vocal repertoire for navigating their environment and conveying social information. A handful of species also show evidence for vocal learning, an extremely rare ability shared only with humans and few other animals. However, despite their potential for the study of vocal communication, bats remain severely understudied at a molecular level. To address this fundamental gap we performed the first transcriptome profiling and genetic interrogation of molecular networks in the brain of a highly vocal bat species, Phyllostomus discolor. RESULTS Gene network analysis typically needs large sample sizes for correct clustering, this can be prohibitive where samples are limited, such as in this study. To overcome this, we developed a novel bioinformatics methodology for identifying robust co-expression gene networks using few samples (N=6). Using this approach, we identified tissue-specific functional gene networks from the bat PAG, a brain region fundamental for mammalian vocalisation. The most highly connected network identified represented a cluster of genes involved in glutamatergic synaptic transmission. Glutamatergic receptors play a significant role in vocalisation from the PAG, suggesting that this gene network may be mechanistically important for vocal-motor control in mammals. CONCLUSION We have developed an innovative approach to cluster co-expressing gene networks and show that it is highly effective in detecting robust functional gene networks with limited sample sizes. Moreover, this work represents the first gene network analysis performed in a bat brain and establishes bats as a novel, tractable model system for understanding the genetics of vocal mammalian communication.
Collapse
Affiliation(s)
- Pedro Rodenas-Cuadrado
- Max Planck Institute for Psycholinguistics, Wundtlaan 1, Nijmegen, 6525 XD, The Netherlands.
| | - Xiaowei Sylvia Chen
- Max Planck Institute for Psycholinguistics, Wundtlaan 1, Nijmegen, 6525 XD, The Netherlands.
| | - Lutz Wiegrebe
- Ludwig-Maximilians-Universität, Division of Neurobiology, Department Biology II, Großhaderner Straße 2, Planegg-Martinsried, Munich, D-82152, Germany.
| | - Uwe Firzlaff
- Lehrstuhl für Zoologie, TU München, Liesel-Beckmann-Str. 4, Freising-Weihenstephan, Munich, 85350, Germany.
| | - Sonja C Vernes
- Max Planck Institute for Psycholinguistics, Wundtlaan 1, Nijmegen, 6525 XD, The Netherlands. .,Donders Centre for Cognitive Neuroimaging, Kapittelweg 29, Nijmegen, 6525 EN, The Netherlands.
| |
Collapse
|
22
|
Pluta T, Bernardo R, Shin HW, Bernardo DR. Unsupervised learning of electrocorticography motifs with binary descriptors of wavelet features and hierarchical clustering. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2015; 2014:2657-60. [PMID: 25570537 DOI: 10.1109/embc.2014.6944169] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
We describe a novel method for data mining spectro-spatiotemporal network motifs from electrocorticographic (ECoG) data. The method utilizes wavelet feature extraction from ECoG data, generation of compact binary vectors from these features, and binary vector hierarchical clustering. The potential utility of this method in the discovery of recurring neural patterns is demonstrated in an example showing clustering of ictal and post-ictal gamma activity patterns. The method allows for the efficient and scalable retrieval and clustering of neural motifs occurring in massive amounts of neural data, such as in prolonged EEG/ECoG recordings and in brain computer interfaces.
Collapse
|
23
|
Iancu OD, Colville A, Oberbeck D, Darakjian P, McWeeney SK, Hitzemann R. Cosplicing network analysis of mammalian brain RNA-Seq data utilizing WGCNA and Mantel correlations. Front Genet 2015; 6:174. [PMID: 26029240 PMCID: PMC4429622 DOI: 10.3389/fgene.2015.00174] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2014] [Accepted: 04/21/2015] [Indexed: 01/06/2023] Open
Abstract
Across species and tissues and especially in the mammalian brain, production of gene isoforms is widespread. While gene expression coordination has been previously described as a scale-free coexpression network, the properties of transcriptome-wide isoform production coordination have been less studied. Here we evaluate the system-level properties of cosplicing in mouse, macaque, and human brain gene expression data using a novel network inference procedure. Genes are represented as vectors/lists of exon counts and distance measures sensitive to exon inclusion rates quantifies differences across samples. For all gene pairs, distance matrices are correlated across samples, resulting in cosplicing or cotranscriptional network matrices. We show that networks including cosplicing information are scale-free and distinct from coexpression. In the networks capturing cosplicing we find a set of novel hubs with unique characteristics distinguishing them from coexpression hubs: heavy representation in neurobiological functional pathways, strong overlap with markers of neurons and neuroglia, long coding lengths, and high number of both exons and annotated transcripts. Further, the cosplicing hubs are enriched in genes associated with autism spectrum disorders. Cosplicing hub homologs across eukaryotes show dramatically increasing intronic lengths but stable coding region lengths. Shared transcription factor binding sites increase coexpression but not cosplicing; the reverse is true for splicing-factor binding sites. Genes with protein-protein interactions have strong coexpression and cosplicing. Additional factors affecting the networks include shared microRNA binding sites, spatial colocalization within the striatum, and sharing a chromosomal folding domain. Cosplicing network patterns remain relatively stable across species.
Collapse
Affiliation(s)
- Ovidiu D Iancu
- Department of Behavioral Neuroscience, Oregon Health & Science University Portland, OR, USA
| | - Alexandre Colville
- Department of Behavioral Neuroscience, Oregon Health & Science University Portland, OR, USA
| | - Denesa Oberbeck
- Department of Behavioral Neuroscience, Oregon Health & Science University Portland, OR, USA
| | - Priscila Darakjian
- Department of Behavioral Neuroscience, Oregon Health & Science University Portland, OR, USA
| | - Shannon K McWeeney
- Division of Biostatistics, Public Health and Preventative Medicine, Oregon Health & Science University Portland, OR, USA
| | - Robert Hitzemann
- Department of Behavioral Neuroscience, Oregon Health & Science University Portland, OR, USA ; Research Service, Veterans Affairs Medical Center Portland, OR, USA
| |
Collapse
|
24
|
McKinney GJ, Hale MC, Goetz G, Gribskov M, Thrower FP, Nichols KM. Ontogenetic changes in embryonic and brain gene expression in progeny produced from migratory and residentOncorhynchus mykiss. Mol Ecol 2015; 24:1792-809. [DOI: 10.1111/mec.13143] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2015] [Revised: 02/27/2015] [Accepted: 03/02/2015] [Indexed: 01/17/2023]
Affiliation(s)
- Garrett J. McKinney
- Department of Biological Sciences; Purdue University; West Lafayette IN 47907 USA
- School of Aquatic and Fishery Sciences; University of Washington; Seattle WA 98195-5020 USA
| | - Matthew C. Hale
- Department of Biological Sciences; Purdue University; West Lafayette IN 47907 USA
| | - Giles Goetz
- Conservation Biology Division; Northwest Fisheries Science Center; National Marine Fisheries Service; National Oceanic and Atmospheric Administration; Seattle WA 98112 USA
| | - Michael Gribskov
- Department of Biological Sciences; Purdue University; West Lafayette IN 47907 USA
| | - Frank P. Thrower
- Ted Stevens Marine Research Institute; Alaska Fisheries Science Center; National Marine Fisheries Service; National Oceanic and Atmospheric Administration; Juneau AK 99801 USA
| | - Krista M. Nichols
- Department of Biological Sciences; Purdue University; West Lafayette IN 47907 USA
- Conservation Biology Division; Northwest Fisheries Science Center; National Marine Fisheries Service; National Oceanic and Atmospheric Administration; Seattle WA 98112 USA
| |
Collapse
|
25
|
Barter RL, Schramm SJ, Mann GJ, Yang YH. Network-based biomarkers enhance classical approaches to prognostic gene expression signatures. BMC SYSTEMS BIOLOGY 2014; 8 Suppl 4:S5. [PMID: 25521200 PMCID: PMC4290694 DOI: 10.1186/1752-0509-8-s4-s5] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
BACKGROUND Classical approaches to predicting patient clinical outcome via gene expression information are primarily based on differential expression of unrelated genes (single-gene approaches) or genes related by, for example, biologic pathway or function (gene-sets). Recently, network-based approaches utilising interaction information between genes have emerged. An open problem is whether such approaches add value to the more traditional methods of signature modelling. We explored this question via comparison of the most widely employed single-gene, gene-set, and network-based methods, using gene expression microarray data from two different cancers: melanoma and ovarian. We considered two kinds of network approaches. The first of these identifies informative genes using gene expression and network connectivity information combined, the latter drawn from prior knowledge of protein-protein interactions. The second approach focuses on identification of informative sub-networks (small networks of interacting proteins, again from prior knowledge networks). For all methods we performed 100 rounds of 5-fold cross-validation under 3 different classifiers. For network-based approaches, we considered two different protein-protein interaction networks. We quantified resulting patterns of misclassification and discussed the relative value of each relative to ongoing development of prognostic biomarkers. RESULTS We found that single-gene, gene-set and network methods yielded similar error rates in melanoma and ovarian cancer data. Crucially, however, our novel and detailed patient-level analyses revealed that the different methods were correctly classifying alternate subsets of patients in each cohort. We also found that the network-based NetRank feature selection method was the most stable. CONCLUSIONS Next-generation methods of gene expression signature modelling harness data from external networks and are foreshadowed as a standard mode of analysis. But what do they add to traditional approaches? Our findings indicate there is value in the way in which different subspaces of the patient sample are captured differently among the various methods, highlighting the possibility of 'combination' classifiers capable of identifying which patients will be more accurately classified by one particular method over another. We have seen this clearly for the first time because of our in-depth analysis at the level of individual patients.
Collapse
Affiliation(s)
- Rebecca L Barter
- School of Mathematics and Statistics at The University of Sydney, F07, The University of Sydney, NSW, 2006, Australia
| | - Sarah-Jane Schramm
- Westmead Millennium Institute at The University of Sydney, 176 Hawkesbury Road, Westmead, NSW, 2145, Australia
- Melanoma Institute Australia, 40 Rocklands Rd, North Sydney, NSW, 2060, Australia
| | - Graham J Mann
- Westmead Millennium Institute at The University of Sydney, 176 Hawkesbury Road, Westmead, NSW, 2145, Australia
- Melanoma Institute Australia, 40 Rocklands Rd, North Sydney, NSW, 2060, Australia
| | - Yee Hwa Yang
- School of Mathematics and Statistics at The University of Sydney, F07, The University of Sydney, NSW, 2006, Australia
- Melanoma Institute Australia, 40 Rocklands Rd, North Sydney, NSW, 2060, Australia
| |
Collapse
|
26
|
Exploration of preterm birth rates using the public health exposome database and computational analysis methods. INTERNATIONAL JOURNAL OF ENVIRONMENTAL RESEARCH AND PUBLIC HEALTH 2014; 11:12346-66. [PMID: 25464130 PMCID: PMC4276617 DOI: 10.3390/ijerph111212346] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/08/2014] [Revised: 11/12/2014] [Accepted: 11/19/2014] [Indexed: 01/19/2023]
Abstract
Recent advances in informatics technology has made it possible to integrate, manipulate, and analyze variables from a wide range of scientific disciplines allowing for the examination of complex social problems such as health disparities. This study used 589 county-level variables to identify and compare geographical variation of high and low preterm birth rates. Data were collected from a number of publically available sources, bringing together natality outcomes with attributes of the natural, built, social, and policy environments. Singleton early premature county birth rate, in counties with population size over 100,000 persons provided the dependent variable. Graph theoretical techniques were used to identify a wide range of predictor variables from various domains, including black proportion, obesity and diabetes, sexually transmitted infection rates, mother’s age, income, marriage rates, pollution and temperature among others. Dense subgraphs (paracliques) representing groups of highly correlated variables were resolved into latent factors, which were then used to build a regression model explaining prematurity (R-squared = 76.7%). Two lists of counties with large positive and large negative residuals, indicating unusual prematurity rates given their circumstances, may serve as a starting point for ways to intervene and reduce health disparities for preterm births.
Collapse
|
27
|
Ha T, Swanson D, Larouche M, Glenn R, Weeden D, Zhang P, Hamre K, Langston M, Phillips C, Song M, Ouyang Z, Chesler E, Duvvurru S, Yordanova R, Cui Y, Campbell K, Ricker G, Phillips C, Homayouni R, Goldowitz D. CbGRiTS: cerebellar gene regulation in time and space. Dev Biol 2014; 397:18-30. [PMID: 25446528 DOI: 10.1016/j.ydbio.2014.09.032] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2014] [Revised: 08/23/2014] [Accepted: 09/27/2014] [Indexed: 01/09/2023]
Abstract
The mammalian CNS is one of the most complex biological systems to understand at the molecular level. The temporal information from time series transcriptome analysis can serve as a potent source of associative information between developmental processes and regulatory genes. Here, we introduce a new transcriptome database called, Cerebellar Gene Regulation in Time and Space (CbGRiTS). This dataset is populated with transcriptome data across embryonic and postnatal development from two standard mouse strains, C57BL/6J and DBA/2J, several recombinant inbred lines and cerebellar mutant strains. Users can evaluate expression profiles across cerebellar development in a deep time series with graphical interfaces for data exploration and link-out to anatomical expression databases. We present three analytical approaches that take advantage of specific aspects of the time series for transcriptome analysis. We demonstrate the use of CbGRiTS dataset as a community resource to explore patterns of gene expression and develop hypotheses concerning gene regulatory networks in brain development.
Collapse
Affiliation(s)
- Thomas Ha
- Centre for Molecular Medicine and Therapeutics, Child and Family Research Institute, Department of Medical Genetics, University of British Columbia, 950 West 28th Avenue, Vancouver, BC, Canada V5Z 4H4
| | - Douglas Swanson
- Centre for Molecular Medicine and Therapeutics, Child and Family Research Institute, Department of Medical Genetics, University of British Columbia, 950 West 28th Avenue, Vancouver, BC, Canada V5Z 4H4
| | - Matt Larouche
- Centre for Molecular Medicine and Therapeutics, Child and Family Research Institute, Department of Medical Genetics, University of British Columbia, 950 West 28th Avenue, Vancouver, BC, Canada V5Z 4H4
| | - Randy Glenn
- Centre for Molecular Medicine and Therapeutics, Child and Family Research Institute, Department of Medical Genetics, University of British Columbia, 950 West 28th Avenue, Vancouver, BC, Canada V5Z 4H4
| | - Dave Weeden
- Centre for Molecular Medicine and Therapeutics, Child and Family Research Institute, Department of Medical Genetics, University of British Columbia, 950 West 28th Avenue, Vancouver, BC, Canada V5Z 4H4
| | - Peter Zhang
- Centre for Molecular Medicine and Therapeutics, Child and Family Research Institute, Department of Medical Genetics, University of British Columbia, 950 West 28th Avenue, Vancouver, BC, Canada V5Z 4H4
| | - Kristin Hamre
- Department of Anatomy and Neurobiology, University of Tennessee Health Science Center, Memphis, TN, USA
| | - Michael Langston
- Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, TN, USA
| | - Charles Phillips
- Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, TN, USA
| | - Mingzhou Song
- Department of Computer Science, New Mexico State University, Las Cruces, NM, USA
| | - Zhengyu Ouyang
- Department of Computer Science, New Mexico State University, Las Cruces, NM, USA
| | | | | | | | - Yan Cui
- Department of Molecular Science, University of Tennessee Health Science Center, Memphis, TN, USA
| | - Kate Campbell
- Centre for Molecular Medicine and Therapeutics, Child and Family Research Institute, Department of Medical Genetics, University of British Columbia, 950 West 28th Avenue, Vancouver, BC, Canada V5Z 4H4
| | - Greg Ricker
- Department of Biology, Bowdoin College, Brunswick, ME, USA
| | - Carey Phillips
- Department of Biology, Bowdoin College, Brunswick, ME, USA
| | - Ramin Homayouni
- Bioinformatics Program, Department of Biology, University of Memphis, Memphis, TN, USA
| | - Dan Goldowitz
- Centre for Molecular Medicine and Therapeutics, Child and Family Research Institute, Department of Medical Genetics, University of British Columbia, 950 West 28th Avenue, Vancouver, BC, Canada V5Z 4H4.
| |
Collapse
|
28
|
Abstract
RNA-seq has emerged as the technology of choice to quantify gene expression. This technology is a convenient accurate tool to quantify diurnal changes in gene expression, gene discovery, differential use of promoters, and splice variants for all genes expressed in a single tissue. Thus, RNA-seq experiments provide sequence information and absolute expression values about transcripts in addition to relative quantification available with microarrays or qRT-PCR. The depth of information by sequencing requires careful assessment of RNA intactness and DNA contamination. Although the RNA-seq is comparatively recent, a standard analysis framework has emerged with the packages of Bowtie2, TopHat, and Cufflinks. With rising popularity of RNA-seq tools have become manageable for researchers without much bioinformatical knowledge or programming skills. Here, we present a workflow for a RNA-seq experiment from experimental planning to biological data extraction.
Collapse
Affiliation(s)
- Canan Külahoglu
- iGrad Plant, Heinrich Heine University Düsseldorf, Universitätsstrasse 1, 40225, Düsseldorf, Germany
| | | |
Collapse
|
29
|
Gaiteri C, Ding Y, French B, Tseng GC, Sibille E. Beyond modules and hubs: the potential of gene coexpression networks for investigating molecular mechanisms of complex brain disorders. GENES, BRAIN, AND BEHAVIOR 2014; 13:13-24. [PMID: 24320616 PMCID: PMC3896950 DOI: 10.1111/gbb.12106] [Citation(s) in RCA: 187] [Impact Index Per Article: 18.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/29/2013] [Revised: 09/25/2013] [Accepted: 11/10/2013] [Indexed: 12/12/2022]
Abstract
In a research environment dominated by reductionist approaches to brain disease mechanisms, gene network analysis provides a complementary framework in which to tackle the complex dysregulations that occur in neuropsychiatric and other neurological disorders. Gene-gene expression correlations are a common source of molecular networks because they can be extracted from high-dimensional disease data and encapsulate the activity of multiple regulatory systems. However, the analysis of gene coexpression patterns is often treated as a mechanistic black box, in which looming 'hub genes' direct cellular networks, and where other features are obscured. By examining the biophysical bases of coexpression and gene regulatory changes that occur in disease, recent studies suggest it is possible to use coexpression networks as a multi-omic screening procedure to generate novel hypotheses for disease mechanisms. Because technical processing steps can affect the outcome and interpretation of coexpression networks, we examine the assumptions and alternatives to common patterns of coexpression analysis and discuss additional topics such as acceptable datasets for coexpression analysis, the robust identification of modules, disease-related prioritization of genes and molecular systems and network meta-analysis. To accelerate coexpression research beyond modules and hubs, we highlight some emerging directions for coexpression network research that are especially relevant to complex brain disease, including the centrality-lethality relationship, integration with machine learning approaches and network pharmacology.
Collapse
Affiliation(s)
- Chris Gaiteri
- . Modeling, Analysis and Theory Group, Allen Institute for Brain Science, Seattle WA, USA
| | - Ying Ding
- . Carnegie Mellon-University of Pittsburgh PhD Program in Computational Biology, Pittsburgh, PA, USA
| | - Beverly French
- . Department of Psychiatry, University of Pittsburgh, Pittsburgh, PA, USA
| | - George C. Tseng
- . Department of Biostatistics, University of Pittsburgh, Pittsburgh, PA, USA
| | - Etienne Sibille
- . Department of Psychiatry, University of Pittsburgh, Pittsburgh, PA, USA
| |
Collapse
|
30
|
Iancu OD, Colville A, Darakjian P, Hitzemann R. Coexpression and cosplicing network approaches for the study of mammalian brain transcriptomes. INTERNATIONAL REVIEW OF NEUROBIOLOGY 2014; 116:73-93. [PMID: 25172472 DOI: 10.1016/b978-0-12-801105-8.00004-7] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Next-generation sequencing experiments have demonstrated great potential for transcriptome profiling. While transcriptome sequencing greatly increases the level of biological detail, system-level analysis of these high-dimensional datasets is becoming essential. We illustrate gene network approaches to the analysis of transcriptional data, with particular focus on the advantage of RNA-Seq technology compared to microarray platforms. We introduce a novel methodology for constructing cosplicing networks, based on distance measures combined with matrix correlations. We find that the cosplicing network is distinct and complementary to the coexpression network, although it shares the scale-free properties. In the cosplicing network, we find a set of novel hubs that have unique characteristics distinguishing them from coexpression hubs: they are heavily represented in neurobiological functional pathways and have strong overlap with markers of neurons and neuroglia, long-coding lengths, and high number of both exons and annotated transcripts. We also find that gene networks are plastic in the face of genetic and environmental pressures.
Collapse
Affiliation(s)
- Ovidiu Dan Iancu
- Department of Behavioral Neuroscience, Oregon Health & Science University, Portland, Oregon, USA.
| | - Alexandre Colville
- Department of Behavioral Neuroscience, Oregon Health & Science University, Portland, Oregon, USA
| | - Priscila Darakjian
- Department of Behavioral Neuroscience, Oregon Health & Science University, Portland, Oregon, USA
| | - Robert Hitzemann
- Department of Behavioral Neuroscience, Oregon Health & Science University, Portland, Oregon, USA; Research Service, Veterans Affairs Medical Center, Portland, Oregon, USA
| |
Collapse
|
31
|
Tong M, Li X, Wegener Parfrey L, Roth B, Ippoliti A, Wei B, Borneman J, McGovern DPB, Frank DN, Li E, Horvath S, Knight R, Braun J. A modular organization of the human intestinal mucosal microbiota and its association with inflammatory bowel disease. PLoS One 2013; 8:e80702. [PMID: 24260458 PMCID: PMC3834335 DOI: 10.1371/journal.pone.0080702] [Citation(s) in RCA: 127] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2013] [Accepted: 10/07/2013] [Indexed: 02/08/2023] Open
Abstract
Abnormalities of the intestinal microbiota are implicated in the pathogenesis of Crohn's disease (CD) and ulcerative colitis (UC), two spectra of inflammatory bowel disease (IBD). However, the high complexity and low inter-individual overlap of intestinal microbial composition are formidable barriers to identifying microbial taxa representing this dysbiosis. These difficulties might be overcome by an ecologic analytic strategy to identify modules of interacting bacteria (rather than individual bacteria) as quantitative reproducible features of microbial composition in normal and IBD mucosa. We sequenced 16S ribosomal RNA genes from 179 endoscopic lavage samples from different intestinal regions in 64 subjects (32 controls, 16 CD and 16 UC patients in clinical remission). CD and UC patients showed a reduction in phylogenetic diversity and shifts in microbial composition, comparable to previous studies using conventional mucosal biopsies. Analysis of weighted co-occurrence network revealed 5 microbial modules. These modules were unprecedented, as they were detectable in all individuals, and their composition and abundance was recapitulated in an independent, biopsy-based mucosal dataset 2 modules were associated with healthy, CD, or UC disease states. Imputed metagenome analysis indicated that these modules displayed distinct metabolic functionality, specifically the enrichment of oxidative response and glycan metabolism pathways relevant to host-pathogen interaction in the disease-associated modules. The highly preserved microbial modules accurately classified IBD status of individual patients during disease quiescence, suggesting that microbial dysbiosis in IBD may be an underlying disorder independent of disease activity. Microbial modules thus provide an integrative view of microbial ecology relevant to IBD.
Collapse
Affiliation(s)
- Maomeng Tong
- Department of Molecular and Medical Pharmacology, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, California, United States of America
| | - Xiaoxiao Li
- Cedars-Sinai F. Widjaja Inflammatory Bowel and Immunobiology Research Institute, Los Angeles, California, United States of America
| | - Laura Wegener Parfrey
- Department of Chemistry & Biochemistry, University of Colorado, Boulder, Colorado, United States of America
| | - Bennett Roth
- Department of Medicine, Division of Digestive Disease, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, California, United States of America
| | - Andrew Ippoliti
- Cedars-Sinai F. Widjaja Inflammatory Bowel and Immunobiology Research Institute, Los Angeles, California, United States of America
| | - Bo Wei
- Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, California, United States of America
| | - James Borneman
- Department of Plant Pathology and Microbiology, University of California Riverside, Riverside, California, United States of America
| | - Dermot P. B. McGovern
- Cedars-Sinai F. Widjaja Inflammatory Bowel and Immunobiology Research Institute, Los Angeles, California, United States of America
| | - Daniel N. Frank
- Division of Infectious Diseases, University of Colorado, School of Medicine, Aurora, Colorado, United States of America
- Union Council, Denver Microbiome Research Consortium (MiRC), University of Colorado, School of Medicine, Aurora, Colorado, United States of America
| | - Ellen Li
- Department of Medicine, Stony Brook University, Stony Brook, New York, United States of America
| | - Steve Horvath
- Department of Human Genetics and Biostatistics, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, California, United States of America
| | - Rob Knight
- Department of Chemistry & Biochemistry, University of Colorado, Boulder, Colorado, United States of America
- Howard Hughes Medical Institute, University of Colorado, Boulder, Colorado, United States of America;
| | - Jonathan Braun
- Department of Molecular and Medical Pharmacology, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, California, United States of America
- Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, California, United States of America
| |
Collapse
|
32
|
Giampieri R, Scartozzi M, Loretelli C, Piva F, Mandolesi A, Lezoche G, Prete MD, Bittoni A, Faloppi L, Bianconi M, Cecchini L, Guerrieri M, Bearzi I, Cascinu S. Cancer stem cell gene profile as predictor of relapse in high risk stage II and stage III, radically resected colon cancer patients. PLoS One 2013; 8:e72843. [PMID: 24023782 PMCID: PMC3762853 DOI: 10.1371/journal.pone.0072843] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2013] [Accepted: 07/16/2013] [Indexed: 12/11/2022] Open
Abstract
Clinical data indicate that prognostic stratification of radically resected colorectal cancer based on disease stage only may not be always be adequate. Preclinical findings suggest that cancer stem cells may influence the biological behaviour of colorectal cancer independently from stage: objective of the study was to assess whether a panel of stemness markers were correlated with clinical outcome in resected stage II and III colon cancer patients. A panel of 66 markers of stemness were analysed and thus patients were divided into two groups (A and B) with most patients clustering in a manner consistent with different time to relapse by using a statistical algorithm. A total of 62 patients were analysed. Thirty-six (58%) relapsed during the follow-up period (range 1.63–86.5 months). Twelve (19%) and 50 (81%) patients were allocated into group A and B, respectively. A significantly different median relapse-free survival was observed between the 2 groups (22.18 vs 42.85 months, p = 0.0296). Among of all genes tested, those with the higher “weight” in determining different prognosis were CD44, ALCAM, DTX2, HSPA9, CCNA2, PDX1, MYST1, COL1A1 and ABCG2. This analysis supports the idea that, other than stage, biological variables, such as expression levels of colon cancer stem cell genes, may be relevant in determining an increased risk of relapse in resected colorectal cancer patients.
Collapse
Affiliation(s)
- Riccardo Giampieri
- Department of Medical Oncology, AO Ospedali Riuniti-UNIVPM, Ancona, Italy
| | - Mario Scartozzi
- Department of Medical Oncology, AO Ospedali Riuniti-UNIVPM, Ancona, Italy
- * E-mail:
| | - Cristian Loretelli
- Department of Medical Oncology, AO Ospedali Riuniti-UNIVPM, Ancona, Italy
| | | | | | | | - Michela Del Prete
- Department of Medical Oncology, AO Ospedali Riuniti-UNIVPM, Ancona, Italy
| | - Alessandro Bittoni
- Department of Medical Oncology, AO Ospedali Riuniti-UNIVPM, Ancona, Italy
| | - Luca Faloppi
- Department of Medical Oncology, AO Ospedali Riuniti-UNIVPM, Ancona, Italy
| | | | - Luca Cecchini
- Department of Medical Oncology, AO Ospedali Riuniti-UNIVPM, Ancona, Italy
| | - Mario Guerrieri
- Institute of Surgery, AO ospedali Riuniti-UNIVPM, Ancona, Italy
| | - Italo Bearzi
- Institute of Pathology, AO ospedali Riuniti-UNIVPM, Ancona, Italy
| | - Stefano Cascinu
- Department of Medical Oncology, AO Ospedali Riuniti-UNIVPM, Ancona, Italy
| |
Collapse
|
33
|
McHardy IH, Goudarzi M, Tong M, Ruegger PM, Schwager E, Weger JR, Graeber TG, Sonnenburg JL, Horvath S, Huttenhower C, McGovern DPB, Fornace AJ, Borneman J, Braun J. Integrative analysis of the microbiome and metabolome of the human intestinal mucosal surface reveals exquisite inter-relationships. MICROBIOME 2013; 1:17. [PMID: 24450808 PMCID: PMC3971612 DOI: 10.1186/2049-2618-1-17] [Citation(s) in RCA: 198] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/31/2013] [Accepted: 05/12/2013] [Indexed: 05/10/2023]
Abstract
BACKGROUND Consistent compositional shifts in the gut microbiota are observed in IBD and other chronic intestinal disorders and may contribute to pathogenesis. The identities of microbial biomolecular mechanisms and metabolic products responsible for disease phenotypes remain to be determined, as do the means by which such microbial functions may be therapeutically modified. RESULTS The composition of the microbiota and metabolites in gut microbiome samples in 47 subjects were determined. Samples were obtained by endoscopic mucosal lavage from the cecum and sigmoid colon regions, and each sample was sequenced using the 16S rRNA gene V4 region (Illumina-HiSeq 2000 platform) and assessed by UPLC mass spectroscopy. Spearman correlations were used to identify widespread, statistically significant microbial-metabolite relationships. Metagenomes for identified microbial OTUs were imputed using PICRUSt, and KEGG metabolic pathway modules for imputed genes were assigned using HUMAnN. The resulting metabolic pathway abundances were mostly concordant with metabolite data. Analysis of the metabolome-driven distribution of OTU phylogeny and function revealed clusters of clades that were both metabolically and metagenomically similar. CONCLUSIONS The results suggest that microbes are syntropic with mucosal metabolome composition and therefore may be the source of and/or dependent upon gut epithelial metabolites. The consistent relationship between inferred metagenomic function and assayed metabolites suggests that metagenomic composition is predictive to a reasonable degree of microbial community metabolite pools. The finding that certain metabolites strongly correlate with microbial community structure raises the possibility of targeting metabolites for monitoring and/or therapeutically manipulating microbial community function in IBD and other chronic diseases.
Collapse
Affiliation(s)
- Ian H McHardy
- Pathology and Laboratory Medicine UCLA, Los Angeles, CA, USA
| | - Maryam Goudarzi
- Biochemistry and Molecular and Cellular Biology, Georgetown University, Washington, DC, USA
| | - Maomeng Tong
- Molecular and Medical Pharmacology, UCLA, Los Angeles, CA, USA
| | | | | | - John R Weger
- Plant Pathology, UC Riverside, Riverside, CA, USA
| | | | | | | | | | - Dermot PB McGovern
- The F. Widjaja Family Foundation Inflammatory Bowel and Immunobiology Research Institute, Cedar's Sinai Medical Center, Los Angeles, CA, USA
| | - Albert J Fornace
- Biochemistry and Molecular and Cellular Biology, Georgetown University, Washington, DC, USA
| | | | - Jonathan Braun
- Pathology and Laboratory Medicine UCLA, Los Angeles, CA, USA
| |
Collapse
|
34
|
Genetic analysis of Mps3 SUN domain mutants in Saccharomyces cerevisiae reveals an interaction with the SUN-like protein Slp1. G3-GENES GENOMES GENETICS 2012; 2:1703-18. [PMID: 23275891 PMCID: PMC3516490 DOI: 10.1534/g3.112.004614] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/04/2012] [Accepted: 10/26/2012] [Indexed: 01/16/2023]
Abstract
In virtually all eukaryotic cells, protein bridges formed by the conserved inner nuclear membrane SUN (for Sad1-UNC-84) domain-containing proteins and their outer nuclear membrane binding partners span the nuclear envelope (NE) to connect the nucleoplasm and cytoplasm. These linkages are important for chromosome movements within the nucleus during meiotic prophase and are essential for nuclear migration and centrosome attachment to the NE. In Saccharomyces cerevisiae, MPS3 encodes the sole SUN protein. Deletion of MPS3 or the conserved SUN domain is lethal in three different genetic backgrounds. Mutations in the SUN domain result in defects in duplication of the spindle pole body, the yeast centrosome-equivalent organelle. A genome-wide screen for mutants that exhibited synthetic fitness defects in combination with mps3 SUN domain mutants yielded a large number of hits in components of the spindle apparatus and the spindle checkpoint. Mutants in lipid metabolic processes and membrane organization also exacerbated the growth defects of mps3 SUN domain mutants, pointing to a role for Mps3 in nuclear membrane organization. Deletion of SLP1 or YER140W/EMP65 (for ER membrane protein of 65 kDa) aggravated growth of mps3 SUN domain mutants. Slp1 and Emp65 form an ER-membrane associated protein complex that is not required directly for spindle pole body duplication or spindle assembly. Rather, Slp1 is involved in Mps3 localization to the NE.
Collapse
|