Transcriptome-Wide Association Analysis Identifies Novel Candidate Susceptibility Genes for Prostate-Specific Antigen Levels in Men Without Prostate Cancer
===========================================================================================================================================================

* Dorothy M. Chen
* Ruocheng Dong
* Linda Kachuri
* Thomas Hoffmann
* Yu Jiang
* Sonja I. Berndt
* John P. Shelley
* Kerry R. Schaffer
* Mitchell J. Machiela
* Neal D. Freedman
* Wen-Yi Huang
* Shengchao A. Li
* Hans Lilja
* Stephen K. Van Den Eeden
* Stephen Chanock
* Christopher A. Haiman
* David V. Conti
* Robert J. Klein
* Jonathan D. Mosley
* John S. Witte
* Rebecca E. Graff

## Abstract

Deciphering the genetic basis of prostate-specific antigen (PSA) levels may improve their utility to screen for prostate cancer (PCa). We thus conducted a transcriptome-wide association study (TWAS) of PSA levels using genome-wide summary statistics from 95,768 PCa-free men, the MetaXcan framework, and gene prediction models trained in Genotype-Tissue Expression (GTEx) project data. Tissue-specific analyses identified 41 statistically significant (p < 0.05/12,192 = 4.10e-6) associations in whole blood and 39 statistically significant (p < 0.05/13,844 = 3.61e-6) associations in prostate tissue, with 18 genes associated in both tissues. Cross-tissue analyses that combined associations across 45 tissues identified 155 genes that were statistically significantly (p < 0.05/22,249 = 2.25e-6) associated with PSA levels. Based on conditional analyses that assessed whether TWAS associations were attributable to a lead GWAS variant, we found 20 novel genes (11 single-tissue, 9 cross-tissue) that were associated with PSA levels in the TWAS. Of these novel genes, five showed evidence of colocalization (colocalization probability > 0.5): *EXOSC9, CCNA2, HIST1H2BN, RP11-182L21.6*, and *RP11-327J17.2*. Six of the 20 novel genes are not known to impact PCa risk. These findings yield new hypotheses for genetic factors underlying PSA levels that should be further explored toward improving our understanding of PSA biology.

## Introduction

Prostate-specific antigen (PSA) is a serine protease of the human tissue kallikrein-related (KLK) peptidase family. Serum levels are commonly used as a biomarker for detection, monitoring, and risk stratification of prostate cancer (PCa).1–3 The process by which a small portion of the highly abundant PSA in the prostate is released into the blood and elevated in men prostate cancer is incompletely understood, but disruption of the prostate gland architecture by neoplastic transformation has been speculated as a possible mechanism.1, 3, 4 PSA levels can additionally be influenced by age, race, body mass index (BMI), infection, prostate volume, benign prostate hyperplasia (BPH), and germline genetics.5, 6

PSA screening has been used for over 25 years for the detection of PCa, which is the second leading cause of cancer death among men in the United States.1, 7–10 However, low test specificity and discrimination have complicated the use and interpretation of PSA as a screening tool.9, 11 Long-term outcomes data from large population-based randomized PSA screening trials show that it significantly reduces death from prostate cancer but also results in considerable overdiagnosis of low risk disease.12, 13 Screening might be improved were it to account for variation in PSA levels that is attributable to germline genetics rather than PCa. Twin studies and genome-wide association studies (GWAS) have estimated that 30-45% of serum PSA variation is influenced by hereditary factors.14, 15 A recent GWAS meta-analysis from our group identified 128 independent variants associated with serum PSA levels, which explained approximately 9% of genetic variation in PSA.16, 17

Additional genetic variation in PSA levels may be determined by analyses of genes, as opposed to individual variants. Such work could lead to the identification of regions or biological pathways that affect PSA levels and provide clarity on mechanisms underlying constitutional increases in PSA levels in the absence of carcinogenesis. Transcriptome-wide association studies (TWAS) using expression quantitative trait loci (eQTL) allow for identifying genes whose cis-regulated expression is associated with complex polygenic traits.18–21 We thus performed a TWAS of PSA levels based on summary statistics from a GWAS meta-analysis of 95,768 men without PCa to identify genes associated with PSA levels and prioritize them for functional investigation.15, 22

## Methods

### Discovery Populations

A total of 95,768 PCa-free men from the following five study populations were included: UK Biobank, Kaiser Permanente’s Genetic Epidemiology Research on Adult Health and Aging, Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial, Vanderbilt University Medical Center’s BioVU, and Malmö Diet and Cancer Study. The men included in these analyses were restricted to PCa-free individuals (i.e., without a PCa diagnosis or history of prostate resection, where information was available). Additionally, analyses were restricted to men with PSA values between 0.01 ng/mL and 10 ng/mL. Additional details about the study cohorts are described in detail in Kachuri et al. and Table S1. Median PSA levels were used for individuals with multiple PSA measurements available, with the exception of PLCO, which used PSA levels at the start of follow-up of the trial.

### Quality Control and GWAS Meta-analysis

Genotyping, imputation, study-specific quality control of genetic data, and the GWAS meta-analysis have been previously described for Kachuri et al.17 Briefly, ancestry and study-specific GWAS analyses used linear regression of log(PSA) on genetic variants, age, and the first 10 genetic ancestry principal components. After results across studies were meta-analyzed within ancestral groups, ancestry-specific summary statistics from individuals of European (n=85,824), African (n=3,509), East Asian (n=3,337), and Hispanic/Latino (n= 3,098) ancestry were meta-analyzed to generate multi-ancestry summary statistics.

### MetaXcan Transcriptome-wide Gene-based Analysis

