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.

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.

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.

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.

Modeling Notch-Induced Tumor Cell Survival in the Drosophila Ovary Identifies Cellular and Transcriptional Response to Nuclear NICD Accumulation.

Notch is a conserved developmental signaling pathway that is dysregulated in many cancer types, most often through constitutive activation. Tumor cells with nuclear accumulation of the active Notch receptor, NICD, generally exhibit enhanced survival while patients experience poorer outcomes. To understand the impact of NICD accumulation during tumorigenesis, we developed a tumor model using the Drosophila ovarian follicular epithelium. Using this system we demonstrated that NICD accumulation contributed to larger tumor growth, reduced apoptosis, increased nuclear size, and fewer incidents of DNA damage without altering ploidy. Using bulk RNA sequencing we identified key genes involved in both a pre- and post- tumor response to NICD accumulation. Among these are genes involved in regulating double-strand break repair, chromosome organization, metabolism, like raptor, which we experimentally validated contributes to early Notch-induced tumor growth. Finally, using single-cell RNA sequencing we identified follicle cell-specific targets in NICD-overexpressing cells which contribute to DNA repair and negative regulation of apoptosis. This valuable tumor model for nuclear NICD accumulation in adult Drosophila follicle cells has allowed us to better understand the specific contribution of nuclear NICD accumulation to cell survival in tumorigenesis and tumor 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.

A single-cell atlas of adult Drosophila ovary identifies transcriptional programs and somatic cell lineage regulating oogenesis.

Oogenesis is a complex developmental process that involves spatiotemporally regulated coordination between the germline and supporting, somatic cell populations. This process has been modeled extensively using the Drosophila ovary. Although different ovarian cell types have been identified through traditional means, the large-scale expression profiles underlying each cell type remain unknown. Using single-cell RNA sequencing technology, we have built a transcriptomic data set for the adult Drosophila ovary and connected tissues. Using this data set, we identified the transcriptional trajectory of the entire follicle-cell population over the course of their development from stem cells to the oogenesis-to-ovulation transition. We further identify expression patterns during essential developmental events that take place in somatic and germline cell types such as differentiation, cell-cycle switching, migration, symmetry breaking, nurse-cell engulfment, egg-shell formation, and corpus luteum signaling. Extensive experimental validation of unique expression patterns in both ovarian and nearby, nonovarian cells also led to the identification of many new cell type-and stage-specific markers. The inclusion of several nearby tissue types in this data set also led to our identification of functional convergence in expression between distantly related cell types such as the immune-related genes that were similarly expressed in immune cells (hemocytes) and ovarian somatic cells (stretched cells) during their brief phagocytic role in nurse-cell engulfment. Taken together, these findings provide new insight into the temporal regulation of genes in a cell-type specific manner during oogenesis and begin to reveal the relatedness in expression between cell and tissues types.

CRL4Mahj E3 ubiquitin ligase promotes neural stem cell reactivation.

The ability of neural stem cells (NSCs) to transit between quiescence and proliferation is crucial for brain development and homeostasis. Drosophila Hippo pathway maintains NSC quiescence, but its regulation during brain development remains unknown. Here, we show that CRL4Mahj, an evolutionarily conserved E3 ubiquitin ligase, is essential for NSC reactivation (exit from quiescence). We demonstrate that damaged DNA-binding protein 1 (DDB1) and Cullin4, two core components of Cullin4-RING ligase (CRL4), are intrinsically required for NSC reactivation. We have identified a substrate receptor of CRL4, Mahjong (Mahj), which is necessary and sufficient for NSC reactivation. Moreover, we show that CRL4Mahj forms a protein complex with Warts (Wts/large tumor suppressor [Lats]), a kinase of the Hippo signaling pathway, and Mahj promotes the ubiquitination of Wts. Our genetic analyses further support the conclusion that CRL4Mahj triggers NSC reactivation by inhibition of Wts. Given that Cullin4B mutations cause mental retardation and cerebral malformation, similar regulatory mechanisms may be applied to the human brain.

The Hox proteins Ubx and AbdA collaborate with the transcription pausing factor M1BP to regulate gene transcription.

