Genome-wide pQTL analysis of protein expression regulatory networks in the human liver

Identifying genetic variants that influence the abundance of cell states in single-cell data

Disease risk alleles influence the composition of cells present in the body, but modeling genetic effects on the cell states revealed by single-cell profiling is difficult because variant-associated states may reflect diverse combinations of the profiled cell features that are challenging to predefine. We introduce Genotype-Neighborhood Associations (GeNA), a statistical tool to identify cell-state abundance quantitative trait loci (csaQTLs) in high-dimensional single-cell datasets. Instead of testing associations to predefined cell states, GeNA flexibly identifies the cell states whose abundance is most associated with genetic variants. In a genome-wide survey of single-cell RNA sequencing peripheral blood profiling from 969 individuals, GeNA identifies five independent loci associated with shifts in the relative abundance of immune cell states. For example, rs3003-T (P = 1.96 × 10-11) associates with increased abundance of natural killer cells expressing tumor necrosis factor response programs. This csaQTL colocalizes with increased risk for psoriasis, an autoimmune disease that responds to anti-tumor necrosis factor treatments. Flexibly characterizing csaQTLs for granular cell states may help illuminate how genetic background alters cellular composition to confer disease risk.

Statistically and functionally fine-mapped blood eQTLs and pQTLs from 1,405 humans reveal distinct regulation patterns and disease relevance

Studying the genetic regulation of protein expression (through protein quantitative trait loci (pQTLs)) offers a deeper understanding of regulatory variants uncharacterized by mRNA expression regulation (expression QTLs (eQTLs)) studies. Here we report cis-eQTL and cis-pQTL statistical fine-mapping from 1,405 genotyped samples with blood mRNA and 2,932 plasma samples of protein expression, as part of the Japan COVID-19 Task Force (JCTF). Fine-mapped eQTLs (n = 3,464) were enriched for 932 variants validated with a massively parallel reporter assay. Fine-mapped pQTLs (n = 582) were enriched for missense variations on structured and extracellular domains, although the possibility of epitope-binding artifacts remains. Trans-eQTL and trans-pQTL analysis highlighted associations of class I HLA allele variation with KIR genes. We contrast the multi-tissue origin of plasma protein with blood mRNA, contributing to the limited colocalization level, distinct regulatory mechanisms and trait relevance of eQTLs and pQTLs. We report a negative correlation between ABO mRNA and protein expression because of linkage disequilibrium between distinct nearby eQTLs and pQTLs.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk postmortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA editing and alternative polyadenylation (APA) within a cell-type-specific population of human neural progenitors and neurons. More RNA editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting that genetically mediated post-transcriptional regulation during brain development leads to differences in brain function.

Multiomics Analyses Identify AKR1A1 as a Biomarker for Diabetic Kidney Disease

Diabetic kidney disease (DKD) is the leading cause of end-stage kidney disease. Because many genes associate with DKD, multiomics approaches were used to narrow the list of functional genes, gene products, and related pathways providing insights into the pathophysiological mechanisms of DKD. The Kidney Precision Medicine Project human kidney single-cell RNA-sequencing (scRNA-seq) data set and Mendeley Data on human kidney cortex biopsy proteomics were used. The R package Seurat was used to analyze scRNA-seq data and data from a subset of proximal tubule cells. PathfindR was applied for pathway analysis in cell type-specific differentially expressed genes and the R limma package was used to analyze differential protein expression in kidney cortex. A total of 790 differentially expressed genes were identified in proximal tubule cells, including 530 upregulated and 260 downregulated transcripts. Compared with differentially expressed proteins, 24 genes or proteins were in common. An integrated analysis combining protein quantitative trait loci, genome-wide association study hits (namely, estimated glomerular filtration rate), and a plasma metabolomics analysis was performed using baseline metabolites predictive of DKD progression in our longitudinal Diabetes Heart Study samples. The aldo-keto reductase family 1 member A1 gene (AKR1A1) was revealed as a potential molecular hub for DKD cellular dysfunction in several cross-linked pathways featured by deficiency of this enzyme.

ARTICLE HIGHLIGHTS

Promises and Challenges of populational Proteomics in Health and Disease

Advances in proteomic assay technologies have significantly increased coverage and throughput, enabling recent increases in the number of large-scale population-based proteomic studies of human plasma and serum. Improvements in multiplexed protein assays have facilitated the quantification of thousands of proteins over a large dynamic range, a key requirement for detecting the lowest-ranging, and potentially the most disease-relevant, blood-circulating proteins. In this perspective, we examine how populational proteomic datasets in conjunction with other concurrent omic measures can be leveraged to better understand the genomic and non-genomic correlates of the soluble proteome, constructing biomarker panels for disease prediction, among others. Mass spectrometry workflows are discussed as they are becoming increasingly competitive with affinity-based array platforms in terms of speed, cost, and proteome coverage due to advances in both instrumentation and workflows. Despite much success, there remain considerable challenges such as orthogonal validation and absolute quantification. We also highlight emergent challenges associated with study design, analytical considerations, and data integration as population-scale studies are run in batches and may involve longitudinal samples collated over many years. Lastly, we take a look at the future of what the nascent next-generation proteomic technologies might provide to the analysis of large sets of blood samples, as well as the difficulties in designing large-scale studies that will likely require participation from multiple and complex funding sources and where data sharing, study designs, and financing must be solved.

Genetic determinants of plasma protein levels in the Estonian population

The proteome holds great potential as an intermediate layer between the genome and phenome. Previous protein quantitative trait locus studies have focused mainly on describing the effects of common genetic variations on the proteome. Here, we assessed the impact of the common and rare genetic variations as well as the copy number variants (CNVs) on 326 plasma proteins measured in up to 500 individuals. We identified 184 cis and 94 trans signals for 157 protein traits, which were further fine-mapped to credible sets for 101 cis and 87 trans signals for 151 proteins. Rare genetic variation contributed to the levels of 7 proteins, with 5 cis and 14 trans associations. CNVs were associated with the levels of 11 proteins (7 cis and 5 trans), examples including a 3q12.1 deletion acting as a hub for multiple trans associations; and a CNV overlapping NAIP, a sensor component of the NAIP-NLRC4 inflammasome which is affecting pro-inflammatory cytokine interleukin 18 levels. In summary, this work presents a comprehensive resource of genetic variation affecting the plasma protein levels and provides the interpretation of identified effects.

Secreted Protein Profiling of Human Aortic Smooth Muscle Cells Identifies Vascular Disease Associations

BACKGROUND