We undertook a TWAS using the MetaXcan approach, which directly estimates Z-scores for associations between gene expression levels and PSA levels using meta-analyzed GWAS summary statistics with tissue-specific prediction models. Multivariate adaptive shrinkage (mashr) prediction models were trained on GTEx (version 8) eQTL data for whole blood and, separately, prostate tissue. The covariances and weights for models trained on individual tissue types were from the PredictDB Data Repository ([https://predictdb.org/](https://predictdb.org/)).23 We also undertook a cross-tissue analysis across 45 GTEx tissues that excluded tissues found primarily or exclusively in women to identify potential associations with genes that may have been missed in single-tissue analyses. Tissue types found primarily or exclusively in females (mammary breast, ovary, uterus, vagina) were excluded from the multiple tissue analysis. A Bonferroni correction was implemented based on the number of genes tested (whole blood: 0.05/12,192 =4.10e-6; prostate: 0.05/13,844 =3.61e-6; cross-tissue: 0.05/22,249 =2.25e-6).

### Conditional Analysis Ascertaining Independence of TWAS and Previous GWAS Findings

To identify the set of TWAS genes for PSA that were novel in comparison to previous GWAS results, we first limited consideration to significant TWAS genes that do not contain genome-wide significant (p < 5e-8) variants from the Kachuri et al GWAS (i.e., within exact gene boundaries).17 Second, for the remaining genes, we performed conditional analyses using COJO24 that simultaneously modeled eQTLs for a given TWAS-identified gene and GWAS results. For individual tissue analyses, eQTLs used in the prediction models of TWAS-significant genes were selected from the tissue-specific PredictDB TWAS eQTL weight files for each gene. For the cross-tissue analyses, the maximum noncollinear eQTLs (variant inflation factor < 8 within a window size of 50kb) were extracted from all 45 GTEx tissues (i.e., to remove variants in linkage disequilibrium). COJO was performed using LD reference panels from 10,000 European ancestry UK Biobank participants (as 89.6% of our population was European ancestry).24 GWAS summary statistics were then conditioned on the eQTLs for one gene at a time. If the set of genome-wide significant variants was reduced after conditioning on eQTLs used to predict gene expression, then the eQTLs and GWAS hits were assumed to be in LD and thus not independent of one another. (Methods that assess the conditional fit of eQTLs on GWAS hits are lacking.) However, if the set of genome-wide significant variants remained unchanged after conditional analyses, then we considered the TWAS gene to be independent of GWAS variants from Kachuri et al.17

### Colocalization Analyses

To further investigate whether genes and variants exhibited shared genetic signals, we used fastENLOC software (version 2) to undertake colocalization analyses of eQTLs for all genes with prediction models across TWAS analyses and all variants investigated in the Kachuri et al. GWAS.17, 25, 26 For these analyses, LD blocks were based on 1000 Genome Project European ancestry individuals.27 Z-scores were calculated by dividing the betas from GWAS summary statistics by their standard errors and converting them to posterior inclusion probabilities.25 Finally, the GWAS summary statistics were colocalized for 49 GTEx version 8 tissues using pre-computed GTEx multiple-tissue annotations. Signal-level results returned regional-level colocalization probabilities (RCP) between eQTL and GWAS signals, which sum up the variant-level colocalization probabilities of correlated variants within an LD block that harbors a single GWAS association signal. Gene-results were based on gene-level regional colocalization probabilities (GRCP), which represent the probability that a candidate gene contains at least one colocalized variant.26 RCP and GRCP values > 0.5 indicate strong evidence of shared genetic signals between eQTLs and GWAS variants. 26

### Pathway Enrichment Analysis

To explore the potential biological relevance of PSA-associated genes, we applied the Enrichr tool to all significant genes identified in the whole blood, prostate, and cross-tissue analyses to assess enrichment against three gene set libraries: Kyoto Encyclopedia of Genes and Genomes (KEGG) 2021 human; gene ontology (GO) biological process 2021; and GO molecular function 2021. Enrichment was assessed by multiplying the p-value from a Fisher’s exact test with the z-score of the deviation from the expected rank.28, 29 Pathways with a Benjamini-Hochberg-corrected p-value < 0.05 were considered significantly enriched.

## Results

Using MetaXcan and GWAS summary statistics based on 95,769 individuals (Table S1), TWAS analyses were conducted for whole blood, prostate, and cross-tissue matrices. The sample was primarily European (n=85,824), although individuals of African (n=3,509), East Asian (n=3,337), and Hispanic/Latino (n= 3,098) ancestry were also included in analyses. The median PSA value across all possible values across different cohorts was 2.35 ng/mL.

### Whole Blood TWAS

In the TWAS based on whole blood gene expression models, 41 out of 12,192 genes were associated with PSA levels at the Bonferroni corrected threshold (Figure 1, Table S2). Expression of 20 genes was positively associated with elevated PSA, and expression of 21 genes was inversely associated. Two of the significant genes are located at 19q13.33, which contains *KLK3*, the gene encoding serum PSA.30 While a prediction model for *KLK3* was not available in whole blood, a member of the same gene family at 19q13.33, *KLK2*, was the gene most strongly associated with PSA levels (p = 1.28e-62). Twenty-four of the 41 genes did not contain genome-wide significant variants and had not been annotated in previous GWAS. Among them, 7 genes were independent of previously published GWAS findings (Table 1).17 Of these, increased expression of 6 genes was associated with elevated PSA levels in PCa-free men: *GPBP1L1* (1p34.1, p = 2.21e-7); *TMEM69* (1p34.1, p = 2.21e-7); *UQCRH* (1p33, p = 4.76e-7); *ACTRT3* (3q26.2, p = 5.90e-7); *EXOSC9* (4q27, p = 1.80e-6); and *CCNA2* (4q27, p = 7.80e-7). Decreased expression of *ITH4* was associated with increased PSA levels (3p21.1, p = 3.65e-6).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/05/05/2023.05.04.23289526/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2023/05/05/2023.05.04.23289526/F1)

Figure 1. 
**TWAS results by chromosome**, as depicted by Miami plot for whole blood analyses. Each dot represents a gene, and the red dashed lines represent the significance thresholds after Bonferroni correction (whole blood: 4.10e-6, prostate: 3.61e-6, cross-tissue: 2.25e-6). For the Miami plots in figures 1 and 2, all statistically significant genes are annotated, and the genes independent of results from the largest prior GWAS are highlighted in yellow; the genes in the upper half of the plots were positively associated with PSA levels, and the genes in the lower half of the plots were inversely associated with PSA levels. For the Manhattan plot in figure 3, the 30 genes with the smallest association p-values are annotated and colored blue, and significant genes from the whole blood and prostate tissue analyses are annotated and colored yellow. Directions of association between genes and PSA levels are not available for cross-tissue analyses due to the joint tissue inference methodology used in S-MultiXcan.

View this table:
[Table 1.](http://medrxiv.org/content/early/2023/05/05/2023.05.04.23289526/T1)

Table 1. 
Novel genes associated with prostate-specific antigen (PSA) levels in transcriptome-wide association studies (TWAS) based on whole blood, prostate, and cross-tissue predicted gene expression*

### Prostate Tissue TWAS

In TWAS based on prostate tissue gene expression models, 39 out of 13,884 genes were statistically significantly associated with PSA levels (Figure 2, Table S3). Among them, increased expression of 18 and decreased expression of 21 was associated with elevated PSA levels. No prediction models were available for genes in the KLK family in prostate tissue. Instead, *MMP7* on 11q22 was most strongly associated with PSA levels (p= 2.78e-18). *MMP7* was among 18 PSA-associated TWAS genes containing known genome-wide significant variants. Of the remaining 21 significant genes, six were conditionally independent of previous GWAS results after conditioning on the eQTL variants used in gene prediction models (Table 1). Of these, increased expression of 3 genes was associated with elevated PSA levels: *LRRC41* (1p33, p = 4.76e-7), *HIST1H2BN* (6p22.1, p = 8.86e-8), and *HEXIM1* (17q21.31, p = 4.95e-7). Decreased expression of the remaining 3 genes was associated with increased PSA levels: *ITIH4* (3p21.1, p = 2.69e-6), *EXOSC9* (4q27 p = 1.18e-6), and *OTOA* (16p12.2, p = 1.44e-6).

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/05/05/2023.05.04.23289526/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2023/05/05/2023.05.04.23289526/F2)

Figure 2: 
**TWAS results by chromosome**, as depicted by Miami plot for prostate tissue analyses. Each dot represents a gene, and the red dashed lines represent the significance thresholds after Bonferroni correction (whole blood: 4.10e-6, prostate: 3.61e-6, cross-tissue: 2.25e-6). For the Miami plots in figures 1 and 2, all statistically significant genes are annotated, and the genes independent of results from the largest prior GWAS are highlighted in yellow; the genes in the upper half of the plots were positively associated with PSA levels, and the genes in the lower half of the plots were inversely associated with PSA levels. For the Manhattan plot in figure 3, the 30 genes with the smallest association p-values are annotated and colored blue, and significant genes from the whole blood and prostate tissue analyses are annotated and colored yellow. Directions of association between genes and PSA levels are not available for cross-tissue analyses due to the joint tissue inference methodology used in S-MultiXcan.

### Cross-Tissue TWAS

Using S-MultiXcan to integrate signals from 45 tissues excluding tissues primarily or exclusively present in women, we identified 155 out of 22,249 genes whose predicted expression was associated with elevated PSA levels in cross-tissue matrices (Figure 3). Among the 155 statistically significant genes, 8 are located on 19q13.33 and part of the expanded human *KLK* gene family. We found 71 genes that do not contain genome-wide significant loci within their gene boundaries. Conditional analyses determined that 13 genes from this set were independent of known genome-wide significant variants: *CDC42* (1p36.12, p = 6.24e-07), *GPBP1L1* (1p34.1, p = 2.08e-08), *SLBP* (4p16.3, p = 8.40e-08), *EXOSC9* (4q27, p = 8.60e-07), *CCNA2* (4q27, p = 1.02e-06), *ANXA8L1* (10q11.22, p = 2.17e-07), *RP11-182L21.6* (10q22.3, p = 1.54e-07), *E2F7* (12q21.2, p = 8.54e-09), *C15orf62* (15q15.1, p = 2.01e-07), *RP11-327J17.2* (15q26.2, p = 5.82e-07), *C16orf47* (16q22.3, p = 1.11e-07), *MYOCD* (17p12, p=1.20E-07), and *ABCD4* (17q21.31, p = 2.17e-06) (Table 1).

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/05/05/2023.05.04.23289526/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2023/05/05/2023.05.04.23289526/F3)

