CHARACTERISTICS AND GENE EXPRESSION OF SINGLE-CELL

CROSS-REFERENCE OF RELATED APPLICATION [0001] The present application claims the priority of U.S. Provisional Patent Application No. 63/197331 filed on June 5, 2021, which is incorporated herein by reference in Patent Application.

FIELD

[0002] The present application relates to biotechnologies, and more particularly, to a method and a system for integrating morphological characteristics and gene expression of individual cells.

BACKGROUND

[0003] Single-cell RNA sequencing (scRNA-seq) technologies have revolutionized the way for transcriptomic analysis of multicellular tissues. Via investigation of thousands of cells at single-particle resolution, scRNA-seq can provide quantitative expression profiles of individual cells with valuable insights into cellular differences, such as cell types and cell states, which are usually elusive in the traditional bulk RNA-seq analysis. The distinct advantages of single-cell transcriptomics enable multi-dimensional investigation of individual cells to, for example, decipher tumor heterogeneity, reveal complex and rare cell populations, and uncover regulatory relationships between genes, offering a strong basis for designing precision medicine and targeted therapy.

[0004] On the other hand, due to the rapid evolution of computational technologies (e.g., deep learning and neural networks), cellular morphological profiling based on high-throughput and high-content image processing has become an emerging tool for various biological and medical applications. For example, imaging flow cytometry (IFC), the representative technology in this field, can extract hundreds of morphological features (e.g., shape, size, intensity, and texture) from each individual cell with AI-driven algorithms, allowing classification of complex cell phenotypes, identification of rare cells, and discovery of useful targets for disease diagnosis, personalized medicine, and drug development.

[0005] Since the morphological characteristics are highly correlated with gene expression patterns, researchers have shown substantial interests in the integration of phenomics and transcriptomics data for novel biological insights. However, currently there is no such an approach that could link gene expression profiles to morphological phenotypes at single-cell level with a high-throughput manner. Some studies are carried out based on bulk RNA-seq analysis while ignore cellular heterogeneity. The others rely on manual collection of each target cell using pipette, which have a limited throughput.

SUMMARY

[0006] In view of this, the present application provides a method and a system for integrating morphological characteristics and gene expression of individual cells, so as to link gene expression information with morphological characteristics at single-cell level in a high-throughput and high-efficiency manner.

[0007] An embodiment of the present application provides a method for integrating the morphological characteristics and gene expression of individual cells, which includes the following steps:

[0008] provide a microfluidic device, the microfluidic device includes a microwell array composed of a plurality of microwells and an interdigital electrode, each microwell includes a plurality of capture oligonucleotides of a certain sequence, and each capture oligonucleotide includes a known cell barcode sequence, a unique molecular identifier sequence, and a capture sequence, the cell barcode sequence has a unique corresponding relationship with the microwell in which the cell barcode sequence is located, the cell barcode is used to mark a cell in the microwell, the unique molecular identifier sequence is used to mark mRNA captured by the capture oligonucleotide;

[0009] inject cells into the micro wells, and use the interdigital electrode to capture a single cell above each of the microwells, and record the morphological characteristics of the cell in the micro well;

[0010] lyse the cell so that the mRNA released by the cell is captured by the capture oligonucleotide in the microwell where the cell is located;

[0011] reverse transcribe the captured mRNA to obtain cDNA, each cDNA includes a capture oligonucleotide and a nucleotide sequence complementary to the captured mRNA; [0012] perform a PCR amplification reaction on the cDNA to obtain a cDNA library, and sequence the cDNA library;

[0013] read the cell barcode sequence and the unique molecular identifier sequence according to a sequencing result, and the morphological characteristics and gene expression information of the cell in the microwell are integrated together.

[0014] In one embodiment, the cell barcode sequence includes 6 to 10 bases.

[0015] In one embodiment, the unique molecular identifier sequence includes 7 to 15 bases.

[0016] In one embodiment, the capture oligonucleotide further includes a PCR handle sequence, and the PCR handle sequence includes 30 to 35 bases.

[0017] In one embodiment, the morphological characteristics of each cell include bright-field images and fluorescent images.

