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.

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.

Prophylactic Anticoagulation With Intermediate-Dose Certoparin in Vascular-Risk Pregnancies-The PACER-VARP Registry.

The management of pregnant women at increased risk of thromboembolic/other vascular events is still a matter of debate. In a single-center, retrospective, observational trial, we analyzed the safety and efficacy of prophylactic anticoagulation with certoparin in pregnant women at intermediate- or high-risk by EThIG criteria of thromboembolic/other vascular events. Subcutaneous certoparin 8,000 IU once daily was administered immediately after pregnancy confirmation and continued for 6 weeks postpartum. We investigated 74 pregnancies (49 women; mean age 31.8 years; weight 77.3 kg). Most prevalent risk factors were factor V Leiden mutation (40.5%), thrombogenic factor II mutation (12.2%) and protein S deficiency (8.1%). In 76 control pregnancies prior to registry inclusion/without anticoagulation there were 14 cases [18.4%] of venous thromboembolism (between week 7 gestation and week 8 postpartum); 63.2% pregnancies resulted in abortion (median week 8.6 gestation). With certoparin anticoagulation, thromboembolism was 1.4%, exclusively non-major bleeding was 4.1% and abortion was 10.8%. One case of pre-eclampsia necessitating obstetric intervention occurred. Prophylactic anticoagulation with intermediate-dose certoparin throughout pregnancies at increased venous vascular risk was safe and effective.

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.

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.

The CHEK2 Variant C.349A>G Is Associated with Prostate Cancer Risk and Carriers Share a Common Ancestor.

The identification of recurrent founder variants in cancer predisposing genes may have important implications for implementing cost-effective targeted genetic screening strategies. In this study, we evaluated the prevalence and relative risk of the CHEK2 recurrent variant c.349A>G in a series of 462 Portuguese patients with early-onset and/or familial/hereditary prostate cancer (PrCa), as well as in the large multicentre PRACTICAL case-control study comprising 55,162 prostate cancer cases and 36,147 controls. Additionally, we investigated the potential shared ancestry of the carriers by performing identity-by-descent, haplotype and age estimation analyses using high-density SNP data from 70 variant carriers belonging to 11 different populations included in the PRACTICAL consortium. The CHEK2 missense variant c.349A>G was found significantly associated with an increased risk for PrCa (OR 1.9; 95% CI: 1.1-3.2). A shared haplotype flanking the variant in all carriers was identified, strongly suggesting a common founder of European origin. Additionally, using two independent statistical algorithms, implemented by DMLE+2.3 and ESTIAGE, we were able to estimate the age of the variant between 2300 and 3125 years. By extending the haplotype analysis to 14 additional carrier families, a shared core haplotype was revealed among all carriers matching the conserved region previously identified in the high-density SNP analysis. These findings are consistent with CHEK2 c.349A>G being a founder variant associated with increased PrCa risk, suggesting its potential usefulness for cost-effective targeted genetic screening in PrCa families.

A Genetic Risk Score to Personalize Prostate Cancer Screening, Applied to Population Data.

BACKGROUND: A polygenic hazard score (PHS), the weighted sum of 54 SNP genotypes, was previously validated for association with clinically significant prostate cancer and for improved prostate cancer screening accuracy. Here, we assess the potential impact of PHS-informed screening. METHODS: United Kingdom population incidence data (Cancer Research United Kingdom) and data from the Cluster Randomized Trial of PSA Testing for Prostate Cancer were combined to estimate age-specific clinically significant prostate cancer incidence (Gleason score ≥7, stage T3-T4, PSA ≥10, or nodal/distant metastases). Using HRs estimated from the ProtecT prostate cancer trial, age-specific incidence rates were calculated for various PHS risk percentiles. Risk-equivalent age, when someone with a given PHS percentile has prostate cancer risk equivalent to an average 50-year-old man (50-year-standard risk), was derived from PHS and incidence data. Positive predictive value (PPV) of PSA testing for clinically significant prostate cancer was calculated using PHS-adjusted age groups. RESULTS: The expected age at diagnosis of clinically significant prostate cancer differs by 19 years between the 1st and 99th PHS percentiles: men with PHS in the 1st and 99th percentiles reach the 50-year-standard risk level at ages 60 and 41, respectively. PPV of PSA was higher for men with higher PHS-adjusted age. CONCLUSIONS: PHS provides individualized estimates of risk-equivalent age for clinically significant prostate cancer. Screening initiation could be adjusted by a man's PHS. IMPACT: Personalized genetic risk assessments could inform prostate cancer screening decisions.

The effect of sample size on polygenic hazard models for prostate cancer.

