Search

Extendend analysis of PARP7 expression in prostate cancer

FACS analysis of NCI-H660 cells for detection of necrosis and apoptosis

Supplementary Table 1: RBN2397 growth effects in prostate cencer cell lines

RBN2397 inhibition of prostate cancer cell growth requires PARP7 over-expression

Cell cycle distributions of cell lines treated with RBN2397 plus androgen or plus AHR agonist

BBQ and FICZ treatment activates AHR signaling in prostate cancer cells

The ADP-ribosyltransferase PARP7 modulates protein function by conjugating ADP-ribose to the side chains of acceptor amino acids. PARP7 has been shown to affect gene expression in prostate cancer cells and certain other cell types by mechanisms that include transcription factor ADP-ribosylation. Here, we use a recently developed catalytic inhibitor to PARP7, RBN2397, to study the effects of PARP7 inhibition in androgen receptor (AR)-positive and AR-negative prostate cancer cells. We find that RBN2397 has nanomolar potency for inhibiting androgen-induced ADP-ribosylation of the AR. RBN2397 inhibits the growth of prostate cancer cells in culture when cells are treated with ligands that activate the AR, or the aryl hydrocarbon receptor, and induce PARP7 expression. We show that the growth-inhibitory effects of RBN2397 are distinct from its enhancement of IFN signaling recently shown to promote tumor immunogenicity. RBN2397 treatment also induces trapping of PARP7 in a detergent-resistant fraction within the nucleus, which is reminiscent of how inhibitors such as talazoparib affect PARP1 compartmentalization. Because PARP7 is expressed in AR-negative metastatic tumors and RBN2397 can affect cancer cells through multiple mechanisms, PARP7 may be an actionable target in advanced prostate cancer.Significance:RBN2397 is a potent and selective inhibitor of PARP7 that reduces the growth of prostate cancer cells, including a model for treatment-emergent neuroendocrine prostate cancer. RBN2397 induces PARP7 trapping on chromatin, suggesting its mechanism of action might be similar to clinically used PARP1 inhibitors.

S1. A boxplot representing the percentage of variants called in the whole exome sequenced (WXS) tumor sample that is likely somatic mutations. S2. Correlation between the variant allele frequencies calculated from the four variant sets and the distribution of allele frequencies. S3. Principal components calculated from germline variants from the whole exome sequencing data from the non-tumor samples. S4. Kaplan-Meier plot for the germline variant rs28672782 in BRSK2. Table S1. Quality control checks reveal that somatic mutations and RNA editing did not affect the results of our analysis. Table S2. Correlation between the allele frequencies calculated in our four variant sets and the allele frequencies reported by gnomAD. Table S3. Variants genetically linked to rs61757955 in the European population, the population that is most similar to the TCGA lower grade glioma patient population. Table S4. Results from testing for an association between the germline variant rs34988193 and genomic and histological variables.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hematologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Background: Plasma proteins are critical mediators of cardiovascular processes and are the targets of many drugs. Previous efforts to characterize the genetic architecture of the plasma proteome have been limited by a focus on individuals of European descent and leveraged genotyping arrays and imputation. Here we describe whole genome sequence analysis of the plasma proteome in individuals with greater African ancestry, increasing our power to identify novel genetic determinants. Methods: Proteomic profiling of 1301 proteins was performed in 1852 Black adults from the Jackson Heart Study using aptamer-based proteomics (SomaScan). Whole genome sequencing association analysis was ascertained for all variants with minor allele count ≥5. Results were validated using an alternative, antibody-based, proteomic platform (Olink) as well as replicated in the Multi-Ethnic Study of Atherosclerosis and the HERITAGE Family Study (Health, Risk Factors, Exercise Training and Genetics). Results: We identify 569 genetic associations between 479 proteins and 438 unique genetic regions at a Bonferroni-adjusted significance level of 3.8×10 -11 . These associations include 114 novel locus-protein relationships and an additional 217 novel sentinel variant-protein relationships. Novel cardiovascular findings include new protein associations at the APOE gene locus including ZAP70 (sentinel single nucleotide polymorphism [SNP] rs7412-T, β=0.61±0.05, P =3.27×10 -30 ) and MMP-3 (β=-0.60±0.05, P =1.67×10 -32 ), as well as a completely novel pleiotropic locus at the HPX gene, associated with 9 proteins. Further, the associations suggest new mechanisms of genetically mediated cardiovascular disease linked to African ancestry; we identify a novel association between variants linked to APOL1-associated chronic kidney and heart disease and the protein CKAP2 (rs73885319-G, β=0.34±0.04, P =1.34×10 -17 ) as well as an association between ATTR amyloidosis and RBP4 levels in community-dwelling individuals without heart failure. Conclusions: Taken together, these results provide evidence for the functional importance of variants in non-European populations, and suggest new biological mechanisms for ancestry-specific determinants of lipids, coagulation, and myocardial function.

