Gene Expression Profile of Benign, Intermediate, and Malignant Spitz and Spitzoid Melanocytic Lesions.

Spitz and Spitzoid lesions represent one of the most challenging melanocytic neoplasms in dermatopathology. Nosologic classification has been more recently improved by the discovery of novel molecular drivers, particularly translocations. In the current study, we aimed to use an unbiased approach to explore the gene expression profile of a group of melanocytic Spitz and Spitzoid melanocytic lesions ranging from benign lesions to melanoma, including intermediate lesions such as SPARK nevi and atypical Spitz tumors/melanocytomas. Using unsupervised analysis of gene expression data, we found some distinct hierarchical clusters of lesions, including groups characterized by ALK and NTRK translocations. Few non-ALK translocated tumors demonstrated increased ALK expression, confirmed by immunohistochemistry. Spitz tumors with overlapping features of dysplastic nevi, so-called SPARK nevi, appear to have a common gene expression profile by hierarchical clustering. Finally, weighted gene correlation network analysis identified gene modules variably regulated in subtypes of these cases. Thus, gene expression profiling of Spitz and Spitzoid lesions represents a viable instrument for the characterization of these lesions.

Co-evolution of AR gene copy number and structural complexity in endocrine therapy resistant prostate cancer.

Androgen receptor (AR) inhibition is standard of care for advanced prostate cancer (PC). However, efficacy is limited by progression to castration-resistant PC (CRPC), usually due to AR re-activation via mechanisms that include AR amplification and structural rearrangement. These two classes of AR alterations often co-occur in CRPC tumors, but it is unclear whether this reflects intercellular or intracellular heterogeneity of AR. Resolving this is important for developing new therapies and predictive biomarkers. Here, we analyzed 41 CRPC tumors and 6 patient-derived xenografts (PDXs) using linked-read DNA-sequencing, and identified 7 tumors that developed complex, multiply-rearranged AR gene structures in conjunction with very high AR copy number. Analysis of PDX models by optical genome mapping and fluorescence in situ hybridization showed that AR residing on extrachromosomal DNA (ecDNA) was an underlying mechanism, and was associated with elevated levels and diversity of AR expression. This study identifies co-evolution of AR gene copy number and structural complexity via ecDNA as a mechanism associated with endocrine therapy resistance.

Nuclear corepressors NCOR1/NCOR2 regulate B cell development, maintain genomic integrity and prevent transformation.

The nuclear corepressors NCOR1 and NCOR2 interact with transcription factors involved in B cell development and potentially link these factors to alterations in chromatin structure and gene expression. Herein, we demonstrate that Ncor1/2 deletion limits B cell differentiation via impaired recombination, attenuates pre-BCR signaling and enhances STAT5-dependent transcription. Furthermore, NCOR1/2-deficient B cells exhibited derepression of EZH2-repressed gene modules, including the p53 pathway. These alterations resulted in aberrant Rag1 and Rag2 expression and accessibility. Whole-genome sequencing of Ncor1/2 DKO B cells identified increased number of structural variants with cryptic recombination signal sequences. Finally, deletion of Ncor1 alleles in mice facilitated leukemic transformation, whereas human leukemias with less NCOR1 correlated with worse survival. NCOR1/2 mutations in human leukemia correlated with increased RAG expression and number of structural variants. These studies illuminate how the corepressors NCOR1/2 regulate B cell differentiation and provide insights into how NCOR1/2 mutations may promote B cell transformation.

Tracking Regulatory T Cell Development in the Thymus Using Single-Cell RNA Sequencing/TCR Sequencing.

