1
|
Heaton EC, Seo EH, Butkovich LM, Yount ST, Gourley SL. Control of goal-directed and inflexible actions by dorsal striatal melanocortin systems, in coordination with the central nucleus of the amygdala. Prog Neurobiol 2024; 238:102629. [PMID: 38763506 PMCID: PMC11198735 DOI: 10.1016/j.pneurobio.2024.102629] [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: 12/19/2023] [Revised: 04/15/2024] [Accepted: 05/10/2024] [Indexed: 05/21/2024]
Abstract
The dorsomedial striatum (DMS) is associated with flexible goal seeking, as opposed to routinized habits. Whether local mechanisms brake this function, for instance when habits may be adaptive, is incompletely understood. We find that a sub-population of dopamine D1 receptor-containing striatal neurons express the melanocortin-4 receptor (MC4R) for α-melanocyte stimulating hormone. These neurons within the DMS are necessary and sufficient for controlling the capacity of mice to flexibly adjust actions based on the likelihood that they will be rewarded. In investigating MC4R function, we found that it suppresses immediate-early gene levels in the DMS and concurrently, flexible goal seeking. MC4R+ neurons receive input from the central nucleus of the amygdala, and behavioral experiments indicate that they are functionally integrated into an amygdalo-striatal circuit that suppresses action flexibility in favor of routine. Publicly available spatial transcriptomics datasets were analyzed for gene transcript correlates of Mc4r expression across the striatal subregions, revealing considerable co-variation in dorsal structures. This insight led to the discovery that the function of MC4R in the dorsolateral striatum complements that in the DMS, in this case suppressing habit-like behavior. Altogether, our findings suggest that striatal MC4R controls the capacity for goal-directed and inflexible actions alike.
Collapse
Affiliation(s)
- Elizabeth C Heaton
- Graduate Program in Neuroscience, Emory University, United States; Emory National Primate Research Center, Departments of Pediatrics and Psychiatry and Behavioral Sciences, Emory University School of Medicine, United States
| | - Esther H Seo
- Emory National Primate Research Center, Departments of Pediatrics and Psychiatry and Behavioral Sciences, Emory University School of Medicine, United States
| | - Laura M Butkovich
- Emory National Primate Research Center, Departments of Pediatrics and Psychiatry and Behavioral Sciences, Emory University School of Medicine, United States
| | - Sophie T Yount
- Emory National Primate Research Center, Departments of Pediatrics and Psychiatry and Behavioral Sciences, Emory University School of Medicine, United States; Graduate Program in Molecular and Systems Pharmacology, Emory University, United States
| | - Shannon L Gourley
- Graduate Program in Neuroscience, Emory University, United States; Emory National Primate Research Center, Departments of Pediatrics and Psychiatry and Behavioral Sciences, Emory University School of Medicine, United States; Graduate Program in Molecular and Systems Pharmacology, Emory University, United States; Children's Healthcare of Atlanta, Atlanta, GA, United States.
| |
Collapse
|
2
|
Branciamore S, Gogoshin G, Rodin AS, Myers AJ. Changes in expression of VGF, SPECC1L, HLA-DRA and RANBP3L act with APOE E4 to alter risk for late onset Alzheimer's disease. Sci Rep 2024; 14:14954. [PMID: 38942763 PMCID: PMC11213882 DOI: 10.1038/s41598-024-65010-7] [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: 11/28/2023] [Accepted: 06/16/2024] [Indexed: 06/30/2024] Open
Abstract
While there are currently over 40 replicated genes with mapped risk alleles for Late Onset Alzheimer's disease (LOAD), the Apolipoprotein E locus E4 haplotype is still the biggest driver of risk, with odds ratios for neuropathologically confirmed E44 carriers exceeding 30 (95% confidence interval 16.59-58.75). We sought to address whether the APOE E4 haplotype modifies expression globally through networks of expression to increase LOAD risk. We have used the Human Brainome data to build expression networks comparing APOE E4 carriers to non-carriers using scalable mixed-datatypes Bayesian network (BN) modeling. We have found that VGF had the greatest explanatory weight. High expression of VGF is a protective signal, even on the background of APOE E4 alleles. LOAD risk signals, considering an APOE background, include high levels of SPECC1L, HLA-DRA and RANBP3L. Our findings nominate several new transcripts, taking a combined approach to network building including known LOAD risk loci.
Collapse
Affiliation(s)
- Sergio Branciamore
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, CA, 91010, USA
| | - Grigoriy Gogoshin
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, CA, 91010, USA
| | - Andrei S Rodin
- Department of Computational and Quantitative Medicine, Beckman Research Institute of the City of Hope, Duarte, CA, 91010, USA.
| | - Amanda J Myers
- Department of Cell Biology, University of Miami Miller School of Medicine, Miami, FL, 33136, USA.
- Institute for Data Science and Computing, University of Miami Miller School of Medicine, Miami, FL, 33136, USA.
- Interdepartmental Program in Neuroscience, University of Miami Miller School of Medicine, Miami, FL, 33136, USA.
- Interdepartmental Program in Human Genetics and Genomics, University of Miami Miller School of Medicine, Miami, FL, 33136, USA.
| |
Collapse
|
3
|
Tanei T, Seno S, Sota Y, Hatano T, Kitahara Y, Abe K, Masunaga N, Tsukabe M, Yoshinami T, Miyake T, Shimoda M, Matsuda H, Shimazu K. High HER2 Intratumoral Heterogeneity Is a Predictive Factor for Poor Prognosis in Early-Stage and Locally Advanced HER2-Positive Breast Cancer. Cancers (Basel) 2024; 16:1062. [PMID: 38473420 DOI: 10.3390/cancers16051062] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Revised: 02/20/2024] [Accepted: 02/26/2024] [Indexed: 03/14/2024] Open
Abstract
PURPOSE Breast cancer tumors frequently have intratumoral heterogeneity (ITH). Tumors with high ITH cause therapeutic resistance and have human epidermal growth factor receptor 2 (HER2) heterogeneity in response to HER2-targeted therapies. This study aimed to investigate whether high HER2 heterogeneity levels were clinically related to a poor prognosis for HER2-targeted adjuvant therapy resistance in primary breast cancers. METHODS This study included patients with primary breast cancer (n = 251) treated with adjuvant HER2-targeted therapies. HER2 heterogeneity was manifested by the shape of HER2 fluorescence in situ hybridization amplification (FISH) distributed histograms with the HER2 gene copy number within a tumor sample. Each tumor was classified into a biphasic grade graph (high heterogeneity [HH]) group or a monophasic grade graph (low heterogeneity [LH]) group based on heterogeneity. Both groups were evaluated for disease-free survival (DFS) and overall survival (OS) for a median of ten years of annual follow-up. RESULTS Of 251 patients with HER2-positive breast cancer, 46 (18.3%) and 205 (81.7%) were classified into the HH and LH groups, respectively. The HH group had more distant metastases and a poorer prognosis than the LH group (DFS: p < 0.001 (HH:63% vs. LH:91% at 10 years) and for the OS: p = 0.012 (HH:78% vs. LH:95% at 10 years). CONCLUSIONS High HER2 heterogeneity is a poor prognostic factor in patients with HER2-positive breast cancer. A novel approach to heterogeneity, which is manifested by the shape of HER2 FISH distributions, might be clinically useful in the prognosis prediction of patients after HER2 adjuvant therapy.
Collapse
Affiliation(s)
- Tomonori Tanei
- Department of Breast and Endocrine Surgery, Graduate School of Medicine, Osaka University, 2-2-E10 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Shigeto Seno
- Department of Bioinformatic Engineering, Graduate School of Information Science and Technology, Osaka University, 1-5 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Yoshiaki Sota
- Department of Breast and Endocrine Surgery, Graduate School of Medicine, Osaka University, 2-2-E10 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Takaaki Hatano
- Department of Breast and Endocrine Surgery, Graduate School of Medicine, Osaka University, 2-2-E10 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Yuri Kitahara
- Department of Breast and Endocrine Surgery, Graduate School of Medicine, Osaka University, 2-2-E10 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Kaori Abe
- Department of Breast and Endocrine Surgery, Graduate School of Medicine, Osaka University, 2-2-E10 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Nanae Masunaga
- Department of Breast and Endocrine Surgery, Graduate School of Medicine, Osaka University, 2-2-E10 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Masami Tsukabe
- Department of Breast and Endocrine Surgery, Graduate School of Medicine, Osaka University, 2-2-E10 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Tetsuhiro Yoshinami
- Department of Breast and Endocrine Surgery, Graduate School of Medicine, Osaka University, 2-2-E10 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Tomohiro Miyake
- Department of Breast and Endocrine Surgery, Graduate School of Medicine, Osaka University, 2-2-E10 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Masafumi Shimoda
- Department of Breast and Endocrine Surgery, Graduate School of Medicine, Osaka University, 2-2-E10 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Hideo Matsuda
- Department of Bioinformatic Engineering, Graduate School of Information Science and Technology, Osaka University, 1-5 Yamadaoka, Suita 565-0871, Osaka, Japan
| | - Kenzo Shimazu
- Department of Breast and Endocrine Surgery, Graduate School of Medicine, Osaka University, 2-2-E10 Yamadaoka, Suita 565-0871, Osaka, Japan
| |
Collapse
|
4
|
Xu D, Li P, Zhang C, Shen Y, Cai J, Wei Q, Cao M, Xu Z, Wu D, Wang H, Bi X, Wang B, Li K. Development of an m6A-Related lncRNAs Signature Predicts Tumor Stemness and Prognosis for Low-Grade Glioma Patients. Stem Cells Int 2024; 2024:2062283. [PMID: 38229597 PMCID: PMC10791469 DOI: 10.1155/2024/2062283] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2023] [Revised: 10/25/2023] [Accepted: 11/22/2023] [Indexed: 01/18/2024] Open
Abstract
Background Growing evidence has revealed that m6A modification of long noncoding RNAs (lncRNAs) dynamically controls tumor stemness and tumorigenesis-related processes. However, the prognostic significance of m6A-related lncRNAs and their associations with stemness in low-grade glioma (LGG) remain to be clarified. Methods A multicenter transcriptome analysis of lncRNA expression in 1,247 LGG samples was performed in this study. The stemness landscape of LGG tumors was presented and associations with clinical features were revealed. The m6A-related lncRNAs were identified between stemness groups and were further prioritized via least absolute shrinkage and selection operator Cox regression analysis. A risk score model based on m6A-related lncRNAs was constructed and validated in external LGG datasets. Results Based on the expression of LINC02984, PFKP-DT, and CRNDE, a risk model and nomogram were constructed; they successfully predicted the survival of patients and were extended to external datasets. Significant correlations were observed between the risk score and tumor stemness. Moreover, patients in different risk groups exhibited distinct tumor immune microenvironments and immune signatures. We finally provided several potential compounds suitable for specific risk groups, which may aid in LGG treatment. Conclusions This novel signature presents noteworthy value in the prediction of prognosis and stemness status for LGG patients and will foster future research on the development of clinical regimens.
Collapse
Affiliation(s)
- Dahua Xu
- Key Laboratory of Tropical Translational Medicine of Ministry of Education, College of Biomedical Information and Engineering, Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 571199, China
| | - Peihu Li
- Key Laboratory of Tropical Translational Medicine of Ministry of Education, College of Biomedical Information and Engineering, Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 571199, China
| | - Chunrui Zhang
- Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100020, China
| | - Yutong Shen
- Key Laboratory of Tropical Translational Medicine of Ministry of Education, College of Biomedical Information and Engineering, Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 571199, China
| | - Jiale Cai
- Key Laboratory of Tropical Translational Medicine of Ministry of Education, College of Biomedical Information and Engineering, Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 571199, China
| | - Qingchen Wei
- Key Laboratory of Tropical Translational Medicine of Ministry of Education, College of Biomedical Information and Engineering, Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 571199, China
| | - Meng Cao
- Key Laboratory of Tropical Translational Medicine of Ministry of Education, College of Biomedical Information and Engineering, Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 571199, China
| | - Zhizhou Xu
- Key Laboratory of Tropical Translational Medicine of Ministry of Education, College of Biomedical Information and Engineering, Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 571199, China
| | - Deng Wu
- School of Life Sciences, Faculty of Science, The Chinese University of Hong Kong, Hong Kong 999077, China
| | - Hong Wang
- Key Laboratory of Tropical Translational Medicine of Ministry of Education, College of Biomedical Information and Engineering, Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 571199, China
| | - Xiaoman Bi
- Key Laboratory of Tropical Translational Medicine of Ministry of Education, College of Biomedical Information and Engineering, Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 571199, China
| | - Bo Wang
- Key Laboratory of Tropical Translational Medicine of Ministry of Education, College of Biomedical Information and Engineering, Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 571199, China
| | - Kongning Li
- Key Laboratory of Tropical Translational Medicine of Ministry of Education, College of Biomedical Information and Engineering, Hainan General Hospital, Hainan Affiliated Hospital of Hainan Medical University, Haikou 571199, China
| |
Collapse
|
5
|
Danino R, Nachman I, Sharan R. Batch correction of single-cell sequencing data via an autoencoder architecture. BIOINFORMATICS ADVANCES 2023; 4:vbad186. [PMID: 38213820 PMCID: PMC10781938 DOI: 10.1093/bioadv/vbad186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/29/2023] [Revised: 12/09/2023] [Accepted: 12/17/2023] [Indexed: 01/13/2024]
Abstract
Motivation Technical differences between gene expression sequencing experiments can cause variations in the data in the form of batch effect biases. These do not represent true biological variations between samples and can lead to false conclusions or hinder the ability to integrate multiple datasets. Since there is a growing need for the joint analysis of single-cell sequencing datasets from different sources, there is also a need to correct the resulting batch effects while maintaining the true biological variations in the data. Results We developed a semi-supervised deep learning architecture called Autoencoder-based Batch Correction (ABC) for integrating single-cell sequencing datasets. Our method removes batch effects through a guided process of data compression using supervised cell type classifier branches for biological signal retention. It aligns the different batches using an adversarial training approach. We comprehensively evaluate the performance of our method using four single-cell sequencing datasets and multiple measures for batch effect removal and biological variation conservation. ABC outperforms 10 state-of-the-art methods for this task including Seurat, scGen, ComBat, scanorama, scVI, scANVI, AutoClass, Harmony, scDREAMER, and CLEAR, correcting various types of batch effects while preserving intricate biological variations.
Collapse
Affiliation(s)
- Reut Danino
- Blavatnik School of Computer Science, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Iftach Nachman
- School of Neurobiology, Biochemistry and Biophysics, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Roded Sharan
- Blavatnik School of Computer Science, Tel Aviv University, Tel Aviv 6997801, Israel
| |
Collapse
|
6
|
Branciamore S, Gogoshin G, Rodin AS, Myers AJ. The Human Brainome: changes in expression of VGF, SPECC1L, HLA-DRA and RANBP3L act with APOE E4 to alter risk for late onset Alzheimer's disease. RESEARCH SQUARE 2023:rs.3.rs-3678057. [PMID: 38168398 PMCID: PMC10760217 DOI: 10.21203/rs.3.rs-3678057/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
While there are currently over 40 replicated genes with mapped risk alleles for Late Onset Alzheimer's disease (LOAD), the Apolipoprotein E locus E4 haplotype is still the biggest driver of risk, with odds ratios for neuropathologically confirmed E44 carriers exceeding 30 (95% confidence interval 16.59-58.75). We sought to address whether the APOE E4 haplotype modifies expression globally through networks of expression to increase LOAD risk. We have used the Human Brainome data to build expression networks comparing APOE E4 carriers to non-carriers using scalable mixed-datatypes Bayesian network (BN) modeling. We have found that VGF had the greatest explanatory weight. High expression of VGF is a protective signal, even on the background of APOE E4 alleles. LOAD risk signals, considering an APOE background, include high levels of SPECC1L, HLA-DRA and RANBP3L. Our findings nominate several new transcripts, taking a combined approach to network building including known LOAD risk loci.
Collapse
|
7
|
Kim D, Tran A, Kim HJ, Lin Y, Yang JYH, Yang P. Gene regulatory network reconstruction: harnessing the power of single-cell multi-omic data. NPJ Syst Biol Appl 2023; 9:51. [PMID: 37857632 PMCID: PMC10587078 DOI: 10.1038/s41540-023-00312-6] [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: 08/14/2023] [Accepted: 10/02/2023] [Indexed: 10/21/2023] Open
Abstract
Inferring gene regulatory networks (GRNs) is a fundamental challenge in biology that aims to unravel the complex relationships between genes and their regulators. Deciphering these networks plays a critical role in understanding the underlying regulatory crosstalk that drives many cellular processes and diseases. Recent advances in sequencing technology have led to the development of state-of-the-art GRN inference methods that exploit matched single-cell multi-omic data. By employing diverse mathematical and statistical methodologies, these methods aim to reconstruct more comprehensive and precise gene regulatory networks. In this review, we give a brief overview on the statistical and methodological foundations commonly used in GRN inference methods. We then compare and contrast the latest state-of-the-art GRN inference methods for single-cell matched multi-omics data, and discuss their assumptions, limitations and opportunities. Finally, we discuss the challenges and future directions that hold promise for further advancements in this rapidly developing field.
Collapse
Affiliation(s)
- Daniel Kim
- School of Mathematics and Statistics, University of Sydney, Camperdown, NSW, Australia
- Computational Systems Biology Unit, Children's Medical Research Institute, University of Sydney, Camperdown, NSW, Australia
- Sydney Precision Data Science Centre, University of Sydney, Camperdown, NSW, Australia
| | - Andy Tran
- School of Mathematics and Statistics, University of Sydney, Camperdown, NSW, Australia
- Sydney Precision Data Science Centre, University of Sydney, Camperdown, NSW, Australia
- Charles Perkins Centre, University of Sydney, Camperdown, NSW, Australia
| | - Hani Jieun Kim
- Computational Systems Biology Unit, Children's Medical Research Institute, University of Sydney, Camperdown, NSW, Australia
- Sydney Precision Data Science Centre, University of Sydney, Camperdown, NSW, Australia
| | - Yingxin Lin
- School of Mathematics and Statistics, University of Sydney, Camperdown, NSW, Australia
- Sydney Precision Data Science Centre, University of Sydney, Camperdown, NSW, Australia
- Charles Perkins Centre, University of Sydney, Camperdown, NSW, Australia
| | - Jean Yee Hwa Yang
- School of Mathematics and Statistics, University of Sydney, Camperdown, NSW, Australia.
- Sydney Precision Data Science Centre, University of Sydney, Camperdown, NSW, Australia.
- Charles Perkins Centre, University of Sydney, Camperdown, NSW, Australia.
| | - Pengyi Yang
- School of Mathematics and Statistics, University of Sydney, Camperdown, NSW, Australia.
- Computational Systems Biology Unit, Children's Medical Research Institute, University of Sydney, Camperdown, NSW, Australia.
- Sydney Precision Data Science Centre, University of Sydney, Camperdown, NSW, Australia.
- Charles Perkins Centre, University of Sydney, Camperdown, NSW, Australia.
| |
Collapse
|
8
|
Tschernoster N, Erger F, Kohl S, Reusch B, Wenzel A, Walsh S, Thiele H, Becker C, Franitza M, Bartram MP, Kömhoff M, Schumacher L, Kukat C, Borodina T, Quedenau C, Nürnberg P, Rinschen MM, Driller JH, Pedersen BP, Schlingmann KP, Hüttel B, Bockenhauer D, Beck B, Altmüller J. Long-read sequencing identifies a common transposition haplotype predisposing for CLCNKB deletions. Genome Med 2023; 15:62. [PMID: 37612755 PMCID: PMC10464140 DOI: 10.1186/s13073-023-01215-1] [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: 03/15/2023] [Accepted: 07/27/2023] [Indexed: 08/25/2023] Open
Abstract
BACKGROUND Long-read sequencing is increasingly used to uncover structural variants in the human genome, both functionally neutral and deleterious. Structural variants occur more frequently in regions with a high homology or repetitive segments, and one rearrangement may predispose to additional events. Bartter syndrome type 3 (BS 3) is a monogenic tubulopathy caused by deleterious variants in the chloride channel gene CLCNKB, a high proportion of these being large gene deletions. Multiplex ligation-dependent probe amplification, the current diagnostic gold standard for this type of mutation, will indicate a simple homozygous gene deletion in biallelic deletion carriers. However, since the phenotypic spectrum of BS 3 is broad even among biallelic deletion carriers, we undertook a more detailed analysis of precise breakpoint regions and genomic structure. METHODS Structural variants in 32 BS 3 patients from 29 families and one BS4b patient with CLCNKB deletions were investigated using long-read and synthetic long-read sequencing, as well as targeted long-read sequencing approaches. RESULTS We report a ~3 kb duplication of 3'-UTR CLCNKB material transposed to the corresponding locus of the neighbouring CLCNKA gene, also found on ~50 % of alleles in healthy control individuals. This previously unknown common haplotype is significantly enriched in our cohort of patients with CLCNKB deletions (45 of 51 alleles with haplotype information, 2.2 kb and 3.0 kb transposition taken together, p=9.16×10-9). Breakpoint coordinates for the CLCNKB deletion were identifiable in 28 patients, with three being compound heterozygous. In total, eight different alleles were found, one of them a complex rearrangement with three breakpoint regions. Two patients had different CLCNKA/CLCNKB hybrid genes encoding a predicted CLCNKA/CLCNKB hybrid protein with likely residual function. CONCLUSIONS The presence of multiple different deletion alleles in our cohort suggests that large CLCNKB gene deletions originated from many independently recurring genomic events clustered in a few hot spots. The uncovered associated sequence transposition haplotype apparently predisposes to these additional events. The spectrum of CLCNKB deletion alleles is broader than expected and likely still incomplete, but represents an obvious candidate for future genotype/phenotype association studies. We suggest a sensitive and cost-efficient approach, consisting of indirect sequence capture and long-read sequencing, to analyse disease-relevant structural variant hotspots in general.
Collapse
Affiliation(s)
- Nikolai Tschernoster
- Cologne Center for Genomics (CCG), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany
- Institute of Human Genetics, Faculty of Medicine and University Hospital Cologne, University of Cologne, Kerpener Str. 34, 50931, Cologne, Germany
- Center for Molecular Medicine Cologne (CMMC), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany
| | - Florian Erger
- Institute of Human Genetics, Faculty of Medicine and University Hospital Cologne, University of Cologne, Kerpener Str. 34, 50931, Cologne, Germany
- Center for Molecular Medicine Cologne (CMMC), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany
| | - Stefan Kohl
- Department of Pediatrics, Cologne Children's Hospital, Cologne, Germany
| | - Björn Reusch
- Institute of Human Genetics, Faculty of Medicine and University Hospital Cologne, University of Cologne, Kerpener Str. 34, 50931, Cologne, Germany
- Center for Molecular Medicine Cologne (CMMC), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany
| | - Andrea Wenzel
- Institute of Human Genetics, Faculty of Medicine and University Hospital Cologne, University of Cologne, Kerpener Str. 34, 50931, Cologne, Germany
- Center for Molecular Medicine Cologne (CMMC), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany
| | - Stephen Walsh
- Department of Renal Medicine, UCL, University College London, London, UK
| | - Holger Thiele
- Cologne Center for Genomics (CCG), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany
| | - Christian Becker
- Cologne Center for Genomics (CCG), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany
| | - Marek Franitza
- Cologne Center for Genomics (CCG), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany
| | - Malte P Bartram
- Center for Molecular Medicine Cologne (CMMC), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany
- Department II of Internal Medicine, University of Cologne, Cologne, Germany
| | - Martin Kömhoff
- Department of Pediatrics, University Marburg, Marburg, Germany
| | - Lena Schumacher
- FACS & Imaging Core Facility, Max Planck Institute for Biology of Ageing, Cologne, Germany
| | - Christian Kukat
- FACS & Imaging Core Facility, Max Planck Institute for Biology of Ageing, Cologne, Germany
| | - Tatiana Borodina
- Max Delbrück Center for Molecular Medicine in the Helmholtz Association (MDC), Hannoversche Straße 28, 10115, Berlin, Germany
| | - Claudia Quedenau
- Max Delbrück Center for Molecular Medicine in the Helmholtz Association (MDC), Hannoversche Straße 28, 10115, Berlin, Germany
| | - Peter Nürnberg
- Cologne Center for Genomics (CCG), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany
- Center for Molecular Medicine Cologne (CMMC), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany
| | - Markus M Rinschen
- Department of Biomedicine, Aarhus University, Aarhus, Denmark
- Aarhus Institute of Advanced Studies, Aarhus University, Aarhus, Denmark
- Department III of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Jan H Driller
- Department of Molecular Biology and Genetics, Aarhus University, Universitetsbyen 81, DK-8000, Aarhus C, Denmark
| | - Bjørn P Pedersen
- Department of Molecular Biology and Genetics, Aarhus University, Universitetsbyen 81, DK-8000, Aarhus C, Denmark
| | - Karl P Schlingmann
- Department of General Pediatrics, University Children's Hospital, Münster, Germany
| | - Bruno Hüttel
- Max Planck Genome-Centre Cologne, Max Planck Institute for Plant Breeding Research, Cologne, Germany
| | - Detlef Bockenhauer
- Department of Renal Medicine, UCL, University College London, London, UK
- Great Ormond Street Hospital for Children, NHS Foundation Trust, London, UK
| | - Bodo Beck
- Institute of Human Genetics, Faculty of Medicine and University Hospital Cologne, University of Cologne, Kerpener Str. 34, 50931, Cologne, Germany.
- Center for Molecular Medicine Cologne (CMMC), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany.
| | - Janine Altmüller
- Center for Molecular Medicine Cologne (CMMC), University of Cologne, Faculty of Medicine and University Hospital Cologne, Cologne, Germany.
- Max Delbrück Center for Molecular Medicine in the Helmholtz Association (MDC), Hannoversche Straße 28, 10115, Berlin, Germany.
- Berlin Institute of Health at Charité - Universitätsmedizin Berlin, Core Facility Genomics, Berlin, Germany.
| |
Collapse
|
9
|
Matsumoto Y, Nakamura S. Rapid and Comprehensive Identification of Nontuberculous Mycobacteria. Methods Mol Biol 2023; 2632:247-255. [PMID: 36781733 DOI: 10.1007/978-1-0716-2996-3_17] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023]
Abstract
Next-generation sequencing is a powerful tool to accurately identify pathogens. The MinION sequencer is best suited for the rapid identification of bacterial species due to its real-time sequence output. In this chapter, we introduce a method to identify nontuberculous mycobacteria (NTM) in one sequencing analysis from culture isolates using the MinION sequencer. NTM disease is now recognized as a growing global health concern due to its increasing incidence and prevalence. There are over 200 NTM species, of which the major pathogens are further classified into many subspecies showing different antibiotic susceptibilities. Therefore, identifying the pathogens at the subspecies level of NTM is necessary to select an appropriate treatment regimen. The protocol described here includes DNA extraction by lysis using silica beads, library preparation, sequencing by the MinION sequencer, and analysis of multilocus sequence typing using the software "mlstverse" and enables rapid and comprehensive identification of 175 species of NTM at the subspecies level with high sensitivity and accuracy.
Collapse
Affiliation(s)
- Yuki Matsumoto
- Department of Infection Metagenomics, Research Institute for Microbial Diseases, Osaka University, Osaka, Japan
| | - Shota Nakamura
- Department of Infection Metagenomics, Research Institute for Microbial Diseases, Osaka University, Osaka, Japan.
- Center for Infectious Disease Education and Research, Osaka University, Osaka, Japan.
| |
Collapse
|
10
|
Shchepina OA, Menshanov PN. Neuron-Glia-Ratio-Like Approach Evidenced for Limited Variability and In-Aggregate Circadian Shifts in Cortical Cell-Specific Transcriptomes. J Mol Neurosci 2023; 73:159-170. [PMID: 36745298 DOI: 10.1007/s12031-023-02103-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2022] [Accepted: 01/22/2023] [Indexed: 02/07/2023]
Abstract
Regardless of shifts in levels of individual transcripts, it remains elusive whether natural variability in cell-specific transcriptomes within the cerebral cortex is limited in aggregate. It is also unclear whether cortical cell-specific transcriptomes might change dynamically in absence of cell number changes. Total variation in neuron- and glia-specific in-aggregate transcriptomes could be identified in a model-free way via glia-neuron ratio approach, by univariate median-to-median ratios comparing integral levels of cell-specific transcripts within a tissue sample. When deleterious, regenerative or developmental events affecting cortical cell numbers were subtle, median-to-median ratios demonstrated within-group variability not exceeding <20-25% in most cases. These levels of total variability might be explained in part by limited (~5-10%) circadian and stress-induced shifts in cell-specific cortical transcriptomes. Relevant in-aggregate transcriptomic alterations were identified after shifts in cell numbers induced by well-validated deleterious events including ischemia, traumatic injury, microglia's activation/depletion or specific mutations. Cortical median-to-median ratios also follow naturally occurring changes in the numbers of excitatory, inhibitory neurons and glial cells during perinatal brain development. These findings characterize cortical cell-specific transcriptomes as subjects to circadian shifts and lifetime events, urging the importance of reporting full details on an origin of any transcriptomic sample collected in vivo.
Collapse
Affiliation(s)
- Olesya A Shchepina
- Ermine Educational Center, Novosibirsk State University, Novosibirsk, Novosibirsk Region, 630117, Russian Federation.,Higher College of Informatics, Novosibirsk State University, Novosibirsk, Novosibirsk Region, 630058, Russian Federation
| | - Petr N Menshanov
- Physiology Department, Novosibirsk State University, Novosibirsk, Novosibirsk Region, 630090, Russian Federation. .,Laser Systems Department, Novosibirsk State Technical University, Novosibirsk, Novosibirsk Region, 630073, Russian Federation. .,AI Tech Department, Novosibirsk State University, Novosibirsk, Novosibirsk Region, 630090, Russian Federation.
| |
Collapse
|
11
|
Expression of type 1 cannabinoid receptor gene in bipolar disorder. J Psychiatr Res 2022; 156:406-413. [PMID: 36323143 DOI: 10.1016/j.jpsychires.2022.10.006] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/23/2021] [Revised: 08/22/2022] [Accepted: 10/03/2022] [Indexed: 11/21/2022]
Abstract
BACKGROUND The Endocannabinoid System (ECBs) may have a crucial role in bipolar disorder (BD). Previous reports have not detected abnormalities in the expression of the cannabinoid receptor gene CNR1, encoding for CB1. However, we hypothesized that differentiating between mania and depression may uncover differences in CNR1 expression levels. METHODS We recruited 44 subjects with BD type I (BD-I), in mania (n = 22) and depression (n = 22) and 25 Healthy Controls (HC). CNR1 gene expression was analyzed using a quantitative real-time polymerase chain reaction from peripheral blood mononuclear cells. Data were analyzed using frequentist non-parametric and Bayesian approaches (generalized location-scale model based on lognormal and gamma distributions). RESULTS Using the frequentist non-parametric approach, the depression group had lower CNR1 expression compared to the mania group (p = 0.004). In addition, there was a negative correlation between CNR1 expression and Hamilton Depression Scale score (rho = -0.37; p = 0.007). Bayesian analyses further revealed that CNR1 expression in the mania group was higher and less variable than among HC (>95% probability), while CNR1 expression in the depression group was lower and more variable than among HC (100% probability). LIMITATIONS Lack of participants with bipolar disorder in the euthymic phase, lack of toxicology screening and evaluation of CNR1 variants. CONCLUSION CNR1 expression is higher and less variable in mania than in depression. It is highly probable that these differences also distinguish individuals in different illness phases from healthy controls. Future studies are needed to clarify the role of the endocannabinoid system in bipolar disorder.
Collapse
|
12
|
Huang ST, Lederer J. DeepMoM: Robust Deep Learning With Median-of-Means. J Comput Graph Stat 2022. [DOI: 10.1080/10618600.2022.2090947] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
13
|
Graph-Based Integration of Histone Modification Profiles. MATHEMATICS 2022. [DOI: 10.3390/math10111842] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
In this work, we introduce a similarity-network-based approach to explore the role of interacting single-cell histone modification signals in haematopoiesis—the process of differentiation of blood cells. Histones are proteins that provide structural support to chromosomes. They are subject to chemical modifications—acetylation or methylation—that affect the degree of accessibility of genes and, in turn, the formation of different phenotypes. The concentration of histone modifications can be modelled as a continuous signal, which can be used to build single-cell profiles. In the present work, the profiles of cell types involved in haematopoiesis are built based on all the major histone modifications (i.e., H3K27ac, H3K27me3, H3K36me3, H3K4me1, H3K4me3, H3K9me3) by counting the number of peaks in the modification signals; then, the profiles are used to compute modification-specific similarity networks among the considered phenotypes. As histone modifications come as interacting signals, we applied a similarity network fusion technique to integrate these networks in a unique graph, with the aim of studying the simultaneous effect of all the modifications for the determination of different phenotypes. The networks permit defining of a graph-cut-based separation score for evaluating the homogeneity of subgroups of cell types corresponding to the myeloid and lymphoid phenotypes in the classical representation of the haematopoietic tree. Resulting scores show that separation into myeloid and lymphoid phenotypes reflects the actual process of haematopoiesis.
Collapse
|
14
|
Marchese S, Cancelmo L, Diab O, Cahn L, Aaronson C, Daskalakis NP, Schaffer J, Horn SR, Johnson JS, Schechter C, Desarnaud F, Bierer LM, Makotkine I, Flory JD, Crane M, Moline JM, Udasin IG, Harrison DJ, Roussos P, Charney DS, Koenen KC, Southwick SM, Yehuda R, Pietrzak RH, Huckins LM, Feder A. Altered gene expression and PTSD symptom dimensions in World Trade Center responders. Mol Psychiatry 2022; 27:2225-2246. [PMID: 35177824 DOI: 10.1038/s41380-022-01457-2] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/11/2021] [Revised: 11/16/2021] [Accepted: 01/18/2022] [Indexed: 11/09/2022]
Abstract
Despite experiencing a significant trauma, only a subset of World Trade Center (WTC) rescue and recovery workers developed posttraumatic stress disorder (PTSD). Identification of biomarkers is critical to the development of targeted interventions for treating disaster responders and potentially preventing the development of PTSD in this population. Analysis of gene expression from these individuals can help in identifying biomarkers of PTSD. We established a well-phenotyped sample of 371 WTC responders, recruited from a longitudinal WTC responder cohort using stratified random sampling, by obtaining blood, self-reported and clinical interview data. Using bulk RNA-sequencing from whole blood, we examined the association between gene expression and WTC-related PTSD symptom severity on (i) highest lifetime Clinician-Administered PTSD Scale (CAPS) score, (ii) past-month CAPS score, and (iii) PTSD symptom dimensions using a 5-factor model of re-experiencing, avoidance, emotional numbing, dysphoric arousal and anxious arousal symptoms. We corrected for sex, age, genotype-derived principal components and surrogate variables. Finally, we performed a meta-analysis with existing PTSD studies (total N = 1016), using case/control status as the predictor and correcting for these variables. We identified 66 genes significantly associated with total highest lifetime CAPS score (FDR-corrected p < 0.05), and 31 genes associated with total past-month CAPS score. Our more granular analyses of PTSD symptom dimensions identified additional genes that did not reach statistical significance in our analyses with total CAPS scores. In particular, we identified 82 genes significantly associated with lifetime anxious arousal symptoms. Several genes significantly associated with multiple PTSD symptom dimensions and total lifetime CAPS score (SERPINA1, RPS6KA1, and STAT3) have been previously associated with PTSD. Geneset enrichment of these findings has identified pathways significant in metabolism, immune signaling, other psychiatric disorders, neurological signaling, and cellular structure. Our meta-analysis revealed 10 genes that reached genome-wide significance, all of which were downregulated in cases compared to controls (CIRBP, TMSB10, FCGRT, CLIC1, RPS6KB2, HNRNPUL1, ALDOA, NACA, ZNF429 and COPE). Additionally, cellular deconvolution highlighted an enrichment in CD4 T cells and eosinophils in responders with PTSD compared to controls. The distinction in significant genes between total lifetime CAPS score and the anxious arousal symptom dimension of PTSD highlights a potential biological difference in the mechanism underlying the heterogeneity of the PTSD phenotype. Future studies should be clear about methods used to analyze PTSD status, as phenotypes based on PTSD symptom dimensions may yield different gene sets than combined CAPS score analysis. Potential biomarkers implicated from our meta-analysis may help improve therapeutic target development for PTSD.
Collapse
Affiliation(s)
- Shelby Marchese
- Pamela Sklar Division of Psychiatric Genomics, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Leo Cancelmo
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Olivia Diab
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Leah Cahn
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Cindy Aaronson
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Nikolaos P Daskalakis
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.,Department of Psychiatry, McLean Hospital, Harvard Medical School, Belmont, MA, USA
| | - Jamie Schaffer
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Sarah R Horn
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Jessica S Johnson
- Pamela Sklar Division of Psychiatric Genomics, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Clyde Schechter
- Department of Family and Social Medicine, Albert Einstein College of Medicine, Bronx, NY, USA
| | - Frank Desarnaud
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.,Department of Psychiatry, James J. Peters Veterans Affairs Medical Center, Bronx, NY, USA
| | - Linda M Bierer
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.,Department of Psychiatry, James J. Peters Veterans Affairs Medical Center, Bronx, NY, USA
| | - Iouri Makotkine
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.,Department of Psychiatry, James J. Peters Veterans Affairs Medical Center, Bronx, NY, USA
| | - Janine D Flory
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.,Department of Psychiatry, James J. Peters Veterans Affairs Medical Center, Bronx, NY, USA
| | - Michael Crane
- Department of Environmental Medicine and Public Health, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Jacqueline M Moline
- Department of Occupational Medicine, Epidemiology and Prevention, Zucker School of Medicine at Hofstra/Northwell, Great Neck, NY, USA
| | - Iris G Udasin
- Environmental and Occupational Health Sciences Institute, School of Public Health, Rutgers University, Piscataway, NJ, USA
| | - Denise J Harrison
- Department of Medicine, Division of Pulmonary Critical Care and Sleep Medicine, NYU School of Medicine, New York, NY, USA
| | - Panos Roussos
- Pamela Sklar Division of Psychiatric Genomics, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.,Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.,Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.,Icahn Institute for Genomics and Multiscale Biology, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.,Mental Illness Research, Education and Clinical Centers, James J. Peters Department of Veterans Affairs Medical Center, Bronx, NY, 14068, USA
| | - Dennis S Charney
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.,Department of Neuroscience, Icahn School of Medicine at Mount Sinai, New York, NY, USA.,Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Karestan C Koenen
- Massachusetts General Hospital, Psychiatric and Neurodevelopmental Genetics Unit (PNGU), Boston, MA, USA.,Broad Institute of MIT and Harvard, Stanley Center for Psychiatric Research, Cambridge, MA, USA.,Harvard School of Public Health, Department of Epidemiology, Boston, MA, USA
| | - Steven M Southwick
- Department of Psychiatry, Yale University School of Medicine, New Haven, CT, USA.,Department of Veterans Affairs National Center for Posttraumatic Stress Disorder, VA Connecticut Healthcare System, West Haven, CT, USA
| | - Rachel Yehuda
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.,Department of Psychiatry, James J. Peters Veterans Affairs Medical Center, Bronx, NY, USA
| | - Robert H Pietrzak
- Department of Psychiatry, Yale University School of Medicine, New Haven, CT, USA.,Department of Veterans Affairs National Center for Posttraumatic Stress Disorder, VA Connecticut Healthcare System, West Haven, CT, USA
| | - Laura M Huckins
- Pamela Sklar Division of Psychiatric Genomics, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA. .,Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA. .,Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA. .,Icahn Institute for Genomics and Multiscale Biology, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA. .,Mental Illness Research, Education and Clinical Centers, James J. Peters Department of Veterans Affairs Medical Center, Bronx, NY, 14068, USA. .,Seaver Autism Center for Research and Treatment, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.
| | - Adriana Feder
- Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA.
| |
Collapse
|
15
|
Qayed A, Han D. Homogeneity test of several high-dimensional covariance matrices for stationary processes under non-normality. COMMUN STAT-THEOR M 2021. [DOI: 10.1080/03610926.2021.1960375] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Abdullah Qayed
- School of Mathematical Sciences, Department of Statistics, Shanghai Jiao Tong University, Shanghai, China
| | - Dong Han
- School of Mathematical Sciences, Department of Statistics, Shanghai Jiao Tong University, Shanghai, China
| |
Collapse
|
16
|
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
|
17
|
Wu M, Xu J, Ding T, Gao J. Mixed Distribution Models Based on Single-Cell RNA Sequencing Data. Interdiscip Sci 2021; 13:362-370. [PMID: 33754241 DOI: 10.1007/s12539-021-00427-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2020] [Revised: 03/02/2021] [Accepted: 03/13/2021] [Indexed: 10/21/2022]
Abstract
Progress in single-cell RNA sequencing (scRNA-seq) has yielded a lot of valuable data. Analysis of these data can provide a new perspective for studying the intratumoral heterogeneity and identifying gene markers. In this paper, the scRNA-seq data of colorectal cancer (CRC) are analyzed, and it is found that the shape of the gene expression difference (GED) data shows certain distribution regularity. To study the distribution regularity, mixed stable-normal distribution (MSND) model and mixed stable-exponential distribution (MSED) model are constructed to fit the GED data. And the estimated parameters of MSND and MSED are used to describe some characteristics of their distribution. Through the comparison of root mean square error and the chi-squared goodness of fit test, it is found that the fitting effect of MSED and MSND are both better than that of stable distribution and Cauchy distribution. Considering the given quantile thresholds, MSND and MSED can be used to identify tumor-related genes. The results of functional analysis indicate that the selected genes are highly correlated with CRC. In addition, the parameters of MSND and MSED exhibit a certain trend with the development of CRC. To explore the association, Gene-set enrichment analysis (GSEA) is performed. The results of GSEA reveal that the trend can well characterize the intratumoral heterogeneity of CRC. In addition, the application of MSED model on hepatocellular carcinoma shows that our model can analyze other cancers. Overall, MSND model and MSED model can well fit the GED data in different disease stages, the parameters of the two models can characterize the heterogeneity of CRC tumor cells, and the two models can be used to identify genes highly correlated with tumors.
Collapse
Affiliation(s)
- Min Wu
- School of Science, Jiangnan University, Wuxi, 214122, China
| | - Junhua Xu
- School of Science, Jiangnan University, Wuxi, 214122, China
| | - Tao Ding
- School of Mathematics Statistics and Physics, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK
| | - Jie Gao
- School of Science, Jiangnan University, Wuxi, 214122, China.
| |
Collapse
|
18
|
Guo Y, Shen L, Shi X, Wang K, Dai Y, Zhao Z. Accelerating bioinformatics research with International Conference on Intelligent Biology and Medicine 2020. BMC Bioinformatics 2020; 21:563. [PMID: 33371868 PMCID: PMC7767910 DOI: 10.1186/s12859-020-03890-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/18/2020] [Indexed: 11/10/2022] Open
Abstract
The International Association for Intelligent Biology and Medicine (IAIBM) is a nonprofit organization that promotes intelligent biology and medical science. It hosts an annual International Conference on Intelligent Biology and Medicine (ICIBM), which was initially established in 2012. Due to the coronavirus (COVID-19) pandemic, the ICIBM 2020 was held for the first time as a virtual online conference on August 9 to 10. The virtual conference had ~ 300 registered participants and featured 41 online real-time presentations. ICIBM 2020 received a total of 75 manuscript submissions, and 12 were selected to be published in this special issue of BMC Bioinformatics. These 12 manuscripts cover a wide range of bioinformatics topics including network analysis, imaging analysis, machine learning, gene expression analysis, and sequence analysis.
Collapse
Affiliation(s)
- Yan Guo
- Department of Internal Medicine, Comprehensive Cancer Center, University of New Mexico, Albuquerque, NM 87131 USA
| | - Li Shen
- Department of Biostatistics, Epidemiology and Informatics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104 USA
| | - Xinghua Shi
- Department of Computer and Information Sciences, College of Science and Technology, Temple University, Philadelphia, PA 19122 USA
| | - Kai Wang
- Raymond G. Perelman Center for Cellular and Molecular Therapeutics, Children’s Hospital of Philadelphia, Philadelphia, PA 19104 USA
| | - Yulin Dai
- Center for Precision Health, School of Biomedical Informatics, The University of Texas Health Science Center at Houston, Houston, TX 77030 USA
| | - Zhongming Zhao
- Center for Precision Health, School of Biomedical Informatics, The University of Texas Health Science Center at Houston, Houston, TX 77030 USA
| |
Collapse
|