CN106591443A - Method for analytical processing of crop fluorescence capillary electrophoretogram - Google Patents

Method for analytical processing of crop fluorescence capillary electrophoretogram Download PDF

Info

Publication number
CN106591443A
CN106591443A CN201611102078.XA CN201611102078A CN106591443A CN 106591443 A CN106591443 A CN 106591443A CN 201611102078 A CN201611102078 A CN 201611102078A CN 106591443 A CN106591443 A CN 106591443A
Authority
CN
China
Prior art keywords
peak
pull
iii
pulls
coincidence
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
CN201611102078.XA
Other languages
Chinese (zh)
Other versions
CN106591443B (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.)
Beijing Academy of Agricultural Sciences Seed Industry Technology Co.,Ltd.
Original Assignee
BEIJING HUASHENG HENGYE TECHNOLOGY Co Ltd
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 BEIJING HUASHENG HENGYE TECHNOLOGY Co Ltd filed Critical BEIJING HUASHENG HENGYE TECHNOLOGY Co Ltd
Priority to CN201611102078.XA priority Critical patent/CN106591443B/en
Publication of CN106591443A publication Critical patent/CN106591443A/en
Application granted granted Critical
Publication of CN106591443B publication Critical patent/CN106591443B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6844Nucleic acid amplification reactions
    • C12Q1/6858Allele-specific amplification
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B45/00ICT specially adapted for bioinformatics-related data visualisation, e.g. displaying of maps or networks

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Organic Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Wood Science & Technology (AREA)
  • Biophysics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Zoology (AREA)
  • General Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biotechnology (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Biology (AREA)
  • Immunology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Analytical Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Microbiology (AREA)
  • Molecular Biology (AREA)
  • Medical Informatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Theoretical Computer Science (AREA)
  • Biochemistry (AREA)
  • General Engineering & Computer Science (AREA)
  • Genetics & Genomics (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

The invention relates to a method and device for analytical processing of a crop fluorescence capillary electrophoretogram, wherein the method includes 3 major processing steps of curve processing, internal standard calibration, and fragment analysis and genotyping, and is characterized in that curve processing includes pull-up peak processing, the pull-up peak processing includes the steps of overlapped peak positioning, overlapped peak grouping, overlapped peak processing, quasi pull-up peak grouping, quasi pull-up peak re-grouping, pull-up peak III judging and pull-up peak III processing. The device includes three modules of curve processing, internal standard calibration, and fragment analysis and genotyping and is characterized in that the curve processing module includes a pull-up peak processing sub module. The problem of internal standard deviation caused because a peak at an internal standard peak position is deleted can be effectively avoided, the sequencing accuracy rate of the crop fluorescence capillary electrophoretogram is improved significantly, and multiple laboratories are facilitated to work together to construct a crop genetic information fingerprint database.

Description

A kind of method of analyzing and processing crops fluorescent capillary electrophoresis tube collection of illustrative plates
Technical field
This invention relates generally to biogenetics and genomics field, and more specifically it relates at a kind of analysis The method and its device of reason crops fluorescent capillary electrophoresis tube collection of illustrative plates.
Background technology
Fluorescent capillary electrophoresis tube is the conventional method of Treatment Analysis bio-genetic material.The genetic fragments of different length can be with It is separated by electrophoresis, and the genetic fragments of different but identical length of originating then can be made a distinction by different fluorescence colors. When fluorescent dye primer is used, the presence or absence of pcr amplification product or how much by fluorescence intensity judging, but it is different There is overlap between the fluorescing frequencies of color.If the overflow of a certain color, has exceeded the process range of program, or Also there is the small peak of another kind of color in power of test of the person beyond instrument, the position that superelevation peak will occur in a kind of color, This small peak is and pulls up (pull up) peak.It is this as color is diffused into " pulling up existing for another from spectrum channel As " it is widely present in fluorescence electrophoresis experimental situation.
Relative to the research of mankind's inheritance resources, the research of crops genetic resourcess falls behind relatively.Currently used for crops The design of primers of genetic marker is relatively deficient, the test kit for also not having maturation, and this causes Jing PCR amplifications and fluorescent capillary electrophoresis tube There is the more peak that pulls up in the crops gene test collection of illustrative plates of generation, and these pull up peak peak height it is higher, it is not easy to other Distinguish at normal peak.Also, due to reference material produce internal standard colour band peak height it is relatively low, if now sample produce common color Band has the peak for overlapping with internal standard colour band, and normal internal standard base peak will likely be taken as pulling up peak and being deleted, so as to meeting Cause internal standard calibration poor or failure.Such situation is pulled up in crops can be tighter under the higher environment of the universal peak height in peak Weight.In addition, the gene test collection of illustrative plates produced after crops electrophoresis easily produces continuous multimodal, i.e., a series of equidistant, spacing are not More than 2.5bp and peak number more than 2 peak, if pulling up peak occurs in continuous multimodal region, it is likely that continuous multimodal can be taken as By and be ignored, cause more to read peaks or misread the situation at peak.
It is existing process fluorescent capillary electrophoresis tube collection of illustrative plates in pull up peak method be generally only by adjacent bands collection of illustrative plates contrast with Examination goes out to pull up peak, and the peak that pulls up that examination is gone out simply is deleted.But this method cannot effective process farming analyte detection collection of illustrative plates In have with internal standard colour band and overlap or pull up peak in continuous many crest segments, often so that internal standard calibration is poor in addition failure or Cause many reading peaks or misread the situation at peak, failure is even detected so as to cause testing result relatively large deviation occur.Therefore, how Peak can will be pulled up in the case where internal standard calibration and loss peak is not affected and do effective process, so as to improve crops fluorescent capillary pipe Electrophoresis pattern gene test accuracy is the key technical problem that those skilled in the art need to solve.
The content of the invention
Asking for peak cannot be pulled up in effective process crops fluorescent labeling detection collection of illustrative plates to solve to exist in prior art Topic, it is it is an object of the invention to provide a kind of method that analyzing and processing crops fluorescent labeling detects collection of illustrative plates including curve processing, interior Calibration standard and 3 main processing steps of fragment analysis and gene type, it is characterised in that also include pulling up peak in curve processing Process, it is described pull up peak process comprise the following steps:
1) coincidence peak positioning:A colour bands are contrasted with B colour bands, in finding B colour bands, position is overlapped with the position at A colour bands peak Peak, these peaks are set to into coincidence peak, wherein A colors have Chong Die and A colour bands that B colour bands can be caused to pull up with B color fluorescence spectruies;
2) coincidence peak packet:
When step 1) position at coincidence peak that positions when being overlapped with the position of internal standard base peak, these coincidence peaks is set to and pulled out Peak I is played,
When coincidence peak is located at continuous many crest segments, these coincidence peaks for being located at continuous many crest segments is set to into standard and pulls up peak;
When the position at coincidence peak is not neither overlapped with internal standard base peak in continuous many crest segments again, these coincidence peaks are set to Peak II is pulled up,
3) coincidence peak is processed:
Record each position for pulling up peak I peak I is pulled up so that follow-up internal standard calibration steps is referred to and eliminated,
Elimination pulls up peak II,
Record standard pulls up the position X that total number N and each standard at peak pulls up peaki B, wherein, i is no less than 0 and is not more than The integer of N,
When N is 0, end pulls up peak process step;
4) standard pulls up peak packet:It is X by positioni BPeak in all B colour bands of ± 32frame is divided into i-th group;
5) standard pulls up peak and is grouped again:N number of group in step 4 is grouped again, that is, is found per the summit in group, will most Peak and its two peaks of left and right are divided into subgroup 1i of this group, and other peaks in group are divided into another subgroup 2i;
6) peak III judgements are pulled up:In each ai subgroup that step 5 determines,
A) peak height at A colour bands peak in the range of subgroup position is sorted and is set to the position of the summit in the range of this Xai A, the position of summit in subgroup is set to into X in additionai B,
B) by Xai BWith Xai BContrast, if Xai B≠Xai B, all peaks in subgroup ai are all judged as not being to pull up peak III, Otherwise next step judgement is then carried out,
C) when position is Xai BPeak peak height Hai BIt is X with positionai APeak peak height Hai ABetween ratio be more than c when, Position is Xai BPeak be judged to be defined and pull up peak III, the wherein value of c is adjusted according to tested species,
D) sort from high to low according to peak height, judge that other peaks are in subgroup ai according to the height ratio of the position and peak at peak No being defined pulls up peak III,
E) when the ratio of the number at all peaks in the number and subgroup ai that standard in subgroup ai pulls up peak III is less than 0.75 When, all peaks in subgroup ai are all judged as not being to pull up peak III, otherwise are all then to pull up peak III,
Wherein a is subgroup number, and a values are 1 or 2;
7) pull up peak III process:Peak III eliminations will be pulled up, and terminates behind to pull up peak process step.
Further, step 7 pull up peak III process described in elimination be:The adjacent situations of peak III are pulled up when existing When, adjacent two pull up peak III and adopt curve difference approximating method to be combined process and pull up peak III to eliminate, otherwise then Carry out unimodal process peak III is pulled up to eliminate.Can be made Jing pulling up after peak III Processing for removing using curve difference approximating method Picture is more smoothed, it is to avoid the figures such as spine is occurred and is produced impact to follow-up atlas analysis.
Further, the curve difference approximating method is linear interpolation or cubic spline interpolation.Imitate in same process Under fruit, linear interpolation and cubic spline interpolation graphics calculations amount are less, and graphics process speed is faster.
Further, the curve processing also include nail peak process step, nail peak process step be by due to During electrophoresis, electric current is unstable and caused nail-like spike is removed.
Further, the curve processing also include saturation peak process step, saturation peak process step be by due to Peak intensity has exceeded the threshold value of sequenator and caused tack or concave head type saturation peak and the pulling out of occurring of correspondence under saturation peak Play peak superimposed so as to form a complete peak.
It is another object of the present invention to be to provide the device that a kind of analyzing and processing crops fluorescent labeling detects collection of illustrative plates, institute Stating device includes curve processing, internal standard calibration and 3 modules of fragment analysis and gene type, it is characterised in that curve processing Also include in module that pulling up peak processes submodule, the peak processing module that pulls up is included with lower unit:
Coincidence peak positioning unit, for A colour bands are contrasted with B colour bands, finds the position at position and A colour bands peak in B colour bands The peak of coincidence is put, these peaks are set to into coincidence peak, wherein A colors have Chong Die with B color fluorescence spectruies and A colour bands can cause B colour bands Pull up;
Coincidence peak grouped element, for being grouped to coincidence peak, specially:
When the position for being set to the coincidence peak positioned by unit is overlapped with internal standard standard peak position, these are overlapped by this unit Peak is set to and pulls up peak I,
When the position at coincidence peak is not neither overlapped with internal standard base peak in continuous many crest segments again, these are overlapped by this unit Peak is set to and pulls up peak II,
When coincidence peak is located at continuous many crest segments, these coincidence peaks for being located at continuous many crest segments are set to standard and are pulled up by this unit Peak;
Coincidence peak processing unit, for processing the Chong Hefeng after packet respectively, specially:
Record each position for pulling up peak I peak I is pulled up so that follow-up internal standard calibration module is referred to and eliminated,
Elimination pulls up peak II,
Record standard pulls up the position X that total number N and each standard at peak pulls up peaki B, wherein, i is no less than 0 and is not more than The integer of N,
When N is 0, end pulls up peak process;
Standard pulls up peak grouped element, pulls up peak for alignment and is grouped, specially:It is X by positioni B± 32frame's Peak in all B colour bands is divided into i-th group;
Standard pulls up peak grouped element again, pulls up peak for alignment and is further grouped, specially:It is single that standard pulls up peak packet N number of group after unit is processed is grouped again, that is, find per the summit in group, and summit and its two peaks of left and right are divided into this group Subgroup 1i, other peaks in group are divided into another subgroup 2i;
Peak III judging units are pulled up, and the judgement at peak are pulled up for alignment, specially:Standard pulls up peak, and grouped element determines again Each ai subgroup in,
A) peak height at the peak of A colour bands in the range of subgroup position is sorted and is set to the position of the summit in the range of this Xai A, the position of summit in subgroup is set to into X in additionai B,
B) by Xai BWith Xai BContrast, if Xai B≠Xai B, all peaks in subgroup ai are all judged as not being to pull up peak III, Otherwise next step judgement is then carried out,
C) when position is Xai BPeak peak height Hai BIt is X with positionai APeak peak height Hai ABetween ratio be more than c when, Position is Xai BPeak be judged as pulling up peak III, the wherein value of c is adjusted according to tested species,
D) judge whether other peaks are to pull up peak III in subgroup ai according to the height of the position and peak at peak,
E) when the ratio of the number at all peaks in number and subgroup ai that peak III is pulled up in subgroup ai is less than 0.75, All peaks in subgroup ai are all judged as not being to pull up peak III, otherwise are all then to pull up peak III,
Wherein a is subgroup number, and a values are 1 or 2;
Peak III processing units are pulled up, for peak III eliminations will be pulled up, and terminates behind to pull up peak process.
Further, the curve processing module also includes that nail peak processes submodule, and the nail peak processes submodule For will due to during electrophoresis electric current it is unstable and caused nail-like spike is removed.
Further, the curve processing module also includes that saturation peak processes submodule, and the saturation peak processes submodule For will exceed the threshold value of sequenator due to peak intensity and caused tack or concave head type saturation peak and under saturation peak it is right Should occur to pull up peak superimposed so as to form a complete peak.
It is a further object of the present invention to provide a kind of electric terminal, it is characterised in that the electric terminal includes above-mentioned Any one analyzing and processing fluorescent labeling detects the device of collection of illustrative plates.
Further, the electric terminal also includes:
Communication module, for by way of wiredly and/or wirelessly by the analyzing and processing crops SSR marker collection of illustrative plates Collection of illustrative plates or the corresponding data transfer of collection of illustrative plates after device process is to outside other equipment.
Curve processing of the present invention may also include the means such as baseline calibration, curve smoothing process.
Due to the interference that there is fluorescent material in electrophoresis system, the baseline of each colour band in collection of illustrative plates there may be deviation. Baseline calibration of the present invention is that the minimum peak of each colour band is set to 0, makes each colour band on same datum line.
Smoothing processing is to be smoothed zigzag fashion present in collection of illustrative plates, becomes apparent from original peaks readable.
Due to the interference that there is operating error and other materials in PCR amplifications and electrophoretic process, SSR collection of illustrative plates there may be The situation of horizontal displacement.Internal standard calibration steps of the present invention is will to detect internal standard substance in obtained collection of illustrative plates (also known as reference material) The position at corresponding peak is compared with the position at each peak in internal standard substance standard diagram, to calibrate the abscissa of collection of illustrative plates, that is, is migrated Distance (frame).
Fragment analysis (GeneScan) of the present invention are that the internal standard substance of known length is done molecular weight standard first to obtain Standard curve of the mrna length (bp) to migration distance (frame), the peak of the unknown sample of detection and standard curve are compared, Calculate the corresponding actual gene fragment length (bp) in each peak.
Gene type (Genotyping) of the present invention is by the peak and allelic ladder of unknown sample (allelic ladder) compares, to confirm the corresponding allele in peak.
The peak number mentioned in the present invention refers to the number at peak.
Method compared with prior art excellent of analyzing and processing crops fluorescent capillary electrophoresis tube collection of illustrative plates of the present invention Put and be:As the peak that pulls up to diverse location is made and being handled differently, peak can be pulled up in effectively elimination and spectrum recognition is caused Impact while can also be prevented effectively from continuous multimodal be directly deleted so as to cause hereditary information omit problem, can also have Effect avoids the internal standard calibration offset issue for being deleted and being caused due to the peak of internal standard standard peak position, makes crops fluorescent capillary pipe The sequencing accuracy rate of electrophoresis pattern is obviously improved, and is conducive to many laboratorys to cooperate jointly and is built crops hereditary information fingerprint Storehouse.
Description of the drawings
Fig. 1 is the part detection collection of illustrative plates contrast of an embodiment of the method through the present invention and a method comparative example process Figure, wherein upper half figure are the part detection collection of illustrative plates processed through the comparative example, and lower half figure is the same part that the embodiment is processed Detection collection of illustrative plates;
Fig. 2 is the functional block diagram of a device embodiment of the present invention;
Fig. 3 is the functional block diagram of another device embodiment of the present invention.
Specific embodiment
Embodiment of the method 1
Be sequenced using the SSR marker of P40 primer pair Semen Maydiss, and corresponding detection is formed using PCR amplifications and fluorescent labeling Collection of illustrative plates.
Process is analyzed to the detection collection of illustrative plates, the analyzing and processing includes curve processing, internal standard calibration and fragment point Analysis and 3 main processing steps of gene type, it is characterised in that also include pulling up peak process in curve processing, it is described to pull up peak Process is comprised the following steps:
1) coincidence peak positioning:A colour bands are contrasted with B colour bands, in finding B colour bands, position is overlapped with the position at A colour bands peak Peak, these peaks are set to into coincidence peak, wherein A colors have Chong Die and A colour bands that B colour bands can be caused to pull up with B color fluorescence spectruies;
2) coincidence peak packet:
When step 1) position at coincidence peak that positions when being overlapped with the position of internal standard base peak, these coincidence peaks is set to and pulled out Peak I is played,
When coincidence peak is located at continuous many crest segments, these coincidence peaks for being located at continuous many crest segments is set to into standard and pulls up peak;
When the position at coincidence peak is not neither overlapped with internal standard base peak in continuous many crest segments again, these coincidence peaks are set to Peak II is pulled up,
3) coincidence peak is processed:
Record each position for pulling up peak I peak I is pulled up so that follow-up internal standard calibration steps is referred to and eliminated,
Elimination pulls up peak II,
Record standard pulls up the position X that total number N and each standard at peak pulls up peaki B, wherein, i is no less than 0 and is not more than The integer of N,
When N is 0, end pulls up peak process step;
4) standard pulls up peak packet:It is X by positioni BPeak in all B colour bands of ± 32frame is divided into i-th group;
5) standard pulls up peak and is grouped again:N number of group in step 4 is grouped again, that is, is found per the summit in group, will most Peak and its two peaks of left and right are divided into subgroup 1i of this group, and other peaks in group are divided into another subgroup 2i;
6) peak III judgements are pulled up:In each ai subgroup that step 5 determines,
A) peak height at A colour bands peak in the range of subgroup position is sorted and is set to the position of the summit in the range of this Xai A, the position of summit in subgroup is set to into X in additionai B,
B) by Xai BWith Xai BContrast, if Xai B≠Xai B, all peaks in subgroup ai are all judged as not being to pull up peak III, Otherwise next step judgement is then carried out,
C) when position is Xai BPeak peak height Hai BIt is X with positionai APeak peak height Hai ABetween ratio be more than c when, Position is Xai BPeak be judged to be defined and pull up peak III, the wherein value of c is adjusted according to tested species,
D) sort from high to low according to peak height, judge that other peaks are in subgroup ai according to the height ratio of the position and peak at peak No being defined pulls up peak III,
E) when the ratio of the number at all peaks in the number and subgroup ai that standard in subgroup ai pulls up peak III is less than 0.75 When, all peaks in subgroup ai are all judged as not being to pull up peak III, otherwise are all then to pull up peak III,
Wherein a is subgroup number, and a values are 1 or 2;
7) pull up peak III process:Peak III eliminations will be pulled up, and terminates behind to pull up peak process step.
Detection collection of illustrative plates after process see accompanying drawing 1 (under).
Method comparative example 1
Be sequenced using the SSR marker of P40 primer pair Semen Maydiss, and corresponding detection is formed using PCR amplifications and fluorescent labeling Collection of illustrative plates.
Be analyzed process to the detection collection of illustrative plates, the analyzing and processing in addition to not pulling up peak and processing, remaining and enforcement Treatment Analysis described in example 1 are identical.
Detection collection of illustrative plates after process see accompanying drawing 1 (on).
Can be seen by accompanying drawing 1, carrying out pulling up after peak processes, B colour bands (lower section colour band) 298frame nearby with 310frame is nearby eliminated by the peak that pulls up that the continuous multimodal of A colour bands takes up, in addition, the peak of the corresponding continuous multimodal of A colour bands Height is corrected, and is conducive to follow-up other image procossings.
Experimental example 1
Sequencing 3 times are repeated with sample sequencing point to 200 known arrays of corn and soybean, Semen Tritici aestivi, it is same to every kind of crops The sequencing procedure of one known array is identical.Sequencing uses SSR marker, PCR amplifications and fluorescent labeling to form detection collection of illustrative plates, Mei Zhongnong Crop respectively obtains 300 detection collection of illustrative plates.By these collection of illustrative plates respectively through process (i.e. similar to Example 1) and discontinuous multimodal Processing method (i.e. similar to comparative example 1) is analyzed process, and the result after analyzing and processing and known array are contrasted, and obtains To the average detected error rate of each method, concrete outcome see the table below.
After the process of continuous multimodal, detect that the average detected error rate of collection of illustrative plates is substantially reduced, substantially improve farming The sequencing accuracy of thing fluorescent capillary electrophoresis tube collection of illustrative plates, is conducive to many laboratorys to cooperate jointly and builds crops hereditary information fingerprint Storehouse.
Fig. 2 is the functional block diagram of a device embodiment of the present invention.As shown in Fig. 2 the device 100 includes:
Including curve processing module 110, internal standard calibration module 120, fragment analysis and gene type module 130, at curve Also include in reason that pulling up peak processes submodule 111, the peak process submodule that pulls up is included with lower unit:
Coincidence peak positioning unit 1111, for A colour bands are contrasted with B colour bands, finds position and A colour bands peak in B colour bands Position overlap peak;
Coincidence peak grouped element 1112, for being grouped to coincidence peak;
Coincidence peak processing unit 1113, for processing the Chong Hefeng after packet respectively;
Standard pulls up peak grouped element 1114, pulls up peak for alignment and is grouped;
Standard pulls up peak grouped element 1115 again, pulls up peak for alignment and is further grouped;
Peak III judging unit 1116 is pulled up, and the judgement at peak is pulled up for alignment;
Peak III processing unit 1117 is pulled up, for peak III eliminations will be pulled up.
Fig. 3 is the functional block diagram of a device embodiment of the present invention.As shown in figure 3, the curve processing in the device 100 Module 110 can also include that nail peak processes submodule 112, for will as during electrophoresis, electric current is unstable caused nail-like point Peak is removed;And saturation peak processes submodule 113, for the threshold value of sequenator will have been exceeded due to peak intensity and caused tack Or concave head type saturation peak and under saturation peak correspondence occur to pull up peak superimposed so as to form a complete peak.
The method of work of the device of the embodiment of the present invention please refer to above method embodiment, therefore will not be described in detail herein.
The embodiment of the present invention also provides a kind of electric terminal, and the electric terminal includes that the analyzing and processing crops of the above are glimmering The device of light capillary electrophoresis collection of illustrative plates.Further, the electric terminal can also include:Communication unit, for by wired And/or the collection of illustrative plates or collection of illustrative plates after wirelessly processed the device of the analyzing and processing crops fluorescent capillary electrophoresis tube collection of illustrative plates Corresponding data transfer is to outside other equipment.
Through the above description of the embodiments, those skilled in the art can be understood that the present invention can be by Software adds the mode of required general hardware platform to realize, naturally it is also possible to by hardware, or the combination of the two is implementing. Based on such understanding, the part that technical scheme is substantially contributed to prior art in other words can be with software The form of product is embodied, and the software module or computer software product can be stored in a storage medium, if including Dry instruction is used so that a computer equipment (can be personal computer, server, or network equipment etc.) performs this Method described in bright each embodiment.Storage medium can be random access memory (RAM), internal memory, read only memory (ROM), electricity Well known in programming ROM, electrically erasable ROM, depositor, hard disk, moveable magnetic disc, CD-ROM or technical field Any other form of storage medium.
Although the principle of the present invention has had specific descriptions herein, it will be appreciated by those of skill in the art that this Individual description only by way of example carrying out, and not as the restriction with regard to the scope of the present invention.Except shown herein Outside the exemplary embodiment of description, within the scope of the invention by the modification and replacement of one of ordinary skill in the art Other embodiment, be regarded as within the scope of the invention.

