Lineage-specific canonical and non-canonical activity of EZH2 in advanced prostate cancer subtypes.

Enhancer of zeste homolog 2 (EZH2) is a histone methyltransferase and emerging therapeutic target that is overexpressed in most castration-resistant prostate cancers and implicated as a driver of disease progression and resistance to hormonal therapies. Here we define the lineage-specific action and differential activity of EZH2 in both prostate adenocarcinoma (PRAD) and neuroendocrine prostate cancer (NEPC) subtypes of advanced prostate cancer to better understand the role of EZH2 in modulating differentiation, lineage plasticity, and to identify mediators of response and resistance to EZH2 inhibitor therapy. Mechanistically, EZH2 modulates bivalent genes that results in upregulation of NEPC-associated transcriptional drivers (e.g., ASCL1) and neuronal gene programs, and leads to forward differentiation after targeting EZH2 in NEPC. Subtype-specific downstream effects of EZH2 inhibition on cell cycle genes support the potential rationale for co-targeting cyclin/CDK to overcome resistance to EZH2 inhibition.

Heterogeneity and transcriptional drivers of triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is a heterogeneous disease with limited treatment options. To characterize TNBC heterogeneity, we defined transcriptional, epigenetic, and metabolic subtypes and subtype-driving super-enhancers and transcription factors by combining functional and molecular profiling with computational analyses. Single-cell RNA sequencing revealed relative homogeneity of the major transcriptional subtypes (luminal, basal, and mesenchymal) within samples. We found that mesenchymal TNBCs share features with mesenchymal neuroblastoma and rhabdoid tumors and that the PRRX1 transcription factor is a key driver of these tumors. PRRX1 is sufficient for inducing mesenchymal features in basal but not in luminal TNBC cells via reprogramming super-enhancer landscapes, but it is not required for mesenchymal state maintenance or for cellular viability. Our comprehensive, large-scale, multiplatform, multiomics study of both experimental and clinical TNBC is an important resource for the scientific and clinical research communities and opens venues for future investigation.

Androgen receptor inhibition suppresses anti-tumor neutrophil response against bone metastatic prostate cancer via regulation of TßRI expression.

Bone metastatic disease of prostate cancer (PCa) is incurable and progression in bone is largely dictated by tumor-stromal interactions in the bone microenvironment. We showed previously that bone neutrophils initially inhibit bone metastatic PCa growth yet metastatic PCa becomes resistant to neutrophil response. Further, neutrophils isolated from tumor-bone lost their ability to suppress tumor growth through unknown mechanisms. With this study, our goal was to define the impact of metastatic PCa on neutrophil function throughout tumor progression and to determine the potential of neutrophils as predictive biomarkers of metastatic disease. Using patient peripheral blood polymorphonuclear neutrophils (PMNs), we identified that PCa progression dictates PMN cell surface markers and gene expression, but not cytotoxicity against PCa. Importantly, we also identified a novel phenomenon in which second generation androgen deprivation therapy (ADT) suppresses PMN cytotoxicity via increased transforming growth factor beta receptor I (TßRI). High dose testosterone and genetic or pharmacologic TßRI inhibition rescued androgen receptor-mediated neutrophil suppression and restored neutrophil anti-tumor immune response. These studies highlight the ability to leverage standard-care ADT to generate neutrophil anti-tumor responses against bone metastatic PCa.

Proteogenomic analysis of chemo-refractory high-grade serous ovarian cancer.

To improve the understanding of chemo-refractory high-grade serous ovarian cancers (HGSOCs), we characterized the proteogenomic landscape of 242 (refractory and sensitive) HGSOCs, representing one discovery and two validation cohorts across two biospecimen types (formalin-fixed paraffin-embedded and frozen). We identified a 64-protein signature that predicts with high specificity a subset of HGSOCs refractory to initial platinum-based therapy and is validated in two independent patient cohorts. We detected significant association between lack of Ch17 loss of heterozygosity (LOH) and chemo-refractoriness. Based on pathway protein expression, we identified 5 clusters of HGSOC, which validated across two independent patient cohorts and patient-derived xenograft (PDX) models. These clusters may represent different mechanisms of refractoriness and implicate putative therapeutic vulnerabilities.

Telomerase Upregulation Induces Progression of Mouse BrafV600E-Driven Thyroid Cancers and Triggers Nontelomeric Effects.

