Abstract

Spatial transcriptome analysis of formalin-fixed paraffin-embedded (FFPE) tissues using RNA-sequencing (RNA-seq) provides interactive information on morphology and gene expression, which is useful for clinical applications. However, despite the advantages of long-term storage at room temperature, FFPE tissues may be severely damaged by methylene crosslinking and provide less gene information than fresh-frozen tissues. In this study, we proposed a sensitive FFPE micro-tissue RNA-seq method that combines the punching of tissue sections (diameter: 100 μm) and the direct construction of RNA-seq libraries. We evaluated a method using mouse liver tissues at 2 years after fixation and embedding and detected approximately 7,000 genes in micro-punched tissue-spots (thickness: 10 μm), similar to that detected with purified total RNA (2.5 ng) equivalent to the several dozen cells in the spot. We applied this method to clinical FFPE specimens of lung cancer that had been fixed and embedded 6 years prior, and found that it was possible to determine characteristic gene expression in the microenvironment containing tumor and non-tumor cells of different morphologies. This result indicates that spatial gene expression analysis of the tumor microenvironment is feasible using FFPE tissue sections stored for extensive periods in medical facilities.

Abstract

Gene expression analysis at the single-cell level by next generation sequencing has revealed the existence of clonal dissemination and microheterogeneity in cancer metastasis. The current spatial analysis technologies can elucidate the heterogeneity of cell–cell interactions in situ. To reveal the regional and expressional heterogeneity in primary tumors and metastases, we performed transcriptomic analysis of microtissues dissected from a triple-negative breast cancer (TNBC) cell line MDA-MB-231 xenograft model with our automated tissue microdissection punching technology. This multiple-microtissue transcriptome analysis revealed three cancer cell-type clusters in the primary tumor and axillary lymph node metastasis, two of which were cancer stem cell (CSC)-like clusters (CD44/MYC-high, HMGA1-high). Reanalysis of public single-cell RNA-seq (scRNA-seq) datasets confirmed that the two CSC-like populations existed both in TNBC xenograft models and TNBC patients. In addition, the gene signature of the HMGA1-high CSC-like cluster has the potential to serve as a novel biomarker for diagnosis. The diversity of these multiple CSC-like populations may cause differential anticancer drug resistance, increasing the difficulty of curing this cancer.

The central nucleus of the amygdala (CeA) and the lateral division of the bed nucleus of the stria terminalis (BNST) are the two major nuclei of the central extended amygdala that plays essential roles in threat processing, responsible for emotional states such as fear and anxiety. While some studies suggested functional differences between these nuclei, others showed anatomical and neurochemical similarities. Despite their complex subnuclear organization, subnuclei-specific functional impact on behavior and their underlying molecular profiles remain obscure. We here constitutively inhibited neurotransmission of protein kinase C-δ-positive (PKCδ+) neurons—a major cell type of the lateral subdivision of the CeA (CeL) and the oval nucleus of the BNST (BNSTov)—and found striking subnuclei-specific effects on fear- and anxiety-related behaviors, respectively. To obtain molecular clues for this dissociation, we conducted RNA sequencing in subnuclei-targeted micropunch samples. The CeL and the BNSTov displayed similar gene expression profiles at the basal level; however, both displayed differential gene expression when animals were exposed to fear-related stimuli, with a more robust expression change in the CeL. These findings provide novel insights into the molecular makeup and differential engagement of distinct subnuclei of the extended amygdala, critical for regulation of threat processing.

