CRISPR-Cas9 for selective targeting of somatic mutations in pancreatic cancers.

Somatic mutations are desirable targets for selective elimination of cancer, yet most are found within noncoding regions. We have adapted the CRISPR-Cas9 gene editing tool as a novel, cancer-specific killing strategy by targeting the subset of somatic mutations that create protospacer adjacent motifs (PAMs), which have evolutionally allowed bacterial cells to distinguish between self and non-self DNA for Cas9-induced double strand breaks. Whole genome sequencing (WGS) of paired tumor minus normal (T-N) samples from three pancreatic cancer patients (Panc480, Panc504, and Panc1002) showed an average of 417 somatic PAMs per tumor produced from single base substitutions. Further analyses of 591 paired T-N samples from The International Cancer Genome Consortium found medians of ∼455 somatic PAMs per tumor in pancreatic, ∼2800 in lung, and ∼3200 in esophageal cancer cohorts. Finally, we demonstrated 69-99% selective cell death of three targeted pancreatic cancer cell lines using 4-9 sgRNAs designed using the somatic PAM discovery approach. We also showed no off-target activity from these tumor-specific sgRNAs in either the patient's normal cells or an irrelevant cancer using WGS. This study demonstrates the potential of CRISPR-Cas9 as a novel and selective anti-cancer strategy, and supports the genetic targeting of adult cancers.

Islands of genomic stability in the face of genetically unstable metastatic cancer.

Introduction: Metastatic cancer affects millions of people worldwide annually and is the leading cause of cancer-related deaths. Most patients with metastatic disease are not eligible for surgical resection, and current therapeutic regimens have varying success rates, some with 5-year survival rates below 5%. Here we test the hypothesis that metastatic cancer can be genetically targeted by exploiting single base substitution mutations unique to individual cells that occur as part of normal aging prior to transformation. These mutations are targetable because ~10% of them form novel tumor-specific "NGG" protospacer adjacent motif (PAM) sites targetable by CRISPR-Cas9. Methods: Whole genome sequencing was performed on five rapid autopsy cases of patient-matched primary tumor, normal and metastatic tissue from pancreatic ductal adenocarcinoma decedents. CRISPR-Cas9 PAM targets were determined by bioinformatic tumor-normal subtraction for each patient and verified in metastatic samples by high-depth capture-based sequencing. Results: We found that 90% of PAM targets were maintained between primary carcinomas and metastases overall. We identified rules that predict PAM loss or retention, where PAMs located in heterozygous regions in the primary tumor can be lost in metastases (private LOH), but PAMs occurring in regions of loss of heterozygosity (LOH) in the primary tumor were universally conserved in metastases. Conclusions: Regions of truncal LOH are strongly retained in the presence of genetic instability, and therefore represent genetic vulnerabilities in pancreatic adenocarcinomas. A CRISPR-based gene therapy approach targeting these regions may be a novel way to genetically target metastatic cancer.

Genomic data resources of the Brain Somatic Mosaicism Network for neuropsychiatric diseases.

Somatic mosaicism is defined as an occurrence of two or more populations of cells having genomic sequences differing at given loci in an individual who is derived from a single zygote. It is a characteristic of multicellular organisms that plays a crucial role in normal development and disease. To study the nature and extent of somatic mosaicism in autism spectrum disorder, bipolar disorder, focal cortical dysplasia, schizophrenia, and Tourette syndrome, a multi-institutional consortium called the Brain Somatic Mosaicism Network (BSMN) was formed through the National Institute of Mental Health (NIMH). In addition to genomic data of affected and neurotypical brains, the BSMN also developed and validated a best practices somatic single nucleotide variant calling workflow through the analysis of reference brain tissue. These resources, which include >400 terabytes of data from 1087 subjects, are now available to the research community via the NIMH Data Archive (NDA) and are described here.

Somatic mutations reveal hyperactive Notch signaling and racial disparities in prurigo nodularis.