Smooth muscle cells (SMCs), which make up the medial layer of arteries, are key cell types involved in cardiovascular disease, the leading cause of mortality and morbidity worldwide. In response to microenvironment alterations, SMCs dedifferentiate from a contractile to a synthetic phenotype characterized by an increased proliferation, migration, production of ECM (extracellular matrix) components, and decreased expression of SMC-specific contractile markers. These phenotypic changes result in vascular remodeling and contribute to the pathogenesis of cardiovascular disease, including coronary artery disease, stroke, hypertension, and aortic aneurysms. Here, we aim to identify the genetic variants that regulate ECM secretion in SMCs and predict the causal proteins associated with vascular disease-related loci identified in genome-wide association studies.

METHODS

Using human aortic SMCs from 123 multiancestry healthy heart transplant donors, we collected the serum-free media in which the cells were cultured for 24 hours and conducted liquid chromatography-tandem mass spectrometry-based proteomic analysis of the conditioned media.

RESULTS

We measured the abundance of 270 ECM and related proteins. Next, we performed protein quantitative trait locus mapping and identified 20 loci associated with secreted protein abundance in SMCs. We functionally annotated these loci using a colocalization approach. This approach prioritized the genetic variant rs6739323-A at the 2p22.3 locus, which is associated with lower expression of LTBP1 (latent-transforming growth factor beta-binding protein 1) in SMCs and atherosclerosis-prone areas of the aorta, and increased risk for SMC calcification. We found that LTBP1 expression is abundant in SMCs, and its expression at mRNA and protein levels was reduced in unstable and advanced atherosclerotic plaque lesions.

CONCLUSIONS

Our results unravel the SMC proteome signature associated with vascular disorders, which may help identify potential therapeutic targets to accelerate the pathway to translation.

Genomic insights into the comorbidity between type 2 diabetes and schizophrenia

Multimorbidity represents an increasingly important public health challenge with far-reaching implications for health management and policy. Mental health and metabolic diseases have a well-established epidemiological association. In this study, we investigate the genetic intersection between type 2 diabetes and schizophrenia. We use Mendelian randomization to examine potential causal relationships between the two conditions and related endophenotypes. We report no compelling evidence that type 2 diabetes genetic liability potentially causally influences schizophrenia risk and vice versa. Our findings show that increased body mass index (BMI) has a protective effect against schizophrenia, in contrast to the well-known risk-increasing effect of BMI on type 2 diabetes risk. We identify evidence of colocalization of association signals for these two conditions at 11 genomic loci, six of which have opposing directions of effect for type 2 diabetes and schizophrenia. To elucidate these colocalizing signals, we integrate multi-omics data from bulk and single-cell gene expression studies, along with functional information. We identify putative effector genes and find that they are enriched for homeostasis and lipid-related pathways. We also highlight drug repurposing opportunities including N-methyl-D-aspartate (NMDA) receptor antagonists. Our findings provide insights into shared biological mechanisms for type 2 diabetes and schizophrenia, highlighting common factors that influence the risk of the two conditions in opposite directions and shedding light on the complex nature of this comorbidity.

Novel Independent Trans- and Cis-Genetic Variants Associated with CYP2D6 Expression and Activity in Human Livers

Cytochrome P450 2D6 (CYP2D6) is a critical hepatic drug-metabolizing enzyme in humans, responsible for metabolizing approximately 20%-25% of commonly used medications such as codeine, desipramine, fluvoxamine, paroxetine, and tamoxifen. The CYP2D6 gene is highly polymorphic, resulting in substantial interindividual variability in its catalytic function and the pharmacokinetics and therapeutic outcomes of its substrate drugs. Although many functional CYP2D6 variants have been discovered and validated, a significant portion of the variability in the expression and activity of CYP2D6 remains unexplained. In this study, we performed a genome-wide association study (GWAS) to identify novel variants associated with CYP2D6 protein expression in individual human livers, followed by a conditional analysis to control for the effect of functional CYP2D6 star alleles. We also examined their impact on hepatic CYP2D6 activity. Genotyping on a genome-wide scale was achieved using the Illumina Multi-Ethnic Genotyping Array (MEGA). A data-independent acquisition (DIA)-based proteomics method was used to quantify CYP2D6 protein concentrations. CYP2D6 activity was determined by measuring the dextromethorphan O-demethylation in individual human liver s9 fractions. The GWAS identified 44 single nuclear polymorphisms (SNPs) that are significantly associated with CYP2D6 protein expressions with a P value threshold of 5.0 × 10-7 After the conditional analysis, five SNPs, including the cis-variants rs1807493 and rs1062753 and the trans-variants rs4073010, rs729559, and rs80274432, emerged as independent variants significantly correlated with hepatic CYP2D6 protein expressions. Notably, four of these SNPs, except for rs80274432, also exhibited a significant association with CYP2D6 activities in human livers, suggesting their potential as novel and independent cis- and trans-variants regulating CYP2D6. SIGNIFICANT STATEMENT: Using individual human livers, we identified four novel cis- and trans-pQTLs/aQTLs (protein quantitative trait loci/activity quantitative trait loci) of Cytochrome P450 2D6 (CYP2D6) that are independent from known functional CYP2D6 star alleles. This study connects the CYP2D6 gene expression and activity, enhancing our understanding of the genetic variants associated with CYP2D6 protein expression and activity, potentially advancing our insight into the interindividual variability in CYP2D6 substrate medication response.

Identifying genetic variants that influence the abundance of cell states in single-cell data

To understand genetic mechanisms driving disease, it is essential but difficult to map how risk alleles affect the composition of cells present in the body. Single-cell profiling quantifies granular information about tissues, but variant-associated cell states may reflect diverse combinations of the profiled cell features that are challenging to predefine. We introduce GeNA (Genotype-Neighborhood Associations), a statistical tool to identify cell state abundance quantitative trait loci (csaQTLs) in high-dimensional single-cell datasets. Instead of testing associations to predefined cell states, GeNA flexibly identifies the cell states whose abundance is most associated with genetic variants. In a genome-wide survey of scRNA-seq peripheral blood profiling from 969 individuals,1GeNA identifies five independent loci associated with shifts in the relative abundance of immune cell states. For example, rs3003-T (p=1.96×10-11) associates with increased abundance of NK cells expressing TNF-α response programs. This csaQTL colocalizes with increased risk for psoriasis, an autoimmune disease that responds to anti-TNF treatments. Flexibly characterizing csaQTLs for granular cell states may help illuminate how genetic background alters cellular composition to confer disease risk.

Secreted protein profiling of human aortic smooth muscle cells identifies vascular disease associations

Background

Smooth muscle cells (SMCs), which make up the medial layer of arteries, are key cell types involved in cardiovascular diseases (CVD), the leading cause of mortality and morbidity worldwide. In response to microenvironment alterations, SMCs dedifferentiate from a "contractile" to a "synthetic" phenotype characterized by an increased proliferation, migration, production of extracellular matrix (ECM) components, and decreased expression of SMC-specific contractile markers. These phenotypic changes result in vascular remodeling and contribute to the pathogenesis of CVD, including coronary artery disease (CAD), stroke, hypertension, and aortic aneurysms. Here, we aim to identify the genetic variants that regulate ECM secretion in SMCs and predict the causal proteins associated with vascular disease-related loci identified in genome-wide association studies (GWAS).