Figure 3: 
**TWAS results by chromosome**, as depicted by Manhattan plot for cross-tissue analyses. Each dot represents a gene, and the red dashed lines represent the significance thresholds after Bonferroni correction (whole blood: 4.10e-6, prostate: 3.61e-6, cross-tissue: 2.25e-6). For the Miami plots in figures 1 and 2, all statistically significant genes are annotated, and the genes independent of results from the largest prior GWAS are highlighted in yellow; the genes in the upper half of the plots were positively associated with PSA levels, and the genes in the lower half of the plots were inversely associated with PSA levels. For the Manhattan plot in figure 3, the 30 genes with the smallest association p-values are annotated and colored blue, and significant genes from the whole blood and prostate tissue analyses are annotated and colored yellow. Directions of association between genes and PSA levels are not available for cross-tissue analyses due to the joint tissue inference methodology used in S-MultiXcan.

### Significant Genes Overlapping Across TWAS

Out of the 41 genes that were significantly associated with PSA levels in the whole blood TWAS, 34 were also imputed in the prostate tissue TWAS. Thirty (30) out of 39 significantly associated genes from the prostate tissue TWAS were also evaluated in whole blood analyses. Of the 155 significantly associated genes detected in the cross-tissue analysis, 89 genes and 105 genes had transcriptome prediction models available for whole blood and prostate tissues, respectively.

Across all TWAS in whole blood, prostate tissue, and cross-tissue models, we identified 13 genes whose predicted expression was associated with elevated PSA levels: *RAB29* (1q32.1)*, LDAH* (2p24.1)*, ZFP36L2* (2p21)*, SIDT1* (3q13.2)*, TRIM59* (3q25.33)*, EXOSC9* (4q27)*, HIBADH* (7p15.2), *TAX1BP1* (7p15.2)*, MEGF9* (9q33.2)*, ADD3* (10q25.1)*, APIP* (11p13)*, LINC01569* (16p13.3), and *EIF2S2* (20q11.22) (Figure 4). Of these genes, only *EXOSC9* was conditionally independent after accounting for prior loci and demonstrated evidence of colocalization (described below). Five (5) associations were not statistically significant in the cross-tissue analysis, but had significant tissue-specific signals in whole blood and/or prostate. Seventeen significant genes were shared by only the whole blood and cross-tissue TWAS, six of which did not have prediction models in prostate tissue. Similarly, the prostate tissue and cross-tissue TWAS shared 14 significant genes; of these, seven genes did not have prediction models in whole blood.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/05/05/2023.05.04.23289526/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2023/05/05/2023.05.04.23289526/F4)

Figure 4. Statistically significant genes overlapping across TWAS analyses.
Of the significant gene sets from whole blood (41 genes), prostate tissue (39 genes), and cross-tissue analyses (155 genes), 13 genes were significant across all three analyses, 36 genes were significant in exactly two analyses, and 124 genes were significant in a single analysis.

### Colocalization Analysis

Colocalization identified 372 signals with RCP > 0.5 at the variant level (Table S6) and 307 genes with GRCP > 0.5 at the gene level (Table S7). The top 10 genes with the largest GRCP were *OTX1* (2p15)*, MAFF* (22q13.1)*, FUT2* (19q13.33)*, EMSY* (11q13.5)*, IFT80* (3q25.33)*, EXOSC9* (4q27)*, SMC4* (3q25.33)*, RCAN3* (1p36.11)*, PBXIP1* (1q21.3), and *PMVK* (1q21.3). Of these genes, only *EXOSC9* was determined to be independent of previously discovered GWAS loci.17

