The EN-TEx resource of multi-tissue personal epigenomes & variant-impact models.

Understanding how genetic variants impact molecular phenotypes is a key goal of functional genomics, currently hindered by reliance on a single haploid reference genome. Here, we present the EN-TEx resource of 1,635 open-access datasets from four donors (∼30 tissues × âˆ¼15 assays). The datasets are mapped to matched, diploid genomes with long-read phasing and structural variants, instantiating a catalog of >1 million allele-specific loci. These loci exhibit coordinated activity along haplotypes and are less conserved than corresponding, non-allele-specific ones. Surprisingly, a deep-learning transformer model can predict the allele-specific activity based only on local nucleotide-sequence context, highlighting the importance of transcription-factor-binding motifs particularly sensitive to variants. Furthermore, combining EN-TEx with existing genome annotations reveals strong associations between allele-specific and GWAS loci. It also enables models for transferring known eQTLs to difficult-to-profile tissues (e.g., from skin to heart). Overall, EN-TEx provides rich data and generalizable models for more accurate personal functional genomics.

Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease.

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

A pan-cancer transcriptome analysis of exitron splicing identifies novel cancer driver genes and neoepitopes.

Exitron splicing (EIS) creates a cryptic intron (called an exitron) within a protein-coding exon to increase proteome diversity. EIS is poorly characterized, but emerging evidence suggests a role for EIS in cancer. Through a systematic investigation of EIS across 33 cancers from 9,599 tumor transcriptomes, we discovered that EIS affected 63% of human coding genes and that 95% of those events were tumor specific. Notably, we observed a mutually exclusive pattern between EIS and somatic mutations in their affected genes. Functionally, we discovered that EIS altered known and novel cancer driver genes for causing gain- or loss-of-function, which promotes tumor progression. Importantly, we identified EIS-derived neoepitopes that bind to major histocompatibility complex (MHC) class I or II. Analysis of clinical data from a clear cell renal cell carcinoma cohort revealed an association between EIS-derived neoantigen load and checkpoint inhibitor response. Our findings establish the importance of considering EIS alterations when nominating cancer driver events and neoantigens.

ATG7(2) Interacts With Metabolic Proteins and Regulates Central Energy Metabolism.

Macroautophagy/autophagy is an essential catabolic process that targets a wide variety of cellular components including proteins, organelles, and pathogens. ATG7, a protein involved in the autophagy process, plays a crucial role in maintaining cellular homeostasis and can contribute to the development of diseases such as cancer. ATG7 initiates autophagy by facilitating the lipidation of the ATG8 proteins in the growing autophagosome membrane. The noncanonical isoform ATG7(2) is unable to perform ATG8 lipidation; however, its cellular regulation and function are unknown. Here, we uncovered a distinct regulation and function of ATG7(2) in contrast with ATG7(1), the canonical isoform. First, affinity-purification mass spectrometry analysis revealed that ATG7(2) establishes direct protein-protein interactions (PPIs) with metabolic proteins, whereas ATG7(1) primarily interacts with autophagy machinery proteins. Furthermore, we identified that ATG7(2) mediates a decrease in metabolic activity, highlighting a novel splice-dependent function of this important autophagy protein. Then, we found a divergent expression pattern of ATG7(1) and ATG7(2) across human tissues. Conclusively, our work uncovers the divergent patterns of expression, protein interactions, and function of ATG7(2) in contrast to ATG7(1). These findings suggest a molecular switch between main catabolic processes through isoform-dependent expression of a key autophagy gene.

A Massive Proteogenomic Screen Identifies Thousands of Novel Peptides From the Human "Dark" Proteome.