Recent studies have demonstrated that regulatory T cells (Tregs) develop in the thymus via two pathways involving distinct Treg progenitors (TregP): CD25+FOXP3- (CD25+ TregP) and CD25-FOXP3lo (FOXP3lo TregP) Treg progenitors. To examine this process in more detail, we carried out single-cell RNA sequencing (scRNA-Seq) and TCR-Seq on sorted murine CD4+CD8+ double-positive (DP) thymocytes, CD4+ single-positive (CD4SP) thymocytes, CD25+FOXP3-CD73- TregP, CD25-FOXP3loCD73- TregP, newly generated mature CD25+FOXP3+CD73- Tregs, and FOXP3+CD73+ recirculating/long-term resident Tregs (RT-Tregs). Sorted populations were individually hashtagged and then combined into one scRNA-Seq/TCR-Seq library before sequencing and subsequent analysis. We found that both CD25+ TregP and FOXP3lo TregP arise via an initial agonist-activated state that gives rise to a second transitional stage before differentiating into mature Tregs Using both scRNA-Seq and bulk RNA-Seq on sorted thymocyte subsets, we demonstrate that CD25+ TregP are significantly enriched for Il2 production, suggesting that they are the major source of IL-2 needed to convert TregP into mature Tregs Using TCR-Seq, we found that several TCRs were clearly biased in favor of the conventional or Treg lineages, but that a large fraction of TCRs were found in both these lineages. Finally, we found that RT-Tregs in the thymus are not monomorphic but are composed of multiple distinct subsets and that these RT-Tregs express the most diverse TCR repertoire of all CD4SP thymocytes. Thus, our studies define multiple stages of Treg differentiation within the murine thymus and serve as a resource for future studies on CD4+ thymocyte development and Treg differentiation.

Molecular characterization of biphenotypic epithelioid and plexiform melanoma with deep penetrating nevus-like features.

Tumor heterogeneity is a relevant hallmark of melanoma due to the high mutation burden and immunogenicity commonly encountered. Heterogeneity at the histologic level frequently corresponds to heterogeneity at the molecular level. A better understanding of this feature of malignancy can help refine the development of predictive biomarkers and to define more effective targeted therapies. Here, we describe a case of melanoma displaying a dual phenotype: a DPN-like/plexiform portion in conjunction with a conventional epithelioid morphology. Molecular studies revealed shared BRAF and PTEN mutations in both components but a CTNNB1 mutation was exclusively found in the DPN-like area of the tumor, consistent with the distinct morphology observed. There was considerable heterogeneity in sequence variants identified in the two regions. Gene expression analysis highlighted differentially regulated genes between the two histologies, including a relevant cluster of genes in the receptor tyrosine kinase (RTK) family and related signaling pathways upregulated in the DPN-like/plexiform area.

Opposing transcriptional programs of KLF5 and AR emerge during therapy for advanced prostate cancer.

Endocrine therapies for prostate cancer inhibit the androgen receptor (AR) transcription factor. In most cases, AR activity resumes during therapy and drives progression to castration-resistant prostate cancer (CRPC). However, therapy can also promote lineage plasticity and select for AR-independent phenotypes that are uniformly lethal. Here, we demonstrate the stem cell transcription factor Krüppel-like factor 5 (KLF5) is low or absent in prostate cancers prior to endocrine therapy, but induced in a subset of CRPC, including CRPC displaying lineage plasticity. KLF5 and AR physically interact on chromatin and drive opposing transcriptional programs, with KLF5 promoting cellular migration, anchorage-independent growth, and basal epithelial cell phenotypes. We identify ERBB2 as a point of transcriptional convergence displaying activation by KLF5 and repression by AR. ERBB2 inhibitors preferentially block KLF5-driven oncogenic phenotypes. These findings implicate KLF5 as an oncogene that can be upregulated in CRPC to oppose AR activities and promote lineage plasticity.

Single-cell analysis identifies dynamic gene expression networks that govern B cell development and transformation.