In metazoans, the pausing of RNA polymerase II at the promoter (paused Pol II) has emerged as a widespread and conserved mechanism in the regulation of gene transcription. While critical in recruiting Pol II to the promoter, the role transcription factors play in transitioning paused Pol II into productive Pol II is, however, little known. By studying how Drosophila Hox transcription factors control transcription, we uncovered a molecular mechanism that increases productive transcription. We found that the Hox proteins AbdA and Ubx target gene promoters previously bound by the transcription pausing factor M1BP, containing paused Pol II and enriched with promoter-proximal Polycomb Group (PcG) proteins, yet lacking the classical H3K27me3 PcG signature. We found that AbdA binding to M1BP-regulated genes results in reduction in PcG binding, the release of paused Pol II, increases in promoter H3K4me3 histone marks and increased gene transcription. Linking transcription factors, PcG proteins and paused Pol II states, these data identify a two-step mechanism of Hox-driven transcription, with M1BP binding leading to Pol II recruitment followed by AbdA targeting, which results in a change in the chromatin landscape and enhanced transcription.

Systematic analysis reveals tumor-enhancing and -suppressing microRNAs in Drosophila epithelial tumors.

Despite their emergence as an important class of noncoding RNAs involved in cancer cell transformation, invasion, and migration, the precise role of microRNAs (miRNAs) in tumorigenesis remains elusive. To gain insights into how miRNAs contribute to primary tumor formation, we conducted an RNA sequencing (RNA-Seq) analysis of Drosophila wing disc epithelial tumors induced by knockdown of a neoplastic tumor-suppressor gene (nTSG) lethal giant larvae (lgl), combined with overexpression of an active form of oncogene Ras (RasV12 ), and identified 51 mature miRNAs that changed significantly in tumorous discs. Followed by in vivo tumor enhancer and suppressor screens in sensitized genetic backgrounds, we identified 10 tumor-enhancing (TE) miRNAs and 11 tumor-suppressing (TS) miRNAs that contributed to the nTSG defect-induced tumorigenesis. Among these, four TE and three TS miRNAs have human homologs. From this study, we also identified 29 miRNAs that individually had no obvious role in enhancing or alleviating tumorigenesis despite their changed expression levels in nTSG tumors. This systematic analysis, which includes both RNA-Seq and in vivo functional studies, helps to categorize miRNAs into different groups based on their expression profile and functional relevance in epithelial tumorigenesis, whereas the evolutionarily conserved TE and TS miRNAs provide potential therapeutic targets for epithelial tumor treatment.

The SWI/SNF Complex Protein Snr1 Is a Tumor Suppressor in Drosophila Imaginal Tissues.

Components of the SWI/SNF chromatin-remodeling complex are among the most frequently mutated genes in various human cancers, yet only SMARCB1/hSNF5, a core member of the SWI/SNF complex, is mutated in malignant rhabdoid tumors (MRT). How SMARCB1/hSNF5 functions differently from other members of the SWI/SNF complex remains unclear. Here, we use Drosophila imaginal epithelial tissues to demonstrate that Snr1, the conserved homolog of human SMARCB1/hSNF5, prevents tumorigenesis by maintaining normal endosomal trafficking-mediated signaling cascades. Removal of Snr1 resulted in neoplastic tumorigenic overgrowth in imaginal epithelial tissues, whereas depletion of any other members of the SWI/SNF complex did not induce similar phenotypes. Unlike other components of the SWI/SNF complex that were detected only in the nucleus, Snr1 was observed in both the nucleus and the cytoplasm. Aberrant regulation of multiple signaling pathways, including Notch, JNK, and JAK/STAT, was responsible for tumor progression upon snr1-depletion. Our results suggest that the cytoplasmic Snr1 may play a tumor suppressive role in Drosophila imaginal tissues, offering a foundation for understanding the pivotal role of SMARCB1/hSNF5 in suppressing MRT during early childhood. Cancer Res; 77(4); 862-73. ©2017 AACR.

Maternal AP determinants in the Drosophila oocyte and embryo.