Among 155 significant genes from the cross-tissue TWAS, locus-level results identified 50 genes with strong evidence of colocalization (Table 2). Variant-level results identified 66 regions in 51 genes (Table S5). When comparing the overlap between transcriptome-wide significant genes and results from colocalization analyses, we identified 5 novel genes (*EXOSC9, CCNA2, HIST1H2BN, RP11-182L21.6*, and *RP11-327J17.2*) with GRCP/RCP > 0.5. Two of these genes are long non-coding RNA while the remaining three are protein-coding genes.

View this table:
[Table 2.](http://medrxiv.org/content/early/2023/05/05/2023.05.04.23289526/T2)

Table 2. 
Gene level colocalization results for significant genes from the cross-tissue transcriptome-wide association study*

### Pathway Enrichment Analysis

TWAS genes were significantly enriched in 23 pathways underlying molecular activities and biological processes present in KEGG and GO gene set databases (Table 3, Figure 5). Enrichment analyses in the whole blood gene set identified the renin-angiotensin-aldosterone system in KEGG, which works to regulate arterial blood pressure. No pathways were significantly enriched by prostate and cross-tissue gene sets in the KEGG catalog.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/05/05/2023.05.04.23289526/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2023/05/05/2023.05.04.23289526/F5)

Figure 5. GO term enrichment results from the cross-tissue TWAS.
The GO molecular function and biological process data repositories were queried for pathway enrichment analyses. The color gradient represents the magnitude of p-values, with darker colors corresponding to smaller values. The size of the circles represents the number of genes in the pathway.

