1
|
Rudra P, Zhou YH, Nobel A, Wright FA. Control of false discoveries in grouped hypothesis testing for eQTL data. BMC Bioinformatics 2024; 25:147. [PMID: 38605284 PMCID: PMC11007981 DOI: 10.1186/s12859-024-05736-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Accepted: 03/08/2024] [Indexed: 04/13/2024] Open
Abstract
BACKGROUND Expression quantitative trait locus (eQTL) analysis aims to detect the genetic variants that influence the expression of one or more genes. Gene-level eQTL testing forms a natural grouped-hypothesis testing strategy with clear biological importance. Methods to control family-wise error rate or false discovery rate for group testing have been proposed earlier, but may not be powerful or easily apply to eQTL data, for which certain structured alternatives may be defensible and may enable the researcher to avoid overly conservative approaches. RESULTS In an empirical Bayesian setting, we propose a new method to control the false discovery rate (FDR) for grouped hypotheses. Here, each gene forms a group, with SNPs annotated to the gene corresponding to individual hypotheses. The heterogeneity of effect sizes in different groups is considered by the introduction of a random effects component. Our method, entitled Random Effects model and testing procedure for Group-level FDR control (REG-FDR), assumes a model for alternative hypotheses for the eQTL data and controls the FDR by adaptive thresholding. As a convenient alternate approach, we also propose Z-REG-FDR, an approximate version of REG-FDR, that uses only Z-statistics of association between genotype and expression for each gene-SNP pair. The performance of Z-REG-FDR is evaluated using both simulated and real data. Simulations demonstrate that Z-REG-FDR performs similarly to REG-FDR, but with much improved computational speed. CONCLUSION Our results demonstrate that the Z-REG-FDR method performs favorably compared to other methods in terms of statistical power and control of FDR. It can be of great practical use for grouped hypothesis testing for eQTL analysis or similar problems in statistical genomics due to its fast computation and ability to be fit using only summary data.
Collapse
Affiliation(s)
- Pratyaydipta Rudra
- Department of Statistics, Oklahoma State University, Stillwater, OK, USA.
| | - Yi-Hui Zhou
- Bioinformatics Research Center, Departments of Statistics and Biological Sciences, North Carolina State University, Raleigh, NC, USA
| | - Andrew Nobel
- Department of Statistics and Operations Research, University of North Carolina, Chapel Hill, NC, USA
| | - Fred A Wright
- Bioinformatics Research Center, Departments of Statistics and Biological Sciences, North Carolina State University, Raleigh, NC, USA.
| |
Collapse
|
2
|
Pandey GK, Vadlamudi S, Currin KW, Moxley AH, Nicholas JC, McAfee JC, Broadaway KA, Mohlke KL. Liver regulatory mechanisms of noncoding variants at lipid and metabolic trait loci. HGG Adv 2024; 5:100275. [PMID: 38297830 PMCID: PMC10881423 DOI: 10.1016/j.xhgg.2024.100275] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2023] [Revised: 01/23/2024] [Accepted: 01/23/2024] [Indexed: 02/02/2024] Open
Abstract
Genome-wide association studies (GWASs) have identified hundreds of risk loci for liver disease and lipid-related metabolic traits, although identifying their target genes and molecular mechanisms remains challenging. We predicted target genes at GWAS signals by integrating them with molecular quantitative trait loci for liver gene expression (eQTL) and liver chromatin accessibility QTL (caQTL). We predicted specific regulatory caQTL variants at four GWAS signals located near EFHD1, LITAF, ZNF329, and GPR180. Using transcriptional reporter assays, we determined that caQTL variants rs13395911, rs11644920, rs34003091, and rs9556404 exhibit allelic differences in regulatory activity. We also performed a protein binding assay for rs13395911 and found that FOXA2 differentially interacts with the alleles of rs13395911. For variants rs13395911 and rs11644920 in putative enhancer regulatory elements, we used CRISPRi to demonstrate that repression of the enhancers altered the expression of the predicted target and/or nearby genes. Repression of the element at rs13395911 reduced the expression of EFHD1, and repression of the element at rs11644920 reduced the expression of LITAF, SNN, and TXNDC11. Finally, we showed that EFHD1 is a metabolically active gene in HepG2 cells. Together, these results provide key steps to connect genetic variants with cellular mechanisms and help elucidate the causes of liver disease.
Collapse
Affiliation(s)
- Gautam K Pandey
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27599, USA
| | | | - Kevin W Currin
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27599, USA
| | - Anne H Moxley
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27599, USA
| | - Jayna C Nicholas
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27599, USA
| | - Jessica C McAfee
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27599, USA; UNC Neuroscience Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - K Alaine Broadaway
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27599, USA
| | - Karen L Mohlke
- Department of Genetics, University of North Carolina, Chapel Hill, NC 27599, USA.
| |
Collapse
|
3
|
Tsouris A, Brach G, Friedrich A, Hou J, Schacherer J. Diallel panel reveals a significant impact of low-frequency genetic variants on gene expression variation in yeast. Mol Syst Biol 2024; 20:362-373. [PMID: 38355920 PMCID: PMC10987670 DOI: 10.1038/s44320-024-00021-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Revised: 01/29/2024] [Accepted: 01/30/2024] [Indexed: 02/16/2024] Open
Abstract
Unraveling the genetic sources of gene expression variation is essential to better understand the origins of phenotypic diversity in natural populations. Genome-wide association studies identified thousands of variants involved in gene expression variation, however, variants detected only explain part of the heritability. In fact, variants such as low-frequency and structural variants (SVs) are poorly captured in association studies. To assess the impact of these variants on gene expression variation, we explored a half-diallel panel composed of 323 hybrids originated from pairwise crosses of 26 natural Saccharomyces cerevisiae isolates. Using short- and long-read sequencing strategies, we established an exhaustive catalog of single nucleotide polymorphisms (SNPs) and SVs for this panel. Combining this dataset with the transcriptomes of all hybrids, we comprehensively mapped SNPs and SVs associated with gene expression variation. While SVs impact gene expression variation, SNPs exhibit a higher effect size with an overrepresentation of low-frequency variants compared to common ones. These results reinforce the importance of dissecting the heritability of complex traits with a comprehensive catalog of genetic variants at the population level.
Collapse
Affiliation(s)
- Andreas Tsouris
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Gauthier Brach
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Anne Friedrich
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Jing Hou
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France.
| | - Joseph Schacherer
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France.
- Institut Universitaire de France (IUF), Paris, France.
| |
Collapse
|
4
|
Zachariou M, Loizidou EM, Spyrou GM. Immediate-Early Genes as Influencers in Genetic Networks and their Role in Alzheimer's Disease. bioRxiv 2024:2024.03.29.586739. [PMID: 38585978 PMCID: PMC10996630 DOI: 10.1101/2024.03.29.586739] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/09/2024]
Abstract
Immediate-early genes (IEGs) are a class of activity-regulated genes (ARGs) that are transiently and rapidly activated in the absence of de novo protein synthesis in response to neuronal activity. We explored the role of IEGs in genetic networks to pinpoint potential drug targets for Alzheimer's disease (AD). Using a combination of network analysis and genome-wide association study (GWAS) summary statistics we show that (1) IEGs exert greater topological influence across different human and mouse gene networks compared to other ARGs, (2) ARGs are sparsely involved in diseases and significantly more mutational constrained compared to non-ARGs, (3) Many AD-linked variants are in ARGs gene regions, mainly in MARK4 near FOSB, with an AD risk eQTL that increases MARK4 expression in cortical areas, (4) MARK4 holds an influential place in a dense AD multi-omic network and a high AD druggability score. Our work on IEGs' influential network role is a valuable contribution to guiding interventions for diseases marked by dysregulation of their downstream targets and highlights MARK4 as a promising underexplored AD-target.
Collapse
Affiliation(s)
| | - Eleni M Loizidou
- biobank.cy, Center of Excellence in Biobanking and Biomedical Research, University of Cyprus
| | - George M Spyrou
- Bioinformatics Department, The Cyprus Institute of Neurology and Genetics
| |
Collapse
|
5
|
Chen F, Zhang Y, Sedlazeck FJ, Creighton CJ. Germline structural variation globally impacts the cancer transcriptome including disease-relevant genes. Cell Rep Med 2024; 5:101446. [PMID: 38442712 PMCID: PMC10983041 DOI: 10.1016/j.xcrm.2024.101446] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Revised: 01/01/2024] [Accepted: 02/06/2024] [Indexed: 03/07/2024]
Abstract
Germline variation and somatic alterations contribute to the molecular profile of cancers. We combine RNA with whole genome sequencing across 1,218 cancer patients to determine the extent germline structural variants (SVs) impact expression of nearby genes. For hundreds of genes, recurrent and common germline SV breakpoints within 100 kb associate with increased or decreased expression in tumors spanning various tissues of origin. A significant fraction of germline SV expression associations involves duplication of intergenic enhancers or 3' UTR disruption. Genes altered by both somatic and germline SVs include ATRX and CEBPA. Genes essential in cancer cell lines include BARD1 and IRS2. Genes with both expression and germline SV breakpoint patterns associated with patient survival include GCLM. Our results capture a class of phenotypic variation at work in the disease setting, including genes with cancer roles. Specific germline SVs represent potential cancer risk variants for genetic testing, including those involving genes with targeting implications.
Collapse
Affiliation(s)
- Fengju Chen
- Dan L. Duncan Comprehensive Cancer Center, Baylor College of Medicine, Houston, TX 77030, USA
| | - Yiqun Zhang
- Dan L. Duncan Comprehensive Cancer Center, Baylor College of Medicine, Houston, TX 77030, USA
| | - Fritz J Sedlazeck
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX 77030, USA; Department of Computer Science, Rice University, Houston, TX 77005, USA
| | - Chad J Creighton
- Dan L. Duncan Comprehensive Cancer Center, Baylor College of Medicine, Houston, TX 77030, USA; Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX 77030, USA; Department of Medicine, Baylor College of Medicine, Houston, TX 77030, USA.
| |
Collapse
|
6
|
Zhao S, Xie Y, Ding X, Zheng C, Chen J, Zhao N, Ji Y, Wang Q, Liu Y, Cheng C. Exploring the causal relationship between antihypertensive drugs and glioblastoma by combining drug target Mendelian randomization study, eQTL colocalization, and single-cell RNA sequencing. Environ Toxicol 2024. [PMID: 38450887 DOI: 10.1002/tox.24210] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/09/2024] [Revised: 02/20/2024] [Accepted: 02/25/2024] [Indexed: 03/08/2024]
Abstract
Recent reports indicate a potential oncogenic role of antihypertensive drugs in common cancers. However, it remains uncertain whether this phenomenon influences the risk of glioblastoma multiforme (GBM). This study aimed to assess the potential causal effects of blood pressure (BP) and antihypertensive drugs on GBM. Genome-wide association study (GWAS) summary statistics for systolic blood pressure (SBP), diastolic blood pressure (DBP), and GBM in Europeans were downloaded. To represent the effects of antihypertensive drugs, we utilized single nucleotide polymorphisms (SNPs) associated with SBP/DBP adjacent to the coding regions of different antihypertensive drugs as instrumental variables to model five antihypertensive drugs, including angiotensin-converting enzyme inhibitors, angiotensin receptor blockers, calcium channel blockers, β-receptor blockers (BBs), and thiazide diuretics. Positive control studies were performed using GWAS data in chronic heart failure. The primary method for causality estimation was the inverse-variance-weighted method. Mendelian randomization analysis showed that BBs with the β1-adrenergic receptor (ADRB1) as a therapeutic target could significantly reduce the risk of GBM by mediating DBP (OR = 0.431, 95% CI: 0.267-0.697, p < .001) and that they could also significantly reduce the risk of GBM by mediating SBP (OR = 0.595, 95% CI: 0.422-0.837, p = .003). Sensitivity analysis and colocalization analysis reinforced the robustness of these findings. Finally, the low expression of the ADRB1 gene in malignant gliomas was found by GBM data from TCGA and single-cell RNA sequencing, which most likely contributed to the poor prognosis of GBM patients. In summary, our study provides preliminary evidence of some causal relationship between ADRB1-targeted BBs and glioblastoma development. However, more studies are needed to validate these findings and further reveal the complex relationship between BP and GBM.
Collapse
Affiliation(s)
- Songyun Zhao
- Department of Neurosurgery, The Affiliated Wuxi People's Hospital of Nanjing Medical University, Wuxi, China
- Wuxi Medical Center of Nanjing Medical University, Wuxi, China
| | - Yi Xie
- Department of Burns and Plastic Surgery, Nanjing Drum Tower Hospital, The Affiliated Hospital of Nanjing University Medical School, Nanjing, China
| | - Xu Ding
- Wuxi Medical Center of Nanjing Medical University, Wuxi, China
| | - Chuanhua Zheng
- Department of Neurosurgery, The Affiliated Hospital of Youjiang Medical University for Nationalities, Baise, China
| | - Jiaxing Chen
- Department of Neurosurgery, The Affiliated Hospital of Youjiang Medical University for Nationalities, Baise, China
| | - Ning Zhao
- Wuxi Medical Center of Nanjing Medical University, Wuxi, China
| | - Yi Ji
- Wuxi Medical Center of Nanjing Medical University, Wuxi, China
| | - Qi Wang
- Department of Gastroenterology, The Affiliated Hospital of Jiangsu University, Zhenjiang, China
| | - Yuankun Liu
- Department of Neurosurgery, The Affiliated Wuxi People's Hospital of Nanjing Medical University, Wuxi, China
- Wuxi Medical Center of Nanjing Medical University, Wuxi, China
| | - Chao Cheng
- Department of Neurosurgery, The Affiliated Wuxi People's Hospital of Nanjing Medical University, Wuxi, China
- Wuxi Medical Center of Nanjing Medical University, Wuxi, China
| |
Collapse
|
7
|
Johnson M, Chelysheva I, Öner D, McGinley J, Lin GL, O'Connor D, Robinson H, Drysdale SB, Gammin E, Vernon S, Muller J, Wolfenden H, Westcar S, Anguvaa L, Thwaites RS, Bont L, Wildenbeest J, Martinón-Torres F, Aerssens J, Openshaw PJM, Pollard AJ. A Genome-Wide Association Study of Respiratory Syncytial Virus Infection Severity in Infants. J Infect Dis 2024; 229:S112-S119. [PMID: 38271230 DOI: 10.1093/infdis/jiae029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2023] [Revised: 01/16/2024] [Accepted: 01/20/2024] [Indexed: 01/27/2024] Open
Abstract
BACKGROUND Respiratory syncytial virus (RSV) is a significant cause of infant morbidity and mortality worldwide. Most children experience at least one 1 RSV infection by the age of two 2 years, but not all develop severe disease. However, the understanding of genetic risk factors for severe RSV is incomplete. Consequently, we conducted a genome-wide association study of RSV severity. METHODS Disease severity was assessed by the ReSVinet scale, in a cohort of 251 infants aged 1 week to 1 year. Genotyping data were collected from multiple European study sites as part of the RESCEU Consortium. Linear regression models were used to assess the impact of genotype on RSV severity and gene expression as measured by microarray. RESULTS While no SNPs reached the genome-wide statistical significance threshold (P < 5 × 10-8), we identified 816 candidate SNPs with a P-value of <1 × 10-4. Functional annotation of candidate SNPs highlighted genes relevant to neutrophil trafficking and cytoskeletal functions, including LSP1 and RAB27A. Moreover, SNPs within the RAB27A locus significantly altered gene expression (false discovery rate, FDR P < .05). CONCLUSIONS These findings may provide insights into genetic mechanisms driving severe RSV infection, offering biologically relevant information for future investigations.
Collapse
Affiliation(s)
- Mari Johnson
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford
- NIHR Oxford Biomedical Research Centre and Oxford University Hospitals NHS Foundation Trust, United Kingdom
| | - Irina Chelysheva
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford
- NIHR Oxford Biomedical Research Centre and Oxford University Hospitals NHS Foundation Trust, United Kingdom
| | - Deniz Öner
- Biomarkers Infectious Diseases, Janssen Pharmaceutica NV, Beerse, Belgium
| | - Joseph McGinley
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford
- NIHR Oxford Biomedical Research Centre and Oxford University Hospitals NHS Foundation Trust, United Kingdom
| | - Gu-Lung Lin
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford
- NIHR Oxford Biomedical Research Centre and Oxford University Hospitals NHS Foundation Trust, United Kingdom
| | - Daniel O'Connor
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford
- NIHR Oxford Biomedical Research Centre and Oxford University Hospitals NHS Foundation Trust, United Kingdom
| | - Hannah Robinson
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford
- NIHR Oxford Biomedical Research Centre and Oxford University Hospitals NHS Foundation Trust, United Kingdom
| | - Simon B Drysdale
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford
- NIHR Oxford Biomedical Research Centre and Oxford University Hospitals NHS Foundation Trust, United Kingdom
| | - Emma Gammin
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford
- NIHR Oxford Biomedical Research Centre and Oxford University Hospitals NHS Foundation Trust, United Kingdom
| | - Sophie Vernon
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford
- NIHR Oxford Biomedical Research Centre and Oxford University Hospitals NHS Foundation Trust, United Kingdom
| | - Jill Muller
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford
- NIHR Oxford Biomedical Research Centre and Oxford University Hospitals NHS Foundation Trust, United Kingdom
| | | | | | | | - Ryan S Thwaites
- National Heart and Lung Institute, Imperial College London, United Kingdom
| | - Louis Bont
- Department of Paediatric Infectious Diseases and Immunology, Wilhelmina Children's Hospital, University Medical Centre Utrecht, Netherlands
| | - Joanne Wildenbeest
- Department of Paediatric Infectious Diseases and Immunology, Wilhelmina Children's Hospital, University Medical Centre Utrecht, Netherlands
| | - Federico Martinón-Torres
- Translational Pediatrics and Infectious Diseases, Department of Pediatrics, Hospital Clínico Universitario de Santiago de Compostela
- Genetics, Vaccines and Infections Research Group, Instituto de Investigación Sanitaria de Santiago, Universidade de Santiago de Compostela
- Centro de Investigación Biomédica en Red de Enfermedades Respiratorias, Instituto de Salud Carlos III, Madrid, Spain
| | - Jeroen Aerssens
- Biomarkers Infectious Diseases, Janssen Pharmaceutica NV, Beerse, Belgium
| | - Peter J M Openshaw
- National Heart and Lung Institute, Imperial College London, United Kingdom
| | - Andrew J Pollard
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford
- NIHR Oxford Biomedical Research Centre and Oxford University Hospitals NHS Foundation Trust, United Kingdom
| |
Collapse
|
8
|
Qiao L, Lv S, Meng K, Yang J. Genetically proxied therapeutic inhibition of lipid-lowering drug targets and risk of rheumatoid arthritis disease: a Mendelian randomization study. Clin Rheumatol 2024; 43:939-947. [PMID: 38198113 DOI: 10.1007/s10067-023-06837-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2023] [Revised: 11/12/2023] [Accepted: 11/27/2023] [Indexed: 01/11/2024]
Abstract
OBJECTIVE To evaluate the potential impact of consistent use of similar treatments over a long period; it is essential to investigate the potential correlation between genetic variations that influence the expression or function of pharmacological targets for reducing lipid levels and the risk of developing rheumatoid arthritis. METHODS We used variants in the following genes to conduct Mendelian randomization analyses: HMGCR (encoding the target for statins), PCSK9 (encoding the target for PCSK9 inhibitors, such as evolocumab and alirocumab), and NPC1L1 (encoding the target for ezetimibe). Data from lipid genetics consortia (173,082 sample size) were used to weight variations according to their correlations with low-density lipoprotein cholesterol (LDL-C). In two large datasets (total n = 19,562 cases, 501,655 controls). We conducted a meta-analysis of Mendelian randomization estimates, weighted by LDL-C levels, on the regional differences in the risk of rheumatoid arthritis using data from two large databases. RESULTS We approached SMR and IVW-MR analyses to examine the relationship between target gene expression (including HMGCR, PCSK9, and NPC1L1) and LDL-C levels mediated by these genes with RA. The IVW-MR analysis revealed no significant association between genetically predicted LDL-C concentration and the risk of RA (OR = 0.88, 95% CI = 0.59-1.29; OR = 0.91, 95% CI = 0.67-1.23; OR = 0.81, 95% CI = 0.49-1.36; all p > 0.05). Similarly, our findings from the SMR approach provided no evidence to suggest that gene expression of HMGCR, PCSK9, and NPC1L1 was associated with the risk of RA (OR = 0.91, 95% CI = 0.79-1.05, p = 0.207; OR = 0.96, 95% CI = 0.85-1.09, p = 0.493). CONCLUSIONS Our results do not provide evidence to support the hypothesis that reducing LDL-C levels with statins, alirocumab, or ezetimibe effectively prevents the risk of developing RA. However, our study provides valuable insights into the assessment of lipid-lowering agents in RA, which can enhance our understanding of the condition and assist in clinical practice by aiding in the determination and monitoring of RA status to clinical response.
Collapse
Affiliation(s)
- Liang Qiao
- Shanghai Xuhui District Central Hospital, Shanghai, China
| | - Shun Lv
- Shanghai Xuhui District Central Hospital, Shanghai, China
| | - Kai Meng
- Shanghai Xuhui District Central Hospital, Shanghai, China
| | - Jianmei Yang
- Shanghai Xuhui District Central Hospital, Shanghai, China.
| |
Collapse
|
9
|
Niciura SCM, Cardoso TF, Ibelli AMG, Okino CH, Andrade BG, Benavides MV, Chagas ACDS, Esteves SN, Minho AP, Regitano LCDA, Gondro C. Multi-omics data elucidate parasite-host-microbiota interactions and resistance to Haemonchus contortus in sheep. Parasit Vectors 2024; 17:102. [PMID: 38429820 PMCID: PMC10908167 DOI: 10.1186/s13071-024-06205-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2023] [Accepted: 02/18/2024] [Indexed: 03/03/2024] Open
Abstract
BACKGROUND The integration of molecular data from hosts, parasites, and microbiota can enhance our understanding of the complex biological interactions underlying the resistance of hosts to parasites. Haemonchus contortus, the predominant sheep gastrointestinal parasite species in the tropics, causes significant production and economic losses, which are further compounded by the diminishing efficiency of chemical control owing to anthelmintic resistance. Knowledge of how the host responds to infection and how the parasite, in combination with microbiota, modulates host immunity can guide selection decisions to breed animals with improved parasite resistance. This understanding will help refine management practices and advance the development of new therapeutics for long-term helminth control. METHODS Eggs per gram (EPG) of feces were obtained from Morada Nova sheep subjected to two artificial infections with H. contortus and used as a proxy to select animals with high resistance or susceptibility for transcriptome sequencing (RNA-seq) of the abomasum and 50 K single-nucleotide genotyping. Additionally, RNA-seq data for H. contortus were generated, and amplicon sequence variants (ASV) were obtained using polymerase chain reaction amplification and sequencing of bacterial and archaeal 16S ribosomal RNA genes from sheep feces and rumen content. RESULTS The heritability estimate for EPG was 0.12. GAST, GNLY, IL13, MGRN1, FGF14, and RORC genes and transcripts were differentially expressed between resistant and susceptible animals. A genome-wide association study identified regions on chromosomes 2 and 11 that harbor candidate genes for resistance, immune response, body weight, and adaptation. Trans-expression quantitative trait loci were found between significant variants and differentially expressed transcripts. Functional co-expression modules based on sheep genes and ASVs correlated with resistance to H. contortus, showing enrichment in pathways of response to bacteria, immune and inflammatory responses, and hub features of the Christensenellaceae, Bacteroides, and Methanobrevibacter genera; Prevotellaceae family; and Verrucomicrobiota phylum. In H. contortus, some mitochondrial, collagen-, and cuticle-related genes were expressed only in parasites isolated from susceptible sheep. CONCLUSIONS The present study identified chromosome regions, genes, transcripts, and pathways involved in the elaborate interactions between the sheep host, its gastrointestinal microbiota, and the H. contortus parasite. These findings will assist in the development of animal selection strategies for parasite resistance and interdisciplinary approaches to control H. contortus infection in sheep.
Collapse
|
10
|
Stone K, Platig J, Quackenbush J, Fagny M. Complex Traits Heritability is Highly Clustered in the eQTL Bipartite Network. bioRxiv 2024:2024.02.27.582063. [PMID: 38464142 PMCID: PMC10925220 DOI: 10.1101/2024.02.27.582063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/12/2024]
Abstract
Single Nucleotide Polymorphisms (SNPs) associated with traits typically explain a small part of the trait genetic heritability-with the remainder thought to be distributed throughout the genome. Such SNPs are likely to alter expression levels of biologically relevant genes. Expression Quantitative Trait Locus (eQTL) networks analysis has helped to functionally characterize such variants. We systematically analyze the distribution of SNP heritability for ten traits across 29 tissue-specific eQTL networks. We find that heritability is clustered in a small number or tissue-specific, functionally relevant SNP-gene modules and that the greatest occurs in local "hubs" that are both the cornerstone of the network's modules and tissue-specific regulatory elements. The network structure could thus both amplify the genotype-phenotype connection and buffer the deleterious effect of the genetic variations on other traits. Together, these results define a conceptual framework for understanding complex trait architecture and identifying key mutations carrying most of the heritability.
Collapse
Affiliation(s)
- Katherine Stone
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States
- Department of Data Science and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts, USA
| | - John Platig
- Center for Public Health Genomics, University of Virginia, Charlottesville, Virginia, USA
- Department of Public Health Sciences, University of Virginia, Charlottesville, Virginia, USA
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
| | - John Quackenbush
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States
- Department of Data Science and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts, USA
- Channing Division of Network Medicine, Brigham and Women's Hospital, Boston, Massachusetts, United States
| | - Maud Fagny
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States
- Department of Data Science and Center for Cancer Computational Biology, Dana-Farber Cancer Institute, Boston, Massachusetts, USA
- Université Paris-Saclay, INRAE, CNRS, AgroParisTech, Genetique Quantitative et Evolution - Le Moulon, Gif-sur-Yvette 91190 France
| |
Collapse
|
11
|
Gomez-Cano F, Rodriguez J, Zhou P, Chu YH, Magnusson E, Gomez-Cano L, Krishnan A, Springer NM, de Leon N, Grotewold E. Prioritizing Metabolic Gene Regulators through Multi-Omic Network Integration in Maize. bioRxiv 2024:2024.02.26.582075. [PMID: 38464086 PMCID: PMC10925184 DOI: 10.1101/2024.02.26.582075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/12/2024]
Abstract
Elucidating gene regulatory networks (GRNs) is a major area of study within plant systems biology. Phenotypic traits are intricately linked to specific gene expression profiles. These expression patterns arise primarily from regulatory connections between sets of transcription factors (TFs) and their target genes. In this study, we integrated publicly available co-expression networks derived from more than 6,000 RNA-seq samples, 283 protein-DNA interaction assays, and 16 million of SNPs used to identify expression quantitative loci (eQTL), to construct TF-target networks. In total, we analyzed ~4.6M interactions to generate four distinct types of TF-target networks: co-expression, protein-DNA interaction (PDI), trans-expression quantitative loci (trans-eQTL), and cis-eQTL combined with PDIs. To improve the functional annotation of TFs based on its target genes, we implemented three different strategies to integrate these four types of networks. We subsequently evaluated the effectiveness of our method through loss-of function mutant and random networks. The multi-network integration allowed us to identify transcriptional regulators of hormone-, metabolic- and development-related processes. Finally, using the topological properties of the fully integrated network, we identified potentially functional redundant TF paralogs. Our findings retrieved functions previously documented for numerous TFs and revealed novel functions that are crucial for informing the design of future experiments. The approach here-described lays the foundation for the integration of multi-omic datasets in maize and other plant systems.
Collapse
Affiliation(s)
- Fabio Gomez-Cano
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824-6473, USA
- Current address: Department of Molecular, Cellular, and Development Biology, University of Michigan, Ann Arbor, MI, 48109, USA
| | - Jonas Rodriguez
- Department of Plant and Agroecosystem Sciences, University of Wisconsin Madison, Madison, WI 53706, USA
| | - Peng Zhou
- Department of Plant and Microbial Biology, University of Minnesota, Saint Paul, MN 55108
| | - Yi-Hsuan Chu
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824-6473, USA
| | - Erika Magnusson
- Department of Plant and Microbial Biology, University of Minnesota, Saint Paul, MN 55108
| | - Lina Gomez-Cano
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824-6473, USA
| | - Arjun Krishnan
- Department of Biomedical Informatics, University of Colorado Anschutz Medical Campus, Aurora, CO 80045, USA
| | - Nathan M Springer
- Department of Plant and Microbial Biology, University of Minnesota, Saint Paul, MN 55108
- Current address: Global Breeding, Bayer Crop Sciences, Chesterfield MO 63017, USA
| | - Natalia de Leon
- Department of Plant and Agroecosystem Sciences, University of Wisconsin Madison, Madison, WI 53706, USA
| | - Erich Grotewold
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824-6473, USA
| |
Collapse
|
12
|
Zhang C, Huang W, Niu W, Yang H, Zheng Y, Gao X, Qiu X. Five genes identified as prognostic markers for colorectal cancer through the integration of genome-wide association study and expression quantitative trait loci data. Per Med 2024. [PMID: 38380524 DOI: 10.2217/pme-2023-0103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2024]
Abstract
Background: Colorectal cancer (CRC) is a prominent form of cancer globally, ranking second in terms of prevalence and serving as a leading cause of cancer-related deaths, but the underlying biological interpretation remains largely unknown. Methods: We used the summary data-based Mendelian randomization method to integrate CRC genome-wide association studies (ncase = 7062; ncontrol = 195,745) and expression quantitative trait loci summary data in peripheral whole blood (Consortium for Architecture of Gene Expression: n = 2765; Genotype-Tissue Expression [v8]: n = 755) and colon tissue (colon-transverse: n = 406; colon-sigmoid: n = 373) and identified related genes. Results: Genes ABTB1, CYP21A2, NLRP1, PHKG1 and PIP5K1C have emerged as significant prognostic markers for CRC patient survival. Functional analysis revealed their involvement in cancer cell migration and invasion mechanisms, providing valuable insights for the development of future anti-CRC drugs. Conclusion: We successfully identified five CRC risk genes, providing new insights and research directions for the effective mechanisms of CRC.
Collapse
Affiliation(s)
- Cuizhen Zhang
- Department of Pharmacy, Huashan Hospital, Fudan University, Shanghai, 200040, China
| | - Wenjie Huang
- Department of Pharmacy, Huashan Hospital, Fudan University, Shanghai, 200040, China
| | - Wanjie Niu
- Department of Pharmacy, Huashan Hospital, Fudan University, Shanghai, 200040, China
| | - Huiying Yang
- Department of Pharmacy, Huashan Hospital, Fudan University, Shanghai, 200040, China
| | - Yingyi Zheng
- Department of Pharmacy, Huashan Hospital, Fudan University, Shanghai, 200040, China
| | - Xuan Gao
- Outpatient and Emergency Management Office, National Children's Medical Center, Children's Hospital of Fudan University, Shanghai, 201102, China
| | - Xiaoyan Qiu
- Department of Pharmacy, Huashan Hospital, Fudan University, Shanghai, 200040, China
| |
Collapse
|
13
|
Yang G, Pan Y, Pan W, Song Q, Zhang R, Tong W, Cui L, Ji W, Song W, Song B, Deng P, Nie X. Combined GWAS and eGWAS reveals the genetic basis underlying drought tolerance in emmer wheat (Triticum turgidum L.). New Phytol 2024. [PMID: 38358006 DOI: 10.1111/nph.19589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Accepted: 01/25/2024] [Indexed: 02/16/2024]
Abstract
Drought is one of the major environmental constraints for wheat production world-wide. As the progenitor and genetic reservoir of common wheat, emmer wheat is considered as an invaluable gene pool for breeding drought-tolerant wheat. Combining GWAS and eGWAS analysis of 107 accessions, we identified 86 QTLs, 105 462 eQTLs as well as 68 eQTL hotspots associating with drought tolerance (DT) in emmer wheat. A complex regulatory network composed of 185 upstream regulator and 2432 downstream drought-responsive candidates was developed, of which TtOTS1 was found to play a negative effect in determining DT through affecting root development. This study sheds light on revealing the genetic basis underlying DT, which will provide the indispensable genes and germplasm resources for elite drought tolerance wheat improvement and breeding.
Collapse
Affiliation(s)
- Guang Yang
- State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China
- National Key Laboratory of Wheat Improvement, Peking University Institute of Advanced Agricultural Sciences, Shandong Laboratory of Advanced Agriculture Sciences in Weifang, Weifang, 261325, Shandong, China
| | - Yan Pan
- State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China
| | - Wenqiu Pan
- State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China
| | - Qingting Song
- State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China
| | - Ruoyu Zhang
- State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China
| | - Wei Tong
- State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China
| | - Licao Cui
- College of Biological Science and Engineering, Jiangxi Agricultural University, Nanchang, 330045, Jiangxi, China
| | - Wanquan Ji
- State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China
| | - Weining Song
- State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China
| | - Baoxing Song
- State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China
- National Key Laboratory of Wheat Improvement, Peking University Institute of Advanced Agricultural Sciences, Shandong Laboratory of Advanced Agriculture Sciences in Weifang, Weifang, 261325, Shandong, China
| | - Pingchuan Deng
- State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China
| | - Xiaojun Nie
- State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China
- Pioneering Innovation Center for Wheat Stress Tolerance Improvement, State Key Laboratory of Crop Stress Biology for Arid Areas, Yangling, 712100, Shaanxi, China
| |
Collapse
|
14
|
Abbas M, Diallo A, Goodney G, Gaye A. Leveraging the transcriptome to further our understanding of GWAS findings: eQTLs associated with genes related to LDL and LDL subclasses, in a cohort of African Americans. Front Genet 2024; 15:1345541. [PMID: 38384714 PMCID: PMC10879560 DOI: 10.3389/fgene.2024.1345541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Accepted: 01/16/2024] [Indexed: 02/23/2024] Open
Abstract
Background: GWAS discoveries often pose a significant challenge in terms of understanding their underlying mechanisms. Further research, such as an integration with expression quantitative trait locus (eQTL) analyses, are required to decipher the mechanisms connecting GWAS variants to phenotypes. An eQTL analysis was conducted on genes associated with low-density lipoprotein (LDL) cholesterol and its subclasses, with the aim of pinpointing genetic variants previously implicated in GWAS studies focused on lipid-related traits. Notably, the study cohort consisted of African Americans, a population characterized by a heightened prevalence of hypercholesterolemia. Methods: A comprehensive differential expression (DE) analysis was undertaken, with a dataset of 17,948 protein-coding mRNA transcripts extracted from the whole-blood transcriptomes of 416 samples to identify mRNA transcripts associated with LDL, with further granularity delineated between small LDL and large LDL subclasses. Subsequently, eQTL analysis was conducted with a subset of 242 samples for which whole-genome sequencing data were available to identify single-nucleotide polymorphisms (SNPs) associated with the LDL-related mRNA transcripts. Lastly, plausible functional connections were established between the identified eQTLs and genetic variants reported in the GWAS catalogue. Results: DE analysis revealed 1,048, 284, and 94 mRNA transcripts that exhibited differential expression in response to LDL, small LDL, and large LDL, respectively. The eQTL analysis identified a total of 9,950 significant SNP-mRNA associations involving 6,955 SNPs including a subset 101 SNPs previously documented in GWAS of LDL and LDL-related traits. Conclusion: Through comprehensive differential expression analysis, we identified numerous mRNA transcripts responsive to LDL, small LDL, and large LDL. Subsequent eQTL analysis revealed a rich landscape of eQTL-mRNA associations, including a subset of eQTL reported in GWAS studies of LDL and related traits. The study serves as a testament to the important role of integrative genomics in unraveling the enigmatic GWAS relationships between genetic variants and the complex fabric of human traits and diseases.
Collapse
Affiliation(s)
- Malak Abbas
- National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, United States
| | - Ana Diallo
- School of Nursing, Virginia Commonwealth University, Richmond, VA, United States
| | - Gabriel Goodney
- National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, United States
| | - Amadou Gaye
- National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, United States
| |
Collapse
|
15
|
Ban M, Bredikhin D, Huang Y, Bonder MJ, Katarzyna K, Oliver AJ, Wilson NK, Coupland P, Hadfield J, Göttgens B, Madissoon E, Stegle O, Sawcer S. Expression profiling of cerebrospinal fluid identifies dysregulated antiviral mechanisms in multiple sclerosis. Brain 2024; 147:554-565. [PMID: 38038362 PMCID: PMC10834244 DOI: 10.1093/brain/awad404] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2023] [Revised: 11/06/2023] [Accepted: 11/18/2023] [Indexed: 12/02/2023] Open
Abstract
Despite the overwhelming evidence that multiple sclerosis is an autoimmune disease, relatively little is known about the precise nature of the immune dysregulation underlying the development of the disease. Reasoning that the CSF from patients might be enriched for cells relevant in pathogenesis, we have completed a high-resolution single-cell analysis of 96 732 CSF cells collected from 33 patients with multiple sclerosis (n = 48 675) and 48 patients with other neurological diseases (n = 48 057). Completing comprehensive cell type annotation, we identified a rare population of CD8+ T cells, characterized by the upregulation of inhibitory receptors, increased in patients with multiple sclerosis. Applying a Multi-Omics Factor Analysis to these single-cell data further revealed that activity in pathways responsible for controlling inflammatory and type 1 interferon responses are altered in multiple sclerosis in both T cells and myeloid cells. We also undertook a systematic search for expression quantitative trait loci in the CSF cells. Of particular interest were two expression quantitative trait loci in CD8+ T cells that were fine mapped to multiple sclerosis susceptibility variants in the viral control genes ZC3HAV1 (rs10271373) and IFITM2 (rs1059091). Further analysis suggests that these associations likely reflect genetic effects on RNA splicing and cell-type specific gene expression respectively. Collectively, our study suggests that alterations in viral control mechanisms might be important in the development of multiple sclerosis.
Collapse
Affiliation(s)
- Maria Ban
- Department of Clinical Neurosciences, University of Cambridge, Cambridge CB2 0QQ, UK
| | - Danila Bredikhin
- European Molecular Biology Laboratory, Genome Biology Unit, 69117 Heidelberg, Germany
- Division of Computational Genomics and Systems Genetics, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Yuanhua Huang
- Department of Clinical Neurosciences, University of Cambridge, Cambridge CB2 0QQ, UK
- European Molecular Biology Laboratory, European Bioinformatics Institute, Cambridge CB10 1SD, UK
| | - Marc Jan Bonder
- European Molecular Biology Laboratory, Genome Biology Unit, 69117 Heidelberg, Germany
- Division of Computational Genomics and Systems Genetics, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Kania Katarzyna
- University of Cambridge, CRUK Cambridge Institute, Cambridge CB2 0RE, UK
| | - Amanda J Oliver
- Wellcome Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK
| | - Nicola K Wilson
- Department of Haematology, University of Cambridge, Cambridge CB2 0AW, UK
- Wellcome-Medical Research Council Cambridge Stem Cell Institute, University of Cambridge, Cambridge CB2 0AW, UK
| | - Paul Coupland
- University of Cambridge, CRUK Cambridge Institute, Cambridge CB2 0RE, UK
| | - James Hadfield
- University of Cambridge, CRUK Cambridge Institute, Cambridge CB2 0RE, UK
| | - Berthold Göttgens
- Department of Haematology, University of Cambridge, Cambridge CB2 0AW, UK
- Wellcome-Medical Research Council Cambridge Stem Cell Institute, University of Cambridge, Cambridge CB2 0AW, UK
| | - Elo Madissoon
- European Molecular Biology Laboratory, European Bioinformatics Institute, Cambridge CB10 1SD, UK
- Wellcome Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK
| | - Oliver Stegle
- European Molecular Biology Laboratory, Genome Biology Unit, 69117 Heidelberg, Germany
- Division of Computational Genomics and Systems Genetics, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
- European Molecular Biology Laboratory, European Bioinformatics Institute, Cambridge CB10 1SD, UK
| | - Stephen Sawcer
- Department of Clinical Neurosciences, University of Cambridge, Cambridge CB2 0QQ, UK
| |
Collapse
|
16
|
Boltz T, Schwarz T, Bot M, Hou K, Caggiano C, Lapinska S, Duan C, Boks MP, Kahn RS, Zaitlen N, Pasaniuc B, Ophoff R. Cell-type deconvolution of bulk-blood RNA-seq reveals biological insights into neuropsychiatric disorders. Am J Hum Genet 2024; 111:323-337. [PMID: 38306997 PMCID: PMC10870131 DOI: 10.1016/j.ajhg.2023.12.018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2023] [Revised: 12/12/2023] [Accepted: 12/13/2023] [Indexed: 02/04/2024] Open
Abstract
Genome-wide association studies (GWASs) have uncovered susceptibility loci associated with psychiatric disorders such as bipolar disorder (BP) and schizophrenia (SCZ). However, most of these loci are in non-coding regions of the genome, and the causal mechanisms of the link between genetic variation and disease risk is unknown. Expression quantitative trait locus (eQTL) analysis of bulk tissue is a common approach used for deciphering underlying mechanisms, although this can obscure cell-type-specific signals and thus mask trait-relevant mechanisms. Although single-cell sequencing can be prohibitively expensive in large cohorts, computationally inferred cell-type proportions and cell-type gene expression estimates have the potential to overcome these problems and advance mechanistic studies. Using bulk RNA-seq from 1,730 samples derived from whole blood in a cohort ascertained from individuals with BP and SCZ, this study estimated cell-type proportions and their relation with disease status and medication. For each cell type, we found between 2,875 and 4,629 eGenes (genes with an associated eQTL), including 1,211 that are not found on the basis of bulk expression alone. We performed a colocalization test between cell-type eQTLs and various traits and identified hundreds of associations that occur between cell-type eQTLs and GWASs but that are not detected in bulk eQTLs. Finally, we investigated the effects of lithium use on the regulation of cell-type expression loci and found examples of genes that are differentially regulated according to lithium use. Our study suggests that applying computational methods to large bulk RNA-seq datasets of non-brain tissue can identify disease-relevant, cell-type-specific biology of psychiatric disorders and psychiatric medication.
Collapse
Affiliation(s)
- Toni Boltz
- Department of Human Genetics, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, USA.
| | - Tommer Schwarz
- Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, CA, USA
| | - Merel Bot
- Center for Neurobehavioral Genetics, Semel Institute for Neuroscience and Human Behavior, University of California, Los Angeles, Los Angeles, CA, USA
| | - Kangcheng Hou
- Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, CA, USA
| | - Christa Caggiano
- Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, CA, USA
| | - Sandra Lapinska
- Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, CA, USA
| | - Chenda Duan
- Department of Computer Science, University of California, Los Angeles, Los Angeles, CA, USA
| | - Marco P Boks
- Department of Psychiatry, Brain Center, University Medical Center Utrecht, University Utrecht, Utrecht, the Netherlands
| | - Rene S Kahn
- Department of Psychiatry, Brain Center, University Medical Center Utrecht, University Utrecht, Utrecht, the Netherlands; Department of Psychiatry, Icahn School of Medicine, Mount Sinai, NY, USA
| | - Noah Zaitlen
- Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, CA, USA; Department of Neurology, University of California Los Angeles, Los Angeles, Los Angeles, CA, USA
| | - Bogdan Pasaniuc
- Department of Human Genetics, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, USA; Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, CA, USA; Department of Computational Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, USA; Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, USA
| | - Roel Ophoff
- Department of Human Genetics, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, USA; Bioinformatics Interdepartmental Program, University of California Los Angeles, Los Angeles, CA, USA; Center for Neurobehavioral Genetics, Semel Institute for Neuroscience and Human Behavior, University of California, Los Angeles, Los Angeles, CA, USA; Department of Psychiatry, Erasmus University Medical Center, Rotterdam, the Netherlands.
| |
Collapse
|
17
|
Diniz WJS, Afonso J, Kertz NC, Dyce PW, Banerjee P. Mapping Expression Quantitative Trait Loci Targeting Candidate Genes for Pregnancy in Beef Cows. Biomolecules 2024; 14:150. [PMID: 38397387 PMCID: PMC10886872 DOI: 10.3390/biom14020150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2023] [Revised: 01/23/2024] [Accepted: 01/24/2024] [Indexed: 02/25/2024] Open
Abstract
Despite collective efforts to understand the complex regulation of reproductive traits, no causative genes and/or mutations have been reported yet. By integrating genomics and transcriptomics data, potential regulatory mechanisms may be unveiled, providing opportunities to dissect the genetic factors governing fertility. Herein, we identified regulatory variants from RNA-Seq data associated with gene expression regulation in the uterine luminal epithelial cells of beef cows. We identified 4676 cis and 7682 trans eQTLs (expression quantitative trait loci) affecting the expression of 1120 and 2503 genes, respectively (FDR < 0.05). These variants affected the expression of transcription factor coding genes (71 cis and 193 trans eQTLs) and genes previously reported as differentially expressed between pregnant and nonpregnant cows. Functional over-representation analysis highlighted pathways related to metabolism, immune response, and hormone signaling (estrogen and GnRH) affected by eQTL-regulated genes (p-value ≤ 0.01). Furthermore, eQTLs were enriched in QTL regions for 13 reproduction-related traits from the CattleQTLdb (FDR ≤ 0.05). Our study provides novel insights into the genetic basis of reproductive processes in cattle. The underlying causal mechanisms modulating the expression of uterine genes warrant further investigation.
Collapse
Affiliation(s)
- Wellison J. S. Diniz
- Department of Animal Sciences, Auburn University, Auburn, AL 36849, USA; (N.C.K.); (P.W.D.); (P.B.)
| | - Juliana Afonso
- Embrapa Pecuária Sudeste, Rodovia Washington Luiz, Km 234, s/n, Fazenda Canchim, São Carlos 13560-970, SP, Brazil;
| | - Nicholas C. Kertz
- Department of Animal Sciences, Auburn University, Auburn, AL 36849, USA; (N.C.K.); (P.W.D.); (P.B.)
| | - Paul W. Dyce
- Department of Animal Sciences, Auburn University, Auburn, AL 36849, USA; (N.C.K.); (P.W.D.); (P.B.)
| | - Priyanka Banerjee
- Department of Animal Sciences, Auburn University, Auburn, AL 36849, USA; (N.C.K.); (P.W.D.); (P.B.)
| |
Collapse
|
18
|
Strober BJ, Tayeb K, Popp J, Qi G, Gordon MG, Perez R, Ye CJ, Battle A. SURGE: uncovering context-specific genetic-regulation of gene expression from single-cell RNA sequencing using latent-factor models. Genome Biol 2024; 25:28. [PMID: 38254214 PMCID: PMC10801966 DOI: 10.1186/s13059-023-03152-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Accepted: 12/20/2023] [Indexed: 01/24/2024] Open
Abstract
Genetic regulation of gene expression is a complex process, with genetic effects known to vary across cellular contexts such as cell types and environmental conditions. We developed SURGE, a method for unsupervised discovery of context-specific expression quantitative trait loci (eQTLs) from single-cell transcriptomic data. This allows discovery of the contexts or cell types modulating genetic regulation without prior knowledge. Applied to peripheral blood single-cell eQTL data, SURGE contexts capture continuous representations of distinct cell types and groupings of biologically related cell types. We demonstrate the disease-relevance of SURGE context-specific eQTLs using colocalization analysis and stratified LD-score regression.
Collapse
Affiliation(s)
- Benjamin J Strober
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Karl Tayeb
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
| | - Joshua Popp
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Guanghao Qi
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - M Grace Gordon
- Biological and Medical Informatics Graduate Program, University of California, San Francisco, CA, USA
- Division of Rheumatology, Department of Medicine, University of California, San Francisco, CA, USA
- Institute for Human Genetics, University of California, San Francisco, CA, USA
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
| | - Richard Perez
- Institute for Human Genetics, University of California, San Francisco, CA, USA
| | - Chun Jimmie Ye
- Division of Rheumatology, Department of Medicine, University of California, San Francisco, CA, USA
- Institute for Human Genetics, University of California, San Francisco, CA, USA
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
- Institute for Computational Health Sciences, University of California, San Francisco, San Francisco, CA, USA
- Chan-Zuckerberg Biohub, San Francisco, CA, USA
| | - Alexis Battle
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA.
- Department of Computer Science, Johns Hopkins University, Baltimore, MD, USA.
- Department of Genetic Medicine, Johns Hopkins University, Baltimore, MD, USA.
| |
Collapse
|
19
|
Venema WJ, Hiddingh S, van Loosdregt J, Bowes J, Balliu B, de Boer JH, Ossewaarde-van Norel J, Thompson SD, Langefeld CD, de Ligt A, van der Veken LT, Krijger PHL, de Laat W, Kuiper JJW. A cis-regulatory element regulates ERAP2 expression through autoimmune disease risk SNPs. Cell Genom 2024; 4:100460. [PMID: 38190099 PMCID: PMC10794781 DOI: 10.1016/j.xgen.2023.100460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/16/2023] [Revised: 10/04/2023] [Accepted: 11/09/2023] [Indexed: 01/09/2024]
Abstract
Single-nucleotide polymorphisms (SNPs) near the ERAP2 gene are associated with various autoimmune conditions, as well as protection against lethal infections. Due to high linkage disequilibrium, numerous trait-associated SNPs are correlated with ERAP2 expression; however, their functional mechanisms remain unidentified. We show by reciprocal allelic replacement that ERAP2 expression is directly controlled by the splice region variant rs2248374. However, disease-associated variants in the downstream LNPEP gene promoter are independently associated with ERAP2 expression. Allele-specific conformation capture assays revealed long-range chromatin contacts between the gene promoters of LNPEP and ERAP2 and showed that interactions were stronger in patients carrying the alleles that increase susceptibility to autoimmune diseases. Replacing the SNPs in the LNPEP promoter by reference sequences lowered ERAP2 expression. These findings show that multiple SNPs act in concert to regulate ERAP2 expression and that disease-associated variants can convert a gene promoter region into a potent enhancer of a distal gene.
Collapse
Affiliation(s)
- Wouter J Venema
- Department of Ophthalmology, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands; Center for Translational Immunology, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands
| | - Sanne Hiddingh
- Department of Ophthalmology, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands; Center for Translational Immunology, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands
| | - Jorg van Loosdregt
- Center for Translational Immunology, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands
| | - John Bowes
- Centre for Genetics and Genomics Versus Arthritis, Centre for Musculoskeletal Research, Manchester Academic Health Science Centre, The University of Manchester, Manchester, UK
| | - Brunilda Balliu
- Department of Computational Medicine, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA, USA
| | - Joke H de Boer
- Department of Ophthalmology, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands
| | | | - Susan D Thompson
- Department of Pediatrics, University of Cincinnati College of Medicine, Division of Human Genetics, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA
| | - Carl D Langefeld
- Department of Biostatistics and Data Science, and Center for Precision Medicine, Wake Forest University School of Medicine, Winston-Salem, NC, USA
| | - Aafke de Ligt
- Department of Ophthalmology, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands; Center for Translational Immunology, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands
| | - Lars T van der Veken
- Department of Genetics, Division Laboratories, Pharmacy and Biomedical Genetics, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands
| | - Peter H L Krijger
- Oncode Institute, Hubrecht Institute-KNAW and University Medical Center Utrecht, 3584 CT Utrecht, the Netherlands
| | - Wouter de Laat
- Oncode Institute, Hubrecht Institute-KNAW and University Medical Center Utrecht, 3584 CT Utrecht, the Netherlands
| | - Jonas J W Kuiper
- Department of Ophthalmology, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands; Center for Translational Immunology, University Medical Center Utrecht, Utrecht University, Utrecht, the Netherlands.
| |
Collapse
|
20
|
Hodonsky CJ, Turner AW, Khan MD, Barrientos NB, Methorst R, Ma L, Lopez NG, Mosquera JV, Auguste G, Farber E, Ma WF, Wong D, Onengut-Gumuscu S, Kavousi M, Peyser PA, van der Laan SW, Leeper NJ, Kovacic JC, Björkegren JLM, Miller CL. Multi-ancestry genetic analysis of gene regulation in coronary arteries prioritizes disease risk loci. Cell Genom 2024; 4:100465. [PMID: 38190101 PMCID: PMC10794848 DOI: 10.1016/j.xgen.2023.100465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Revised: 09/07/2023] [Accepted: 11/19/2023] [Indexed: 01/09/2024]
Abstract
Genome-wide association studies (GWASs) have identified hundreds of risk loci for coronary artery disease (CAD). However, non-European populations are underrepresented in GWASs, and the causal gene-regulatory mechanisms of these risk loci during atherosclerosis remain unclear. We incorporated local ancestry and haplotypes to identify quantitative trait loci for expression (eQTLs) and splicing (sQTLs) in coronary arteries from 138 ancestrally diverse Americans. Of 2,132 eQTL-associated genes (eGenes), 47% were previously unreported in coronary artery; 19% exhibited cell-type-specific expression. Colocalization revealed subgroups of eGenes unique to CAD and blood pressure GWAS. Fine-mapping highlighted additional eGenes, including TBX20 and IL5. We also identified sQTLs for 1,690 genes, among which TOR1AIP1 and ULK3 sQTLs demonstrated the importance of evaluating splicing to accurately identify disease-relevant isoform expression. Our work provides a patient-derived coronary artery eQTL resource and exemplifies the need for diverse study populations and multifaceted approaches to characterize gene regulation in disease processes.
Collapse
Affiliation(s)
- Chani J Hodonsky
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA
| | - Adam W Turner
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA
| | - Mohammad Daud Khan
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA
| | - Nelson B Barrientos
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA; Department of Genetic Medicine, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Ruben Methorst
- Central Diagnostics Laboratory, Division Laboratories, Pharmacy, and Biomedical Genetics, University Medical Center Utrecht, Utrecht University, 3584 CX Utrecht, the Netherlands
| | - Lijiang Ma
- Department of Genetics and Genomic Sciences, Icahn Institute for Genomics and Multiscale Biology, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Nicolas G Lopez
- Division of Vascular Surgery, Department of Surgery, Stanford University, Stanford, CA 94305, USA
| | - Jose Verdezoto Mosquera
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA; Department of Biochemistry and Molecular Genetics, University of Virginia, Charlottesville, VA 22908, USA
| | - Gaëlle Auguste
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA
| | - Emily Farber
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA
| | - Wei Feng Ma
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA; Medical Scientist Training Program, Department of Pathology, University of Virginia, Charlottesville, VA 22908, USA
| | - Doris Wong
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA; Department of Biochemistry and Molecular Genetics, University of Virginia, Charlottesville, VA 22908, USA
| | - Suna Onengut-Gumuscu
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA
| | - Maryam Kavousi
- Department of Epidemiology, Erasmus University Medical Center, 3000 CA Rotterdam, the Netherlands
| | - Patricia A Peyser
- Department of Epidemiology, University of Michigan, Ann Arbor, MI 48019, USA
| | - Sander W van der Laan
- Central Diagnostics Laboratory, Division Laboratories, Pharmacy, and Biomedical Genetics, University Medical Center Utrecht, Utrecht University, 3584 CX Utrecht, the Netherlands
| | - Nicholas J Leeper
- Division of Vascular Surgery, Department of Surgery, Stanford University, Stanford, CA 94305, USA
| | - Jason C Kovacic
- Cardiovascular Research Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA; Victor Chang Cardiac Research Institute, Darlinghurst, NSW 2010, Australia; St. Vincent's Clinical School, University of New South Wales, Sydney, NSW 2052, Australia
| | - Johan L M Björkegren
- Department of Genetics and Genomic Sciences, Icahn Institute for Genomics and Multiscale Biology, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA; Department of Medicine, Huddinge, Karolinska Institutet, 141 52 Huddinge, Sweden
| | - Clint L Miller
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA 22908, USA; Division of Vascular Surgery, Department of Surgery, Stanford University, Stanford, CA 94305, USA; Department of Public Health Sciences, University of Virginia, Charlottesville, VA 22908, USA.
| |
Collapse
|
21
|
Chen L, Liu L, Yang G, Li X, Dai X, Xue L, Yin T. Expression Quantitative Trait Locus of Wood Formation-Related Genes in Salix suchowensis. Int J Mol Sci 2023; 25:247. [PMID: 38203430 PMCID: PMC10778782 DOI: 10.3390/ijms25010247] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Revised: 12/20/2023] [Accepted: 12/21/2023] [Indexed: 01/12/2024] Open
Abstract
Shrub willows are widely planted for landscaping, soil remediation, and biomass production, due to their rapid growth rates. Identification of regulatory genes in wood formation would provide clues for genetic engineering of willows for improved growth traits on marginal lands. Here, we conducted an expression quantitative trait locus (eQTL) analysis, using a full sibling F1 population of Salix suchowensis, to explore the genetic mechanisms underlying wood formation. Based on variants identified from simplified genome sequencing and gene expression data from RNA sequencing, 16,487 eQTL blocks controlling 5505 genes were identified, including 2148 cis-eQTLs and 16,480 trans-eQTLs. eQTL hotspots were identified, based on eQTL frequency in genomic windows, revealing one hotspot controlling genes involved in wood formation regulation. Regulatory networks were further constructed, resulting in the identification of key regulatory genes, including three transcription factors (JAZ1, HAT22, MYB36) and CLV1, BAM1, CYCB2;4, CDKB2;1, associated with the proliferation and differentiation activity of cambium cells. The enrichment of genes in plant hormone pathways indicates their critical roles in the regulation of wood formation. Our analyses provide a significant groundwork for a comprehensive understanding of the regulatory network of wood formation in S. suchowensis.
Collapse
Affiliation(s)
| | | | | | | | | | - Liangjiao Xue
- State Key Laboratory of Tree Genetics and Breeding, Jiangsu Key Laboratory for Poplar Germplasm Enhancement and Variety Improvement, Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037, China
| | - Tongming Yin
- State Key Laboratory of Tree Genetics and Breeding, Jiangsu Key Laboratory for Poplar Germplasm Enhancement and Variety Improvement, Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037, China
| |
Collapse
|
22
|
Zeng Y, Jain R, Lam M, Ahmed M, Guo H, Xu W, Zhong Y, Wei GH, Xu W, He HH. DNA methylation modulated genetic variant effect on gene transcriptional regulation. Genome Biol 2023; 24:285. [PMID: 38066556 PMCID: PMC10709945 DOI: 10.1186/s13059-023-03130-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2023] [Accepted: 11/28/2023] [Indexed: 12/18/2023] Open
Abstract
BACKGROUND Expression quantitative trait locus (eQTL) analysis has emerged as an important tool in elucidating the link between genetic variants and gene expression, thereby bridging the gap between risk SNPs and associated diseases. We recently identified and validated a specific case where the methylation of a CpG site influences the relationship between the genetic variant and gene expression. RESULTS Here, to systematically evaluate this regulatory mechanism, we develop an extended eQTL mapping method, termed DNA methylation modulated eQTL (memo-eQTL). Applying this memo-eQTL mapping method to 128 normal prostate samples enables identification of 1063 memo-eQTLs, the majority of which are not recognized as conventional eQTLs in the same cohort. We observe that the methylation of the memo-eQTL CpG sites can either enhance or insulate the interaction between SNP and gene expression by altering CTCF-based chromatin 3D structure. CONCLUSIONS This study demonstrates the prevalence of memo-eQTLs paving the way to identify novel causal genes for traits or diseases associated with genetic variations.
Collapse
Affiliation(s)
- Yong Zeng
- Princess Margaret Cancer Centre, University Health Network, Toronto, Canada.
| | - Rahi Jain
- Princess Margaret Cancer Centre, University Health Network, Toronto, Canada
| | - Magnus Lam
- Princess Margaret Cancer Centre, University Health Network, Toronto, Canada
- Department of Medical Biophysics, University of Toronto, Toronto, Canada
| | - Musaddeque Ahmed
- Princess Margaret Cancer Centre, University Health Network, Toronto, Canada
| | - Haiyang Guo
- Department of Clinical Laboratory, the Second Hospital, Cheeloo College of Medicine, Shandong University, Jinan, 250033, Shandong, China
| | - Wenjie Xu
- MOE Key Laboratory of Metabolism and Molecular Medicine and Department of Biochemistry and Molecular Biology of School of Basic Medical Sciences, Fudan University Shanghai Cancer Center, Shanghai Medical College of Fudan University, Shanghai, China
| | - Yuan Zhong
- Princess Margaret Cancer Centre, University Health Network, Toronto, Canada
| | - Gong-Hong Wei
- MOE Key Laboratory of Metabolism and Molecular Medicine and Department of Biochemistry and Molecular Biology of School of Basic Medical Sciences, Fudan University Shanghai Cancer Center, Shanghai Medical College of Fudan University, Shanghai, China
- Biocenter Oulu & Faculty of Biochemistry and Molecular Medicine, University of Oulu, Oulu, Finland
| | - Wei Xu
- Dalla Lana School of Public Health, University of Toronto, Toronto, ON, Canada.
| | - Housheng Hansen He
- Princess Margaret Cancer Centre, University Health Network, Toronto, Canada.
- Department of Medical Biophysics, University of Toronto, Toronto, Canada.
| |
Collapse
|
23
|
Seifert F, Eisenblätter R, Beckmann J, Schürmann P, Hanel P, Jentschke M, Böhmer G, Strauß HG, Hirchenhain C, Schmidmayr M, Müller F, Fasching P, Luyten A, Häfner N, Dürst M, Runnebaum IB, Hillemanns P, Dörk T, Ramachandran D. Association of two genomic variants with HPV type-specific risk of cervical cancer. Tumour Virus Res 2023; 16:200269. [PMID: 37499979 PMCID: PMC10415783 DOI: 10.1016/j.tvr.2023.200269] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2023] [Revised: 07/22/2023] [Accepted: 07/24/2023] [Indexed: 07/29/2023] Open
Abstract
PROBLEM Human papillomavirus infection is integral to developing invasive cervical cancer in the majority of patients. In a recent genome-wide association study, rs9357152 and rs4243652 have been associated with seropositivity for HPV16 or HPV18, respectively. It is unknown whether these variants also associate with cervical cancer triggered by either HPV16 or HPV18. METHODS We investigate whether the two HPV susceptibility variants show association with type-specific cervical cancer in a genetic case-control study with cases stratified by HPV16 or HPV18, respectively. We further tested whether rs9357152 modulates gene expression of any of 36 genes at the human leukocyte antigen locus in 256 cervical tissues. RESULTS rs9357152 was associated with invasive HPV16-positive cervical cancer (OR 1.33, 95%CI 1.03-1.70, p = 0.03), and rs4243652 was associated with HPV18-positive adenocarcinomas (OR 2.96, 95%CI 1.18-7.41, p = 0.02). These associations remained borderline significant after testing against different sets of controls. rs9357152 was found to be an eQTL for HLA-DRB1 in HPV-positive cervical tissues (pANOVA = 0.0009), with the risk allele lowering mRNA levels. CONCLUSIONS We find evidence that HPV seropositivity variants at chromosome 6 and 14 may modulate type-specific cervical cancer risk. rs9357152 may exert its effect through regulating HLA-DRB1 induction in the presence of HPV. In regard of multiple testing, these results need to be confirmed in larger studies.
Collapse
Affiliation(s)
- Finja Seifert
- Gynaecology Research Unit, Comprehensive Cancer Center, Hannover Medical School, D-30625, Hannover, Germany
| | - Rieke Eisenblätter
- Gynaecology Research Unit, Comprehensive Cancer Center, Hannover Medical School, D-30625, Hannover, Germany
| | - Julia Beckmann
- Gynaecology Research Unit, Comprehensive Cancer Center, Hannover Medical School, D-30625, Hannover, Germany
| | - Peter Schürmann
- Gynaecology Research Unit, Comprehensive Cancer Center, Hannover Medical School, D-30625, Hannover, Germany
| | - Patricia Hanel
- Gynaecology Research Unit, Comprehensive Cancer Center, Hannover Medical School, D-30625, Hannover, Germany
| | - Matthias Jentschke
- Clinics of Gynaecology and Obstetrics, Hannover Medical School, D-30625, Hannover, Germany
| | | | - Hans-Georg Strauß
- Department of Gynaecology, University Clinics, Martin-Luther University, Halle-Wittenberg, Germany
| | - Christine Hirchenhain
- Department of Gynaecology, Clinics Carl Gustav Carus, University of Dresden, Dresden, Germany
| | - Monika Schmidmayr
- Department of Gynaecology, Technische Universität München, Munich, Germany
| | - Florian Müller
- Martin-Luther Hospital, Charite University, Berlin, Germany
| | - Peter Fasching
- Department of Gynaecology and Obstetrics, Erlangen University Hospital, Friedrich-Alexander University of Erlangen-Nuremberg, Erlangen, Germany
| | - Alexander Luyten
- Dysplasia Unit, Department of Gynaecology and Obstetrics, Mare Klinikum, Kronshagen, Germany; Department of Gynaecology, Wolfsburg Hospital, Wolfsburg, Germany
| | - Norman Häfner
- Department of Gynaecology, Jena University Hospital, Friedrich -Schiller-University Jena, Jena, Germany
| | - Matthias Dürst
- Department of Gynaecology, Jena University Hospital, Friedrich -Schiller-University Jena, Jena, Germany
| | - Ingo B Runnebaum
- Department of Gynaecology, Jena University Hospital, Friedrich -Schiller-University Jena, Jena, Germany
| | - Peter Hillemanns
- Clinics of Gynaecology and Obstetrics, Hannover Medical School, D-30625, Hannover, Germany
| | - Thilo Dörk
- Gynaecology Research Unit, Comprehensive Cancer Center, Hannover Medical School, D-30625, Hannover, Germany
| | - Dhanya Ramachandran
- Gynaecology Research Unit, Comprehensive Cancer Center, Hannover Medical School, D-30625, Hannover, Germany.
| |
Collapse
|
24
|
Zhang J, Zhao H. eQTL studies: from bulk tissues to single cells. J Genet Genomics 2023; 50:925-933. [PMID: 37207929 PMCID: PMC10656365 DOI: 10.1016/j.jgg.2023.05.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2023] [Revised: 05/02/2023] [Accepted: 05/04/2023] [Indexed: 05/21/2023]
Abstract
An expression quantitative trait locus (eQTL) is a chromosomal region where genetic variants are associated with the expression levels of specific genes that can be both nearby or distant. The identifications of eQTLs for different tissues, cell types, and contexts have led to a better understanding of the dynamic regulations of gene expressions and implications of functional genes and variants for complex traits and diseases. Although most eQTL studies have been performed on data collected from bulk tissues, recent studies have demonstrated the importance of cell-type-specific and context-dependent gene regulations in biological processes and disease mechanisms. In this review, we discuss statistical methods that have been developed to enable the detection of cell-type-specific and context-dependent eQTLs from bulk tissues, purified cell types, and single cells. We also discuss the limitations of the current methods and future research opportunities.
Collapse
Affiliation(s)
- Jingfei Zhang
- Information Systems and Operations Management, Emory University, Atlanta, GA 30322, USA
| | - Hongyu Zhao
- Department of Biostatistics, Yale School of Public Health, New Haven, CT 208034, USA.
| |
Collapse
|
25
|
Kulski JK, Pfaff AL, Marney LD, Fröhlich A, Bubb VJ, Quinn JP, Koks S. Regulation of expression quantitative trait loci by SVA retrotransposons within the major histocompatibility complex. Exp Biol Med (Maywood) 2023; 248:2304-2318. [PMID: 38031415 PMCID: PMC10903234 DOI: 10.1177/15353702231209411] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2023] [Accepted: 10/04/2023] [Indexed: 12/01/2023] Open
Abstract
Genomic and transcriptomic studies of expression quantitative trait loci (eQTL) revealed that SINE-VNTR-Alu (SVA) retrotransposon insertion polymorphisms (RIPs) within human genomes markedly affect the co-expression of many coding and noncoding genes by coordinated regulatory processes. This study examined the polymorphic SVA modulation of gene co-expression within the major histocompatibility complex (MHC) genomic region where more than 160 coding genes are involved in innate and adaptive immunity. We characterized the modulation of SVA RIPs utilizing the genomic and transcriptomic sequencing data obtained from whole blood of 1266 individuals in the Parkinson's Progression Markers Initiative (PPMI) cohort that included an analysis of human leukocyte antigen (HLA)-A regulation in a subpopulation of the cohort. The regulatory properties of eight SVAs located within the class I and class II MHC regions were associated with differential co-expression of 71 different genes within and 75 genes outside the MHC region. Some of the same genes were affected by two or more different SVA. Five SVA are annotated in the human genomic reference sequence GRCh38.p14/hg38, whereas the other three were novel insertions within individuals. We also examined and found distinct structural effects (long and short variants and the CT internal variants) for one of the SVA (R_SVA_24) insertions on the differential expression of the HLA-A gene within a subpopulation (550 individuals) of the PPMI cohort. This is the first time that many HLA and non-HLA genes (multilocus expression units) and splicing mechanisms have been shown to be regulated by eight structurally polymorphic SVA within the MHC genomic region by applying precise statistical analysis of RNA data derived from the blood samples of a human cohort population. This study shows that SVA within the MHC region are important regulators or rheostats of gene co-expression that might have potential roles in diversity, health, and disease.
Collapse
Affiliation(s)
- Jerzy K Kulski
- Department of Molecular Life Sciences, School of Medicine, Tokai University, Isehara, Kanagawa 259–1193, Japan
- Health and Medical Science. Division of Immunology and Microbiology, School of Biomedical Sciences, The University of Western Australia, Nedlands, WA 6009, Australia
| | - Abigail L Pfaff
- Perron Institute for Neurological and Translational Science, Perth, WA 6009, Australia
- Centre for Molecular Medicine and Innovative Therapeutics, Murdoch University, Perth, WA 6150, Australia
| | - Luke D Marney
- Department of Pharmacology and Therapeutics, Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, L69 3BX, UK
| | - Alexander Fröhlich
- Department of Pharmacology and Therapeutics, Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, L69 3BX, UK
| | - Vivien J Bubb
- Department of Pharmacology and Therapeutics, Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, L69 3BX, UK
| | - John P Quinn
- Department of Pharmacology and Therapeutics, Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, L69 3BX, UK
| | - Sulev Koks
- Perron Institute for Neurological and Translational Science, Perth, WA 6009, Australia
- Centre for Molecular Medicine and Innovative Therapeutics, Murdoch University, Perth, WA 6150, Australia
| |
Collapse
|
26
|
Wang J, Cheng X, Liang Q, Owen LA, Lu J, Zheng Y, Wang M, Chen S, DeAngelis MM, Li Y, Chen R. Single-cell multiomics of the human retina reveals hierarchical transcription factor collaboration in mediating cell type-specific effects of genetic variants on gene regulation. Genome Biol 2023; 24:269. [PMID: 38012720 PMCID: PMC10680294 DOI: 10.1186/s13059-023-03111-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2022] [Accepted: 11/15/2023] [Indexed: 11/29/2023] Open
Abstract
BACKGROUND Systematic characterization of how genetic variation modulates gene regulation in a cell type-specific context is essential for understanding complex traits. To address this question, we profile gene expression and chromatin accessibility in cells from healthy retinae of 20 human donors through single-cell multiomics and genomic sequencing. RESULTS We map eQTL, caQTL, allelic-specific expression, and allelic-specific chromatin accessibility in major retinal cell types. By integrating these results, we identify and characterize regulatory elements and genetic variants effective on gene regulation in individual cell types. The majority of identified sc-eQTLs and sc-caQTLs display cell type-specific effects, while the cis-elements containing genetic variants with cell type-specific effects are often accessible in multiple cell types. Furthermore, the transcription factors whose binding sites are perturbed by genetic variants tend to have higher expression levels in the cell types where the variants exert their effects, compared to the cell types where the variants have no impact. We further validate our findings with high-throughput reporter assays. Lastly, we identify the enriched cell types, candidate causal variants and genes, and cell type-specific regulatory mechanism underlying GWAS loci. CONCLUSIONS Overall, genetic effects on gene regulation are highly context dependent. Our results suggest that cell type-dependent genetic effect is driven by precise modulation of both trans-factor expression and chromatin accessibility of cis-elements. Our findings indicate hierarchical collaboration among transcription factors plays a crucial role in mediating cell type-specific effects of genetic variants on gene regulation.
Collapse
Affiliation(s)
- Jun Wang
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, USA
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
| | - Xuesen Cheng
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, USA
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
| | - Qingnan Liang
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, USA
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
- Verna and Marrs McLean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, TX, USA
| | - Leah A Owen
- Department of Ophthalmology and Visual Sciences, John A. Moran Eye Center, University of Utah, Salt Lake City, UT, USA
| | - Jiaxiong Lu
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, USA
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
- Verna and Marrs McLean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, TX, USA
| | - Yiqiao Zheng
- Department of Ophthalmology and Visual Sciences, Washington University in St Louis, Saint Louis, MO, USA
| | - Meng Wang
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, USA
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
| | - Shiming Chen
- Department of Ophthalmology and Visual Sciences, Washington University in St Louis, Saint Louis, MO, USA
- Department of Developmental Biology, Washington University in St Louis, Saint Louis, MO, USA
| | - Margaret M DeAngelis
- Department of Ophthalmology, University at Buffalo the State University of New York, Buffalo, NY, USA
| | - Yumei Li
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, USA
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
| | - Rui Chen
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, USA.
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA.
| |
Collapse
|
27
|
Ding R, Zou X, Qin Y, Gong L, Chen H, Ma X, Guang S, Yu C, Wang G, Li L. xQTLbiolinks: a comprehensive and scalable tool for integrative analysis of molecular QTLs. Brief Bioinform 2023; 25:bbad440. [PMID: 38058186 PMCID: PMC10701093 DOI: 10.1093/bib/bbad440] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Revised: 10/23/2023] [Accepted: 11/11/2023] [Indexed: 12/08/2023] Open
Abstract
Genome-wide association studies (GWAS) have identified thousands of disease-associated non-coding variants, posing urgent needs for functional interpretation. Molecular Quantitative Trait Loci (xQTLs) such as eQTLs serve as an essential intermediate link between these non-coding variants and disease phenotypes and have been widely used to discover disease-risk genes from many population-scale studies. However, mining and analyzing the xQTLs data presents several significant bioinformatics challenges, particularly when it comes to integration with GWAS data. Here, we developed xQTLbiolinks as the first comprehensive and scalable tool for bulk and single-cell xQTLs data retrieval, quality control and pre-processing from public repositories and our integrated resource. In addition, xQTLbiolinks provided a robust colocalization module through integration with GWAS summary statistics. The result generated by xQTLbiolinks can be flexibly visualized or stored in standard R objects that can easily be integrated with other R packages and custom pipelines. We applied xQTLbiolinks to cancer GWAS summary statistics as case studies and demonstrated its robust utility and reproducibility. xQTLbiolinks will profoundly accelerate the interpretation of disease-associated variants, thus promoting a better understanding of disease etiologies. xQTLbiolinks is available at https://github.com/lilab-bioinfo/xQTLbiolinks.
Collapse
Affiliation(s)
- Ruofan Ding
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Xudong Zou
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Yangmei Qin
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Lihai Gong
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Hui Chen
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Xuelian Ma
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Shouhong Guang
- Department of Obstetrics and Gynecology, The First Affiliated Hospital of USTC, The USTC RNA Institute, Ministry of Education Key Laboratory for Membraneless Organelles & Cellular Dynamics, School of Life Sciences, Division of Life Sciences and Medicine, Biomedical Sciences and Health Laboratory of Anhui Province, University of Science and Technology of China, Hefei, Anhui 230027, China
| | - Chen Yu
- Institute of Cancer Research, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Gao Wang
- The Gertrude H. Sergievsky Center and the Department of Neurology, Columbia University, NY 10032, USA
| | - Lei Li
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| |
Collapse
|
28
|
Zhao J, Wang R, Song L, Han H, Wang P, Zhao Y, Zhang Y, Zhang H. Causal association between lipid-lowering drugs and female reproductive endocrine diseases: a drug-targeted Mendelian randomization study. Front Endocrinol (Lausanne) 2023; 14:1295412. [PMID: 38027179 PMCID: PMC10668027 DOI: 10.3389/fendo.2023.1295412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/16/2023] [Accepted: 10/30/2023] [Indexed: 12/01/2023] Open
Abstract
Purpose The relationship between dyslipidemia and female reproductive endocrine diseases has been increasingly studied. The use of lipid-lowering drugs in treating various related diseases, including coronary heart disease, may affect female reproductive endocrine diseases. Therefore, our study aims to investigate the effects of lipid-lowering drugs on female reproductive endocrine diseases and provide a basis for the appropriate selection of drugs. Methods In this study, we focused on three drug targets of statins, namely HMG-CoA reductase (HMGCR) inhibitors, proprotein convertase kexin 9 (PCSK9) inhibitors, and Niemann-Pick C1-Like 1 (NPC1L1) inhibitors. To identify potential inhibitors for these targets, we collected single nucleotide polymorphisms (SNPs) associated with HMGCR, PCSK9, and NPC1L1 from published genome-wide association study statistics. Subsequently, we conducted a drug target Mendelian randomization (MR) analysis to investigate the effects of these inhibitors on reproductive endocrine diseases mediated by low-density lipoprotein cholesterol (LDL-C) levels. Alongside coronary heart disease as a positive control, our main outcomes of interest included the risk of polycystic ovary syndrome (PCOS), premature ovarian insufficiency (POI), premenstrual syndrome (PMS), abnormal uterine bleeding (including menorrhagia and oligomenorrhea), and infertility. Results PCSK9 inhibitors significantly increased the risk of infertility in patients (OR [95%CI] = 1.14 [1.06, 1.23], p<0.05). In contrast, HMGCR inhibitors significantly reduced the risk of menorrhagia in female patients (OR [95%CI] = 0.85 [0.75, 0.97], p<0.05), but had no statistical impact on patients with oligomenorrhea. Conclusion The findings suggest that PCSK9 inhibitors may significantly increase the risk of infertility in patients. On the other hand, HMGCR inhibitors could potentially offer protection against menorrhagia in women. However, no effects of lipid-lowering drugs have been observed on other reproductive endocrine disorders, such as PCOS, POF, PMS and oligomenorrhea.
Collapse
Affiliation(s)
- Jing Zhao
- Department of Gynecology, Hebei General Hospital, Shijiazhuang, China
| | - Runfang Wang
- Department of Obstetrics, Hebei General Hospital, Shijiazhuang, China
| | - Liyun Song
- Department of Gynecology, Hebei General Hospital, Shijiazhuang, China
| | - Hua Han
- Department of Gynecology, Hebei General Hospital, Shijiazhuang, China
| | - Pei Wang
- Department of Gynecology, Hebei General Hospital, Shijiazhuang, China
| | - Yuan Zhao
- Department of Clinical Laboratories, Kunhua Affiliated Hospital, Kunming University of Science and Technology, Kunming, China
| | - Yunxia Zhang
- Department of Gynecology, Hebei General Hospital, Shijiazhuang, China
| | - Hongzhen Zhang
- Department of Obstetrics and Gynecology, The First Hospital of Hebei Medical University, Shijiazhuang, China
| |
Collapse
|
29
|
Million CR, Wijeratne S, Karhoff S, Cassone BJ, McHale LK, Dorrance AE. Molecular mechanisms underpinning quantitative resistance to Phytophthora sojae in Glycine max using a systems genomics approach. Front Plant Sci 2023; 14:1277585. [PMID: 38023885 PMCID: PMC10662313 DOI: 10.3389/fpls.2023.1277585] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Accepted: 10/16/2023] [Indexed: 12/01/2023]
Abstract
Expression of quantitative disease resistance in many host-pathogen systems is controlled by genes at multiple loci, each contributing a small effect to the overall response. We used a systems genomics approach to study the molecular underpinnings of quantitative disease resistance in the soybean-Phytophthora sojae pathosystem, incorporating expression quantitative trait loci (eQTL) mapping and gene co-expression network analysis to identify the genes putatively regulating transcriptional changes in response to inoculation. These findings were compared to previously mapped phenotypic (phQTL) to identify the molecular mechanisms contributing to the expression of this resistance. A subset of 93 recombinant inbred lines (RILs) from a Conrad × Sloan population were inoculated with P. sojae isolate 1.S.1.1 using the tray-test method; RNA was extracted, sequenced, and the normalized read counts were genetically mapped from tissue collected at the inoculation site 24 h after inoculation from both mock and inoculated samples. In total, more than 100,000 eQTLs were mapped. There was a switch from predominantly cis-eQTLs in the mock treatment to an almost entirely nonoverlapping set of predominantly trans-eQTLs in the inoculated treatment, where greater than 100-fold more eQTLs were mapped relative to mock, indicating vast transcriptional reprogramming due to P. sojae infection occurred. The eQTLs were organized into 36 hotspots, with the four largest hotspots from the inoculated treatment corresponding to more than 70% of the eQTLs, each enriched for genes within plant-pathogen interaction pathways. Genetic regulation of trans-eQTLs in response to the pathogen was predicted to occur through transcription factors and signaling molecules involved in plant-pathogen interactions, plant hormone signal transduction, and MAPK pathways. Network analysis identified three co-expression modules that were correlated with susceptibility to P. sojae and associated with three eQTL hotspots. Among the eQTLs co-localized with phQTLs, two cis-eQTLs with putative functions in the regulation of root architecture or jasmonic acid, as well as the putative master regulators of an eQTL hotspot nearby a phQTL, represent candidates potentially underpinning the molecular control of these phQTLs for resistance.
Collapse
Affiliation(s)
- Cassidy R. Million
- Department of Plant Pathology, The Ohio State University, Wooster, OH, United States
- Center for Soybean Research and Center for Applied Plant Sciences, The Ohio State University, Columbus, OH, United States
| | - Saranga Wijeratne
- Molecular and Cellular Imaging Center, The Ohio State University, Wooster, OH, United States
| | - Stephanie Karhoff
- Center for Soybean Research and Center for Applied Plant Sciences, The Ohio State University, Columbus, OH, United States
- Translational Plant Sciences Graduate Program, The Ohio State University, Columbus, OH, United States
| | - Bryan J. Cassone
- Center for Soybean Research and Center for Applied Plant Sciences, The Ohio State University, Columbus, OH, United States
- Department of Biology, Brandon University, Brandon, Manitoba, MB, Canada
| | - Leah K. McHale
- Center for Soybean Research and Center for Applied Plant Sciences, The Ohio State University, Columbus, OH, United States
- Department of Horticulture and Crop Science, The Ohio State University, Columbus, OH, United States
| | - Anne E. Dorrance
- Department of Plant Pathology, The Ohio State University, Wooster, OH, United States
- Center for Soybean Research and Center for Applied Plant Sciences, The Ohio State University, Columbus, OH, United States
| |
Collapse
|
30
|
Lin C, Liu W, Jiang W, Zhao H. Robustness of quantifying mediating effects of genetically regulated expression on complex traits with mediated expression score regression. Biol Methods Protoc 2023; 8:bpad024. [PMID: 37901453 PMCID: PMC10599978 DOI: 10.1093/biomethods/bpad024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Revised: 10/08/2023] [Accepted: 10/16/2023] [Indexed: 10/31/2023] Open
Abstract
Genetic association signals have been mostly found in noncoding regions through genome-wide association studies (GWAS), suggesting the roles of gene expression regulation in human diseases and traits. However, there has been limited success in colocalizing expression quantitative trait locus (eQTL) with disease-associated variants. Mediated expression score regression (MESC) is a recently proposed method to quantify the proportion of trait heritability mediated by genetically regulated gene expressions (GReX). Applications of MESC to GWAS results have yielded low estimation of mediated heritability for many traits. As MESC relies on stringent independence assumptions between cis-eQTL effects, gene effects, and nonmediated SNP effects, it may fail to characterize the true relationships between those effect sizes, which leads to biased results. Here, we consider the robustness of MESC to investigate whether the low fraction of mediated heritability inferred by MESC reflects biological reality for complex traits or is an underestimation caused by model misspecifications. Our results suggest that MESC may lead to biased estimates of mediated heritability with misspecification of gene annotations leading to underestimation, whereas misspecification of SNP annotations may lead to overestimation. Furthermore, errors in eQTL effect estimates may lead to underestimation of mediated heritability.
Collapse
Affiliation(s)
- Chen Lin
- Department of Biostatistics, Yale University, New Haven, CT 06510, United States
| | - Wei Liu
- Program of Computational Biology and Bioinformatics, Yale University, New Haven, CT 06510, United States
| | - Wei Jiang
- Department of Biostatistics, Yale University, New Haven, CT 06510, United States
| | - Hongyu Zhao
- Department of Biostatistics, Yale University, New Haven, CT 06510, United States
- Program of Computational Biology and Bioinformatics, Yale University, New Haven, CT 06510, United States
| |
Collapse
|
31
|
Gong R, Qiu M, Cao J, Zhou Z, Wei Y, Wen Q, Lin Q, Wei X, Liang X, Jiang Y, Chen P, Wei J, Zhan S, Liu Y, Yu H. Potentially Functional Genetic Variants in the NRF2 Signaling Pathway Genes are Associated With HBV-related Hepatocellular Carcinoma Survival. J Cancer 2023; 14:3387-3396. [PMID: 38021150 PMCID: PMC10647191 DOI: 10.7150/jca.88561] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Accepted: 09/17/2023] [Indexed: 12/01/2023] Open
Abstract
The nuclear factor E2-related factor 2 (NRF2) signaling pathway is one of the most important cell defense pathways. However, it is unclear whether genetic variants in NRF2 signaling pathway genes are associated with the survival of hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC). In the present study, we utilized a new hypothesis-driven approach based on biological pathways to investigate the associations between 17919 single nucleotide polymorphisms (SNPs) in 137 NRF2 signaling pathway genes and the overall survival (OS) of 866 patients with HBV-related HCC. As a result, two independent SNPs with potential biological function were identified to be significantly associated with HBV-related HCC OS: [SLC2A9 rs28643326 T>C: hazard ratio (HR) = 0.74, 95% confidence interval (95% CI) = 0.62-0.89, P < 0.001 and SLC5A10 rs2472711 G>T: HR = 0.81, 95% CI = 0.71-0.93, P = 0.003, respectively]. The expression quantitative trait loci (eQTL) analysis further revealed that the rs28643326 C allele was significantly associated with increased levels of SLC2A9 mRNA expression (P < 0.001), and higher mRNA expression levels of SLC2A9 in adjacent normal liver tissues were associated with better survival. Although the association between the rs2472711 T allele and the mRNA expression of SLC5A10 was not statistically significant (P = 0.200), the fact that rs2472711 is located at the DNase I hypersensitivity site and is a marker for promoter and enhancer histones also suggests that it may have the function of regulating its corresponding gene expression. In conclusion, genetic variants of NRF2 signaling pathway genes may serve as potential prognostic biomarkers for HBV-related HCC and also provide a solid basis for further mechanistic exploration.
Collapse
Affiliation(s)
- Rongbin Gong
- Department of Experimental Research, Guangxi Medical University Cancer Hospital, Nanning 530000, China
- Department of Epidemiology and Health Statistics, School of Public Health, Guangxi Medical University, Nanning 530000, China
| | - Moqin Qiu
- Department of Respiratory Oncology, Guangxi Medical University Cancer Hospital, Nanning 530000, China
| | - Ji Cao
- Department of Cancer Prevention and Control, Guangxi Medical University Cancer Hospital, Nanning 530000, China
| | - Zihan Zhou
- Department of Cancer Prevention and Control, Guangxi Medical University Cancer Hospital, Nanning 530000, China
| | - Yuying Wei
- Department of Experimental Research, Guangxi Medical University Cancer Hospital, Nanning 530000, China
| | - Qiuping Wen
- Department of Experimental Research, Guangxi Medical University Cancer Hospital, Nanning 530000, China
| | - Qiuling Lin
- Department of Clinical Research, Guangxi Medical University Cancer Hospital, Nanning 530000, China
| | - Xiaoxia Wei
- Department of Clinical Research, Guangxi Medical University Cancer Hospital, Nanning 530000, China
| | - Xiumei Liang
- Department of Disease Process Management, Guangxi Medical University Cancer Hospital, Nanning 530000, China
| | - Yanji Jiang
- Department of Scientific Research Dept, Guangxi Medical University Cancer Hospital, Nanning 530000, China
| | - Peiqin Chen
- Department of Experimental Research, Guangxi Medical University Cancer Hospital, Nanning 530000, China
| | - Junjie Wei
- Department of Epidemiology and Health Statistics, School of Public Health, Guangxi Medical University, Nanning 530000, China
| | - Shicheng Zhan
- Department of Epidemiology and Health Statistics, School of Public Health, Guangxi Medical University, Nanning 530000, China
| | - Yingchun Liu
- Department of Experimental Research, Guangxi Medical University Cancer Hospital, Nanning 530000, China
- Key Cultivated Laboratory of Cancer Molecular Medicine of Guangxi Health Commission, Guangxi Medical University Cancer Hospital, Nanning 530000, China
| | - Hongping Yu
- Department of Experimental Research, Guangxi Medical University Cancer Hospital, Nanning 530000, China
- Key Cultivated Laboratory of Cancer Molecular Medicine of Guangxi Health Commission, Guangxi Medical University Cancer Hospital, Nanning 530000, China
- Key Laboratory of Early Prevention and Treatment for Regional High Frequency Tumor (Guangxi Medical University), Ministry of Education, Nanning 530000, China
| |
Collapse
|
32
|
Badia-Bringué G, Canive M, Fernandez-Jimenez N, Lavín JL, Casais R, Blanco-Vázquez C, Vázquez P, Fernández A, Bilbao JR, Garrido JM, Juste RA, González-Recio O, Alonso-Hearn M. Summary-data based Mendelian randomization identifies gene expression regulatory polymorphisms associated with bovine paratuberculosis by modulation of the nuclear factor Kappa β (NF-κß)-mediated inflammatory response. BMC Genomics 2023; 24:605. [PMID: 37821814 PMCID: PMC10568764 DOI: 10.1186/s12864-023-09710-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2023] [Accepted: 10/02/2023] [Indexed: 10/13/2023] Open
Abstract
Genome-wide association studies (GWAS) have identified host genetic variants associated with paratuberculosis (PTB) susceptibility. Most of the GWAS-identified SNPs are in non-coding regions. Connecting these non-coding variants and downstream affected genes is a challenge and, up to date, only a few functional mutations or expression quantitative loci (cis-eQTLs) associated with PTB susceptibility have been identified. In the current study, the associations between imputed whole-genome sequence genotypes and whole RNA-Sequencing data from peripheral blood (PB) and ileocecal valve (ICV) samples of Spanish Holstein cows (N = 16) were analyzed with TensorQTL. This approach allowed the identification of 88 and 37 cis-eQTLs regulating the expression levels of 90 and 37 genes in PB and ICV samples, respectively (False discorey rate, FDR ≤ 0.05). Next, we applied summary-based data Mendelian randomization (SMR) to integrate the cis-eQTL dataset with GWAS data obtained from a cohort of 813 culled cattle that were classified according to the presence or absence of PTB-associated histopathological lesions in gut tissues. After multiple testing corrections (FDR ≤ 0.05), we identified two novel cis-eQTLs affecting the expression of the early growth response factor 4 (EGR4) and the bovine neuroblastoma breakpoint family member 6-like protein isoform 2 (MGC134040) that showed pleiotropic associations with the presence of multifocal and diffuse lesions in gut tissues; P = 0.002 and P = 0.017, respectively. While EGR4 acts as a brake on T-cell proliferation and cytokine production through interaction with the nuclear factor Kappa β (NF-κß), MGC134040 is a target gene of NF-κß. Our findings provide a better understanding of the genetic factors influencing PTB outcomes, confirm that the multifocal lesions are localized/confined lesions that have different underlying host genetics than the diffuse lesions, and highlight regulatory SNPs and regulated-gene targets to design future functional studies.
Collapse
Affiliation(s)
- Gerard Badia-Bringué
- Department of Animal Health, NEIKER- Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain
- Doctoral Program in Molecular Biology and Biomedicine, Universidad del País Vasco/Euskal Herriko Unibertsitatea (UPV/EHU), Leioa, Bizkaia, Spain
| | - Maria Canive
- Department of Animal Health, NEIKER- Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain
| | - Nora Fernandez-Jimenez
- Department of Genetics, Physical Anthropology and Animal Physiology, University of the Basque Country (UPV/EHU), Biocruces-Bizkaia HRI, Leioa, Bizkaia, Spain
| | - José Luis Lavín
- Department of Applied Mathematics, NEIKER- Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain
| | - Rosa Casais
- Center of Animal Biotechnology, SERIDA, Servicio Regional de Investigación y Desarrollo Agroalimentario, Deva, Asturias, Spain
| | - Cristina Blanco-Vázquez
- Center of Animal Biotechnology, SERIDA, Servicio Regional de Investigación y Desarrollo Agroalimentario, Deva, Asturias, Spain
| | - Patricia Vázquez
- Department of Animal Health, NEIKER- Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain
| | - Almudena Fernández
- Departamento de Mejora Genética Animal, Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria, CSIC, Madrid, Spain
| | - Jose Ramón Bilbao
- Department of Genetics, Physical Anthropology and Animal Physiology, University of the Basque Country (UPV/EHU), Biocruces-Bizkaia HRI, Leioa, Bizkaia, Spain
| | - Joseba M Garrido
- Department of Animal Health, NEIKER- Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain
| | - Ramón A Juste
- Department of Animal Health, NEIKER- Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain
| | - Oscar González-Recio
- Departamento de Mejora Genética Animal, Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria, CSIC, Madrid, Spain
| | - Marta Alonso-Hearn
- Department of Animal Health, NEIKER- Basque Institute for Agricultural Research and Development, Basque Research and Technology Alliance (BRTA), Derio, Bizkaia, Spain.
| |
Collapse
|
33
|
McAfee JC, Lee S, Lee J, Bell JL, Krupa O, Davis J, Insigne K, Bond ML, Zhao N, Boyle AP, Phanstiel DH, Love MI, Stein JL, Ruzicka WB, Davila-Velderrain J, Kosuri S, Won H. Systematic investigation of allelic regulatory activity of schizophrenia-associated common variants. Cell Genomics 2023; 3:100404. [PMID: 37868037 PMCID: PMC10589626 DOI: 10.1016/j.xgen.2023.100404] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/22/2022] [Revised: 02/23/2023] [Accepted: 08/21/2023] [Indexed: 10/24/2023]
Abstract
Genome-wide association studies (GWASs) have successfully identified 145 genomic regions that contribute to schizophrenia risk, but linkage disequilibrium makes it challenging to discern causal variants. We performed a massively parallel reporter assay (MPRA) on 5,173 fine-mapped schizophrenia GWAS variants in primary human neural progenitors and identified 439 variants with allelic regulatory effects (MPRA-positive variants). Transcription factor binding had modest predictive power, while fine-map posterior probability, enhancer overlap, and evolutionary conservation failed to predict MPRA-positive variants. Furthermore, 64% of MPRA-positive variants did not exhibit expressive quantitative trait loci signature, suggesting that MPRA could identify yet unexplored variants with regulatory potentials. To predict the combinatorial effect of MPRA-positive variants on gene regulation, we propose an accessibility-by-contact model that combines MPRA-measured allelic activity with neuronal chromatin architecture.
Collapse
Affiliation(s)
- Jessica C McAfee
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Neuroscience Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Curriculum in Genetics and Molecular Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Sool Lee
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Neuroscience Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Jiseok Lee
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Neuroscience Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Jessica L Bell
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Neuroscience Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Oleh Krupa
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Neuroscience Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Jessica Davis
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, CA 90095, USA
- UCLA-DOE Institute for Genomics and Proteomics, University of California, Los Angeles, Los Angeles, CA 90095, USA
- Molecular Biology Institute, University of California, Los Angeles, Los Angeles, CA 90095, USA
- Quantitative and Computational Biology Institute, University of California, Los Angeles, Los Angeles, CA 90095, USA
- Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Kimberly Insigne
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, CA 90095, USA
- UCLA-DOE Institute for Genomics and Proteomics, University of California, Los Angeles, Los Angeles, CA 90095, USA
- Molecular Biology Institute, University of California, Los Angeles, Los Angeles, CA 90095, USA
- Quantitative and Computational Biology Institute, University of California, Los Angeles, Los Angeles, CA 90095, USA
- Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Marielle L Bond
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Curriculum in Genetics and Molecular Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Nanxiang Zhao
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48109, USA
| | - Alan P Boyle
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48109, USA
- Department of Human Genetics, University of Michigan, Ann Arbor, MI 48109, USA
| | - Douglas H Phanstiel
- Curriculum in Genetics and Molecular Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Thurston Arthritis Research Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Department of Cell Biology and Physiology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Michael I Love
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Department of Biostatistics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - Jason L Stein
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Neuroscience Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | - W Brad Ruzicka
- Laboratory for Epigenomics in Human Psychopathology, McLean Hospital, Belmont, MA 02141, USA
- Harvard Medical School, Boston, MA 02115, USA
- Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | | | - Sriram Kosuri
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, CA 90095, USA
- UCLA-DOE Institute for Genomics and Proteomics, University of California, Los Angeles, Los Angeles, CA 90095, USA
- Molecular Biology Institute, University of California, Los Angeles, Los Angeles, CA 90095, USA
- Quantitative and Computational Biology Institute, University of California, Los Angeles, Los Angeles, CA 90095, USA
- Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Hyejung Won
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Neuroscience Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| |
Collapse
|
34
|
Xiang R, Fang L, Liu S, Macleod IM, Liu Z, Breen EJ, Gao Y, Liu GE, Tenesa A, Mason BA, Chamberlain AJ, Wray NR, Goddard ME. Gene expression and RNA splicing explain large proportions of the heritability for complex traits in cattle. Cell Genom 2023; 3:100385. [PMID: 37868035 PMCID: PMC10589627 DOI: 10.1016/j.xgen.2023.100385] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/08/2022] [Revised: 08/10/2022] [Accepted: 07/26/2023] [Indexed: 10/24/2023]
Abstract
Many quantitative trait loci (QTLs) are in non-coding regions. Therefore, QTLs are assumed to affect gene regulation. Gene expression and RNA splicing are primary steps of transcription, so DNA variants changing gene expression (eVariants) or RNA splicing (sVariants) are expected to significantly affect phenotypes. We quantify the contribution of eVariants and sVariants detected from 16 tissues (n = 4,725) to 37 traits of ∼120,000 cattle (average magnitude of genetic correlation between traits = 0.13). Analyzed in Bayesian mixture models, averaged across 37 traits, cis and trans eVariants and sVariants detected from 16 tissues jointly explain 69.2% (SE = 0.5%) of heritability, 44% more than expected from the same number of random variants. This 69.2% includes an average of 24% from trans e-/sVariants (14% more than expected). Averaged across 56 lipidomic traits, multi-tissue cis and trans e-/sVariants also explain 71.5% (SE = 0.3%) of heritability, demonstrating the essential role of proximal and distal regulatory variants in shaping mammalian phenotypes.
Collapse
Affiliation(s)
- Ruidong Xiang
- Faculty of Veterinary & Agricultural Science, the University of Melbourne, Parkville, VIC 3052, Australia
- Agriculture Victoria, AgriBio, Centre for AgriBiosciences, Bundoora, VIC 3083, Australia
- Cambridge-Baker Systems Genomics Initiative, Baker Heart and Diabetes Institute, Melbourne, VIC 3004, Australia
| | - Lingzhao Fang
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, the University of Edinburgh, Edinburgh, UK
- Center for Quantitative Genetics and Genomics, Aarhus University, Aarhus, Denmark
| | - Shuli Liu
- Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
| | - Iona M. Macleod
- Agriculture Victoria, AgriBio, Centre for AgriBiosciences, Bundoora, VIC 3083, Australia
| | - Zhiqian Liu
- Agriculture Victoria, AgriBio, Centre for AgriBiosciences, Bundoora, VIC 3083, Australia
| | - Edmond J. Breen
- Agriculture Victoria, AgriBio, Centre for AgriBiosciences, Bundoora, VIC 3083, Australia
| | - Yahui Gao
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
| | - George E. Liu
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
| | - Albert Tenesa
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, the University of Edinburgh, Edinburgh, UK
- The Roslin Institute, Royal (Dick) School of Veterinary Studies, the University of Edinburgh, Midlothian EH25 9RG, UK
| | - CattleGTEx Consortium
- Faculty of Veterinary & Agricultural Science, the University of Melbourne, Parkville, VIC 3052, Australia
- Agriculture Victoria, AgriBio, Centre for AgriBiosciences, Bundoora, VIC 3083, Australia
- Cambridge-Baker Systems Genomics Initiative, Baker Heart and Diabetes Institute, Melbourne, VIC 3004, Australia
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, the University of Edinburgh, Edinburgh, UK
- Center for Quantitative Genetics and Genomics, Aarhus University, Aarhus, Denmark
- Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
- The Roslin Institute, Royal (Dick) School of Veterinary Studies, the University of Edinburgh, Midlothian EH25 9RG, UK
- Institute for Molecular Bioscience, the University of Queensland, Brisbane, QLD 4072, Australia
- Queensland Brain Institute, the University of Queensland, Brisbane, QLD 4072, Australia
- School of Applied Systems Biology, La Trobe University, Bundoora, VIC 3083, Australia
| | - Brett A. Mason
- Agriculture Victoria, AgriBio, Centre for AgriBiosciences, Bundoora, VIC 3083, Australia
| | - Amanda J. Chamberlain
- Agriculture Victoria, AgriBio, Centre for AgriBiosciences, Bundoora, VIC 3083, Australia
- School of Applied Systems Biology, La Trobe University, Bundoora, VIC 3083, Australia
| | - Naomi R. Wray
- Institute for Molecular Bioscience, the University of Queensland, Brisbane, QLD 4072, Australia
- Queensland Brain Institute, the University of Queensland, Brisbane, QLD 4072, Australia
| | - Michael E. Goddard
- Faculty of Veterinary & Agricultural Science, the University of Melbourne, Parkville, VIC 3052, Australia
- Agriculture Victoria, AgriBio, Centre for AgriBiosciences, Bundoora, VIC 3083, Australia
| |
Collapse
|
35
|
van Duijvenboden S, Ramírez J, Young WJ, Olczak KJ, Ahmed F, Alhammadi MJAY, Bell CG, Morris AP, Munroe PB. Integration of genetic fine-mapping and multi-omics data reveals candidate effector genes for hypertension. Am J Hum Genet 2023; 110:1718-1734. [PMID: 37683633 PMCID: PMC10577090 DOI: 10.1016/j.ajhg.2023.08.009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Revised: 08/11/2023] [Accepted: 08/11/2023] [Indexed: 09/10/2023] Open
Abstract
Genome-wide association studies of blood pressure (BP) have identified >1,000 loci, but the effector genes and biological pathways at these loci are mostly unknown. Using published association summary statistics, we conducted annotation-informed fine-mapping incorporating tissue-specific chromatin segmentation and colocalization to identify causal variants and candidate effector genes for systolic BP, diastolic BP, and pulse pressure. We observed 532 distinct signals associated with ≥2 BP traits and 84 with all three. For >20% of signals, a single variant accounted for >75% posterior probability, 65 were missense variants in known (SLC39A8, ADRB2, and DBH) and previously unreported BP candidate genes (NRIP1 and MMP14). In disease-relevant tissues, we colocalized >80 and >400 distinct signals for each BP trait with cis-eQTLs and regulatory regions from promoter capture Hi-C, respectively. Integrating mouse, human disorder, gene expression and tissue abundance data, and literature review, we provide consolidated evidence for 436 BP candidate genes for future functional validation and discover several potential drug targets.
Collapse
Affiliation(s)
- Stefan van Duijvenboden
- William Harvey Research Institute, Barts and the London Faculty of Medicine and Dentistry, Queen Mary University of London, EC1M 6BQ London, UK; Institute of Cardiovascular Science, University College London, London, UK; Nuffield Department of Population Health, University of Oxford, Oxford, UK
| | - Julia Ramírez
- William Harvey Research Institute, Barts and the London Faculty of Medicine and Dentistry, Queen Mary University of London, EC1M 6BQ London, UK; Aragon Institute of Engineering Research, University of Zaragoza, Zaragoza, Spain; Centro de Investigación Biomédica en Red - Bioingeniería, Biomateriales y Nanomedicina, Zaragoza, Spain
| | - William J Young
- William Harvey Research Institute, Barts and the London Faculty of Medicine and Dentistry, Queen Mary University of London, EC1M 6BQ London, UK; Barts Heart Centre, St Bartholomew's Hospital, EC1A 7BE London, UK
| | - Kaya J Olczak
- William Harvey Research Institute, Barts and the London Faculty of Medicine and Dentistry, Queen Mary University of London, EC1M 6BQ London, UK
| | - Farah Ahmed
- William Harvey Research Institute, Barts and the London Faculty of Medicine and Dentistry, Queen Mary University of London, EC1M 6BQ London, UK
| | | | - Christopher G Bell
- William Harvey Research Institute, Barts and the London Faculty of Medicine and Dentistry, Queen Mary University of London, EC1M 6BQ London, UK
| | - Andrew P Morris
- Centre for Genetics and Genomics Versus Arthritis, Centre for Musculoskeletal Research, The University of Manchester, Manchester, UK; National Institute of Health and Care Research, Manchester Biomedical Research Centre, Manchester University NHS Foundation Trust, Manchester Academic Health Science Centre, Manchester, UK.
| | - Patricia B Munroe
- William Harvey Research Institute, Barts and the London Faculty of Medicine and Dentistry, Queen Mary University of London, EC1M 6BQ London, UK; National Institute of Health and Care Research, Barts Cardiovascular Biomedical Research Centre, Queen Mary University of London, EC1M 6BQ London, UK.
| |
Collapse
|
36
|
Chepelev I, Harley IT, Harley JB. Modeling of horizontal pleiotropy identifies possible causal gene expression in systemic lupus erythematosus. Front Lupus 2023; 1:1234578. [PMID: 37799268 PMCID: PMC10554754 DOI: 10.3389/flupu.2023.1234578] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/07/2023]
Abstract
Background Systemic lupus erythematosus (SLE) is a chronic autoimmune condition with complex causes involving genetic and environmental factors. While genome-wide association studies (GWASs) have identified genetic loci associated with SLE, the functional genomic elements responsible for disease development remain largely unknown. Mendelian Randomization (MR) is an instrumental variable approach to causal inference based on data from observational studies, where genetic variants are employed as instrumental variables (IVs). Methods This study utilized a two-step strategy to identify causal genes for SLE. In the first step, the classical MR method was employed, assuming the absence of horizontal pleiotropy, to estimate the causal effect of gene expression on SLE. In the second step, advanced probabilistic MR methods (PMR-Egger, MRAID, and MR-MtRobin) were applied to the genes identified in the first step, considering horizontal pleiotropy, to filter out false positives. PMR-Egger and MRAID analyses utilized whole blood expression quantitative trait loci (eQTL) and SLE GWAS summary data, while MR-MtRobin analysis used an independent eQTL dataset from multiple immune cell types along with the same SLE GWAS data. Results The initial MR analysis identified 142 genes, including 43 outside of chromosome 6. Subsequently, applying the advanced MR methods reduced the number of genes with significant causal effects on SLE to 66. PMR-Egger, MRAID, and MR-MtRobin, respectively, identified 13, 7, and 16 non-chromosome 6 genes with significant causal effects. All methods identified expression of PHRF1 gene as causal for SLE. A comprehensive literature review was conducted to enhance understanding of the functional roles and mechanisms of the identified genes in SLE development. Conclusions The findings from the three MR methods exhibited overlapping genes with causal effects on SLE, demonstrating consistent results. However, each method also uncovered unique genes due to different modelling assumptions and technical factors, highlighting the complementary nature of the approaches. Importantly, MRAID demonstrated a reduced percentage of causal genes from the Major Histocompatibility complex (MHC) region on chromosome 6, indicating its potential in minimizing false positive findings. This study contributes to unraveling the mechanisms underlying SLE by employing advanced probabilistic MR methods to identify causal genes, thereby enhancing our understanding of SLE pathogenesis.
Collapse
Affiliation(s)
- Iouri Chepelev
- Research Service, US Department of Veterans Affairs Medical Center, Cincinnati, OH, United States
- Cincinnati Education and Research for Veterans Foundation, Cincinnati, OH, United States
| | - Isaac T.W. Harley
- US Department of Veterans Affairs Medical Center, Aurora, CO, United States
- Department of Immunology and Microbiology, University of Colorado School of Medicine, Aurora, CO, United States
- Division of Rheumatology, Department of Medicine, University of Colorado School of Medicine, Aurora, CO, United States
| | - John B. Harley
- Research Service, US Department of Veterans Affairs Medical Center, Cincinnati, OH, United States
- Cincinnati Education and Research for Veterans Foundation, Cincinnati, OH, United States
| |
Collapse
|
37
|
Li X, Li H, Christenson SA, Castro M, Denlinger LC, Erzurum SC, Fahy JV, Gaston BM, Israel E, Jarjour NN, Levy BD, Mauger DT, Moore WC, Zein J, Kaminski N, Wenzel SE, Woodruff PG, Bleecker ER, Meyers DA. Genetic analyses of chr11p15.5 region identify MUC5AC- MUC5B associated with asthma-related phenotypes. J Asthma 2023; 60:1824-1835. [PMID: 36946148 PMCID: PMC10524756 DOI: 10.1080/02770903.2023.2193631] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Accepted: 03/16/2023] [Indexed: 03/23/2023]
Abstract
OBJECTIVE Genome-wide association studies (GWASs) have identified single nucleotide polymorphisms (SNPs) in chr11p15.5 region associated with asthma and idiopathic interstitial pneumonias (IIPs). We sought to identify functional genes for asthma by combining SNPs and mRNA expression in bronchial epithelial cells (BEC) in the Severe Asthma Research Program (SARP). METHODS Correlation analyses of mRNA expression of six candidate genes (AP2A2, MUC6, MUC2, MUC5AC, MUC5B, and TOLLIP) and asthma phenotypes were performed in the longitudinal cohort (n = 156) with RNAseq in BEC, and replicated in the cross-sectional cohort (n = 155). eQTL (n = 114) and genetic association analysis of asthma severity (426 severe vs. 531 non-severe asthma) were performed, and compared with previously published GWASs of IIPs and asthma. RESULTS Higher expression of AP2A2 and MUC5AC and lower expression of MUC5B in BEC were correlated with asthma, asthma exacerbations, and T2 biomarkers (P < 0.01). SNPs associated with asthma and IIPs in previous GWASs were eQTL SNPs for MUC5AC, MUC5B, or TOLLIP, however, they were not in strong linkage disequilibrium. The risk alleles for asthma or protective alleles for IIPs were associated with higher expression of MUC5AC and lower expression of MUC5B. rs11603634, rs12788104, and rs28415845 associated with moderate-to-severe asthma or adult onset asthma in previous GWASs were not associated with asthma severity (P > 0.8). CONCLUSIONS SNPs associated with asthma in chr11p15.5 region are not associated with asthma severity neither with IIPs. Higher expression of MUC5AC and lower expression of MUC5B are risk for asthma but protective for IIPs.
Collapse
Affiliation(s)
- Xingnan Li
- Division of Genetics, Genomics and Precision Medicine, Department of Medicine, University of Arizona, Tucson, Arizona, USA
| | - Huashi Li
- Division of Genetics, Genomics and Precision Medicine, Department of Medicine, University of Arizona, Tucson, Arizona, USA
| | - Stephanie A. Christenson
- Division of Pulmonary, Critical Care, Sleep and Allergy, Department of Medicine, University of California at San Francisco, San Francisco, California, USA
| | - Mario Castro
- Division of Pulmonary, Critical Care and Sleep Medicine, University of Kansas School of Medicine, Kansas City, Kansas, USA
| | - Loren C. Denlinger
- Department of Medicine, University of Wisconsin School of Medicine & Public Health, Madison, Wisconsin, USA
| | - Serpil C. Erzurum
- Lerner Research Institute and the Respiratory Institute, Cleveland Clinic, Cleveland, Ohio, USA
| | - John V. Fahy
- Division of Pulmonary, Critical Care, Sleep and Allergy, Department of Medicine, University of California at San Francisco, San Francisco, California, USA
| | - Benjamin M. Gaston
- Wells Center for Pediatric Research and Riley Hospital for Children, Indiana University, Indianapolis, Indiana, USA
| | - Elliot Israel
- Division of Pulmonary and Critical Care Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, Massachusetts, USA
| | - Nizar N. Jarjour
- Department of Medicine, University of Wisconsin School of Medicine & Public Health, Madison, Wisconsin, USA
| | - Bruce D. Levy
- Division of Pulmonary and Critical Care Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, Massachusetts, USA
| | - David T Mauger
- Department of Public Health Sciences, College of Medicine, Penn State University, Hershey, Pennsylvania, USA
| | - Wendy C. Moore
- Department of Medicine, Wake Forest School of Medicine, Winston-Salem, North Carolina, USA
| | - Joe Zein
- Lerner Research Institute and the Respiratory Institute, Cleveland Clinic, Cleveland, Ohio, USA
| | - Naftali Kaminski
- Division of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, Yale School of Medicine, New Haven, Connecticut, USA
| | - Sally E. Wenzel
- Department of Environmental and Occupational Health, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| | - Prescott G. Woodruff
- Division of Pulmonary, Critical Care, Sleep and Allergy, Department of Medicine, University of California at San Francisco, San Francisco, California, USA
| | - Eugene R. Bleecker
- Division of Genetics, Genomics and Precision Medicine, Department of Medicine, University of Arizona, Tucson, Arizona, USA
| | - Deborah A. Meyers
- Division of Genetics, Genomics and Precision Medicine, Department of Medicine, University of Arizona, Tucson, Arizona, USA
| | | |
Collapse
|
38
|
Zhao T, Wu H, Wang X, Zhao Y, Wang L, Pan J, Mei H, Han J, Wang S, Lu K, Li M, Gao M, Cao Z, Zhang H, Wan K, Li J, Fang L, Zhang T, Guan X. Integration of eQTL and machine learning to dissect causal genes with pleiotropic effects in genetic regulation networks of seed cotton yield. Cell Rep 2023; 42:113111. [PMID: 37676770 DOI: 10.1016/j.celrep.2023.113111] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 06/19/2023] [Accepted: 08/24/2023] [Indexed: 09/09/2023] Open
Abstract
The dissection of a gene regulatory network (GRN) that complements the genome-wide association study (GWAS) locus and the crosstalk underlying multiple agronomical traits remains a major challenge. In this study, we generate 558 transcriptional profiles of lint-bearing ovules at one day post-anthesis from a selective core cotton germplasm, from which 12,207 expression quantitative trait loci (eQTLs) are identified. Sixty-six known phenotypic GWAS loci are colocalized with 1,090 eQTLs, forming 38 functional GRNs associated predominantly with seed yield. Of the eGenes, 34 exhibit pleiotropic effects. Combining the eQTLs within the seed yield GRNs significantly increases the portion of narrow-sense heritability. The extreme gradient boosting (XGBoost) machine learning approach is applied to predict seed cotton yield phenotypes on the basis of gene expression. Top-ranking eGenes (NF-YB3, FLA2, and GRDP1) derived with pleiotropic effects on yield traits are validated, along with their potential roles by correlation analysis, domestication selection analysis, and transgenic plants.
Collapse
Affiliation(s)
- Ting Zhao
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China; Hainan Institute of Zhejiang University, Building 11, Yonyou Industrial Park, Yazhou Bay Science and Technology City, Yazhou District, Sanya 572025, China
| | - Hongyu Wu
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China
| | - Xutong Wang
- Hubei Hongshan Laboratory, Wuhan 430070, China
| | - Yongyan Zhao
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China; Hainan Institute of Zhejiang University, Building 11, Yonyou Industrial Park, Yazhou Bay Science and Technology City, Yazhou District, Sanya 572025, China
| | - Luyao Wang
- Hainan Institute of Zhejiang University, Building 11, Yonyou Industrial Park, Yazhou Bay Science and Technology City, Yazhou District, Sanya 572025, China
| | - Jiaying Pan
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China; Hainan Institute of Zhejiang University, Building 11, Yonyou Industrial Park, Yazhou Bay Science and Technology City, Yazhou District, Sanya 572025, China
| | - Huan Mei
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China
| | - Jin Han
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China
| | - Siyuan Wang
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China
| | - Kening Lu
- State Key Laboratory of Crop Genetics and Germplasm Enhancement, Cotton Hybrid R & D Engineering Center (the Ministry of Education), College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China
| | - Menglin Li
- State Key Laboratory of Crop Genetics and Germplasm Enhancement, Cotton Hybrid R & D Engineering Center (the Ministry of Education), College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China
| | - Mengtao Gao
- State Key Laboratory of Crop Genetics and Germplasm Enhancement, Cotton Hybrid R & D Engineering Center (the Ministry of Education), College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China
| | - Zeyi Cao
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China
| | - Hailin Zhang
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China
| | - Ke Wan
- State Key Laboratory of Crop Genetics and Germplasm Enhancement, Cotton Hybrid R & D Engineering Center (the Ministry of Education), College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China
| | - Jie Li
- State Key Laboratory of Crop Genetics and Germplasm Enhancement, Cotton Hybrid R & D Engineering Center (the Ministry of Education), College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China
| | - Lei Fang
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China; Hainan Institute of Zhejiang University, Building 11, Yonyou Industrial Park, Yazhou Bay Science and Technology City, Yazhou District, Sanya 572025, China
| | - Tianzhen Zhang
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China; Hainan Institute of Zhejiang University, Building 11, Yonyou Industrial Park, Yazhou Bay Science and Technology City, Yazhou District, Sanya 572025, China
| | - Xueying Guan
- Zhejiang Provincial Key Laboratory of Crop Genetic Resources, The Advanced Seed Institute, Plant Precision Breeding Academy, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 300058, China; Hainan Institute of Zhejiang University, Building 11, Yonyou Industrial Park, Yazhou Bay Science and Technology City, Yazhou District, Sanya 572025, China.
| |
Collapse
|
39
|
Kim MS, Song M, Kim B, Shim I, Kim DS, Natarajan P, Do R, Won HH. Prioritization of therapeutic targets for dyslipidemia using integrative multi-omics and multi-trait analysis. Cell Rep Med 2023; 4:101112. [PMID: 37582372 PMCID: PMC10518515 DOI: 10.1016/j.xcrm.2023.101112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2022] [Revised: 12/22/2022] [Accepted: 06/19/2023] [Indexed: 08/17/2023]
Abstract
Drug targets with genetic support are several-fold more likely to succeed in clinical trials. We introduce a genetic-driven approach based on causal inferences that can inform drug target prioritization, repurposing, and adverse effects of using lipid-lowering agents. Given that a multi-trait approach increases the power to detect meaningful variants/genes, we conduct multi-omics and multi-trait analyses, followed by network connectivity investigations, and prioritize 30 potential therapeutic targets for dyslipidemia, including SORT1, PSRC1, CELSR2, PCSK9, HMGCR, APOB, GRN, HFE2, FJX1, C1QTNF1, and SLC5A8. 20% (6/30) of prioritized targets from our hypothesis-free drug target search are either approved or under investigation for dyslipidemia. The prioritized targets are 22-fold higher in likelihood of being approved or under investigation in clinical trials than genome-wide association study (GWAS)-curated targets. Our results demonstrate that the genetic-driven approach used in this study is a promising strategy for prioritizing targets while informing about the potential adverse effects and repurposing opportunities.
Collapse
Affiliation(s)
- Min Seo Kim
- Samsung Advanced Institute for Health Sciences & Technology (SAIHST), Sungkyunkwan University, Samsung Medical Center, Seoul, Republic of Korea
| | - Minku Song
- Samsung Advanced Institute for Health Sciences & Technology (SAIHST), Sungkyunkwan University, Samsung Medical Center, Seoul, Republic of Korea
| | - Beomsu Kim
- Samsung Advanced Institute for Health Sciences & Technology (SAIHST), Sungkyunkwan University, Samsung Medical Center, Seoul, Republic of Korea
| | - Injeong Shim
- Samsung Advanced Institute for Health Sciences & Technology (SAIHST), Sungkyunkwan University, Samsung Medical Center, Seoul, Republic of Korea
| | - Dan Say Kim
- Samsung Advanced Institute for Health Sciences & Technology (SAIHST), Sungkyunkwan University, Samsung Medical Center, Seoul, Republic of Korea
| | - Pradeep Natarajan
- Medical and Population Genetics and Cardiovascular Disease Initiative, Broad Institute of MIT and Harvard, Cambridge, MA, USA; Cardiovascular Research Center, Massachusetts General Hospital, Boston, MA, USA; Department of Medicine, Harvard Medical School, Boston, MA, USA
| | - Ron Do
- The Charles Bronfman Institute for Personalized Medicine, Icahn School of Medicine at Mount Sinai, New York, NY, USA; Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Hong-Hee Won
- Samsung Advanced Institute for Health Sciences & Technology (SAIHST), Sungkyunkwan University, Samsung Medical Center, Seoul, Republic of Korea; Samsung Genome Institute, Samsung Medical Center, Seoul, Republic of Korea.
| |
Collapse
|
40
|
Li X, Shen A, Zhao Y, Xia J. Mendelian Randomization Using the Druggable Genome Reveals Genetically Supported Drug Targets for Psychiatric Disorders. Schizophr Bull 2023; 49:1305-1315. [PMID: 37418754 PMCID: PMC10483453 DOI: 10.1093/schbul/sbad100] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 07/09/2023]
Abstract
BACKGROUND AND HYPOTHESIS Psychiatric disorders impose a huge health and economic burden on modern society. However, there is currently no proven completely effective treatment available, partly owing to the inefficiency of drug target identification and validation. We aim to identify therapeutic targets relevant to psychiatric disorders by conducting Mendelian randomization (MR) analysis. STUDY DESIGN We performed genome-wide MR analysis by integrating expression quantitative trait loci (eQTL) of 4479 actionable genes that encode druggable proteins and genetic summary statistics from genome-wide association studies of psychiatric disorders. After conducting colocalization analysis on the brain MR findings, we employed protein quantitative trait loci (pQTL) data as genetic proposed instruments for intersecting the colocalized genes to provide further genetic evidence. STUDY RESULTS By performing MR and colocalization analysis with eQTL genetic instruments, we obtained 31 promising drug targets for psychiatric disorders, including 21 significant genes for schizophrenia, 7 for bipolar disorder, 2 for depression, 1 for attention deficit and hyperactivity (ADHD) and none for autism spectrum disorder. Combining MR results using pQTL genetic instruments, we finally proposed 8 drug-targeting genes supported by the strongest MR evidence, including gene ACE, BTN3A3, HAPLN4, MAPK3 and NEK4 for schizophrenia, gene NEK4 and HAPLN4 for bipolar disorder, and gene TIE1 for ADHD. CONCLUSIONS Our findings with genetic support were more likely to be to succeed in clinical trials. In addition, our study prioritizes approved drug targets for the development of new therapies and provides critical drug reuse opportunities for psychiatric disorders.
Collapse
Affiliation(s)
- Xiaoyan Li
- Key Laboratory of Intelligent Computing and Signal Processing of Ministry of Education and Information Materials and Intelligent Sensing Laboratory of Anhui Province, Institutes of Physical Science and Information Technology, Anhui University, Hefei, Anhui 230601, China
| | - Aotian Shen
- Key Laboratory of Intelligent Computing and Signal Processing of Ministry of Education and Information Materials and Intelligent Sensing Laboratory of Anhui Province, Institutes of Physical Science and Information Technology, Anhui University, Hefei, Anhui 230601, China
| | - Yiran Zhao
- Key Laboratory of Intelligent Computing and Signal Processing of Ministry of Education and Information Materials and Intelligent Sensing Laboratory of Anhui Province, Institutes of Physical Science and Information Technology, Anhui University, Hefei, Anhui 230601, China
| | - Junfeng Xia
- Key Laboratory of Intelligent Computing and Signal Processing of Ministry of Education and Information Materials and Intelligent Sensing Laboratory of Anhui Province, Institutes of Physical Science and Information Technology, Anhui University, Hefei, Anhui 230601, China
| |
Collapse
|
41
|
Xu C, Song LY, Zhou Y, Ma DN, Ding QS, Guo ZJ, Li J, Song SW, Zhang LD, Zheng HL. Integration of eQTL and GWAS analysis uncovers a genetic regulation of natural ionomic variation in Arabidopsis. Plant Cell Rep 2023; 42:1473-1485. [PMID: 37516984 DOI: 10.1007/s00299-023-03042-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/27/2022] [Accepted: 06/12/2023] [Indexed: 08/01/2023]
Abstract
KEY MESSAGE This study provided important insights into the genetic architecture of variations in A. thaliana leaf ionome in a cell-type-specific manner. The functional interpretation of traits associated variants by expression quantitative trait loci (eQTL) analysis is usually performed in bulk tissue samples. While the regulation of gene expression is context-dependent, such as cell-type-specific manner. In this study, we estimated cell-type abundances from 728 bulk tissue samples using single-cell RNA-sequencing dataset, and performed cis-eQTL mapping to identify cell-type-interaction eQTL (cis-eQTLs(ci)) in A. thaliana. Also, we performed Genome-wide association studies (GWAS) analyses for 999 accessions to identify the genetic basis of variations in A. thaliana leaf ionome. As a result, a total of 5,664 unique eQTL genes and 15,038 unique cis-eQTLs(ci) were significant. The majority (62.83%) of cis-eQTLs(ci) were cell-type-specific eQTLs. Using colocalization, we uncovered one interested gene AT2G25590 in Phloem cell, encoding a kind of plant Tudor-like protein with possible chromatin-associated functions, which colocalized with the most significant cis-eQTL(ci) of a Mo-related locus (Chr2:10,908,806:A:C; P = 3.27 × 10-27). Furthermore, we prioritized eight target genes associated with AT2G25590, which were previously reported in regulating the concentration of Mo element in A. thaliana. This study revealed the genetic regulation of ionomic variations and provided a foundation for further studies on molecular mechanisms of genetic variants controlling the A. thaliana ionome.
Collapse
Affiliation(s)
- Chaoqun Xu
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, 361104, China
| | - Ling-Yu Song
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, 361104, China
| | - Ying Zhou
- School of Medicine, National Institute for Data Science in Health and Medicine, Xiamen University, Xiamen, 361102, China
| | - Dong-Na Ma
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, 361104, China
- National Engineering Research Center of Cereal Fermentation and Food Biomanufacturing, School of Food Science and Technology, Jiangnan University, Wuxi, 214122, Jiangsu, China
| | - Qian-Su Ding
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, 361104, China
| | - Ze-Jun Guo
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, 361104, China
| | - Jing Li
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, 361104, China
| | - Shi-Wei Song
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, 361104, China
| | - Lu-Dan Zhang
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, 361104, China
| | - Hai-Lei Zheng
- Key Laboratory of the Ministry of Education for Coastal and Wetland Ecosystems, College of the Environment and Ecology, Xiamen University, Xiamen, 361104, China.
| |
Collapse
|
42
|
Abstract
Recent advancements in single-cell technologies have enabled expression quantitative trait locus (eQTL) analysis across many individuals at single-cell resolution. Compared with bulk RNA sequencing, which averages gene expression across cell types and cell states, single-cell assays capture the transcriptional states of individual cells, including fine-grained, transient, and difficult-to-isolate populations at unprecedented scale and resolution. Single-cell eQTL (sc-eQTL) mapping can identify context-dependent eQTLs that vary with cell states, including some that colocalize with disease variants identified in genome-wide association studies. By uncovering the precise contexts in which these eQTLs act, single-cell approaches can unveil previously hidden regulatory effects and pinpoint important cell states underlying molecular mechanisms of disease. Here, we present an overview of recently deployed experimental designs in sc-eQTL studies. In the process, we consider the influence of study design choices such as cohort, cell states, and ex vivo perturbations. We then discuss current methodologies, modeling approaches, and technical challenges as well as future opportunities and applications.
Collapse
Affiliation(s)
- Joyce B Kang
- Center for Data Sciences and Divisions of Genetics and Rheumatology, Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts, USA; ,
- Department of Biomedical Informatics, Harvard Medical School, Boston, Massachusetts, USA
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA;
| | | | - Aparna Nathan
- Center for Data Sciences and Divisions of Genetics and Rheumatology, Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts, USA; ,
- Department of Biomedical Informatics, Harvard Medical School, Boston, Massachusetts, USA
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA;
| | - Nicole Soranzo
- Human Technopole, Milan, Italy; ,
- Department of Human Genetics, Wellcome Sanger Institute, Hinxton, United Kingdom
- British Heart Foundation Centre of Research Excellence and Department of Haematology, University of Cambridge, Cambridge, United Kingdom
| | - Soumya Raychaudhuri
- Center for Data Sciences and Divisions of Genetics and Rheumatology, Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts, USA; ,
- Department of Biomedical Informatics, Harvard Medical School, Boston, Massachusetts, USA
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA;
- Centre for Genetics and Genomics Versus Arthritis, University of Manchester, Manchester, United Kingdom
| |
Collapse
|
43
|
Brookes KJ. Evaluating the Classification Accuracy of Expression Quantitative Trait Loci Calculated Polygenic Risk Scores in Alzheimer's Disease. Int J Mol Sci 2023; 24:12799. [PMID: 37628980 PMCID: PMC10454324 DOI: 10.3390/ijms241612799] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Revised: 08/09/2023] [Accepted: 08/12/2023] [Indexed: 08/27/2023] Open
Abstract
Polygenic risk scores (PRS) hold promise for the early identification of those at risk from neurodegenerative disorders such as Alzheimer's Disease (AD), allowing for intervention to occur prior to neuronal damage. The current selection of informative single nucleotide polymorphisms (SNPs) to generate the risk scores is based on the modelling of large genome-wide association data using significance thresholds. However, the biological relevance of these SNPs is largely unknown. This study, in contrast, aims to identify SNPs with biological relevance to AD and then assess them for their ability to accurately classify cases and controls. Samples selected from the Brains for Dementia Research (BDR) were used to produce gene expression data to identify potential expression quantitative trait loci (eQTLs) relevant to AD. These SNPs were then incorporated into a PRS model to classify AD and controls in the full BDR cohort. Models derived from these eQTLs demonstrate modest classification potential with an accuracy between 61% and 67%. Although the model accuracy is not as high as some values in the literature based on significance thresholds from genome-wide association studies, these models may reflect a more biologically relevant model, which may provide novel targets for therapeutic intervention.
Collapse
Affiliation(s)
- Keeley J Brookes
- Department of Biosciences, School of Science & Technology, Nottingham Trent University, Nottingham NG11 8NS, UK
| |
Collapse
|
44
|
Tian Y, Dong D, Wang Z, Wu L, Park JY, Wei GH, Wang L. Combined CRISPRi and proteomics screening reveal a cohesin-CTCF-bound allele contributing to increased expression of RUVBL1 and prostate cancer progression. Am J Hum Genet 2023; 110:1289-1303. [PMID: 37541187 PMCID: PMC10432188 DOI: 10.1016/j.ajhg.2023.07.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Revised: 07/06/2023] [Accepted: 07/06/2023] [Indexed: 08/06/2023] Open
Abstract
Genome-wide association studies along with expression quantitative trait locus (eQTL) mapping have identified hundreds of single-nucleotide polymorphisms (SNPs) and their target genes in prostate cancer (PCa), yet functional characterization of these risk loci remains challenging. To screen for potential regulatory SNPs, we designed a CRISPRi library containing 9,133 guide RNAs (gRNAs) to cover 2,166 candidate SNP loci implicated in PCa and identified 117 SNPs that could regulate 90 genes for PCa cell growth advantage. Among these, rs60464856 was covered by multiple gRNAs significantly depleted in screening (FDR < 0.05). Pooled SNP association analysis in the PRACTICAL and FinnGen cohorts showed significantly higher PCa risk for the rs60464856 G allele (p value = 1.2 × 10-16 and 3.2 × 10-7, respectively). Subsequent eQTL analysis revealed that the G allele is associated with increased RUVBL1 expression in multiple datasets. Further CRISPRi and xCas9 base editing confirmed that the rs60464856 G allele leads to elevated RUVBL1 expression. Furthermore, SILAC-based proteomic analysis demonstrated allelic binding of cohesin subunits at the rs60464856 region, where the HiC dataset showed consistent chromatin interactions in prostate cell lines. RUVBL1 depletion inhibited PCa cell proliferation and tumor growth in a xenograft mouse model. Gene-set enrichment analysis suggested an association of RUVBL1 expression with cell-cycle-related pathways. Increased expression of RUVBL1 and activation of cell-cycle pathways were correlated with poor PCa survival in TCGA datasets. Our CRISPRi screening prioritized about one hundred regulatory SNPs essential for prostate cell proliferation. In combination with proteomics and functional studies, we characterized the mechanistic role of rs60464856 and RUVBL1 in PCa progression.
Collapse
Affiliation(s)
- Yijun Tian
- Department of Tumor Biology, H. Lee Moffitt Cancer Center and Research Institute, 12902 Magnolia Drive, Tampa, FL 33612, USA
| | - Dandan Dong
- MOE Key Laboratory of Metabolism and Molecular Medicine, Shanghai Medical College of Fudan University, Shanghai, China
| | - Zixian Wang
- MOE Key Laboratory of Metabolism and Molecular Medicine, Shanghai Medical College of Fudan University, Shanghai, China; Department of Biochemistry and Molecular Biology of School of Basic Medical Sciences, Shanghai Medical College of Fudan University, Shanghai, China; Fudan University Shanghai Cancer Center, Shanghai Medical College of Fudan University, Shanghai, China
| | - Lang Wu
- Population Sciences in the Pacific Program, University of Hawaii Cancer Center, University of Hawaii at Manoa, Honolulu, HI 96813, USA
| | - Jong Y Park
- Department of Cancer Epidemiology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA
| | - Gong-Hong Wei
- MOE Key Laboratory of Metabolism and Molecular Medicine, Shanghai Medical College of Fudan University, Shanghai, China; Department of Biochemistry and Molecular Biology of School of Basic Medical Sciences, Shanghai Medical College of Fudan University, Shanghai, China; Fudan University Shanghai Cancer Center, Shanghai Medical College of Fudan University, Shanghai, China; Disease Networks Research Unit, Biocenter Oulu, University of Oulu, Oulu, Finland; Faculty of Biochemistry and Molecular Medicine, University of Oulu, Oulu, Finland.
| | - Liang Wang
- Department of Tumor Biology, H. Lee Moffitt Cancer Center and Research Institute, 12902 Magnolia Drive, Tampa, FL 33612, USA.
| |
Collapse
|
45
|
Swart PC, Du Plessis M, Rust C, Womersley JS, van den Heuvel LL, Seedat S, Hemmings SMJ. Identifying genetic loci that are associated with changes in gene expression in PTSD in a South African cohort. J Neurochem 2023; 166:705-719. [PMID: 37522158 PMCID: PMC10953375 DOI: 10.1111/jnc.15919] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2023] [Revised: 06/30/2023] [Accepted: 07/05/2023] [Indexed: 08/01/2023]
Abstract
The molecular mechanisms underlying posttraumatic stress disorder (PTSD) are yet to be fully elucidated, especially in underrepresented population groups. Expression quantitative trait loci (eQTLs) are DNA sequence variants that influence gene expression, in a local (cis-) or distal (trans-) manner, and subsequently impact cellular, tissue, and system physiology. This study aims to identify genetic loci associated with gene expression changes in a South African PTSD cohort. Genome-wide genotype and RNA-sequencing data were obtained from 32 trauma-exposed controls and 35 PTSD cases of mixed-ancestry, as part of the SHARED ROOTS project. The first approach utilised 108 937 single-nucleotide polymorphisms (SNPs) (MAF > 10%) and 11 312 genes with Matrix eQTL to map potential eQTLs, while controlling for covariates as appropriate. The second analysis was focused on 5638 SNPs related to a previously calculated PTSD polygenic risk score for this cohort. SNP-gene pairs were considered eQTLs if they surpassed Bonferroni correction and had a false discovery rate <0.05. We did not identify eQTLs that significantly influenced gene expression in a PTSD-dependent manner. However, several known cis-eQTLs, independent of PTSD diagnosis, were observed. rs8521 (C > T) was associated with TAGLN and SIDT2 expression, and rs11085906 (C > T) was associated with ZNF333 expression. This exploratory study provides insight into the molecular mechanisms associated with PTSD in a non-European, admixed sample population. This study was limited by the cross-sectional design and insufficient statistical power. Overall, this study should encourage further multi-omics approaches towards investigating PTSD in diverse populations.
Collapse
Affiliation(s)
- Patricia C. Swart
- Department of Psychiatry, Faculty of Medicine and Health SciencesStellenbosch UniversityCape TownSouth Africa
- South African Medical Research Council/Stellenbosch University Genomics of Brain Disorders UnitCape TownSouth Africa
| | - Morne Du Plessis
- Department of Psychiatry, Faculty of Medicine and Health SciencesStellenbosch UniversityCape TownSouth Africa
- South African Medical Research Council/Stellenbosch University Genomics of Brain Disorders UnitCape TownSouth Africa
| | - Carlien Rust
- Department of Psychiatry, Faculty of Medicine and Health SciencesStellenbosch UniversityCape TownSouth Africa
- South African Medical Research Council/Stellenbosch University Genomics of Brain Disorders UnitCape TownSouth Africa
| | - Jacqueline S. Womersley
- Department of Psychiatry, Faculty of Medicine and Health SciencesStellenbosch UniversityCape TownSouth Africa
- South African Medical Research Council/Stellenbosch University Genomics of Brain Disorders UnitCape TownSouth Africa
| | - Leigh L. van den Heuvel
- Department of Psychiatry, Faculty of Medicine and Health SciencesStellenbosch UniversityCape TownSouth Africa
- South African Medical Research Council/Stellenbosch University Genomics of Brain Disorders UnitCape TownSouth Africa
| | - Soraya Seedat
- Department of Psychiatry, Faculty of Medicine and Health SciencesStellenbosch UniversityCape TownSouth Africa
- South African Medical Research Council/Stellenbosch University Genomics of Brain Disorders UnitCape TownSouth Africa
| | - Sian M. J. Hemmings
- Department of Psychiatry, Faculty of Medicine and Health SciencesStellenbosch UniversityCape TownSouth Africa
- South African Medical Research Council/Stellenbosch University Genomics of Brain Disorders UnitCape TownSouth Africa
| |
Collapse
|
46
|
Tsouris A, Brach G, Friedrich A, Hou J, Schacherer J. Diallel panel reveals a significant impact of low-frequency genetic variants on gene expression variation in yeast. bioRxiv 2023:2023.07.21.550015. [PMID: 37503053 PMCID: PMC10370210 DOI: 10.1101/2023.07.21.550015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
Unraveling the genetic sources of gene expression variation is essential to better understand the origins of phenotypic diversity in natural populations. Genome-wide association studies identified thousands of variants involved in gene expression variation, however, variants detected only explain part of the heritability. In fact, variants such as low-frequency and structural variants (SVs) are poorly captured in association studies. To assess the impact of these variants on gene expression variation, we explored a half-diallel panel composed of 323 hybrids originated from pairwise crosses of 26 natural Saccharomyces cerevisiae isolates. Using short- and long-read sequencing strategies, we established an exhaustive catalog of single nucleotide polymorphisms (SNPs) and SVs for this panel. Combining this dataset with the transcriptomes of all hybrids, we comprehensively mapped SNPs and SVs associated with gene expression variation. While SVs impact gene expression variation, SNPs exhibit a higher effect size with an overrepresentation of low-frequency variants compared to common ones. These results reinforce the importance of dissecting the heritability of complex traits with a comprehensive catalog of genetic variants at the population level.
Collapse
Affiliation(s)
- Andreas Tsouris
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Gauthier Brach
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Anne Friedrich
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Jing Hou
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Joseph Schacherer
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
- Institut Universitaire de France (IUF), Paris, France
| |
Collapse
|
47
|
Tsouris A, Brach G, Schacherer J, Hou J. Non-additive genetic components contribute significantly to population-wide gene expression variation. bioRxiv 2023:2023.07.21.550013. [PMID: 37546809 PMCID: PMC10401925 DOI: 10.1101/2023.07.21.550013] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/08/2023]
Abstract
Gene expression variation, an essential step between genomic variation and phenotypic landscape, is collectively controlled by local (cis) and distant (trans) regulatory changes. Nevertheless, how these regulatory elements differentially influence the heritability of expression traits remains unclear. Here, we bridge this gap by analyzing the transcriptomes of a large diallel panel consisting of 323 unique hybrids originated from genetically divergent yeast isolates. We estimated the broad- and narrow-sense heritability across 5,087 transcript abundance traits and showed that non-additive components account for 36% of the phenotypic variance on average. By comparing allelic expression ratios in the hybrid and the corresponding parental pair, we identified regulatory changes in 25% of all cases, with a majority acting in trans. We further showed that trans-regulation could underlie coordinated expression variation across highly connected genes, resulting in significantly higher non-additive variance and most likely in some of the missing heritability of gene expression traits.
Collapse
Affiliation(s)
- Andreas Tsouris
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Gauthier Brach
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Joseph Schacherer
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
- Institut Universitaire de France (IUF), Paris, France
| | - Jing Hou
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| |
Collapse
|
48
|
Castanera R, Morales-Díaz N, Gupta S, Purugganan M, Casacuberta JM. Transposons are important contributors to gene expression variability under selection in rice populations. eLife 2023; 12:RP86324. [PMID: 37467142 PMCID: PMC10393045 DOI: 10.7554/elife.86324] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/21/2023] Open
Abstract
Transposable elements (TEs) are an important source of genome variability. Here, we analyze their contribution to gene expression variability in rice by performing a TE insertion polymorphism expression quantitative trait locus mapping using expression data from 208 varieties from the Oryza sativa ssp. indica and O. sativa ssp. japonica subspecies. Our data show that TE insertions are associated with changes of expression of many genes known to be targets of rice domestication and breeding. An important fraction of these insertions were already present in the rice wild ancestors, and have been differentially selected in indica and japonica rice populations. Taken together, our results show that small changes of expression in signal transduction genes induced by TE insertions accompany the domestication and adaptation of rice populations.
Collapse
Affiliation(s)
- Raúl Castanera
- Centre for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus UAB, Cerdanyola del Vallès, Barcelona, Spain
| | - Noemia Morales-Díaz
- Centre for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus UAB, Cerdanyola del Vallès, Barcelona, Spain
| | - Sonal Gupta
- Center for Genomics and Systems Biology, New York University, New York, United States
| | - Michael Purugganan
- Center for Genomics and Systems Biology, New York University, New York, United States
- Center for Genomics and Systems Biology, New York University Abu Dhabi, Saadiyat Island, Abu Dhabi, United Arab Emirates
| | - Josep M Casacuberta
- Centre for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus UAB, Cerdanyola del Vallès, Barcelona, Spain
| |
Collapse
|
49
|
Sadler MC, Auwerx C, Deelen P, Kutalik Z. Multi-layered genetic approaches to identify approved drug targets. Cell Genom 2023; 3:100341. [PMID: 37492104 PMCID: PMC10363916 DOI: 10.1016/j.xgen.2023.100341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 04/25/2023] [Accepted: 05/16/2023] [Indexed: 07/27/2023]
Abstract
Drugs targeting genes linked to disease via evidence from human genetics have increased odds of approval. Approaches to prioritize such genes include genome-wide association studies (GWASs), rare variant burden tests in exome sequencing studies (Exome), or integration of a GWAS with expression/protein quantitative trait loci (eQTL/pQTL-GWAS). Here, we compare gene-prioritization approaches on 30 clinically relevant traits and benchmark their ability to recover drug targets. Across traits, prioritized genes were enriched for drug targets with odds ratios (ORs) of 2.17, 2.04, 1.81, and 1.31 for the GWAS, eQTL-GWAS, Exome, and pQTL-GWAS methods, respectively. Adjusting for differences in testable genes and sample sizes, GWAS outperforms e/pQTL-GWAS, but not the Exome approach. Furthermore, performance increased through gene network diffusion, although the node degree, being the best predictor (OR = 8.7), revealed strong bias in literature-curated networks. In conclusion, we systematically assessed strategies to prioritize drug target genes, highlighting the promises and pitfalls of current approaches.
Collapse
Affiliation(s)
- Marie C. Sadler
- University Center for Primary Care and Public Health, Route de Berne 113, 1010 Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Quartier Sorge, 1015 Lausanne, Switzerland
- Department of Computational Biology, University of Lausanne, 1015 Lausanne, Switzerland
| | - Chiara Auwerx
- University Center for Primary Care and Public Health, Route de Berne 113, 1010 Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Quartier Sorge, 1015 Lausanne, Switzerland
- Department of Computational Biology, University of Lausanne, 1015 Lausanne, Switzerland
- Center for Integrative Genomics, University of Lausanne, 1015 Lausanne, Switzerland
| | - Patrick Deelen
- Department of Genetics, University of Groningen, 9700 Groningen, the Netherlands
- Oncode Institute, 3521 Utrecht, the Netherlands
| | - Zoltán Kutalik
- University Center for Primary Care and Public Health, Route de Berne 113, 1010 Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Quartier Sorge, 1015 Lausanne, Switzerland
- Department of Computational Biology, University of Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
50
|
Drown MK, Oleksiak MF, Crawford DL. Trans-Acting Genotypes Associted with mRNA Expression Affect Metabolic and Thermal Tolerance Traits. Genome Biol Evol 2023:evad123. [PMID: 37392472 PMCID: PMC10370451 DOI: 10.1093/gbe/evad123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Revised: 06/23/2023] [Accepted: 06/28/2023] [Indexed: 07/03/2023] Open
Abstract
Evolutionary processes driving physiological trait variation depend on the underlying genomic mechanisms. Evolution of these mechanisms depends on the genetic complexity (involving many genes) and how gene expression impacting the traits is converted to phenotype. Yet, genomic mechanisms that impact physiological traits are diverse and context dependent (e.g., vary by environment, tissues), making them difficult to discern. We examine the relationships between genotype, mRNA expression, and physiological traits to discern the genetic complexity and whether the gene expression affecting the physiological traits is primarily cis or trans-acting. We use low-coverage whole genome sequencing and heart- or brain-specific mRNA expression to identify polymorphisms directly associated with physiological traits and expressed quantitative trait loci (eQTL) indirectly associated with variation in six temperature specific physiological traits (standard metabolic rate, thermal tolerance, and four substrate specific cardiac metabolic rates). Focusing on a select set of mRNAs belonging to co-expression modules that explain up to 82% of temperature specific traits, we identified hundreds of significant eQTL for mRNA whose expression affects physiological traits. Surprisingly, most eQTL (97.4% for heart and 96.7% for brain) were trans-acting. This could be due to higher effect size of trans versus cis-acting eQTL for mRNAs that are central to co-expression modules. That is, we may have enhanced the identification of trans-acting factors by looking for SNPs associated with mRNAs in co-expression modules that broadly influence gene expression patterns. Overall, these data indicate that the genomic mechanism driving physiological variation across environments is driven by trans-acting heart- or brain-specific mRNA expression.
Collapse
Affiliation(s)
- Melissa K Drown
- Department of Marine Biology and Ecology, The Rosenstiel School, University of Miami, Miami, FL, USA
| | - Marjorie F Oleksiak
- Department of Marine Biology and Ecology, The Rosenstiel School, University of Miami, Miami, FL, USA
| | - Douglas L Crawford
- Department of Marine Biology and Ecology, The Rosenstiel School, University of Miami, Miami, FL, USA
| |
Collapse
|