Non-biased and efficient global amplification of a single-cell cDNA library.

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.

Position-Specific Gene Expression Analysis Using a Microgram Dissection Method Combined with On-Bead cDNA Library Construction.

Gene expression analysis is a key technology that is used to understand living systems. Multicellular organisms, including plants, are composed of various tissues and cell types, each of which exhibits a unique gene expression pattern. However, because of their rigid cell walls, plant cells are difficult to isolate from the whole plant. Although laser dissection has been used to circumvent this problem, the plant sample needs to be fixed beforehand, which presents several problems. In the present study, we developed an alternative method to conduct highly reliable gene expression profiling. First, we assembled a dissection apparatus that used a narrow, sharpened needle to dissect out a microsample of fresh plant tissue (0.1-0.2 mm on each side) automatically from a target site within a short time frame. Then, we optimized a protocol to synthesize a high-quality cDNA library on magnetic beads using a single microsample. The cDNA library was amplified and subjected to high-throughput sequencing. In this way, a stable and reliable system was developed to conduct gene expression profiling in small regions of a plant. The system was used to analyze the gene expression patterns at successive 50 µm intervals in the shoot apex of a 4-day-old Arabidopsis seedling. Clustering analysis of the data demonstrated that two small, adjacent domains, the shoot apical meristem and the leaf primordia, were clearly distinguishable. This system should be broadly applicable in the investigation of the spatial organization of gene expression in various contexts.

Spatial Regulation of the Gene Expression Response to Shade in Arabidopsis Seedlings.

The shade avoidance response, which allows plants to escape from nearby competitors, is triggered by a reduction in the PFR form of phytochrome in response to shade. Classic physiological experiments have demonstrated that the shade signal perceived by the leaves is transmitted to the other parts of the plant. Recently, a simple method was developed to analyze the transcriptome in a single microgram tissue sample. In the present study, we adopted this method to conduct organ-specific transcriptomic analysis of the shade avoidance response in Arabidopsis seedlings. The shoot apical samples, which contained the meristem, basal parts of leaf primordia and short fragments of vasculature, were collected from the topmost part of the hypocotyl and subjected to RNA sequencing analysis. Unexpectedly, many more genes were up-regulated in the shoot apical region than in the cotyledons. Spotlight irradiation demonstrated that the apex-responsive genes were mainly controlled by phytochrome in the cotyledons. In accordance with the involvement of many auxin-responsive genes in this category, auxin biosynthesis was genetically shown to be essential for this response. In contrast, organ-autonomous regulation was more important for the genes that were up-regulated preferentially either in the cotyledons or in both the cotyledons and the apical region. Their responses to shade depended variously on auxin and PIFs (phytochrome-interacting factors), indicating the mechanistic diversity of the organ-autonomous response. Finally, we examined the expression of the auxin synthesis genes, the YUC genes, and found that three YUC genes, which were differently spatially regulated, co-ordinately elevated the auxin level within the shoot apical region.

Involvement of Ca2+ in Vacuole Degradation Caused by a Rapid Temperature Decrease in Saintpaulia Palisade Cells: A Case of Gene Expression Analysis in a Specialized Small Tissue.

Saintpaulia (African violet) leaves are known to be damaged by a rapid temperature decrease when cold water is applied to the leaf surface; the injury is ascribed to the chloroplast damage caused by the cytosolic pH decrease following the degradation of the vacuolar membrane in the palisade cells. In this report, we present evidence for the involvement of Ca(2+) in facilitating the collapse of the vacuolar membrane and in turn in the temperature sensitivity of Saintpaulia leaves. In the presence of a Ca(2+) chelator (EGTA) or certain Ca(2+) channel inhibitors (Gd(3+) or La(3+)) but not others (verapamil or nifedipine), the pH of the vacuole, monitored through BCECF (2',7'-bis(carboxyethyl)-4 or 5-carboxyfluorescein) fluorescence, did not increase in response to a rapid temperature drop. These pharmacological observations are consistent with the involvement of mechanosensitive Ca(2+) channels in the collapse of the vacuolar membrane. The high level of expression of an MCA- (Arabidopsis mechanosensitive Ca(2+) channel) like gene, a likely candidate for a mechanosensitive Ca(2+) channel(s) in plant cells, was confirmed in the palisade tissue in Saintpaulia leaves by using a newly developed method of gene expression analysis for the specialized small tissues.

A highly sensitive and accurate gene expression analysis by sequencing ("bead-seq") for a single cell.

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.

Self-oriented immobilization of DNA polymerase tagged by titanium-binding peptide motif.