View this table:
[Table 3.](http://medrxiv.org/content/early/2023/05/05/2023.05.04.23289526/T3)

Table 3. 
Results from pathway analyses based on all significant genes identified in the whole blood, prostate tissue, and cross-tissue transcriptome-wide association studies

Enrichment analyses within the GO biological process and molecular function gene set libraries identified 23 pathways that were significantly enriched in the cross-tissue gene set. No pathways were significantly enriched in GO for whole blood nor prostate tissue TWAS sets. Of note, the biological process underlying the regulation of cyclin-dependent protein serine/threonine kinase (CDK) activity in cell cycle progression is particularly relevant for the function of *CCNA2*. We also observed enrichment in 13 pathways involving both single and double stranded DNA binding mechanisms and regulation by transcription factors, including RNA polymerase II (pol II), a core component of the DNA transcription machinery. Moreover, two pathways underlying extracellular matrix and cellular component disassembly were enriched in the cross-tissue gene set. Analogous to the findings in KEGG for the whole blood TWAS, biological pathways underlying the regulation of systemic arterial blood pressure were also identified in the GO biological process catalog.

## Discussion

This TWAS of 95,768 PCa-free men identified 173 genes whose predicted expression levels were significantly associated with PSA. Conditional analyses identified 20 novel candidate susceptibility genes for PSA, and colocalization analyses highlighted five of them: *EXOSC9* and *CCNA2* at 4q27, *HIST1H2BN* at 6p22.1, *RP11-182L21.6* at 10q22.3, and *RP11-327J17.2* at 15q26.2. Pathway enrichment analyses across three gene ontology catalogs implicated regulatory pathways related to transcription, cell signaling, and disassembly of cellular and noncellular components such as the extracellular matrix.

The five colocalized genes introduce new hypotheses and insights regarding the genetic mechanisms regulating PSA production. Of particular interest are the signals at the 4q27 locus (*EXOSC9* and *CCNA2*), a genomic region that was not detected in previous GWAS of PSA.16, 17 Researchers have previously investigated an autoimmune related block on the 4q27 locus and did not find an association with PCa risk. 32–35 *EXOSC9*, which encodes a core protein involved in the RNA degradation machinery in humans, was identified across all three TWAS analyses. Currently, there is no literature supporting an association between *EXCOSC9* and PSA levels or PCa risk. Genetically predicted expression of *EXOSC9* in whole blood was positively associated with PSA levels, while the opposite direction of effect was observed for expression in prostate tissue. There exists limited tissue-specific knowledge of the physiological and pathological function of *EXOSC9* in non-cancer cell lines.34, 35 Literature documenting the relationship between *CCNA2* and PSA is also limited, though *CCNA2* was previously tested as a candidate autoantibody signature marker for distinguishing PCa from BPH in patients with elevated serum PSA;36 it was not determined to be a top antibody signature of any specific PCa targets. Gene co-expression network analyses have linked *CCNA2* to biochemical recurrence and survival in men with PCa.37, 38 Further exploration of *CCNA2* expression in individuals without PCa may unveil biological pathways that influence PSA levels.

Pathway enrichment results also highlighted a number of interesting pathways involving Pol II, which plays a critical role in the regulated synthesis of both protein-coding and noncoding messenger RNA in eukaryotic genomes.39–41 In addition to its numerous roles in DNA transcription, Pol II is known to associate with androgen receptor binding regions of the PSA enhancer and promoter to initiate androgen-dependent transcription.42 The two colocalized lncRNA genes, *RP11-182L21.6* and *RP11-327J17.2*, could feasibly regulate epigenetic modification through histone or DNA acetylation or methylation.43 That histone deacetylase binding was enriched in our pathway analyses suggests a possible mechanism through which lncRNAs could affect constitutive PSA levels. Aberrant expression of various lncRNA products has also been observed for various cancers, including PCa.44

Because PSA is used for both PCa detection and monitoring of PCa progression, it is difficult to disentangle the mechanisms underlying our observed associations. PSA-related screening bias may account for a portion of the observed relationship between gene expression and PSA levels. Genes may also exert pleiotropic effects on PCa and PSA through overlapping biological mechanisms.1, 9 Nevertheless, our restriction to PCa-free men for all analyses minimizes the potential for reverse causation and bolsters confidence that the observed genetic signals inform constitutive PSA levels. Fourteen of the 20 novel PSA-related genes identified from conditional analyses, spanning 12 genomic regions not implicated by prior PSA GWAS, have been associated with prostatic malignancies: *CDC42* (1p36.12), *GPBP1L1* and *TMEM69* (1p34.1), *LRRC41* (1p33), *ITIH4* (3p21.1), *SLBP* (4p16.3), *CCNA2* (4q27), *HIST1H2BN* (6p22.1), ANXA8L1 (10q11.2), *E2F7* (12q21.2), *C15orf62* (15q15.1), *OTOA* (16p12.2), *C16orf47* (16q22.3), and *HEXIM1*(17q21.31).45–55 Colocalized gene *HIST1H2BN*, which encodes a component of a core nucleosome histone, has been linked to PCa cell growth and epithelial-mesenchymal transition through upregulated NF-kB/Rel expression.45 Activation of the NF-kB pathway can induce activation of the PSA promoter-enhancer, even in the absence of androgens, and NF-kB can directly bind to the PSA enhancer in prostate cancer cell lines.56 To our knowledge, no experiments have been conducted to investigate the relationship between NF-kB and PSA in PCa-free populations.

Conditional analyses further identified 6 genes that have not been implicated in PCa susceptibility in gene function, experimental, or human population research: *UQCRH* (1p33)*, ACTRT3* (3q26.2)*, EXOSC9* (4q27)*, RP11-182L21.6* (10q22.3), *RP11-327J17.2* (15q26.2), and *ACBD4* (17q21.31)*. ACTRT3* is of particular interest, as it is critical in regulating sperm nucleus cytomorphology upstream of the processing of spermatid into mature motile sperm.57 Increased expression of *ACTRT3*, which forms a testis-specific profilin III-*ACTRT3* complex that facilitates male germ cell head cytomorphology and maintains sperm motility in animal models, was associated with elevated PSA levels in the TWAS based on whole blood.57, 58 Though there is no documented link between *ACTRT3* and PSA levels or PCa risk, the pronounced role of *ACTRT3* in mediating conformational changes in sperm nuclei suggests possible shared biological pathways between PSA and the production and processing of male germ cells.5, 57 Loci that map to the multigenic region that contains *ACTRT3* have been linked to melanoma, colorectal cancer, and lung cancer susceptibility.59–61 The functional role of *ACTRT3* is poorly characterized and warrants further analysis.

Our study had several limitations. First, nearly 90% of men in the study population were of primarily European ancestry, and 85% of the tissue samples used to train the TWAS models were derived from European populations.62 Therefore, our TWAS may have missed ancestry-specific signals, and our findings may not be generalizable to broader ancestral populations. To better characterize the genetic mechanisms underlying circulating PSA levels, it will be imperative for future analyses to expand to multi-ethnic populations; diversity in genetic studies of PSA levels is critical for equitably improving PSA screening.63 Second, our study assessed only cis-eQTLs, so any trans-eQTL effects are not incorporated. Third, colocalization analyses have a high type 2 error rate and may be underpowered to detect shared association signals.26 Fourth, although we restricted our analyses to men who had not been diagnosed with PCa, we cannot rule out the possibility of latent, undiagnosed disease or disease diagnosed at a later timepoint. However, the prevalence of undiagnosed PCa in our population was likely to be low on account of low prevalence at the time of the first PSA measurement or increased monitoring and surveillance in longitudinal cohorts.7, 9 Finally, many of the novel genes that we identified were clustered at multi-gene loci, in part due to co-regulation by a shared set of eQTLs. Sentinel genes at these loci should be interpreted with caution, as there may be correlated predicted expression.64

Our study also had several key strengths. The use of GWAS summary statistics from a cohort of 95,768 men provided us with high statistical power to quantify PSA-associated genes. In addition to conducting tissue-specific association analyses in tissues that are biologically meaningful for PSA, we integrated association signals across 45 GTEx tissues to improve power for genes with similar regulatory mechanisms across tissues. This method allowed us to capture PSA-related genes with expression patterns that were significant for specific tissues, as well as genes with expression that was similar across multiple tissues. Our study also used COJO analysis to find novel genes conditionally independent from known GWAS variants.

In summary, our TWAS identified gene expression profiles associated with PSA levels in men without PCa. These findings provide several novel hypotheses for genes that affect constitutive PSA. Further exploration of these results, including functional analyses of these genes in in-vivo settings, will augment our understanding of the genetic etiology of PSA variation. Transcriptomic analyses might also be vertically integrated with downstream -omic approaches to uncover complete mechanisms through which genetics influence circulating PSA levels. In addition, TWAS findings may be used to develop polygenic transcriptome risk scores65 for PSA levels, which could be leveraged for improving PSA as a tool for PCa screening.

## Supporting information

Supplemental Tables [[supplements/289526_file03.xlsx]](pending:yes)

## Data Availability

GWAS summary statistics used in this analysis are available for download from Kachuri et al are available from the following Zenodo repository: https://doi.org/10.5281/zenodo.7460134. Transcriptome prediction weights and models used for MetaXcan analyses are available from: https://predictdb.org.

[https://doi.org/10.5281/zenodo.7460134](https://doi.org/10.5281/zenodo.7460134) 

[https://predictdb.org](https://predictdb.org) 

## Declaration of Interests

JSW is a non-employee, cofounder of Avail Bio. HL is named on a patent for assays to measure intact prostate-specific antigen and a patent for a statistical method to detect prostate cancer commercialized by OPKO Health (4KScore). HL receives royalties from sales of the assay and has stock in OPKO Health. HL serves on the Scientific Advisory Board for Fujirebio Diagnostics Inc and owns stock in Diaprost AB and Acousort AB.

## Data Availability

GWAS summary statistics used in this analysis are available for download from Kachuri et al are available from the following Zenodo repository: [https://doi.org/10.5281/zenodo.7460134](https://doi.org/10.5281/zenodo.7460134). Transcriptome prediction weights and models used for MetaXcan analyses are available from: [https://predictdb.org](https://predictdb.org).

## Funding sources

The Precision PSA study is supported by funding from the National Institutes of Health (NIH) and National Cancer Institute (NCI) under award number R01CA241410 (PI: JSW). REG is supported by a Prostate Cancer Foundation Young Investigator Award. LK is supported by funding from National Cancer Institute (R00CA246076). JPS is supported by funding from the National Institute of General Medical Sciences (T32GM007347). HL is supported in part by funding from NIH/NCI (P30-CA008748, U01-CA199338, R01CA244948) and the Swedish Cancer Society (Cancerfonden 20 1354 PjF). RK is supported by funding from the NIH (R01 CA244948 and R01 CA175491).

The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

## Footnotes

*   † These authors jointly supervised this work

*   Received May 4, 2023.
*   Revision received May 4, 2023.
*   Accepted May 5, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Lilja H, Ulmert D, Vickers AJ. Prostate-specific antigen and prostate cancer: prediction, detection and monitoring. Nat Rev Cancer. 2008;8(4):268–278. doi:10.1038/nrc2351
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrc2351&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18337732&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000254133500015&link_type=ISI) 

2.  2.Cucchiara V, Cooperberg MR, Dall’Era M, et al. Genomic Markers in Prostate Cancer Decision Making. Eur Urol. 2018;73(4):572–582. doi:10.1016/j.eururo.2017.10.036
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.eururo.2017.10.036&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29129398&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

3.  3.LeBeau AM, Kostova M, Craik CS, Denmeade SR. Prostate-specific antigen: an overlooked candidate for the targeted treatment and selective imaging of prostate cancer. Biol Chem. 2010;391(4):333–343. doi:10.1515/BC.2010.044
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1515/BC.2010.044&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20180648&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

4.  4.Qiu SD, Young CY, Bilhartz DL, et al. In situ hybridization of prostate-specific antigen mRNA in human prostate. J Urol. 1990;144(6):1550–1556. doi:10.1016/s0022-5347(17)39797-5
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0022-5347(17)39797-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1700164&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990EK57900073&link_type=ISI) 

5.  5.Lilja H, Oldbring J, Rannevik G, Laurell CB. Seminal vesicle-secreted proteins and their reactions during gelation and liquefaction of human semen. J Clin Invest. 1987;80(2):281–285.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI113070&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3611349&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1987J417200001&link_type=ISI) 

6.  6.Pinsky PF, Kramer BS, Crawford ED, et al. Prostate volume and prostate-specific antigen levels in men enrolled in a large screening trial. Urology. 2006;68(2):352–356. doi:10.1016/j.urology.2006.02.026
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.urology.2006.02.026&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16904451&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

7.  7.Ilic D, Neuberger MM, Djulbegovic M, Dahm P. Screening for prostate cancer. Cochrane Database Syst Rev. 2013;2013(1):CD004720. doi: 10.1002/14651858.CD004720.pub3.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/14651858.CD004720.pub3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23440794&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

8.  8.Fenton JJ, Weyrich MS, Durbin S, Liu Y, Bang H, Melnikow J. Prostate-Specific Antigen-Based Screening for Prostate Cancer: Evidence Report and Systematic Review for the US Preventive Services Task Force. JAMA. 2018;319(18):1914–1931. doi:10.1001/jama.2018.3712
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2018.3712&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29801018&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

9.  9.Ilic D, Djulbegovic M, Jung JH, et al. Prostate cancer screening with prostate-specific antigen (PSA) test: a systematic review and meta-analysis. The BMJ. 2018;362:k3519. doi:10.1136/bmj.k3519
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNjIvc2VwMDVfNy9rMzUxOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA1LzA1LzIwMjMuMDUuMDQuMjMyODk1MjYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

10. 10.Potosky AL, Feuer EJ, Levin DL. Impact of screening on incidence and mortality of prostate cancer in the United States. Epidemiol Rev. 2001;23(1):181–186. doi:10.1093/oxfordjournals.epirev.a000787
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/oxfordjournals.epirev.a000787&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11588846&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171238600029&link_type=ISI) 