Mutations in the promoter of the telomerase reverse transcriptase (TERT) gene are the paradigm of a cross-cancer alteration in a noncoding region. TERT promoter mutations (TPM) are biomarkers of poor prognosis in cancer, including thyroid tumors. TPMs enhance TERT transcription, which is otherwise silenced in adult tissues, thus reactivating a bona fide oncoprotein. To study TERT deregulation and its downstream consequences, we generated a Tert mutant promoter mouse model via CRISPR/Cas9 engineering of the murine equivalent locus (Tert-123C>T) and crossed it with thyroid-specific BrafV600E-mutant mice. We also employed an alternative model of Tert overexpression (K5-Tert). Whereas all BrafV600E animals developed well-differentiated papillary thyroid tumors, 29% and 36% of BrafV600E+Tert-123C>T and BrafV600E+K5-Tert mice progressed to poorly differentiated cancers at week 20, respectively. Tert-upregulated tumors showed increased mitosis and necrosis in areas of solid growth, and older animals displayed anaplastic-like features, that is, spindle cells and macrophage infiltration. Murine TPM increased Tert transcription in vitro and in vivo, but temporal and intratumoral heterogeneity was observed. RNA-sequencing of thyroid tumor cells showed that processes other than the canonical Tert-mediated telomere maintenance role operate in these specimens. Pathway analysis showed that MAPK and PI3K/AKT signaling, as well as processes not previously associated with this tumor etiology, involving cytokine, and chemokine signaling, were overactivated. These models constitute useful preclinical tools to understand the cell-autonomous and microenvironment-related consequences of Tert-mediated progression in advanced thyroid cancers and other aggressive tumors carrying TPMs. IMPLICATIONS: Telomerase-driven cancer progression activates pathways that can be dissected and perhaps therapeutically exploited.

Uncovering transcriptomic landscape alterations of CAN-2409 in in vitro and in vivo glioma models.

Rationale: CAN-2409 is a locally delivered oncolytic therapy, which results in vaccination against the injected tumor. CAN-2409 consists of a non-replicating adenovirus armed with the Herpes virus thymidine kinase, which metabolizes ganciclovir into a phosphorylated nucleotide that is incorporated into the tumor cell's genome, thereby inflicting immunogenic cancer cell death. While CAN-2409's immunological impact has been well characterized, its effects on the tumor cells transcriptome remains unknown. We compared the transcriptomic landscape after treatment of glioblastoma models with CAN-2409 in vitro and in vivo to assess how the interplay with the tumor microenvironment influences CAN-2409-mediated transcriptome alterations. Methods: We performed RNA-Seq with CAN-2409 treated patient-derived glioma stem-like cells and tumors of C57/BL6 mice and compared KEGG pathway usage and differential gene expression focusing on immune cell and cytokine profiles. T-cell -killing assays were performed to assess candidate effectors. Results: PCA analysis showed distinct clustering of control and CAN-2409 samples under both conditions. KEGG pathway analysis revealed significant enrichment for p53 signaling and cell cycle pathway, with similar dynamics for key regulators of both pathways in vitro and in vivo, including MYC, CCNB1, PLK1 and CDC20. Selected alterations (PLK1 and CCNB1) were validated at the protein level. Cytokine expression analysis revealed upregulation of pro-inflammatory IL12a under both conditions; immune cell gene profiling showed reduction of myeloid associated genes. T-cell-killing assays showed increased killing in the presence of IL-12. Conclusion: CAN-2409 significantly alters the transcriptome both in vitro and in vivo. Comparison of pathway enrichment revealed mutual and differential utilization of pathways under both conditions, suggesting a modulating influence on the cell cycle in tumor cells, and of the tumor microenvironment on the transcriptome in vivo. IL-12 synthesis likely depends on interactions with the tumor microenvironment, and it facilitates CAN-2409 cell killing. This dataset provides potential to understand resistance mechanisms and identify potential biomarkers for future studies.

Comprehensive genetic and phenotype analysis of 95 individuals with mosaic tuberous sclerosis complex.