Prurigo nodularis (PN) is a chronic inflammatory skin disease that disproportionately affects African Americans and is characterized by pruritic skin nodules of unknown etiology. Little is known about genetic alterations in PN pathogenesis, especially relating to somatic events which are often implicated in inflammatory conditions. We thus performed whole-exome sequencing on 54 lesional and nonlesional skin biopsies from 17 PN patients and 10 atopic dermatitis (AD) patients for comparison. Somatic mutational analysis revealed that PN lesional skin harbors pervasive somatic mutations in fibrotic, neurotropic, and cancer-associated genes. Nonsynonymous mutations were most frequent in NOTCH1 and the Notch signaling pathway, a regulator of cellular proliferation and tissue fibrosis, and NOTCH1 mutations were absent in AD. Somatic copy-number analysis, combined with expression data, showed that recurrently deleted and downregulated genes in PN lesional skin are associated with axonal guidance and extension. Follow-up immunofluorescence validation demonstrated increased NOTCH1 expression in PN lesional skin fibroblasts and increased Notch signaling in PN lesional dermis. Finally, multi-center data revealed a significantly increased risk of NOTCH1-associated diseases in PN patients. In characterizing the somatic landscape of PN, we uncover novel insights into its pathophysiology and identify a role for dysregulated Notch signaling in PN.

Corrigendum: Convergent antibody responses are associated with broad neutralization of hepatitis C virus.

[This corrects the article DOI: 10.3389/fimmu.2023.1135841.].

CRISPR-Cas9 for selective targeting of somatic mutations in pancreatic cancers.

Somatic mutations are desirable targets for selective elimination of cancer, yet most are found within the noncoding regions. We propose a novel, cancer-specific killing approach using CRISPR-Cas9 which exploits the requirement of a protospacer adjacent motif (PAM) for Cas9 activity. Through whole genome sequencing (WGS) of paired tumor minus normal (T-N) samples from three pancreatic cancer patients (Panc480, Panc504, and Panc1002), we identified an average of 417 somatic PAMs per tumor produced from single base substitutions. We analyzed 591 paired T-N samples from The International Cancer Genome Consortium and discovered medians of ~455 somatic PAMs per tumor in pancreatic, ~2800 in lung, and ~3200 in esophageal cancer cohorts. Finally, we demonstrated >80% selective cell death of two targeted pancreatic cancer cell lines in co-cultures using 4-9 sgRNAs, targeting noncoding regions, designed from the somatic PAM discovery approach. We also showed no off-target activity from these tumor-specific sgRNAs through WGS.

Convergent antibody responses are associated with broad neutralization of hepatitis C virus.

Introduction: Early development of broadly neutralizing antibodies (bNAbs) targeting the hepatitis C virus (HCV) envelope glycoprotein E2 is associated with spontaneous clearance of infection, so induction of bNAbs is a major goal of HCV vaccine development. However, the molecular antibody features important for broad neutralization are not known. Methods: To identify B cell repertoire features associated with broad neutralization, we performed RNA sequencing of the B cell receptors (BCRs) of HCV E2-reactive B cells of HCV-infected individuals with either high or low plasma neutralizing breadth. We then produced a monoclonal antibody (mAb) expressed by pairing the most abundant heavy and light chains from public clonotypes identified among clearance, high neutralization subjects. Results: We found distinctive BCR features associated with broad neutralization of HCV, including long heavy chain complementarity determining region 3 (CDRH3) regions, specific VH gene usage, increased frequencies of somatic hypermutation, and particular VH gene mutations. Most intriguing, we identified many E2-reactive public BCR clonotypes (heavy and light chain clones with the same V and J-genes and identical CDR3 sequences) present only in subjects who produced highly neutralizing plasma. The majority of these public clonotypes were shared by two subjects who cleared infection. A mAb expressing the most abundant public heavy and light chains from these clearance, high neutralization subjects had features enriched in high neutralization clonotypes, such as increased somatic hypermutation frequency and usage of IGHV1-69, and was cross-neutralizing. Discussion: Together, these results demonstrate distinct BCR repertoires associated with high plasma neutralizing capacity. Further characterization of the molecular features and function of these antibodies can inform HCV vaccine development.

Single-cell atlas of epithelial and stromal cell heterogeneity by lobe and strain in the mouse prostate.