[0018] The present application also provides a system for integrating the morphological characteristics and gene expression of single-cell. The system includes a microfluidic device, an imaging device, a PCR machine, a sequencer, and a controller. The microfluidic device is used for dispersing a cell population into individual cells and marking the individual cells respectively, and the microfluidic device includes a microwell array composed of a plurality of microwells and an interdigital electrode. Each microwell includes a capture oligonucleotide of a known sequence, and each capture oligonucleotide includes a cell barcode sequence, a unique molecular identifier sequence, and a capture sequence. The cell barcode sequence has a unique corresponding relationship with the microwell in which it is located, the cell barcode sequence is used to mark the cells in the microwells, and the unique molecular identifier sequence is used to mark the captured mRNA. The imaging device is used to capture images the cells and record the morphological characteristics of the cells for morphological analysis. The PCR machine is used for a PCR amplification reaction on cDNA to obtain a cDNA library. The sequencer is used to sequence the cDNA library. The controller is used to integrate the morphological characteristics and gene expression of the cells.

[0019] In one embodiment, the diameter of the microwell is 20 pm, the depth of the micro well is 20 pm, and the distance between two adjacent micro wells is 85 pm.

[0020] In one embodiment, the number of the microwells is 10,000 to 100,000.

[0021] In one embodiment, the interdigital electrode is arranged beneath the microwell array.

[0022] In one embodiment, the imaging device includes a CCD camera and a microscope connected to the CCD camera.

[0023] In this application, a single cell is placed in a single microwell in the microfluidic device, and from the beginning (imaging) to the end (sequencing), each cell is assigned a unique known cell barcode sequence to observe their phenotype before processing for sequencing. The cell barcode sequence in the capture oligonucleotide can be "read" in the micro wells and also can be "read" from the sequence reads obtained from the cDNA library, thereby the genome/transcriptome data (mRNA sequence information) is linked to the observed phenotype of single-cell, so that the morphological phenotype is directly related to gene expression. This method focuses on integrating the morphological characteristics and gene expression profiles of isolated single cells. It has the characteristics of high efficiency, single-cell resolution, high recovery rate and sensitivity, and low RNA contamination. It can facilitate fundamental biological studies, develop multi-dimensional biomarker signatures for diseases, and accelerate drug discovery and development.

BRIEF DESCRIPTION OF THE DRAWINGS [0024] Implementations of the present disclosure will now be described, by way of embodiments only, with reference to the attached figures.

[0025] FIG. 1 is a diagrammatic view of a microfluidic device according to an embodiment of the present application.

[0026] FIG. 2 is a cross-sectional view taken along A-A in FIG. 1.

[0027] FIG. 3 is a diagrammatic view of a capture oligonucleotide according to an embodiment of the present application.

[0028] FIG. 4 is a diagrammatic view of a process of synthesizing capture oligonucleotides via inkjet printing.

[0029] FIG. 5 is a diagrammatic view of the capture oligonucleotide synthesized in FIG. 4.

[0030] FIG. 6 is a diagrammatic view of a process for performing transcription analysis on cells according to an embodiment of the application.

[0031] FIG. 7 is a diagrammatic view illustrating morphological features linking to gene expression profiles at single-cell level with a high throughput. [0032] FIG. 8 is a diagram of a system for integrating morphological characteristics and gene expression of a single cell according to an embodiment of the application.

[0033] FIG. 9 is an image of the cells separated from the microwells in Example 1 of the present application.

[0034] FIG. 10 is a diagram of cDNA library analysis via a bioanalyzer in Example 1 of the present application.

[0035] FIG. 11 is a diagram of the analysis of mouse cells and human cells in Example 1 of the present application.

Symbol description of main components microfluidic device 10 micro well 101 interdigital electrode 103

SU-8 photoresist layer 105 glass substrate 107 glass channel 109 inlet 102 outlet 104 capture oligonucleotide 1010 cell barcode sequence 1011 unique molecular identifier 1012 capture sequence 1013

PCR handle sequence 1014 cell 20 system 200 imaging device 30

PCR machine 50 sequencer 70 controller 90

[0036] The following detailed description will further illustrate the present application in conjunction with the above-mentioned drawings.

DETAILED DESCRIPTION

[0037] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by those skilled persons in the art. The terms used herein is only for the purpose of describing specific implementation manners, and is not intended to limit the embodiments of the present application.

[0038] In the case of no conflict, the following embodiments and features in the embodiments can be combined with each other.

[0039] An embodiment of the present application provides a method for integrating morphological characteristics and gene expression of individual cells, which comprises the following steps:

[0040] SI: a microfluidic device 10 is provided.