Although the human gene annotation has been continuously improved over the past 2 decades, numerous studies demonstrated the existence of a "dark proteome", consisting of proteins that were critical for biological processes but not included in widely used gene catalogs. The Genotype-Tissue Expression project generated more than 15,000 RNA-seq datasets from multiple tissues, which modeled 30 million transcripts in the human genome. To provide a resource of high-confidence novel proteins from the dark proteome, we screened 50,000 mass spectrometry runs from over 900 projects to identify proteins translated from the Genotype-Tissue Expression transcript model with proteomic support. We also integrated 3.8 million common genetic variants from the gnomAD database to improve peptide identification. As a result, we identified 170,529 novel peptides with proteomic evidence, of which 6048 passed the strictest standard we defined and were supported by PepQuery. We provided a user-friendly website (https://ncorf.genes.fun/) for researchers to check the evidence of novel peptides from their studies. The findings will improve our understanding of coding genes and facilitate genomic data interpretation in biomedical research.

Role of autophagy-related genes in liver cancer prognosis.

Autophagy, a highly conserved process of protein and organelle degradation, has emerged as a critical regulator in various diseases, including cancer progression. In the context of liver cancer, the predictive value of autophagy-related genes remains ambiguous. Leveraging chip datasets from the TCGA and GTEx databases, we identified 23 differentially expressed autophagy-related genes in liver cancer. Notably, five key autophagy genes, PRKAA2, BIRC5, MAPT, IGF1, and SPNS1, were highlighted as potential prognostic markers, with MAPT showing significant overexpression in clinical samples. In vitro cellular assays further demonstrated that MAPT promotes liver cancer cell proliferation, migration, and invasion by inhibiting autophagy and suppressing apoptosis. Subsequent in vivo studies further corroborated the pro-tumorigenic role of MAPT by suppressing autophagy. Collectively, our model based on the five key genes provides a promising tool for predicting liver cancer prognosis, with MAPT emerging as a pivotal factor in tumor progression through autophagy modulation.

Patterns of Unwanted Biological and Technical Expression Variation Among 49 Human Tissues.

Tissue gene expression studies are impacted by biological and technical sources of variation, which can be broadly classified into wanted and unwanted variation. The latter, if not addressed, results in misleading biological conclusions. Methods have been proposed to reduce unwanted variation, such as normalization and batch correction. A more accurate understanding of all causes of variation could significantly improve the ability of these methods to remove unwanted variation while retaining variation corresponding to the biological question of interest. We used 17,282 samples from 49 human tissues in the Genotype-Tissue Expression data set (v8) to investigate patterns and causes of expression variation. Transcript expression was transformed to z-scores, and only the most variable 2% of transcripts were evaluated and clustered based on coexpression patterns. Clustered gene sets were assigned to different biological or technical causes based on histologic appearances and metadata elements. We identified 522 variable transcript clusters (median: 11 per tissue) among the samples. Of these, 63% were confidently explained, 16% were likely explained, 7% were low confidence explanations, and 14% had no clear cause. Histologic analysis annotated 46 clusters. Other common causes of variability included sex, sequencing contamination, immunoglobulin diversity, and compositional tissue differences. Less common biological causes included death interval (Hardy score), disease status, and age. Technical causes included blood draw timing and harvesting differences. Many of the causes of variation in bulk tissue expression were identifiable in the Tabula Sapiens data set of single-cell expression. This is among the largest explorations of the underlying sources of tissue expression variation. It uncovered expected and unexpected causes of variable gene expression and demonstrated the utility of matched histologic specimens. It further demonstrated the value of acquiring meaningful tissue harvesting metadata elements to use for improved normalization, batch correction, and analysis of both bulk and single-cell RNA-seq data.

Where Are the Disease-Associated eQTLs?

Most disease-associated variants, although located in putatively regulatory regions, do not have detectable effects on gene expression. One explanation could be that we have not examined gene expression in the cell types or conditions that are most relevant for disease. Even large-scale efforts to study gene expression across tissues are limited to human samples obtained opportunistically or postmortem, mostly from adults. In this review we evaluate recent findings and suggest an alternative strategy, drawing on the dynamic and highly context-specific nature of gene regulation. We discuss new technologies that can extend the standard regulatory mapping framework to more diverse, disease-relevant cell types and states.

Balancing the transcriptome: leveraging sample similarity to improve measures of gene specificity.

The spatial and temporal domain of a gene's expression can range from ubiquitous to highly specific. Quantifying the degree to which this expression is unique to a specific tissue or developmental timepoint can provide insight into the etiology of genetic diseases. However, quantifying specificity remains challenging as measures of specificity are sensitive to similarity between samples in the sample set. For example, in the Gene-Tissue Expression project (GTEx), brain subregions are overrepresented at 13 of 54 (24%) unique tissues sampled. In this dataset, existing specificity measures have a decreased ability to identify genes specific to the brain relative to other organs. To solve this problem, we leverage sample similarity information to weight samples such that overrepresented tissues do not have an outsized effect on specificity estimates. We test this reweighting procedure on 4 measures of specificity, Z-score, Tau, Tsi and Gini, in the GTEx data and in single cell datasets for zebrafish and mouse. For all of these measures, incorporating sample similarity information to weight samples results in greater stability of sets of genes called as specific and decreases the overall variance in the change of specificity estimates as sample sets become more unbalanced. Furthermore, the genes with the largest improvement in their specificity estimate's stability are those with functions related to the overrepresented sample types. Our results demonstrate that incorporating similarity information improves specificity estimates' stability to the choice of the sample set used to define the transcriptome, providing more robust and reproducible measures of specificity for downstream analyses.

Effect of ZIC2 on immune infiltration and ceRNA axis regulation in lung adenocarcinoma via bioinformatics and experimental studies.

OBJECTIVE: This study aimed to conclude the effect and mechanism of ZIC2 on immune infiltration in lung adenocarcinoma (LUAD). METHODS: Expression of ZIC2 in several kinds of normal tissues of TCGA data was analyzed and its correlation with the baseline characteristic of LUAD patients were analyzed. The immune infiltration analysis of LUAD patients was performed by CIBERSORT algorithm. The correlation analysis between ZIC2 and immune cell composition was performed. Additionally, the potential upstream regulatory mechanisms of ZIC2 were predicted to identify the possible miRNAs and lncRNAs that regulated ZIC2 in LUAD. In vitro and in vivo experiments were also conducted to confirm the potential effect of ZIC2 on cell proliferation and invasion ability of LUAD cells. RESULTS: ZIC2 expression was decreased in various normal tissues, but increased in multiple tumors, including LUAD, and correlated with the prognosis of LUAD patients. Enrichment by GO and KEGG suggested the possible association of ZIC2 with cell cycle and p53 signal pathway. ZIC2 expression was significantly correlated with T cells CD4 memory resting, Macrophages M1, and plasma cells, indicating that dysregulated ZIC2 expression in LUAD may directly influence immune infiltration. ZIC2 might be regulated by several different lncRNA-mediated ceRNA mechanisms. In vitro experiments validated the promotive effect of ZIC2 on cell viability and invasion ability of LUAD cells. In vivo experiments validated ZIC2 can accelerate tumor growth in nude mouse. CONCLUSION: ZIC2 regulated by different lncRNA-mediated ceRNA mechanisms may play a critical regulatory role in LUAD through mediating the composition of immune cells in tumor microenvironment.

Analysis of genetically determined gene expression suggests role of inflammatory processes in exfoliation syndrome.

BACKGROUND: Exfoliation syndrome (XFS) is an age-related systemic disorder characterized by excessive production and progressive accumulation of abnormal extracellular material, with pathognomonic ocular manifestations. It is the most common cause of secondary glaucoma, resulting in widespread global blindness. The largest global meta-analysis of XFS in 123,457 multi-ethnic individuals from 24 countries identified seven loci with the strongest association signal in chr15q22-25 region near LOXL1. Expression analysis have so far correlated coding and a few non-coding variants in the region with LOXL1 expression levels, but functional effects of these variants is unclear. We hypothesize that analysis of the contribution of the genetically determined component of gene expression to XFS risk can provide a powerful method to elucidate potential roles of additional genes and clarify biology that underlie XFS. RESULTS: Transcriptomic Wide Association Studies (TWAS) using PrediXcan models trained in 48 GTEx tissues leveraging on results from the multi-ethnic and European ancestry GWAS were performed. To eliminate the possibility of false-positive results due to Linkage Disequilibrium (LD) contamination, we i) performed PrediXcan analysis in reduced models removing variants in LD with LOXL1 missense variants associated with XFS, and variants in LOXL1 models in both multiethnic and European ancestry individuals, ii) conducted conditional analysis of the significant signals in European ancestry individuals, and iii) filtered signals based on correlated gene expression, LD and shared eQTLs, iv) conducted expression validation analysis in human iris tissues. We observed twenty-eight genes in chr15q22-25 region that showed statistically significant associations, which were whittled down to ten genes after statistical validations. In experimental analysis, mRNA transcript levels for ARID3B, CD276, LOXL1, NEO1, SCAMP2, and UBL7 were significantly decreased in iris tissues from XFS patients compared to control samples. TWAS genes for XFS were significantly enriched for genes associated with inflammatory conditions. We also observed a higher incidence of XFS comorbidity with inflammatory and connective tissue diseases. CONCLUSION: Our results implicate a role for connective tissues and inflammation pathways in the etiology of XFS. Targeting the inflammatory pathway may be a potential therapeutic option to reduce progression in XFS.