11. 11.Han PKJ, Kobrin S, Breen N, et al. National Evidence on the Use of Shared Decision Making in Prostate-Specific Antigen Screening. Ann Fam Med. 2013;11(4):306–314. doi:10.1370/afm.1539
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiYW5uYWxzZm0iO3M6NToicmVzaWQiO3M6ODoiMTEvNC8zMDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNS8wNS8yMDIzLjA1LjA0LjIzMjg5NTI2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

12. 12.Hugosson J, Roobol MJ, Månsson M, et al. A 16-yr Follow-up of the European Randomized study of Screening for Prostate Cancer. Eur Urol. 2019;76(1):43–51. doi:10.1016/j.eururo.2019.02.009
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.eururo.2019.02.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30824296&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

13. 13.de Koning HJ, Gulati R, Moss SM, et al. The efficacy of prostate-specific antigen screening: Impact of key components in the ERSPC and PLCO trials. Cancer. 2018;124(6):1197–1206. doi:10.1002/cncr.31178
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cncr.31178&link_type=DOI) 

14. 14.Bansal A, Murray DK, Wu JT, Stephenson RA, Middleton RG, Meikle AW. Heritability of Prostate-Specific Antigen and Relationship with Zonal Prostate Volumes in Aging Twins*. J Clin Endocrinol Metab. 2000;85(3):1272–1276. doi:10.1210/jcem.85.3.6399
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.85.3.1272&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10720075&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

15. 15.Pilia G, Chen WM, Scuteri A, et al. Heritability of Cardiovascular and Personality Traits in 6,148 Sardinians. PLoS Genet. 2006;2(8):e132. doi:10.1371/journal.pgen.0020132
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.0020132&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16934002&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

16. 16.Hoffmann TJ, Passarelli MN, Graff RE, et al. Genome-wide association study of prostate-specific antigen levels identifies novel loci independent of prostate cancer. Nat Commun. 2017;8:14248. doi:10.1038/ncomms14248
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms14248&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28139693&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

17. 17.Kachuri L, Hoffmann TJ, Jiang Y, et al. Incorporating Genetic Determinants of Prostate-Specific Antigen Levels Improves Prostate Cancer Screening. MedRvix. doi:10.1101/2022.04.18.22273850
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMi4wNC4xOC4yMjI3Mzg1MHYzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDUvMDUvMjAyMy4wNS4wNC4yMzI4OTUyNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

18. 18.Nicolae DL, Gamazon E, Zhang W, Duan S, Dolan ME, Cox NJ. Trait-associated SNPs are more likely to be eQTLs: annotation to enhance discovery from GWAS. PLoS Genet. 2010;6(4):e1000888. doi:10.1371/journal.pgen.1000888
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1000888&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20369019&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

19. 19.Gamazon ER, Wheeler HE, Shah KP, et al. A gene-based association method for mapping traits using reference transcriptome data. Nat Genet. 2015;47(9):1091–1098. doi:10.1038/ng.3367
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3367&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26258848&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

20. 20.Gusev A, Ko A, Shi H, et al. Integrative approaches for large-scale transcriptome-wide association studies. Nat Genet. 2016;48(3):245–252. doi:10.1038/ng.3506
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3506&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26854917&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

21. 21.Pasaniuc B, Price AL. Dissecting the genetics of complex traits using summary association statistics. Nat Rev Genet. 2017;18(2):117–127. doi:10.1038/nrg.2016.142
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg.2016.142&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

22. 22.Barbeira AN, Dickinson SP, Bonazzola R, et al. Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics. Nat Commun. 2018;9(1):1825. doi:10.1038/s41467-018-03621-1
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-03621-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29739930&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

23. 23.Barbeira, A.N., Bonazzola, R., Gamazon, E.R. et al. Exploiting the GTEx resources to decipher the mechanisms at GWAS loci. Genome Biol. 2021;22(49). doi: 10.1186/s13059-020-02252-4
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-020-02252-4&link_type=DOI) 

24. 24.Yang J, Ferreira T, Morris AP, et al. Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. Nat Genet. 2012;44(4):369–S3. doi:10.1038/ng.2213
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2213&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22426310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