We determined the effect of sample size on performance of polygenic hazard score (PHS) models in prostate cancer. Age and genotypes were obtained for 40,861 men from the PRACTICAL consortium. The dataset included 201,590 SNPs per subject, and was split into training and testing sets. Established-SNP models considered 65 SNPs that had been previously associated with prostate cancer. Discovery-SNP models used stepwise selection to identify new SNPs. The performance of each PHS model was calculated for random sizes of the training set. The performance of a representative Established-SNP model was estimated for random sizes of the testing set. Mean HR98/50 (hazard ratio of top 2% to average in test set) of the Established-SNP model increased from 1.73 [95% CI: 1.69-1.77] to 2.41 [2.40-2.43] when the number of training samples was increased from 1 thousand to 30 thousand. Corresponding HR98/50 of the Discovery-SNP model increased from 1.05 [0.93-1.18] to 2.19 [2.16-2.23]. HR98/50 of a representative Established-SNP model using testing set sample sizes of 0.6 thousand and 6 thousand observations were 1.78 [1.70-1.85] and 1.73 [1.71-1.76], respectively. We estimate that a study population of 20 thousand men is required to develop Discovery-SNP PHS models while 10 thousand men should be sufficient for Established-SNP models.

Detecting predictive androgen receptor modifications in circulating prostate cancer cells.

Molecular modifications of the androgen receptor (AR) can cause resistance to androgen deprivation therapy (ADT) in prostate cancer patients. Since lack of representative tumor samples hinders therapy adjustments according to emerging AR-modifications, we evaluated simultaneous detection of the two most common AR modifications (AR-V7 splice variant and AR point mutations) in circulating tumor cells (CTCs). We devised a single-tube assay to detect AR-V7 splice variants and AR point mutations in CTCs using immunomagnetic cell isolation, followed by quantitative real-time PCR and DNA pyrosequencing. We prospectively investigated 47 patients with PSA progression awaiting therapy switch. Comparison of response to newly administered therapy and CTC-AR-status allowed effect size estimation. Nineteen (51%) of 37 patients with detectable CTCs carried AR-modifications. Seventeen patients carried the AR-V7 splice variant, one harbored a p.T878A point mutation and one harbored both AR-V7 and a p.H875Y mutation. We estimated a positive predictive value for response and non-response to therapy by AR status in CTCs of ~94%. Based on a conservative calculation, we estimated the effect size for molecularly-informed therapy switches for prospective clinical trial planning to ~27%. In summary, the ability to determine key resistance-mediating AR modifications in CTCs has the potential to considerably improve prostate cancer treatment.

PU.1 controls fibroblast polarization and tissue fibrosis.

Fibroblasts are polymorphic cells with pleiotropic roles in organ morphogenesis, tissue homeostasis and immune responses. In fibrotic diseases, fibroblasts synthesize abundant amounts of extracellular matrix, which induces scarring and organ failure. By contrast, a hallmark feature of fibroblasts in arthritis is degradation of the extracellular matrix because of the release of metalloproteinases and degrading enzymes, and subsequent tissue destruction. The mechanisms that drive these functionally opposing pro-fibrotic and pro-inflammatory phenotypes of fibroblasts remain unknown. Here we identify the transcription factor PU.1 as an essential regulator of the pro-fibrotic gene expression program. The interplay between transcriptional and post-transcriptional mechanisms that normally control the expression of PU.1 expression is perturbed in various fibrotic diseases, resulting in the upregulation of PU.1, induction of fibrosis-associated gene sets and a phenotypic switch in extracellular matrix-producing pro-fibrotic fibroblasts. By contrast, pharmacological and genetic inactivation of PU.1 disrupts the fibrotic network and enables reprogramming of fibrotic fibroblasts into resting fibroblasts, leading to regression of fibrosis in several organs.

AA9int: SNP interaction pattern search using non-hierarchical additive model set.

Motivation: The use of single nucleotide polymorphism (SNP) interactions to predict complex diseases is getting more attention during the past decade, but related statistical methods are still immature. We previously proposed the SNP Interaction Pattern Identifier (SIPI) approach to evaluate 45 SNP interaction patterns/patterns. SIPI is statistically powerful but suffers from a large computation burden. For large-scale studies, it is necessary to use a powerful and computation-efficient method. The objective of this study is to develop an evidence-based mini-version of SIPI as the screening tool or solitary use and to evaluate the impact of inheritance mode and model structure on detecting SNP-SNP interactions. Results: We tested two candidate approaches: the 'Five-Full' and 'AA9int' method. The Five-Full approach is composed of the five full interaction models considering three inheritance modes (additive, dominant and recessive). The AA9int approach is composed of nine interaction models by considering non-hierarchical model structure and the additive mode. Our simulation results show that AA9int has similar statistical power compared to SIPI and is superior to the Five-Full approach, and the impact of the non-hierarchical model structure is greater than that of the inheritance mode in detecting SNP-SNP interactions. In summary, it is recommended that AA9int is a powerful tool to be used either alone or as the screening stage of a two-stage approach (AA9int+SIPI) for detecting SNP-SNP interactions in large-scale studies. Availability and implementation: The 'AA9int' and 'parAA9int' functions (standard and parallel computing version) are added in the SIPI R package, which is freely available at https://linhuiyi.github.io/LinHY_Software/. Supplementary information: Supplementary data are available at Bioinformatics online.