We developed a titanium-binding-peptide-1 (TBP-1)-tagged DNA polymerase, for self-oriented immobilization onto a titanium oxide (TiO2) substrate. The enzymatic function of a polymerase immobilized on a solid state device is strongly dependent on the orientation of the enzyme. The TBP-tagged DNA polymerase, which was derived from a hyperthermophilic archaeon, was designed to incorporate the RKLPDA peptide at the N-terminus, and synthesized by translation processes in Escherichia coli (E. coli). The specific binding of the TBP-tagged DNA polymerase onto a TiO2 substrate was clearly monitored by surface plasmon resonance spectroscopy (SPR) and by surface potential detection with an extended-gate field effect transistor (FET). In the SPR analyses, constant quantities of the DNA polymerase were stably immobilized on the titanium substrate under flow conditions, regardless of the concentration of the DNA polymerase, and could be completely removed by a 4 M MgCl2 wash after measurement. The FET signal showed the contribution of the molecular charge in the TBP motif to the binding with TiO2. In addition, the TBP-tagged DNA polymerase-tethered TiO2 gate electrode enabled the effective detection of the positive charges of hydrogen ions produced by the DNA extension reaction, according to the FET principle. Therefore, the self-oriented immobilization platform based on the motif-inserted enzyme is suitable for the quick and stable immobilization of functional enzymes on biosensing devices.

Hedgehog-Gli activators direct osteo-chondrogenic function of bone morphogenetic protein toward osteogenesis in the perichondrium.

Specification of progenitors into the osteoblast lineage is an essential event for skeletogenesis. During endochondral ossification, cells in the perichondrium give rise to osteoblast precursors. Hedgehog (Hh) and bone morphogenetic protein (BMP) are suggested to regulate the commitment of these cells. However, properties of perichondrial cells and regulatory mechanisms of the specification process are still poorly understood. Here, we investigated the machineries by combining a novel organ culture system and single-cell expression analysis with mouse genetics and biochemical analyses. In a metatarsal organ culture reproducing bone collar formation, activation of BMP signaling enhanced the bone collar formation cooperatively with Hh input, whereas the signaling induced ectopic chondrocyte formation in the perichondrium without Hh input. Similar phenotypes were also observed in compound mutant mice, where signaling activities of Hh and BMP were genetically manipulated. Single-cell quantitative RT-PCR analyses showed heterogeneity of perichondrial cells in terms of natural characteristics and responsiveness to Hh input. In vitro analyses revealed that Hh signaling suppressed BMP-induced chondrogenic differentiation; Gli1 inhibited the expression of Sox5, Sox6, and Sox9 (SRY box-containing gene 9) as well as transactivation by Sox9. Indeed, ectopic expression of chondrocyte maker genes were observed in the perichondrium of metatarsals in Gli1(-/-) fetuses, and the phenotype was more severe in Gli1(-/-);Gli2(-/-) newborns. These data suggest that Hh-Gli activators alter the function of BMP to specify perichondrial cells into osteoblasts; the timing of Hh input and its target populations are critical for BMP function.

Improvement of LATE-PCR to allow single-cell analysis by pyrosequencing.

Nucleic acid analysis in a single cell is very important, but the extremely small amount of template in a single cell requires a detection method more sensitive than the conventional method. In this paper, we describe a novel assay allowing a single cell genotyping by coupling improved linear-after-the-exponential-PCR (imLATE-PCR) on a modified glass slide with highly sensitive pyrosequencing. Due to the significantly increased yield of ssDNA in imLATE-PCR amplicons, it is possible to employ pyrosequencing to sequence the products from 1 µL chip PCR which directly used a single cell as the starting material. As a proof-of-concept, the 1555A>G mutation (related to inherited deafness) on mitochondrial DNA and the SNP 2731C>T of the BRCA1 gene on genomic DNA from a single cell were successfully detected, indicating that our single-cell-pyrosequencing method has high sensitivity, simple operation and is low cost. The approach has promise to be of efficient usage in the fields of diagnosis of genetic disease from a single cell, for example, preimplantation genetic diagnosis (PGD).

Digital detection of multiple minority mutants in stool DNA for noninvasive colorectal cancer diagnosis.