An animal embryo cannot initiate its journey of forming a new life on its own. It must rely on maternally provided resources and inputs to kick-start its developmental process. In Drosophila, the initial polarities of the embryo along both the anterior-posterior (AP) and dorsal-ventral (DV) axes are also specified by maternal determinants. Over the past several decades, genetic and molecular studies have identified and characterized such determinants, as well as the zygotic genetic regulatory networks that control patterning in the early embryo. Extensive studies of oogenesis have also led to a detailed knowledge of the cellular and molecular interactions that control the formation of a mature egg. Despite these efforts, oogenesis and embryogenesis have been studied largely as separate problems, except for qualitative aspects with regard to maternal regulation of the asymmetric localization of maternal determinants. Can oogenesis and embryogenesis be viewed from a unified perspective at a quantitative level, and can that improve our understanding of how robust embryonic patterning is achieved? Here, we discuss the basic knowledge of the regulatory mechanisms controlling oogenesis and embryonic patterning along the AP axis. We explore properties of the maternal Bicoid gradient in relation to embryo size in search for a unified framework for robust AP patterning. WIREs Dev Biol 2016, 5:562-581. doi: 10.1002/wdev.235 For further resources related to this article, please visit the WIREs website.

RNA helicase Belle/DDX3 regulates transgene expression in Drosophila.

Belle (Bel), the Drosophila homolog of the yeast DEAD-box RNA helicase DED1 and human DDX3, has been shown to be required for oogenesis and female fertility. Here we report a novel role of Bel in regulating the expression of transgenes. Abrogation of Bel by mutations or RNAi induces silencing of a variety of P-element-derived transgenes. This silencing effect depends on downregulation of their RNA levels. Our genetic studies have revealed that the RNA helicase Spindle-E (Spn-E), a nuage RNA helicase that plays a crucial role in regulating RNA processing and PIWI-interacting RNA (piRNA) biogenesis in germline cells, is required for loss-of-bel-induced transgene silencing. Conversely, Bel abrogation alleviates the nuage-protein mislocalization phenotype in spn-E mutants, suggesting a competitive relationship between these two RNA helicases. Additionally, disruption of the chromatin remodeling factor Mod(mdg4) or the microRNA biogenesis enzyme Dicer-1 (Dcr-1) also alleviates the transgene-silencing phenotypes in bel mutants, suggesting the involvement of chromatin remodeling and microRNA biogenesis in loss-of-bel-induced transgene silencing. Finally we show that genetic inhibition of Bel function leads to de novo generation of piRNAs from the transgene region inserted in the genome, suggesting a potential piRNA-dependent mechanism that may mediate transgene silencing as Bel function is inhibited.

E(y)1/TAF9 mediates the transcriptional output of Notch signaling in Drosophila.

Transcriptional activation of Notch signaling targets requires the formation of a ternary complex that involves the intracellular domain of the Notch receptor (NICD), DNA-binding protein Suppressor of Hairless [Su(H), RPBJ in mammals] and coactivator Mastermind (Mam). Here, we report that E(y)1/TAF9, a component of the transcription factor TFIID complex, interacts specifically with the NICD-Su(H)-Mam complex to facilitate the transcriptional output of Notch signaling. We identified E(y)1/TAF9 in a large-scale in vivo RNA interference (RNAi) screen for genes that are involved in a Notch-dependent mitotic-to-endocycle transition in Drosophila follicle cells. Knockdown of e(y)1/TAF9 displayed Notch-mutant-like phenotypes and defects in target gene and activity reporter expression in both the follicle cells and wing imaginal discs. Epistatic analyses in these two tissues indicated that E(y)1/TAF9 functions downstream of Notch cleavage. Biochemical studies in S2 cells demonstrated that E(y)1/TAF9 physically interacts with the transcriptional effectors of Notch signaling Su(H) and NICD. Taken together, our data suggest that the association of the NICD-Su(H)-Mastermind complex with E(y)1/TAF9 in response to Notch activation recruits the transcription initiation complex to induce Notch target genes, coupling Notch signaling with the transcription machinery.

Uif, a large transmembrane protein with EGF-like repeats, can antagonize Notch signaling in Drosophila.