A Multi-tissue Transcriptome Analysis of Human Metabolites Guides Interpretability of Associations Based on Multi-SNP Models for Gene Expression.

There is particular interest in transcriptome-wide association studies (TWAS) gene-level tests based on multi-SNP predictive models of gene expression-for identifying causal genes at loci associated with complex traits. However, interpretation of TWAS associations may be complicated by divergent effects of model SNPs on phenotype and gene expression. We developed an iterative modeling scheme for obtaining multi-SNP models of gene expression and applied this framework to generate expression models for 43 human tissues from the Genotype-Tissue Expression (GTEx) Project. We characterized the performance of single- and multi-SNP models for identifying causal genes in GWAS data for 46 circulating metabolites. We show that: (A) multi-SNP models captured more variation in expression than did the top cis-eQTL (median 2-fold improvement); (B) predicted expression based on multi-SNP models was associated (false discovery rate < 0.01) with metabolite levels for 826 unique gene-metabolite pairs, but, after stepwise conditional analyses, 90% were dominated by a single eQTL SNP; (C) among the 35% of associations where a SNP in the expression model was a significant cis-eQTL and metabolomic-QTL (met-QTL), 92% demonstrated colocalization between these signals, but interpretation was often complicated by incomplete overlap of QTLs in multi-SNP models; and (D) using a "truth" set of causal genes at 61 met-QTLs, the sensitivity was high (67%), but the positive predictive value was low, as only 8% of TWAS associations (19% when restricted to colocalized associations at met-QTLs) involved true causal genes. These results guide the interpretation of TWAS and highlight the need for corroborative data to provide confident assignment of causality.