Methods

Using human aortic SMCs from 123 multi-ancestry healthy heart transplant donors, we collected the serum-free media in which the cells were cultured for 24 hours and conducted Liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based proteomic analysis of the conditioned media.

Results

We measured the abundance of 270 ECM and related proteins. Next, we performed protein quantitative trait locus mapping (pQTL) and identified 20 loci associated with secreted protein abundance in SMCs. We functionally annotated these loci using a colocalization approach. This approach prioritized the genetic variant rs6739323-A at the 2p22.3 locus, which is associated with lower expression of LTBP1 in SMCs and atherosclerosis-prone areas of the aorta, and increased risk for SMC calcification. We found that LTBP1 expression is abundant in SMCs, and its expression at mRNA and protein levels was reduced in unstable and advanced atherosclerotic plaque lesions.

Conclusions

Our results unravel the SMC proteome signature associated with vascular disorders, which may help identify potential therapeutic targets to accelerate the pathway to translation.

Genome-Wide Association Study for the Genetic Determinants of Thiopurine Methyltransferase Protein Expression in Human Livers and Racial Differences

INTRODUCTION

Polymorphisms in the Thiopurine S-Methyltransferase (TPMT) gene are associated with decreased TPMT activity, but little is known about their impact on TPMT protein expression in the liver. This project is to conduct a genome-wide association study (GWAS) to identify single nucleotide polymorphisms (SNPs) associated with altered TPMT protein expression in human livers and to determine if demographics affect hepatic TPMT protein expression.

METHODS

Human liver samples (n = 287) were genotyped using a whole genome genotyping panel and quantified for TPMT protein expression using a Data-Independent Acquisition proteomics approach.

RESULTS AND DISCUSSION

Thirty-one SNPs were found to be associated with differential expression of TPMT protein in the human livers. Subsequent analysis, conditioning on rs1142345, a SNP associated with the TPMT*3A and TPMT*3C alleles, showed no additional independent signals. Mean TPMT expression is significantly higher in wildtype donors compared to those carrying the known TPMT alleles, including TPMT*3A, TPMT*3C, and TPMT*24 (0.107 ± 0.028 vs. 0.052 ± 0.014 pmol/mg total protein, P = 2.2 × 10-16). After removing samples carrying the known TPMT variants, European ancestry donors exhibited significantly higher expression than African ancestry donors (0.109 ± 0.026 vs. 0.090 ± 0.041 pmol/mg total protein, P = 0.020).

CONCLUSION

The GWAS identified 31 SNPs associated with TPMT protein expression in human livers. Hepatic TPMT protein expression was significantly lower in subjects carrying the TPMT*3A, TPMT*3C, and TPMT*24 alleles compared to non-carriers. European ancestry was associated with significantly higher hepatic TPMT protein expression than African ancestry, independent of known TPMT variants.

Analysis of liver miRNA in Hu sheep with different residual feed intake

Feed efficiency (FE), an important economic trait in sheep production, is indirectly assessed by residual feed intake (RFI). However, RFI in sheep is varied, and the molecular processes that regulate RFI are unclear. It is thus vital to investigate the molecular mechanism of RFI to developing a feed-efficient sheep. The miRNA-sequencing (RNA-Seq) was utilized to investigate miRNAs in liver tissue of 6 out of 137 sheep with extreme RFI phenotypic values. In these animals, as a typical metric of FE, RFI was used to distinguish differentially expressed miRNAs (DE_miRNAs) between animals with high (n = 3) and low (n = 3) phenotypic values. A total of 247 miRNAs were discovered in sheep, with four differentially expressed miRNAs (DE_miRNAs) detected. Among these DE_miRNAs, three were found to be upregulated and one was downregulated in animals with low residual feed intake (Low_RFI) compared to those with high residual feed intake (High_RFI). The target genes of DE_miRNAs were primarily associated with metabolic processes and biosynthetic process regulation. Furthermore, they were also considerably enriched in the FE related to glycolysis, protein synthesis and degradation, and amino acid biosynthesis pathways. Six genes were identified by co-expression analysis of DE_miRNAs target with DE_mRNAs. These results provide a theoretical basis for us to understand the sheep liver miRNAs in RFI molecular regulation.

Genomic insights into the comorbidity between type 2 diabetes and schizophrenia

Multimorbidity represents an increasingly important public health challenge with far-reaching implications for health management and policy. Mental health and metabolic diseases have a well-established epidemiological association. In this study, we investigate the genetic intersection between type 2 diabetes and schizophrenia. We use Mendelian randomization to examine potential causal relationships between the two conditions and related endophenotypes. We report no compelling evidence that type 2 diabetes genetic liability potentially causally influences schizophrenia risk and vice versa. Our findings show that increased body mass index (BMI) has a protective effect against schizophrenia, in contrast to the well-known risk-increasing effect of BMI on type 2 diabetes risk. We identify evidence of colocalization of association signals for these two conditions at 11 genomic loci, six of which have opposing directions of effect for type 2 diabetes and schizophrenia. To elucidate these colocalizing signals, we integrate multi-omics data from bulk and single-cell gene expression studies, along with functional information. We identify high-confidence effector genes and find that they are enriched for homeostasis and lipid-related pathways. We also highlight drug repurposing opportunities including N-methyl-D-aspartate (NMDA) receptor antagonists. Our findings provide insights into shared biological mechanisms for type 2 diabetes and schizophrenia, highlighting common factors that influence the risk of the two conditions in opposite directions and shedding light on the complex nature of this comorbidity.

Circulating Copper and Liver Cancer

The association between circulating copper and the risk of liver cancer has been investigated by previous studies, while the findings were inconsistent. Thus, we aimed to evaluate the association between circulating copper and liver cancer by using meta-analysis and Mendelian randomization (MR). For meta-analysis, PubMed and Web of Science were searched to identify eligible studies published before April 4, 2022. Standardized mean difference (SMD) with 95% confidence interval (CI) in circulating copper level between liver cancer patients and controls were pooled. Furthermore, we selected genetic instruments for circulating copper from a genome-wide association study (GWAS) to conduct MR analysis. The summary statistics related to liver cancer were obtained from two large independent cohorts, UKBB and FinnGen, respectively. MR analysis was performed mainly by inverse-variance weighted (IVW) approach, followed by maximum-likelihood method as sensitivity analysis. In meta-analysis of eight studies, circulating copper was found to be higher in liver cancer patients (SMD: 1.65; 95% CI: 0.65 to 2.65) with high heterogeneity (I2 = 96.40%, P = 0.001). However, inconsistent findings were observed among subgroups with high evidence. In MR analysis, genetically predicted circulating copper was not significantly associated with the risk of liver cancer by IVW in UKBB (OR: 1.38; 95% CI: 0.72 to 2.65) and FinnGen (OR: 1.10; 95% CI: 0.69 to 1.73) separately, and the pooled results produced similar results (OR: 1.18, 95% CI: 0.81 to 1.72). Moreover, non-significant finding was confirmed by using maximum-likelihood method. There is no sufficient evidence to demonstrate that high levels of circulating copper increase the risks of liver cancer.

