1
|
Wang Y, Wei Z, Su J, Coenen F, Meng J. RgnTX: Colocalization analysis of transcriptome elements in the presence of isoform heterogeneity and ambiguity. Comput Struct Biotechnol J 2023; 21:4110-4117. [PMID: 37671241 PMCID: PMC10475473 DOI: 10.1016/j.csbj.2023.08.021] [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: 03/14/2023] [Revised: 08/13/2023] [Accepted: 08/23/2023] [Indexed: 09/07/2023] Open
Abstract
Colocalization analysis of genomic region sets has been widely adopted to unveil potential functional interactions between corresponding biological attributes, which often serves as the basis for further investigation. A number of methods have been developed for colocalization analysis of genomic elements. However, none of them explicitly considered the transcriptome heterogeneity and isoform ambiguity, making them less appropriate for analyzing transcriptome elements. Here, we developed RgnTX, an R/Bioconductor tool for the colocalization analysis of transcriptome elements with permutation tests. Different from existing approaches, RgnTX directly takes advantage of transcriptome annotation, and offers high flexibility in the null model to simulate realistic transcriptome-wide background, such as the complex alternative splicing patterns. Importantly, it supports the testing of transcriptome elements without clear isoform association, which is often the real scenario due to technical limitations. Proposed package offers a wide selection of pre-defined functions, easy to be utilized by users for visualizing permutation results, calculating shifted z-scores and conducting multiple hypothesis testing under Benjamini-Hochberg correction. Moreover, with synthetic and real datasets, we show that RgnTX novel testing modes return distinct and more significant results compared to existing genome-based methods. We believe RgnTX should make a useful tool to characterize the randomness of the transcriptome, and for conducting statistical association analysis for genomic region sets within the heterogeneous transcriptome. The package now has been accepted by Bioconductor and is freely available at: https://bioconductor.org/packages/RgnTX.
Collapse
Affiliation(s)
- Yue Wang
- Department of Mathematical Sciences, Xi’an Jiaotong-Liverpool University, Suzhou, Jiangsu 215123, China
- Department of Computer Science, University of Liverpool, L69 7ZB Liverpool, United Kingdom
| | - Zhen Wei
- Department of Biological Sciences, Xi’an Jiaotong-Liverpool University, Suzhou, Jiangsu 215123, China
- Institute of Systems, Molecular and Integrative Biology, University of Liverpool, L69 7ZB Liverpool, United Kingdom
| | - Jionglong Su
- School of AI and Advanced Computing, Xi’an Jiaotong-Liverpool University, Suzhou, Jiangsu 215123, China
| | - Frans Coenen
- Department of Computer Science, University of Liverpool, L69 7ZB Liverpool, United Kingdom
| | - Jia Meng
- Department of Biological Sciences, Xi’an Jiaotong-Liverpool University, Suzhou, Jiangsu 215123, China
- AI University Research Centre, Xi’an Jiaotong-Liverpool University, Suzhou, Jiangsu 215123, China
- Institute of Systems, Molecular and Integrative Biology, University of Liverpool, L69 7ZB Liverpool, United Kingdom
| |
Collapse
|
2
|
Kuksa PP, Leung Y, Gangadharan P, Katanic Z, Kleidermacher L, Amlie-Wolf A, Lee CY, Qu L, Greenfest-Allen E, Valladares O, Wang LS. FILER: a framework for harmonizing and querying large-scale functional genomics knowledge. NAR Genom Bioinform 2022; 4:lqab123. [PMID: 35047815 PMCID: PMC8759563 DOI: 10.1093/nargab/lqab123] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2021] [Revised: 11/02/2021] [Accepted: 12/15/2021] [Indexed: 12/30/2022] Open
Abstract
Querying massive functional genomic and annotation data collections, linking and summarizing the query results across data sources/data types are important steps in high-throughput genomic and genetic analytical workflows. However, these steps are made difficult by the heterogeneity and breadth of data sources, experimental assays, biological conditions/tissues/cell types and file formats. FILER (FunctIonaL gEnomics Repository) is a framework for querying large-scale genomics knowledge with a large, curated integrated catalog of harmonized functional genomic and annotation data coupled with a scalable genomic search and querying interface. FILER uniquely provides: (i) streamlined access to >50 000 harmonized, annotated genomic datasets across >20 integrated data sources, >1100 tissues/cell types and >20 experimental assays; (ii) a scalable genomic querying interface; and (iii) ability to analyze and annotate user's experimental data. This rich resource spans >17 billion GRCh37/hg19 and GRCh38/hg38 genomic records. Our benchmark querying 7 × 109 hg19 FILER records shows FILER is highly scalable, with a sub-linear 32-fold increase in querying time when increasing the number of queries 1000-fold from 1000 to 1 000 000 intervals. Together, these features facilitate reproducible research and streamline integrating/querying large-scale genomic data within analyses/workflows. FILER can be deployed on cloud or local servers (https://bitbucket.org/wanglab-upenn/FILER) for integration with custom pipelines and is freely available (https://lisanwanglab.org/FILER).
Collapse
Affiliation(s)
| | - Yuk Yee Leung
- Correspondence may also be addressed to Yuk Yee Leung. Tel: +1 215 573 3729; Fax: +1 215 573 3111;
| | | | | | | | | | | | | | | | | | - Li-San Wang
- To whom correspondence should be addressed. Tel: +1 215 746 7015; Fax: +1 215 573 3111;
| |
Collapse
|
3
|
Van Remmen H, Freeman WM, Miller BF, Kinter M, Wren JD, Chiao A, Towner RA, Snider TA, Sonntag WE, Richardson A. Oklahoma Nathan Shock Aging Center - assessing the basic biology of aging from genetics to protein and function. GeroScience 2021; 43:2183-2203. [PMID: 34606039 PMCID: PMC8599778 DOI: 10.1007/s11357-021-00454-7] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Accepted: 09/03/2021] [Indexed: 01/21/2023] Open
Abstract
The Oklahoma Shock Nathan Shock Center is designed to deliver unique, innovative services that are not currently available at most institutions. The focus of the Center is on geroscience and the development of careers of young investigators. Pilot grants are provided through the Research Development Core to junior investigators studying aging/geroscience throughout the USA. However, the services of our Center are available to the entire research community studying aging and geroscience. The Oklahoma Nathan Shock Center provides researchers with unique services through four research cores. The Multiplexing Protein Analysis Core uses the latest mass spectrometry technology to simultaneously measure the levels, synthesis, and turnover of hundreds of proteins associated with pathways of importance to aging, e.g., metabolism, antioxidant defense system, proteostasis, and mitochondria function. The Genomic Sciences Core uses novel next-generation sequencing that allows investigators to study the effect of age, or anti-aging manipulations, on DNA methylation, mitochondrial genome heteroplasmy, and the transcriptome of single cells. The Geroscience Redox Biology Core provides investigators with a comprehensive state-of-the-art assessment of the oxidative stress status of a cell, e.g., measures of oxidative damage and redox couples, which are important in aging as well as many major age-related diseases as well as assays of mitochondrial function. The GeroInformatics Core provides investigators assistance with data analysis, which includes both statistical support as well as analysis of large datasets. The Core also has developed number of unique software packages to help with interpretation of results and discovery of new leads relevant to aging. In addition, the Geropathology Research Resource in the Program Enhancement Core provides investigators with pathological assessments of mice using the recently developed Geropathology Grading Platform.
Collapse
Affiliation(s)
- Holly Van Remmen
- Aging & Metabolism Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, 73104, USA.
- Biochemistry & Molecular Biology, The University of Oklahoma Health Sciences Center, Oklahoma, City, OK, USA.
| | - Willard M Freeman
- Biochemistry & Molecular Biology, The University of Oklahoma Health Sciences Center, Oklahoma, City, OK, USA
- Genes & Human Disease Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, USA
| | - Benjamin F Miller
- Aging & Metabolism Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, 73104, USA
- Biochemistry & Molecular Biology, The University of Oklahoma Health Sciences Center, Oklahoma, City, OK, USA
| | - Michael Kinter
- Aging & Metabolism Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, 73104, USA
| | - Jonathan D Wren
- Genes & Human Disease Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, USA
| | - Ann Chiao
- Aging & Metabolism Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, 73104, USA
| | - Rheal A Towner
- Advanced Magnetic Resonance Center, Oklahoma Medical Research Foundation, Oklahoma City, OK, USA
| | - Timothy A Snider
- Department of Veterinary Pathology, Center for Veterinary Health Sciences at, Oklahoma State University, Stillwater, OK, USA
| | - William E Sonntag
- Biochemistry & Molecular Biology, The University of Oklahoma Health Sciences Center, Oklahoma, City, OK, USA
| | - Arlan Richardson
- Biochemistry & Molecular Biology, The University of Oklahoma Health Sciences Center, Oklahoma, City, OK, USA
- Oklahoma City VA Medical Center, Oklahoma City, OK, USA
| |
Collapse
|
4
|
Gu A, Cho HJ, Sheffield NC. Bedshift: perturbation of genomic interval sets. Genome Biol 2021; 22:238. [PMID: 34416909 PMCID: PMC8379854 DOI: 10.1186/s13059-021-02440-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2020] [Accepted: 07/26/2021] [Indexed: 12/25/2022] Open
Abstract
Functional genomics experiments, like ChIP-Seq or ATAC-Seq, produce results that are summarized as a region set. There is no way to objectively evaluate the effectiveness of region set similarity metrics. We present Bedshift, a tool for perturbing BED files by randomly shifting, adding, and dropping regions from a reference file. The perturbed files can be used to benchmark similarity metrics, as well as for other applications. We highlight differences in behavior between metrics, such as that the Jaccard score is most sensitive to added or dropped regions, while coverage score is most sensitive to shifted regions.
Collapse
Affiliation(s)
- Aaron Gu
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA, USA
- Department of Computer Science, University of Virginia School of Engineering, Charlottesville, VA, USA
| | - Hyun Jae Cho
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA, USA
- Department of Computer Science, University of Virginia School of Engineering, Charlottesville, VA, USA
| | - Nathan C Sheffield
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA, USA.
- Department of Public Health Sciences, University of Virginia, Charlottesville, VA, USA.
- Department of Biomedical Engineering, University of Virginia, Charlottesville, VA, USA.
- Department of Biochemistry and Molecular Genetics, University of Virginia, Charlottesville, VA, USA.
| |
Collapse
|
5
|
Fuemmeler BF, Dozmorov MG, Do EK, Zhang J(J, Grenier C, Huang Z, Maguire RL, Kollins SH, Hoyo C, Murphy SK. DNA Methylation in Babies Born to Nonsmoking Mothers Exposed to Secondhand Smoke during Pregnancy: An Epigenome-Wide Association Study. ENVIRONMENTAL HEALTH PERSPECTIVES 2021; 129:57010. [PMID: 34009014 PMCID: PMC8132610 DOI: 10.1289/ehp8099] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/17/2020] [Revised: 02/09/2021] [Accepted: 04/19/2021] [Indexed: 06/12/2023]
Abstract
BACKGROUND Maternal smoking during pregnancy is related to altered DNA methylation in infant umbilical cord blood. The extent to which low levels of smoke exposure among nonsmoking pregnant women relates to offspring DNA methylation is unknown. OBJECTIVE This study sought to evaluate relationships between maternal prenatal plasma cotinine levels and DNA methylation in umbilical cord blood in newborns using the Infinium HumanMethylation 450K BeadChip. METHODS Participants from the Newborn Epigenetics Study cohort who reported not smoking during pregnancy had verified low levels of cotinine from maternal prenatal plasma (0 ng / mL to < 4 ng / mL ), and offspring epigenetic data from umbilical cord blood were included in this study (n = 79 ). Multivariable linear regression models were fit to the data, controlling for cell proportions, age, race, education, and parity. Estimates represent changes in response to any 1 -ng / mL unit increase in exposure. RESULTS Multivariable linear regression models yielded 29,049 CpGs that were differentially methylated in relation to increases in cotinine at a 5% false discovery rate. Top CpGs were within or near genes involved in neuronal functioning (PRKG1, DLGAP2, BSG), carcinogenesis (FHIT, HSPC157) and inflammation (AGER). Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggest cotinine was related to methylation of gene pathways controlling neuronal signaling, metabolic regulation, cell signaling and regulation, and cancer. Further, enhancers associated with transcription start sites were enriched in altered CpGs. Using an independent sample from the same study population (n = 115 ), bisulfite pyrosequencing was performed with infant cord blood DNA for two genes within our top 20 hits (AGER and PRKG1). Results from pyrosequencing replicated epigenome results for PRKG1 (cg17079497, estimate = - 1.09 , standard error ( SE ) = 0.45 , p = 0.018 ) but not for AGER (cg09199225; estimate = - 0.16 , SE = 0.21 , p = 0.44 ). DISCUSSION Secondhand smoke exposure among nonsmoking women may alter DNA methylation in regions involved in development, carcinogenesis, and neuronal functioning. These novel findings suggest that even low levels of smoke exposure during pregnancy may be sufficient to alter DNA methylation in distinct sites of mixed umbilical cord blood leukocytes in pathways that are known to be altered in cord blood from pregnant active smokers. https://doi.org/10.1289/EHP8099.
Collapse
Affiliation(s)
- Bernard F. Fuemmeler
- Department of Health Behavior and Policy, Virginia Commonwealth University, Richmond, Virginia, USA
- Massey Cancer Center, Virginia Commonwealth University, Richmond, Virginia, USA
| | - Mikhail G. Dozmorov
- Department of Biostatistics, Virginia Commonwealth University, Richmond, Virginia, USA
- Department of Pathology, Virginia Commonwealth University, Richmond, Virginia, USA
| | - Elizabeth K. Do
- Department of Health Behavior and Policy, Virginia Commonwealth University, Richmond, Virginia, USA
- Massey Cancer Center, Virginia Commonwealth University, Richmond, Virginia, USA
| | - Junfeng (Jim) Zhang
- Nicholas School of the Environment and Duke Global Health Institute, Duke University, Durham, North Carolina, USA
| | - Carole Grenier
- Department of Obstetrics and Gynecology, Duke University Medical Center, Durham, North Carolina, USA
| | - Zhiqing Huang
- Department of Obstetrics and Gynecology, Duke University Medical Center, Durham, North Carolina, USA
| | - Rachel L. Maguire
- Department of Obstetrics and Gynecology, Duke University Medical Center, Durham, North Carolina, USA
- Department of Biological Sciences, Center for Human Health and the Environment North Carolina State University, Raleigh, North Carolina, USA
| | - Scott H. Kollins
- Department of Psychiatry and Behavioral Sciences, Duke University, Durham, North Carolina, USA
| | - Cathrine Hoyo
- Department of Biological Sciences, Center for Human Health and the Environment North Carolina State University, Raleigh, North Carolina, USA
| | - Susan K. Murphy
- Department of Obstetrics and Gynecology, Duke University Medical Center, Durham, North Carolina, USA
| |
Collapse
|
6
|
Dozmorov MG, Tyc KM, Sheffield NC, Boyd DC, Olex AL, Reed J, Harrell JC. Chromatin conformation capture (Hi-C) sequencing of patient-derived xenografts: analysis guidelines. Gigascience 2021; 10:6242846. [PMID: 33880552 DOI: 10.1093/gigascience/giab022] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2020] [Revised: 01/14/2021] [Accepted: 03/09/2021] [Indexed: 12/14/2022] Open
Abstract
BACKGROUND Sequencing of patient-derived xenograft (PDX) mouse models allows investigation of the molecular mechanisms of human tumor samples engrafted in a mouse host. Thus, both human and mouse genetic material is sequenced. Several methods have been developed to remove mouse sequencing reads from RNA-seq or exome sequencing PDX data and improve the downstream signal. However, for more recent chromatin conformation capture technologies (Hi-C), the effect of mouse reads remains undefined. RESULTS We evaluated the effect of mouse read removal on the quality of Hi-C data using in silico created PDX Hi-C data with 10% and 30% mouse reads. Additionally, we generated 2 experimental PDX Hi-C datasets using different library preparation strategies. We evaluated 3 alignment strategies (Direct, Xenome, Combined) and 3 pipelines (Juicer, HiC-Pro, HiCExplorer) on Hi-C data quality. CONCLUSIONS Removal of mouse reads had little-to-no effect on data quality as compared with the results obtained with the Direct alignment strategy. Juicer extracted more valid chromatin interactions for Hi-C matrices, regardless of the mouse read removal strategy. However, the pipeline effect was minimal, while the library preparation strategy had the largest effect on all quality metrics. Together, our study presents comprehensive guidelines on PDX Hi-C data processing.
Collapse
Affiliation(s)
- Mikhail G Dozmorov
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA 23298, USA.,Department of Pathology, Virginia Commonwealth University, Richmond, VA 23284, USA
| | - Katarzyna M Tyc
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA 23298, USA.,Department of Pharmacology and Toxicology, Virginia Commonwealth University, Richmond, VA, 23298, USA
| | - Nathan C Sheffield
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA
| | - David C Boyd
- Department of Pathology, Virginia Commonwealth University, Richmond, VA 23284, USA.,Integrative Life Sciences Doctoral Program, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Amy L Olex
- C. Kenneth and Dianne Wright Center for Clinical and Translational Research, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Jason Reed
- Virginia Commonwealth University, Massey Cancer Center, Richmond, VA, 23298, USA.,Department of Physics, Virginia Commonwealth University, Richmond, VA 23220, USA
| | - J Chuck Harrell
- Department of Pathology, Virginia Commonwealth University, Richmond, VA 23284, USA
| |
Collapse
|
7
|
Gundersen S, Boddu S, Capella-Gutierrez S, Drabløs F, Fernández JM, Kompova R, Taylor K, Titov D, Zerbino D, Hovig E. Recommendations for the FAIRification of genomic track metadata. F1000Res 2021; 10. [PMID: 34249331 PMCID: PMC8226415 DOI: 10.12688/f1000research.28449.1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 03/17/2021] [Indexed: 01/25/2023] Open
Abstract
Background: Many types of data from genomic analyses can be represented as genomic tracks,
i.e. features linked to the genomic coordinates of a reference genome. Examples of such data are epigenetic DNA methylation data, ChIP-seq peaks, germline or somatic DNA variants, as well as RNA-seq expression levels. Researchers often face difficulties in locating, accessing and combining relevant tracks from external sources, as well as locating the raw data, reducing the value of the generated information. Description of work: We propose to advance the application of FAIR data principles (Findable, Accessible, Interoperable, and Reusable) to produce searchable metadata for genomic tracks. Findability and Accessibility of metadata can then be ensured by a track search service that integrates globally identifiable metadata from various track hubs in the Track Hub Registry and other relevant repositories. Interoperability and Reusability need to be ensured by the specification and implementation of a basic set of recommendations for metadata. We have tested this concept by developing such a specification in a JSON Schema, called FAIRtracks, and have integrated it into a novel track search service, called TrackFind. We demonstrate practical usage by importing datasets through TrackFind into existing examples of relevant analytical tools for genomic tracks: EPICO and the GSuite HyperBrowser. Conclusion: We here provide a first iteration of a draft standard for genomic track metadata, as well as the accompanying software ecosystem. It can easily be adapted or extended to future needs of the research community regarding data, methods and tools, balancing the requirements of both data submitters and analytical end-users.
Collapse
Affiliation(s)
| | - Sanjay Boddu
- European Molecular Biology Laboratory, European Bioinformatics Institute, Hinxton, UK
| | | | - Finn Drabløs
- Department of Clinical and Molecular Medicine, NTNU - Norwegian University of Science and Technology, Trondheim, Norway
| | - José M Fernández
- Life Sciences Department, Barcelona Supercomputing Center (BSC), Barcelona, Spain
| | - Radmila Kompova
- Center for Bioinformatics, University of Oslo (UiO), Oslo, Norway
| | - Kieron Taylor
- European Molecular Biology Laboratory, European Bioinformatics Institute, Hinxton, UK
| | - Dmytro Titov
- Center for Bioinformatics, University of Oslo (UiO), Oslo, Norway
| | - Daniel Zerbino
- European Molecular Biology Laboratory, European Bioinformatics Institute, Hinxton, UK
| | - Eivind Hovig
- Center for Bioinformatics, University of Oslo (UiO), Oslo, Norway.,Department of Tumor Biology, Institute for Cancer Research, Oslo University Hospital (OUH), Oslo, Norway
| |
Collapse
|
8
|
Ortiz-Fernández L, Sawalha AH. Genetics of Behçet's Disease: Functional Genetic Analysis and Estimating Disease Heritability. Front Med (Lausanne) 2021; 8:625710. [PMID: 33644100 PMCID: PMC7907152 DOI: 10.3389/fmed.2021.625710] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2020] [Accepted: 01/14/2021] [Indexed: 01/17/2023] Open
Abstract
Behçet's disease is a chronic multisystemic inflammatory disorder characterized by recurrent oral and genital ulcers. Although its etiology remains unclear, it is thought that both genetic and environmental factors contribute to the onset and progression of Behçet's disease. Here, we provide an updated view of the genetic landscape and architecture of Behçet's disease. Large-scale genetic studies performed to date revealed 21 genetic susceptibility loci associated with the disease at a GWAS level of significance (p-value = 5 × 10−8). We performed epigenetic pattern enrichment analysis in Behçet's disease associated loci, providing new insights into the molecular mechanisms underlying its pathophysiology. Our data suggest the crucial involvement of several immune cell types, including natural killer cells, monocytes, and B cells in the pathogenesis of the disease. Pathway enrichment analysis identified important biological processes involved. Using large-scale genetic data available from ~200 immune-related loci (Immunochip), we estimate Behçet's disease heritability to be at least 16%. We further used the same approach to estimate the heritability explained by the known Behçet's disease-associated loci, suggesting that they explain ~ 60% of the genetic component underlying Behçet's disease. These results indicate a significant role of non-genetic factors in causing Behçet's disease and that additional genetic variation influencing the risk of Behçet's disease remains to be identified. Finally, we calculated a cumulative genetic risk score across populations reinforcing the link between geographic variations in disease prevalence with its genetic component.
Collapse
Affiliation(s)
- Lourdes Ortiz-Fernández
- Division of Rheumatology, Department of Pediatrics, University of Pittsburgh School of Medicine, Pittsburgh, PA, United States
| | - Amr H Sawalha
- Division of Rheumatology, Department of Pediatrics, University of Pittsburgh School of Medicine, Pittsburgh, PA, United States.,Division of Rheumatology and Clinical Immunology, Department of Medicine, University of Pittsburgh School of Medicine, Pittsburgh, PA, United States.,Lupus Center of Excellence, University of Pittsburgh School of Medicine, Pittsburgh, PA, United States.,Department of Immunology, University of Pittsburgh School of Medicine, Pittsburgh, PA, United States
| |
Collapse
|
9
|
Ortiz-Fernández L, Saruhan-Direskeneli G, Alibaz-Oner F, Kaymaz-Tahra S, Coit P, Kong X, Kiprianos AP, Maughan RT, Aydin SZ, Aksu K, Keser G, Kamali S, Inanc M, Springer J, Akar S, Onen F, Akkoc N, Khalidi NA, Koening C, Karadag O, Kiraz S, Forbess L, Langford CA, McAlear CA, Ozbalkan Z, Yavuz S, Çetin GY, Alpay-Kanitez N, Chung S, Ates A, Karaaslan Y, McKinnon-Maksimowicz K, Monach PA, Ozer HT, Seyahi E, Fresko I, Cefle A, Seo P, Warrington KJ, Ozturk MA, Ytterberg SR, Cobankara V, Onat AM, Duzgun N, Bıcakcıgil M, Yentür SP, Lally L, Manfredi AA, Baldissera E, Erken E, Yazici A, Kısacık B, Kaşifoğlu T, Dalkilic E, Cuthbertson D, Pagnoux C, Sreih A, Reales G, Wallace C, Wren JD, Cunninghame-Graham DS, Vyse TJ, Sun Y, Chen H, Grayson PC, Tombetti E, Jiang L, Mason JC, Merkel PA, Direskeneli H, Sawalha AH. Identification of susceptibility loci for Takayasu arteritis through a large multi-ancestral genome-wide association study. Am J Hum Genet 2021; 108:84-99. [PMID: 33308445 DOI: 10.1016/j.ajhg.2020.11.014] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2020] [Accepted: 11/17/2020] [Indexed: 12/21/2022] Open
Abstract
Takayasu arteritis is a rare inflammatory disease of large arteries. We performed a genetic study in Takayasu arteritis comprising 6,670 individuals (1,226 affected individuals) from five different populations. We discovered HLA risk factors and four non-HLA susceptibility loci in VPS8, SVEP1, CFL2, and chr13q21 and reinforced IL12B, PTK2B, and chr21q22 as robust susceptibility loci shared across ancestries. Functional analysis proposed plausible underlying disease mechanisms and pinpointed ETS2 as a potential causal gene for chr21q22 association. We also identified >60 candidate loci with suggestive association (p < 5 × 10-5) and devised a genetic risk score for Takayasu arteritis. Takayasu arteritis was compared to hundreds of other traits, revealing the closest genetic relatedness to inflammatory bowel disease. Epigenetic patterns within risk loci suggest roles for monocytes and B cells in Takayasu arteritis. This work enhances understanding of the genetic basis and pathophysiology of Takayasu arteritis and provides clues for potential new therapeutic targets.
Collapse
|
10
|
Cresswell KG, Stansfield JC, Dozmorov MG. SpectralTAD: an R package for defining a hierarchy of topologically associated domains using spectral clustering. BMC Bioinformatics 2020; 21:319. [PMID: 32689928 PMCID: PMC7372752 DOI: 10.1186/s12859-020-03652-w] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2019] [Accepted: 07/10/2020] [Indexed: 01/08/2023] Open
Abstract
BACKGROUND The three-dimensional (3D) structure of the genome plays a crucial role in gene expression regulation. Chromatin conformation capture technologies (Hi-C) have revealed that the genome is organized in a hierarchy of topologically associated domains (TADs), sub-TADs, and chromatin loops. Identifying such hierarchical structures is a critical step in understanding genome regulation. Existing tools for TAD calling are frequently sensitive to biases in Hi-C data, depend on tunable parameters, and are computationally inefficient. METHODS To address these challenges, we developed a novel sliding window-based spectral clustering framework that uses gaps between consecutive eigenvectors for TAD boundary identification. RESULTS Our method, implemented in an R package, SpectralTAD, detects hierarchical, biologically relevant TADs, has automatic parameter selection, is robust to sequencing depth, resolution, and sparsity of Hi-C data. SpectralTAD outperforms four state-of-the-art TAD callers in simulated and experimental settings. We demonstrate that TAD boundaries shared among multiple levels of the TAD hierarchy were more enriched in classical boundary marks and more conserved across cell lines and tissues. In contrast, boundaries of TADs that cannot be split into sub-TADs showed less enrichment and conservation, suggesting their more dynamic role in genome regulation. CONCLUSION SpectralTAD is available on Bioconductor, http://bioconductor.org/packages/SpectralTAD/ .
Collapse
Affiliation(s)
- Kellen G Cresswell
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA
| | - John C Stansfield
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA
| | - Mikhail G Dozmorov
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA.
| |
Collapse
|
11
|
Smith JP, Sheffield NC. Analytical Approaches for ATAC-seq Data Analysis. CURRENT PROTOCOLS IN HUMAN GENETICS 2020; 106:e101. [PMID: 32543102 PMCID: PMC8191135 DOI: 10.1002/cphg.101] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
ATAC-seq, the assay for transposase-accessible chromatin using sequencing, is a quick and efficient approach to investigating the chromatin accessibility landscape. Investigating chromatin accessibility has broad utility for answering many biological questions, such as mapping nucleosomes, identifying transcription factor binding sites, and measuring differential activity of DNA regulatory elements. Because the ATAC-seq protocol is both simple and relatively inexpensive, there has been a rapid increase in the availability of chromatin accessibility data. Furthermore, advances in ATAC-seq protocols are rapidly extending its breadth to additional experimental conditions, cell types, and species. Accompanying the increase in data, there has also been an explosion of new tools and analytical approaches for analyzing it. Here, we explain the fundamentals of ATAC-seq data processing, summarize common analysis approaches, and review computational tools to provide recommendations for different research questions. This primer provides a starting point and a reference for analysis of ATAC-seq data. © 2020 Wiley Periodicals LLC.
Collapse
Affiliation(s)
- Jason P. Smith
- Center for Public Health Genomics, University of Virginia, Charlottesville, Virginia
- Department of Biochemistry and Molecular Genetics, University of Virginia, Charlottesville, Virginia
| | - Nathan C. Sheffield
- Center for Public Health Genomics, University of Virginia, Charlottesville, Virginia
- Department of Biochemistry and Molecular Genetics, University of Virginia, Charlottesville, Virginia
- Department of Public Health Sciences, University of Virginia, Charlottesville, Virginia
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia
| |
Collapse
|
12
|
Cresswell KG, Dozmorov MG. TADCompare: An R Package for Differential and Temporal Analysis of Topologically Associated Domains. Front Genet 2020; 11:158. [PMID: 32211023 PMCID: PMC7076128 DOI: 10.3389/fgene.2020.00158] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2019] [Accepted: 02/11/2020] [Indexed: 12/02/2022] Open
Abstract
Recent research using chromatin conformation capture technologies, such as Hi-C, has demonstrated the importance of topologically associated domains (TADs) and smaller chromatin loops, collectively referred hereafter as "interacting domains." Many such domains change during development or disease, and exhibit cell- and condition-specific differences. Quantification of the dynamic behavior of interacting domains will help to better understand genome regulation. Methods for comparing interacting domains between cells and conditions are highly limited. We developed TADCompare, a method for differential analysis of boundaries of interacting domains between two or more Hi-C datasets. TADCompare is based on a spectral clustering-derived measure called the eigenvector gap, which enables a loci-by-loci comparison of boundary differences. Using this measure, we introduce methods for identifying differential and consensus boundaries of interacting domains and tracking boundary changes over time. We further propose a novel framework for the systematic classification of boundary changes. Colocalization- and gene enrichment analysis of different types of boundary changes demonstrated distinct biological functionality associated with them. TADCompare is available on https://github.com/dozmorovlab/TADCompare and Bioconductor (submitted).
Collapse
|
13
|
Zhou Y, Sun Y, Huang D, Li MJ. epiCOLOC: Integrating Large-Scale and Context-Dependent Epigenomics Features for Comprehensive Colocalization Analysis. Front Genet 2020; 11:53. [PMID: 32117461 PMCID: PMC7029718 DOI: 10.3389/fgene.2020.00053] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2019] [Accepted: 01/17/2020] [Indexed: 12/18/2022] Open
Abstract
High-throughput genome-wide epigenomic assays, such as ChIP-seq, DNase-seq and ATAC-seq, have profiled a huge number of functional elements across numerous human tissues/cell types, which provide an unprecedented opportunity to interpret human genome and disease in context-dependent manner. Colocalization analysis determines whether genomic features are functionally related to a given search and will facilitate identifying the underlying biological functions characterizing intricate relationships with queries for genomic regions. Existing colocalization methods leveraged diverse assumptions and background models to assess the significance of enrichment, however, they only provided limited and predefined sets of epigenomic features. Here, we comprehensively collected and integrated over 44,385 bulk or single-cell epigenomic assays across 53 human tissues/cell types, such as transcription factor binding, histone modification, open chromatin and transcriptional event. By classifying these profiles into hierarchy of tissue/cell type, we developed a web portal, epiCOLOC (http://mulinlab.org/epicoloc or http://mulinlab.tmu.edu.cn/epicoloc), for users to perform context-dependent colocalization analysis in a convenient way.
Collapse
Affiliation(s)
- Yao Zhou
- Department of Pharmacology, Tianjin Key Laboratory of Inflammation Biology, School of Basic Medical Sciences, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, Tianjin Medical University, Tianjin, China
| | - Yongzheng Sun
- Department of Pharmacology, Tianjin Key Laboratory of Inflammation Biology, School of Basic Medical Sciences, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, Tianjin Medical University, Tianjin, China
| | - Dandan Huang
- Department of Pharmacology, Tianjin Key Laboratory of Inflammation Biology, School of Basic Medical Sciences, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, Tianjin Medical University, Tianjin, China
| | - Mulin Jun Li
- Department of Pharmacology, Tianjin Key Laboratory of Inflammation Biology, School of Basic Medical Sciences, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, Tianjin Medical University, Tianjin, China.,Collaborative Innovation Center of Tianjin for Medical Epigenetics, Tianjin Key Laboratory of Medical Epigenetics, Tianjin Medical University, Tianjin, China
| |
Collapse
|
14
|
Nagraj VP, Magee NE, Sheffield NC. LOLAweb: a containerized web server for interactive genomic locus overlap enrichment analysis. Nucleic Acids Res 2019; 46:W194-W199. [PMID: 29878235 PMCID: PMC6030814 DOI: 10.1093/nar/gky464] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2018] [Accepted: 05/21/2018] [Indexed: 12/21/2022] Open
Abstract
The past few years have seen an explosion of interest in understanding the role of regulatory DNA. This interest has driven large-scale production of functional genomics data and analytical methods. One popular analysis is to test for enrichment of overlaps between a query set of genomic regions and a database of region sets. In this way, new genomic data can be easily connected to annotations from external data sources. Here, we present an interactive interface for enrichment analysis of genomic locus overlaps using a web server called LOLAweb. LOLAweb accepts a set of genomic ranges from the user and tests it for enrichment against a database of region sets. LOLAweb renders results in an R Shiny application to provide interactive visualization features, enabling users to filter, sort, and explore enrichment results dynamically. LOLAweb is built and deployed in a Linux container, making it scalable to many concurrent users on our servers and also enabling users to download and run LOLAweb locally.
Collapse
Affiliation(s)
- V P Nagraj
- School of Medicine Research Computing, University of Virginia, USA
| | - Neal E Magee
- School of Medicine Research Computing, University of Virginia, USA
| | - Nathan C Sheffield
- Center for Public Health Genomics, University of Virginia, USA.,Departments of Public Health Sciences, Biomedical Engineering, and Biochemistry and Molecular Genetics, University of Virginia, USA
| |
Collapse
|
15
|
Dozmorov MG. Disease classification: from phenotypic similarity to integrative genomics and beyond. Brief Bioinform 2019; 20:1769-1780. [DOI: 10.1093/bib/bby049] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2018] [Revised: 05/01/2018] [Indexed: 02/06/2023] Open
Abstract
Abstract
A fundamental challenge of modern biomedical research is understanding how diseases that are similar on the phenotypic level are similar on the molecular level. Integration of various genomic data sets with the traditionally used phenotypic disease similarity revealed novel genetic and molecular mechanisms and blurred the distinction between monogenic (Mendelian) and complex diseases. Network-based medicine has emerged as a complementary approach for identifying disease-causing genes, genetic mediators, disruptions in the underlying cellular functions and for drug repositioning. The recent development of machine and deep learning methods allow for leveraging real-life information about diseases to refine genetic and phenotypic disease relationships. This review describes the historical development and recent methodological advancements for studying disease classification (nosology).
Collapse
Affiliation(s)
- Mikhail G Dozmorov
- Department of Biostatistics, Virginia Commonwealth University, 830 East Main Street, Richmond, VA, USA
| |
Collapse
|
16
|
Stansfield JC, Tran D, Nguyen T, Dozmorov MG. R Tutorial: Detection of Differentially Interacting Chromatin Regions From Multiple Hi-C Datasets. CURRENT PROTOCOLS IN BIOINFORMATICS 2019; 66:e76. [PMID: 31125519 PMCID: PMC6588411 DOI: 10.1002/cpbi.76] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The three-dimensional (3D) interactions of chromatin regulate cell-type-specific gene expression, recombination, X-chromosome inactivation, and many other genomic processes. High-throughput chromatin conformation capture (Hi-C) technologies capture the structure of the chromatin on a global scale by measuring all-vs.-all interactions and can provide new insights into genomic regulation. The workflow presented here describes how to analyze and interpret a comparative Hi-C experiment. We describe the process of obtaining Hi-C data from public repositories and give suggestions for pre-processing pipelines for users who intend to analyze their own raw data. We then describe the data normalization and comparative analysis process. We present three protocols describing the use of the multiHiCcompare, diffHic, and FIND R packages, respectively, to perform a comparative analysis of Hi-C experiments. Finally, visualization of the results and downstream interpretation of the differentially interacting regions are discussed. The bulk of this tutorial uses the R programming environment, and the processes described can be performed with most operating systems and a single computer. © 2019 by John Wiley & Sons, Inc.
Collapse
Affiliation(s)
- John C. Stansfield
- Dept. of Biostatistics, Virginia Commonwealth University, Richmond, VA, 23298, USA
| | - Duc Tran
- Dept. of Computer Science & Engineering, University of Nevada, Reno, NV, 89557, USA
| | - Tin Nguyen
- Dept. of Computer Science & Engineering, University of Nevada, Reno, NV, 89557, USA
| | - Mikhail G. Dozmorov
- Dept. of Biostatistics, Virginia Commonwealth University, Richmond, VA, 23298, USA
| |
Collapse
|
17
|
Oulas A, Minadakis G, Zachariou M, Sokratous K, Bourdakou MM, Spyrou GM. Systems Bioinformatics: increasing precision of computational diagnostics and therapeutics through network-based approaches. Brief Bioinform 2019; 20:806-824. [PMID: 29186305 PMCID: PMC6585387 DOI: 10.1093/bib/bbx151] [Citation(s) in RCA: 69] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2017] [Revised: 10/17/2017] [Indexed: 02/01/2023] Open
Abstract
Systems Bioinformatics is a relatively new approach, which lies in the intersection of systems biology and classical bioinformatics. It focuses on integrating information across different levels using a bottom-up approach as in systems biology with a data-driven top-down approach as in bioinformatics. The advent of omics technologies has provided the stepping-stone for the emergence of Systems Bioinformatics. These technologies provide a spectrum of information ranging from genomics, transcriptomics and proteomics to epigenomics, pharmacogenomics, metagenomics and metabolomics. Systems Bioinformatics is the framework in which systems approaches are applied to such data, setting the level of resolution as well as the boundary of the system of interest and studying the emerging properties of the system as a whole rather than the sum of the properties derived from the system's individual components. A key approach in Systems Bioinformatics is the construction of multiple networks representing each level of the omics spectrum and their integration in a layered network that exchanges information within and between layers. Here, we provide evidence on how Systems Bioinformatics enhances computational therapeutics and diagnostics, hence paving the way to precision medicine. The aim of this review is to familiarize the reader with the emerging field of Systems Bioinformatics and to provide a comprehensive overview of its current state-of-the-art methods and technologies. Moreover, we provide examples of success stories and case studies that utilize such methods and tools to significantly advance research in the fields of systems biology and systems medicine.
Collapse
Affiliation(s)
- Anastasis Oulas
- Bioinformatics European Research Area Chair, The Cyprus Institute of Neurology and Genetics, Nicosia, Cyprus
| | - George Minadakis
- Bioinformatics European Research Area Chair, The Cyprus Institute of Neurology and Genetics, Nicosia, Cyprus
| | - Margarita Zachariou
- Bioinformatics European Research Area Chair, The Cyprus Institute of Neurology and Genetics, Nicosia, Cyprus
| | - Kleitos Sokratous
- Bioinformatics European Research Area Chair, The Cyprus Institute of Neurology and Genetics, Nicosia, Cyprus
| | - Marilena M Bourdakou
- Bioinformatics European Research Area Chair, The Cyprus Institute of Neurology and Genetics, Nicosia, Cyprus
| | - George M Spyrou
- Bioinformatics European Research Area Chair, The Cyprus Institute of Neurology and Genetics, Nicosia, Cyprus
| |
Collapse
|
18
|
Kanduri C, Bock C, Gundersen S, Hovig E, Sandve GK. Colocalization analyses of genomic elements: approaches, recommendations and challenges. Bioinformatics 2019; 35:1615-1624. [PMID: 30307532 PMCID: PMC6499241 DOI: 10.1093/bioinformatics/bty835] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2018] [Revised: 09/03/2018] [Accepted: 10/10/2018] [Indexed: 12/23/2022] Open
Abstract
MOTIVATION Many high-throughput methods produce sets of genomic regions as one of their main outputs. Scientists often use genomic colocalization analysis to interpret such region sets, for example to identify interesting enrichments and to understand the interplay between the underlying biological processes. Although widely used, there is little standardization in how these analyses are performed. Different practices can substantially affect the conclusions of colocalization analyses. RESULTS Here, we describe the different approaches and provide recommendations for performing genomic colocalization analysis, while also discussing common methodological challenges that may influence the conclusions. As illustrated by concrete example cases, careful attention to analysis details is needed in order to meet these challenges and to obtain a robust and biologically meaningful interpretation of genomic region set data. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Chakravarthi Kanduri
- Department of Informatics, University of Oslo, Oslo, Norway
- K. G. Jebsen Coeliac Disease Research Centre, Oslo, Norway
| | - Christoph Bock
- CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria
- Department of Laboratory Medicine, Medical University of Vienna, Vienna, Austria
- Max Planck Institute for Informatics, Saarbrücken, Germany
| | - Sveinung Gundersen
- Department of Informatics, University of Oslo, Oslo, Norway
- Elixir Norway, Oslo Node, University of Oslo, Oslo, Norway
| | - Eivind Hovig
- Department of Informatics, University of Oslo, Oslo, Norway
- Elixir Norway, Oslo Node, University of Oslo, Oslo, Norway
- Department of Tumor Biology, Institute for Cancer Research, Oslo, Norway
- Institute for Cancer Genetics and Informatics, The Norwegian Radium Hospital, Oslo, Norway, UK
| | - Geir Kjetil Sandve
- Department of Informatics, University of Oslo, Oslo, Norway
- K. G. Jebsen Coeliac Disease Research Centre, Oslo, Norway
| |
Collapse
|
19
|
Márquez A, Kerick M, Zhernakova A, Gutierrez-Achury J, Chen WM, Onengut-Gumuscu S, González-Álvaro I, Rodriguez-Rodriguez L, Rios-Fernández R, González-Gay MA, Mayes MD, Raychaudhuri S, Rich SS, Wijmenga C, Martín J. Meta-analysis of Immunochip data of four autoimmune diseases reveals novel single-disease and cross-phenotype associations. Genome Med 2018; 10:97. [PMID: 30572963 PMCID: PMC6302306 DOI: 10.1186/s13073-018-0604-8] [Citation(s) in RCA: 65] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2018] [Accepted: 11/22/2018] [Indexed: 12/15/2022] Open
Abstract
BACKGROUND In recent years, research has consistently proven the occurrence of genetic overlap across autoimmune diseases, which supports the existence of common pathogenic mechanisms in autoimmunity. The objective of this study was to further investigate this shared genetic component. METHODS For this purpose, we performed a cross-disease meta-analysis of Immunochip data from 37,159 patients diagnosed with a seropositive autoimmune disease (11,489 celiac disease (CeD), 15,523 rheumatoid arthritis (RA), 3477 systemic sclerosis (SSc), and 6670 type 1 diabetes (T1D)) and 22,308 healthy controls of European origin using the R package ASSET. RESULTS We identified 38 risk variants shared by at least two of the conditions analyzed, five of which represent new pleiotropic loci in autoimmunity. We also identified six novel genome-wide associations for the diseases studied. Cell-specific functional annotations and biological pathway enrichment analyses suggested that pleiotropic variants may act by deregulating gene expression in different subsets of T cells, especially Th17 and regulatory T cells. Finally, drug repositioning analysis evidenced several drugs that could represent promising candidates for CeD, RA, SSc, and T1D treatment. CONCLUSIONS In this study, we have been able to advance in the knowledge of the genetic overlap existing in autoimmunity, thus shedding light on common molecular mechanisms of disease and suggesting novel drug targets that could be explored for the treatment of the autoimmune diseases studied.
Collapse
Affiliation(s)
- Ana Márquez
- Instituto de Parasitología y Biomedicina “López-Neyra”, CSIC, PTS Granada, Granada, Spain
- Systemic Autoimmune Disease Unit, Instituto de Investigación Biosanitaria de Granada, Granada, Spain
| | - Martin Kerick
- Instituto de Parasitología y Biomedicina “López-Neyra”, CSIC, PTS Granada, Granada, Spain
| | - Alexandra Zhernakova
- Department of Genetics, University of Groningen, University Medical Centre Groningen, Groningen, The Netherlands
| | | | - Wei-Min Chen
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA USA
| | - Suna Onengut-Gumuscu
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA USA
| | | | | | - Raquel Rios-Fernández
- Systemic Autoimmune Diseases Unit, Complejo Hospitalario de Granada, Hospital Campus de la Salud, Granada, Spain
| | - Miguel A. González-Gay
- Epidemiology, Genetics and Atherosclerosis Research Group on Systemic Inflammatory Diseases, IDIVAL, Santander, Spain
| | - Maureen D. Mayes
- Division of Rheumatology and Clinical Immunogenetics, The University of Texas Health Science Center-Houston, Houston, USA
| | - Soumya Raychaudhuri
- Division of Rheumatology, Immunology, and Allergy, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA USA
- Division of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA USA
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA USA
| | - Stephen S. Rich
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA USA
| | - Cisca Wijmenga
- Department of Genetics, University of Groningen, University Medical Centre Groningen, Groningen, The Netherlands
| | - Javier Martín
- Instituto de Parasitología y Biomedicina “López-Neyra”, CSIC, PTS Granada, Granada, Spain
| |
Collapse
|
20
|
Dozmorov MG. Epigenomic annotation-based interpretation of genomic data: from enrichment analysis to machine learning. Bioinformatics 2018; 33:3323-3330. [PMID: 29028263 DOI: 10.1093/bioinformatics/btx414] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2017] [Accepted: 06/22/2017] [Indexed: 12/12/2022] Open
Abstract
Motivation One of the goals of functional genomics is to understand the regulatory implications of experimentally obtained genomic regions of interest (ROIs). Most sequencing technologies now generate ROIs distributed across the whole genome. The interpretation of these genome-wide ROIs represents a challenge as the majority of them lie outside of functionally well-defined protein coding regions. Recent efforts by the members of the International Human Epigenome Consortium have generated volumes of functional/regulatory data (reference epigenomic datasets), effectively annotating the genome with epigenomic properties. Consequently, a wide variety of computational tools has been developed utilizing these epigenomic datasets for the interpretation of genomic data. Results The purpose of this review is to provide a structured overview of practical solutions for the interpretation of ROIs with the help of epigenomic data. Starting with epigenomic enrichment analysis, we discuss leading tools and machine learning methods utilizing epigenomic and 3D genome structure data. The hierarchy of tools and methods reviewed here presents a practical guide for the interpretation of genome-wide ROIs within an epigenomic context. Contact mikhail.dozmorov@vcuhealth.org. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Mikhail G Dozmorov
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA 23298, USA
| |
Collapse
|
21
|
Gensterblum-Miller E, Wu W, Sawalha AH. Novel Transcriptional Activity and Extensive Allelic Imbalance in the Human MHC Region. THE JOURNAL OF IMMUNOLOGY 2018; 200:1496-1503. [PMID: 29311362 DOI: 10.4049/jimmunol.1701061] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/02/2017] [Accepted: 12/07/2017] [Indexed: 12/19/2022]
Abstract
The MHC region encodes HLA genes and is the most complex region in the human genome. The extensively polymorphic nature of the HLA hinders accurate localization and functional assessment of disease risk loci within this region. Using targeted capture sequencing and constructing individualized genomes for transcriptome alignment, we identified 908 novel transcripts within the human MHC region. These include 593 novel isoforms of known genes, 137 antisense strand RNAs, 119 novel long intergenic noncoding RNAs, and 5 transcripts of 3 novel putative protein-coding human endogenous retrovirus genes. We revealed allele-dependent expression imbalance involving 88% of all heterozygous transcribed single nucleotide polymorphisms throughout the MHC transcriptome. Among these variants, the genetic variant associated with Behçet's disease in the HLA-B/MICA region, which tags HLA-B*51, is within novel long intergenic noncoding RNA transcripts that are exclusively expressed from the haplotype with the protective but not the disease risk allele. Further, the transcriptome within the MHC region can be defined by 14 distinct coexpression clusters, with evidence of coregulation by unique transcription factors in at least 9 of these clusters. Our data suggest a very complex regulatory map of the human MHC, and can help uncover functional consequences of disease risk loci in this region.
Collapse
Affiliation(s)
| | - Weisheng Wu
- Biomedical Research Core Facilities, Bioinformatics Core, University of Michigan, Ann Arbor, MI 48109; and
| | - Amr H Sawalha
- Division of Rheumatology, Department of Internal Medicine, University of Michigan, Ann Arbor, MI 48109; .,Center for Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48109
| |
Collapse
|
22
|
Dozmorov MG, Coit P, Maksimowicz-McKinnon K, Sawalha AH. Age-associated DNA methylation changes in naive CD4 + T cells suggest an evolving autoimmune epigenotype in aging T cells. Epigenomics 2017; 9:429-445. [PMID: 28322571 DOI: 10.2217/epi-2016-0143] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
AIM We sought to define age-associated DNA methylation changes in naive CD4+ T cells. MATERIALS & METHODS Naive CD4+ T cells were collected from 74 healthy individuals (age 19-66 years), and age-related DNA methylation changes were characterized. RESULTS We identified 11,431 age-associated CpG sites, 57% of which were hypermethylated with age. Hypermethylated sites were enriched in CpG islands and repressive transcription factor binding sites, while hypomethylated sites showed T cell specific enrichment in active enhancers marked by H3K27ac and H3K4me1. Our data emphasize cancer-related DNA methylation changes with age, and also reveal age-associated hypomethylation in immune-related pathways, such as T cell receptor signaling, FCγR-mediated phagocytosis, apoptosis and the mammalian target of rapamycin signaling pathway. The MAPK signaling pathway was hypermethylated with age, consistent with a defective MAPK signaling in aging T cells. CONCLUSION Age-associated DNA methylation changes may alter regulatory mechanisms and signaling pathways that predispose to autoimmunity.
Collapse
Affiliation(s)
- Mikhail G Dozmorov
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Patrick Coit
- Division of Rheumatology, Department of Internal Medicine, University of Michigan, Ann Arbor, MI 48109, USA
| | | | - Amr H Sawalha
- Division of Rheumatology, Department of Internal Medicine, University of Michigan, Ann Arbor, MI 48109, USA.,Center for Computational Medicine & Bioinformatics, University of Michigan, Ann Arbor, MI 48109, USA
| |
Collapse
|
23
|
Candidate SNP Markers of Chronopathologies Are Predicted by a Significant Change in the Affinity of TATA-Binding Protein for Human Gene Promoters. BIOMED RESEARCH INTERNATIONAL 2016; 2016:8642703. [PMID: 27635400 PMCID: PMC5011241 DOI: 10.1155/2016/8642703] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/04/2016] [Revised: 06/25/2016] [Accepted: 06/28/2016] [Indexed: 01/14/2023]
Abstract
Variations in human genome (e.g., single nucleotide polymorphisms, SNPs) may be associated with hereditary diseases, their complications, comorbidities, and drug responses. Using Web service SNP_TATA_Comparator presented in our previous paper, here we analyzed immediate surroundings of known SNP markers of diseases and identified several candidate SNP markers that can significantly change the affinity of TATA-binding protein for human gene promoters, with circadian consequences. For example, rs572527200 may be related to asthma, where symptoms are circadian (worse at night), and rs367732974 may be associated with heart attacks that are characterized by a circadian preference (early morning). By the same method, we analyzed the 90 bp proximal promoter region of each protein-coding transcript of each human gene of the circadian clock core. This analysis yielded 53 candidate SNP markers, such as rs181985043 (susceptibility to acute Q fever in male patients), rs192518038 (higher risk of a heart attack in patients with diabetes), and rs374778785 (emphysema and lung cancer in smokers). If they are properly validated according to clinical standards, these candidate SNP markers may turn out to be useful for physicians (to select optimal treatment for each patient) and for the general population (to choose a lifestyle preventing possible circadian complications of diseases).
Collapse
|