Integration of external signals and B-lymphoid transcription factor activities organise B cell lineage commitment through alternating cycles of proliferation and differentiation, producing a diverse repertoire of mature B cells. We use single-cell transcriptomics/proteomics to identify differentially expressed gene networks across B cell development and correlate these networks with subtypes of B cell leukemia. Here we show unique transcriptional signatures that refine the pre-B cell expansion stages into pre-BCR-dependent and pre-BCR-independent proliferative phases. These changes correlate with reciprocal changes in expression of the transcription factor EBF1 and the RNA binding protein YBX3, that are defining features of the pre-BCR-dependent stage. Using pseudotime analysis, we further characterize the expression kinetics of different biological modalities across B cell development, including transcription factors, cytokines, chemokines, and their associated receptors. Our findings demonstrate the underlying heterogeneity of developing B cells and characterise developmental nodes linked to B cell transformation.

Identification of mutations that cooperate with defects in B cell transcription factors to initiate leukemia.

The transcription factors PAX5, IKZF1, and EBF1 are frequently mutated in B cell acute lymphoblastic leukemia (B-ALL). We demonstrate that compound heterozygous loss of multiple genes critical for B and T cell development drives transformation, including Pax5+/-xEbf1+/-, Pax5+/-xIkzf1+/-, and Ebf1+/-xIkzf1+/- mice for B-ALL, or Tcf7+/-xIkzf1+/- mice for T-ALL. To identify genetic defects that cooperate with Pax5 and Ebf1 compound heterozygosity to initiate leukemia, we performed a Sleeping Beauty (SB) transposon screen that identified cooperating partners including gain-of-function mutations in Stat5b (~65%) and Jak1 (~68%), or loss-of-function mutations in Cblb (61%) and Myb (32%). These findings underscore the role of JAK/STAT5B signaling in B cell transformation and demonstrate roles for loss-of-function mutations in Cblb and Myb in transformation. RNA-Seq studies demonstrated upregulation of a PDK1>SGK3>MYC pathway; treatment of Pax5+/-xEbf1+/- leukemia cells with PDK1 inhibitors blocked proliferation in vitro. In addition, we identified a conserved transcriptional gene signature between human and murine leukemias characterized by upregulation of myeloid genes, most notably involving the GM-CSF pathway, that resemble a B cell/myeloid mixed-lineage leukemia. Thus, our findings identify multiple mechanisms that cooperate with defects in B cell transcription factors to generate either progenitor B cell or mixed B/myeloid-like leukemias.

Gene Expression and Mutational Profile in BAP-1 Inactivated Melanocytic Lesions of Progressive Malignancy from a Patient with Multiple Lesions.

BAP-1 (BRCA1-associated protein 1) inactivated melanocytic lesions are a group of familial or sporadic lesions with unique histology and molecular features. They are of great clinical interest, at least in part due to the potential for malignant transformation and association with a familial cancer predisposition syndrome. Here, we describe a patient with multiple spatially and temporally distinct melanocytic lesions with loss of BAP1 expression by immunohistochemistry. RNA sequencing was performed on three independent lesions spanning the morphologic spectrum: a benign nevus, an atypical tumor, and a melanoma arising from a pre-existing BAP1-inactivated nevus. The three lesions demonstrated largely distinct gene expression and mutational profiles. Gene expression analysis revealed that genes involved in receptor protein kinase pathways were progressively upregulated from nevus to melanoma. Moreover, a clear enrichment of genes regulated in response to UV radiation was found in the melanoma from this patient, as well as upregulation of MAPK pathway-related genes and several transcription factors related to melanomagenesis.

Erratum: Witschen, P.M., et al. Tumor Cell Associated Hyaluronan-CD44 Signaling Promotes Pro-Tumor Inflammation in Breast Cancer. Cancers 2020, 12, 1325.

The authors wish to make the following corrections to this paper [...].

Multiomic Analysis of Subtype Evolution and Heterogeneity in High-Grade Serous Ovarian Carcinoma.

