Phenotypic, genomic, and transcriptomic heterogeneity in a pancreatic cancer cell line.

OBJECTIVES: To evaluate the suitability of the MIA PaCa-2 cell line for studying pancreatic cancer intratumor heterogeneity, we aim to further characterize the nature of MIA PaCa-2 cells' phenotypic, genomic, and transcriptomic heterogeneity. METHODS: MIA PaCa-2 single-cell clones were established through flow cytometry. For the phenotypic study, we quantified the cellular morphology, proliferation rate, migration potential, and drug sensitivity of the clones. The chromosome copy number and transcriptomic profiles were quantified using SNPa and RNA-seq, respectively. RESULTS: Four MIA PaCa-2 clones showed distinctive phenotypes, with differences in cellular morphology, proliferation rate, migration potential, and drug sensitivity. We also observed a degree of genomic variations between these clones in form of chromosome copy number alterations and single nucleotide variations, suggesting the genomic heterogeneity of the population, and the intrinsic genomic instability of MIA PaCa-2 cells. Lastly, transcriptomic analysis of the clones also revealed gene expression profile differences between the clones, including the uniquely regulated ITGAV, which dictates the morphology of MIA PaCa-2 clones. CONCLUSIONS: MIA PaCa-2 is comprised of cells with distinctive phenotypes, heterogeneous genomes, and differential transcriptomic profiles, suggesting its suitability as a model to study the underlying mechanisms behind pancreatic cancer heterogeneity.

NKX2-2 based nuclei sorting on frozen human archival pancreas enables the enrichment of islet endocrine populations for single-nucleus RNA sequencing.

BACKGROUND: Current approaches to profile the single-cell transcriptomics of human pancreatic endocrine cells almost exclusively rely on freshly isolated islets. However, human islets are limited in availability. Furthermore, the extensive processing steps during islet isolation and subsequent single cell dissolution might alter gene expressions. In this work, we report the development of a single-nucleus RNA sequencing (snRNA-seq) approach with targeted islet cell enrichment for endocrine-population focused transcriptomic profiling using frozen archival pancreatic tissues without islet isolation. RESULTS: We cross-compared five nuclei isolation protocols and selected the citric acid method as the best strategy to isolate nuclei with high RNA integrity and low cytoplasmic contamination from frozen archival human pancreata. We innovated fluorescence-activated nuclei sorting based on the positive signal of NKX2-2 antibody to enrich nuclei of the endocrine population from the entire nuclei pool of the pancreas. Our sample preparation procedure generated high-quality single-nucleus gene-expression libraries while preserving the endocrine population diversity. In comparison with single-cell RNA sequencing (scRNA-seq) library generated with live cells from freshly isolated human islets, the snRNA-seq library displayed comparable endocrine cellular composition and cell type signature gene expression. However, between these two types of libraries, differential enrichments of transcripts belonging to different functional classes could be observed. CONCLUSIONS: Our work fills a technological gap and helps to unleash frozen archival pancreatic tissues for molecular profiling targeting the endocrine population. This study opens doors to retrospective mappings of endocrine cell dynamics in pancreatic tissues of complex histopathology. We expect that our protocol is applicable to enrich nuclei for transcriptomics studies from various populations in different types of frozen archival tissues.

Integrating SNVs and CNAs on a phylogenetic tree from single-cell DNA sequencing data.

Single-cell DNA sequencing enables the construction of evolutionary trees that can reveal how tumors gain mutations and grow. Different whole-genome amplification procedures render genomic materials of different characteristics, often suitable for the detection of either single-nucleotide variation or copy number aberration, but not ideally for both. Consequently, this hinders the inference of a comprehensive phylogenetic tree and limits opportunities to investigate the interplay of SNVs and CNAs. Existing methods such as SCARLET and COMPASS require that the SNVs and CNAs are detected from the same sets of cells, which is technically challenging. Here we present a novel computational tool, SCsnvcna, that places SNVs on a tree inferred from CNA signals, whereas the sets of cells rendering the SNVs and CNAs are independent, offering a more practical solution in terms of the technical challenges. SCsnvcna is a Bayesian probabilistic model using both the genotype constraints on the tree and the cellular prevalence to search the optimal solution. Comprehensive simulations and comparison with seven state-of-the-art methods show that SCsnvcna is robust and accurate in a variety of circumstances. Particularly, SCsnvcna most frequently produces the lowest error rates, with ability to scale to a wide range of numerical values for leaf nodes in the tree, SNVs, and SNV cells. The application of SCsnvcna to two published colorectal cancer data sets shows highly consistent placement of SNV cells and SNVs with the original study while also supporting a refined placement of ATP7B, illustrating SCsnvcna's value in analyzing complex multitumor samples.