Supervised structural learning of semiparametric regression on high-dimensional correlated covariates with applications to eQTL studies.

Expression quantitative trait loci (eQTL) studies utilize regression models to explain the variance of gene expressions with genetic loci or single nucleotide polymorphisms (SNPs). However, regression models for eQTL are challenged by the presence of high dimensional non-sparse and correlated SNPs with small effects, and nonlinear relationships between responses and SNPs. Principal component analyses are commonly conducted for dimension reduction without considering responses. Because of that, this non-supervised learning method often does not work well when the focus is on discovery of the response-covariate relationship. We propose a new supervised structural dimensional reduction method for semiparametric regression models with high dimensional and correlated covariates; we extract low-dimensional latent features from a vast number of correlated SNPs while accounting for their relationships, possibly nonlinear, with gene expressions. Our model identifies important SNPs associated with gene expressions and estimates the association parameters via a likelihood-based algorithm. A GTEx data application on a cancer related gene is presented with 18 novel eQTLs detected by our method. In addition, extensive simulations show that our method outperforms the other competing methods in bias, efficiency, and computational cost.

Cell-Type Heterogeneity in Adipose Tissue Is Associated with Complex Traits and Reveals Disease-Relevant Cell-Specific eQTLs.

Adipose tissue is an important endocrine organ with a role in many cardiometabolic diseases. It is comprised of a heterogeneous collection of cell types that can differentially impact disease phenotypes. Cellular heterogeneity can also confound -omic analyses but is rarely taken into account in analysis of solid-tissue transcriptomes. Here, we investigate cell-type heterogeneity in two population-level subcutaneous adipose-tissue RNA-seq datasets (TwinsUK, n = 766 and the Genotype-Tissue Expression project [GTEx], n = 326) by estimating the relative proportions of four distinct cell types (adipocytes, macrophages, CD4+ T cells, and micro-vascular endothelial cells). We find significant cellular heterogeneity within and between the TwinsUK and GTEx adipose datasets. We find that adipose cell-type composition is heritable and confirm the positive association between adipose-resident macrophage proportion and obesity (high BMI), but we find a stronger BMI-independent association with dual-energy X-ray absorptiometry (DXA) derived body-fat distribution traits. We benchmark the impact of adipose-tissue cell composition on a range of standard analyses, including phenotype-gene expression association, co-expression networks, and cis-eQTL discovery. Our results indicate that it is critical to account for cell-type composition when combining adipose transcriptome datasets in co-expression analysis and in differential expression analysis with obesity-related traits. We applied gene expression by cell-type proportion interaction models (G × Cell) to identify 26 cell-type-specific expression quantitative trait loci (eQTLs) in 20 genes, including four autoimmune disease genome-wide association study (GWAS) loci. These results identify cell-specific eQTLs and demonstrate the potential of in silico deconvolution of bulk tissue to identify cell-type-restricted regulatory variants.