BACKGROUND: Evaluating the complex interplay of cell types in the tissue microenvironment is critical to understanding the origin and progression of diseases in the prostate and potential opportunities for intervention. Mouse models are an essential tool to investigate the molecular and cell-type-specific contributions of prostate disease at an organismal level. While there are well-documented differences in the extent, timing, and nature of disease development in various genetically engineered and exposure-based mouse models in different mouse strains and prostate lobes within each mouse strain, the underlying molecular phenotypic differences in cell types across mouse strains and prostate lobes are incompletely understood. METHODS: In this study, we used single-cell RNA-sequencing (scRNA-seq) methods to assess the single-cell transcriptomes of 6-month-old mouse prostates from two commonly used mouse strains, friend virus B/NIH jackson (FVB/NJ) (N = 2) and C57BL/6J (N = 3). For each mouse, the lobes of the prostate were dissected (anterior, dorsal, lateral, and ventral), and individual scRNA-seq libraries were generated. In situ and pathological analyses were used to explore the spatial and anatomical distributions of novel cell types and molecular markers defining these cell types. RESULTS: Data dimensionality reduction and clustering analysis of scRNA-seq data revealed that basal and luminal cells possessed strain-specific transcriptomic differences, with luminal cells also displaying marked lobe-specific differences. Gene set enrichment analysis comparing luminal cells by strain showed enrichment of proto-Oncogene targets in FVB/NJ mice. Additionally, three rare populations of epithelial cells clustered independently of strain and lobe: one population of luminal cells expressing Foxi1 and components of the vacuolar ATPase proton pump (Atp6v0d2 and Atp6v1g3), another population expressing Psca and other stem cell-associated genes (Ly6a/Sca-1, Tacstd2/Trop-2), and a neuroendocrine population expressing Chga, Chgb, and Syp. In contrast, stromal cell clusters, including fibroblasts, smooth muscle cells, endothelial cells, pericytes, and immune cell types, were conserved across strain and lobe, clustering largely by cell type and not by strain or lobe. One notable exception to this was the identification of two distinct fibroblast populations that we term subglandular fibroblasts and interstitial fibroblasts based on their strikingly distinct spatial distribution in the mouse prostate. CONCLUSIONS: Altogether, these data provide a practical reference of the transcriptional profiles of mouse prostate from two commonly used mouse strains and across all four prostate lobes.

Identifying Phased Mutations and Complex Rearrangements in Human Prostate Cancer Cell Lines through Linked-Read Whole-Genome Sequencing.

A limited number of cell lines have fueled the majority of preclinical prostate cancer research, but their genomes remain incompletely characterized. Here, we utilized whole-genome linked-read sequencing for comprehensive characterization of phased mutations and rearrangements in the most commonly used cell lines in prostate cancer research including PC3, LNCaP, DU145, CWR22Rv1, VCaP, LAPC4, MDA-PCa-2b, RWPE-1, and four derivative castrate-resistant (CR) cell lines LNCaP_Abl, LNCaP_C42b, VCaP-CR, and LAPC4-CR. Phasing of mutations allowed determination of "gene-level haplotype" to assess whether genes harbored heterozygous mutations in one or both alleles. Phased structural variant analysis allowed identification of complex rearrangement chains consistent with chromothripsis and chromoplexy. In addition, comparison of parental and derivative CR lines revealed previously known and novel genomic alterations associated with the CR phenotype. IMPLICATIONS: This study therefore comprehensively characterized phased genomic alterations in the commonly used prostate cancer cell lines, providing a useful resource for future prostate cancer research.

Entinostat Decreases Immune Suppression to Promote Antitumor Responses in a HER2+ Breast Tumor Microenvironment.