Differential modulation of cellular phenotype and drug sensitivity by extracellular matrix proteins in primary and metastatic pancreatic cancer cells.

Pancreatic ductal adenocarcinoma (PDAC) is reported to be the third highest cause of cancer-related deaths in the United States. PDAC is known for its high proportion of stroma, which accounts for 90% of the tumor mass. The stroma is made up of extracellular matrix (ECM) and nonmalignant cells such as inflammatory cells, cancer-associated fibroblasts, and lymphatic and blood vessels. Here, we decoupled the effects of the ECM on PDAC cell lines by culturing cells on surfaces coated with different ECM proteins. Our data show that the primary tumor-derived cell lines have different morphology depending on the ECM proteins on which they are cultured, while metastatic lesion-derived PDAC lines' morphology does not change with respect to the different ECM proteins. Similarly, ECM proteins modulate the proliferation rate and the gemcitabine sensitivity of the primary tumor PDAC cell lines, but not the metastatic PDAC lines. Lastly, transcriptomics analysis of the primary tumor PDAC cells cultured on different ECM proteins reveals the regulation of various pathways, such as cell cycle, cell-adhesion molecules, and focal adhesion, including the regulation of several integrin genes that are essential for ECM recognition.

Differences in cell shape, motility, and growth reflect chromosomal number variations that can be visualized with live-cell ChReporters.

Chromosome numbers often change dynamically in tumors and cultured cells, which complicates therapy as well as understanding genotype-mechanotype relationships. Here we use a live-cell "ChReporter" method to identify cells with a single chromosomal loss in efforts to better understand differences in cell shape, motility, and growth. We focus on a standard cancer line and first show clonal populations that retain the ChReporter exhibit large differences in cell and nuclear morphology as well as motility. Phenotype metrics follow simple rules, including migratory persistence scaling with speed, and cytoskeletal differences are evident from drug responses, imaging, and single-cell RNA sequencing. However, mechanotype-genotype relationships between fluorescent ChReporter-positive clones proved complex and motivated comparisons of clones that differ only in loss or retention of a Chromosome-5 ChReporter. When lost, fluorescence-null cells show low expression of Chromosome-5 genes, including a key tumor suppressor APC that regulates microtubules and proliferation. Colonies are compact, nuclei are rounded, and cells proliferate more, with drug results implicating APC, and patient survival data indicating an association in multiple tumor-types. Visual identification of genotype with ChReporters can thus help clarify mechanotype and mechano-evolution.

Nucleus Mechanosensing in Cardiomyocytes.

Cardiac muscle contraction is distinct from the contraction of other muscle types. The heart continuously undergoes contraction-relaxation cycles throughout an animal's lifespan. It must respond to constantly varying physical and energetic burdens over the short term on a beat-to-beat basis and relies on different mechanisms over the long term. Muscle contractility is based on actin and myosin interactions that are regulated by cytoplasmic calcium ions. Genetic variants of sarcomeric proteins can lead to the pathophysiological development of cardiac dysfunction. The sarcomere is physically connected to other cytoskeletal components. Actin filaments, microtubules and desmin proteins are responsible for these interactions. Therefore, mechanical as well as biochemical signals from sarcomeric contractions are transmitted to and sensed by other parts of the cardiomyocyte, particularly the nucleus which can respond to these stimuli. Proteins anchored to the nuclear envelope display a broad response which remodels the structure of the nucleus. In this review, we examine the central aspects of mechanotransduction in the cardiomyocyte where the transmission of mechanical signals to the nucleus can result in changes in gene expression and nucleus morphology. The correlation of nucleus sensing and dysfunction of sarcomeric proteins may assist the understanding of a wide range of functional responses in the progress of cardiomyopathic diseases.