Long non-coding RNA transcriptome of uncharacterized samples can be accurately imputed using protein-coding genes.

Long non-coding RNAs (lncRNAs) play an important role in gene regulation and are increasingly being recognized as crucial mediators of disease pathogenesis. However, the vast majority of published transcriptome datasets lack high-quality lncRNA profiles compared to protein-coding genes (PCGs). Here we propose a framework to harnesses the correlative expression patterns between lncRNA and PCGs to impute unknown lncRNA profiles. The lncRNA expression imputation (LEXI) framework enables characterization of lncRNA transcriptome of samples lacking any lncRNA data using only their PCG profiles. We compare various machine learning and missing value imputation algorithms to implement LEXI and demonstrate the feasibility of this approach to impute lncRNA transcriptome of normal and cancer tissues. Additionally, we determine the factors that influence imputation accuracy and provide guidelines for implementing this approach.

MCAM abnormal expression and clinical outcome associations are highly cancer dependent as revealed through pan-cancer analysis.

MCAM (CD146) is a cell surface adhesion molecule that has been reported to promote cancer development, progression and metastasis and is considered as a potential tumor biomarker and therapeutic target. However, inconsistent reports exist, and its clinical value is yet to be confirmed. Here we took advantage of several large genomic data collections (Genotype-Tissue Expression, The Cancer Genome Atlas and Cancer Cell Line Encyclopedia) and comprehensively analyzed MCAM expression in thousands of normal and cancer samples and cell lines along with their clinical phenotypes and drug response information. Our results show that MCAM is very highly expressed in large vessel tissues while majority of tissues have low or minimal expression. Its expression is dramatically increased in a few tumors but significantly decreased in most other tumors relative to their pairing normal tissues. Increased MCAM expression is associated with a higher tumor stage and worse patient survival for some less common tumors but not for major ones. Higher MCAM expression in primary tumors may be complicated by tumor-associated or normal stromal blood vessels yet its significance may differ from the one from cancer cells. MCAM expression is weakly associated with the response to a few small molecular drugs and the association with targeted anti-BRAF agents suggests its involvement in that pathway which warrants further investigation.