Claims (10)

1. a kind of method that analyzing and processing crops fluorescent labeling detects collection of illustrative plates, including the calibration of curve processing, internal standard and fragment Analysis and 3 main processing steps of gene type, it is characterised in that also include pulling up peak process in curve processing, it is described to pull up Peak process is comprised the following steps:
1) coincidence peak positioning:A colour bands are contrasted with B colour bands, the peak that position in B colour bands is overlapped with the position at A colour bands peak is found, These peaks are set to into coincidence peak, wherein A colors there are Chong Die and A colour bands that B colour bands can be caused to pull up with B color fluorescence spectruies;
2) coincidence peak packet:
When step 1) position at coincidence peak that positions when being overlapped with the position of internal standard base peak, these coincidence peaks is set to and pulls up peak I,
When coincidence peak is located at continuous many crest segments, these coincidence peaks for being located at continuous many crest segments is set to into standard and pulls up peak;
When the position at coincidence peak is not neither overlapped with internal standard base peak in continuous many crest segments again, these coincidence peaks are set to and are pulled up Peak II,
3) coincidence peak is processed:
Record each position for pulling up peak I peak I is pulled up so that follow-up internal standard calibration steps is referred to and eliminated,
Elimination pulls up peak II,
Record standard pulls up the position X that total number N and each standard at peak pulls up peaki B, wherein, i is no less than 0 and no more than N Integer,
When N is 0, end pulls up peak process step;
4) standard pulls up peak packet:It is X by positioni BPeak in all B colour bands of ± 32frame is divided into i-th group;
5) standard pulls up peak and is grouped again:N number of group in step 4 is grouped again, that is, is found per the summit in group, by summit And its two peaks of left and right are divided into subgroup 1i of this group, other peaks in group are divided into another subgroup 2i;
6) peak III judgements are pulled up:In each ai subgroup that step 5 determines,
A) peak height at A colour bands peak in the range of subgroup position is sorted and the position of the summit in the range of this is set to into Xai A, separately The outer position by summit in subgroup is set to Xai B,
B) by Xai BWith Xai BContrast, if Xai B≠Xai B, all peaks in subgroup ai are all judged as not being to pull up peak III, otherwise Next step judgement is carried out then,
C) when position is Xai BPeak peak height Hai BIt is X with positionai APeak peak height Hai ABetween ratio be more than c when, position For Xai BPeak be judged to be defined and pull up peak III, the wherein value of c is adjusted according to tested species,
D) sort from high to low according to peak height, judge in subgroup ai that whether other peaks are according to the height ratio of the position and peak at peak Standard pulls up peak III,
E) when the ratio of the number at all peaks in number and subgroup ai that standard in subgroup ai pulls up peak III is less than 0.75, son Group ai in all peaks be all judged as not being to pull up peak III, otherwise be all then to pull up peak III,
Wherein a is subgroup number, and a values are 1 or 2;
7) pull up peak III process:Peak III eliminations will be pulled up, and terminates behind to pull up peak process step.
2. method according to claim 1, it is characterised in that step 7 pull up peak III process described in elimination be:When When presence pulls up peak III adjacent situations, adjacent two pull up peak III and are combined process using curve difference approximating method Peak III is pulled up to eliminate, otherwise is then carried out unimodal process and peak III is pulled up to eliminate.
3. method according to claim 2, it is characterised in that:The curve difference approximating method is linear interpolation or three times Spline interpolation.
4. the method according to any one of claim 1-3, it is characterised in that:The curve processing is also included at nail peak Reason step, nail peak process step are by due to electric current shakiness during electrophoresis, caused nail-like spike is removed.
5. the method according to any one of claim 1-4, it is characterised in that:The curve processing is also included at saturation peak Reason step, saturation peak process step are will to have exceeded the threshold value of sequenator due to peak intensity and caused tack or concave head type Saturation peak and under saturation peak correspondence occur to pull up peak superimposed so as to form a complete peak.
6. a kind of analyzing and processing crops fluorescent labeling detects the device of collection of illustrative plates, and described device includes curve processing, internal standard calibration And 3 modules of fragment analysis and gene type, it is characterised in that also include in curve processing module that pulling up peak processes submodule Block, the peak processing module that pulls up are included with lower unit:
Coincidence peak positioning unit, for A colour bands are contrasted with B colour bands, in finding B colour bands, position is heavy with the position at A colour bands peak These peaks are set to coincidence peak by the peak of conjunction, and wherein A colors have Chong Die and A colour bands that B colour bands can be caused to pull up with B color fluorescence spectruies;
Coincidence peak grouped element, for being grouped to coincidence peak, specially:
When the position for being set to the coincidence peak positioned by unit is overlapped with internal standard standard peak position, this unit is fixed by these coincidence peaks To pull up peak I,
When the position at coincidence peak is not neither overlapped with internal standard base peak in continuous many crest segments again, this unit is fixed by these coincidence peaks To pull up peak II,
When coincidence peak is located at continuous many crest segments, these coincidence peaks for being located at continuous many crest segments are set to standard and pull up peak by this unit;
Coincidence peak processing unit, for processing the Chong Hefeng after packet respectively, specially:
Record each position for pulling up peak I peak I is pulled up so that follow-up internal standard calibration module is referred to and eliminated,
Elimination pulls up peak II,
Record standard pulls up the position X that total number N and each standard at peak pulls up peaki B, wherein, i is no less than 0 and no more than N Integer,
When N is 0, end pulls up peak process;
Standard pulls up peak grouped element, pulls up peak for alignment and is grouped, specially:It is X by positioni BAll B of ± 32frame Peak in colour band is divided into i-th group;
Standard pulls up peak grouped element again, pulls up peak for alignment and is further grouped, specially:Standard is pulled up at peak grouped element N number of group after reason is grouped again, that is, find per the summit in group, and summit and its two peaks of left and right are divided into the one of this group Individual subgroup 1i, other peaks in group are divided into another subgroup 2i;
Peak III judging units are pulled up, and the judgement at peak are pulled up for alignment, specially:Standard pull up peak again grouped element determine it is every In individual ai subgroups,
A) peak height at the peak of A colour bands in the range of subgroup position is sorted and the position of the summit in the range of this is set to into Xai A, The position of summit in subgroup is set to into X in additionai B,
B) by Xai BWith Xai BContrast, if Xai B≠Xai B, all peaks in subgroup ai are all judged as not being to pull up peak III, otherwise Next step judgement is carried out then,
C) when position is Xai BPeak peak height Hai BIt is X with positionai APeak peak height Hai ABetween ratio be more than c when, position For Xai BPeak be judged as pulling up peak III, the wherein value of c is adjusted according to tested species,
D) judge whether other peaks are to pull up peak III in subgroup ai according to the height of the position and peak at peak,
E) when the ratio of the number at all peaks in number and subgroup ai that peak III is pulled up in subgroup ai is less than 0.75, subgroup All peaks in ai are all judged as not being to pull up peak III, otherwise are all then to pull up peak III,
Wherein a is subgroup number, and a values are 1 or 2;
Peak III processing units are pulled up, for peak III eliminations will be pulled up, and terminates behind to pull up peak process.
7. device according to claim 6, it is characterised in that:The curve processing module also includes that nail peak processes submodule Block, the nail peak process submodule is used for will the caused nail-like spike removal due to electric current shakiness during electrophoresis.
8. device according to claim 6, it is characterised in that:The curve processing module also includes that saturation peak processes submodule Block, the saturation peak process submodule for having exceeded the threshold value of sequenator due to peak intensity and caused tack or concave head type Saturation peak and under saturation peak correspondence occur to pull up peak superimposed so as to form a complete peak.
9. a kind of electric terminal, it is characterised in that the electric terminal is included at the analysis any one of claim 4-6 Reason crops fluorescent labeling detects the device of collection of illustrative plates.
10. electric terminal according to claim 7, it is characterised in that the electric terminal also includes:
Communication module, for the analyzing and processing crops fluorescent labeling is detected collection of illustrative plates by way of wiredly and/or wirelessly Device process after collection of illustrative plates or the corresponding data transfer of collection of illustrative plates to outside other equipment.
CN201611102078.XA 2016-12-05 2016-12-05 A kind of method for analyzing and processing crops fluorescent capillary electrophoresis tube collection of illustrative plates Active CN106591443B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611102078.XA CN106591443B (en) 2016-12-05 2016-12-05 A kind of method for analyzing and processing crops fluorescent capillary electrophoresis tube collection of illustrative plates

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611102078.XA CN106591443B (en) 2016-12-05 2016-12-05 A kind of method for analyzing and processing crops fluorescent capillary electrophoresis tube collection of illustrative plates