Confinement plus myosin-II suppression maximizes heritable loss of chromosomes, as revealed by live-cell ChReporters.

The mechanical environment of a cell can have many effects, but whether it impacts the DNA sequence of a cell has remained unexamined. To investigate this, we developed a live-cell method to measure changes in chromosome numbers. We edited constitutive genes with GFP or RFP tags on single alleles and discovered that cells that lose Chromosome reporters (ChReporters) become non-fluorescent. We applied our new tools to confined mitosis and to inhibition of the putative tumor suppressor myosin-II. We quantified compression of mitotic chromatin in vivo and demonstrated that similar compression in vitro resulted in cell death, but also rare and heritable ChReptorter loss. Myosin-II suppression rescued lethal multipolar divisions and maximized ChReporter loss during three-dimensional (3D) compression and two-dimensional (2D) lateral confinement, but not in standard 2D culture. ChReporter loss was associated with chromosome mis-segregation, rather than just the number of divisions, and loss in vitro and in mice was selected against in subsequent 2D cultures. Inhibition of the spindle assembly checkpoint (SAC) caused ChReporter loss in 2D culture, as expected, but not during 3D compression, suggesting a SAC perturbation. Thus, ChReporters enable diverse studies of viable genetic changes, and show that confinement and myosin-II affect DNA sequence and mechano-evolution.

Cooperative phagocytosis of solid tumours by macrophages triggers durable anti-tumour responses.

In solid tumours, the abundance of macrophages is typically associated with a poor prognosis. However, macrophage clusters in tumour-cell nests have been associated with survival in some tumour types. Here, by using tumour organoids comprising macrophages and cancer cells opsonized via a monoclonal antibody, we show that highly ordered clusters of macrophages cooperatively phagocytose cancer cells to suppress tumour growth. In mice with poorly immunogenic tumours, the systemic delivery of macrophages with signal-regulatory protein alpha (SIRPα) genetically knocked out or else with blockade of the CD47-SIRPα macrophage checkpoint was combined with the monoclonal antibody and subsequently triggered the production of endogenous tumour-opsonizing immunoglobulin G, substantially increased the survival of the animals and helped confer durable protection from tumour re-challenge and metastasis. Maximizing phagocytic potency by increasing macrophage numbers, by tumour-cell opsonization and by disrupting the phagocytic checkpoint CD47-SIRPα may lead to durable anti-tumour responses in solid cancers.

Genomic heterogeneity in pancreatic cancer organoids and its stability with culture.

The establishment of patient-derived pancreatic cancer organoid culture in recent years creates an exciting opportunity for researchers to perform a wide range of in vitro studies on a model that closely recapitulates the tumor. One of the outstanding question in pancreatic cancer biology is the causes and consequences of genomic heterogeneity observed in the disease. However, to use pancreatic cancer organoids as a model to study genomic variations, we need to first understand the degree of genomic heterogeneity and its stability within organoids. Here, we used single-cell whole-genome sequencing to investigate the genomic heterogeneity of two independent pancreatic cancer organoid lines, as well as their genomic stability with extended culture. Clonal populations with similar copy number profiles were observed within the organoids, and the proportion of these clones was shifted with extended culture, suggesting the growth advantage of some clones. However, sub-clonal genomic heterogeneity was also observed within each clonal population, indicating the genomic instability of the pancreatic cancer cells themselves. Furthermore, our transcriptomic analysis also revealed a positive correlation between copy number alterations and gene expression regulation, suggesting the "gene dosage" effect of these copy number alterations that translates to gene expression regulation.

CD47-SIRPα Checkpoint Disruption in Metastases Requires Tumor-Targeting Antibody for Molecular and Engineered Macrophage Therapies.

