Purpose

Pathologic staging of bladder cancer patients remains a challenge. Standard-of-care histology exhibits limited sensitivity in detection of micrometastases, which can increase risk of cancer progression and delay potential adjuvant therapies. Here, we sought to develop a proof of concept novel molecular approach to improve detection of cancer micrometastasis.

Experimental design

We combined fluorescence activated cell sorting and next-generation sequencing and performed whole-exome sequencing of total cancer cells and cancer cell subpopulations in multiple tumor specimens and regional lymph nodes in a single patient with muscle-invasive urothelial carcinoma of the bladder following radical cystectomy.

Results

Mean allele frequency analysis demonstrated a significant correlation between primary tumor cancer cells and cancer cells isolated from the lymph nodes, confirming lymph node disease despite negative pathologic staging. RNA-sequencing revealed intratumoral heterogeneity as well as enrichment for immune system and lipid metabolism gene sets in the micrometastatic cancer cell subpopulations.

Conclusions

Our analysis illustrates how next-generation sequencing of cancer cell subpopulations may be utilized to enrich for cancer cell markers and enhance detection of bladder cancer micrometastases to improve pathologic staging and provide insight into cancer cell biology.

RNA sequencing has allowed high-throughput screening of differential gene expression in many tissues and organisms. Xenopus laevis is a classical embryological and cell-free extract model system, but its genomic sequence had been lacking due to difficulties arising from allotetraploidy. There is currently much excitement surrounding the release of the completed X. laevis genome (version 9.1) by the Joint Genome Institute (JGI), which provides a platform for genome-wide studies. Here we present a deep RNA-seq dataset of transcripts expressed in dorsal and ventral lips of the early Xenopus gastrula embryo using the new genomic information, which was further annotated by blast searches against the human proteome. Overall, our findings confirm previous results from differential screenings using other methods that uncovered classical dorsal genes such as Chordin, Noggin and Cerberus, as well as ventral genes such as Sizzled, Ventx, Wnt8 and Bambi. Complete transcriptome-wide tables of mRNAs suitable for data mining are presented, which include many novel dorsal- and ventral-specific genes. RNA-seq was very quantitative and reproducible, and allowed us to define dorsal and ventral signatures useful for gene set expression analyses (GSEA). As an example of a new gene, we present here data on an organizer-specific secreted protein tyrosine kinase known as Pkdcc (protein kinase domain containing, cytoplasmic) or Vlk (vertebrate lonesome kinase). Overexpression experiments indicate that Pkdcc can act as a negative regulator of Wnt/ β-catenin signaling independently of its kinase activity. We conclude that RNA-Seq in combination with the X. laevis complete genome now available provides a powerful tool for unraveling cell-cell signaling pathways during embryonic induction.

Global DNA demethylation in humans is a fundamental process that occurs in pre-implantation embryos and reversion to naive ground state pluripotent stem cells (PSCs). However, the extent of DNA methylation reprogramming in human germline cells is unknown. Here, we performed whole-genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq) of human prenatal germline cells from 53 to 137 days of development. We discovered that the transcriptome and methylome of human germline is distinct from both human PSCs and the inner cell mass (ICM) of human blastocysts. Using this resource to monitor the outcome of global DNA demethylation with reversion of primed PSCs to the naive ground state, we uncovered hotspots of ultralow methylation at transposons that are protected from demethylation in the germline and ICM. Taken together, the human germline serves as a valuable in vivo tool for monitoring the epigenome of cells that have emerged from a global DNA demethylation event.

Screens for genes that orchestrate neural circuit formation in mammals have been hindered by practical constraints of germline mutagenesis. To overcome these limitations, we combined RNA-seq with somatic CRISPR mutagenesis to study synapse development in the mouse retina. Here synapses occur between cellular layers, forming two multilayered neuropils. The outer neuropil, the outer plexiform layer (OPL), contains synapses made by rod and cone photoreceptor axons on rod and cone bipolar dendrites, respectively. We used RNA-seq to identify selectively expressed genes encoding cell surface and secreted proteins and CRISPR-Cas9 electroporation with cell-specific promoters to assess their roles in OPL development. Among the genes identified in this way are Wnt5a and Wnt5b. They are produced by rod bipolars and activate a non-canonical signaling pathway in rods to regulate early OPL patterning. The approach we use here can be applied to other parts of the brain.

Abstract B cell development in bone marrow is followed by specification into spleen subsets, including marginal zone (MZ) cells. MZ require elaboration of distinct gene expression programs for development. Given their role in gene regulation, its not surprising that microRNAs (miRNAs) influence cell development. Recent work demonstrated that deficiency of NF-κB feedback regulator, Mir146 (miR-146a), led to a range of hematopoietic phenotypes, but B cells have not been extensively characterized. Here, we found miR-146a deficient mice demonstrate a reduction in MZ B cells, likely from a T cell independent developmental block. Utilizing comparative analysis of developmental stage-specific transcriptomes, we show MZ cell differentiation was impaired due to decreases in Notch2 signaling. Further, we discovered that the cell-fate regulatory protein, Numb, is a direct target of miR-146a, and its derepression in miR-146a deficient B cells underlies the decreases in Notch2. Our studies reveal miR-146a-dependent B cell phenotypes regulated by the Numb-Notch2 pathway. Disclosures No relevant conflicts of interest to declare.