Oxidized Soybean Oil Evoked Hepatic Fatty Acid Metabolism Disturbance in Rats and their Offspring

The oxidation of fats and oils is an undisputed subject of science, given the effect of oxidized fats and oils on food quality and safety. This study aimed to determine whether maternal exposure to oxidized soybean oil (OSO) causes lipid metabolism disorders in the liver and whether this lipid metabolism disorder can be transmitted to offspring or even worsened. A total of 60 female Sprague-Dawley (SD) rats were divided randomly into four groups in this study. Treatment groups received a pure diet of OSO with a peroxide value of 200, 400, or 800 mEqO2/kg, while the control group received fresh soybean oil (FSO). As for our results, OSO affected serum biochemical parameters in the maternal generation (F0) and induced liver histopathology changes, inflammation, and oxidative stress. Moreover, the expression of genes related to the liver X receptor α (LXRα)─sterol regulatory element binding protein-1c (SREBP-1c) signaling pathway was changed. Similar trends were found in the livers of offspring on postnatal days 21 and 56. In conclusion, exposure to OSO during gestation and lactation can affect liver lipid synthesis. Additionally, it is detrimental to the development of the offspring's liver, affecting fatty acid metabolism and causing liver damage.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk post-mortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA-editing and alternative polyadenylation (APA), within a cell-type-specific population of human neural progenitors and neurons. More RNA-editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting genetically mediated post-transcriptional regulation during brain development lead to differences in brain function.

Genetic architecture of heart mitochondrial proteome influencing cardiac hypertrophy

Mitochondria play an important role in both normal heart function and disease etiology. We report analysis of common genetic variations contributing to mitochondrial and heart functions using an integrative proteomics approach in a panel of inbred mouse strains called the Hybrid Mouse Diversity Panel (HMDP). We performed a whole heart proteome study in the HMDP (72 strains, n=2-3 mice) and retrieved 848 mitochondrial proteins (quantified in ≥50 strains). High-resolution association mapping on their relative abundance levels revealed three trans-acting genetic loci on chromosomes (chr) 7, 13 and 17 that regulate distinct classes of mitochondrial proteins as well as cardiac hypertrophy. DAVID enrichment analyses of genes regulated by each of the loci revealed that the chr13 locus was highly enriched for complex-I proteins (24 proteins, P=2.2E-61), the chr17 locus for mitochondrial ribonucleoprotein complex (17 proteins, P=3.1E-25) and the chr7 locus for ubiquinone biosynthesis (3 proteins, P=6.9E-05). Follow-up high resolution regional mapping identified NDUFS4, LRPPRC and COQ7 as the candidate genes for chr13, chr17 and chr7 loci, respectively, and both experimental and statistical analyses supported their causal roles. Furthermore, a large cohort of Diversity Outbred mice was used to corroborate Lrpprc gene as a driver of mitochondrial DNA (mtDNA)-encoded gene regulation, and to show that the chr17 locus is specific to heart. Variations in all three loci were associated with heart mass in at least one of two independent heart stress models, namely, isoproterenol-induced heart failure and diet-induced obesity. These findings suggest that common variations in certain mitochondrial proteins can act in trans to influence tissue-specific mitochondrial functions and contribute to heart hypertrophy, elucidating mechanisms that may underlie genetic susceptibility to heart failure in human populations.

Functional characterization of human genomic variation linked to polygenic diseases

The burden of human disease lies predominantly in polygenic diseases. Since the early 2000s, genome-wide association studies (GWAS) have identified genetic variants and loci associated with complex traits. These have ranged from variants in coding sequences to mutations in regulatory regions, such as promoters and enhancers, as well as mutations affecting mediators of mRNA stability and other downstream regulators, such as 5' and 3'-untranslated regions (UTRs), long noncoding RNA (lncRNA), and miRNA. Recent research advances in genetics have utilized a combination of computational techniques, high-throughput in vitro and in vivo screening modalities, and precise genome editing to impute the function of diverse classes of genetic variants identified through GWAS. In this review, we highlight the vastness of genomic variants associated with polygenic disease risk and address recent advances in how genetic tools can be used to functionally characterize them.

Proteome and genome integration analysis of obesity

ABSTRACT

The prevalence of obesity has increased worldwide in recent decades. Genetic factors are now known to play a substantial role in the predisposition to obesity and may contribute up to 70% of the risk for obesity. Technological advancements during the last decades have allowed the identification of many hundreds of genetic markers associated with obesity. However, the transformation of current genetic variant-obesity associations into biological knowledge has been proven challenging. Genomics and proteomics are complementary fields, as proteomics extends functional analyses. Integrating genomic and proteomic data can help to bridge a gap in knowledge regarding genetic variant-obesity associations and to identify new drug targets for the treatment of obesity. We provide an overview of the published papers on the integrated analysis of proteomic and genomic data in obesity and summarize four mainstream strategies: overlap, colocalization, Mendelian randomization, and proteome-wide association studies. The integrated analyses identified many obesity-associated proteins, such as leptin, follistatin, and adenylate cyclase 3. Despite great progress, integrative studies focusing on obesity are still limited. There is an increased demand for large prospective cohort studies to identify and validate findings, and further apply these findings to the prevention, intervention, and treatment of obesity. In addition, we also discuss several other potential integration methods.

Mass spectrometry-based high-throughput proteomics and its role in biomedical studies and systems biology

There are multiple reasons why the next generation of biological and medical studies require increasing numbers of samples. Biological systems are dynamic, and the effect of a perturbation depends on the genetic background and environment. As a consequence, many conditions need to be considered to reach generalizable conclusions. Moreover, human population and clinical studies only reach sufficient statistical power if conducted at scale and with precise measurement methods. Finally, many proteins remain without sufficient functional annotations, because they have not been systematically studied under a broad range of conditions. In this review, we discuss the latest technical developments in mass spectrometry (MS)-based proteomics that facilitate large-scale studies by fast and efficient chromatography, fast scanning mass spectrometers, data-independent acquisition (DIA), and new software. We further highlight recent studies which demonstrate how high-throughput (HT) proteomics can be applied to capture biological diversity, to annotate gene functions or to generate predictive and prognostic models for human diseases.

Multi-Omics Studies in Historically Excluded Populations: The Road to Equity