The macrophage checkpoint interaction CD47-SIRPα is an emerging target for cancer therapy, but clinical trials of monoclonal anti-CD47 show efficacy only in liquid tumors when combined with tumor-opsonizing IgG. Here, in challenging metastatic solid tumors, CD47 deletion shows no effect on tumor growth unless combined with otherwise ineffective tumor-opsonization, and we likewise show wild-type metastases are suppressed by SIRPα-blocked macrophages plus tumor-opsonization. Lung tumor nodules of syngeneic B16F10 melanoma cells with CD47 deletion show opsonization drives macrophage phagocytosis of B16F10s, consistent with growth versus phagocytosis calculus for exponential suppression of cancer. Wild-type CD47 levels on metastases in lungs of immunocompetent mice and on human metastases in livers of immunodeficient mice show that systemic injection of antibody-engineered macrophages also suppresses growth. Such in vivo functionality can be modulated by particle pre-loading of the macrophages. Thus, even though CD47-SIRPα disruption and tumor-opsonizing IgG are separately ineffective against established metastatic solid tumors, their combination in molecular and cellular therapies prolongs survival.

Scaling concepts in 'omics: Nuclear lamin-B scales with tumor growth and often predicts poor prognosis, unlike fibrosis.

Physicochemical principles such as stoichiometry and fractal assembly can give rise to characteristic scaling between components that potentially include coexpressed transcripts. For key structural factors within the nucleus and extracellular matrix, we discover specific gene-gene scaling exponents across many of the 32 tumor types in The Cancer Genome Atlas, and we demonstrate utility in predicting patient survival as well as scaling-informed machine learning (SIML). All tumors with adjacent tissue data show cancer-elevated proliferation genes, with some genes scaling with the nuclear filament LMNB1, including the transcription factor FOXM1 that we show directly regulates LMNB1 SIML shows that such regulated cancers cluster together with longer overall survival than dysregulated cancers, but high LMNB1 and FOXM1 in half of regulated cancers surprisingly predict poor survival, including for liver cancer. COL1A1 is also studied because it too increases in tumors, and a pan-cancer set of fibrosis genes shows substoichiometric scaling with COL1A1 but predicts patient outcome only for liver cancer-unexpectedly being prosurvival. Single-cell RNA-seq data show nontrivial scaling consistent with power laws from bulk RNA and protein analyses, and SIML segregates synthetic from contractile cancer fibroblasts. Our scaling approach thus yields fundamentals-based power laws relatable to survival, gene function, and experiments.

Gene Signatures of NEUROGENIN3+ Endocrine Progenitor Cells in the Human Pancreas.

NEUROGENIN3+ (NEUROG3+) cells are considered to be pancreatic endocrine progenitors. Our current knowledge on the molecular program of NEUROG3+ cells in humans is largely extrapolated from studies in mice. We hypothesized that single-cell RNA-seq enables in-depth exploration of the rare NEUROG3+ cells directly in humans. We aligned four large single-cell RNA-seq datasets from postnatal human pancreas. Our integrated analysis revealed 10 NEUROG3+ epithelial cells from a total of 11,174 pancreatic cells. Noticeably, human NEUROG3+ cells clustered with mature pancreatic cells and epsilon cells displayed the highest frequency of NEUROG3 positivity. We confirmed the co-expression of NEUROG3 with endocrine markers and the high percentage of NEUROG3+ cells among epsilon cells at the protein level based on immunostaining on pancreatic tissue sections. We further identified unique genetic signatures of the NEUROG3+ cells. Regulatory network inference revealed novel transcription factors including Prospero homeobox protein 1 (PROX1) may act jointly with NEUROG3. As NEUROG3 plays a central role in endocrine differentiation, knowledge gained from our study will accelerate the development of beta cell regeneration therapies to treat diabetes.

Polyploid mitosis and depolyploidization promote chromosomal instability and tumor progression in a Notch-induced tumor model.

Ploidy variation is a cancer hallmark and is frequently associated with poor prognosis in high-grade cancers. Using a Drosophila solid-tumor model where oncogenic Notch drives tumorigenesis in a transition-zone microenvironment in the salivary gland imaginal ring, we find that the tumor-initiating cells normally undergo endoreplication to become polyploid. Upregulation of Notch signaling, however, induces these polyploid transition-zone cells to re-enter mitosis and undergo tumorigenesis. Growth and progression of the transition-zone tumor are fueled by a combination of polyploid mitosis, endoreplication, and depolyploidization. Both polyploid mitosis and depolyploidization are error prone, resulting in chromosomal copy-number variation and polyaneuploidy. Comparative RNA-seq and epistasis analysis reveal that the DNA-damage response genes, also active during meiosis, are upregulated in these tumors and are required for the ploidy-reduction division. Together, these findings suggest that polyploidy and associated cell-cycle variants are critical for increased tumor-cell heterogeneity and genome instability during cancer progression.