Therapeutic combinations to alter immunosuppressive, solid tumor microenvironments (TME), such as in breast cancer, are essential to improve responses to immune checkpoint inhibitors (ICI). Entinostat, an oral histone deacetylase inhibitor, has been shown to improve responses to ICIs in various tumor models with immunosuppressive TMEs. The precise and comprehensive alterations to the TME induced by entinostat remain unknown. Here, we employed single-cell RNA sequencing on HER2-overexpressing breast tumors from mice treated with entinostat and ICIs to fully characterize changes across multiple cell types within the TME. This analysis demonstrates that treatment with entinostat induced a shift from a protumor to an antitumor TME signature, characterized predominantly by changes in myeloid cells. We confirmed myeloid-derived suppressor cells (MDSC) within entinostat-treated tumors associated with a less suppressive granulocytic (G)-MDSC phenotype and exhibited altered suppressive signaling that involved the NFκB and STAT3 pathways. In addition to MDSCs, tumor-associated macrophages were epigenetically reprogrammed from a protumor M2-like phenotype toward an antitumor M1-like phenotype, which may be contributing to a more sensitized TME. Overall, our in-depth analysis suggests that entinostat-induced changes on multiple myeloid cell types reduce immunosuppression and increase antitumor responses, which, in turn, improve sensitivity to ICIs. Sensitization of the TME by entinostat could ultimately broaden the population of patients with breast cancer who could benefit from ICIs.

Comparative genomics of tadpole shrimps (Crustacea, Branchiopoda, Notostraca): Dynamic genome evolution against the backdrop of morphological stasis.

This analysis presents five genome assemblies of four Notostraca taxa. Notostraca origin dates to the Permian/Upper Devonian and the extant forms show a striking morphological similarity to fossil taxa. The comparison of sequenced genomes with other Branchiopoda genomes shows that, despite the morphological stasis, Notostraca share a dynamic genome evolution with high turnover for gene families' expansion/contraction and a transposable elements content comparable to other branchiopods. While Notostraca substitutions rate appears similar or lower in comparison to other branchiopods, a subset of genes shows a faster evolutionary pace, highlighting the difficulty of generalizing about genomic stasis versus dynamism. Moreover, we found that the variation of Triops cancriformis transposable elements content appeared linked to reproductive strategies, in line with theoretical expectations. Overall, besides providing new genomic resources for the study of these organisms, which appear relevant for their ecology and evolution, we also confirmed the decoupling of morphological and molecular evolution.

Decreased Activated CD4+ T Cell Repertoire Diversity After Antiretroviral Therapy in HIV-1/HCV Coinfection Correlates with CD4+ T Cell Recovery.

Dysfunctional immune activation accumulates during chronic viral infection and contributes to disease pathogenesis. In HIV-1, immune activation is exacerbated by concurrent infection with hepatitis C virus (HCV), accelerating depletion of CD4+ T cells. HIV-1 suppression with antiretroviral therapy (ART) generally reconstitutes CD4+ T cell counts, while also reducing the proportion that is activated. Whether this immune reconstitution also reduces the complexity of the CD4+ T cell population is unknown. We sought to characterize the relationship between activated CD4+ T cell repertoire diversity and immune reconstitution following ART in HIV-1/HCV coinfection. We extracted T cell receptor (TCR) sequences from RNA sequencing data obtained from activated CD4+ T cells of HIV-1/HCV coinfected individuals before and after treatment with ART (clinical trial NCT01285050). There was notable heterogeneity in both the extent of CD4+ T cell reconstitution and in the change in activated CD4+ TCR repertoire diversity following ART. Decreases in activated CD4+ TCR repertoire diversity following ART were predictive of the degree of CD4+ T cell reconstitution. The association of decreased activated CD4+ TCR repertoire diversity and improved CD4+ T cell reconstitution may represent loss of nonspecifically activated TCR clonotypes, and possibly selective expansion of specifically activated CD4+ clones. These results provide insight into the dynamic relationship between activated CD4+ TCR diversity and CD4+ T cell recovery of HIV-1/HCV coinfected individuals after suppression of HIV-1 viremia.

Transcriptional programs of neoantigen-specific TIL in anti-PD-1-treated lung cancers.