Recent studies have shown that mutation at Ser522 causes inhibition of collapsin response mediator protein 2 (CRMP2) phosphorylation and induces axon elongation and partial recovery of the lost sensorimotor function after spinal cord injury (SCI). We aimed to reveal the intracellular mechanism in axotomized neurons in the CRMP2 knock-in (CRMP2KI) mouse model by performing transcriptome analysis in mouse sensorimotor cortex using micro-dissection punching system. Prior to that, we analyzed the structural pathophysiology in axotomized or neighboring neurons after SCI and found that somatic atrophy and dendritic spine reduction in sensorimotor cortex were suppressed in CRMP2KI mice. Further analysis of the transcriptome has aided in the identification of four hemoglobin genes Hba-a1, Hba-a2, Hbb-bs, and Hbb-bt that are significantly upregulated in wild-type mice with concomitant upregulation of genes involved in the oxidative phosphorylation and ribosomal pathways after SCI. However, we observed substantial upregulation in channel activity genes and downregulation of genes regulating vesicles, synaptic function, glial cell differentiation in CRMP2KI mice. Moreover, the transcriptome profile of CRMP2KI mice has been discussed wherein energy metabolism and neuronal pathways were found to be differentially regulated. Our results showed that CRMP2KI mice displayed improved SCI pathophysiology not only via microtubule stabilization in neurons, but also possibly via the whole metabolic system in the central nervous system, response changes in glial cells, and synapses. Taken together, we reveal new insights on SCI pathophysiology and the regenerative mechanism of central nervous system by the inhibition of CRMP2 phosphorylation at Ser522. All these experiments were performed in accordance with the guidelines of the Institutional Animal Care and Use Committee at Waseda University, Japan (2017-A027 approved on March 21, 2017; 2018-A003 approved on March 25, 2018; 2019-A026 approved on March 25, 2019).

Humanized mouse models are attractive experimental models for analyzing the development and functions of human dendritic cells (DCs) in vivo. Although various types of DC subsets, including DC type 3 (DC3s), have been identified in humans, it remains unclear whether humanized mice can reproduce heterogeneous DC subsets. CD14, classically known as a monocyte/macrophage marker, is reported as an indicator of DC3s. We previously observed that some CD14+ myeloid cells expressed CD1c, a pan marker for bona fide conventional DC2 (cDC2s), in humanized mouse models in which human FLT3L and GM-CSF genes were transiently expressed using in vivo transfection (IVT). Here, we aimed to elucidate the identity of CD14+CD1c+ DC-like cells in humanized mouse models. We found that CD14+CD1c+ cells were phenotypically different from cDC2s; CD14+CD1c+ cells expressed CD163 but not CD5, whereas cDC2s expressed CD5 but not CD163. Furthermore, CD14+CD1c+ cells primed and polarized naïve CD4+ T cells toward IFN-γ+ Th1 cells more profoundly than cDC2s. Transcriptional analysis revealed that CD14+CD1c+ cells expressed several DC3-specific transcripts, such as CD163, S100A8, and S100A9, and were clearly segregated from cDC2s and monocytes. When lipopolysaccharide was administered to the humanized mice, the frequency of CD14+CD1c+ cells producing IL-6 and TNF-α was elevated, indicating a pro-inflammatory signature. Thus, humanized mice are able to sustain development of functional CD14+CD1c+ DCs, which are equivalent to DC3 subset observed in humans, and they could be useful for analyzing the development and function of DC3s in vivo.

Spatial transcriptomics is useful for understanding the molecular organization of a tissue and providing insights into cellular function in a morphological context. In order to obtain reproducible results in spatial transcriptomics, we have to maintain tissue morphology and RNA molecule stability during the image acquisition and biomolecule collection processes. Here, we developed a tissue processing method for robust and reproducible RNA-seq from tissue microdissection samples. In this method, we suppressed RNA degradation in fresh-frozen tissue specimens by dehydration fixation and effectively collected a small amount of RNA molecules from microdissection samples by magnetic beads. We demonstrated the spatial transcriptome analysis of the mouse liver and brain in serial microdissection samples (100 μm in a diameter and 10 μm in thickness) produced by a microdissection punching system. Using our method, we could prevent RNA degradation at room temperature and effectively produce a sequencing library with Smart-seq2. This resulted in reproducible sequence read mapping in exon regions and the detection of more than 2000 genes compared to non-fixed samples in the RNA-seq analysis. Our method would be applied to various transcriptome analyses, providing the information for region specific gene expression in tissue specimens.

A simple and robust ultra-small four-color-fluorescence detection system was developed by integrating its components, namely, a four-capillary array, an injection-molded-plastic four-lens array, a four-dichroic-mirror array, and a CMOS sensor, as one device. The developed system was applied to a high-dynamic-range capillary-array electrophoresis (HiDy CE) to quantify a rare EGFR mutant (MT) of exon 19 deletion in a large excess of EGFR wild type (WT). Samples with serially diluted MT and constant-concentration WT were co-amplified by competitive PCR and subjected to HiDy CE. The MT peak in each electropherogram was then compared to the WT peak. As a result, MT was quantified with high-sensitivity (LOD of 0.004% MT/WT) and four-orders-of-magnitude dynamic range (0.01-100% MT/WT) by HiDy CE. Moreover, compared with existing methods, HiDy CE achieves higher speed, higher sample throughput, and lower consumable cost per sample. It has therefore great potential to be used in clinical practice.