Nuclear mechanosensing: mechanism and consequences of a nuclear rupture.

The physical connections between the cytoskeletal system and the nucleus provide a route for the nucleus to sense the mechanical stress both inside and outside of the cell. Failure to withstand such stress leads to nuclear rupture, which is observed in human diseases. In this review, we will go through the recent findings and our current understandings of nuclear rupture. Starting with the triggers of nuclear rupture, including the aberrant nuclear lamina composition and the elevated actomyosin contractility. We will also discuss the role of ESCRT-III in nuclear rupture repair and the biological consequences of nuclear rupture, including the negative impacts on cellular compartmentalization, DNA damage, and cellular differentiation. Recent studies on nuclear rupture provide further insights into the direct mechanistic link between nuclear rupture and several pathological conditions. Such knowledge can guide us in developing potential therapeutic solutions for the patients.

Aggregation-induced integrated stress response rejuvenates culture-expanded human mesenchymal stem cells.

Protein homeostasis is critical for cellular function, as loss of homeostasis is attributed to aging and the accumulation of unwanted proteins. Human mesenchymal stem cells (MSCs) have shown promising therapeutic potential due to their impressive abilities to secrete inflammatory modulators, angiogenic, and regenerative cytokines. However, there exists the problem of human MSC expansion with compromised therapeutic quality. Duringin vitro expansion, human MSCs are plated on stiff plastics and undergo culture adaptation, which results in aberrant proliferation, shifts in metabolism, and decreased autophagic activity. It has previously been shown that three-dimensional (3D) aggregation can reverse some of these alterations by heightening autophagy and recovering the metabolic state back to a naïve phenotype. To further understand the proteostasis in human MSC culture, this study investigated the effects of 3D aggregation on the human MSC proteome to determine the specific pathways altered by aggregation. The 3D aggregates and 2D cultures of human MSCs derived from bone marrow (bMSC) and adipose tissue (ASC) were analyzed along with differentiated human dermal fibroblasts (FB). The proteomics analysis showed the elevated eukaryotic initiation factor 2 pathway and the upregulated activity of the integrated stress response (ISR) in 3D aggregates. Specific protein quantification further determined that bMSC and ASC responded to ISR, while FB did not. 3D aggregation significantly increased the ischemic survival of bMSCs and ASCs. Perturbation of ISR with small molecules salubrinal and GSK2606414 resulted in differential responses of bMSC, ASC, and FB. This study indicates that aggregation-based preconditioning culture holds the potential for improving the therapeutic efficacy of expanded human MSCs via the establishment of ISR and homeostasis.

Nuclear Mechanics and Cancer Cell Migration.

As a cancer cell invades adjacent tissue, penetrates a basement membrane barrier, or squeezes into a blood capillary, its nucleus can be greatly constricted. Here, we examine: (1) the passive and active deformation of the nucleus during 3D migration; (2) the nuclear structures-namely, the lamina and chromatin-that govern nuclear deformability; (3) the effect of large nuclear deformation on DNA and nuclear factors; and (4) the downstream consequences of mechanically stressing the nucleus. We focus especially on recent studies showing that constricted migration causes nuclear envelope rupture and excess DNA damage, leading to cell cycle suppression, possibly cell death, and ultimately it seems to heritable genomic variation. We first review the latest understanding of nuclear dynamics during cell migration, and then explore the functional effects of nuclear deformation, especially in relation to genome integrity and potentially cancerous mutations.

Rescue of DNA damage after constricted migration reveals a mechano-regulated threshold for cell cycle.