PD-1 blockade unleashes CD8 T cells1, including those specific for mutation-associated neoantigens (MANA), but factors in the tumour microenvironment can inhibit these T cell responses. Single-cell transcriptomics have revealed global T cell dysfunction programs in tumour-infiltrating lymphocytes (TIL). However, the majority of TIL do not recognize tumour antigens2, and little is known about transcriptional programs of MANA-specific TIL. Here, we identify MANA-specific T cell clones using the MANA functional expansion of specific T cells assay3 in neoadjuvant anti-PD-1-treated non-small cell lung cancers (NSCLC). We use their T cell receptors as a 'barcode' to track and analyse their transcriptional programs in the tumour microenvironment using coupled single-cell RNA sequencing and T cell receptor sequencing. We find both MANA- and virus-specific clones in TIL, regardless of response, and MANA-, influenza- and Epstein-Barr virus-specific TIL each have unique transcriptional programs. Despite exposure to cognate antigen, MANA-specific TIL express an incompletely activated cytolytic program. MANA-specific CD8 T cells have hallmark transcriptional programs of tissue-resident memory (TRM) cells, but low levels of interleukin-7 receptor (IL-7R) and are functionally less responsive to interleukin-7 (IL-7) compared with influenza-specific TRM cells. Compared with those from responding tumours, MANA-specific clones from non-responding tumours express T cell receptors with markedly lower ligand-dependent signalling, are largely confined to HOBIThigh TRM subsets, and coordinately upregulate checkpoints, killer inhibitory receptors and inhibitors of T cell activation. These findings provide important insights for overcoming resistance to PD-1 blockade.

Cell wall protein variation, break-induced replication, and subtelomere dynamics in Candida glabrata.

Candida glabrata is an opportunistic pathogen of humans, responsible for up to 30% of disseminated candidiasis. Adherence of C. glabrata to host cells is mediated by adhesin-like proteins (ALPs), about half of which are encoded in the subtelomeres. We performed a de novo assembly of two C. glabrata strains, BG2 and BG3993, using long single-molecule real-time (SMRT) reads, and constructed high-quality telomere-to-telomere assemblies of all 13 chromosomes to assess differences between C. glabrata strains. We documented variation between strains, and in agreement with earlier studies, found high (~0.5%-1%) frequencies of SNVs across the genome, including within subtelomeric regions. We documented changes in ALP gene structure and complement: there are large length differences in ALP genes in different strains, resulting from copy number variation in tandem repeats. We compared strains to characterize chromosome rearrangement events including within the poorly characterized subtelomeric regions. We show that rearrangements within the subtelomere regions all affect ALP-encoding genes, and 14/16 involve just the most terminal ALP gene. We present evidence that these rearrangements are mediated by break-induced replication. This study highlights the constrained nature of subtelomeric changes impacting ALP gene complement and subtelomere structure.

Hobotnica: exploring molecular signature quality.

A Molecular Features Set (MFS), is a result of a vast diversity of bioinformatics pipelines. The lack of a "gold standard" for most experimental data modalities makes it difficult to provide valid estimation for a particular MFS's quality. Yet, this goal can partially be achieved by analyzing inner-sample Distance Matrices (DM) and their power to distinguish between phenotypes. The quality of a DM can be assessed by summarizing its power to quantify the differences of inner-phenotype and outer-phenotype distances. This estimation of the DM quality can be construed as a measure of the MFS's quality.  Here we propose Hobotnica, an approach to estimate MFSs quality by their ability to stratify data, and assign them significance scores, that allow for collating various signatures and comparing their quality for contrasting groups.

Multi-ancestry fine mapping of interferon lambda and the outcome of acute hepatitis C virus infection.

Clearance of acute infection with hepatitis C virus (HCV) is associated with the chr19q13.13 region containing the rs368234815 (TT/ΔG) polymorphism. We fine-mapped this region to detect possible causal variants that may contribute to HCV clearance. First, we performed sequencing of IFNL1-IFNL4 region in 64 individuals sampled according to rs368234815 genotype: TT/clearance (N = 16) and ΔG/persistent (N = 15) (genotype-outcome concordant) or TT/persistent (N = 19) and ΔG/clearance (N = 14) (discordant). 25 SNPs had a difference in counts of alternative allele >5 between clearance and persistence individuals. Then, we evaluated those markers in an association analysis of HCV clearance conditioning on rs368234815 in two groups of European (692 clearance/1 025 persistence) and African ancestry (320 clearance/1 515 persistence) individuals. 10/25 variants were associated (P < 0.05) in the conditioned analysis leaded by rs4803221 (P value = 4.9 × 10-04) and rs8099917 (P value = 5.5 × 10-04). In the European ancestry group, individuals with the haplotype rs368234815ΔG/rs4803221C were 1.7× more likely to clear than those with the rs368234815ΔG/rs4803221G haplotype (P value = 3.6 × 10-05). For another nearby SNP, the haplotype of rs368234815ΔG/rs8099917T was associated with HCV clearance compared to rs368234815ΔG/rs8099917G (OR: 1.6, P value = 1.8 × 10-04). We identified four possible causal variants: rs368234815, rs12982533, rs10612351 and rs4803221. Our results suggest a main signal of association represented by rs368234815, with contributions from rs4803221, and/or nearby SNPs including rs8099917.