Somatic mutations in stool DNA are quite specific to colorectal cancer (CRC), but a method being able to detect the extraordinarily low amounts of mutants is challengeable in sensitivity. We proposed a hydrogel bead-array to digitally count CRC-specific mutants in stool at a low cost. At first, multiplex amplification of targets containing multiple mutation loci of interest is carried out by a target enriched multiplex PCR (Tem-PCR), yielding the templates qualified for emulsion PCR (emPCR). Then, after immobilizing the beads from emPCR on a glass surface, the incorporation of Cy3-dUTP into the mutant-specific probes, which are specifically hybridized with the amplified beads from emPCR, is used to color the beads coated with mutants. As all amplified beads are hybridized with the Cy5-labeled universal probe, a mutation rate is readily obtained by digitally counting the beads with different colors (yellow and red). A high specificity of the method is achieved by removing the mismatched probes in a bead-array with electrophoresis. The approach has been used to simultaneously detect 8 mutation loci within the APC, TP53, and KRAS genes in stools from eight CRC patients, and 50% of CRC patients were positively diagnosed; therefore, our method can be a potential tool for the noninvasive diagnosis of CRC.

Multiplex loop-mediated isothermal amplification detection by sequence-based barcodes coupled with nicking endonuclease-mediated pyrosequencing.

The loop-mediated isothermal amplification (LAMP) is a well-developed method for replicating a targeted DNA sequence with a high specificity, but multiplex LAMP detection is difficult because LAMP amplicons are very complicated in structure. To allow simultaneous detection of multiple LAMP products, a series of target-specific barcodes were designed and tagged in LAMP amplicons by FIP primers. The targeted barcodes were decoded by pyrosequencing on nicked LAMP amplicons. To enable the nicking reaction to occur just near the barcode regions, the recognition sequence of the nicking endonuclease (NEase) was also introduced into the FIP primer. After the nicking reaction, pyrosequencing started at the nicked 3' end when the added deoxyribonucleoside triphosphate (dNTP) was complementary to the non-nicked strand. To efficiently encode multiple targets, the barcodes were designed with a reporter base and two stuffer bases, so that the decoding of a target-specific barcode only required a single peak in a pyrogram. We have successfully detected the four kinds of pathogens including hepatitis B virus (HBV), hepatitis C virus (HCV), human immunodeficiency virus (HIV), and Treponema pallidum (TP), which are easily infected in blood, by a 4-plex LAMP in a single tube, indicating that barcoded LAMP coupled with NEase-mediated pyrosequencing is a simple, rapid, and reliable way in multiple target identification.

Quantitative analysis of gene expression in a single cell by qPCR.

We developed a quantitative PCR method featuring a reusable single-cell cDNA library immobilized on beads for measuring the expression of multiple genes in a single cell. We used this method to analyze multiple cDNA targets (from several copies to several hundred thousand copies) with an experimental error of 15.9% or less. This method is sufficiently accurate to investigate the heterogeneity of single cells.

Integrated spatial analysis of gene mutation and gene expression for understanding tumor diversity in formalin-fixed paraffin-embedded lung adenocarcinoma.

Introduction: A deeper understanding of intratumoral heterogeneity is essential for prognosis prediction or accurate treatment plan decisions in clinical practice. However, due to the cross-links and degradation of biomolecules within formalin-fixed paraffin-embedded (FFPE) specimens, it is challenging to analyze them. In this study, we aimed to optimize the simultaneous extraction of mRNA and DNA from microdissected FFPE tissues (φ = 100 µm) and apply the method to analyze tumor diversity in lung adenocarcinoma before and after erlotinib administration. Method: Two magnetic beads were used for the simultaneous extraction of mRNA and DNA. The decross-linking conditions were evaluated for gene mutation and gene expression analyses of microdissected FFPE tissues. Lung lymph nodes before treatment and lung adenocarcinoma after erlotinib administration were collected from the same patient and were preserved as FFPE specimens for 4 years. Gene expression and gene mutations between histologically classified regions of lung adenocarcinoma (pre-treatment tumor in lung lymph node biopsies and post-treatment tumor, normal lung, tumor stroma, and remission stroma, in resected lung tissue) were compared in a microdissection-based approach. Results: Using the optimized simultaneous extraction of DNA and mRNA and whole-genome amplification, we detected approximately 4,000-10,000 expressed genes and the epidermal growth factor receptor (EGFR) driver gene mutations from microdissected FFPE tissues. We found the differences in the highly expressed cancer-associated genes and the positive rate of EGFR exon 19 deletions among the tumor before and after treatment and tumor stroma, even though they were collected from tumors of the same patient or close regions of the same specimen. Conclusion: Our integrated spatial analysis method would be applied to various FFPE pathology specimens providing area-specific gene expression and gene mutation information.

Reproducible and sensitive micro-tissue RNA sequencing from formalin-fixed paraffin-embedded tissues for spatial gene expression analysis.

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 two years after fixation and embedding and detected approximately 7000 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.