Migration through 3D constrictions can cause nuclear rupture and mislocalization of nuclear proteins, but damage to DNA remains uncertain, as does any effect on cell cycle. Here, myosin II inhibition rescues rupture and partially rescues the DNA damage marker γH2AX, but an apparent block in cell cycle appears unaffected. Co-overexpression of multiple DNA repair factors or antioxidant inhibition of break formation also exert partial effects, independently of rupture. Combined treatments completely rescue cell cycle suppression by DNA damage, revealing a sigmoidal dependence of cell cycle on excess DNA damage. Migration through custom-etched pores yields the same damage threshold, with ∼4-µm pores causing intermediate levels of both damage and cell cycle suppression. High curvature imposed rapidly by pores or probes or else by small micronuclei consistently associates nuclear rupture with dilution of stiff lamin-B filaments, loss of repair factors, and entry from cytoplasm of chromatin-binding cGAS (cyclic GMP-AMP synthase). The cell cycle block caused by constricted migration is nonetheless reversible, with a potential for DNA misrepair and genome variation.

Constricted migration modulates stem cell differentiation.

Tissue regeneration at an injured site depends on proliferation, migration, and differentiation of resident stem or progenitor cells, but solid tissues are often sufficiently dense and constricting that nuclei are highly stressed by migration. In this study, constricted migration of myoblastic cell types and mesenchymal stem cells (MSCs) increases nuclear rupture, increases DNA damage, and modulates differentiation. Fewer myoblasts fuse into regenerating muscle in vivo after constricted migration in vitro, and myodifferentiation in vitro is likewise suppressed. Myosin II inhibition rescues rupture and DNA damage, implicating nuclear forces, while mitosis and the cell cycle are suppressed by constricted migration, consistent with a checkpoint. Although perturbed proliferation fails to explain defective differentiation, nuclear rupture mislocalizes differentiation-relevant MyoD and KU80 (a DNA repair factor), with nuclear entry of the DNA-binding factor cGAS. Human MSCs exhibit similar damage, but osteogenesis increases-which is relevant to bone and to calcified fibrotic tissues, including diseased muscle. Tissue repair can thus be modulated up or down by the curvature of pores through which stem cells squeeze.

Mechanosensing by the Lamina Protects against Nuclear Rupture, DNA Damage, and Cell-Cycle Arrest.

Whether cell forces or extracellular matrix (ECM) can impact genome integrity is largely unclear. Here, acute perturbations (∼1 h) to actomyosin stress or ECM elasticity cause rapid and reversible changes in lamin-A, DNA damage, and cell cycle. The findings are especially relevant to organs such as the heart because DNA damage permanently arrests cardiomyocyte proliferation shortly after birth and thereby eliminates regeneration after injury including heart attack. Embryonic hearts, cardiac-differentiated iPS cells (induced pluripotent stem cells), and various nonmuscle cell types all show that actomyosin-driven nuclear rupture causes cytoplasmic mis-localization of DNA repair factors and excess DNA damage. Binucleation and micronuclei increase as telomeres shorten, which all favor cell-cycle arrest. Deficiencies in lamin-A and repair factors exacerbate these effects, but lamin-A-associated defects are rescued by repair factor overexpression and also by contractility modulators in clinical trials. Contractile cells on stiff ECM normally exhibit low phosphorylation and slow degradation of lamin-A by matrix-metalloprotease-2 (MMP2), and inhibition of this lamin-A turnover and also actomyosin contractility are seen to minimize DNA damage. Lamin-A is thus stress stabilized to mechano-protect the genome.

Nuclear rupture at sites of high curvature compromises retention of DNA repair factors.

The nucleus is physically linked to the cytoskeleton, adhesions, and extracellular matrix-all of which sustain forces, but their relationships to DNA damage are obscure. We show that nuclear rupture with cytoplasmic mislocalization of multiple DNA repair factors correlates with high nuclear curvature imposed by an external probe or by cell attachment to either aligned collagen fibers or stiff matrix. Mislocalization is greatly enhanced by lamin A depletion, requires hours for nuclear reentry, and correlates with an increase in pan-nucleoplasmic foci of the DNA damage marker γH2AX. Excess DNA damage is rescued in ruptured nuclei by cooverexpression of multiple DNA repair factors as well as by soft matrix or inhibition of actomyosin tension. Increased contractility has the opposite effect, and stiff tumors with low lamin A indeed exhibit increased nuclear curvature, more frequent nuclear rupture, and excess DNA damage. Additional stresses likely play a role, but the data suggest high curvature promotes nuclear rupture, which compromises retention of DNA repair factors and favors sustained damage.