The ADP-ribosyltransferase PARP7 modulates protein function by conjugating ADP-ribose to the side chains of acceptor amino acids. PARP7 has been shown to affect gene expression in prostate cancer cells and certain other cell types by mechanisms that include transcription factor ADP-ribosylation. Here, we use a recently developed catalytic inhibitor to PARP7, RBN2397, to study the effects of PARP7 inhibition in androgen receptor-positive and androgen receptor-negative prostate cancer cells. We find that RBN2397 has nanomolar potency for inhibiting androgen-induced ADP-ribosylation of the androgen receptor. RBN2397 inhibits the growth of prostate cancer cells in culture when cells are treated with ligands that activate the androgen receptor, or the aryl hydrocarbon receptor, and induce PARP7 expression. We show that the growth inhibitory effects of RBN2397 are distinct from its enhancement of interferon signaling recently shown to promote tumor immunogenicity. RBN2397 treatment also induces trapping of PARP7 in a detergentresistant fraction within the nucleus, which is reminiscent of how inhibitors such as Talazoparib affect PARP1 fractionation. Because PARP7 is expressed in AR negative metastatic tumors and RBN2397 can affect cancer cells through multiple mechanisms, PARP7 may be an actionable target in advanced prostate cancer.SignificanceRBN2397 is a potent and selective inhibitor of PARP7 that reduces the growth of prostate cancer cells, including a model for treatment-emergent neuroendocrine prostate cancer. RBN2397 induces PARP7 trapping on chromatin, suggesting its mechanism of action might be similar to clinically-used PARP1 inhibitors.

Species across the tree of life can switch between asexual and sexual reproduction. In facultatively sexual species, the ability to switch between reproductive modes is often environmentally dependent and subject to local adaptation. However, the ecological and evolutionary factors that influence the maintenance and turnover of polymorphism associated with facultative sex remain unclear. We studied the ecological and evolutionary dynamics of reproductive investment in the facultatively sexual model species, Daphnia pulex. We found that patterns of clonal diversity, but not genetic diversity varied among ponds consistent with the predicted relationship between ephemerality and clonal structure. Reconstruction of a multi-year pedigree demonstrated the coexistence of clones that differ in their investment into male production. Mapping of quantitative variation in male production using lab-generated and field-collected individuals identified multiple putative quantitative trait loci (QTL) underlying this trait, and we identified a plausible candidate gene. The evolutionary history of these QTL suggests that they are relatively young, and male limitation in this system is a rapidly evolving trait. Our work highlights the dynamic nature of the genetic structure and composition of facultative sex across space and time and suggests that quantitative genetic variation in reproductive strategy can undergo rapid evolutionary turnover.