Tuberous sclerosis complex (TSC) is a neurogenetic disorder due to loss-of-function TSC1 or TSC2 variants, characterized by tumors affecting multiple organs, including skin, brain, heart, lung, and kidney. Mosaicism for TSC1 or TSC2 variants occurs in 10%-15% of individuals diagnosed with TSC. Here, we report comprehensive characterization of TSC mosaicism by using massively parallel sequencing (MPS) of 330 TSC samples from a variety of tissues and fluids from a cohort of 95 individuals with mosaic TSC. TSC1 variants in individuals with mosaic TSC are much less common (9%) than in germline TSC overall (26%) (p < 0.0001). The mosaic variant allele frequency (VAF) is significantly higher in TSC1 than in TSC2, in both blood and saliva (median VAF: TSC1, 4.91%; TSC2, 1.93%; p = 0.036) and facial angiofibromas (median VAF: TSC1, 7.7%; TSC2 3.7%; p = 0.004), while the number of TSC clinical features in individuals with TSC1 and TSC2 mosaicism was similar. The distribution of mosaic variants across TSC1 and TSC2 is similar to that for pathogenic germline variants in general TSC. The systemic mosaic variant was not present in blood in 14 of 76 (18%) individuals with TSC, highlighting the value of analysis of multiple samples from each individual. A detailed comparison revealed that nearly all TSC clinical features are less common in individuals with mosaic versus germline TSC. A large number of previously unreported TSC1 and TSC2 variants, including intronic and large rearrangements (n = 11), were also identified.

Telomerase reactivation induces progression of mouse Braf V600E -driven thyroid cancers without telomere lengthening.

Mutations in the promoter of the telomerase reverse transcriptase ( TERT ) gene are the paradigm of a cross-cancer alteration in a non-coding region. TERT promoter mutations (TPMs) are biomarkers of poor prognosis in several tumors, including thyroid cancers. TPMs enhance TERT transcription, which is otherwise silenced in adult tissues, thus reactivating a bona fide oncoprotein. To study TERT deregulation and its downstream consequences, we generated a Tert mutant promoter mouse model via CRISPR/Cas9 engineering of the murine equivalent locus (Tert -123C>T ) and crossed it with thyroid-specific Braf V600E -mutant mice. We also employed an alternative model of Tert overexpression (K5-Tert). Whereas all Braf V600E animals developed well-differentiated papillary thyroid tumors, 29% and 36% of Braf V600E +Tert -123C>T and Braf V600E +K5-Tert mice progressed to poorly differentiated thyroid cancers at week 20, respectively. Braf+Tert tumors showed increased mitosis and necrosis in areas of solid growth, and older animals from these cohorts displayed anaplastic-like features, i.e., spindle cells and macrophage infiltration. Murine Tert promoter mutation increased Tert transcription in vitro and in vivo , but temporal and intra-tumoral heterogeneity was observed. RNA-sequencing of thyroid tumor cells showed that processes other than the canonical Tert-mediated telomere maintenance role operate in these specimens. Pathway analysis showed that MAPK and PI3K/AKT signaling, as well as processes not previously associated with this tumor etiology, involving cytokine and chemokine signaling, were overactivated. Braf+Tert animals remained responsive to MAPK pathway inhibitors. These models constitute useful pre-clinical tools to understand the cell-autonomous and microenvironment-related consequences of Tert-mediated progression in advanced thyroid cancers and other aggressive tumors carrying TPMs.

Splicing factor SRSF1 controls autoimmune-related molecular pathways in regulatory T cells distinct from FoxP3.

