CN113960001A - AutoQT-FRET method based on primary imaging measurement system correction factor and application thereof - Google Patents

AutoQT-FRET method based on primary imaging measurement system correction factor and application thereof Download PDF

Info

Publication number
CN113960001A
CN113960001A CN202111211348.1A CN202111211348A CN113960001A CN 113960001 A CN113960001 A CN 113960001A CN 202111211348 A CN202111211348 A CN 202111211348A CN 113960001 A CN113960001 A CN 113960001A
Authority
CN
China
Prior art keywords
fret
cells
plasmid
cell
equal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111211348.1A
Other languages
Chinese (zh)
Other versions
CN113960001B (en
Inventor
陈同生
高璐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China Normal University
Original Assignee
South China Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by South China Normal University filed Critical South China Normal University
Priority to CN202111211348.1A priority Critical patent/CN113960001B/en
Publication of CN113960001A publication Critical patent/CN113960001A/en
Application granted granted Critical
Publication of CN113960001B publication Critical patent/CN113960001B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/645Specially adapted constructive features of fluorimeters
    • G01N21/6456Spatial resolved fluorescence measurements; Imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/6428Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N2021/6417Spectrofluorimetric devices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/6428Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
    • G01N2021/6439Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes" with indicators, stains, dyes, tags, labels, marks

Abstract

The invention discloses an AutoQT-FRET method based on a correction factor of a primary imaging measurement system and application thereof. The method comprises the steps of firstly utilizing a single-transfer plasmid cell sample to obtain the three-channel average fluorescence intensity I of each plasmid through E-FRET imagingDD、IAA、IDAAnd acceptor-sensitized fluorescence intensity FCAnd with IDDAnd FCDrawing a scatter diagram for the horizontal and vertical coordinates, and obtaining the range of the slope k value of each plasmid through least square normal linear fitting; then E-FRET imaging is carried out on the cell sample mixed with the multiple plasmids, background removal, crosstalk correction, region segmentation and region screening are carried out in sequence, and then the cells mixed with the multiple plasmids are distinguished according to the range of the slope k value of each plasmid so as to determine the space standard line of each plasmid; and finally, fitting a plane according to the space standard line to obtain system correction factors G and beta. The system correction factors G and beta determined by the method are accurate and are suitable for different detection systems.

Description

AutoQT-FRET method based on primary imaging measurement system correction factor and application thereof
Technical Field
The invention belongs to the technical field of Fluorescence Resonance Energy Transfer (FRET) detection, and particularly relates to an AutoQT-FRET method based on a correction factor of a primary imaging measurement system and application thereof.
Background
Quantitative Fluorescence Resonance Energy Transfer (FRET) detection is an important technology for studying intermolecular interactions in living cells, analyzing the molecular structure of oligomerized proteins, studying the regulatory mechanisms between proteins in signaling pathways, etc. [ Deal, J., et alT-based sensors of intracellular signals:A biological perspective of the history of FRET."Cellular Signalling(2020):109769.]。Zal[Zal,Tomasz,and Nicholas RJ Gascoigne."Photobleaching-corrected FRET efficiency imaging of live cells."Biophysical journal 86.6(2004):3923-3939.]And the like, a quantitative FRET (E-FRET) measuring method based on 3-cube is provided, and acceptor excitation crosstalk and donor emission crosstalk can be overcome simultaneously. The E-FRET method proceeds through three channels: the combination of the donor detection channel (DD) when excited by donor excitation light, the acceptor detection channel (DA) when excited by donor excitation light, and the acceptor detection channel (AA) when excited by acceptor excitation light eliminates spectral crosstalk, and at the same time, can eliminate spectral crosstalk according to E ═ Fc)/(Fc+G*IDD) The FRET efficiency E was calculated quantitatively. Wherein IDDRepresenting the intensity of donor fluorescence detected by the DD channel; g is a system sensitization quenching conversion factor; fcIs the fluorescence intensity sensitized to the acceptor in the FRET sample.
Accurate measurement of the sensitization-quenching conversion factor G is the key of the quantitative FRET detection technology. Accurate measurement of the correction factor G is therefore particularly important in relation to the intrinsic properties of the instrument and of the fluorescent molecule. In the E-FRET method, zal et al [ Zal, Tomasz, and Nichlas RJ Gascoigne. "Photobleaching-corrected FRET interference of live cells." Biophysical journel 86.6(2004):3923 and 3939.] are measured by Photobleaching with a reference sample identical to the donor and acceptor of the sample to be measured to obtain the sensitization-quenching conversion factor G. However, in the actual measurement, the accurate measurement of the correction factor G is influenced by the defects that the photobleaching method has incomplete photobleaching in the photobleaching process and can only carry out measurement on fixed cells.
Alexis coulomb et al [ Coillomb, Alexis, et al. "QuanTI-FRET. a frame for quantitative FRET measurements in measuring cells." Scientific reports 10.1(2020):1-11 in 2020.]A quantitative triple image FRET method (QuanTi-FRET) is further provided on the basis of the E-FRET method, and a correction factor G and a correction factor beta value related to the concentration ratio of a donor and an acceptor can be simultaneously measured and obtained. The method also obtains an image I related to the fluorescence intensity of the donor and acceptor through the combination of three channels DD, DA and AADD、IDAAnd IAA. Is subdividedIs distinguished byDD、IDAAnd IAAAnd drawing a spatial scatter diagram for the x axis, the y axis and the z axis, and performing plane least square fitting on the spatial scatter diagram to obtain direction vector absolute values of a spatial plane as correction factors G and beta respectively. The QuanTi-FRET method requires the measurement of at least i (i.gtoreq.2) immobilized plasmids individually. However, the method needs to switch multiple cell samples, is complex to operate, and cannot ensure the same contrast and background signal when imaging the multiple cell samples respectively.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides an AutoQT-FRET method based on a correction factor of a primary imaging measurement system, which is a quantitative three-image FRET method further provided on the basis of an E-FRET method.
The invention also aims to provide application of the AutoQT-FRET method based on the correction factor of the primary imaging measurement system.
The purpose of the invention is realized by the following technical scheme:
an AutoQT-FRET method based on a correction factor of a primary imaging measurement system comprises the following steps:
(1) sample preparation
Reference sample: respectively transfecting cells with the I FRET plasmids with different FRET efficiencies, and after the transfected plasmids are successfully expressed, respectively culturing the cells expressing the different FRET plasmids in different sample container carriers until the cells are attached to the walls, and marking as a reference sample 1 and a reference sample 2 … … reference sample i; wherein i is more than or equal to 2 (and is an integer);
② the sample to be tested: respectively transfecting cells with the i FRET plasmids with different FRET efficiencies, merging the cells expressing the different FRET plasmids in the same sample container carrier after the transfected plasmids are successfully expressed, and culturing until the cells adhere to the wall to obtain a sample to be detected; wherein i is more than or equal to 2 (and is an integer);
(2) E-FRET spectroscopic imaging
Thirdly, M fields are selected from the reference sample 1 (the reference sample transfected with the 1 st plasmid) obtained in the step I for E-FRET imaging, and then M cells are selected under the 1 st field to obtainObtaining three channels of average fluorescence intensity IDD1-1、IDD1-2……IDD1-m,IAA1-1、IAA1-2……IAA1-m,IDA1-1、IDA1-2……IDA1-mAnd calculating the acceptor-sensitized fluorescence intensity FC1-1、FC1-2……FC1-m(ii) a Then with IDD1-1、IDD1-2……IDD1-mAs the abscissa, FC1-1、FC1-2……FC1-mA scatter plot is plotted for the ordinate and a least squares linear fit is performed to obtain the slope k at the 1 st field of view of the reference sample 11-1(ii) a By analogy, the slope k under the 2 nd to M th field of the reference sample 1 is obtained1-2、k1-3……k1-MAnd determining the maximum value and the minimum value to obtain the slope k of the reference sample 11Range of values (k)Minimum value≤k1≤kMaximum value) (ii) a Wherein M is more than or equal to 10, and M is more than or equal to 10 (both integers);
replacing the reference sample 1 with the reference sample 2, and repeating the step c to obtain the slope k of the reference sample 22A range of values; and so on until the slope k of the reference sample i is obtainediRange of values (denoted as k)2、k3……ki);
Selecting 1 visual field from the sample to be tested obtained in the step II to carry out E-FRET imaging to obtain three channels I thereofDD、IAA、IDAImage, and the length and width of the area of 100 x 100 pixels are respectively corresponding to IDD、IAA、IDAPerforming region segmentation on the image to obtain n cell region images IDD-S1、IDD-S2……IDD-Sn,IAA-S1、IAA-S2……IAA-Sn,IDA-S1、IDA-S2……IDA-Sn(ii) a Then, starting from the 1 st cell region, the proportion S of the effective data in the cell region is calculated1And measuring the three-channel average fluorescence intensity of all pixel points in the cell region, wherein the total number of the pixel points is L, and the pixel points are marked as IDD-1、IDD-2……IDD-L,IAA-1、IAA-2……IAA-L,IDA-1、IDA-2……IDA-LAnd simultaneously calculating the receptor-sensitized fluorescence intensity FC-1、FC-2……FC-LThen further calculating to obtain the ratio F of the 1 st cell areaC-1/IDD-1、FC-2/IDD-2……FC-L/IDD-LPerforming unimodal Gaussian fitting to obtain unimodal Gaussian function fitting degree value X of the cell region1And according to IDD-1、IDD-2……IDD-LAnd FC-1、FC-2……FC-LDrawing a scatter diagram, and then obtaining a slope k according to least square normal linear fitting1' value; by analogy, the proportion S of the effective data in the 2 nd cell area to the nth cell area is calculated in the same way2、S3……SnDegree of fit value X of single-peak Gaussian function2、X3……XnAnd slope k2′、k3′……kn'; wherein n is more than or equal to 25 (the value of n is related to the image size, for example, the 2048 × 2048 pixel picture is divided by 100 × 100, n is 400, the 1024 × 1024 pixel picture is divided by 100 × 100, n is 100, the 512 × 512 pixel picture is divided by 100 × 100, and n is 25);
s degree1< 30% or X1(ii) < 90%, removing the cell region if S1Not less than 30% and X1More than or equal to 90%, the cell region is retained (excluding the effect on plasmid standard line determination when multiple plasmids are present in a region); in the same way, if S2< 30% or X2(ii) < 90%, removing the cell region if S2Not less than 30% and X2More than or equal to 90 percent, the cell area is reserved, and the like, R cell areas in total are obtained, and the corresponding slopes are marked as k1″、k2″……kR"; wherein R is more than or equal to 10;
(3) cell classification and determination of plasmid standard lines
Seventhly, obtaining k according to the step (c)1″、k2″……kRSize of "and range of slope values obtained in Steps (c) and (dk1、k2……kiCell classification was performed: if k is1At k1In the range of (i), i.e. k1"corresponding cells are cells transfected with plasmid 1; if k is1At k2In the range of (i), i.e. k1"corresponding cells are cells … … f k transfected with plasmid 21At kiIn the range of (i), i.e. k1"corresponding cells are cells transfected with the i plasmid; by analogy, respectively combine k2″、k3″……kR"performing classification; then classifying the cells according to the I of each pixel point in all cell regions of each plasmidDDFc and IAAA spatial scattergram is drawn for the x, y and z axes, and a spatial straight line is fitted by the least square method to serve as a spatial standard line l for each plasmid1、l2……li
(4) Obtaining correction factors for E-FRET systems
Eighty to the step of obtaining the space standard line l of each plasmid1、l2……liPerforming plane least square fitting to obtain a plane equation Ax + By which is z; wherein A and B are constants, and B is a system correction factor beta1A/beta is the system correction factor G1(corresponding to the system correction factors β, G at the 1 st field of view);
ninthly, repeating the steps from five to eight, randomly selecting N different visual fields from the sample to be measured obtained in the step two to carry out E-FRET spectral imaging, and respectively calculating to obtain a system correction factor beta2、β3……βNAnd G2、G3……GNThen, taking the average value of the correction factors to obtain system correction factors G and beta; wherein N is more than or equal to 10.
The types of i in the step (1) are 2-4; preferably 4.
The FRET plasmid in the step (1) is a FRET standard plasmid, and a FRET plasmid with different FRET efficiencies can be constructed by adjusting the distance between a donor and an acceptor (the larger the distance is, the smaller the FRET efficiency is); the efficiency of FRET plasmids can also be measured in advance by methods such as a lifetime measurement method, an E-FRRT measurement method or spectral linear separation and the like (the invention does not need to know the specific efficiency of each plasmid exactly); the FRET plasmid comprises a donor and an acceptor, and the number (concentration) ratio of the donor to the acceptor is 1: n or n: 1 FRET plasmid (standard plasmid of C-V system), n is more than or equal to 1; preferably, the ratio of the number (concentration) of the donor to the number (concentration) of the acceptor is 1: 1, a FRET plasmid; more preferably at least two of C5V, C17V, C32V and CTV.
The sample container carrier in the first step and the second step is a container carrier capable of realizing cell culture and imaging, and comprises but is not limited to a cell culture dish, a pore plate, a glass slide and the like; the invention selects the conventional cell culture dish to culture and image the cells.
The cells in the first step (1) and the second step are preferably Hela cells.
The value range of M in the third step (2) is that M is more than or equal to 20; preferably M is more than or equal to 40; more preferably M.gtoreq.50; the larger the number of cells selected herein, the better.
The selection of the m cells in the third step (2) is preferably realized by the following method: placing a sample container carrier (such as a culture dish) containing cells with transfection plasmids under a wide-field fluorescence microscope, selecting M visual fields, then selecting cells from each visual field, and selecting M cells for measurement; wherein M is more than or equal to 10 (preferably M is more than or equal to 20, more preferably M is more than or equal to 40, and more preferably M is more than or equal to 50), M is more than or equal to 10 (preferably M is more than or equal to 80, more preferably M is more than or equal to 250, and more preferably M is more than or equal to 500); that is, M cells are from M fields of view (all cells in the vector. gtoreq.m cells in all M fields of view).
The number of the selected cells in each visual field is more than 1; preferably 1 to 7; more preferably 1 to 3; more preferably 2 to 3.
The value range of m in the third step (2) is that m is more than or equal to 80; preferably m is more than or equal to 250; more preferably m.gtoreq.500; the larger the number of cells selected herein, the better.
The three-channel average fluorescence intensity in the step (2) is donor fluorescence intensity I detected in the donor fluorescence detection channel when the donor excitation light is excitedDDAcceptor fluorescence detected in the acceptor fluorescence detection channel upon excitation with donor excitation lightLight intensity IDAAcceptor fluorescence intensity I detected in acceptor fluorescence detection channel when excited by acceptor excitation lightAATo 1, pairDD、IDA、IAAThe original image uses Gaussian filtering to determine the range of the background, and background subtraction and crosstalk correction are sequentially carried out on the original image.
The receptor sensitization fluorescence intensity in the third step (2) is the pair IDAImage subtraction FRET channel (I)DA) The fluorescence of the excitation crosstalk of the acceptor and the fluorescence of the emission crosstalk of the donor are corrected to obtain the acceptor-sensitized fluorescence intensity (Fc).
Three channels I described in step (2)DD、IAA、IDAThe image is an image of 2048 × 2048 pixels, 1024 × 1024 pixels, or 512 × 512 pixels; preferably 2048 × 2048 pixels.
The proportion of the effective data in the cell area in the fifth step (2) is calculated by the following method: in the cell region, the background is set to 0 and the non-background is set to 1, and then the proportion of the sum of non-background pixels to the sum of pixels of the entire region is calculated as the proportion of valid data.
The value of R in the step (2) is related to the value of n in the step (2), and the value of n is related to the image size, namely R is more than or equal to 10 and less than or equal to n, namely at least 10 cell areas are ensured to exist, the data of each plasmid can be ensured to exist, and a spatial standard line can be obtained; that is, when the image size is 2048 × 2048 pixels, n is 400, at which time 10 ≦ R ≦ 400; when the image size is 1024 multiplied by 1024 pixels, n is 100, and R is more than or equal to 10 and less than or equal to 100; when the image size is 512 × 512 pixels, n is 25, and at this time, R is 10 ≦ 25.
The cell classification in step (3) can also be achieved by:
(a) selecting N 'cell regions under N' visual fields for E-FRET imaging on the sample to be detected obtained in the step II, and obtaining three-channel average fluorescence intensity I pixel by pixelDDˊ、IDAˊ、IAAAnd calculating the intensity of receptor-sensitized fluorescence F thereofCAnd with I of each pixelDDˊ、FCˊ、IAAIs emptyAn inter-coordinate point;
(b) respectively calculating the 1 st spatial coordinate point to the spatial standard line l obtained in the step (c)1、l2……liDistance d of1、d2……di(ii) a Then d is put1、d2……diSorting by magnitude of value, if d1And if the minimum value is reached, classifying the spatial coordinate point as a spatial standard line l1Namely, the cell region corresponding to the space coordinate point belongs to the cell transfected with the FRET plasmid of the 1 st type; in the same way, if d2And if the minimum value is reached, classifying the spatial coordinate point as a spatial standard line l2That is, the region of the cell corresponding to the spatial coordinate point belongs to the region d of the cell … … transfected with the FRET plasmid of the 2 nd speciesiAnd if the minimum value is reached, classifying the spatial coordinate point as a spatial standard line liNamely, the cell region corresponding to the space coordinate point belongs to the cells transfected with the i-th FRET plasmid; and by analogy, classifying the spatial coordinate point corresponding to each pixel to realize the classification of the cells (plasmids).
The preferable value range of N in the step (5) is as follows: n is more than or equal to 15; more preferably, N.gtoreq.40.
The AutoQT-FRET method based on the correction factor of the primary imaging measurement system is applied to the correction factor of the measurement system.
The system correction factors are system correction factors G and beta.
The basic principle of the invention is as follows:
(1) for FRET samples imaged by a three-channel method, the fluorescence intensity I of the three channels for the acceptor is obtained by imagingDD、IDA、IAA. It can be written as a function of photophysical and instrumental parameters, the number n of donors in the pixel under considerationDAnd the number of acceptors nAFluorophore and FRET probability (efficiency) E:
IAA=nA*LAA AexAAem Adet (1)
IDD=nD*LDD Dex*(1-E)*ΦDDem Ddet (2)
IDA=nD*LDD Dex*E*ΦAAem Adet+nD*LDD Dex*(1-E)*ΦDDem Adet+nA*LDA DexAAem Adet (3)
wherein L isAExcitation intensity at a wavelength selected to excite fluorescent molecule a; l isDExcitation intensity at a wavelength selected to excite fluorescent molecule D; sigmaA AexIs the absorption cross section of A at the excitation wavelength of A; sigmaD DexIs the absorption cross section of D at the excitation wavelength of D; sigmaA DexIs the absorption cross section of D at the excitation wavelength of A; phiAQuantum yield of a; phiDQuantum yield of D; etaAem AdetA detection efficiency for emitting photons in detection channel a; etaDem DdetA detection efficiency for emitting photons in detection channel D; etaDem AdetThe detection efficiency of a for emitting photons in the detection channel D.
To simplify the operation, four correction factors are defined by equation (4) to simplify equations (1) - (3): the direct excitation correction factor is a; the leakage correction factor is d; the correction factor for different detection efficiencies in the two channels is G; the correction factor for the different excitation efficiencies in the two channels is β.
Figure BDA0003308966810000051
(2) By deducting FRET channel (I)DA) Fluorescence of excitation crosstalk of intermediate acceptor and fluorescence of emission crosstalk of donor to obtain fluorescence intensity F sensitized to acceptor in FRET samplec
Fc=IDA-aIAA-dIDD (5)
Thus FRET efficiency can be expressed as:
Figure BDA0003308966810000061
(3) meanwhile, the stoichiometric ratio is defined as the relative number of donor molecules relative to the total number of fluorophores in each pixel:
Figure BDA0003308966810000062
from the equations (1) (2), n can be derivedDAnd nAAnd inserted into equation (7). By simplifying with the use of the correction factor β defined in equation (4), equation (7) reduces to:
Figure BDA0003308966810000063
(4) to accurately obtain the correction factors G and β, equation (8) is rewritten as:
Figure BDA0003308966810000064
this is represented by { IDD、Fc、IAAThe equation of the plane in the defined three-dimensional space. If the stoichiometric ratio S is known, the experimental data { IDD、Fc、IAAThe correction factors G and β are determined by least squares fitting to a plane.
(5) Multiple sets of { I } sets can be obtained by measuring 2 or more samples of FRET tandem structure containing 1 donor and 1 acceptor (number ratio 1: 1)DD、Fc、IAAAnd through the pair { I }DD、Fc、IAAAnd performing least square fitting on a space plane, wherein absolute values of space vectors of the fitting plane are correction factors G and beta in the formula (9) respectively.
(6) K value range for cells transfected with FRET-derived plasmids using single-transfer plasmid samples: the correction factor G is only related to the instrumentAnd the intrinsic properties of the fluorescent molecule, whereas the FRET efficiency E is associated with a different donor FRET pair. Therefore, I shown in the following formula (10)DDSlope of Fc
Figure BDA0003308966810000065
Showing tandem FRET plasmids having different FRET efficiencies.
Figure BDA0003308966810000066
Compared with the prior art, the invention has the following advantages and effects:
(1) the invention belongs to the technical field of FRET quantitative detection, and the method utilizes a cell sample (the sample contains 2 or more than 2 FRET tandem plasmids with the ratio of donor to acceptor being 1: 1) to simultaneously measure a system correction factor G (system sensitization quenching conversion factor) and beta (excitation scale factor of two channels), and specifically comprises the following steps: obtaining the k value range of each plasmid by using E-FRET data of single-transfer plasmid cells; sequentially carrying out processes of background removal, crosstalk correction, region segmentation, region screening, determination of a space standard line of each plasmid and the like on cell E-FRET image data mixed with various plasmids; fitting a plane according to the determined space standard line of each plasmid and obtaining a system correction factor; and the spatial standard lines can be used to classify one dish of images of multiple granulocytes.
(2) The invention provides an AutoQT-FRET method based on a correction factor of a primary imaging measurement system, thereby improving the intelligence and the accuracy of correction parameters G and beta of the measurement system and being the basis for developing an intelligent FRET correction system.
(3) The G and the beta obtained by the method are more accurate and are suitable for different detection systems, so that the method can greatly promote the application range of the (AutoQT-FRET) method, and further improve the application range of the FRET detection technology in cell biology.
(4) Compared with the traditional method, the method can finish the measurement of the correction factor only by one cell sample (each plasmid is not required to be separately measured, and the problems of different contrast and background signals exist in the separate measurement), and in addition, the method can be used for measuring and calculating the correction factor on line (for example, data is processed by a computer), so that the use threshold is reduced, and the efficiency is improved.
Drawings
FIG. 1 is a diagram showing the procedure and results of determining the range of different plasmids k according to the present invention; wherein, (a) is DA channel image for E-FRET spectral imaging of cell sample of single-transfer four plasmids; (b) to select the cell region encircled by the red box in (a), the cell region I is subjected toDDDrawing a scatter diagram by the Fc data, and fitting by a least square method to obtain a linear graph of the slope k of the cell region; (c) to perform a statistical analysis of the k-values of a total of 265 cells in 40 fields, a statistical map of the range of slope k-values for different plasmids was obtained.
FIG. 2 is a graph showing the process and results of determining spatial standard lines of different plasmids according to the present invention; wherein, (a) a petri dish is transfected with E-FRET spectroscopic imaging data of any one field of four different plasmid cell samples; (b) to select the cell region circled by the red box in (a), the cell region was subjected to pixel-by-pixel Fc/IDDCounting the frequency distribution of the values and obtaining a fitting graph by single-peak Gaussian fitting; (c) results are presented as spatial standard lines for three plasmids C5V, C17V, C32V, CTV in this field; (d) a visual field data classification template map was prepared based on (a) the distribution of different plasmids in the visual field (in the figure: blue is data of cell-transfected C5V plasmid, red is data of cell-transfected C17V plasmid, yellow is data of cell-transfected C32V plasmid, and purple is data of cell-transfected CTV plasmid).
FIG. 3 is a graph of the statistical results of correction factors obtained using different methods.
Detailed Description
The present invention will be described in further detail with reference to examples, but the embodiments of the present invention are not limited thereto. Reagents, methods and apparatus used in the present invention are conventional in the art unless otherwise indicated. The test methods in the following examples, in which specific experimental conditions are not specified, are generally performed according to conventional experimental conditions or according to the experimental conditions recommended by the manufacturer. Unless otherwise specified, reagents and starting materials for use in the present invention are commercially available.
Example 1
1. Plasmids
Given donor-acceptor pairs: the donor is gene coding fluorescent protein Cerulean (short for C), and the acceptor is gene coding fluorescent protein Venus (short for V);
FRET tandem plasmid structure C5V: a FRET tandem plasmid structure consisting of a 5 amino acid linker linking C and V;
FRET tandem plasmid structure C17V: a FRET tandem plasmid structure consisting of a 17 amino acid linker linking C and V;
FRET tandem plasmid structure C32V: a FRET tandem plasmid structure consisting of a 32 amino acid linker linking C and V;
FRET tandem plasmid structure CTV (Cerulean-TRAF-Venus): a FRET tandem plasmid structure consisting of a 229 amino acid linker linking C and V; wherein TRAF is a receptor-associated factor domain comprising a long chain tumor necrosis factor of 229 amino acids;
the source is as follows: the donor plasmids Cerulean (C), the acceptor plasmid Venus (V) and the FRET reference plasmids C5V, C17V, C32V and CTV were purchased from the Addge plasmid library, USA.
2. Wide-field spectral microscopic imaging system
Wide field fluorescence microscopes were manufactured by Olympus, Japan, under model IX 73. The light source was a mercury lamp of the olympus HGLGPS series, japan. The objective lens is an oil lens with the magnification of 40 and the numerical aperture of 1.3(40 multiplied by 1.3NA), one exciting rotating wheel provided with four exciting sheets, one rotating wheel provided with eight cubes (each cube can be provided with one exciting sheet, one splitting sheet and one emitting sheet respectively), one electric emitting rotating wheel provided with six emitting sheets, and the electric emitting rotating wheel is externally connected with a CCD camera. The wavelength of the excitation light is selected by rotating the excitation wheel.
3. Cell culture and plasmid transfection
Hela cells were obtained from university of river, Guangzhou, China, and 10% (v/v) of newborn bovine serum was added to DMEM medium and cultured in an incubator at 37 ℃ containing 5% (v/v) carbon dioxide. Trypsin for cellsDigesting, transferring to a cell culture dish, culturing for 24 hours, and then using an in vitro transfection reagent Turbofect when the cells grow to 70-90 percentTMThe plasmid was transiently transferred into cells. The specific steps of transfection are as follows:
(1) taking a sterilized EP tube, firstly adding 100 mu L of serum-free DMEM medium, then adding 1-2 mu L of transfection reagent, then adding 1-2 mu L (500-600 ng/mu L) of plasmid, gently blowing and beating for 6-8 times, and then standing for 20 minutes;
(2) after 20 minutes, adding 100 mu L of serum-free DMEM medium into the EP tube, and gently mixing;
(3) washing cells in the culture dish for 2-3 times by using a serum-free DMEM medium or PBS (phosphate buffered saline), mainly washing off dirt such as dead cells and the like, then transferring the mixture in the step (2) into the culture dish, and putting the culture dish back into the culture box for 4-6 hours;
(4) (ii) for samples of single transfer plasmids C5V, C17V, C32V and CTV in a culture dish, E-FRET imaging can be respectively carried out after the steps (1) to (3) are completed, and donor fluorescence intensity I detected in a donor fluorescence detection channel when excited by donor excitation light is obtainedDDIntensity of acceptor fluorescence I detected in acceptor fluorescence detection channel upon excitation with donor excitation lightDAAnd the intensity of acceptor fluorescence I detected in the acceptor fluorescence detection channel when excited by the acceptor excitation lightAA(i.e., three-channel mean fluorescence intensity IDD、IDA、IAA) And calculating the acceptor-sensitized fluorescence intensity FC(ii) a Wherein, FCTo subtract FRET channel (I)DA) Fluorescence intensity sensitized to acceptor in FRET sample obtained by fluorescence of excitation crosstalk of acceptor and fluorescence of emission crosstalk of donorc
② for a sample containing the multi-plasmid (C5V, C17V, C32V and CTV) in the cells of a culture dish, the preparation method comprises the following steps: after the cells independently transfecting plasmids are subjected to step (3) for 4-6 hours respectively, absorbing transfection liquid, cleaning the cells in the culture dish for 2-3 times by using serum-free DMEM culture medium or PBS, adding 100 mu L of pancreatin cell digestive fluid into the culture dish to digest the cells transfected with different plasmids in different culture dishes, absorbing the digestive fluid, slightly blowing the cells for 2-3 times by using 100 mu L of serum-free DMEM culture medium, putting cell turbid liquid blown from different culture dishes into the same EP tube, uniformly mixing for 2-3 times, and putting the cell turbid liquid into a new culture dish again, wherein the ratio of the DMEM culture medium to fetal bovine serum is 10: 1 (volume ratio), adding fetal bovine serum, culturing for 8-12 hours, and allowing the cells to adhere to the wall for experiment.
4. Multi-plasmid cotransfer sample imaging process
CTV, C5V, C17V and C32V were transfected into Hela cells separately according to the plasmid transfection procedure described above, and then pooled into one dish for expression maturation (cell attachment) followed by E-FRET spectroscopy. The wide field fluorescence microscope was configured as follows: light of 436/20nm and 470/20nm is used as exciting light of Cerulean and Venus respectively, and channels of 510nm and 530nm are used as detection channels of Cerulean and Venus fluorescence respectively. Under 436nm excitation, the detector collects images of all detection channels, and under 470nm excitation, the detector collects images of only 530nm detection channels, and 3 images are recorded as E-FRET data of one field of view.
5. Data processing
5.1 determining the k value range of FRET plasmid transfected cells according to E-FRET data of a multi-plasmid single-transfer sample, and the specific process is as follows:
in the step (4), 40 fields (10 fields for each plasmid sample) are randomly selected from the multi-plasmid single-transfer sample, then cells are selected in each field, and 265 cells are selected for measurement in the experiment. Obtaining an image I after background subtraction and crosstalk correction according to E-FRET spectral imagingDD、Fc、IAA. Drawing I pixel by pixel for each field of viewDDFc dot plot (as I)DDAbscissa and Fc ordinate), and performing least square fitting to calculate a slope k value of each field, which is as follows:
as shown in FIG. 1(a), for the FRET reference plasmid C5V, 1 field was selected for E-FRET spectroscopic imaging, and then the cell region enclosed by the red box in FIG. 1(a) was selected (i.e., the field was selected first, and then the cell region was selected under the field, which has no special requirement and is free to be selectedSelecting and controlling the size of the selected cell region, such as whole cells or part of cells, wherein the more data, the more accurate the measurement, and 265 cells are selected in total), and obtaining the three-channel average fluorescence intensity I of each cell region respectivelyDD、IDA、IAAAnd calculating the acceptor-sensitized fluorescence intensity F according to the calculated resultC. Then based on I of the selected cell regionDDDrawing a scatter diagram by Fc data, and fitting by a least square method to obtain a diagram of a slope k fitting straight line of the cell region, wherein the diagram is shown in a figure 1 (b); by analogy, E-FRET spectroscopic imaging was performed on cell samples of the single-transfer FRET reference plasmids C5V, C17V and C32V, and then slope k was calculated by fitting all cell regions in 10 fields of each plasmid. The slope k values of all cell regions were statistically classified according to different plasmids, resulting in the range of slope k values for different plasmids, as shown in FIG. 1 (c). The finally determined ranges of the slope k values of the different plasmids are respectively: C5V: k is more than or equal to 5.5C5V≤8.5;C17V:3.2≤kC17V≤5;C32V:2≤kC32V≤3;CTV:0.2≤kCTV≤1.3。
TABLE 1 slope k value ranges for different plasmids
Species of plasmid C5V C17V C32V CTV
Slope k [5.5,8.5] [3.2,5] [2,3] [0.2,1.3]
Here, the measurement can be performed on Hela cells transfected with CTV, C5V, C17V and C32V individually (different dishes), and the cells can be selected individually (the number of cells to be selected can be selected as the case may be, preferably 10 or more fields (preferably 20 or more) are selected for each dish transfected with a different plasmid, and 1 or more cells (preferably 3 to 7 cells) are selected for each field). The experiment chose to measure on different dishes, i.e. directly on single-transfer plasmid samples in step 4.
5.2 determining the spatial standard line for each plasmid based on E-FRET data from multiple plasmid cotransformation samples:
5.2.1 Classification of plasmids according to the slope k values and determination of the spatial Standard line for each plasmid
For a total of 223 cells in 40 visual fields of the sample (containing C5V, C17V, C32V and CTV) co-transformed by multiple plasmids in the step 3, calculating a slope k value pixel by pixel for each cell region, classifying according to k value ranges of 4 different plasmids determined in the step 5.1, and then performing space straight line least square fitting to obtain standard lines of different plasmids, wherein the specific steps are as follows:
as shown in fig. 2(a), after performing E-FRET spectral imaging on a cell sample co-transformed with four plasmids, original images DD and AA are obtained; obtaining I after background deduction and crosstalk correctionDD、IDA、IAAAn image; according to the area length and width I of 100 x 100 pixelsDD、IDA、IAAAnd carrying out region segmentation on the image. Selecting the cell region (i.e. 100 x 100 pixel region) encircled by red box in fig. 2(a), calculating the proportion of effective data in the region (in the region, the background is already set to 0 and the non-background is set to 1 in the background subtraction process, so that the proportion of the sum of non-background pixels to the sum of pixels in the whole region is the effective data proportion) pixel by pixel, namely S-10.86. The cell region was then subjected to pixel-by-pixel Fc/IDDOf valueFrequency distribution statistics (i.e. for the first selected cell, in this cell region, starting from the 1 st pixel point, obtaining E-FRET data Fc and I thereofDDThen calculating Fc/IDDA value; then obtaining E-FRET data Fc and I of the 2 nd pixel pointDDThen it calculates Fc/IDDThe value … … is analogized until the last pixel point of the cell region (the size of the map used in the present invention is 2048 × 2048 pixels (pixels)); a unimodal gaussian fit was performed, see fig. 2(b), to obtain the degree of fitting of the unimodal gaussian function for this cell region, i.e., X ═ 0.48 (gaussian fit was performed using matlab software on the blue points of fig. 2(b) obtained statistically, to obtain the gaussian function (red line)). Finally, the slope k value of the least square fitting of the cell region is calculated to be 5.23, and the specific calculation method is the same as the step 5.1 (corresponding to I through each pixel point in the cell region)DDAnd drawing a scatter diagram by the Fc data, and fitting a straight line to obtain a slope k value).
Note: the above is based on Aoyin et al [ Yin, Ao, et al ] "Measuring calibration factors by imaging a dis of cell expression measurements" Cytometry Part A (2021).]The Fc/I of the same plasmid is proposedDDThe frequency histogram should satisfy a monomodal gaussian distribution, and Fc/IDD frequency histograms of two or more plasmids show the presence of multiple gaussian peaks corresponding to the number of plasmids. According to the principle, the Fc/I is calculated region by region for the regions screened out in the stepDDFrequency distribution and fitting with a unimodal gaussian function.
Obtaining effective data proportion S of all cell areas in the visual field, Gaussian function fitting degree X and slope k according to the same method, and respectively recording the effective data proportion S, the Gaussian function fitting degree X and the slope k as S-1、S-2……S-223,X-1、X-2……X-223,k-1、k-2……k-223. And respectively judging whether the effective data proportion in each region exceeds 30%, and if the effective data proportion in each region is less than 30%, removing the cell region. If the cell area is larger than or equal to 30%, the cell area is reserved; and then judging whether the fitting degree of the Gaussian function in each remaining cell area is greater than or equal to 90%, and if the fitting degree of the Gaussian function in each remaining cell area is less than 90%, removing the cell area. If greater thanEqual to 90% the cell area is retained; finally, the slope k values of the remaining cell regions are determined according to the k value range (k is more than or equal to 5.5) of different plasmids determined in the step 5.1C5V≤8.5;3.2≤kC17V≤5;2≤kC32V≤3;0.2≤kCTV≦ 1.3) for classification (the k-value range is determined by the maximum and minimum values of the valid data measured for each plasmid, with a small amount of data error possible, e.g. data with k-values outside the above ranges are obtained, which can be re-measured or removed), i.e. the remaining cell area is divided to determine which plasmid (C5V, C17V, C32V or CTV) it belongs to. Then, the classified cell region set is subjected to fitting of a space straight line least square method (namely I corresponding to the classified cell regionDDFc and IAADrawing a spatial scatter diagram for x, y and z axes, and fitting) to obtain spatial standard lines l of different plasmids1、l2、l3、l4See fig. 2 (c).
5.2.2 Classification of visual field data by spatial Standard lines of different plasmids
In addition to the above classification of plasmids using slope k values, the four spatial standard lines l of plasmids obtained in step 5.2.1 can be used1、l2、l3、l4(corresponding to 4 kinds of plasmids) and the formula (11), the visual field data of FIG. 2(a) was classified, and a template map for classification of the visual field data was prepared:
Figure BDA0003308966810000111
wherein
Figure BDA0003308966810000112
Is a standard line Li(i is 1,2 … … n) and point M is a straight line LiAt any point, point M0Is the coordinates of the spatial point sought.
Firstly, as shown in FIG. 2(a), E-FRET spectral imaging is carried out on a cell sample co-transformed with four plasmids to obtain original images DD, DA and AA, and background subtraction and subtraction are carried outAfter crosstalk correction, I is obtainedDD、Fc、IAA(ii) a Then with IDD、Fc、IAACalculating the coordinate of a space point P corresponding to the 1 st pixel point according to the formula (1) for x, y and z axes as { I }DD、Fc、IAATo the space standard line l1Distance d of1
Secondly, according to the same method, calculating to obtain a space point P { IDD、Fc、IAATo the space standard line l2Distance d of2Space point P { IDD、Fc、IAATo the space standard line l3Distance d of3And spatial point P { IDD、Fc、IAATo the space standard line l3Distance d of4
Comparison d1、d2、d3Size, the spatial point P { I }DD、Fc、IAAClassify as a spatial standard line of closest distance, i.e. d1At the minimum, the spatial point P is classified as a spatial standard line l1In the same way, e.g. d2At the minimum, the spatial point P is classified as a spatial standard line l2E.g. d3At the minimum, the spatial point P is classified as a spatial standard line l3E.g. d4At the minimum, the spatial point P is classified as a spatial standard line l4
And fourthly, classifying all the pixel points in the view to obtain four data set image templates of C5V, C17V, C32V and CTV (namely, marking the background in the view as 0, marking the pixel point classified as C5V as 1, marking the pixel point classified as C17V as 2, marking the pixel point classified as C32V as 3 and marking the pixel point classified as CTV as 4). The four datasets were merged (simple matrix addition, each dataset being an n x n size image; where n is the original image size 1024 x 1024 or 512 x 512) and the visual field data sort template map was made with different color labels, see fig. 2 (d): blue is data of cells transfected with C5V plasmid, red is data of cells transfected with C17V plasmid, yellow is data of cells transfected with C32V plasmid, and purple is data of cells transfected with CTV plasmid.
5.3 calculate the system correction factors G and β:
four plasmid spatial standard lines l obtained in step 5.2 were used1、l2、l3、l4And equation (10), calculating system correction factors G and β; for the four plasmid space standard lines l obtained in step 5.21、l2、l3、l4And performing least square fitting on the space plane, wherein the fitting plane equation should satisfy the form of Ax + By ═ z (the distance from each straight line to the plane is calculated, and the four straight line-to-plane distances and the plane with the minimum distance are the fitting plane). The visual field data is fitted to obtain a plane equation of 2.01IDD+0.59Fc=IAAAccording to formula (10), wherein G β is 2.01, β is 0.59; the fit therefore gives a G value of 3.70 for this field.
5.4 authentication
In order to test whether the data of the correction factors obtained by the method is accurate, the data of 200 cells in 22 fields selected by a co-transformed sample of various plasmids are analyzed by using an AutoQT-FRET method, a two-hybrid-multi-plasmid method and a TP-G method, and the specific steps are as follows:
(1) according to the two-hybrid-multi-plasmid method (reference: Butz, E.S., et al. "Quantifying macromolecular interactions in living cells using FRET two-hybrid assays," Nature Protocols 11.12(2016):2470.), the template image matrixes of the four plasmid datasets C5V, C17V, C32V and CTV classified in step 5.2.2 are respectively multiplied by the template image matrixes of the plasmid datasets obtained in step 5.2.2 to obtain I after background subtraction and crosstalk correctionDD、Fc、IAAImage matrix, obtaining I of each plasmidDD、Fc、IAAAnd (5) an image matrix. And determining each plasmid IDD、Fc、IAAAverage value of image matrix, denoted as IDD(C5V)、Fc(C5V)、IAA(C5V),IDD(C17V)、Fc(C17V)、IAA(C17V),IDD(C32V)、Fc(C32V)、IAA(C32V), and IDD(CTV)、Fc(CTV)、IAA(CTV). The point is then coordinated as
Figure BDA0003308966810000131
Figure BDA0003308966810000132
The four points are drawn in a plane rectangular coordinate system, and the G value obtained according to a two-hybrid-multi-plasmid method in the visual field is obtained by performing least square fitting on the four points. Averaging the G values obtained from 22 visual fields to obtain a correction factor G value of 3.72 +/-0.93 (FIG. 3);
(2) fitting the same data set using the AutoQT-FRET method of the invention (see steps 5.1 to 5.3 for a specific method) resulted in a mean correction factor of 3.51 ± 0.81 for 22 fields (fig. 3).
(3) According to the TP-G method (reference: Chen H, Puhl H L, Koushik S V, et al. measurement of FRET efficiency and ratio of don to acceptor con-centration in living cells [ J)]Biophysical journal,2006,91(5): L39-L41.) also gave values for the correction factor G. The TP-G method obtains the G value by a double-plasmid method, so that the G value can be calculated by combining two plasmids in a system of co-transforming four plasmids. We named the method of obtaining the correction factor G by C5V and C17V plasmids (in step 5.2.2, the four plasmids in the visual field of the cotransformation plasmid sample have been distinguished by C5V, C17V, C32V and CTV, so the correction factor G can be obtained by pairing two plasmids and calculating) as TP-G-1 method, the method of obtaining the correction factor G by C17V and C32V plasmids as TP-G-2 method, and the method of obtaining the correction factor G by C32V and CTV plasmids as TP-G-3 method. Similarly, the template image matrixes of the four plasmid data sets of C5V, C17V, C32V and CTV classified in the step 5.2.2 are respectively subjected to point multiplication by the template image matrixes of the four plasmid data sets obtained in the step 5.2, and subjected to background subtraction and crosstalk correction to obtain IDD、Fc、IAAImage matrix, obtaining I of each plasmidDD、Fc、IAAAnd (5) an image matrix. And determining each plasmid IDD、Fc、IAAAverage value of image matrix, denoted as IDD(C5V)、Fc(C5V)、IAA(C5V),IDD(C17V)、Fc(C17V)、IAA(C17V),IDD(C32V)、Fc(C32V)、IAA(C32V), and IDD(CTV)、Fc(CTV)、IAA(CTV). According to equation (12):
Figure BDA0003308966810000133
wherein 1 and 2 represent two different plasmids. Calculating the mean value of the correction factors to be 5.42 +/-1.81 by a method of obtaining 22 visual fields TP-G-1 according to the formula (12); the mean value of correction factors calculated by a TP-G-2 method is 4.23 +/-1.64; the mean correction factor calculated by the TP-G-3 method was 2.98. + -. 1.33 (FIG. 3).
The error of the correction factor G obtained by comparing the two-hybrid method and the AutoQT-FRET method is 5.64 percent and is within the error allowable range. Therefore, compared with the double-hybridization method, the correction factor obtained by the AutoQT-FRET method has stability. The time and workload for solving the correction factor G can be greatly reduced by using the AutoQT-FRET method.
5.5 the above experiment is to fit the spatial standard line of the visual field in any visual field, and finally obtain the correction factor of the visual field. In order to enable the obtained data of the correction factors to be more accurate, the data can be measured for multiple times and then averaged, for example, multiple visual fields (at least 10-15 visual fields) are selected, G and beta under different visual fields are obtained according to the methods of 5.2 and 5.3, and then the average values are obtained to obtain the correction factors G and beta of the automatic QT-FRET system.
The above embodiments are preferred embodiments of the present invention, but the present invention is not limited to the above embodiments, and any other changes, modifications, substitutions, combinations, and simplifications which do not depart from the spirit and principle of the present invention should be construed as equivalents thereof, and all such changes, modifications, substitutions, combinations, and simplifications are intended to be included in the scope of the present invention.

Claims (10)

1. An AutoQT-FRET method based on a correction factor of a primary imaging measurement system is characterized by comprising the following steps:
(1) sample preparation
Reference sample: respectively transfecting cells with the I FRET plasmids with different FRET efficiencies, and after the transfected plasmids are successfully expressed, respectively culturing the cells expressing the different FRET plasmids in different sample container carriers until the cells are attached to the walls, and marking as a reference sample 1 and a reference sample 2 … … reference sample i; wherein i is more than or equal to 2;
② the sample to be tested: respectively transfecting cells with the i FRET plasmids with different FRET efficiencies, merging the cells expressing the different FRET plasmids in the same sample container carrier after the transfected plasmids are successfully expressed, and culturing until the cells adhere to the wall to obtain a sample to be detected; wherein i is more than or equal to 2;
(2) E-FRET spectroscopic imaging
Thirdly, M visual fields are selected from the reference sample 1 obtained in the step I to carry out E-FRET imaging, and then M cells are selected under the 1 st visual field to obtain three-channel average fluorescence intensity IDD1-1、IDD1-2……IDD1-m,IAA1-1、IAA1-2……IAA1-m,IDA1-1、IDA1-2……IDA1-mAnd calculating the acceptor-sensitized fluorescence intensity FC1-1、FC1-2……FC1-m(ii) a Then with IDD1-1、IDD1-2……IDD1-mAs the abscissa, FC1-1、FC1-2……FC1-mA scatter plot is plotted for the ordinate and a least squares linear fit is performed to obtain the slope k at the 1 st field of view of the reference sample 11-1(ii) a By analogy, the slope k under the 2 nd to M th field of the reference sample 1 is obtained1-2、k1-3……k1-MAnd determining the maximum value and the minimum value to obtain the slope k of the reference sample 11Range of values: k is a radical ofMinimum value≤k1≤kMaximum value(ii) a Wherein M is more than or equal to 10, and M is more than or equal to 10;
replacing the reference sample 1 with the reference sample 2, and repeating the step c to obtain the slope k of the reference sample 22A range of values; and so on until the slope k of the reference sample i is obtainediA range of values;
selecting 1 visual field from the sample to be tested obtained in the step II to carry out E-FRET imaging to obtain three channels I thereofDD、IAA、IDAImage, and the length and width of the area of 100 x 100 pixels are respectively corresponding to IDD、IAA、IDAPerforming region segmentation on the image to obtain n cell region images IDD-S1、IDD-S2……IDD-Sn,IAA-S1、IAA-S2……IAA-Sn,IDA-S1、IDA-S2……IDA-Sn(ii) a Then, starting from the 1 st cell region, the proportion S of the effective data in the cell region is calculated1And measuring the three-channel average fluorescence intensity of all pixel points in the cell region, wherein the total number of the pixel points is L, and the pixel points are marked as IDD-1、IDD-2……IDD-L,IAA-1、IAA-2……IAA-L,IDA-1、IDA-2……IDA-LAnd simultaneously calculating the receptor-sensitized fluorescence intensity FC-1、FC-2……FC-LThen further calculating to obtain the ratio F of the 1 st cell areaC-1/IDD-1、FC-2/IDD-2……FC-L/IDD-LPerforming unimodal Gaussian fitting to obtain unimodal Gaussian function fitting degree value X of the cell region1And according to IDD-1、IDD-2……IDD-LAnd FC-1、FC-2……FC-LDrawing a scatter diagram, and then obtaining a slope k according to least square normal linear fitting1' value; by analogy, the proportion S of the effective data in the 2 nd cell area to the nth cell area is calculated in the same way2、S3……SnDegree of fit value X of single-peak Gaussian function2、X3……XnAnd slope k2′、k3′……kn'; wherein n is more than or equal to 25;
s degree1< 30% or X1(ii) < 90%, removing the cell region if S1Not less than 30% and X1The cell area is reserved if the cell area is more than or equal to 90 percent; in the same way, if S2< 30% or X2(ii) < 90%, removing the cell region if S2Not less than 30% and X2More than or equal to 90 percent, the cell area is reserved, and the like, R cell areas in total are obtained, and the corresponding slopes are marked as k1″、k2″……kR″;
Wherein R is more than or equal to 10;
(3) cell classification and determination of plasmid standard lines
Seventhly, obtaining k according to the step (c)1″、k2″……kR"size and range k of slope values obtained in Steps (c) and (d)1、k2……kiCell classification was performed: if k is1At k1In the range of (i), i.e. k1"corresponding cells are cells transfected with plasmid 1; if k is1At k2In the range of (i), i.e. k1"corresponding cells are cells … … f k transfected with plasmid 21At kiIn the range of (i), i.e. k1"corresponding cells are cells transfected with the i plasmid; by analogy, respectively combine k2″、k3″……kR"performing classification; then classifying the cells according to the I of each pixel point in all cell regions of each plasmidDDFc and IAAA spatial scattergram is drawn for the x, y and z axes, and a spatial straight line is fitted by the least square method to serve as a spatial standard line l for each plasmid1、l2……li
(4) Obtaining correction factors for E-FRET systems
Eighty to the step of obtaining the space standard line l of each plasmid1、l2……liPerforming plane least square fitting to obtain a plane equation Ax + By which is z; wherein A and B are constants, and B is a system correction factor beta1A/beta is the system correction factor G1
Ninthly, repeating the steps from five to eight, randomly selecting N different visual fields from the sample to be measured obtained in the step two to carry out E-FRET spectral imaging, and respectively calculating to obtain a system correction factor beta2、β3……βNAnd G2、G3……GNThen, taking the average value of the correction factors to obtain system correction factors G and beta; wherein N is more than or equal to 10.
2. The AutoQT-FRET method based on primary imaging measurement system correction factors according to claim 1, characterized in that:
the cell classification in step (3) can also be achieved by:
(a) selecting N under N' visual fields for the sample to be detected obtained in the step IIE-FRET imaging is carried out on the cell region of the' cell, and three channels of average fluorescence intensity I are obtained pixel by pixelDDˊ、IDAˊ、IAAAnd calculating the intensity of receptor-sensitized fluorescence F thereofCAnd with I of each pixelDDˊ、FCˊ、IAA' is a spatial coordinate point;
(b) respectively calculating the 1 st spatial coordinate point to the spatial standard line l obtained in the step (c)1、l2……liDistance d of1、d2……di(ii) a Then d is put1、d2……diSorting by magnitude of value, if d1And if the minimum value is reached, classifying the spatial coordinate point as a spatial standard line l1Namely, the cell region corresponding to the space coordinate point belongs to the cell transfected with the FRET plasmid of the 1 st type; in the same way, if d2And if the minimum value is reached, classifying the spatial coordinate point as a spatial standard line l2That is, the region of the cell corresponding to the spatial coordinate point belongs to the region d of the cell … … transfected with the FRET plasmid of the 2 nd speciesiAnd if the minimum value is reached, classifying the spatial coordinate point as a spatial standard line liNamely, the cell region corresponding to the space coordinate point belongs to the cells transfected with the i-th FRET plasmid; and by analogy, classifying the spatial coordinate point corresponding to each pixel to realize the classification of the cells.
3. The AutoQT-FRET method based on primary imaging measurement system correction factors according to claim 1, characterized in that:
the types of i in the step (1) are 2-4;
the number ratio of the FRET plasmid as a donor to an acceptor in the (1) part is 1: n or n: 1, n is more than or equal to 1;
the value range of M in the third step (2) is that M is more than or equal to 20;
the value range of m in the third step (2) is that m is more than or equal to 80;
the value range of R in the step (2) is more than or equal to 10 and less than or equal to 400;
the numerical range of N in the step (5) is as follows: n is more than or equal to 15.
4. The AutoQT-FRET method based on primary imaging measurement system correction factors according to claim 3, characterized in that:
the number ratio of the FRET plasmid as a donor to an acceptor in the (1) part is 1: 1, a FRET plasmid;
the value range of M in the third step (2) is that M is more than or equal to 40;
the value range of m in the third step (2) is that m is more than or equal to 250;
the numerical range of N in the step (5) is as follows: n is more than or equal to 40.
5. The AutoQT-FRET method based on primary imaging measurement system correction factors of claim 4, characterized in that:
the FRET plasmid in the step (1) is at least two of C5V, C17V, C32V and CTV;
the value range of M in the third step (2) is that M is more than or equal to 50;
the value range of m in the third step (2) is that m is more than or equal to 500.
6. The AutoQT-FRET method based on primary imaging measurement system correction factors according to claim 1, characterized in that:
the proportion of the effective data in the cell area in the fifth step (2) is calculated by the following method: in the cell region, the background is set to 0 and the non-background is set to 1, and then the proportion of the sum of non-background pixels to the sum of pixels of the entire region is calculated as the proportion of valid data.
7. The AutoQT-FRET method based on primary imaging measurement system correction factors according to claim 1, characterized in that:
the cells in the first step and the second step are Hela cells.
8. The AutoQT-FRET method based on primary imaging measurement system correction factors according to claim 1, characterized in that:
the selection of the m cells in the third step (2) is realized by the following mode: placing a sample container carrier filled with cells of transfection plasmids under a wide-field fluorescence microscope, firstly selecting M visual fields, then selecting cells from each visual field, and selecting M cells for measurement; wherein M is more than or equal to 10, and M is more than or equal to 10;
the number of the selected cells in each field of view is 1 or more.
9. The AutoQT-FRET method based on primary imaging measurement system correction factors according to claim 1, characterized in that:
the sample container carrier in the first step and the second step is a cell culture dish, a pore plate or a glass slide;
three channels I described in step (2)DD、IAA、IDAThe image is an image of 2048 × 2048 pixels, 1024 × 1024 pixels, or 512 × 512 pixels.
10. The application of the AutoQT-FRET method based on the primary imaging measurement system correction factor in the measurement system correction factor in any one of claims 1 to 9 is characterized in that:
the system correction factors are system correction factors G and beta.
CN202111211348.1A 2021-10-18 2021-10-18 AutoQT-FRET method based on primary imaging measurement system correction factor and application thereof Active CN113960001B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111211348.1A CN113960001B (en) 2021-10-18 2021-10-18 AutoQT-FRET method based on primary imaging measurement system correction factor and application thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111211348.1A CN113960001B (en) 2021-10-18 2021-10-18 AutoQT-FRET method based on primary imaging measurement system correction factor and application thereof

Publications (2)

Publication Number Publication Date
CN113960001A true CN113960001A (en) 2022-01-21
CN113960001B CN113960001B (en) 2022-11-01

Family

ID=79464301

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111211348.1A Active CN113960001B (en) 2021-10-18 2021-10-18 AutoQT-FRET method based on primary imaging measurement system correction factor and application thereof

Country Status (1)

Country Link
CN (1) CN113960001B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114900593A (en) * 2022-04-22 2022-08-12 华南师范大学 Camera mounting device for dual-channel FRET imaging and correction method thereof

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103212074A (en) * 2013-04-16 2013-07-24 中国科学院理化技术研究所 Preparation method of photosensitizer-doped nanoparticles with fluorescence resonance energy transfer and antenna effects and used for photodynamic therapy
CN105466902A (en) * 2016-01-18 2016-04-06 华南师范大学 Method for measuring fluorescence resonance energy transfer sensitized quenching conversion factors
EP3418348A1 (en) * 2017-06-21 2018-12-26 Université de Strasbourg Dye-loaded fluorescent polymeric nanoparticles as nano-antenna
CN112129737A (en) * 2020-09-24 2020-12-25 华南师范大学 Method for simultaneously and automatically measuring FRET system correction parameter and donor/acceptor extinction coefficient ratio and application thereof
CN113049555A (en) * 2021-03-12 2021-06-29 师大瑞利光电科技(清远)有限公司 Method for measuring linear separation and quantification of FRET system correction factor based on same system cell sample and application

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103212074A (en) * 2013-04-16 2013-07-24 中国科学院理化技术研究所 Preparation method of photosensitizer-doped nanoparticles with fluorescence resonance energy transfer and antenna effects and used for photodynamic therapy
CN105466902A (en) * 2016-01-18 2016-04-06 华南师范大学 Method for measuring fluorescence resonance energy transfer sensitized quenching conversion factors
EP3418348A1 (en) * 2017-06-21 2018-12-26 Université de Strasbourg Dye-loaded fluorescent polymeric nanoparticles as nano-antenna
US20200110032A1 (en) * 2017-06-21 2020-04-09 Universite De Strasbourg Dye-loaded fluorescent polymeric nanoparticles as nano-antenna
CN112129737A (en) * 2020-09-24 2020-12-25 华南师范大学 Method for simultaneously and automatically measuring FRET system correction parameter and donor/acceptor extinction coefficient ratio and application thereof
CN113049555A (en) * 2021-03-12 2021-06-29 师大瑞利光电科技(清远)有限公司 Method for measuring linear separation and quantification of FRET system correction factor based on same system cell sample and application

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
C. ZHANG ET AL.: "Automated E-FRET microscope for dynamical live-cell FRET imaging", 《JOURNAL OF MICROSCOPY》 *
魏丽春 等: "基于单波长激发的定量FRET成像", 《激光生物学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114900593A (en) * 2022-04-22 2022-08-12 华南师范大学 Camera mounting device for dual-channel FRET imaging and correction method thereof
CN114900593B (en) * 2022-04-22 2023-06-16 华南师范大学 Camera mounting device for dual-channel FRET imaging and correction method thereof

Also Published As

Publication number Publication date
CN113960001B (en) 2022-11-01

Similar Documents

Publication Publication Date Title
US20230145084A1 (en) Artificial immunohistochemical image systems and methods
Gaigalas et al. The development of fluorescence intensity standards
US9329130B2 (en) Systems and methods for counting cells and biomolecules
JP2023505265A (en) Systems and methods for high-throughput drug screening
CN101243192A (en) Method for detecting and quantitating multiple subcellular components
CN102165489A (en) Reproducible quantification of biomarker expression
EP2630492B1 (en) Internal focus reference beads for imaging cytometry
ES2359449T3 (en) METHODS TO IDENTIFY CONDITIONS THAT AFFECT A CELLULAR STATE.
CN112129737B (en) Method for simultaneously and automatically measuring FRET system correction parameter and donor/acceptor extinction coefficient ratio and application thereof
WO2008018905A2 (en) Method for predicting biological systems responses
EP2171443A1 (en) Method for predicting biological systems responses in hepatocytes
CN113049555B (en) Method for measuring linear separation and quantification of FRET system correction factor based on same system cell sample and application
CN113960001B (en) AutoQT-FRET method based on primary imaging measurement system correction factor and application thereof
JP2008518587A (en) Detection and quantification of subcellular components
US20030100024A1 (en) Method and apparatus for high throughput cell-based assays for screening and diagnostics
US20240125700A1 (en) Systems, methods and assays for outlier clustering unsupervised learning automated report (ocular)
CN108491688B (en) Method for preprocessing FRET (fluorescence resonance energy transfer) double-hybridization detection data based on donor-acceptor concentration ratio
CN104215561A (en) Method for accurately distinguishing cell cycle
US9122904B2 (en) Method for optimization of quantitative video-microscopy and associated system
JP2013526274A (en) Identification and / or characterization of microbial factors using taxonomic hierarchy classification
Siruvallur Murali et al. Evaluating Melanoma Viability and Proliferation in 3D Microenvironments
Keikhosravi et al. HiTIPS: High-Throughput Image Processing Software for the Study of Nuclear Architecture and Gene Expression
CN111596053A (en) Application of TPN (thermoplastic vulcanizate) molecule in preparation of circulating tumor cell detection reagent, detection reagent and kit
Kotrys et al. High-throughput measurement of mitochondrial RNA turnover in human cultured cells
US11908130B2 (en) Apparatuses and methods for digital pathology

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant