Kuura-An automated workflow for analyzing WES and WGS data.

The advent of high-throughput sequencing technologies has revolutionized the field of genomic sciences by cutting down the cost and time associated with standard sequencing methods. This advancement has not only provided the research community with an abundance of data but has also presented the challenge of analyzing it. The paramount challenge in analyzing the copious amount of data is in using the optimal resources in terms of available tools. To address this research gap, we propose "Kuura-An automated workflow for analyzing WES and WGS data", which is optimized for both whole exome and whole genome sequencing data. This workflow is based on the nextflow pipeline scripting language and uses docker to manage and deploy the workflow. The workflow consists of four analysis stages-quality control, mapping to reference genome & quality score recalibration, variant calling & variant recalibration and variant consensus & annotation. An important feature of the DNA-seq workflow is that it uses the combination of multiple variant callers (GATK Haplotypecaller, DeepVariant, VarScan2, Freebayes and Strelka2), generating a list of high-confidence variants in a consensus call file. The workflow is flexible as it integrates the fragmented tools and can be easily extended by adding or updating tools or amending the parameters list. The use of a single parameters file enhances reproducibility of the results. The ease of deployment and usage of the workflow further increases computational reproducibility providing researchers with a standardized tool for the variant calling step in different projects. The source code, instructions for installation and use of the tool are publicly available at our github repository https://github.com/dhanaprakashj/kuura_pipeline.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants.

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

NTHL1 is a recessive cancer susceptibility gene.

In search of novel breast cancer (BC) risk variants, we performed a whole-exome sequencing and variant analysis of 69 Finnish BC patients as well as analysed loss-of-function variants identified in DNA repair genes in the Finns from the Genome Aggregation Database. Additionally, we carried out a validation study of SERPINA3 c.918-1G>C, recently suggested for BC predisposition. We estimated the frequencies of 41 rare candidate variants in 38 genes by genotyping them in 2482-4101 BC patients and in 1273-3985 controls. We further evaluated all coding variants in the candidate genes in a dataset of 18,786 BC patients and 182,927 controls from FinnGen. None of the variants associated significantly with cancer risk in the primary BC series; however, in the FinnGen data, NTHL1 c.244C>T p.(Gln82Ter) associated with BC with a high risk for homozygous (OR = 44.7 [95% CI 6.90-290], P = 6.7 × 10-5) and a low risk for heterozygous women (OR = 1.39 [1.18-1.64], P = 7.8 × 10-5). Furthermore, the results suggested a high risk of colorectal, urinary tract, and basal-cell skin cancer for homozygous individuals, supporting NTHL1 as a recessive multi-tumour susceptibility gene. No significant association with BC risk was detected for SERPINA3 or any other evaluated gene.

Evaluating Approaches for Constructing Polygenic Risk Scores for Prostate Cancer in Men of African and European Ancestry.

Genome-wide polygenic risk scores (GW-PRS) have been reported to have better predictive ability than PRS based on genome-wide significance thresholds across numerous traits. We compared the predictive ability of several GW-PRS approaches to a recently developed PRS of 269 established prostate cancer risk variants from multi-ancestry GWAS and fine-mapping studies (PRS 269 ). GW-PRS models were trained using a large and diverse prostate cancer GWAS of 107,247 cases and 127,006 controls used to develop the multi-ancestry PRS 269 . Resulting models were independently tested in 1,586 cases and 1,047 controls of African ancestry from the California/Uganda Study and 8,046 cases and 191,825 controls of European ancestry from the UK Biobank and further validated in 13,643 cases and 210,214 controls of European ancestry and 6,353 cases and 53,362 controls of African ancestry from the Million Veteran Program. In the testing data, the best performing GW-PRS approach had AUCs of 0.656 (95% CI=0.635-0.677) in African and 0.844 (95% CI=0.840-0.848) in European ancestry men and corresponding prostate cancer OR of 1.83 (95% CI=1.67-2.00) and 2.19 (95% CI=2.14-2.25), respectively, for each SD unit increase in the GW-PRS. However, compared to the GW-PRS, in African and European ancestry men, the PRS 269 had larger or similar AUCs (AUC=0.679, 95% CI=0.659-0.700 and AUC=0.845, 95% CI=0.841-0.849, respectively) and comparable prostate cancer OR (OR=2.05, 95% CI=1.87-2.26 and OR=2.21, 95% CI=2.16-2.26, respectively). Findings were similar in the validation data. This investigation suggests that current GW-PRS approaches may not improve the ability to predict prostate cancer risk compared to the multi-ancestry PRS 269 constructed with fine-mapping.

Evaluating approaches for constructing polygenic risk scores for prostate cancer in men of African and European ancestry.

Genome-wide polygenic risk scores (GW-PRSs) have been reported to have better predictive ability than PRSs based on genome-wide significance thresholds across numerous traits. We compared the predictive ability of several GW-PRS approaches to a recently developed PRS of 269 established prostate cancer-risk variants from multi-ancestry GWASs and fine-mapping studies (PRS269). GW-PRS models were trained with a large and diverse prostate cancer GWAS of 107,247 cases and 127,006 controls that we previously used to develop the multi-ancestry PRS269. Resulting models were independently tested in 1,586 cases and 1,047 controls of African ancestry from the California Uganda Study and 8,046 cases and 191,825 controls of European ancestry from the UK Biobank and further validated in 13,643 cases and 210,214 controls of European ancestry and 6,353 cases and 53,362 controls of African ancestry from the Million Veteran Program. In the testing data, the best performing GW-PRS approach had AUCs of 0.656 (95% CI = 0.635-0.677) in African and 0.844 (95% CI = 0.840-0.848) in European ancestry men and corresponding prostate cancer ORs of 1.83 (95% CI = 1.67-2.00) and 2.19 (95% CI = 2.14-2.25), respectively, for each SD unit increase in the GW-PRS. Compared to the GW-PRS, in African and European ancestry men, the PRS269 had larger or similar AUCs (AUC = 0.679, 95% CI = 0.659-0.700 and AUC = 0.845, 95% CI = 0.841-0.849, respectively) and comparable prostate cancer ORs (OR = 2.05, 95% CI = 1.87-2.26 and OR = 2.21, 95% CI = 2.16-2.26, respectively). Findings were similar in the validation studies. This investigation suggests that current GW-PRS approaches may not improve the ability to predict prostate cancer risk compared to the PRS269 developed from multi-ancestry GWASs and fine-mapping.

Biochemical activity induced by a germline variation in KLK3 (PSA) associates with cellular function and clinical outcome in prostate cancer.

Genetic variation at the 19q13.3 KLK locus is linked with prostate cancer susceptibility. The non-synonymous KLK3 SNP, rs17632542 (c.536T>C; Ile163Thr-substitution in PSA) is associated with reduced prostate cancer risk, however, the functional relevance is unknown. Here, we identify that the SNP variant-induced change in PSA biochemical activity as a previously undescribed function mediating prostate cancer pathogenesis. The 'Thr' PSA variant led to small subcutaneous tumours, supporting reduced prostate cancer risk. However, 'Thr' PSA also displayed higher metastatic potential with pronounced osteolytic activity in an experimental metastasis in-vivo model. Biochemical characterization of this PSA variant demonstrated markedly reduced proteolytic activity that correlated with differences in in-vivo tumour burden. The SNP is associated with increased risk for aggressive disease and prostate cancer-specific mortality in three independent cohorts, highlighting its critical function in mediating metastasis. Carriers of this SNP allele had reduced serum total PSA and a higher free/total PSA ratio that could contribute to late biopsy decisions and delay in diagnosis. Our results provide a molecular explanation for the prominent 19q13.3 KLK locus, rs17632542 SNP, association with a spectrum of prostate cancer clinical outcomes.

Transcripts of the Prostate Cancer-Associated Gene ANO7 Are Retained in the Nuclei of Prostatic Epithelial Cells.

Prostate cancer affects millions of men globally. The prostate cancer-associated gene ANO7 is downregulated in advanced prostate cancer, whereas benign tissue and low-grade cancer display varying expression levels. In this study, we assess the spatial correlation between ANO7 mRNA and protein using fluorescent in situ hybridization and immunohistochemistry for the detection of mRNA and protein in parallel sections of tissue microarrays prepared from radical prostatectomy samples. We show that ANO7 mRNA and protein expression correlate in prostate tissue. Furthermore, we show that ANO7 mRNA is enriched in the nuclei of the luminal cells at 89% in benign ducts and low-grade cancer, and at 78% in high-grade cancer. The nuclear enrichment of ANO7 mRNA was validated in prostate cancer cell lines 22Rv1 and MDA PCa 2b using droplet digital polymerase chain reaction (ddPCR) on RNA isolated from nuclear and cytoplasmic fractions of the cells. The nuclear enrichment of ANO7 mRNA was compared to the nuclearly-enriched lncRNA MALAT1, confirming the surprisingly high nuclear retention of ANO7 mRNA. ANO7 has been suggested to be used as a diagnostic marker and a target for immunotherapy, but a full comprehension of its role in prostate cancer progression is currently lacking. Our results contribute to a better understanding of the dynamics of ANO7 expression in prostatic tissue.

TREX1 p.A129fs and p.Y305C variants in a large multi-ethnic cohort of CADASIL-like unrelated patients.

Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) and retinal vasculopathy with cerebral leukodystrophy and systemic manifestations (RVCL-S) are the most common forms of rare monogenic early-onset cerebral small vessel disease and share clinical, and, to different extents, neuroradiological and neuropathological features. However, whether CADASIL and RVCL-S overlapping phenotype may be explained by shared genetic risk or causative factors such as TREX1 coding variants remains poorly understood. To investigate this intriguing hypothesis, we used exome sequencing to screen TREX1 protein-coding variability in a large multi-ethnic cohort of 180 early-onset independent familial and apparently sporadic CADASIL-like Caucasian patients from the USA, Portugal, Finland, Serbia and Turkey. We report 2 very rare and likely pathogenic TREX1 mutations: a loss of function mutation (p.Ala129fs) clustering in the catalytic domain, in an apparently sporadic 46-year-old patient from the USA and a missense mutation (p.Tyr305Cys) in the well conserved C-terminal region, in a 57-year-old patient with positive family history from Serbia. In concert with recent findings, our study expands the clinical spectrum of diseases associated with TREX1 mutations.

Clinical features and spectrum of NOTCH3 variants in Finnish patients with cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL).

OBJECTIVES: Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is a cerebral small vessel disease caused by pathogenic variants in the NOTCH3 gene. In Finland, the majority of CADASIL patients carry the pathogenic founder variant c.397C>T, (p.Arg133Cys), but the spectrum of other NOTCH3 variants has not been investigated previously. The aim of the study was to investigate the spectrum and prevalence of NOTCH3 variants Finnish CADASIL patients and to examine the clinical features associated with them. MATERIALS AND METHODS: The spectrum of NOTCH3 variants and the clinical features associated with them were retrospectively examined in 294 Finnish CADASIL patients tested during January 1996 to October 2021 in the Medical Genetics laboratory of Department of Genomics of Turku University Hospital, where practically all samples of patients with suspected CADASIL in Finland are investigated. RESULTS: The most common NOTCH3 variants in the study cohort were c.397C>T, (p.Arg133Cys) (68%) and c.3206A>G p.(Tyr1069Cys) (18%), but other less common NOTCH3 variants were detected in as many as 14% of the patients. Eight of the detected NOTCH3 variants were novel: c.520T>A,p.(Cys174Ser), c.836A>G,p.(Gln279Arg), c.1369T>G,p.(Cys457Gly), c.1338C>G,p.(Cys446Trp), c.1564T>G,p.(Cys522Gly), c.2848T>G,p.(Cys950Gly), c.6102dup,p.(Gly2035Argfs*60), and c.2410+6C>G. Other NOTCH3 variants than p.Arg133Cys and p.Tyr1069Cys were more often associated with more severe clinical features. CONCLUSION: This study revealed the genetic and clinical spectrum of CADASIL in the Finnish population. Sequencing of the whole NOTCH3 gene performing a gene-panel or exome sequencing is recommended when suspecting CADASIL.

Genetic analysis reveals novel variants for vascular cognitive impairment.

OBJECTIVES: The genetic background of vascular cognitive impairment (VCI) is poorly understood compared to other dementia disorders. The aim of the study was to investigate the genetic background of VCI in a well-characterized Finnish cohort. MATERIALS & METHODS: Whole-exome sequencing (WES) was applied in 45 Finnish VCI patients. Copy-number variant (CNV) analysis using a SNP array was performed in 80 VCI patients. This study also examined the prevalence of variants at the miR-29 binding site of COL4A1 in 73 Finnish VCI patients. RESULTS: In 40% (18/45) of the cases, WES detected possibly causative variants in genes associated with cerebral small vessel disease (CSVD) or other neurological or stroke-related disorders. These variants included HTRA1:c.847G>A p.(Gly283Arg), TREX1:c.1079A>G, p.(Tyr360Cys), COLGALT1:c.1411C>T, p.(Arg471Trp), PRNP: c.713C>T, p.(Pro238Leu), and MTHFR:c.1061G>C, p.(Gly354Ala). Additionally, screening of variants in the 3'UTR of COL4A1 gene in a sub-cohort of 73 VCI patients identified a novel variant c.*36T>A. CNV analysis showed that pathogenic CNVs are uncommon in VCI. CONCLUSIONS: These data support pathogenic roles of variants in HTRA1, TREX1 and in the 3'UTR of COL4A1 in CSVD and VCI, and suggest that vascular pathogenic mechanisms are linked to neurodegeneration, expanding the understanding of the genetic background of VCI.

Prostate cancer risk stratification improvement across multiple ancestries with new polygenic hazard score.

BACKGROUND: Prostate cancer risk stratification using single-nucleotide polymorphisms (SNPs) demonstrates considerable promise in men of European, Asian, and African genetic ancestries, but there is still need for increased accuracy. We evaluated whether including additional SNPs in a prostate cancer polygenic hazard score (PHS) would improve associations with clinically significant prostate cancer in multi-ancestry datasets. METHODS: In total, 299 SNPs previously associated with prostate cancer were evaluated for inclusion in a new PHS, using a LASSO-regularized Cox proportional hazards model in a training dataset of 72,181 men from the PRACTICAL Consortium. The PHS model was evaluated in four testing datasets: African ancestry, Asian ancestry, and two of European Ancestry-the Cohort of Swedish Men (COSM) and the ProtecT study. Hazard ratios (HRs) were estimated to compare men with high versus low PHS for association with clinically significant, with any, and with fatal prostate cancer. The impact of genetic risk stratification on the positive predictive value (PPV) of PSA testing for clinically significant prostate cancer was also measured. RESULTS: The final model (PHS290) had 290 SNPs with non-zero coefficients. Comparing, for example, the highest and lowest quintiles of PHS290, the hazard ratios (HRs) for clinically significant prostate cancer were 13.73 [95% CI: 12.43-15.16] in ProtecT, 7.07 [6.58-7.60] in African ancestry, 10.31 [9.58-11.11] in Asian ancestry, and 11.18 [10.34-12.09] in COSM. Similar results were seen for association with any and fatal prostate cancer. Without PHS stratification, the PPV of PSA testing for clinically significant prostate cancer in ProtecT was 0.12 (0.11-0.14). For the top 20% and top 5% of PHS290, the PPV of PSA testing was 0.19 (0.15-0.22) and 0.26 (0.19-0.33), respectively. CONCLUSIONS: We demonstrate better genetic risk stratification for clinically significant prostate cancer than prior versions of PHS in multi-ancestry datasets. This is promising for implementing precision-medicine approaches to prostate cancer screening decisions in diverse populations.

The variant rs77559646 associated with aggressive prostate cancer disrupts ANO7 mRNA splicing and protein expression.

Prostate cancer is among the most common cancers in men, with a large fraction of the individual risk attributable to heritable factors. A majority of the diagnosed cases does not lead to a lethal disease, and hence biological markers that can distinguish between indolent and fatal forms of the disease are of great importance for guiding treatment decisions. Although over 300 genetic variants are known to be associated with prostate cancer risk, few have been associated with the risk of an aggressive disease. One such variant is rs77559646 located in ANO7. This variant has a dual function. It constitutes a missense mutation in the short isoform of ANO7 and a splice region mutation in full-length ANO7. In this study, we have analyzed the impact of the variant allele of rs77559646 on ANO7 mRNA splicing using a minigene splicing assay and by performing splicing analysis with the tools IRFinder (intron retention finder), rMATS (replicate multivariate analysis of transcript splicing) and LeafCutter on RNA sequencing data from prostate tissue of six rs77559646 variant allele carriers and 43 non-carriers. The results revealed a severe disruption of ANO7 mRNA splicing in rs77559646 variant allele carriers. Immunohistochemical analysis of prostate samples from patients homozygous for the rs77559646 variant allele demonstrated a loss of apically localized ANO7 protein. Our study is the first to provide a mechanistic explanation for the impact of a prostate cancer risk SNP on ANO7 protein production. Furthermore, the rs77559646 variant is the first known germline loss-of-function mutation described for ANO7. We suggest that loss of ANO7 contributes to prostate cancer progression.

A Rare Variant in ERF (rs144812092) Predisposes to Prostate and Bladder Cancers in an Extended Pedigree.

Pairs of related bladder cancer cases who belong to pedigrees with an excess of bladder cancer were sequenced to identify rare, shared variants as candidate predisposition variants. Candidate variants were tested for association with bladder cancer risk. A validated variant was assayed for segregation to other related cancer cases, and the predicted protein structure of this variant was analyzed. This study of affected bladder cancer relative pairs from high-risk pedigrees identified 152 bladder cancer predisposition candidate variants. One variant in ERF (ETS Repressing Factor) was significantly associated with bladder cancer risk in an independent population, was observed to segregate with bladder and prostate cancer in relatives, and showed evidence for altering the function of the associated protein. This finding of a rare variant in ERF that is strongly associated with bladder and prostate cancer risk in an extended pedigree both validates ERF as a cancer predisposition gene and shows the continuing value of analyzing affected members of high-risk pedigrees to identify and validate rare cancer predisposition variants.

Novel prostate cancer susceptibility gene SP6 predisposes patients to aggressive disease.

Prostate cancer (PrCa) is one of the most common cancers in men, but little is known about factors affecting its clinical outcomes. Genome-wide association studies have identified more than 170 germline susceptibility loci, but most of them are not associated with aggressive disease. We performed a genome-wide analysis of 185,478 SNPs in Finnish samples (2738 cases, 2400 controls) from the international Collaborative Oncological Gene-Environment Study (iCOGS) to find underlying PrCa risk variants. We identified a total of 21 common, low-penetrance susceptibility loci, including 10 novel variants independently associated with PrCa risk. Novel risk loci were located in the 8q24 (CASC8 rs16902147, OR 1.86, padj = 3.53 × 10-8 and rs58809953, OR 1.71, padj = 4.00 × 10-6; intergenic rs79012498, OR 1.81, padj = 4.26 × 10-8), 17q21 (SP6 rs2074187, OR 1.66, padj = 3.75 × 10-5), 11q13 (rs12795301, OR 1.42, padj = 2.89 × 10-5) and 8p21 (rs995432, OR 1.38, padj = 3.00 × 10-11) regions. Here, we describe SP6, a transcription factor gene, as a new, potentially high-risk gene for PrCa. The intronic variant rs2074187 in SP6 was associated not only with overall susceptibility to PrCa (OR 1.66) but also with a higher odds ratio for aggressive PrCa (OR 1.89) and lower odds for non-aggressive PrCa (OR 1.43). Furthermore, the new intergenic variant rs79012498 at 8q24 conferred risk for aggressive PrCa. Our findings highlighted the power of a population-stratified approach to identify novel, clinically actionable germline PrCa risk loci and strongly suggested SP6 as a new PrCa candidate gene that may be involved in the pathogenesis of PrCa.