Over the past few decades, genomewide association studies (GWASs) have identified the specific genetics variants contributing to many complex diseases by testing millions of genetic variations across the human genome against a variety of phenotypes. However, GWASs are limited in their ability to uncover mechanistic insight given that most significant associations are found in non-coding region of the genome. Furthermore, the lack of diversity in studies has stymied the advance of precision medicine for many historically excluded populations. In this review, we summarize most popular multi-omics approaches (genomics, transcriptomics, proteomics, and metabolomics) related to precision medicine and highlight if diverse populations have been included and how their findings have advance biological understanding of disease and drug response. New methods that incorporate local ancestry have been to improve the power of GWASs for admixed populations (such as African Americans and Latinx). Because most signals from GWAS are in the non-coding region, other machine learning and omics approaches have been developed to identify the potential causative single-nucleotide polymorphisms and genes that explain these phenotypes. These include polygenic risk scores, expression quantitative trait locus mapping, and transcriptome-wide association studies. Analogous protein methods, such as proteins quantitative trait locus mapping, proteome-wide association studies, and metabolomic approaches provide insight into the consequences of genetic variation on protein abundance. Whereas, integrated multi-omics studies have improved our understanding of the mechanisms for genetic association, we still lack the datasets and cohorts for historically excluded populations to provide equity in precision medicine and pharmacogenomics.

Exome-wide screening identifies novel rare risk variants for bone mineral density

Bone mineral density (BMD) is an independent risk factor of osteoporosis-related fractures. We performed gene-based burden tests to assess the association between rare variants and BMD, and identified several BMD candidate genes.

PURPOSE

BMD is highly heritable and a major predictor of osteoporotic fractures, but its genetic basis remains unclear. We aimed to identify rare risk variants contributing to BMD.

METHODS

Utilizing the newly released UK Biobank 200,643 exome dataset, we conducted a gene-based exome-wide association study in males and females, respectively. First, 100,639 males and 117,338 females with BMD values were included in the polygenic risk scores (PRS) analysis. Among individuals with lower 30% PRS, cases were individuals with top 10% BMD, and individuals with bottom 10% BMD were the controls. Considering the effects of vitamin D (VD), individuals with the highest 30% VD concentration were selected for VD-BMD analysis. After quality control, 741 males and 697 females were included in the BMD analysis, and 717 males and 708 females were included in the VD-BMD analysis. The variants were annotated by ANNOVAR software, then BMD and VD-BMD qualified variants were imported into the SKAT R-package to perform gene-based burden tests, respectively.

RESULTS

The gene-based burden test of the exonic variants identified genome-wide candidate associations in ANKRD18A (P = 1.60 × 10^-5, P_{Bonferroni adjust} = 2.11 × 10^-3), C22orf31 (P = 3.49 × 10^-4, P_{Bonferroni adjust} = 3.17 × 10^-2), and SPATC1L (P = 1.09 × 10^-5, P_{Bonferroni adjust} = 8.80 × 10^-3). For VD-BMD analysis, three genes were associated with BMD, such as NIPAL1 (P = 1.06 × 10^-3, P_{Bonferroni adjust} = 3.91 × 10^-2).

CONCLUSIONS

Our study suggested that rare variants contribute to BMD, providing new sights for broadening the genetic structure of BMD.

Integration of Candida albicans-induced single-cell gene expression data and secretory protein concentrations reveal genetic regulators of inflammation

Both gene expression and protein concentrations are regulated by genetic variants. Exploring the regulation of both eQTLs and pQTLs simultaneously in a context- and cell-type dependent manner may help to unravel mechanistic basis for genetic regulation of pQTLs. Here, we performed meta-analysis of Candida albicans-induced pQTLs from two population-based cohorts and intersected the results with Candida-induced cell-type specific expression association data (eQTL). This revealed systematic differences between the pQTLs and eQTL, where only 35% of the pQTLs significantly correlated with mRNA expressions at single cell level, indicating the limitation of eQTLs use as a proxy for pQTLs. By taking advantage of the tightly co-regulated pattern of the proteins, we also identified SNPs affecting protein network upon Candida stimulations. Colocalization of pQTLs and eQTLs signals implicated several genomic loci including MMP-1 and AMZ1. Analysis of Candida-induced single cell gene expression data implicated specific cell types that exhibit significant expression QTLs upon stimulation. By highlighting the role of trans-regulatory networks in determining the abundance of secretory proteins, our study serve as a framework to gain insights into the mechanisms of genetic regulation of protein levels in a context-dependent manner.

Integrative Post-Genome-Wide Association Study Analyses Relevant to Psychiatric Disorders: Imputing Transcriptome and Proteome Signals

Background

The genome-wide association study (GWAS) is a common tool to identify genetic variants associated with complex traits, including psychiatric disorders (PDs). However, post-GWAS analyses are needed to extend the statistical inference to biologically relevant entities, e.g., genes, proteins, and pathways. To achieve this goal, researchers developed methods that incorporate biologically relevant intermediate molecular phenotypes, such as gene expression and protein abundance, which are posited to mediate the variant-trait association. Transcriptome-wide association study (TWAS) and proteome-wide association study (PWAS) are commonly used methods to test the association between these molecular mediators and the trait.

Summary

In this review, we discuss the most recent developments in TWAS and PWAS. These methods integrate existing "omic" information with the GWAS summary statistics for trait(s) of interest. Specifically, they impute transcript/protein data and test the association between imputed gene expression/protein level with phenotype of interest by using (i) GWAS summary statistics and (ii) reference transcriptomic/proteomic/genomic datasets. TWAS and PWAS are suitable as analysis tools for (i) primary association scan and (ii) fine-mapping to identify potentially causal genes for PDs.

Key Messages

As post-GWAS analyses, TWAS and PWAS have the potential to highlight causal genes for PDs. These prioritized genes could indicate targets for the development of novel drug therapies. For researchers attempting such analyses, we recommend Mendelian randomization tools that use GWAS statistics for both trait and reference datasets, e.g., summary Mendelian randomization (SMR). We base our recommendation on (i) being able to use the same tool for both TWAS and PWAS, (ii) not requiring the pre-computed weights (and thus easier to update for larger reference datasets), and (iii) most larger transcriptome reference datasets are publicly available and easy to transform into a compatible format for SMR analysis.

Transcriptomic and proteomic retinal pigment epithelium signatures of age-related macular degeneration