High expression of syndecan-4 is related to clinicopathological features and poor prognosis of pancreatic adenocarcinoma.

OBJECTIVE: Pancreatic adenocarcinoma (PAAD) is a leading cause of cancer-related mortality in adults. Syndecan-4 (SDC4) is involved in cancer pathogenesis. Therefore, this study aimed to explore the expression and clinical significance of SDC4 in PAAD. METHODS: Differentially expressed genes (DEGs) between PAAD and normal pancreas were screened from the GTEx and TCGA databases, and the correlationship between the DEGs and prognosis were analyzed. The prognostic value of the screened SDC4, SERPINE1, and SLC2A1 was evaluated using the Kaplan-Meier curve and SDC4 was subsequently selected as the better candidate. Also, SDC4 expression was analyzed in PAAD tissues, the other risk factors affecting postoperative survival were analyzed using Cox regression analysis, and SDC4-mediated pathways enrichment was identified by GSVA and GSEA. SDC4 expression in PAAD tissues and adjacent normal tissues of selected PAAD patients was detected by RT-qPCR and immunohistochemistry. The correlation between SDC4 and clinical features was evaluated by the χ2 test. RESULTS: SDC4 was highly expressed in PAAD tissues. Elevated SDC4 was correlated with reduced overall survival. SDC4 enrichment pathways included spliceosome function, proteasome activity, pentose phosphate pathway, base excision repair, mismatch repair, DNA replication, oxidative phosphorylation, mitotic spindle formation, epithelial-mesenchymal transition, and G2M checkpoints. SDC4 was elevated in PAAD tissues of PAAD patients compared with adjacent normal tissues. High SDC4 expression was related to metastatic differentiation, TNM stage, lymphatic metastasis, and lower 3-year survival rate. SDC4 was an independent risk factor affecting postoperative survival. CONCLUSION: SDC4 was highly expressed in PAAD and was related to clinicopathological features and poor prognosis, which might be an important index for PAAD early diagnosis and prognosis.

Colocalization analysis of pancreas eQTLs with risk loci from alcoholic and novel non-alcoholic chronic pancreatitis GWAS suggests potential disease causing mechanisms.

BACKGROUND: Previous genome-wide association studies (GWAS) identified genome-wide significant risk loci in chronic pancreatitis and investigated underlying disease causing mechanisms by simple overlaps with expression quantitative trait loci (eQTLs), a procedure which may often result in false positive conclusions. METHODS: We conducted a GWAS in 584 non-alcoholic chronic pancreatitis (NACP) patients and 6040 healthy controls. Next, we applied Bayesian colocalization analysis of identified genome-wide significant risk loci from both, our recently published alcoholic chronic pancreatitis (ACP) and the novel NACP dataset, with pancreas eQTLs from the GTEx V8 European cohort to prioritize candidate causal genes and extracted credible sets of shared causal variants. RESULTS: Variants at the CTRC (p = 1.22 × 10-21) and SPINK1 (p = 6.59 × 10-47) risk loci reached genome-wide significance in NACP. CTRC risk variants colocalized with CTRC eQTLs in ACP (PP4 = 0.99, PP4/PP3 = 95.51) and NACP (PP4 = 0.99, PP4/PP3 = 95.46). For both diseases, the 95% credible set of shared causal variants consisted of rs497078 and rs545634. CLDN2-MORC4 risk variants colocalized with CLDN2 eQTLs in ACP (PP4 = 0.98, PP4/PP3 = 42.20) and NACP (PP4 = 0.67, PP4/PP3 = 7.18), probably driven by the shared causal variant rs12688220. CONCLUSIONS: A shared causal CTRC risk variant might unfold its pathogenic effect in ACP and NACP by reducing CTRC expression, while the CLDN2-MORC4 shared causal variant rs12688220 may modify ACP and NACP risk by increasing CLDN2 expression.