[0041] Referring to Fig. 1, in some embodiments, the microfluidic device 10 comprises a micro well array composed of a plurality of micro wells 101 and an interdigital electrode 103. Referring to Fig. 2, in some embodiments, the microwell array can be fabricated on a SU-8 photoresist layer 105 by photolithography, and the SU-8 photoresist layer 105 is disposed on a glass substrate 107. In order to improve cell capturing efficiency, an interdigital electrode (IDE) 103 is patterned on the glass substrate 107 and under the micro well array, that is, the interdigital electrode 103 is arranged between the SU-8 photoresist layer 105 and the glass substrate 107. As shown in Fig. 1, the glass channel 109 is bonded onto the top of the SU-8 photoresist layerl05 via an adhesive to form a sealed flow cell with a single inlet 102 and a single outlet 104. During the injection of a cell sample, a dielectrophoresis force (DEP force) generated by the interdigital electrode 103 will trap cells and guide them into the micro wells 101.

[0042] Referring to Fig. 2 and Fig. 3, in some embodiments, each microwell 101 comprises a plurality of capture oligonucleotides 1010, and each capture oligonucleotide 1010 comprises a known cell barcode sequence 1011, a unique molecule identifier (UMI) sequence 1012, and a capture sequence 1013. Furthermore, the capture oligonucleotide 1010 also comprises a PCR handle sequence 1014. There is a unique correspondence between the cell barcode sequence 1011 and the microwell 101 where the cell barcode sequence 1011 is located, and the cell barcode sequence 1011 is used to mark the micro well 101 where the cell barcode sequence 1011 is located, so as to mark the cell in the microwell 1011. The unique molecular identifier sequence 1012 is used to label the captured mRNA to avoid repeated counting after PCR. In some embodiments, the capture sequence 1013 may be a Poly dT sequence with a size of 30 bases, which is used to hybridize with the Poly dA segment of the mRNA released after cell lysis to capture the mRNA. In some embodiments, the cell barcode sequence 1011 comprises 6 to 10 bases. In some embodiments, the unique molecular identifier sequence 1012 comprises 7 to 15 bases. In some embodiments, the PCR handle sequence 1014 comprises 30 to 35 bases. Furthermore, the PCR handle sequence 1014 may include 32 bases.

[0043] In some embodiments, the capture oligonucleotide 1010 can be synthesized by an inkjet printing technology, and the printing process is shown in Fig. 4. The epoxide groups on the SU-8 surface are available for direct conjugation of amine-modified oligonucleotides. First, the PCR handle sequence 1014 with 3’-dimethoxytrityl (DMT) protection is printed to the substrate, with the amine-modified 5’ end being linked to the surface of the SU-8 photoresist layer 105. Then, in order to generate the cell barcode sequence 1011, DMTs are selectively removed by chemicals, followed by a coupling reaction with specific nucleotides which only occurs on the spots that are de-protected. By repeating this process, a pre-designed 7-mers unique cell barcode sequence 1011 is developed for each microwell. Next, in order to generate a unique molecular identifier (UMI) sequence, a mixture of 4 types of DMT-protected nucleotides (A, G, C, T) is repeatedly added to all spots to generate 10-mers random unique sequences within each spot. The function of UMIs is to assign each captured mRNA transcript an identifier to avoid repeated counting after PCR. Finally, a poly dT tail is linked to each oligonucleotide, and the capture oligonucleotide 1010 as shown in FIG. 5 is obtained. The cell barcodes (7-mers) and UMIs (10-mers) have a maximum of 16,384 (4 ⁷ ) and 1,048,576 (4 ¹⁰ ) different sequences, that is, up to 16,384 different cells can be distinguished at a time. For samples with more than 16384 cells, the length of cell barcode sequence can be further increased.

[0044] To validate that the cell barcode sequence 1011 are well synthesized, 100 spots with known oligonucleotide sequence are distributed at different locations of the microwell array for hybridization with fluorescent probes. The hybridization rate of the 100 spots and the CV (coefficient of variation) of fluorescent signal intensity are characterized by fluorescence imaging. In the present application, the probe hybridization rate is greater than 95%, and CV of fluorescent signal intensity is less thanl0%, indicating that the cell barcode sequence 1011 synthesized by this application has a high accuracy rate.

[0045] S2: cells are injected into the microwells 101, and the interdigital electrode

103 is used to capture a single cell in the microwell 101. The morphological characteristics of the cells in the microwells 101 are recorded for morphological analysis. [0046] During operations, cells are injected through the inlet 102, which flow into the microwell array through the glass channel 109. The dielectrophoretic (DEP) force generated by the interdigital electrode 103 can trap single cells above the microwells. After the cells stop flowing, the interdigital electrode 103 is turned off, allowing the cells to be precipitated into the microwells. The excess cells outside the microwells are washed away from the outlet 104. The ideal situation is that every microwell traps only one cell. In order to achieve the ideal situation, the key parameters, including but not limited to microwell dimension (diameter and depth), DEP force intensity, channel height, input cell concentration, and flow rate, are optimized to maximize the single-cell purity and microwell occupancy rate. In one embodiment, the diameter of the microwell is 25 mih, the depth of the microwell is 20 pm, the distance between two adjacent microwells (the distance from the center to the center) is 85 pm, and the number of microwells is 10,000 to 100,000.