Highly sensitive pyrosequencing based on the capture of free adenosine 5' phosphosulfate with adenosine triphosphate sulfurylase.

In pyrosequencing chemistry, four cascade enzymatic reactions with the catalysis of polymerase, adenosine triphosphate (ATP) sulfurylase, luciferase, and apyrase are employed. The sensitivity of pyrosequencing mainly depends on the concentration of luciferase which catalyzes a photoemission reaction. However, the side-reaction of adenosine 5' phosphosulfate (APS, an analogue of ATP) with luciferase resulted in an unavoidable background signal; hence, the sensitivity cannot be much higher due to the simultaneous increase of the background signal when a larger amount of luciferase is used. In this study, we demonstrated a sensitive pyrosequencing using a large amount of ATP sulfurylase to lower the concentration of free APS in the pyrosequencing mixture. As the complex of ATP sulfurylase and APS does not react with luciferase, a large amount of luciferase can be used to achieve a sensitive pyrosequencing reaction. This sensitivity-improving pyrosequencing chemistry allows the use of an inexpensive light sensor photodiode array for constructing a portable pyrosequencer, a potential tool in a point-of-care test (POCT).

Highly sensitive pyrosequencing system with polymer-supported enzymes for high-throughput DNA analysis.

A highly sensitive massively parallel pyrosequencing system employing a gel matrix to immobilize enzymes at high density in microreaction chambers is demonstrated. Reducing the size of microreaction chambers in a DNA analyzer is important to achieve a high throughput utilizing a commercially available detection device or camera. A high-performance system can be attained by detecting signals from one reaction chamber with one photopixel of around several micrometers by utilizing a 1:1 image magnification. However, the use of small beads immobilizing DNA has a disadvantage in detecting luminescence because only small amounts of DNA can be immobilized on the bead surfaces for sequencing. As luminescence intensity could be enhanced by increasing the luciferase density in the chambers, we overcame this difficulty by using a gel matrix to immobilize luciferase at a high concentration in the microreaction chambers. Luminescence 1 order of magnitude higher could be observed with the new method compared to the conventional method. Consequently, the chamber size and bead size immobilizing DNA could be reduced to as small as 6.5 and 4 µm, respectively. This can be successfully applied to achieving small, inexpensive, pyrosequencing systems with high throughput.

Optimization of pyrosequencing reads by superior successive incorporation efficiency of improved 2'-deoxyadenosine-5'-triphosphate analogs.

Conventional pyrosequencing using 2'-deoxyadenosine-5'-O-(1-thiotriphosphate) (dATPαS) is problematic due to the high cost of the substrate (dATPαS) and deterioration in the accuracy of incorporation to read through poly(T) regions. One reason for these problems is that dATPαS has a sulfur on the α-phosphate and also has isomers (Sp and Rp). To solve these problems, 11 nucleotide substrates, which could replace dATPαS in pyrosequencing, were newly synthesized. All substrates were modified on the seventh or eighth position of the adenine base from normal dATP. We found that the substrate that had an ethenyl-linked modified group on the seventh position of the adenine base had low activity in the luciferase reaction and high incorporation efficiency with the thymine base. One substrate in particular had 10-fold better incorporation efficiency than dATPαS. The new nucleotide substrate satisfied all conditions as a replacement of dATPαS.

Digital analysis of the expression levels of multiple colorectal cancer-related genes by multiplexed digital-PCR coupled with hydrogel bead-array.

To digitally analyze expression levels of multiple genes in one reaction, we proposed a method termed as 'MDHB' (Multiplexed Digital-PCR coupled with Hydrogel Bead-array). The template for bead-based emulsion PCR (emPCR) was prepared by reverse transcription using sequence-tagged primers. The beads recovered from emPCR were immobilized with hydrogel to form a single-bead layer on a chip, and then decoded by gene-specific probe hybridization and Cy3-dUTP based primer extension reaction. The specificity of probe hybridization was improved by using electrophoresis to remove mismatched probes on the bead's surface. The number of positive beads reflects the abundance of expressed genes; the expression levels of target genes were normalized to a housekeeping gene and expressed as the number ratio of green beads to red beads. The discrimination limit of MDHB is 0.1% (i.e., one target molecule from 1000 background molecules), and the sensitivity of the method is below 100 cells when using the ß-actin gene as the detection target. We have successfully employed MDHB to detect the relative expression levels of four colorectal cancer (CRC)-related genes (c-myc, COX-2, MMP7, and DPEP1) in 8 tissue samples and 9 stool samples from CRC patients, giving the detection rates of 100% and 77%, respectively. The results suggest that MDHB could be a potential tool for early non-invasive diagnosis of CRC.