There are currently no treatments for geographic atrophy, the advanced form of age-related macular degeneration. Hence, innovative studies are needed to model this condition and prevent or delay its progression. Induced pluripotent stem cells generated from patients with geographic atrophy and healthy individuals were differentiated to retinal pigment epithelium. Integrating transcriptional profiles of 127,659 retinal pigment epithelium cells generated from 43 individuals with geographic atrophy and 36 controls with genotype data, we identify 445 expression quantitative trait loci in cis that are asssociated with disease status and specific to retinal pigment epithelium subpopulations. Transcriptomics and proteomics approaches identify molecular pathways significantly upregulated in geographic atrophy, including in mitochondrial functions, metabolic pathways and extracellular cellular matrix reorganization. Five significant protein quantitative trait loci that regulate protein expression in the retinal pigment epithelium and in geographic atrophy are identified - two of which share variants with cis- expression quantitative trait loci, including proteins involved in mitochondrial biology and neurodegeneration. Investigation of mitochondrial metabolism confirms mitochondrial dysfunction as a core constitutive difference of the retinal pigment epithelium from patients with geographic atrophy. This study uncovers important differences in retinal pigment epithelium homeostasis associated with geographic atrophy.

Heritability Estimation of Multiple Sclerosis Related Plasma Protein Levels in Sardinian Families with Immunochip Genotyping Data

This work aimed at estimating narrow-sense heritability, defined as the proportion of the phenotypic variance explained by the sum of additive genetic effects, via Haseman–Elston regression for a subset of 56 plasma protein levels related to Multiple Sclerosis (MS). These were measured in 212 related individuals (with 69 MS cases and 143 healthy controls) obtained from 20 Sardinian families with MS history. Using pedigree information, we found seven statistically significant heritable plasma protein levels (after multiple testing correction), i.e., Gc (h2 = 0.77; 95%CI: 0.36, 1.00), Plat (h2 = 0.70; 95%CI: 0.27, 0.95), Anxa1 (h2 = 0.68; 95%CI: 0.27, 1.00), Sod1 (h2 = 0.58; 95%CI: 0.18, 0.96), Irf8 (h2 = 0.56; 95%CI: 0.19, 0.99), Ptger4 (h2 = 0.45; 95%CI: 0.10, 0.96), and Fadd (h2 = 0.41; 95%CI: 0.06, 0.84). A subsequent analysis was performed on these statistically significant heritable plasma protein levels employing Immunochip genotyping data obtained in 155 healthy controls (92 related and 63 unrelated); we found a meaningful proportion of heritable plasma protein levels’ variability explained by a small set of SNPs. Overall, the results obtained, for these seven MS-related proteins, emphasized a high additive genetic variance component explaining plasma levels’ variability.

Mendelian randomization study for the roles of IL-18 and IL-1 receptor antagonist in the development of inflammatory bowel disease

BACKGROUND AND AIMS

IL-1 and IL-18 play important roles in intestine barrier integrity maintenance and inflammatory response. However, their net effects on the risk of IBD are still inconclusive. Here, we used Mendelian randomization (MR) approaches to investigate the causal associations of IL-18 and IL-1Ra (receptor antagonist) on the risks of IBD and subtypes.

METHODS

For IL-18, both three-sample and two-sample MR approaches were used for the causal inferences. In three-sample MR, three single nucleotide polymorphisms (SNPs) and the effect values were extracted from two quantitative trait loci (pQTL) datasets with non-overlapping populations. In two-sample MR, we extracted genetic instruments information from the same larger pQTL dataset. For IL-1Ra, we applied the two-sample MR method with summary-statistics from the larger pQTL dataset. Summary-level results of three large IBD/CD/UC genome-wide association studies in European ancestry were employed. Inverse-variance weighted method, various sensitivity analyses and meta-analysis were performed to give causal estimates, detect heterogeneity and correct for outliers.

RESULTS

We observed consistent positive causal effects of IL-18 on all three major outcomes using three-sample MR, with meta-analyses odds ratios (ORs) equal to 1.240 (IBD), 1.199 (CD) and 1.274 (UC) respectively. The two-sample MR demonstrated similar results. Moreover, genetically predicted IL-1Ra is inversely associated with the risk of IBD/UC/CD with ORs equal to 0.915 (IBD), 0.902 (CD) and 0.899 (UC) respectively in meta-analyses.

CONCLUSIONS

This study suggested genetically predicted IL-18 and IL-1Ra level are causally associated with an increased and decreased risk of IBD and subtypes.

Genetically regulated gene expression and proteins revealed discordant effects

Background

Although gene-expression (GE) and protein levels are typically strongly genetically regulated, their correlation is known to be low. Here we investigate this phenomenon by focusing on the genetic background of this correlation in order to understand the similarities and differences in the genetic regulation of these omics layers.

Methods and results

We performed locus-wide association studies of 92 protein levels measured in whole blood for 2,014 samples of European ancestry and found that 66 are genetically regulated. Three female- and one male-specific effects were detected. We estimated the genetically regulated GE for all significant genes in 49 GTEx v8 tissues. A total of 7 proteins showed negative correlations with their respective GE across multiple tissues. Finally, we tested for causal links of GE on protein expression via Mendelian Randomization, and confirmed a negative causal effect of GE on protein level for five of these genes in a total of 63 gene-tissue pairs: BLMH, CASP3, CXCL16, IL6R, and SFTPD. For IL6R, we replicated the negative causal effect on coronary-artery disease (CAD), while its GE was positively linked to CAD.

Conclusion

While total GE and protein levels are only weakly correlated, we found high correlations between their genetically regulated components across multiple tissues. Of note, strong negative causal effects of tissue-specific GE on five protein levels were detected. Causal network analyses revealed that GE effects on CAD risks was in general mediated by protein levels.

Detection of common copy number of variation underlying selection pressure in Middle Eastern horse breeds using whole-genome sequence data

Dareshouri, Arabian, and Akhal-Teke are three Middle Eastern horse breeds that have been selected for endurance and adaptation to harsh climates. Deciphering the genetic characteristics of these horses by tracing selection footprints and copy number of variations will be helpful in improving our understanding of equine breeds' development and adaptation. For this purpose, we sequenced the whole-genome of four Dareshouri horses using Illumina Hiseq panels and compared them with publicly available whole-genome sequences of Arabian (n=3) and Akhal-Teke (n=3) horses . Three tests of FLK, hapFLK, and pooled heterozygosity were applied using a sliding window (window size=100kb, step size=50kb) approach to detect putative selection signals. Copy number variation analysis was applied to investigate copy number of variants (CNVs), and the results were used to suggest selection signatures involving CNVs. Whole-genome sequencing demonstrated 8,837,950 single nucleotide polymorphisms (SNPs) in autosomal chromosomes. We suggested 58 genes and three quantitative trait loci (QTLs), including some related to horse gait, insect bite hypersensitivity, and withers height, based on selective signals detected by adjusted p-value of Mahalanobis distance based on the rank-based P-values (Md-rank-P) method. We proposed 12 genomic regions under selection pressure involving CNVs which were previously reported to be associated with metabolism energy (SLC5A8), champagne dilution in horses (SLC36A1), and synthesis of polyunsaturated fatty acids (FAT2). Only 10 Middle Eastern horses were tested in this study; therefore, the conclusions are speculative. Our findings are useful to better understanding the evolution and adaptation of Middle Eastern horse breeds.