Single-cell RNA-seq is a powerful tool for revealing heterogeneity in cancer cells. However, each of the current single-cell RNA-seq platforms has inherent advantages and disadvantages. Here, we show that combining the different single-cell RNA-seq platforms can be an effective approach to obtaining complete information about expression differences and a sufficient cellular population to understand transcriptional heterogeneity in cancers. We demonstrate that it is possible to estimate missing expression information. We further demonstrate that even in the cases where precise information for an individual gene cannot be inferred, the activity of given transcriptional modules can be analyzed. Interestingly, we found that two distinct transcriptional modules, one associated with the Aurora kinase gene and the other with the DUSP gene, are aberrantly regulated in a minor population of cells and may thus contribute to the possible emergence of dormancy or eventual drug resistance within the population.

Analyses of gene expressions in single cells are important for understanding detailed biological phenomena. Here, a highly sensitive and accurate method by sequencing (called "bead-seq") to obtain a whole gene expression profile for a single cell is proposed. A key feature of the method is to use a complementary DNA (cDNA) library on magnetic beads, which enables adding washing steps to remove residual reagents in a sample preparation process. By adding the washing steps, the next steps can be carried out under the optimal conditions without losing cDNAs. Error sources were carefully evaluated to conclude that the first several steps were the key steps. It is demonstrated that bead-seq is superior to the conventional methods for single-cell gene expression analyses in terms of reproducibility, quantitative accuracy, and biases caused during sample preparation and sequencing processes.

Analysis of single-cell gene expression promises a more precise understanding of molecular mechanisms of a living system. Most techniques only allow studies of the expressions for limited numbers of gene species. When amplification of cDNA was carried out for analysing more genes, amplification biases were frequently reported. A non-biased and efficient global-amplification method, which uses a single-cell cDNA library immobilized on beads, was developed for analysing entire gene expressions for single cells. Every step in this analysis from reverse transcription to cDNA amplification was optimized. By removing degrading excess primers, the bias due to the digestion of cDNA was prevented. Since the residual reagents, which affect the efficiency of each subsequent reaction, could be removed by washing beads, the conditions for uniform and maximized amplification of cDNAs were achieved. The differences in the amplification rates for randomly selected eight genes were within 1.5-folds, which could be negligible for most of the applications of single-cell analysis. The global amplification gives a large amount of amplified cDNA (>100 μg) from a single cell (2-pg mRNA), and that amount is enough for downstream analysis. The proposed global-amplification method was used to analyse transcript ratios of multiple cDNA targets (from several copies to several thousand copies) quantitatively.

We studied the fundamental instrumental issues relevant to a capillary-based integrated system to measure expression of a specific gene directly from cells. Samples were introduced into a capillary by use of a syringe pump. All reactions were carried out in a microthermocycler, where a part of the capillary having 1 microL inner volume was used as a reaction vessel. First, cells were lysed by heating to release RNA, followed by deoxyribonuclease (DNase) treatment. Then, reverse transcription-polymerase chain reaction (RT-PCR) was performed to obtain amplified products from the targeted mRNA. Finally, the product was verified by capillary electrophoresis (CE) with laser-induced fluorescence detection. The whole protocol was completed in the system in 3 h. PCR product from beta-actin mRNA in 16 human lymphoblast cells was obtained with a signal-to-noise ratio (S/N) of 3400 +/- 730 (n = 3). Therefore, the system is reproducible and sensitive enough to measure gene expression from a single cell. We show that the amplified fragment from breast cancer-specific mRNA was obtained from cells of breast cancer cell line, but was not obtained from cells of hepatoma cell line. These results therefore lay the foundations for future CE or microchip instrumentation for high-throughput automated gene-expression analysis.