[0047] After individual cells are trapped in the microwell array, the bright-field and fluorescent images of each cell for morphological profiling are recorded via a CCD (Charge Coupled Device) camera connected to a microscope. Meanwhile, the cell barcodes sequence 1011 on the capture oligonucleotide 1010 are assigned to the cells based on their locations in the array. For example, the cell barcode sequence of the microwell in the first row and the first column is known as TACGAGC (TACGAGC is unique among all cell barcode sequences), and TACGAGC is assigned to the cell located in the first row and the first column of the microwells.

[0048] Air bubbles in the flow cell will cause failure to reverser transcription and PCR due to the bubble expansion at high temperature. To remove the air bubbles, the glass channel 107 is treated with 80% ethanol before buffer injection. Meanwhile, by applying a voltage to the interdigital electrode 103, the electro wetting behavior makes the surface hydrophilic, which is helpful to remove the bubbles trapped in the microwells. [0049] S3: cells are lysed so that the mRNA released by the cell is captured by the capture oligonucleotide in the microwell where the cell is located.

[0050] Referring to Fig. 6, after cell imaging, in-site cell lysis is carried out via the injection of lysis buffer. The mRNA (with a poly dA tail) released by cells will be captured by the nearest capture oligonucleotide 1010 due to hybridization between poly dAand poly dT tails (i.e., capture sequence 1013 in the capture oligonucleotide 1010). [0051] S4: the captured mRNA is reverse transcribed to obtain cDNA, each cDNA comprises a capture oligonucleotide sequence 1010 and a nucleotide sequence complementary to the captured mRNA.

[0052] Referring to Fig. 6, after sample wash and injection of reverse transcriptase reagents, the mRNA sequence is transferred to the capture oligonucleotide 1010 via reverse transcription with template switching to obtain a cDNA with a sequence complementary to the mRNA. The cDNA also contains the entire sequence of the capture oligonucleotidelOlO, that is, the cDNA also contains the cell barcode sequence 1011, the unique molecular identifier sequence 1012, the capture sequence 1013, and the PCR handle sequence 1014.

[0053] S5: perform a PCR amplification reaction on the cDNA to obtain a cDNA library, and the cDNA library is sequenced.

[0054] Referring to Fig. 6, after digesting excess barcoded oligonucleotides that did not capture an mRNA using exonuclease enzyme I, on-chip PCR amplification of the cDNA is performed. The reagents and specific steps used for on-chip PCR amplification are common knowledge in the art, and will not be repeated here. After on-chip PCR amplification, high-throughput sequencing of the cDNA library is performed on a DNA sequencer to read the sequence information of the cDNA library. [0055] S6: the cell barcode sequence 1011 and the unique molecular identifier sequence 1012 are read according to a sequencing result, and the morphological characteristics and gene expression information of the cell in the microwell are integrated together. [0056] Specifically, after identifying the genes, reads are organized by their cell barcode sequence 1011 and individual UMIs are counted for each gene in each cell. By reading the cell barcode sequence 1011 in the sequencing result, the micro well 101 can be located, and the morphological characteristics (the bright-field image and the fluorescence image obtained in step S2) of the cell corresponding to the microwell can be obtained. Furthermore, the morphological characteristics and gene expression of the cell are integrated by a controller which includes an analysis software. In some embodiments, the analysis software is t-distributed stochastic neighbor embedding (tSNE), a data visualization tool that can reduce high-dimensional data to two-dimensional or three-dimensional, and then draw it into a graph. Referring to Fig. 7, the single-cell expression profiles are eventually be plotted two-dimensionally (tSNE) for visualized analysis and integrated with their morphological features. It is understandable that other analysis software or methods can also be used to analyze the morphological characteristics and gene expression profiles of single cells.