Regulatory T cells (Tregs) are vital for maintaining immune self-tolerance, and their impaired function leads to autoimmune disease. Mutations in FoxP3, the master transcriptional regulator of Tregs, leads to immune dysregulation, polyendocrinopathy, enteropathy, X-linked (IPEX) syndrome in humans and the early lethal "scurfy" phenotype with multi-organ autoimmune disease in mice. We recently identified serine/arginine-rich splicing factor 1 (SRSF1) as an indispensable regulator of Treg homeostasis and function. Intriguingly, Treg-conditional SRSF1-deficient mice exhibit early lethal systemic autoimmunity with multi-organ inflammation reminiscent of the scurfy mice. Importantly, SRSF1 is decreased in T cells from patients with the autoimmune disease systemic lupus erythematosus (SLE), and low SRSF1 levels inversely correlate with disease severity. Given that the Treg-specific deficiency of SRSF1 causes similarly profound autoimmune disease outcomes in mice as the deficiency/mutation in FoxP3, we aimed to evaluate the genes and molecular pathways controlled by these two indispensable regulatory proteins. We performed comparative bioinformatic analyses of transcriptomic profiles of Tregs from Srsf1-knockout mice and two Foxp3 mutant mice--the FoxP3-deficient ΔFoxp3 and the Foxp3 M370I mutant mice. We identified 132 differentially expressed genes (DEGs) unique to Srsf1-ko Tregs, 503 DEGs unique to Foxp3 M370I Tregs, and 1367 DEGs unique to ΔFoxp3 Tregs. Gene set enrichment and pathway analysis of DEGs unique to Srsf1-ko Tregs indicate that SRSF1 controls cytokine and immune response pathways. Conversely, FoxP3 controls pathways involved in DNA replication and cell cycle. Besides the distinct gene signatures, we identified only 30 shared genes between all three Treg mutants, mostly contributing to cytokine and immune defense pathways. Prominent genes included the chemokines CXCR6 and CCL1 and the checkpoint inhibitors FASLG and PDCD1. Thus, we demonstrate that SRSF1 and FoxP3 control common and distinct molecular pathways implicated in autoimmunity. Our analyses suggest that SRSF1 controls crucial immune functions in Tregs contributing to immune tolerance, and perturbations in its levels lead to systemic autoimmunity via mechanisms that are largely distinct from FoxP3.

Splicing factor SRSF1 is essential for CD8 T cell function and host antigen-specific viral immunity.

Cytotoxic CD8 T cells are crucial for the host antigen-specific immune response to viral pathogens. Here we report the identification of an essential role for the serine/arginine-rich splicing factor (SRSF) 1 in CD8 T cell homeostasis and function. Specifically, SRSF1 is necessary for the maintenance of normal CD8 T lymphocyte numbers in the lymphoid compartment, and for the proliferative capacity and cytotoxic function of CD8 T cells. Furthermore, SRSF1 is required for antigen-specific IFN-γ cytokine responses in a viral infection challenge in mice. Transcriptomics analyses of Srsf1-deficient T cells reveal that SRSF1 controls proliferation, MAP kinase signaling and IFN signaling pathways. Mechanistically, SRSF1 controls the expression and activity of the Mnk2/p38-MAPK axis at the molecular level. Our findings reveal previously unrecognized roles for SRSF1 in the physiology and function of cytotoxic CD8 T lymphocytes and a potential molecular mechanism in viral immunopathogenesis.

Heme oxygenase-1 mitigates liver injury and fibrosis via modulation of LNX1/Notch1 pathway in myeloid cells.

Activation of resident macrophages (Mϕ) and hepatic stellate cells is a key event in chronic liver injury. Mice with heme oxygenase-1 (HO-1; Hmox1)-deficient Mϕ (LysM-Cre:Hmox1 flfl ) exhibit increased inflammation, periportal ductular reaction, and liver fibrosis following bile duct ligation (BDL)-induced liver injury and increased pericellular fibrosis in NASH model. RiboTag-based RNA-sequencing profiling of hepatic HO-1-deficient Mϕ revealed dysregulation of multiple genes involved in lipid and amino acid metabolism, regulation of oxidative stress, and extracellular matrix turnover. Among these genes, ligand of numb-protein X1 (LNX1) expression is strongly suppressed in HO-1-deficient Mϕ. Importantly, HO-1 and LNX1 were expressed by hepatic Mϕ in human biliary and nonbiliary end-stage cirrhosis. We found that Notch1 expression, a downstream target of LNX1, was increased in LysM-Cre:Hmox1 flfl mice. In HO-1-deficient Mϕ treated with heme, transient overexpression of LNX1 drives M2-like Mϕ polarization. In summary, we identified LNX1/Notch1 pathway as a downstream target of HO-1 in liver fibrosis.

Ultrasensitive profiling of UV-induced mutations identifies thousands of subclinical facial tumors in tuberous sclerosis complex.