We have developed a novel high-performance quantitative assay for unamplified nucleic acids that is based on single-molecule imaging. The apparatus is a simple but highly sensitive single-molecule detection system that uses a normal CCD camera instead of an image-intensified CCD camera. After the DNA molecules in a sample were labeled with YOYO-1, they were induced to migrate electrophoretically in a polymer solution and imaged. No chemical or biochemical amplification was required. Direct quantitation of the sample by counting molecules was possible, because the number counted over the measurement period was directly proportional to the concentration of DNA molecules in the sample. Nonspecifically labeled impurities that would degrade the sensitivity of the assay were successfully reduced and discriminated from the DNA molecules by differences in electrophoretic mobility. By using beta-actin DNA (838 bp) as a model sample, we demonstrate that this protocol was fast (10-min measurement period), highly sensitive (limit of quantitation: approximately 10(3) copies/sample, or 3 x 10(-16) M), quantitative, and covered a wide linear dynamic range (approximately 10(4)). This high-performance assay promises to be a powerful technology for the quantitation of specific varieties of mRNA in the study of gene functions and diseases and in the clinical detection of mutant cells.

A blood assay for detection of lung cancer biomarkers could significantly improve cancer patient prognosis and survival rates. Amplified fragment length polymorphism-differential display (AFLP-DD) was used to identify gene transcripts found in lung cancer tissue and the peripheral blood of lung cancer patients. The clones were evaluated for gene expression in lung cancer tissue, peripheral blood of lung cancer patients and healthy volunteers' blood. The isolated gene transcript clones were found to be from the syndecan 1 gene, collagen 1 gene and 2 novel genes. All 4 transcripts were expressed in normal lung tissue, 4 cultured primary lung cells and 6 lung cancer cell lines. RNA was isolated from peripheral blood samples of 69 lung cancer patients. Reverse transcriptase polymerase chain reaction (RT-PCR) was used to test for the presence of cytokeratin 19 and the 4 gene mRNA transcripts in blood RNA. The positive detection rate of at least 1 of the 5 transcripts was 79% for lung adenocarcinoma and 62% for squamous carcinoma. Using RT-PCR, at least 1 of the markers was found in 53% of stage I patients, 100% of stage II, 71% of stage III and 81% of stage IV lung cancer patients. Blood samples from 20 healthy volunteers were also tested, but only 1 of the 5 transcripts was found in 1 patient. These new molecular markers may aid early detection, staging and follow-up of lung cancer patients by RNA isolated from blood.

We have developed a reliable method for eliminating base-mispair amplification in selective polymerase chain reaction (PCR), which is utilized for amplifying unknown sequence fragments produced by restriction enzyme reaction. The proposed procedure applies amplified fragment length polymorphism (AFLP) with high fidelity. Selective PCR utilizes the known polymerase reaction characteristic that the complementary strand extension is strongly affected by matching a template with the 3'-terminus of the primers. However, false positive amplification is frequently observed because the specificity of terminal bases for discrimination of fragments (usually, 1-3 anchor sequences) is not enough to separate each fragment. A protocol for the selective PCR separation of every fragment was therefore investigated. A single-base mismatch was artificially introduced on the 4th base position from the 3' end of the primers to improve the hybridization specificity of anchored 2-bases at the 3' termini of primers. PCR reaction was carried out at 66 degrees C to prevent false positive amplification. The concentration of the primers having anchored-base sequences of AA, AT, TA, and TT must be three times larger than that of other primers because the T-m values for these sequences are lower than the others. As all the fragments can be separated into groups with high fidelity, the improved selective PCR will be applied to gene finding and analyzing differences on genome sequences based on AFLP.

A new method for selecting and amplifying a single DNA fragment from a mixture is proposed. This method is applicable for the rapid classification of DNA fragments from a mixture and for preparation of sequencing templates. DNAs of several to tens of kilobases (kb) are digested with a four-base recognition restriction enzyme to produce smaller fragments. The complementary strand extension reactions are then carried out to produce fluorophore-labeled DNA fragments from the digestion products. These fragments can be rapidly classified according to their terminal-base sequences and their sizes are analyzed by capillary-array gel electrophoresis (CAGE). Electropherograms are used to characterize the fragments and to select polymerase chain reaction (PCR) primers. Any fragment in a digestion mixture can be amplified by PCR with a pair of primers selected from a primer pool by referring to the electropherograms of the fragments. This method was successfully used to compare the electropherograms of two different DNA strands and to sequence a several-kb DNA fragment without subcloning. Combined with CAGE, this method could be used to dramatically simplify DNA fragment analysis.