Multisite Evaluation and Validation of a Sensitive Diagnostic and Screening System for Spinal Muscular Atrophy that Reports SMN1 and SMN2 Copy Number, along with Disease Modifier and Gene Duplication Variants.

Spinal muscular atrophy is a severe autosomal recessive disease caused by disruptions in the SMN1 gene. The nearly identical SMN2 gene copy number is associated with disease severity. SMN1 duplication markers, such as c.∗3+80T>G and c.∗211_∗212del, can assess residual carrier risk. An SMN2 disease modifier (c.859G>C) can help inform prognostic outcomes. The emergence of multiple precision gene therapies for spinal muscular atrophy requires accurate and rapid detection of SMN1 and SMN2 copy numbers to enable early treatment and optimal patient outcomes. We developed and evaluated a single-tube PCR/capillary electrophoresis assay system that quantifies SMN1/2 copy numbers and genotypes three additional clinically relevant variants. Analytical validation was performed with human cell lines and whole blood representing varying SMN1/2 copies on four capillary electrophoresis instrument models. In addition, four independent laboratories used the assay to test 468 residual clinical genomic DNA samples. The results were ≥98.3% concordant with consensus SMN1/2 exon 7 copy numbers, determined using multiplex ligation-dependent probe amplification and droplet digital PCR, and were 100% concordant with Sanger sequencing for the three variants. Furthermore, copy number values were 98.6% (SMN1) and 97.1% (SMN2) concordant to each laboratory's own reference results.

KLK3 SNP-SNP interactions for prediction of prostate cancer aggressiveness.

Risk classification for prostate cancer (PCa) aggressiveness and underlying mechanisms remain inadequate. Interactions between single nucleotide polymorphisms (SNPs) may provide a solution to fill these gaps. To identify SNP-SNP interactions in the four pathways (the angiogenesis-, mitochondria-, miRNA-, and androgen metabolism-related pathways) associated with PCa aggressiveness, we tested 8587 SNPs for 20,729 cases from the PCa consortium. We identified 3 KLK3 SNPs, and 1083 (P < 3.5 × 10-9) and 3145 (P < 1 × 10-5) SNP-SNP interaction pairs significantly associated with PCa aggressiveness. These SNP pairs associated with PCa aggressiveness were more significant than each of their constituent SNP individual effects. The majority (98.6%) of the 3145 pairs involved KLK3. The 3 most common gene-gene interactions were KLK3-COL4A1:COL4A2, KLK3-CDH13, and KLK3-TGFBR3. Predictions from the SNP interaction-based polygenic risk score based on 24 SNP pairs are promising. The prevalence of PCa aggressiveness was 49.8%, 21.9%, and 7.0% for the PCa cases from our cohort with the top 1%, middle 50%, and bottom 1% risk profiles. Potential biological functions of the identified KLK3 SNP-SNP interactions were supported by gene expression and protein-protein interaction results. Our findings suggest KLK3 SNP interactions may play an important role in PCa aggressiveness.

PHACTR1 genetic variability is not critical in small vessel ischemic disease patients and PcomA recruitment in C57BL/6J mice.