Chronic natural killer large granular lymphocyte (NK-LGL) leukemia, also referred to as chronic lymphoproliferative disorder of NK cells, is a rare disorder defined by prolonged expansion of clonal NK cells. Similar prevalence of STAT3 mutations in chronic T-LGL and NK-LGL leukemia is suggestive of common pathogenesis. We undertook whole-genome sequencing to identify mutations unique to NK-LGL leukemia. The results were analyzed to develop a resequencing panel that was applied to 58 patients. Phosphatidylinositol 3-kinase pathway gene mutations (PIK3CD/PIK3AP1) and TNFAIP3 mutations were seen in 5% and 10% of patients, respectively. TET2 was exceptional in that mutations were present in 16 (28%) of 58 patient samples, with evidence that TET2 mutations can be dominant and exclusive to the NK compartment. Reduced-representation bisulfite sequencing revealed that methylation patterns were significantly altered in TET2 mutant samples. The promoter of TET2 and that of PTPRD, a negative regulator of STAT3, were found to be methylated in additional cohort samples, largely confined to the TET2 mutant group. Mutations in STAT3 were observed in 19 (33%) of 58 patient samples, 7 of which had concurrent TET2 mutations. Thrombocytopenia and resistance to immunosuppressive agents were uniquely observed in those patients with only TET2 mutation (Games-Howell post hoc test, P = .0074; Fisher’s exact test, P = .00466). Patients with STAT3 mutation, inclusive of those with TET2 comutation, had lower hematocrit, hemoglobin, and absolute neutrophil count compared with STAT3 wild-type patients (Welch’s t test, P ≤ .015). We present the discovery of TET2 mutations in chronic NK-LGL leukemia and evidence that it identifies a unique molecular subtype.