A medical student scholarly concentrations program: scholarly self-efficacy and impact on future research activities.

BACKGROUND: The Scholarly Concentrations program was established at Johns Hopkins University School of Medicine in 2009 with the aim of instilling passion for scholarship. OBJECTIVE: Our study aimed to determine whether the Scholarly Concentrations program achieves positive changes in medical student self-efficacy in conducting research and, if so, whether this results in future career aspirations toward scholarship. DESIGN: We used the Clinical Research Appraisal Inventory-Short Form (CRAI-SF) to assess changes in self-efficacy among students completing the Scholarly Concentrations program between 2014 and 2017. We calculated composite mean scores of six domains. We included outcomes on whether students published a manuscript, overall program perceptions, and likelihood of future research careers. We analyzed relationships between CRAI-SF scores and outcomes using paired t-tests and multivariable-adjusted logistic regression. RESULTS: A total of 419 students completed the Scholarly Concentrations program. All 6 CRAI domain scores showed significant improvements in self-efficacy between the pre-Scholarly Concentrations and post-Scholarly Concentrations ratings (range of changes 0.76-1.39, p < 0.05 for all). We found significant associations between post-Scholarly Concentrations self-efficacy ratings and course satisfaction (adjusted OR 1.57 [95% CI 1.20, 2.07]) and mentor satisfaction (OR 1.46 [1.15, 1.86]), as well as students' intent to conduct future research (OR 1.46 [1.15, 1.86]). These results were robust to sensitivity analyses, and pronounced in the group of students without prior research experience. CONCLUSIONS: Our findings suggest that a Scholarly Concentrations program is associated with an increased self-efficacy for research, and these changes in self-efficacy are associated with higher satisfaction in the scholarly experience and increased likelihood of pursuing scholarly work. Other medical schools could use such a tool of self-efficacy to both investigate the overall Scholarly Concentrations experience and understand factors that may increase interest in future physician-scientist pathways.

Pervasive promoter hypermethylation of silenced TERT alleles in human cancers.

BACKGROUND: In cancers, maintenance of telomeres often occurs through activation of the catalytic subunit of telomerase, encoded by TERT. Yet, most cancers show only modest levels of TERT gene expression, even in the context of activating hotspot promoter mutations (C228T and C250T). The role of epigenetic mechanisms, including DNA methylation, in regulating TERT gene expression in cancer cells is as yet not fully understood. METHODS: Here, we have carried out the most comprehensive characterization to date of TERT promoter methylation using ultra-deep bisulfite sequencing spanning the CpG island surrounding the core TERT promoter in 96 different human cell lines, including primary, immortalized and cancer cell types, as well as in control and reference samples. RESULTS: In general, we observed that immortalized and cancer cell lines were hypermethylated in a region upstream of the recurrent C228T and C250T TERT promoter mutations, while non-malignant primary cells were comparatively hypomethylated in this region. However, at the allele-level, we generally found that hypermethylation of promoter sequences in cancer cells is associated with repressed expression, and the remaining unmethylated alleles marked with open chromatin are largely responsible for the observed TERT expression in cancer cells. CONCLUSIONS: Our findings suggest that hypermethylation of the TERT promoter alleles signals transcriptional repression of those alleles, leading to attenuation of TERT activation in cancer cells. This type of fine tuning of TERT expression may account for the modest activation of TERT expression in most cancers.