BackgroundTuberous sclerosis complex (TSC) is a neurogenetic syndrome due to loss-of-function mutations in TSC2 or TSC1, characterized by tumors at multiple body sites, including facial angiofibroma (FAF). Here, an ultrasensitive assessment of the extent and range of UV-induced mutations in TSC facial skin was performed.MethodsA multiplex high-sensitivity PCR assay (MHPA) was developed, enabling mutation detection at extremely low (<0.1%) variant allele frequencies (VAFs).ResultsMHPA assays were developed for both TSC2 and TP53, and applied to 81 samples, including 66 skin biopsies. UV-induced second-hit mutation causing inactivation of TSC2 was pervasive in TSC facial skin with an average of 4.8 mutations per 2-mm biopsy at median VAF 0.08%, generating more than 150,000 incipient facial tumors (subclinical "micro-FAFs") in the average TSC subject. The MHPA analysis also led to the identification of a refined UV-related indel signature and a recurrent complex mutation pattern, consisting of both a single-nucleotide or dinucleotide variant and a 1- to 9-nucleotide deletion, in cis.ConclusionTSC facial skin can be viewed as harboring a patchwork of clonal fibroblast proliferations (micro-FAFs) with indolent growth, a small proportion of which develop into clinically observable FAF. Our observations also expand the spectrum of UV-related mutation signatures.FundingThis work was supported by the TSC Alliance; the Engles Family Fund for Research in TSC and LAM; and the NIH, National Heart, Lung, and Blood Institute (U01HL131022-04 and Intramural Research Program).

A distinct core regulatory module enforces oncogene expression in KMT2A-rearranged leukemia.

Acute myeloid leukemia with KMT2A (MLL) rearrangements is characterized by specific patterns of gene expression and enhancer architecture, implying unique core transcriptional regulatory circuitry. Here, we identified the transcription factors MEF2D and IRF8 as selective transcriptional dependencies of KMT2A-rearranged AML, where MEF2D displays partially redundant functions with its paralog, MEF2C. Rapid transcription factor degradation followed by measurements of genome-wide transcription rates and superresolution microscopy revealed that MEF2D and IRF8 form a distinct core regulatory module with a narrow direct transcriptional program that includes activation of the key oncogenes MYC, HOXA9, and BCL2. Our study illustrates a mechanism of context-specific transcriptional addiction whereby a specific AML subclass depends on a highly specialized core regulatory module to directly enforce expression of common leukemia oncogenes.

Astrocyte-neuron crosstalk through Hedgehog signaling mediates cortical synapse development.

Neuron-glia interactions play a critical role in the regulation of synapse formation and circuit assembly. Here we demonstrate that canonical Sonic hedgehog (Shh) pathway signaling in cortical astrocytes acts to coordinate layer-specific synaptic connectivity. We show that the Shh receptor Ptch1 is expressed by cortical astrocytes during development and that Shh signaling is necessary and sufficient to promote the expression of genes involved in regulating synaptic development and layer-enriched astrocyte molecular identity. Loss of Shh in layer V neurons reduces astrocyte complexity and coverage by astrocytic processes in tripartite synapses; conversely, cell-autonomous activation of Shh signaling in astrocytes promotes cortical excitatory synapse formation. Furthermore, Shh-dependent genes Lrig1 and Sparc distinctively contribute to astrocyte morphology and synapse formation. Together, these results suggest that Shh secreted from deep-layer cortical neurons acts to specialize the molecular and functional features of astrocytes during development to shape circuit assembly and function.

Splicing factor SRSF1 controls distinct molecular programs in regulatory and effector T cells implicated in systemic autoimmune disease.