Publications (2)

Publication Number Publication Date
CN106591443A true CN106591443A (en) 2017-04-26
CN106591443B CN106591443B (en) 2017-11-28

Family

ID=58594832

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611102078.XA Active CN106591443B (en) 2016-12-05 2016-12-05 A kind of method for analyzing and processing crops fluorescent capillary electrophoresis tube collection of illustrative plates

Country Status (1)

Country Link
CN (1) CN106591443B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102880812A (en) * 2012-08-02 2013-01-16 北京华生恒业科技有限公司 Method and device for processing solubility curve based on wavelet analysis method, and electronic terminal

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102880812A (en) * 2012-08-02 2013-01-16 北京华生恒业科技有限公司 Method and device for processing solubility curve based on wavelet analysis method, and electronic terminal

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
MITCHELL M. HOLLAND 等: "《GeneMarker_ HID: A Reliable Software Tool for the Analysis of Forensic STR Data》", 《J FORENSIC SCI》 *
孙亚男 等: "《SNaPshot技术检测Y-SNP位点O3-M122单倍群的应用》", 《中国法医学杂志》 *
沈红缨 等: "《DNA数据库样本数据分析中常见问题探讨》", 《中国法医学杂志》 *

Also Published As

Publication number Publication date
CN106591443B (en) 2017-11-28