BACKGROUND: Notch signaling is a highly conserved pathway in multi-cellular organisms ranging from flies to humans. It controls a variety of developmental processes by stimulating the expression of its target genes in a highly specific manner both spatially and temporally. The diversity, specificity and sensitivity of the Notch signaling output are regulated at distinct levels, particularly at the level of ligand-receptor interactions. METHODOLOGY/PRINCIPAL FINDINGS: Here, we report that the Drosophila gene uninflatable (uif), which encodes a large transmembrane protein with eighteen EGF-like repeats in its extracellular domain, can antagonize the canonical Notch signaling pathway. Overexpression of Uif or ectopic expression of a neomorphic form of Uif, Uif*, causes Notch signaling defects in both the wing and the sensory organ precursors. Further experiments suggest that ectopic expression of Uif* inhibits Notch signaling in cis and acts at a step that is dependent on the extracellular domain of Notch. Our results suggest that Uif can alter the accessibility of the Notch extracellular domain to its ligands during Notch activation. CONCLUSIONS/SIGNIFICANCE: Our study shows that Uif can modulate Notch activity, illustrating the importance of a delicate regulation of this signaling pathway for normal patterning.

CAF-1 is essential for Drosophila development and involved in the maintenance of epigenetic memory.

DNA synthesis during S-phase and upon DNA repair is accompanied by chromatin assembly. The chromatin assembly factor CAF-1 has been biochemically well-characterized to deposit histones onto newly synthesized DNA. To gain insights into the in vivo functions of CAF-1 in Drosophila, we generated null mutants of the largest subunit of dCAF-1, dCAF-1-p180. We show that, unlike CAF-1 mutant yeast, dCAF-1-p180 mutant flies are hemizygous lethal. Removal of maternal dCAF-1-p180 activity by germline clones blocks oogenesis. Tissue-specific deletion of dCAF-1-p180 in the eye primordia disrupts eye development. In addition, reduction of dCAF-1-p180 activity suppresses gene silencing at heterochromatin and antagonizes Polycomb-mediated cell fate determination. Furthermore, heterozygous dCAF-1-p180 mutant flies display an increased sensitivity to gamma-irradiation and a reduced efficiency in recombinational double strand break (DSB) repair. Our experiments also show that human hCAF-1-p150 can rescue the dCAF-1-p180 mutant flies, demonstrating a functional conservation of eukaryotic CAF-1 activities in vivo. Together, our results establish that dCAF-1-p180 is an essential gene for Drosophila development and further underscore the importance of dCAF-1 in regulating gene expression and DNA repair in vivo.

Glucose-6-phosphate dehydrogenase plays a pivotal role in nitric oxide-involved defense against oxidative stress under salt stress in red kidney bean roots.

The pivotal role of glucose-6-phosphate dehydrogenase (G-6-PDH)-mediated nitric oxide (NO) production in the tolerance to oxidative stress induced by 100 mM NaCl in red kidney bean (Phaseolus vulgaris) roots was investigated. The results show that the G-6-PDH activity was enhanced rapidly in the presence of NaCl and reached a maximum at 100 mM. Western blot analysis indicated that the increase of G-6-PDH activity in the red kidney bean roots under 100 mM NaCl was mainly due to the increased content of the G-6-PDH protein. NO production and nitrate reductase (NR) activity were also induced by 100 mM NaCl. The NO production was reduced by NaN(3) (an NR inhibitor), but not affected by N(omega)-nitro-L-arginine (L-NNA) (an NOS inhibitor). Application of 2.5 mM Na(3)PO(4), an inhibitor of G-6-PDH, blocked the increase of G-6-PDH and NR activity, as well as NO production in red kidney bean roots under 100 mM NaCl. The activities of antioxidant enzymes in red kidney bean roots increased in the presence of 100 mM NaCl or sodium nitroprusside (SNP), an NO donor. The increased activities of all antioxidant enzymes tested at 100 mM NaCl were completely inhibited by 2.5 mM Na(3)PO(4). Based on these results, we conclude that G-6-PDH plays a pivotal role in NR-dependent NO production, and in establishing tolerance of red kidney bean roots to salt stress.