25. 25.Wen X, Pique-Regi R, Luca F. Integrating molecular QTL data into genome-wide genetic association analysis: Probabilistic assessment of enrichment and colocalization. PLOS Genet. 2017;13(3):e1006646. doi:10.1371/journal.pgen.1006646
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1006646&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28278150&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

26. 26.Hukku A, Sampson MG, Luca F, Pique-Regi R, Wen X. Analyzing and reconciling colocalization and transcriptome-wide association studies from the perspective of inferential reproducibility. Am J Hum Genet. 2022;109(5):825–837. doi:10.1016/j.ajhg.2022.04.005
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2022.04.005&link_type=DOI) 

27. 27.Berisa T, Pickrell JK. Approximately independent linkage disequilibrium blocks in human populations. Bioinformatics. 2016;32(2):283–285. doi:10.1093/bioinformatics/btv546
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btv546&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26395773&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

28. 28.Chen EY, Tan CM, Kou Y, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14(1):128. doi:10.1186/1471-2105-14-128
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-14-128&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23586463&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

29. 29.Xie Z, Bailey A, Kuleshov MV, et al. Gene Set Knowledge Discovery with Enrichr. Curr Protoc. 2021;1(3):e90. doi:10.1002/cpz1.90
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cpz1.90&link_type=DOI) 

30. 30.Balk SP, Ko YJ, Bubley GJ. Biology of Prostate-Specific Antigen. J Clin Oncol. 2003;21(2):383–391. doi:10.1200/JCO.2003.02.083
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjg6IjIxLzIvMzgzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDUvMDUvMjAyMy4wNS4wNC4yMzI4OTUyNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

31. 31.Barbeira AN, Pividori M, Zheng J, Wheeler HE, Nicolae DL, Im HK. Integrating predicted transcriptome from multiple tissues improves association detection. PLOS Genet. 2019;15(1):e1007889. doi:10.1371/journal.pgen.1007889
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1007889&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29739930&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

32. 32.Tindall EA, Hoang HN, Southey MC, et al. The 4q27 locus and prostate cancer risk. BMC Cancer. 2010;10(1):69. doi:10.1186/1471-2407-10-69
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2407-10-69&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20184734&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

33. 33.Yoshino S, Hara T, Weng JS, Takahashi Y, Seiki M, Sakamoto T. Genetic Screening of New Genes Responsible for Cellular Adaptation to Hypoxia Using a Genome-Wide shRNA Library. PLoS ONE. 2012;7(4):e35590. doi:10.1371/journal.pone.0035590
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0035590&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22523603&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

34. 34.Yoshino S, Matsui Y, Fukui Y, et al. EXOSC9 depletion attenuates P-body formation, stress resistance, and tumorigenicity of cancer cells. Sci Rep. 2020;10(1):9275. doi:10.1038/s41598-020-66455-2
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-66455-2&link_type=DOI) 

35. 35.Burns DT, Donkervoort S, Müller JS, et al. Variants in EXOSC9 Disrupt the RNA Exosome and Result in Cerebellar Atrophy with Spinal Motor Neuronopathy. Am J Hum Genet. 2018;102(5):858–873. doi:10.1016/j.ajhg.2018.03.011
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2018.03.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29727687&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

36. 36.O’Rourke DJ, DiJohnson DA, Caiazzo RJ, et al. Autoantibody Signatures as Biomarkers to Distinguish Prostate Cancer from Benign Prostatic Hyperplasia in Patients with Increased Serum Prostate Specific Antigen. Clin Chim Acta Int J Clin Chem. 2012;413(5-6):561–567. doi:10.1016/j.cca.2011.11.027
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cca.2011.11.027&link_type=DOI) 

37. 37.Feng T, Wei D, Li Q, et al. Four Novel Prognostic Genes Related to Prostate Cancer Identified Using Co-expression Structure Network Analysis. Front Genet. 2021;12. doi: 10.3389/fgene.2021.584164
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fgene.2021.584164&link_type=DOI) 

38. 38.Yang R, Du Y, Wang L, Chen Z, Liu X. Weighted gene co-expression network analysis identifies CCNA2 as a treatment target of prostate cancer through inhibiting cell cycle. J Cancer. 2020;11(5):1203–1211. doi:10.7150/jca.38173
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7150/jca.38173&link_type=DOI) 

39. 39.Smale ST, Kadonaga JT. The RNA polymerase II core promoter. Annu Rev Biochem. 2003;72:449–479. doi:10.1146/annurev.biochem.72.121801.161520
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.biochem.72.121801.161520&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12651739&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185092500015&link_type=ISI) 

40. 40.Wagner SD, Yakovchuk P, Gilman B, et al. RNA polymerase II acts as an RNA-dependent RNA polymerase to extend and destabilize a non-coding RNA. EMBO J. 2013;32(6):781–790. doi:10.1038/emboj.2013.18
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZW1ib2pubCI7czo1OiJyZXNpZCI7czo4OiIzMi82Lzc4MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA1LzA1LzIwMjMuMDUuMDQuMjMyODk1MjYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

41. 41.Schier AC, Taatjes DJ. Structure and mechanism of the RNA polymerase II transcription machinery. Genes Dev. 2020;34(7-8):465–488. doi:10.1101/gad.335679.119
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXNkZXYiO3M6NToicmVzaWQiO3M6MTA6IjM0LzctOC80NjUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNS8wNS8yMDIzLjA1LjA0LjIzMjg5NTI2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

42. 42.Louie MC, Yang HQ, Ma AH, et al. Androgen-induced recruitment of RNA polymerase II to a nuclear receptor–p160 coactivator complex. Proc Natl Acad Sci. 2003;100(5):2226–2230. doi:10.1073/pnas.0437824100
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTAwLzUvMjIyNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA1LzA1LzIwMjMuMDUuMDQuMjMyODk1MjYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

43. 43.Zhang X, Wang W, Zhu W, et al. Mechanisms and Functions of Long Non-Coding RNAs at Multiple Regulatory Levels. Int J Mol Sci. 2019;20(22):5573. doi:10.3390/ijms20225573
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijms20225573&link_type=DOI) 

44. 44.Misawa A, Takayama K, Inoue S. Long non-coding RNAs and prostate cancer. Cancer Sci. 2017;108(11):2107–2114. doi:10.1111/cas.13352
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/cas.13352&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28796922&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

45. 45.Zhang, J., Chang, Y., Xia, H. et al. HIST1H2BN induced cell proliferation and EMT phenotype in prostate cancer via NF-κB signal pathway. Genes Genom. 2021;43:1361–1369. doi: 10.1007/s13258-021-01164-2
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s13258-021-01164-2&link_type=DOI) 