Plasma proteome analyses in individuals of European and African ancestry identify cis-pQTLs and models for proteome-wide association studies

Improved understanding of genetic regulation of the proteome can facilitate identification of the causal mechanisms for complex traits. We analyzed data on 4,657 plasma proteins from 7,213 European American (EA) and 1,871 African American (AA) individuals from the Atherosclerosis Risk in Communities study, and further replicated findings on 467 AA individuals from the African American Study of Kidney Disease and Hypertension study. Here, we identified 2,004 proteins in EA and 1,618 in AA, with most overlapping, which showed associations with common variants in cis-regions. Availability of AA samples led to smaller credible sets and notable number of population-specific cis-protein quantitative trait loci. Elastic Net produced powerful models for protein prediction in both populations. An application of proteome-wide association studies to serum urate and gout implicated several proteins, including IL1RN, revealing the promise of the drug anakinra to treat acute gout flares. Our study demonstrates the value of large and diverse ancestry study to investigate the genetic mechanisms of molecular phenotypes and their relationship with complex traits.

Functional Characterization of Genetic Variant Effects on Expression

Thousands of common genetic variants in the human population have been associated with disease risk and phenotypic variation by genome-wide association studies (GWAS). However, the majority of GWAS variants fall into noncoding regions of the genome, complicating our understanding of their regulatory functions, and few molecular mechanisms of GWAS variant effects have been clearly elucidated. Here, we set out to review genetic variant effects, focusing on expression quantitative trait loci (eQTLs), including their utility in interpreting GWAS variant mechanisms. We discuss the interrelated challenges and opportunities for eQTL analysis, covering determining causal variants, elucidating molecular mechanisms of action, and understanding context variability. Addressing these questions can enable better functional characterization of disease-associated loci and provide insights into fundamental biological questions of the noncoding genetic regulatory code and its control of gene expression. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Multi-omics analysis reveals the genetic basis of rice fragrance mediated by betaine aldehyde dehydrogenase 2

INTRODUCTION

Fragrance is an important economic and quality trait in rice. The trait is controlled by the recessive gene betaine aldehyde dehydrogenase 2 (BADH2) via the production of 2-acetyl-1-pyrroline (2AP).

OBJECTIVES

Variation in BADH2 was evaluated at the population, genetic, transcriptional, and metabolic levels to obtain insights into fragrance regulation in rice.

METHODS

Whole-genome resequencing of the Korean World Rice Collection of 475 rice accessions, including 421 breeding lines and 54 wild accessions, was performed. Transcriptome analyses of a subset of 279 accessions, proteome analyses of 64 accessions, and volatile profiling of 421 breeding lines were also performed.

RESULTS

We identified over 3.1 million high-quality single nucleotide polymorphisms (SNPs) in Korean rice collection. Most SNPs were present in intergenic regions (79%), and 190,148 SNPs (6%) were located in the coding sequence, of which 53% were nonsynonymous. In total, 38 haplotypes were identified in the BADH2 coding region, including four novel haplotypes (one in cultivated and three in wild accessions). Tajima's D values suggested that BADH2 was under balancing selection in japonica rice. Furthermore, we identified 316 expression quantitative trait loci (eQTL), including 185 cis-eQTLs and 131 trans-eQTLs, involved in BADH2 regulation. A protein quantitative trait loci (pQTL) analysis revealed the presence of trans-pQTLs; 13 pQTLs were mapped 1 Mbp from the BADH2 region. Based on variable importance in projection (VIP) scores, 15 volatile compounds, including 2AP, discriminated haplotypes and were potential biomarkers for rice fragrance.

CONCLUSION

We generated a catalog of haplotypes based on a resequencing analysis of a large number of rice accessions. eQTLs and pQTLs associated with BADH2 gene expression and protein accumulation are likely involved in the regulation of 2AP variation in fragrant rice. These data improve our understanding of fragrance and provide valuable information for rice breeding.

Pharmaco-proteogenomic profiling of pediatric diffuse midline glioma to inform future treatment strategies

Diffuse midline glioma (DMG) is a deadly pediatric and adolescent central nervous system (CNS) tumor localized along the midline structures of the brain atop the spinal cord. With a median overall survival (OS) of just 9–11-months, DMG is characterized by global hypomethylation of histone H3 at lysine 27 (H3K27me3), driven by recurring somatic mutations in H3 genes including, HIST1H3B/C (H3.1K27M) or H3F3A (H3.3K27M), or through overexpression of EZHIP in patients harboring wildtype H3. The recent World Health Organization’s 5th Classification of CNS Tumors now designates DMG as, ‘H3 K27-altered’, suggesting that global H3K27me3 hypomethylation is a ubiquitous feature of DMG and drives devastating transcriptional programs for which there are no treatments. H3-alterations co-segregate with various other somatic driver mutations, highlighting the high-level of intertumoral heterogeneity of DMG. Furthermore, DMG is also characterized by very high-level intratumoral diversity with tumors harboring multiple subclones within each primary tumor. Each subclone contains their own combinations of driver and passenger lesions that continually evolve, making precision-based medicine challenging to successful execute. Whilst the intertumoral heterogeneity of DMG has been extensively investigated, this is yet to translate to an increase in patient survival. Conversely, our understanding of the non-genomic factors that drive the rapid growth and fatal nature of DMG, including endogenous and exogenous microenvironmental influences, neurological cues, and the posttranscriptional and posttranslational architecture of DMG remains enigmatic or at best, immature. However, these factors are likely to play a significant role in the complex biological sequelae that drives the disease. Here we summarize the heterogeneity of DMG and emphasize how analysis of the posttranslational architecture may improve treatment paradigms. We describe factors that contribute to treatment response and disease progression, as well as highlight the potential for pharmaco-proteogenomics (i.e., the integration of genomics, proteomics and pharmacology) in the management of this uniformly fatal cancer.

A trans-omic Mendelian randomization study of parental lifespan uncovers novel aging biology and therapeutic candidates for chronic diseases

The study of parental lifespan has emerged as an innovative tool to advance aging biology and our understanding of the genetic architecture of human longevity and aging-associated diseases. Here, we leveraged summary statistics of a genome-wide association study including over one million parental lifespans to identify genetically regulated genes from the Genotype-Tissue Expression project. Through a combination of multi-tissue transcriptome-wide association analyses and genetic colocalization, we identified novel genes that may be associated with parental lifespan. Mendelian randomization (MR) analyses also identified circulating proteins and metabolites causally associated with parental lifespan and chronic diseases offering new drug repositioning opportunities such as those targeting apolipoprotein-B-containing lipoproteins. Liver expression of HP, the gene encoding haptoglobin, and plasma haptoglobin levels were causally linked with parental lifespan. Phenome-wide MR analyses were used to map genetically regulated genes, proteins and metabolites with other human traits as well as the disease-related phenome in the FinnGen cohorts (n = 135,638). Altogether, this study identified new candidate genes, circulating proteins and metabolites that may influence human aging as well as potential therapeutic targets for chronic diseases that warrant further investigation.