Gene and pathway level analyses of germline DNA-repair gene variants and prostate cancer susceptibility using the iCOGS-genotyping array.

This corrects the article DOI: 10.1038/bjc.2016.50.

Polygenic hazard score to guide screening for aggressive prostate cancer: development and validation in large scale cohorts.

OBJECTIVES: To develop and validate a genetic tool to predict age of onset of aggressive prostate cancer (PCa) and to guide decisions of who to screen and at what age. DESIGN: Analysis of genotype, PCa status, and age to select single nucleotide polymorphisms (SNPs) associated with diagnosis. These polymorphisms were incorporated into a survival analysis to estimate their effects on age at diagnosis of aggressive PCa (that is, not eligible for surveillance according to National Comprehensive Cancer Network guidelines; any of Gleason score ≥7, stage T3-T4, PSA (prostate specific antigen) concentration ≥10 ng/L, nodal metastasis, distant metastasis). The resulting polygenic hazard score is an assessment of individual genetic risk. The final model was applied to an independent dataset containing genotype and PSA screening data. The hazard score was calculated for these men to test prediction of survival free from PCa. SETTING: Multiple institutions that were members of international PRACTICAL consortium. PARTICIPANTS: All consortium participants of European ancestry with known age, PCa status, and quality assured custom (iCOGS) array genotype data. The development dataset comprised 31 747 men; the validation dataset comprised 6411 men. MAIN OUTCOME MEASURES: Prediction with hazard score of age of onset of aggressive cancer in validation set. RESULTS: In the independent validation set, the hazard score calculated from 54 single nucleotide polymorphisms was a highly significant predictor of age at diagnosis of aggressive cancer (z=11.2, P<10-16). When men in the validation set with high scores (>98th centile) were compared with those with average scores (30th-70th centile), the hazard ratio for aggressive cancer was 2.9 (95% confidence interval 2.4 to 3.4). Inclusion of family history in a combined model did not improve prediction of onset of aggressive PCa (P=0.59), and polygenic hazard score performance remained high when family history was accounted for. Additionally, the positive predictive value of PSA screening for aggressive PCa was increased with increasing polygenic hazard score. CONCLUSIONS: Polygenic hazard scores can be used for personalised genetic risk estimates that can predict for age at onset of aggressive PCa.

A low-frequency haplotype spanning SLX4/FANCP constitutes a new risk locus for early-onset breast cancer (<60 years) and is associated with reduced DNA repair capacity.

Only a fraction of breast cancer (BC) cases can be yet explained by mutations in genes or genomic variants discovered in linkage, genome-wide association and sequencing studies. The known genes entailing medium or high risk for BC are strongly enriched for a function in DNA double strand repair. Thus, aiming at identifying low frequency variants conferring an intermediate risk, we here investigated 17 variants (MAF: 0.01-0.1) in 10 candidate genes involved in DNA repair or cell cycle control. In an exploration cohort of 437 cases and 1189 controls, we show the variant rs3810813 in the SLX4/FANCP gene to be significantly associated with both BC (≤60 years; OR = 2.6(1.6-3.9), p = 1.6E-05) and decreased DNA repair capacity (≤60 years; beta = 37.8(17.9-57.8), p = 5.3E-4). BC association was confirmed in a verification cohort (N = 2441). Both associations were absent from cases diagnosed >60 years and stronger the earlier the diagnosis. By imputation we show that rs3810813 tags a haplotype with 5 additional variants with the same allele frequency (R2 > 0.9), and a pattern of association very similar for both phenotypes (cases <60 years, p < 0.001, the Bonferroni threshold derived from unlinked variants in the region). In young cases (≤60 years) carrying the risk haplotype, micronucleus test results are predictive for BC (AUC > 0.9). Our findings propose a risk variant with high penetrance on the haplotype spanning SLX4/FANCP to be functionally associated to BC predisposition via decreased repair capacity and suggest this variant is carried by a fraction of these haplotypes that is enriched in early onset BC cases.

Nintedanib inhibits macrophage activation and ameliorates vascular and fibrotic manifestations in the Fra2 mouse model of systemic sclerosis.