Multiple studies have identified transcriptome subtypes of high-grade serous ovarian carcinoma (HGSOC), but their interpretation and translation are complicated by tumor evolution and polyclonality accompanied by extensive accumulation of somatic aberrations, varying cell type admixtures, and different tissues of origin. In this study, we examined the chronology of HGSOC subtype evolution in the context of these factors using a novel integrative analysis of absolute copy-number analysis and gene expression in The Cancer Genome Atlas complemented by single-cell analysis of six independent tumors. Tumor purity, ploidy, and subclonality were reliably inferred from different genomic platforms, and these characteristics displayed marked differences between subtypes. Genomic lesions associated with HGSOC subtypes tended to be subclonal, implying subtype divergence at later stages of tumor evolution. Subclonality of recurrent HGSOC alterations was evident for proliferative tumors, characterized by extreme genomic instability, absence of immune infiltration, and greater patient age. In contrast, differentiated tumors were characterized by largely intact genome integrity, high immune infiltration, and younger patient age. Single-cell sequencing of 42,000 tumor cells revealed widespread heterogeneity in tumor cell type composition that drove bulk subtypes but demonstrated a lack of intrinsic subtypes among tumor epithelial cells. Our findings prompt the dismissal of discrete transcriptome subtypes for HGSOC and replacement by a more realistic model of continuous tumor development that includes mixtures of subclones, accumulation of somatic aberrations, infiltration of immune and stromal cells in proportions correlated with tumor stage and tissue of origin, and evolution between properties previously associated with discrete subtypes. SIGNIFICANCE: This study infers whether transcriptome-based groupings of tumors differentiate early in carcinogenesis and are, therefore, appropriate targets for therapy and demonstrates that this is not the case for HGSOC.

Tumor Cell Associated Hyaluronan-CD44 Signaling Promotes Pro-Tumor Inflammation in Breast Cancer.

Cancer has been conceptualized as a chronic wound with a predominance of tumor promoting inflammation. Given the accumulating evidence that the microenvironment supports tumor growth, we investigated hyaluronan (HA)-CD44 interactions within breast cancer cells, to determine whether this axis directly impacts the formation of an inflammatory microenvironment. Our results demonstrate that breast cancer cells synthesize and fragment HA and express CD44 on the cell surface. Using RNA sequencing approaches, we found that loss of CD44 in breast cancer cells altered the expression of cytokine-related genes. Specifically, we found that production of the chemokine CCL2 by breast cancer cells was significantly decreased after depletion of either CD44 or HA. In vivo, we found that CD44 deletion in breast cancer cells resulted in a delay in tumor formation and localized progression. This finding was accompanied by a decrease in infiltrating CD206+ macrophages, which are typically associated with tumor promoting functions. Importantly, our laboratory results were supported by human breast cancer patient data, where increased HAS2 expression was significantly associated with a tumor promoting inflammatory gene signature. Because high levels of HA deposition within many tumor types yields a poorer prognosis, our results emphasize that HA-CD44 interactions potentially have broad implications across multiple cancers.

Isolation of enriched small RNA from cell-lysate using on-chip isotachophoresis.

In spite of the growing interest in the roles and applications of small RNAs (sRNAs), sRNA isolation methods are inconsistent, tedious, and dependent on the starting number of cells. In this work, we employ ITP to isolate sRNAs from the cell-lysate of K562 (chronic myelogenous leukemia) cells in a polydimethylsiloxane (PDMS) mesofluidic device. Our method specifically purifies sRNA of <60 nucleotides from lysate of a wide range of cell number spanning from 100 to 1 000 000 cells. We measured the amount of sRNA using the Agilent Bioanalyzer and further verified the extraction efficiency by reverse transcription quantitative PCR. Our method was shown to be more efficient in sRNA extraction than commercial sRNA isolation kits, especially when using smaller numbers of starting cells. Our assay presents a simple and rapid sRNA extraction method with 20 min assay time and no intermediate transfer steps.

Simultaneous RNA purification and size selection using on-chip isotachophoresis with an ionic spacer.