B-cell development in the bone marrow is followed by specification into functional subsets in the spleen, including marginal zone (MZ) B-cells. MZ B-cells are classically characterized by T-independent antigenic responses and require the elaboration of distinct gene expression programs for development. Given their role in gene regulation, it is not surprising that microRNAs are important factors in B-cell development. Recent work demonstrated that deficiency of the NFκB feedback regulator, miR-146a, led to a range of hematopoietic phenotypes, but B-cell phenotypes have not been extensively characterized. Here, we found that miR-146a-deficient mice demonstrate a reduction in MZ B-cells, likely from a developmental block. Utilizing high-throughput sequencing and comparative analysis of developmental stage-specific transcriptomes, we determined that MZ cell differentiation was impaired due to decreases in Notch2 signaling. Our studies reveal miR-146a-dependent B-cell phenotypes and highlight the complex role of miR-146a in the hematopoietic system.

The earliest event in Xenopus development is the dorsal accumulation of nuclear β-catenin under the influence of cytoplasmic determinants displaced by fertilization. In this study, a genome-wide approach was used to examine transcription of the 43,673 genes annotated in the Xenopus laevis genome under a variety of conditions that inhibit or promote formation of the Spemann organizer signaling center. Loss of function of β-catenin with antisense morpholinos reproducibly reduced the expression of 247 mRNAs at gastrula stage. Interestingly, only 123 β-catenin targets were enriched on the dorsal side and defined an early dorsal β-catenin gene signature. These genes included several previously unrecognized Spemann organizer components. Surprisingly, only 3 of these 123 genes overlapped with the late Wnt signature recently defined by two other groups using inhibition by Dkk1 mRNA or Wnt8 morpholinos, which indicates that the effects of β-catenin/Wnt signaling in early development are exquisitely regulated by stage-dependent mechanisms. We analyzed transcriptome responses to a number of treatments in a total of 46 RNA-seq libraries. These treatments included, in addition to β-catenin depletion, regenerating dorsal and ventral half-embryos, lithium chloride treatment, and the overexpression of Wnt8, Siamois, and Cerberus mRNAs. Only some of the early dorsal β-catenin signature genes were activated at blastula whereas others required the induction of endomesoderm, as indicated by their inhibition by Cerberus overexpression. These comprehensive data provide a rich resource for analyzing how the dorsal and ventral regions of the embryo communicate with each other in a self-organizing vertebrate model embryo.

Lung squamous cell carcinoma (SCC) is thought to arise from premalignant lesions in the airway epithelium; therefore, studying these lesions is critical for understanding lung carcinogenesis. Previous microarray and sequencing studies designed to discover early biomarkers and therapeutic targets for lung SCC had limited success identifying key driver events in lung carcinogenesis, mostly due to the cellular heterogeneity of patient samples examined and the interindividual variability associated with difficult to obtain airway premalignant lesions and appropriate normal control samples within the same patient. We performed RNA sequencing on laser-microdissected representative cell populations along the SCC pathologic continuum of patient-matched normal basal cells, premalignant lesions, and tumor cells. We discovered transcriptomic changes and identified genomic pathways altered with initiation and progression of SCC within individual patients. We used immunofluorescent staining to confirm gene expression changes in premalignant lesions and tumor cells, including increased expression of SLC2A1, CEACAM5, and PTBP3 at the protein level and increased activation of MYC via nuclear translocation. Cytoband enrichment analysis revealed coordinated loss and gain of expression in chromosome 3p and 3q regions, respectively, during carcinogenesis. This is the first gene expression profiling study of airway premalignant lesions with patient-matched SCC tumor samples. Our results provide much needed information about the biology of premalignant lesions and the molecular changes that occur during stepwise carcinogenesis of SCC, and it highlights a novel approach for identifying some of the earliest molecular changes associated with initiation and progression of lung carcinogenesis within individual patients.

Both basal and submucosal gland (SMG) duct stem cells of the airway epithelium are capable of sphere formation in the in vitro sphere assay, although the efficiency at which this occurs is very low. We sought to improve this efficiency of sphere formation by identifying subpopulations of airway basal stem cells (ABSC) and SMG duct cells based on their aldehyde dehydrogenase (ALDH) activity. ALDH(hi) ABSCs and SMG duct cells were highly enriched for the population of cells that could make spheres, while the co-culture of ALDH(hi) differentiated cells with the ALDH(hi) ABSCs increased their sphere-forming efficiency. Specific ALDH agonists and antagonists were used to show that airway specific ALDH isozymes are important for ABSC proliferation. Pathway analysis of gene expression profiling of ALDH(hi) and ALDH(lo) ABSCs revealed a significant upregulation of the arachidonic acid (AA) metabolism pathway in ALDH(hi) ABSCs. We confirmed the importance of this pathway in the metabolism of proliferating ALDH(hi) ABSCs using bioenergetics studies as well as agonists and antagonists of the AA pathway. These studies could lead to the development of novel strategies for altering ABSC proliferation in the airway epithelium.

Information processing relies on precise patterns of synapses between neurons. The cellular recognition mechanisms regulating this specificity are poorly understood. In the medulla of the Drosophila visual system, different neurons form synaptic connections in different layers. Here, we sought to identify candidate cell recognition molecules underlying this specificity. Using RNA sequencing (RNA-seq), we show that neurons with different synaptic specificities express unique combinations of mRNAs encoding hundreds of cell surface and secreted proteins. Using RNA-seq and protein tagging, we demonstrate that 21 paralogs of the Dpr family, a subclass of immunoglobulin (Ig)-domain containing proteins, are expressed in unique combinations in homologous neurons with different layer-specific synaptic connections. Dpr interacting proteins (DIPs), comprising nine paralogs of another subclass of Ig-containing proteins, are expressed in a complementary layer-specific fashion in a subset of synaptic partners. We propose that pairs of Dpr/DIP paralogs contribute to layer-specific patterns of synaptic connectivity.