BACKGROUND: Nintedanib is an inhibitor targeting platelet-derived growth factor receptor, fibroblast growth factor receptor and vascular endothelial growth factor receptor tyrosine kinases that has recently been approved for the treatment of idiopathic pulmonary fibrosis. The aim of this study was to analyse the effects of nintedanib in the fos-related antigen-2 (Fra2) mouse model of systemic sclerosis (SSc). METHODS: The effects of nintedanib on pulmonary arterial hypertension with proliferation of pulmonary vascular smooth muscle cells (PVSMCs) and luminal occlusion, on microvascular disease with apoptosis of microvascular endothelial cells (MVECs) and on fibroblast activation with myofibroblast differentiation and accumulation of extracellular matrix were analysed. We also studied the effects of nintedanib on the levels of key mediators involved in the pathogenesis of SSc and on macrophage polarisation. RESULTS: Nintedanib inhibited proliferation of PVSMCs and prevented thickening of the vessel walls and luminal occlusion of pulmonary arteries. Treatment with nintedanib also inhibited apoptosis of MVECs and blunted the capillary rarefaction in Fra2-transgenic mice. These effects were associated with a normalisation of the serum levels of vascular endothelial growth factor in Fra2 mice on treatment with nintedanib. Nintedanib also effectively blocked myofibroblast differentiation and reduced pulmonary, dermal and myocardial fibrosis in Fra2-transgenic mice. The antifibrotic effects of nintedanib were associated with impaired M2 polarisation of monocytes and reduced numbers of M2 macrophages. CONCLUSION: Nintedanib targets core features of SSc in Fra2-transgenic mice and ameliorates histological features of pulmonary arterial hypertension, destructive microangiopathy and pulmonary and dermal fibrosis. These data might have direct implications for the ongoing phase III clinical trial with nintedanib in SSc-associated interstitial lung disease.

gsSKAT: Rapid gene set analysis and multiple testing correction for rare-variant association studies using weighted linear kernels.

Next-generation sequencing technologies have afforded unprecedented characterization of low-frequency and rare genetic variation. Due to low power for single-variant testing, aggregative methods are commonly used to combine observed rare variation within a single gene. Causal variation may also aggregate across multiple genes within relevant biomolecular pathways. Kernel-machine regression and adaptive testing methods for aggregative rare-variant association testing have been demonstrated to be powerful approaches for pathway-level analysis, although these methods tend to be computationally intensive at high-variant dimensionality and require access to complete data. An additional analytical issue in scans of large pathway definition sets is multiple testing correction. Gene set definitions may exhibit substantial genic overlap, and the impact of the resultant correlation in test statistics on Type I error rate control for large agnostic gene set scans has not been fully explored. Herein, we first outline a statistical strategy for aggregative rare-variant analysis using component gene-level linear kernel score test summary statistics as well as derive simple estimators of the effective number of tests for family-wise error rate control. We then conduct extensive simulation studies to characterize the behavior of our approach relative to direct application of kernel and adaptive methods under a variety of conditions. We also apply our method to two case-control studies, respectively, evaluating rare variation in hereditary prostate cancer and schizophrenia. Finally, we provide open-source R code for public use to facilitate easy application of our methods to existing rare-variant analysis results.

Inhibition of phosphodiesterase 4 (PDE4) reduces dermal fibrosis by interfering with the release of interleukin-6 from M2 macrophages.

OBJECTIVES: To investigate the disease-modifying effects of phosphodiesterase 4 (PDE4) inhibition in preclinical models of systemic sclerosis (SSc). METHODS: We studied the effects of PDE4 inhibition in a prevention and a treatment model of bleomycin-induced skin fibrosis, in the topoisomerase mouse model as well as in a model of sclerodermatous chronic graft-versus-host disease. To better understand the mode of action of PDE4 blockade in preclinical models of SSc, we investigated fibrosis-relevant mediators in fibroblasts and macrophages from healthy individuals and patients suffering from diffuse-cutaneous SSc on blockade of PDE4. RESULTS: Specific inhibition of PDE4 by rolipram and apremilast had potent antifibrotic effects in bleomycin-induced skin fibrosis models, in the topoisomerase I mouse model and in murine sclerodermatous chronic graft-versus-host disease. Fibroblasts were not the direct targets of the antifibrotic effects of PDE4 blockade. Reduced leucocyte infiltration in lesional skin on PDE4 blockade suggested an immune-mediated mechanism. Further analysis revealed that PDE4 inhibition decreased the differentiation of M2 macrophages and the release of several profibrotic cytokines, resulting in reduced fibroblast activation and collagen release. Within these profibrotic mediators, interleukin-6 appeared to play a central role. CONCLUSIONS: PDE4 inhibition reduces inflammatory cell activity and the release of profibrotic cytokines from M2 macrophages, leading to decreased fibroblast activation and collagen release. Importantly, apremilast is already approved for the treatment of psoriasis and psoriatic arthritis. Therefore, PDE4 inhibitors might be further developed as potential antifibrotic therapies for patients with SSc. Our findings suggest that particularly patients with inflammation-driven fibrosis might benefit from PDE4 blockade.