[0057] Referring to Fig. 8, the present application also provides a system 200 that integrates the morphological characteristics and gene expression of single-cell. The system 200 comprises a microfluidic device 10, an imaging device 30, a PCR machine 50, a sequencer 70, and a controller 90. The microfluidic device 10 is used to disperse the cell population into individual cells and mark the individual cell. The microfluidic device 10 comprises a microwell array composed of a plurality of microwells 101 and an interdigital electrode 103. Each micro well 101 comprises a plurality of capture oligonucleotides 1010, and each capture oligonucleotide 1010 comprises a known cell barcode sequence 1011, a unique molecular identifier sequence 1012, and a capture sequence 1013. There is a unique correspondence between the cell barcode sequence

1011 and the microwell 101 where it is located, and the cell barcode sequence 1011 is used to mark the cell in the micro well 1011. The unique molecular identifier sequence

1012 is used to label the captured mRNA. The imaging device 30 is used to capture images the cells 20 and record the morphological characteristics of the cells 20 for morphological analysis. In some embodiments, the imaging device 30 comprises a CCD camera and a microscope connected to the CCD camera. The PCR machine 50 is used for PCR amplification of cDNA to obtain a cDNA library. The sequencer 70 is used to sequence the cDNA library. The controller 90 is used to integrate the morphological characteristics and gene expression of the cells 20. In some embodiments, the controller includes an analysis software, and the analysis software comprises but is not limited to tSNE.

[0058] The application will be further described below in conjunction with specific embodiments.

EXAMPLE 1

[0059] HEK 293 T cells (human embryonic kidney cell line) and mouse 3T3 cell lines were mixed at the same concentration, and the mixture was analyzed for single-cell morphological characteristics and gene expression profiles. Individual cells were isolated in the microwells for imaging, as shown in Fig. 9. After cell lysis, reversed transcription, and on-chip PCR amplification on device, the cDNA library was purified using SPRI beads and measured by Agilent 2100 Bioanalyzer (DNA 7500 Kit) as shown in Fig.10. It can be seen from Fig. 10 that the average size of the cDNA was between 1000-2000 bp. The cDNA library was sequenced using a MiSeq sequencer, and the number of individual UMIs was plotted in Fig. 11. As can be seen from Fig. 11, 960 human-only transcripts, 884 mouse-only transcripts and 9 human-mouse mixed transcripts are collected, indicating a multiplet rate of 0.4%, which demonstrates a good performance of the single-cell isolation.

[0060] If cells from two cell lines are isolated into the same microwell and their mRNAs are captured by the same capture oligonucleotide 1010, the genes from the two species will share the same cell barcode sequence 1011, the proportion of which reveals the single-cell purity. The single-cell purity could be further improved by analyzing the cell images to screen out the cell doublets and multiplets. In this example, the single-cell purity is greater than 95%.

[0061] The recovery rate is calculated as the percentage of recovered cells (number of cell barcode sequences) to the total input cell number. In this example, the recovery rate is greater than 80%.

[0062] The contamination of ambient RNA from original biofluids or cell disruption decreases the accuracy of interpretation. Synthetic RNA controls are spiked into the sample with known concentration and sequence to evaluate the RNA contamination, which is defined as the percentage of recovered spiked RNAs (number of individual UMIs with spike-RNA sequence) to the total input spiked RNA amount. In this example, RNA contamination is less than 5%.

[0063] To characterize the detection sensitivity, mouse cells are spiked into human cells in different ratios of concentration ranging from 1:1 (50%) to 1:99 (1%). The sensitivity is determined as the percentage above which the mouse cells can be detected. In this example, the sensitivity is less than 5%, that is, the lower limit of detection of mouse cell concentration is 5%, and mouse cell can be detected as long as the concentration is greater than 5%.

[0064] In this application, a plurality of individual cells are placed in a plurality of microwells in the microfluidic device, and from the beginning (imaging) to the end (sequencing), each cell will be assigned a unique known cell barcode sequence to observe their phenotype before processing for sequencing. The cell barcode sequence in the capture oligonucleotide can be "read" in the microwells and also can be "read" from the sequence reads obtained from the cDNA library, thereby the genome/transcriptome data (mRNA sequence information) is linked to the observed phenotype of single-cell, so that the morphological phenotype is directly related to gene expression. This method focuses on integrating the morphological characteristics and gene expression profdes of isolated single cells. It has the characteristics of high efficiency, single-cell purity (greater than 95%), recovery rate (greater than 80%) and sensitivity (less than 5%), and low RNA contamination (less than 5%). It can facilitate fundamental biological studies, develop multi-dimensional biomarker signatures for diseases, and accelerate drug discovery and development.

[0065] The above descriptions are some specific implementation manners of the present application, but in the actual application process, they should not be limited to these implementation manners. For persons skilled in the art, other modifications and changes made according to the technical concept of this application should all belong to the protection scope of this application.