Systemic autoimmune diseases are characterized by hyperactive effector T cells (Teffs), aberrant cytokines and chemokines, and dysfunctional regulatory T cells (Tregs). We previously uncovered new roles for serine/arginine-rich splicing factor 1 (SRSF1) in the control of genes involved in T cell signaling and cytokine production in human T cells. SRSF1 levels are decreased in T cells from patients with systemic lupus erythematosus (SLE), and low levels correlate with severe disease. Moreover, T cell-conditional Srsf1-deficient mice recapitulate the autoimmune phenotype, exhibiting CD4 T cell hyperactivity, dysfunctional Tregs, systemic autoimmunity, and tissue inflammation. However, the role of SRSF1 in controlling molecular programs in Teffs and Tregs and how these pathways are implicated in autoimmunity is not known. Here, by comparative bioinformatics analysis, we demonstrate that SRSF1 controls largely distinct gene programs in Tregs and Teffs in vivo. SRSF1 regulates 189 differentially expressed genes (DEGs) unique to Tregs, 582 DEGs unique to Teffs, and 29 DEGs shared between both. Shared genes included IL-17A, IL-17F, CSF1, CXCL10, and CXCR4, and were highly enriched for inflammatory response and cytokine-cytokine receptor interaction pathways. SRSF1 controls distinct pathways in Tregs, which include chemokine signaling and immune cell differentiation, compared with pathways in Teffs, which include cytokine production, T cell homeostasis, and activation. We identified putative mRNA binding targets of SRSF1 which include CSF1, CXCL10, and IL-17F. Finally, comparisons with transcriptomics profiles from lupus-prone MRL/lpr mice reveal that SRSF1 controls genes and pathways implicated in autoimmune disease. The target genes of SRSF1 and putative binding targets we discovered, have known roles in systemic autoimmunity. Our findings suggest that SRSF1 controls distinct molecular pathways in Tregs and Teffs and aberrant SRSF1 levels may contribute to their dysfunction and immunopathogenesis of systemic autoimmune disease.

Performance assessment of DNA sequencing platforms in the ABRF Next-Generation Sequencing Study.

Assessing the reproducibility, accuracy and utility of massively parallel DNA sequencing platforms remains an ongoing challenge. Here the Association of Biomolecular Resource Facilities (ABRF) Next-Generation Sequencing Study benchmarks the performance of a set of sequencing instruments (HiSeq/NovaSeq/paired-end 2 × 250-bp chemistry, Ion S5/Proton, PacBio circular consensus sequencing (CCS), Oxford Nanopore Technologies PromethION/MinION, BGISEQ-500/MGISEQ-2000 and GS111) on human and bacterial reference DNA samples. Among short-read instruments, HiSeq 4000 and X10 provided the most consistent, highest genome coverage, while BGI/MGISEQ provided the lowest sequencing error rates. The long-read instrument PacBio CCS had the highest reference-based mapping rate and lowest non-mapping rate. The two long-read platforms PacBio CCS and PromethION/MinION showed the best sequence mapping in repeat-rich areas and across homopolymers. NovaSeq 6000 using 2 × 250-bp read chemistry was the most robust instrument for capturing known insertion/deletion events. This study serves as a benchmark for current genomics technologies, as well as a resource to inform experimental design and next-generation sequencing variant calling.

Intra- and extra-cellular environments contribute to the fate of HIV-1 infection.

HIV-1 entry into host cells leads to one of the following three alternative fates: (1) HIV-1 elimination by restriction factors, (2) establishment of HIV-1 latency, or (3) active viral replication in target cells. Here, we report the development of an improved system for monitoring HIV-1 fate at single-cell and population levels and show the diverse applications of this system to study specific aspects of HIV-1 fate in different cell types and under different environments. An analysis of the transcriptome of infected, primary CD4+ T cells that support alternative fates of HIV-1 identifies differential gene expression signatures in these cells. Small molecules are able to selectively target cells that support viral replication with no significant effect on viral latency. In addition, HIV-1 fate varies in different tissues following infection of humanized mice in vivo. Altogether, these studies indicate that intra- and extra-cellular environments contribute to the fate of HIV-1 infection.

Multiomic analysis identifies CPT1A as a potential therapeutic target in platinum-refractory, high-grade serous ovarian cancer.

Resistance to platinum compounds is a major determinant of patient survival in high-grade serous ovarian cancer (HGSOC). To understand mechanisms of platinum resistance and identify potential therapeutic targets in resistant HGSOC, we generated a data resource composed of dynamic (±carboplatin) protein, post-translational modification, and RNA sequencing (RNA-seq) profiles from intra-patient cell line pairs derived from 3 HGSOC patients before and after acquiring platinum resistance. These profiles reveal extensive responses to carboplatin that differ between sensitive and resistant cells. Higher fatty acid oxidation (FAO) pathway expression is associated with platinum resistance, and both pharmacologic inhibition and CRISPR knockout of carnitine palmitoyltransferase 1A (CPT1A), which represents a rate limiting step of FAO, sensitize HGSOC cells to platinum. The results are further validated in patient-derived xenograft models, indicating that CPT1A is a candidate therapeutic target to overcome platinum resistance. All multiomic data can be queried via an intuitive gene-query user interface (https://sites.google.com/view/ptrc-cell-line).