Pyrosequencing on nicked dsDNA generated by nicking endonucleases.

Although the pyrosequencing method is simple and fast, the step of ssDNA preparation increases the cost, labor, and cross-contamination risk. In this paper, we proposed a method enabling pyrosequencing directly on dsDNA digested by nicking endonucleases (NEases). Recognition sequence of NEases was introduced using artificially mismatched bases in a PCR primer (in the case of genotyping) or a reverse-transcription primer (in the case of gene expression analysis). PCR products were treated to remove excess amounts of primers, nucleotides, and pyrophosphate (PPi) prior to sequencing. After the nicking reaction, pyrosequencing starts at the nicked 3' end, and extension reaction occurs when the added dNTP is complementary to the non-nicked strand. Although the activity of strand displacement by Klenow is limited, approximately 10 bases are accurately sequenced; this length is long enough for genotyping and SRPP-based differential gene expression analysis. It was observed that the signals of two allele-specific bases in a pyrogram from nicked dsDNA are highly quantitative, enabling quantitative determination of allele-specific templates; thus, Down's Syndrome diagnosis as well as differential gene expression analysis was successfully executed. The results indicate that pyrosequencing using nicked dsDNA as templates is a simple, inexpensive, and reliable way in either quantitative genotyping or gene expression analysis.

Gene expression analysis on a photodiode array-based bioluminescence analyzer by using sensitivity-improved SRPP.

Most methods used for gene expression analysis are based on dye-labeling, which requires costly instruments. Recently a dye-free gene expression analysis method-SRPP (Sequence-tagged reverse-transcription polymerase chain reaction coupled with pyrosequencing) was developed to compare relative gene expression levels in different tissues, but the throughput of the SRPP assay is very limited due to the use of a photomultiplier tube (PMT)-based pyrosequencer for the detection. To increase the throughput of the SRPP assay, an inexpensive photodiode (PD) array-based bioluminescence analyzer (termed as "PD-based pyrosequencer") was coupled to SRPP; however the low sensitivity of PD limited the wide application of SRPP. To enable SRPP analyzing low abundance genes in clinical samples, sequence-tagged gene-specific primers instead of sequence-tagged poly (T)(n) primers were used for reverse-transcription, and the SRPP sensitivity was thus improved more than 10 times. This improvement compensates the sensitivity loss due to the use of PD in a pyrosequencer. The accurate determination of the expression levels of ten prognostic marker genes (AL080059, MMP9, EXT1, ORC6L, AF052162, C9orf30, FBXO31, IGFBP5, ESM1, and RUNDC1) differing between normal tissues and tumor tissues of breast cancer patients demonstrated that SRPP using gene-specific RT primers coupled with the PD array-based bioluminescence analyzer is reliable, inexpensive, and sensitive in gene expression analysis.

Slow-Cycling Cancer Stem Cells Regulate Progression and Chemoresistance in Colon Cancer.

Cancer chemoresistance is often attributed to the presence of cancer stem cell (CSC)-like cells, but whether they are homogeneously chemoresistant remains unclear. We previously showed that in colon tumors, a subpopulation of LGR5+ CSC-like cells driven by TCF1 (TCF7), a Wnt-responsive transcription factor, were responsible for tumorigenicity. Here we demonstrate that the tumorigenic subpopulation of mouse LGR5+ cells exists in a slow-cycling state and identify a unique 22-gene signature that characterizes these slow-cycling CSC. Seven of the signature genes are specifically expressed in slow-cycling LGR5+ cells from xenografted human colon tumors and are upregulated in colon cancer clinical specimens. Among these seven, four genes (APCDD1, NOTUM, PROX1, and SP5) are known to be direct Wnt target genes, and PROX1 was expressed in the invasive fronts of colon tumors. PROX1 was activated by TCF1 to induce CDKN1C and maintain a slow-cycling state in colon cancer organoids. Strikingly, PROX1 was required for recurrent growth after chemotherapeutic treatment, suggesting that inhibition of slow-cycling CSC by targeting the TCF1-PROX1-CDKN1C pathway is an effective strategy to combat refractory colon cancer in combination with conventional chemotherapy. SIGNIFICANCE: These findings illustrate the importance of a slow-cycling CSC subpopulation in colon cancer development and chemoresistance, with potential implications for the identified slow-cycling CSC signatures and the TCF1-PROX1-CDKN1C pathway as therapeutic targets.