We present an on-chip method for the extraction of RNA within a specific size range from low-abundance samples. We use isotachophoresis (ITP) with an ionic spacer and a sieving matrix to enable size-selection with a high yield of RNA in the target size range. The spacer zone separates two concentrated ITP peaks, the first containing unwanted single nucleotides and the second focusing RNA of the target size range (2-35 nt). Our ITP method excludes >90% of single nucleotides and >65% of longer RNAs (>35 nt). Compared to size selection using gel electrophoresis, ITP-based size-selection yields a 2.2-fold increase in the amount of extracted RNAs within the target size range. We also demonstrate compatibility of the ITP-based size-selection with downstream next generation sequencing. On-chip ITP-prepared samples reveal higher reproducibility of transcript-specific measurements compared to samples size-selected by gel electrophoresis. Our method offers an attractive alternative to conventional sample preparation for sequencing with shorter assay time, higher extraction efficiency and reproducibility. Potential applications of ITP-based size-selection include sequencing-based analyses of small RNAs from low-abundance samples such as rare cell types, samples from fluorescence activated cell sorting (FACS), or limited clinical samples.

Programmable mutually exclusive alternative splicing for generating RNA and protein diversity.

Alternative splicing performs a central role in expanding genomic coding capacity and proteomic diversity. However, programming of splicing patterns in engineered biological systems remains underused. Synthetic approaches thus far have predominantly focused on controlling expression of a single protein through alternative splicing. Here, we describe a modular and extensible platform for regulating four programmable exons that undergo a mutually exclusive alternative splicing event to generate multiple functionally-distinct proteins. We present an intron framework that enforces the mutual exclusivity of two internal exons and demonstrate a graded series of consensus sequence elements of varying strengths that set the ratio of two mutually exclusive isoforms. We apply this framework to program the DNA-binding domains of modular transcription factors to differentially control downstream gene activation. This splicing platform advances an approach for generating diverse isoforms and can ultimately be applied to program modular proteins and increase coding capacity of synthetic biological systems.

A minimum information standard for reproducing bench-scale bacterial cell growth and productivity.

Reproducing, exchanging, comparing, and building on each other's work is foundational to technological advances. Advancing biotechnology calls for reliable reuse of engineered organisms. Reliable reuse of engineered organisms requires reproducible growth and productivity. Here, we identify the experimental factors that have the greatest effect on the growth and productivity of our engineered organisms in order to demonstrate reproducibility for biotechnology. We present a draft of a Minimum Information Standard for Engineered Organism Experiments (MIEO) based on this method. We evaluate the effect of 22 factors on Escherichia coli engineered to produce the small molecule lycopene, and 18 factors on E. coli engineered to produce red fluorescent protein. Container geometry and shaking have the greatest effect on product titer and yield. We reproduce our results under two different conditions of reproducibility: conditions of use (different fractional factorial experiments), and time (48 biological replicates performed on 12 different days over 4 months).

Efficient Production of On-Target Reads for Small RNA Sequencing of Single Cells Using Modified Adapters.

Although single-cell mRNA sequencing has been a powerful tool to explore cellular heterogeneity, the sequencing of small RNA at the single-cell level (sc-sRNA-seq) remains a challenge, as these have no consensus sequence, are relatively low abundant, and are difficult to amplify in a bias-free fashion. We present two methods of single-cell-lysis that enable sc-sRNA-seq. The first method is a chemical-based technique with overnight freezing while the second method leverages on-chip electrical lysis of plasma membrane and physical extraction and separation of cytoplasmic RNA via isotachophoresis. We coupled these two methods with off-chip small RNA library preparation using CleanTag modified adapters to prevent the formation of adapter dimers. We then demonstrated sc-sRNA-seq with single K562 human leukemic cells. Our approaches offer a relatively short hands-on time of 6 h and efficient generation of on-target reads. The sc-sRNA-seq with our approaches showed detection of miRNA with various abundances ranging from 16 000 copies/cell to about 10 copies/cell. We anticipate this approach will create a new opportunity to explore cellular heterogeneity through small RNA expression.

Evaluation of the External RNA Controls Consortium (ERCC) reference material using a modified Latin square design.