Recently, several genome-wide association studies identified PHACTR1 as key locus for five diverse vascular disorders: coronary artery disease, migraine, fibromuscular dysplasia, cervical artery dissection and hypertension. Although these represent significant risk factors or comorbidities for ischemic stroke, PHACTR1 role in brain small vessel ischemic disease and ischemic stroke most important survival mechanism, such as the recruitment of brain collateral arteries like posterior communicating arteries (PcomAs), remains unknown. Therefore, we applied exome and genome sequencing in a multi-ethnic cohort of 180 early-onset independent familial and apparently sporadic brain small vessel ischemic disease and CADASIL-like Caucasian patients from US, Portugal, Finland, Serbia and Turkey and in 2 C57BL/6J stroke mouse models (bilateral common carotid artery stenosis [BCCAS] and middle cerebral artery occlusion [MCAO]), characterized by different degrees of PcomAs patency. We report 3 very rare coding variants in the small vessel ischemic disease-CADASIL-like cohort (p.Glu198Gln, p.Arg204Gly, p.Val251Leu) and a stop-gain mutation (p.Gln273*) in one MCAO mouse. These coding variants do not cluster in PHACTR1 known pathogenic domains and are not likely to play a critical role in small vessel ischemic disease or brain collateral circulation. We also exclude the possibility that copy number variants (CNVs) or a variant enrichment in Phactr1 may be associated with PcomA recruitment in BCCAS mice or linked to diverse vascular traits (cerebral blood flow pre-surgery, PcomA size, leptomeningeal microcollateral length and junction density during brain hypoperfusion) in C57BL/6J mice, respectively. Genetic variability in PHACTR1 is not likely to be a common susceptibility factor influencing small vessel ischemic disease in patients and PcomA recruitment in C57BL/6J mice. Nonetheless, rare variants in PHACTR1 RPEL domains may influence the stroke outcome and are worth investigating in a larger cohort of small vessel ischemic disease patients, different ischemic stroke subtypes and with functional studies.

Polygenic hazard score is associated with prostate cancer in multi-ethnic populations.

Genetic models for cancer have been evaluated using almost exclusively European data, which could exacerbate health disparities. A polygenic hazard score (PHS1) is associated with age at prostate cancer diagnosis and improves screening accuracy in Europeans. Here, we evaluate performance of PHS2 (PHS1, adapted for OncoArray) in a multi-ethnic dataset of 80,491 men (49,916 cases, 30,575 controls). PHS2 is associated with age at diagnosis of any and aggressive (Gleason score ≥ 7, stage T3-T4, PSA ≥ 10 ng/mL, or nodal/distant metastasis) cancer and prostate-cancer-specific death. Associations with cancer are significant within European (n = 71,856), Asian (n = 2,382), and African (n = 6,253) genetic ancestries (p < 10-180). Comparing the 80th/20th PHS2 percentiles, hazard ratios for prostate cancer, aggressive cancer, and prostate-cancer-specific death are 5.32, 5.88, and 5.68, respectively. Within European, Asian, and African ancestries, hazard ratios for prostate cancer are: 5.54, 4.49, and 2.54, respectively. PHS2 risk-stratifies men for any, aggressive, and fatal prostate cancer in a multi-ethnic dataset.

Additional SNPs improve risk stratification of a polygenic hazard score for prostate cancer.

BACKGROUND: Polygenic hazard scores (PHS) can identify individuals with increased risk of prostate cancer. We estimated the benefit of additional SNPs on performance of a previously validated PHS (PHS46). MATERIALS AND METHOD: 180 SNPs, shown to be previously associated with prostate cancer, were used to develop a PHS model in men with European ancestry. A machine-learning approach, LASSO-regularized Cox regression, was used to select SNPs and to estimate their coefficients in the training set (75,596 men). Performance of the resulting model was evaluated in the testing/validation set (6,411 men) with two metrics: (1) hazard ratios (HRs) and (2) positive predictive value (PPV) of prostate-specific antigen (PSA) testing. HRs were estimated between individuals with PHS in the top 5% to those in the middle 40% (HR95/50), top 20% to bottom 20% (HR80/20), and bottom 20% to middle 40% (HR20/50). PPV was calculated for the top 20% (PPV80) and top 5% (PPV95) of PHS as the fraction of individuals with elevated PSA that were diagnosed with clinically significant prostate cancer on biopsy. RESULTS: 166 SNPs had non-zero coefficients in the Cox model (PHS166). All HR metrics showed significant improvements for PHS166 compared to PHS46: HR95/50 increased from 3.72 to 5.09, HR80/20 increased from 6.12 to 9.45, and HR20/50 decreased from 0.41 to 0.34. By contrast, no significant differences were observed in PPV of PSA testing for clinically significant prostate cancer. CONCLUSIONS: Incorporating 120 additional SNPs (PHS166 vs PHS46) significantly improved HRs for prostate cancer, while PPV of PSA testing remained the same.