Similar Documents

Publication Publication Date Title
Nelson et al. Whole-genome validation of high-information-content fingerprinting
US6681186B1 (en) System and method for improving the accuracy of DNA sequencing and error probability estimation through application of a mathematical model to the analysis of electropherograms
CN104302781B (en) A kind of method and device detecting chromosomal structural abnormality
Gill et al. Automated short tandem repeat (STR) analysis in forensic casework—a strategy for the future
US20210343367A1 (en) Methods for detecting mutation load from a tumor sample
US20040197925A1 (en) Analysis
EP4004236B1 (en) Method for analysing loss-of-heterozygosity (loh) following deterministic restriction-site whole genome amplification (drs-wga)
JP2005500542A5 (en)
CN105349617A (en) High-throughput RNA sequencing data quality control method and high-throughput RNA sequencing data quality control apparatus
Park et al. GAB2 amplification in squamous cell lung cancer of non-smokers
CN110970091B (en) Label quality control method and device
CN101663408B (en) Melting curve analyzing method and melting curve analyzing device
CN113488105B (en) Microsatellite locus based on amplicon next-generation sequencing MSI detection, screening method and application thereof
CN111292807B (en) Method for analyzing double cells in single-cell transcriptome data
CN106591443B (en) A kind of method for analyzing and processing crops fluorescent capillary electrophoresis tube collection of illustrative plates
CN117743783A (en) Method and system for correcting electrophoretic mobility shift
WO2019014455A4 (en) Pattern recognition system
CN106480215B (en) A kind of method and its device for analyzing and processing crops SSR marker collection of illustrative plates
EP1798651B1 (en) Gene information display method and apparatus
CN115948521A (en) Method for detecting aneuploid missing chromosome information
CN107746884B (en) AFLP primer combination product, kit and method for identifying individual and variety of beef cattle
US20240221954A1 (en) Disease prediction methods and devices, electronic devices, and computer readable storage media
CN111724860B (en) Method and device for identifying chromatin open area based on sequencing data
JP7099759B1 (en) Mechanical detection of candidate break points for variants in the number of copies on the genome sequence
CN108060239B (en) Primer pair combination product, kit and method for distinguishing yak from non-yak

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
TR01 Transfer of patent right

Effective date of registration: 20201028

Address after: Room 107, building 5, Zizhuyuan Road, Haidian District, Beijing 100048

Patentee after: Beijing Academy of Agricultural Sciences Seed Industry Technology Co.,Ltd.

Address before: 100083 Beijing City, Haidian District Zhongguancun Road No. 18 smartfortune International Building C-1511

Patentee before: Beijing Huasheng Hengye Technology Co.,Ltd.

TR01 Transfer of patent right