Introduction: In this study, we investigated the prognostic role of the immune cell state atlas in predicting therapeutic benefits of patients treated with immune checkpoint inhibitors (ICI) within the ORIEN network of 18 collaborating cancer centers under the Total Cancer Care protocol. Methods: We utilized RNA-seq data of 926 samples generated from 875 individuals. Gene expression data were deconvoluted for immune cell states using the Carcinoma EcoTyper software. We then conducted a series of survival analyses to test the association between survival outcomes and predicted cell types and states in five malignant tumors: Genitourinary (GU), Gastrointestinal (GI), Thoracic (THO), Cutaneous (CUT), Head & Neck (H&N). The regularized Cox regression model in R package ‘glmnet’ was then applied to select the complementary pathway signatures (including gene ontology and KEGG pathways) to the immune cell states in predicting survival outcomes. We also explored the immune-related long non-coding RNAs (lncRNA) as potential biomarkers for cell states and patient outcomes. Results: EcoTyper analysis revealed that 692 (~80%) of patients were assigned to the 10 pre-identified Carcinoma Ecotypes (CE1 to CE10) or cell state atlas group. Overall, two immune deficiency ecotype patient groups (CE1 and CE2) pre-identified based on the independent training data were linked to worse survival, while two proinflammatory ecotype groups (CE9 and CE10) were associated with favorable surxvival. Those ecotype groups showed strong prognostic significance in predicting OS in melanoma and H&N. Meanwhile, CE6, a non-neoplastic tissue enriched cell subtype, was also found to be highly associated with longer OS in H&N and GU. CE7, an age-related mutation patient subgroup, contributed to shorter survival in both melanoma and GI. We also found that a subset of activated B cell state and the exhausted/effector CD4 T cell state were significantly associated with patient survival in melanoma and GU, respectively. The penalized Cox regression model revealed that β-catenin signaling pathway, P53 pathway and heme metabolism in the MSigDB Hallmark gene sets are the most complementary pathways to the ecotype scores in multiple cancer types. In additional, multiple pathways in KEGG such as endocytosis were found to jointly contribute to the ecotype-pathway composite prognostic model. In anazlying immune-related lncRNA biomarkers, we highlighted the prognostic role of NKILA in our dataset, which has been studied to promote tumor immune evasion. Conclusion: Our analysis has successfully established the utility of immune cell state atlas in predicting therapeutic benefits with ICIs. We expect that the discovered complementary signatures in the cancer-cell compartment will also lead to a novel spectrum of tumor-based biomarkers to ICI. Citation Format: Tingyi Li, Vineeth Sukrithan, Aakrosh Ratan, Martin McCarter, John Carpten, Howard Colman, Alexandra P. Ikeguchi, Igor Puzanov, Susanne Arnold, Michelle Churchman, Patrick Hwu, Paulo C. Rodriguez, William S. Dalton, George J. Weiner, Ahmad Tarhini, Xuefeng Wang. The immune cell state atlas analysis predicts therapeutic benefits with immune checkpoint inhibitors. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5703.

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the ‘venom-ome’ and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 ‘venom-ome-specific toxins’ (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery., Analysis of a near-chromosomal genome assembly and transcriptome profiling of the Indian cobra identifies genes expressed in the venom glands. These data should help develop a new antivenom.

BACKGROUND: Asthma is a common respiratory disorder with a highly heterogeneous nature that remains poorly understood. The objective was to use whole genome sequencing (WGS) data to identify regions of common genetic variation contributing to lung function in individuals with a diagnosis of asthma. METHODS: WGS data were generated for 1,053 individuals from trios and extended pedigrees participating in the family-based Genetic Epidemiology of Asthma in Costa Rica study. Asthma affection status was defined through a physician’s diagnosis of asthma, and most participants with asthma also had airway hyperresponsiveness (AHR) to methacholine. Family-based association tests for single variants were performed to assess the associations with lung function phenotypes. RESULTS: A genome-wide significant association was identified between baseline FEV(1)/FVC ratio and a single-nucleotide polymorphism in the top hit cysteine-rich secretory protein LCCL domain-containing 2 (CRISPLD2) (rs12051168; P = 3.6 × 10(−8) in the unadjusted model) that retained suggestive significance in the covariate-adjusted model (P = 5.6 × 10(−6)). Rs12051168 was also nominally associated with other related phenotypes: baseline FEV(1) (P = 3.3 × 10(−3)), postbronchodilator (PB) FEV(1) (7.3 × 10(−3)), and PB FEV(1)/FVC ratio (P = 2.7 × 10(−3)). The identified baseline FEV(1)/FVC ratio and rs12051168 association was meta-analyzed and replicated in three independent cohorts in which most participants with asthma also had confirmed AHR (combined weighted z-score P = .015) but not in cohorts without information about AHR. CONCLUSIONS: These findings suggest that using specific asthma characteristics, such as AHR, can help identify more genetically homogeneous asthma subgroups with genotype-phenotype associations that may not be observed in all children with asthma. CRISPLD2 also may be important for baseline lung function in individuals with asthma who also may have AHR.

The underrepresentation of non-Europeans in human genetic studies so far has limited the diversity of individuals in genomic datasets and led to reduced medical relevance for a large proportion of the world’s population. Population-specific reference genome datasets as well as genome-wide association studies in diverse populations are needed to address this issue. Here we describe the pilot phase of the GenomeAsia 100K Project. This includes a whole-genome sequencing reference dataset from 1,739 individuals of 219 population groups and 64 countries across Asia. We catalogue genetic variation, population structure, disease associations and founder effects. We also explore the use of this dataset in imputation, to facilitate genetic studies in populations across Asia and worldwide., Using whole-genome sequencing data from 1,739 individuals, the GenomeAsia 100K Project catalogues genetic variation, population structure and disease associations to facilitate genetic studies in Asian populations and increase representation in genetics studies worldwide.

BACKGROUND: The lymph node metastasis-derived LNCaP, the bone metastasis-derived PC3 (skull), and VCaP (vertebral) cell lines are widely used as preclinical models of human prostate cancer (CaP) and have been described in more than 19,000 publications. Here, we report on short-read whole-genome sequencing and genomic analyses of LNCaP, VCaP, and PC3 cells stably transduced with WT AR (PC3-AR). METHODS: LNCaP, VCaP, and PC3-AR cell lines were sequenced to an average depth of more than 30-fold using Illumina short-read sequencing. Using various computational methods, we identified and compared the single-nucleotide variants, copy-number profiles, and the structural variants observed in the three cell lines. RESULTS: LNCaP cells are composed of multiple subpopulations, which results in nonintegral copy number states and a high mutational load when the data is analyzed in bulk. All three cell lines contain pathogenic mutations and homozygous deletions in genes involved in DNA mismatch repair, along with deleterious mutations in cell-cycle, Wnt signaling, and other critical cellular processes. PC3-AR cells have a truncating mutation in TP53 and do not express the p53 protein. The VCaP cells contain a homozygous gain-of-function mutation in TP53 (p.R248W) that promotes cancer invasion, metastasis, and progression and has also been observed in prostate adenocarcinomas. In addition, we detect the signatures of chromothripsis of the q arms of chromosome 5 in both PC3-AR and VCaP cells, strengthening the association of TP53 inactivation with chromothripsis reported in other systems. CONCLUSIONS: Our work provides a resource for genetic, genomic, and biological studies employing these commonly-used prostate cancer cell lines.