Regulation of protein abundance in genetically diverse mouse populations

Genetically diverse mouse populations are powerful tools for characterizing the regulation of the proteome and its relationship to whole-organism phenotypes. We used mass spectrometry to profile and quantify the abundance of 6,798 proteins in liver tissue from mice of both sexes across 58 Collaborative Cross (CC) inbred strains. We previously collected liver proteomics data from the related Diversity Outbred (DO) mice and their founder strains. We show concordance across the proteomics datasets despite being generated from separate experiments, allowing comparative analysis. We map protein abundance quantitative trait loci (pQTLs), identifying 1,087 local and 285 distal in the CC mice and 1,706 local and 414 distal in the DO mice. We find that regulatory effects on individual proteins are conserved across the mouse populations, in particular for local genetic variation and sex differences. In comparison, proteins that form complexes are often co-regulated, displaying varying genetic architectures, and overall show lower heritability and map fewer pQTLs. We have made this resource publicly available to enable quantitative analyses of the regulation of the proteome.

MRLocus: Identifying causal genes mediating a trait through Bayesian estimation of allelic heterogeneity

Expression quantitative trait loci (eQTL) studies are used to understand the regulatory function of non-coding genome-wide association study (GWAS) risk loci, but colocalization alone does not demonstrate a causal relationship of gene expression affecting a trait. Evidence for mediation, that perturbation of gene expression in a given tissue or developmental context will induce a change in the downstream GWAS trait, can be provided by two-sample Mendelian Randomization (MR). Here, we introduce a new statistical method, MRLocus, for Bayesian estimation of the gene-to-trait effect from eQTL and GWAS summary data for loci with evidence of allelic heterogeneity, that is, containing multiple causal variants. MRLocus makes use of a colocalization step applied to each nearly-LD-independent eQTL, followed by an MR analysis step across eQTLs. Additionally, our method involves estimation of the extent of allelic heterogeneity through a dispersion parameter, indicating variable mediation effects from each individual eQTL on the downstream trait. Our method is evaluated against other state-of-the-art methods for estimation of the gene-to-trait mediation effect, using an existing simulation framework. In simulation, MRLocus often has the highest accuracy among competing methods, and in each case provides more accurate estimation of uncertainty as assessed through interval coverage. MRLocus is then applied to five candidate causal genes for mediation of particular GWAS traits, where gene-to-trait effects are concordant with those previously reported. We find that MRLocus's estimation of the causal effect across eQTLs within a locus provides useful information for determining how perturbation of gene expression or individual regulatory elements will affect downstream traits. The MRLocus method is implemented as an R package available at https://mikelove.github.io/mrlocus.

Epigenetic Element-Based Transcriptome-Wide Association Study Identifies Novel Genes for Bipolar Disorder

Since the bipolar disorder (BD) signals identified by genome-wide association study (GWAS) often reside in the non-coding regions, understanding the biological relevance of these genetic loci has proven to be complicated. Transcriptome-wide association studies (TWAS) providing a powerful approach to identify novel disease risk genes and uncover possible causal genes at loci identified previously by GWAS. However, these methods did not consider the importance of epigenetic regulation in gene expression. Here, we developed a novel epigenetic element-based transcriptome-wide association study (ETWAS) that tested the effects of genetic variants on gene expression levels with the epigenetic features as prior and further mediated the association between predicted expression and BD. We conducted an ETWAS consisting of 20 352 cases and 31 358 controls and identified 44 transcriptome-wide significant hits. We found 14 conditionally independent genes, and 10 genes that did not previously implicate with BD were regarded as novel candidate genes, such as ASB16 in the cerebellar hemisphere (P = 9.29 × 10-8). We demonstrated that several genome-wide significant signals from the BD GWAS driven by genetically regulated expression, and NEK4 explained 90.1% of the GWAS signal. Additionally, ETWAS identified genes could explain heritability beyond that explained by GWAS-associated SNPs (P = 5.60 × 10-66). By querying the SNPs in the final models of identified genes in phenome databases, we identified several phenotypes previously associated with BD, such as schizophrenia and depression. In conclusion, ETWAS is a powerful method, and we identified several novel candidate genes associated with BD.

Simultaneous quantification of mRNA and protein in single cells reveals post-transcriptional effects of genetic variation

Trans-acting DNA variants may specifically affect mRNA or protein levels of genes located throughout the genome. However, prior work compared trans-acting loci mapped in separate studies, many of which had limited statistical power. Here, we developed a CRISPR-based system for simultaneous quantification of mRNA and protein of a given gene via dual fluorescent reporters in single, live cells of the yeast Saccharomyces cerevisiae. In large populations of recombinant cells from a cross between two genetically divergent strains, we mapped 86 trans-acting loci affecting the expression of ten genes. Less than 20% of these loci had concordant effects on mRNA and protein of the same gene. Most loci influenced protein but not mRNA of a given gene. One locus harbored a premature stop variant in the YAK1 kinase gene that had specific effects on protein or mRNA of dozens of genes. These results demonstrate complex, post-transcriptional genetic effects on gene expression.

Genome-wide pQTL analysis of protein expression regulatory networks in the human liver

BACKGROUND

Previous expression quantitative trait loci (eQTL) studies have identified thousands of genetic variants to be associated with gene expression at the mRNA level in the human liver. However, protein expression often correlates poorly with mRNA levels. Thus, protein quantitative trait loci (pQTL) study is required to identify genetic variants that regulate protein expression in human livers.

RESULTS

We conducted a genome-wide pQTL study in 287 normal human liver samples and identified 900 local pQTL variants and 4026 distant pQTL variants. We further discovered 53 genome hotspots of pQTL variants. Transcriptional region mapping analysis showed that 1133 pQTL variants are in transcriptional regulatory regions. Genomic region enrichment analysis of the identified pQTL variants revealed 804 potential regulatory interactions among 595 predicted regulators (e.g., non-coding RNAs) and 394 proteins. Moreover, pQTL variants and trait-variant integration analysis implied several novel mechanisms underlying the relationships between protein expression and liver diseases, such as alcohol dependence. Notably, over 2000 of the identified pQTL variants have not been reported in previous eQTL studies, suggesting extensive involvement of genetic polymorphisms in post-transcriptional regulation of protein expression in human livers.

CONCLUSIONS

We have partially established protein expression regulation networks in human livers and generated a wealth of pQTL data that could serve as a valuable resource for the scientific community.