BACKGROUND: Highly multiplexed assays for quantitation of RNA transcripts are being used in many areas of biology and medicine. Using data generated by these transcriptomic assays requires measurement assurance with appropriate controls. Methods to prototype and evaluate multiple RNA controls were developed as part of the External RNA Controls Consortium (ERCC) assessment process. These approaches included a modified Latin square design to provide a broad dynamic range of relative abundance with known differences between four complex pools of ERCC RNA transcripts spiked into a human liver total RNA background. RESULTS: ERCC pools were analyzed on four different microarray platforms: Agilent 1- and 2-color, Illumina bead, and NIAID lab-made spotted microarrays; and two different second-generation sequencing platforms: the Life Technologies 5500xl and the Illumina HiSeq 2500. Individual ERCC controls were assessed for reproducible performance in signal response to concentration among the platforms. Most demonstrated linear behavior if they were not located near one of the extremes of the dynamic range. Performance issues with any individual ERCC transcript could be attributed to detection limitations, platform-specific target probe issues, or potential mixing errors. Collectively, these pools of spike-in RNA controls were evaluated for suitability as surrogates for endogenous transcripts to interrogate the performance of the RNA measurement process of each platform. The controls were useful for establishing the dynamic range of the assay, as well as delineating the useable region of that range where differential expression measurements, expressed as ratios, would be expected to be accurate. CONCLUSIONS: The modified Latin square design presented here uses a composite testing scheme for the evaluation of multiple performance characteristics: linear performance of individual controls, signal response within dynamic range pools of controls, and ratio detection between pairs of dynamic range pools. This compact design provides an economical sample format for the evaluation of multiple external RNA controls within a single experiment per platform. These results indicate that well-designed pools of RNA controls, spiked into samples, provide measurement assurance for endogenous gene expression studies.

Assessing technical performance in differential gene expression experiments with external spike-in RNA control ratio mixtures.

There is a critical need for standard approaches to assess, report and compare the technical performance of genome-scale differential gene expression experiments. Here we assess technical performance with a proposed standard 'dashboard' of metrics derived from analysis of external spike-in RNA control ratio mixtures. These control ratio mixtures with defined abundance ratios enable assessment of diagnostic performance of differentially expressed transcript lists, limit of detection of ratio (LODR) estimates and expression ratio variability and measurement bias. The performance metrics suite is applicable to analysis of a typical experiment, and here we also apply these metrics to evaluate technical performance among laboratories. An interlaboratory study using identical samples shared among 12 laboratories with three different measurement processes demonstrates generally consistent diagnostic power across 11 laboratories. Ratio measurement variability and bias are also comparable among laboratories for the same measurement process. We observe different biases for measurement processes using different mRNA-enrichment protocols.

Proteomic and physiological experiments to test Thermotoga neapolitana constraint-based model hypotheses of carbon source utilization.

Constraint-based models of biochemical reaction networks require experimental validation to test model-derived hypotheses and iteratively improve the model. Physiological and proteomic analysis of Thermotoga neapolitana growth on cellotetraose was conducted to identify gene products related to growth on cellotetraose to improve a constraint-based model of T. neapolitana central carbon metabolism with incomplete cellotetraose pathways. In physiological experiments comparing cellotetraose to cellobiose and glucose as growth substrates, product formation yields on cellotetraose, cellobiose, and glucose were similar; however cell yields per mol carbon consumed were higher on cellotetraose than on cellobiose or glucose. Proteomic analysis showed increased expression of several proteins from cells grown on cellotetraose compared with glucose cell cultures, including cellobiose phosphorylase (CTN_0783), endo-1,4-ß-glucosidase (CTN_1106), and an ATP-binding protein (CTN_1296). The CTN_1296 gene product should be evaluated further for participation in cellotetraose metabolism and is included as one of two hypothetical gene-protein-reaction associations in the T. neapolitana constraint-based model to reinstate cellotetraose metabolism in model simulations.