46. 46.Liu Q, Reed M, Zhu H, et al. Epigenome-wide DNA methylation and transcriptome profiling of localized and locally advanced prostate cancer: Uncovering new molecular markers. Genomics. 2022;114(5):110474. doi:10.1016/j.ygeno.2022.110474
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ygeno.2022.110474&link_type=DOI) 

47. 47.Pressinotti NC, Klocker H, Schäfer G, et al. Differential expression of apoptotic genes PDIA3 and MAP3K5 distinguishes between low- and high-risk prostate cancer. Mol Cancer. 2009;8:130. doi:10.1186/1476-4598-8-130
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1476-4598-8-130&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20035634&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

48. 48.Mahajan NP, Liu Y, Majumder S, et al. Activated Cdc42-associated kinase Ack1 promotes prostate cancer progression via androgen receptor tyrosine phosphorylation. Proc Natl Acad Sci U S A. 2007;104(20):8438–8443. doi:10.1073/pnas.0700420104
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTA0LzIwLzg0MzgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNS8wNS8yMDIzLjA1LjA0LjIzMjg5NTI2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

49. 49.Davalieva K, Kiprijanovska S, Komina S, Petrusevska G, Zografska NC, Polenakovic M. Proteomics analysis of urine reveals acute phase response proteins as candidate diagnostic biomarkers for prostate cancer. Proteome Sci. 2015;13:2. doi:10.1186/s12953-014-0059-9
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12953-014-0059-9&link_type=DOI) 

50. 50.Krishnan N, Titus MA, Thapar R. The Prolyl Isomerase Pin1 Regulates mRNA Levels of Genes with Short Half-Lives by Targeting Specific RNA Binding Proteins. PLOS ONE. 2014;9(1):e85427. doi:10.1371/journal.pone.0085427
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0085427&link_type=DOI) 

51. 51.Wang Y, Pei X, Xu P, et al. E2F7, regulated by miRl7l30c, inhibits apoptosis and promotes cell cycle of prostate cancer cells. Oncol Rep. 2020;44(3):849–862. doi:10.3892/or.2020.7659
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3892/or.2020.7659&link_type=DOI) 

52. 52.Sheeba JDJ, Hegde S, Tamboli N, Nadig N, Keshavamurthy R, Ranganathan P. Gene expression signature of castrate resistant prostate cancer. BioRxiv. 2022.03.16.484397. doi:10.1101/2022.03.16.484397
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMi4wMy4xNi40ODQzOTd2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA1LzA1LzIwMjMuMDUuMDQuMjMyODk1MjYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

53. 53.Grasso CS, Wu YM, Robinson DR, et al. The mutational landscape of lethal castration-resistant prostate cancer. Nature. 2012;487(7406):239–243. doi:10.1038/nature11125
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature11125&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22722839&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000306278900041&link_type=ISI) 

54. 54.Lew QJ, Chu KL, Chia YL, Cheong N, Chao SH. HEXIM1, a New Player in the p53 Pathway. Cancers. 2013;5(3):838–856. doi:10.3390/cancers5030838
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/cancers5030838&link_type=DOI) 

55. 55.Chen R, Yik JHN, Lew QJ, Chao SH. Brd4 and HEXIM1: Multiple Roles in P-TEFb Regulation and Cancer. BioMed Res Int. 2014;2014:e232870. doi:10.1155/2014/232870
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1155/2014/232870&link_type=DOI) 

56. 56.Chen CD, Sawyers CL. NF-κB Activates Prostate-Specific Antigen Expression and Is Upregulated in Androgen-Independent Prostate Cancer. Mol Cell Biol. 2002;22(8):2862–2870. doi:10.1128/MCB.22.8.2862-2870.2002
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoibWNiIjtzOjU6InJlc2lkIjtzOjk6IjIyLzgvMjg2MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA1LzA1LzIwMjMuMDUuMDQuMjMyODk1MjYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

57. 57.Hara Y, Yamagata K, Oguchi K, Baba T. Nuclear localization of profilin III-ArpM1 complex in mouse spermiogenesis. FEBS Lett. 2008;582(20):2998–3004. doi:10.1016/j.febslet.2008.07.058
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.febslet.2008.07.058&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18692047&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000259123300005&link_type=ISI) 

58. 58.Umer N, Arévalo L, Phadke S, et al. Loss of Profilin3 Impairs Spermiogenesis by Affecting Acrosome Biogenesis, Autophagy, Manchette Development and Mitochondrial Organization. Front Cell Dev Biol. 2021;9 9:749559. doi: 10.3389/fcell.2021.749559
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fcell.2021.749559&link_type=DOI) 

59. 59.Song F, Amos CI, Lee JE, et al. Identification of a melanoma susceptibility locus and somatic mutation in TET2. Carcinogenesis. 2014;35(9):2097–2101. doi:10.1093/carcin/bgu140
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/carcin/bgu140&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24980573&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

60. 60.Law PJ, Timofeeva M, Fernandez-Rozadilla C, et al. Association analyses identify 31 new risk loci for colorectal cancer susceptibility. Nat Commun. 2019;10(1):2154. doi:10.1038/s41467-019-09775-w
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-09775-w&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

61. 61.McKay JD, Hung RJ, Han Y, et al. Large-scale association analysis identifies new lung cancer susceptibility loci and heterogeneity in genetic susceptibility across histological subtypes. Nat Genet. 2017;49(7):1126–1132. doi:10.1038/ng.3892
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3892&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28604730&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

62. 62.THE GTEX CONSORTIUM. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science. 2020;369(6509):1318–1330. doi:10.1126/science.aaz1776
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjkvNjUwOS8xMzE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDUvMDUvMjAyMy4wNS4wNC4yMzI4OTUyNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

63. 63.Martin AR, Kanai M, Kamatani Y, Okada Y, Neale BM, Daly MJ. Clinical use of current polygenic risk scores may exacerbate health disparities. Nat Genet. 2019;51(4):584–591. doi:10.1038/s41588-019-0379-x
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0379-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30926966&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

64. 64.Wainberg M, Sinnott-Armstrong N, Mancuso N, et al. Opportunities and challenges for transcriptome-wide association studies. Nat Genet. 2019;51(4):592–599. doi:10.1038/s41588-019-0385-z
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0385-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30926968&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom) 

65. 65.Liang, Y., Pividori, M., Manichaikul, A. et al. Polygenic transcriptome risk scores (PTRS) can improve portability of polygenic risk scores across ancestries. Genome Biol 2022;23(23). doi: 10.1186/s13059-021-02591-w
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-021-02591-w&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35027082&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F05%2F2023.05.04.23289526.atom)