AdaReg: data adaptive robust estimation in linear regression with application in GTEx gene expressions.

The Genotype-Tissue Expression (GTEx) project provides a valuable resource of large-scale gene expressions across multiple tissue types. Under various technical noise and unknown or unmeasured factors, how to robustly estimate the major tissue effect becomes challenging. Moreover, different genes exhibit heterogeneous expressions across different tissue types. Therefore, we need a robust method which adapts to the heterogeneities of gene expressions to improve the estimation for the tissue effect. We followed the approach of the robust estimation based on γ-density-power-weight in the works of Fujisawa, H. and Eguchi, S. (2008). Robust parameter estimation with a small bias against heavy contamination. J. Multivariate Anal. 99: 2053-2081 and Windham, M.P. (1995). Robustifying model fitting. J. Roy. Stat. Soc. B: 599-609, where γ is the exponent of density weight which controls the balance between bias and variance. As far as we know, our work is the first to propose a procedure to tune the parameter γ to balance the bias-variance trade-off under the mixture models. We constructed a robust likelihood criterion based on weighted densities in the mixture model of Gaussian population distribution mixed with unknown outlier distribution, and developed a data-adaptive γ-selection procedure embedded into the robust estimation. We provided a heuristic analysis on the selection criterion and found that our practical selection trend under various γ's in average performance has similar capability to capture minimizer γ as the inestimable mean squared error (MSE) trend from our simulation studies under a series of settings. Our data-adaptive robustifying procedure in the linear regression problem (AdaReg) showed a significant advantage in both simulation studies and real data application in estimating tissue effect of heart samples from the GTEx project, compared to the fixed γ procedure and other robust methods. At the end, the paper discussed some limitations on this method and future work.

Estimating colocalization probability from limited summary statistics.

BACKGROUND: Colocalization is a statistical method used in genetics to determine whether the same variant is causal for multiple phenotypes, for example, complex traits and gene expression. It provides stronger mechanistic evidence than shared significance, which can be produced through separate causal variants in linkage disequilibrium. Current colocalization methods require full summary statistics for both traits, limiting their use with the majority of reported GWAS associations (e.g. GWAS Catalog). We propose a new approximation to the popular coloc method that can be applied when limited summary statistics are available. Our method (POint EstiMation of Colocalization, POEMColoc) imputes missing summary statistics for one or both traits using LD structure in a reference panel, and performs colocalization using the imputed summary statistics. RESULTS: We evaluate the performance of POEMColoc using real (UK Biobank phenotypes and GTEx eQTL) and simulated datasets. We show good correlation between posterior probabilities of colocalization computed from imputed and observed datasets and similar accuracy in simulation. We evaluate scenarios that might reduce performance and show that multiple independent causal variants in a region and imputation from a limited subset of typed variants have a larger effect while mismatched ancestry in the reference panel has a modest effect. Further, we find that POEMColoc is a better approximation of coloc when the imputed association statistics are from a well powered study (e.g., relatively larger sample size or effect size). Applying POEMColoc to estimate colocalization of GWAS Catalog entries and GTEx eQTL, we find evidence for colocalization of 150,000 trait-gene-tissue triplets. CONCLUSIONS: We find that colocalization analysis performed with full summary statistics can be closely approximated when only the summary statistics of the top SNP are available for one or both traits. When applied to the full GWAS Catalog and GTEx eQTL, we find that colocalized trait-gene pairs are enriched in tissues relevant to disease etiology and for matches to approved drug mechanisms. POEMColoc R package is available at https://github.com/AbbVie-ComputationalGenomics/POEMColoc .