WO2013028699A2 - Cell line discernment using short tandem repeat - Google Patents
Cell line discernment using short tandem repeat Download PDFInfo
- Publication number
- WO2013028699A2 WO2013028699A2 PCT/US2012/051746 US2012051746W WO2013028699A2 WO 2013028699 A2 WO2013028699 A2 WO 2013028699A2 US 2012051746 W US2012051746 W US 2012051746W WO 2013028699 A2 WO2013028699 A2 WO 2013028699A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- cell line
- probable
- target
- str
- line
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/5005—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells
- G01N33/5008—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells for testing or evaluating the effect of chemical or biological compounds, e.g. drugs, cosmetics
- G01N33/5044—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells for testing or evaluating the effect of chemical or biological compounds, e.g. drugs, cosmetics involving specific cell types
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
- G16B20/20—Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6813—Hybridisation assays
- C12Q1/6827—Hybridisation assays for detection of mutation or polymorphism
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6809—Methods for determination or identification of nucleic acids involving differential detection
Definitions
- the purpose of the work is to allow scientists to automatically detect misidentified and unstable lines, oversee the long-term stability of a line for observing general loss of default characteristics in the line, and detect potential of cross-contamination.
- STR Short Tandem Repeat
- Typical STR profile represents the frequency of a class of tetra and/or penta-nucleotides that are short in sequence selected from different regions. These regions are selected in a rigorous and distributed fashion across the genome. These short tetramers are responsible and/or associated with the identification of certain traits that do not necessarily lie on the coding protein regions. Therefore, even with two lines that are found genetically identical at a time, they both may (or may not) confer different protein expressions at that same time. These repeats are called alleles, and are collected across different and variant regions on the genome called loci or reference loci (also known as reference markers). Many of these tetramer repeats can be associated with one or more disease.
- traits can be polygenic (i.e. multiple genes are responsible for the expression of a single physical trait). Identifying a phenotypic class inheritance from a known genotype profile can be a complex routine. As an example, Waardenburg Syndrome is inherited as an autosomal dominant and can cause an individual to have different color eyes among other communication disorders mainly due mutation in PAX3 gene in chromosome 2q35-2q37 [Faivre et al. 1997]. Thus, genes that are responsible for eye color are polygenic and, therefore, different alleles from this trait on a specified locus should be highly unique to an individual [NIH Pub. No. 91-3260. March 1999]. Unlike monogenic traits, polygenic traits do not follow an implicit pattern. Thus, expressing alleles from compound tetra- nucleotides rather than simple tetra-nucleotide may by itself complicate any potential interpretation that may be drawn.
- the tetramer 'GATA' is repeated 16 times on chromosome seven on the reference locus D7S280. Therefore, if different alleles of locus (D7S280) can have any pair combination of alleles from 5 to 16, i.e. 12 homozygous pairs and 66 heterozygous possible pair combinations, a total of 78 possibilities of pairs of alleles (i.e. 78 genotypes, or 78 possible local STR profiles at locus D7S280) can explicitly exist for D7S280. Complete genotyping information and frequency of repeats for the 16 loci are described in Butler, 2005. These alleles are then used to genotype a line. Each locus varies in length from 100 to 350 bp.
- locus like D21 S1 1, have more explicit repeat pattern (i.e. alleles) that appear as a mixture of tetra nucleotides and tri-nucleotides that can vary in size by single base pair. Thus, their repeats are expressed in the allele as x.2 or x.3.
- An example of an allele would be 23.3, showing the tetra pattern repeated in whole 23 times and the tri nucleotides is 3 nucleotides off the motif.
- the TH01 locus shows a typical case of repeat with the designation of the incomplete tetra-nucleotides and it often appears as 9.3 of commonly referenced AATG motif and one repeat with ATG [Puers et al 1993].
- FGA locus is another example of micro- repeats that often repeats as 22.2, that is, 22 repeats with one incomplete tetra-nucleotides formed with 2 base [barber et al 1996].
- TPOX is a simple tetra-nucleotide repeat found in the intron 10 of the human thyroid peroxidase gene near the very end of the short arm of chromosome 2.
- the Motif governing TPOX is AATG and it is considered the least polymorphic locus.
- Short Tandem Repeat (STR) technology has enabled fast DNA informative data for further analysis. These data however, need to be carefully interpreted and analyzed when verifying identities and stability of the lines. Complications arise in cell lines such as instability (microsatellite instability), amplification of a region, deletion, duplication, cell line mixing and non-human cell line contamination. While STR technology generates genomic profiles that offers a coded pattern, the interpretation of the coded pattern often result in more than one possible inference, thus, leaving scientists with great amount of time to consume retrieving and analyzing these data.
- cell line plays as the only available modeling tool to monitor the efficacy of anti-cancer agents.
- model may behave slightly different than it would behave inside its natural host "the human body" for diseased and/or free of disease host as opposed to a line that is cloned and preserved in a media.
- loss of default characteristics of a model line may be carried faster, as opposed to the loss of default characteristics of a similar characterized cell line living in its natural host.
- Cultured cell lines are one of the most used prototypes in pathway discovery, and drug pre-screening research. It is estimated that over fifty percent of the cell lines used do not correctly confer the characteristics assumed for these lines that were published and/or used in research [Rees at al. 1981, Drexler et al. 1999, Drexler et al. 2003]. The direct and/or indirect clinical impact of such consequences is misleading and can be devastating.
- Similarity techniques also known as matching techniques, are measures of identifications that are heavily used in artificial intelligence such as classification networks, and fuzzy logic [Weaver et al. 1949, Miller 1955, Watson et al. 1980, Gibbs et al. 2002].
- the goal is to acquire a similarity score between a target feature vector and a source or reference feature vector.
- most of the matching techniques adopted on genomic data is based on the concept of the probability of shared information to the overall information.
- similarity coefficients that are frequently used with genomic data these are based on either the use of Euclidean, Jaccard, Cosine, or Pearson Correlation among others. Correlation Techniques often assume linear relations among variables.
- ATC American Type Culture Collection Standards Development
- ASN-0002 and other companies such as Genetica, Human Cell line Authentication Laboratory offer services to authenticate Cell Lines using the STR DNA Technology, which is the method mostly used to certify the identity of cell lines. However, this authentication does not include evaluating the stability and cross-contamination of cell lines.
- the process of obtaining a DNA profile using highly polymorphic STR repeats includes the isolation of DNA from a cell line specimen, the amplification of targeted areas of that DNA, and analysis of the amplified product using a genetic analyzer.
- Embodiments of the invention include one or more of the following: 1 -Integrate the short tandem repeat profile of one or more target lines into a DBMS architecture designed to host the high throughput STR profile data generated. 2- Match each target STR profile to other Reference STR profiles, and disseminate if the line is an alias to one or more pre-recognized lines from the public and/or private repositories that are housed in a DBMS.
- An Enterprise In-House Laboratory Information Management System is designed and implemented to allow scientists to periodically observe the long-term stability of the line and to detect miss-identified lines prior to using the line to selectively experiment any targeted molecule interactions.
- LIMS Enterprise In-House Laboratory Information Management System
- An automated system to detect the miss-identified lines and provide decision support for the stability of the line is a critical step to bring any further in-vitro investigation to a solid ground.
- the system described here is deemed successful to deliver solutions to the cell culture dilemma.
- Embodiments of the system detects if a target line is mis identified, detects if the target line is stable or unstable, and checks if there is probable cross-contamination from an extracted pattern. This pattern represents the proportional changes in the complex tetra- nucleotides alleles. If the target line is deemed unique and stable a reference profile of the line is then created thus allowing dynamic platform for future matching.
- the system offers a dynamic, robust, and efficient detection of misidentified, and/or unstable cell lines, and saves a great deal of scientists' time and efforts as it ensures the delivery of eradicated and disseminated information to the scientific community.
- Embodiments of the invention relate to a system, method and a computer readable medium embodied by a computer readable medium, which detects misidentified, unstable and cross-contaminated cell lines using STR profiles.
- An embodiment of the invention is a method comprising a) obtaining a sample comprising a target cell line; b) testing the sample to determine an STR profile of the target cell line; c) receiving information corresponding to a STR profile of a reference cell line; d) comparing the alleles in the STR profile of the target to the alleles in the reference STR profile and calculating a probable identity match; e) calculating probable instability of the target cell line; calculating probable cross-contamination of the target cell line; f) and generating an output comprising the probable identify match of the target cell line with the cross reference cell line, g) the probable instability of the target line, and the probable cross- contamination of the target line.
- the method may additionally comprise identifying the target cell line as the reference cell line or a different cell line from the probable identity match.
- the method may also additionally comprise receiving information corresponding to the instability of the reference cell line and the potential of cross-contamination of the reference cell line.
- calculating the probable instability of the target cell line comprises equation 3.
- calculating the probable cross-contamination of the target line comprises equations 3 and 4.
- comparing the alleles in the STR profile of the target to the alleles in the reference STR profile and calculating the probable identify match comprises equation 1.
- the information corresponding to the STR profile of the target cell may comprise a different number of STR loci than the corresponding number of STR loci that are available in the STR profile of the reference cell line.
- the information corresponding to a STR profile of a reference cell line is received from an electronic database management system of STR profiles or other data set obtained from public STR profiles. Further, the output may be sent to an electronic database management system for storage. In an embodiment of the invention, the method is performed for each reference cell line in the database.
- Another embodiment of the invention is a system for cell line discernment for use by a software application, the system comprising a processor in communication with a memory, where: the memory stores processor-executable program code; and the processor is configured to be operative in conjunction with the processor-executable program code to: a) receive information corresponding to an STR profile of a target cell line; b) receive information corresponding to an STR profile of a reference cell line; c) compare the alleles in the STR profile of the target to the alleles in the reference STR profile and to calculate a probable identity match; d) calculate probable instability of the target cell line; e) calculate probable cross-contamination of the target cell line; and f) generate an output comprising the probable identify match of the target cell line with the cross reference cell line, the probable instability of the target line, and the probable cross-contamination of the target line.
- the system may also include a server comprising an electronic database system that is interfaced with the processor.
- the embodiment may additionally comprise identifying the target cell line as the reference cell line or a different cell line from the probable identity match. Further, additional information may be received corresponding to the instability of the reference cell line and the potential cross-contamination of the reference cell line.
- the calculating instability comprises calculating the probable instability of the target cell line with equation 3.
- calculating the probable cross-contamination of the target line comprises equations 3 and 4.
- the comparing the alleles comprise calculating the probable identity match with equation 1.
- comparing the information corresponding to the STR profile of the target cell may comprise a different number of STR loci than the corresponding number of STR loci that are available in the STR profile of the reference cell line.
- the additional information is received from an electronic database management system of STR profiles or other data set obtained from public STR profiles.
- the output generated may sent to an electronic database management system.
- the embodiment may additionally comprise repeating the executable embodiments of the invention for each reference cell line in the electronic database management system.
- Another embodiment of the invention is a computer readable medium comprising computer-usable program code executable to perform operations comprising: a) receiving information corresponding to a STR profile of a target cell line; b) receiving information corresponding to a STR profile of a reference cell line; c) comparing the alleles in the STR profile of the target to the alleles in the reference STR profile and to calculate a probable identity match; d) calculating probable instability of the target cell line; e) calculating probable cross-contamination of the target cell line; and f) sending the probable identify match of the target cell line with the cross reference cell line, the probable instability of the target line, and the probable cross-contamination of the target line.
- the computer readable medium may additionally comprise identifying the target cell line as the reference cell line or a different cell line from the probable identity match.
- the computer readable medium may also additionally comprise receiving information corresponding to the instability of the reference cell line and the potential cross-contamination of the reference cell line.
- calculating the probable instability of the target cell line comprises equation 3.
- calculating the probable cross-contamination of the target line comprises equations 3 and 4.
- comparing the alleles comprises equation 1.
- the information corresponding to the STR profile of the target cell may comprise a different number of STR loci than the corresponding number of STR loci that are available in the STR profile of the reference cell line.
- the computer readable medium may additional comprise interfacing with an electronic database system.
- the receiving information corresponding to a STR profile of a reference cell line receives information from an electronic database management system of STR profiles or other data set obtained from public STR profiles.
- the data is sent to an electronic database management system.
- the computer readable medium may further comprise executing the method for each reference cell line in the database.
- Fig. 1 shows the association of global hit score, instability score, and potential of cross-contamination score including proper dependencies affecting each score and associated rules.
- Fig. 2 illustrates a 16 locus analysis using hit score, instability score and potential of cross-contamination score showing the alignment STR channels for target line A and an instant reference line B.
- Fig. 3 illustrates the reference and target locus and the hit score, instability score and potential of cross-contamination score, and associated embodiments of mathematical equations.
- Fig. 4 is a schematic flow which illustrates the steps of an embodiment of the algorithm.
- Fig. 5 illustrates STR object models and their relations as used in an embodiment of the system.
- Fig. 6 illustrates a summary of the top-down work flow and example input files.
- Fig. 7 is an embodiment of the system architecture illustrating different architecture layers and elements.
- Fig. 8 is schematic of the layers within an example system architecture showing a) a cross section of the system architecture and b) a layout of an embodiment of the network.
- Fig. 9 is a screenshot of a) a log in screen and b) a STR enterprise web interface.
- Fig. 10 is a screenshot of a STR LIMS submission panel.
- Fig. 1 l is a screenshot of a STR LIMS self-hit results panel.
- Fig. 12 is a screenshot from the STR LIMS of a target/reference preview panel.
- Fig. 13 is a screenshot of a sample STR LIMS instant STR profile search.
- Fig. 14 is a screen shot of the sample queue enterprise LIMS showing sample submission.
- Fig. 15 is the linguistic cluster results of the target line C8161-C9.
- Fig. 16 is a table of selected target lines that are highly matched to other references in the DBMS.
- Fig. 17 is a table of selected target lines highly matched to other references in the DBMS and were marked as unstable.
- Fig. 18 is a table of selected target line (HCT1 16-B) marked as unstable and probable mixed with reference lines SKOV3 and ARTFLOX.
- Fig. 19 is a table of the sensitivity (se) and specificity (sp) analysis for the instability and probable mix detection evaluation that is prepared using pre-selected 77 cultures.
- Fig. 20 is a linguistic cluster of a selected target line AU565.
- Fig. 21 is a table of selected target lines detected as alias to other references and marked as stable.
- Fig. 22 is a table of four selected target lines that were marked as unstable.
- Fig. 23 is a table of selected target lines matched with other reference lines and a) marked as alias and b) partial match, unstable and probable mix.
- Fig. 24 is graph analysis of the HCT116-A and HCT116-B family of cultured cells.
- Fig. 25 is graph analysis for line AU565 hits with other cultured lines.
- Fig. 26 is graph analysis of the CAKI1 family with other cultured lines.
- Fig. 27 is graph analysis of the IGROV1 family with other cultured lines.
- Fig. 28 is a table of the collection of thirty seven cultures that were mixed and not mixed and used for probable mix testing.
- Fig. 29 is a table of a linguistic cluster of a selected target lines CAKI1, HCT116, and TK-10 with selected reference lines.
- genotype STR profile collected is very useful to inspect identity, monitor long term stability of a line, and can be used to infer mutation or mixed lines being observed either across different cultures and/or across large family of passages in an assorted culture.
- the following section is used to outline an embodiment of the method and system flow.
- a 16 -core CODIS system may be used.
- Table 1 illustrates a 16 CODIS generated (i.e. STR profile).
- Other systems may be used such as FFFL, GammaSTR, and AmpF/STR systems.
- the two STR profiles to be matched are produced by the same platform to exclude any potential shifts per locus and that may exist with some loci due to different platforms such as the difference between Power Plex 1.1, and Power Plex 16 for some of the loci, and that all other precautions in handling low DNA are preserved and taken according to all guidelines.
- different systems may be matched.
- Reference lines may come from any STR databases such as ATCC, DSMZ, RIKEN, JCRB and NIST, in addition to institutional databases.
- Target line STR profiles are examined against reference STR line profiles stored as references in the DBMS creating a node in the target linguistic cluster being inspected.
- the first 16 scores represents the local hit scores between the target line (A) and the reference line (B).
- the second 32 scores represent the target line, where 16 scores indicates stability of 16 local locus followed by 16 score inspecting cross-contamination for the target line (A).
- the first 16 local hit scores are then weighted to calculate the global weighted hit scores for the target (A).
- the mathematical notations are described in the following section.
- Matching Algorithm ⁇ , ⁇ , and v serve as three indicators for the global hit score, stability score, and potential of cross-contamination score respectively for a line.
- Fig. 1 shows the association of the three measures used to describe a line during an instant hit between a target line and a reference line. The figure simplifies the inclusions of each score with regard to the other two. (I) Misidentification
- a single amplification routine generates coded repeats pattern, which are used as inputs to the Algorithm.
- the Algorithm uses these sixteen channels, aiming at three objectives namely; 1- Calculate global weighted score between a target line A and a reference line B mathematically, noted as ⁇ .
- 2- State if a line is stable or unstable based on a weighted stability score noted as ⁇ .
- 3- Inspect if a line expresses potential of cross-contamination based on weighted score noted as v, as illustrated in Fig. 2.
- Equation 1 and 2 are given below.
- Every allele in the target domain is inspected against every allele in the corresponding reference allelic domain to calculate ⁇ ; where, a j represents all target instant allele being visited and, thus, domain of all instances of circumstances are accounted
- Equation 1 depicts both absolute and conditional partial match while providing local score between finite sets A (target) and B (reference).
- the a j is the new pattern between instant pair a, and b allele as it defines the concise match without creating anarchy. The later is important to account for small aberrations between the default pair of alleles being matched and ensures they are accounted for and not ignored.
- target allele 9 matches reference allele 9.3 with 0.7 if reference lack instance of 9.
- the notations used denote the cardinality of A and B sets.
- I a and 3 ⁇ 4
- the values I a represent the inferred influence of the relevant profiles known as target and reference, respectively (i.e. number of available informative locus (i.e. fully matched, and partially matched).
- the value I a represents the inferred influence rather than the abstract influence across all loci with available information.
- ⁇ defines the ratio between the inferred to the abstract influences.
- STR allelic profiles may exhibit compound repeats, every individuals is either homozygous or heterozygous i.e. showing inherited alleles from one parent or both parents. Despite the fact that these locus are chosen to represent alleles with high variations in a population, it does not change the reality that one individual cannot possibly represent all incidence of populated events. Therefore, a template filter is designed to extract how variant subsequent repeats are presented per locus STR. The filter is sensitive to the smallest variations and preserves heterozygosis events as the default. Equation 3 gives power not only to the variation in the allelic repeats but also to their spatial order from the first 2 alleles.
- the notation ⁇ depict the per-locus differentiated repeats as a measure of variation to highlight any instability per locus notated as ⁇ ; a mathematical description of ⁇ ; is in equation 3.
- 3 ⁇ 4 used to evaluate the pattern extracted to detect instability as notated in equation (3). Stability is then preserved if for all loci thus preserving heterozygosis events, and maintains the order of the variant alleles where, the index j is the repeats order as its variation is being evaluated.
- STR allelic profiles may exhibit compound repeats, every individual is either homozygous or heterozygous i.e. showing inherited alleles from one parent or both parents. Despite the fact that these locus are chosen to represent alleles with high variations in a population, one individual cannot possibly represent all incidence of populated events. Therefore, in one embodiment, a filter template is used to read how variant subsequent alleles are presented in every local STR given power not only to the variation in the allelic repeats but also to their spatial order from the first 2 alleles. The variation to highlight any instability per locus is notated as Equation 3 is given below.
- a cross-contaminated line can be caused from a variety of reasons.
- One example is if regulations in handling low DNA are not met. This can cause the scientist DNA to show in the STR profile leaving a cross-contamination call to that STR. Further, intra-species cross- contamination is expected particularly in handling cell culture.
- tetra-nucleotide alleles are associated with certain incomplete polymorphisms that appear as x.2, x.3 in the pattern. Some STR alleles are compound tetra-nucleotides. If a line expressed frequent fragments of incomplete polymorphisms that are forming series of compound tetra-nucleotides per locus that raises the possibility that it captures sign of cross- contaminations thus loosing pattern of solid and discrete repeats.
- equation 3 and equation 4 are examined for any indication of probable mix allowing sensing all variations that may be presented with simple and/or compound repeats in a locus.
- probability of having potential mixed lines will increase if a target will show sign of instability from equation 3 where alone with or without
- cross-contamination score can be derived from a stability score i.e. (i.e. v entails ⁇ ).
- a filter to collect the ratio between variations in the irrational domain of an allele to the variation in its rational domain was designed to indicate cross-contamination score as follows in Equation 4:
- Equation 4 extract patterns are associated with the complex alleles rather than the simple alleles in the set of local alleles. Furthermore this pattern examines the proportional change in consecutive complex alleles.
- allelic set as ⁇ 1 1, 12, 13, 19, 9.3, 23.2 ⁇ is equal to (0.0071942446) 2 where, The lower found the higher
- ⁇ depicts the consecutive variations for every 2 consecutive compound alleles in the set that is denoting either simple or complex repeats.
- ⁇ assumes heterozygosis as a default.
- the symbol y j denotes the fraction in the complex allele if there exists complex alleles while x j represents the complex allele.
- Fig. 4 An embodiment of the algorithm or procedure is illustrated in Fig. 4 and given below.
- STR object model has been designed to reproduce all elements needed to describe an instant profile with 16 pre-known loci in a map where each node in the map has a set of alleles (see Fig. 5).
- Line 13-14 lists the instant calculation to scale the severity of cross-contamination if conditions stated are met for the calculation.
- line 16 if there exists more than a pair alleles in the instant locus being examined then line 17 evaluates the differential change that is also being accumulated. However, another parameter ⁇ ° is incremented as ⁇ ; > ⁇ °. Thus, the target is marked as unstable if 1 or more locus is detected meeting ⁇ ° value.
- Line 18-23 is calculated after all alleles in the Loci domain are examined. Specifically, line 22-23 calculates Local hit and updates the global hit ⁇ . The above steps are repeated until all 16 loci are exhausted. Line 25-27 equations establish the weight. Line 28 computes weighted global hit.
- a distance value between target and instant reference is then calculated in line 29.
- line 30-31 pushes the final linguistic cluster between the target A and the instant reference STR object B k while updating the linguistic cluster map for this Target A as the routine is being exhausted over all references collected that ends with a return of the linguistic cluster for a target line A as LC A .
- the system allows public STR profiles generated by other facilities to be matched, which may enable smart logic while mining lines from different laboratories.
- a panel was designed to allow scientists to retrieved matched lines that were subjected to varieties of modification and/or panel of lines collected from across culture i.e. different passages.
- a module is "[a] self- contained hardware or software component that interacts with a larger system.” Alan Freedman, "The Computer Glossary" 268 (8th ed. 1998).
- a module comprises a machine or machines executable instructions.
- a module may be implemented as a hardware circuit comprising custom VLSI circuits or gate arrays, off-the-shelf semiconductors such as logic chips, transistors, or other discrete components.
- a module may also be implemented in programmable hardware devices such as field programmable gate arrays, programmable array logic, programmable logic devices or the like.
- Modules may also include software-defined units or instructions, that when executed by a processing machine or device, transform data stored on a data storage device from a first state to a second state.
- An identified module of executable code may, for instance, comprise one or more physical or logical blocks of computer instructions, which may be organized as an object, procedure, or function. Nevertheless, the executables of an identified module need not be physically located together, but may comprise disparate instructions stored in different locations which, when joined logically together, comprise the module, and when executed by the processor, achieve the stated data transformation. Indeed, a module of executable code may be a single instruction, or many instructions, and may even be distributed over several different code segments, among different programs, and across several memory devices. Similarly, operational data may be identified and illustrated herein within modules, and may be embodied in any suitable form and organized within any suitable type of data structure. The operational data may be collected as a single data set, or may be distributed over different locations including over different storage devices.
- the system for cell line discernment using STR profiles may include a server, a data storage device, a network, and a user interface device.
- the system may include a storage controller, or storage server configured to manage data communications between the data storage device, and the server or other components in communication with the network.
- the storage controller may be coupled to the network.
- the user interface device is referred to broadly and is intended to encompass a suitable processor-based device such as a desktop computer, a laptop computer, a Personal Digital Assistant (PDA), a mobile communication device or organizer device having access to a network.
- the user interface device may access the Internet to access a web application or web services hosted by the server and provide a user interface for enabling a user to enter or receive information. For example, the user may enter text files such as the one shown in Fig. 6 into the system.
- the network may facilitate communications of data between the server and the user interface device, the server and a database, and/or a user and a database.
- the network may include any type of communications network including, but not limited to, a direct PC to PC connection, a local area network (LAN), a wide area network (WAN), a modem to modem connection, the Internet, a combination of the above, or any other communications network now known or later developed within the networking arts which permits two or more computers to communicate, one with another.
- STR profiles against reference cell lines STR profiles to determine the identity, stability, and probable of cross-contamination of the target cell lines.
- the server may access data stored in a data storage device via a Storage Area Network (SAN) connection, a LAN, a data bus, or the like.
- SAN Storage Area Network
- a desktop computer is configured to analyze target cell lines STR profiles against reference cell lines STR profiles to determine the identity, stability and probable of cross-contamination of the target cell lines.
- the data storage device may include a hard disk, including hard disks arranged in an Redundant Array of Independent Disks (RAID) array, a tape storage drive comprising a magnetic tape data storage device, an optical storage device, or the like.
- the data storage device may store reference STR profiles and/or target cell line STR profiles.
- the data may be arranged in an electronic database and accessible through Structured Query Language (SQL) queries, or other database query languages or operations.
- SQL Structured Query Language
- a data storage device may be located on a desktop type computer, on a server, and/or on a dedicated database device.
- the server may host a software application configured for cell line discernment using STR profiles.
- a desktop type computer may host a software application configured for cell line discernment using STR profiles.
- the software application may further include modules for interfacing with the data storage devices, interfacing with a network, interfacing with a user, and the like.
- the server or desktop type computer may host an engine, application plug-in, or application -programming interface (API).
- the server may host a web service or web accessible software application.
- the software application may be embodied on a computer readable medium.
- STR profiles of target cell lines are input into the system through an input file and/or directly from a STR profile-generating machine.
- the STR target cell line profiles are compared to reference cell line STR profiles that are located in one or more database.
- the process of inputting files and comparing the STR profiles may be automated.
- the automation may occur at the system that hosts the database, in the users system, or in a system that bridges the user and the database.
- the input files may include one or more STR profile.
- FIG. 8a illustrates the layers of the architecture.
- the exterior of the circle represents user interaction with the system such as users from lab groups or institutions who interacts with the system through a web browser and/or client APIs, for example.
- the interactions with the system may be through an enterprise Laboratory Information Management System (LIMS) and/or web portal.
- LIMS enterprise Laboratory Information Management System
- Scheduled processes layer which interface with one or more database handles intermediate transactions that do not need an interface with users but rather needs layer of automated scheduled processes.
- Web services layer handles all users and scheduled processes requests in connections with DBMS layers or AI layers.
- Fig. 8b illustrates an example of a layout of a cancer network. Cell profiles may come from any type of cell line in the cancer network.
- Fig. 9 is an example of a system graphical user interface.
- the user can perform various STR related activities such as submitting an STR profile, matching an STR profile or searching for an STR profile.
- the user submits a STR profile and identifies where the profile comes from, the species the cell line belongs to, the technology used to generate the STR profile, the lab ID, the set ID and a description.
- the user can preview self -hit results and preview the STR client sample details.
- Multiple target reference cell STR profiles can be submitted in batch. This is shown in more detail in Fig. 10-12.
- Fig. 10 illustrates submitting STR information for a target cell, as described above.
- Fig. 10 illustrates submitting STR information for a target cell, as described above.
- FIG. 11 shows a view of the preview of STR client self hit results with target cell lines matched against multiple reference cell lines.
- the self hit results show the identity of the target and reference cell lines, the min hit %, the average distance, the hit score, and the probable stability.
- Fig. 12 shows a screen shot of the target/reference preview window and/or panel. This window illustrates the STR profiles of the input target cells STR profiles, in addition to the hit results.
- Fig. 13 is a screenshot showing an instant search of STR samples. The example search is based on any STR with completed analysis.
- Fig. 14 is an illustration showing a multiple sample submission. The sample names, cell lines, tissue type, state, species and STR profiles are shown.
- Fig. 6, as discussed above, is an example STR profile input file generated from GeneMapper and another format of input file that is used or been transformed to, and the system top-down flow to output the match, stability and potential cross-contamination.
- EXAMPLE 1 is an example STR profile input file
- NCI National Cancer Institute
- Table 1 represents the 16 CODIS generated (i.e. STR profile) for two lines.
- DNA was extracted from cell lines using QiaAMP maxi preps (Qiagen) and DNA was quantified by Nanodrop spectrometry.
- Cell line profiles obtained by STR DNA fingerprinting using the AmpFfSTR identifier kit according to manufacturer instructions (Applied Biosystems cat 4322288).
- ATCC American Type Culture Collection Standards Development
- STR DNA Technology which is the method mostly used to certify the identity of cell lines.
- the STR profiles were compared to known ATCC fingerprints, to the Cell Line Integrated Molecular Authentication (CLIMA) database version 0.1.200808, (Nucleic Acids Research 37:D925-D932 PMCID: PMC2686526), and to the complete 16 loci and to the MD Anderson fingerprint database where the STR profiles either matched known DNA fingerprints or were unique.
- the lines HCT-1 16 (NCI-60), and IGROV-1, X, Y were obtained as part of the NCI-60 cell line collection from the National Cancer Institute.
- OVCAR429, OVCAR432, OVACR433, and HCT-1 16 cell lines were provided by several labs. NCI-60 cell lines were grown in RPMI [Lorenzi et al. 2009, Stults et al. 2011], 10% FBS in a Co2 incubator. Same line is introduced in Savas et al. 2011 as it was obtained from the ATCC Organization Workgroup.
- the true positive fraction of the system was tested based on expert biologists and their definition of instability. Further, set of prepared mixtures were designed by expert biologist for the purpose of evaluating the instability and probable mix decisions placed, by the detection schema. A total of 77 cultures were used in this study. Amongst those a panel was prepared with 49 mixture lines which, were used to evaluate the se and sp of cross contamination detection. The culture profiles were presented to two expert biologists for the determination of instability those were compared to the system detection result. Overall, total of 77 lines were used for the se and sp analysis. Fig. 19 shows the se and sp analysis for the instability detection and probable cross contamination.
- the detection schema offers robust solution to the misidentification, instability, and probable cross contamination.
- TPF true positive fraction
- FPF false positive fraction
- ppv positive predictive value
- npv negative predictive value
- the instability detection is influenced by subjective assessment.
- the probable mix detection performance was analyzed.
- the TPF was 0.957 at FPF 3/49.
- the positive predictive value (ppv) for the probable mix detection is 0.937.
- Fig. 15 shows an example of matching results showing selected Target Line (C8161- C9) as it hits several reference lines.
- Fig. 16 shows matched cell lines where the target was misidentified.
- Fig. 17 shows examples of target lines that were highly matched an instant reference and both found unstable in a, and b examples.
- Fig. 18 shows a case of cross contamination.
- the HCT1 16-A and HCT116-B lines were also identified as unstable lines. However, HCT1 16-B was marked as probable cross-contamination in the system while proposing 2 other lines (SKOV3 and ARTFLOX) as potential of being cross contaminated with HCT1 16-B.
- SKOV3 and ARTFLOX 2 other lines
- Unlike HCT116-A the system detected the culture as unstable.
- Fig. 29 shows an example of matching results showing selected target line (CAKI1) as it hits several reference lines.
- a matching technique is presented based on pattern extracted from the STR profiles.
- the system offers a dynamic, and efficient detection of misidentified, and unstable cell lines.
- the system saves a great deal of scientists' time and efforts thus ensuring the welfare of the research as it ensures the delivery of disseminated information to the scientific community.
- the robust performance of the detection techniques suggests that the system provides an efficient and more advanced starting point for the automated cell cultures validation and authentication efforts that may or may not currently, exist for the detection of misidentification.
- cell line plays as the only available modeling tool in monitoring the efficacy of anti-cancer agents.
- good practice is to observe the long-term stability of a line.
- such model may behave slightly different than it would behave inside its natural host "the human body" for diseased and/or free of disease host as opposed to a line that are cloned and preserved in a media.
- the lost of default characteristics of a model line may be carried faster opposed to the lost of default characteristics of a similar characterized cell line living in its natural host.
- the only means to enable further studies to test such hypothesis is to have such a cost effective system in place that deliver fast, effective, and reliable mechanisms to decimate barriers that may otherwise obscure and/or overburden the translation of such large high throughput genomic data.
- Fig. 20 shows the linguistic cluster of the target cell line AU565 matching results.
- the matching results show the reference and target cell lines identification, as well as the min hit, min distance, weighted hit and the stability of the reference cell line.
- Fig. 21 are selected target lines detected as alias to reference lines and marked as stable. The matches are shown to be exact, exact, and exact with difference for the three cultures displayed.
- Fig. 22 shows the results of four target lines that were marked as unstable and
- Fig. 23 shows target lines that were matched with reference lines and marked as alias (a) and (b) partial match, unstable, and probably mix.
- Figs. 24-27 are graph analysis of families of cultured lines namely: HCT116, AU565, CAKI1, and IGROV.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Biotechnology (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Immunology (AREA)
- Genetics & Genomics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Evolutionary Biology (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Organic Chemistry (AREA)
- Urology & Nephrology (AREA)
- Microbiology (AREA)
- Biochemistry (AREA)
- Zoology (AREA)
- Cell Biology (AREA)
- Wood Science & Technology (AREA)
- Hematology (AREA)
- Tropical Medicine & Parasitology (AREA)
- Toxicology (AREA)
- General Engineering & Computer Science (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Design And Manufacture Of Integrated Circuits (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Image Processing (AREA)
Abstract
A method, system and computer readable media for evaluating cell line identification, instability, and cross-contamination using short tandem repeats is disclosed. The method system and computer readable media allow scientists to automatically detect misidentified and unstable lines, oversee the long-term instability of a line, and detect potential of cross-contamination or general loss of default characteristics in the line.
Description
DESCRIPTION
CELL LINE DISCERNMENT USING SHORT TANDEM REPEAT
The application claims priority to U. S. Provisional Patent Application No. 61/525,793 filed August 21, 2011, which is incorporated herein by reference in its entirety.
This invention was made with government support under Grant Number CA 16672 awarded by the National Cancer Institute at the National Institutes of Health. The United States Government has certain rights in the invention. TECHNICAL FIELD
The purpose of the work is to allow scientists to automatically detect misidentified and unstable lines, oversee the long-term stability of a line for observing general loss of default characteristics in the line, and detect potential of cross-contamination.
BACKGROUND OF THE INVENTION Cell lines have been increasingly used in translational medicine for prototyping the tumor tissues, their host, and in the discovery of effective therapeutic agents. A typical prototype should be stable and consistently delivering reliable and inclusive information. Reliability of further inferences radically depends on accurate identification of the cell line being experimentally observed. Human cell lines that are being actively used over the years in major and critical cancer research have been compromised. Cell line validation is a critical issue for both scientific publications and grant applications, as articles are being published with cell lines misidentified [Schweppe et al., 2009; Drexler et al, 2003; MacLeod et al, 1999; Capes-Davis et al, 2010].
The Short Tandem Repeat (STR) technique, also known as flanking microsatellite repeat, was being used in forensic sciences even before cell models existed and it became popular around 1990 [Clayton et al. 1995, Hagelberg et al. 1991]. Typical STR profile represents the frequency of a class of tetra and/or penta-nucleotides that are short in sequence selected from different regions. These regions are selected in a rigorous and distributed fashion across the genome. These short tetramers are responsible and/or associated with the identification of certain traits that do not necessarily lie on the coding protein regions.
Therefore, even with two lines that are found genetically identical at a time, they both may (or may not) confer different protein expressions at that same time. These repeats are called alleles, and are collected across different and variant regions on the genome called loci or reference loci (also known as reference markers). Many of these tetramer repeats can be associated with one or more disease.
Given that traits can be polygenic (i.e. multiple genes are responsible for the expression of a single physical trait). Identifying a phenotypic class inheritance from a known genotype profile can be a complex routine. As an example, Waardenburg Syndrome is inherited as an autosomal dominant and can cause an individual to have different color eyes among other communication disorders mainly due mutation in PAX3 gene in chromosome 2q35-2q37 [Faivre et al. 1997]. Thus, genes that are responsible for eye color are polygenic and, therefore, different alleles from this trait on a specified locus should be highly unique to an individual [NIH Pub. No. 91-3260. March 1999]. Unlike monogenic traits, polygenic traits do not follow an implicit pattern. Thus, expressing alleles from compound tetra- nucleotides rather than simple tetra-nucleotide may by itself complicate any potential interpretation that may be drawn.
There are several advantages of using short sequences rather than long sequences. One advantage is that it is easier to amplify short sequences with a low amount of DNA. Second, even with DNA degradation it is easier to spot PCR stutter caused by DNA degradation with short sequence than with a longer sequence. Third, short tetra or penta- nucleotide are often implicit (i.e. easier to measure), and are more likely to be free from any other artifacts that may often be present when PCR is done on longer sequences, particularly if mutations are present. STR became popular DNA markers because these short repeats are easily amplified by the polymerase chain reaction (PCR) without any differential amplification [Butler 2005]. Thus, these alleles can then be used to identify the genotype profile of instant cell lines conditioned by the fact that these 16 loci are distributed across the genome and uniquely represent different individuals across large and different ethnic populations.
As an example, the tetramer 'GATA' is repeated 16 times on chromosome seven on the reference locus D7S280. Therefore, if different alleles of locus (D7S280) can have any pair combination of alleles from 5 to 16, i.e. 12 homozygous pairs and 66 heterozygous possible pair combinations, a total of 78 possibilities of pairs of alleles (i.e. 78 genotypes, or
78 possible local STR profiles at locus D7S280) can explicitly exist for D7S280. Complete genotyping information and frequency of repeats for the 16 loci are described in Butler, 2005. These alleles are then used to genotype a line. Each locus varies in length from 100 to 350 bp.
Furthermore, some locus, like D21 S1 1, have more explicit repeat pattern (i.e. alleles) that appear as a mixture of tetra nucleotides and tri-nucleotides that can vary in size by single base pair. Thus, their repeats are expressed in the allele as x.2 or x.3. An example of an allele would be 23.3, showing the tetra pattern repeated in whole 23 times and the tri nucleotides is 3 nucleotides off the motif. The TH01 locus shows a typical case of repeat with the designation of the incomplete tetra-nucleotides and it often appears as 9.3 of commonly referenced AATG motif and one repeat with ATG [Puers et al 1993].
FGA locus is another example of micro- repeats that often repeats as 22.2, that is, 22 repeats with one incomplete tetra-nucleotides formed with 2 base [barber et al 1996]. TPOX is a simple tetra-nucleotide repeat found in the intron 10 of the human thyroid peroxidase gene near the very end of the short arm of chromosome 2. The Motif governing TPOX is AATG and it is considered the least polymorphic locus. The same PCR primer for TPOX contained in Power Plex 16 with no shift between Power Plex 1.1 or 2.1 of which careful validation needs to be done before comparing 2 STR profiles, as alleles may be different due to shifting existing between two different platforms and kits used [Butler 2005].
Short Tandem Repeat (STR) technology has enabled fast DNA informative data for further analysis. These data however, need to be carefully interpreted and analyzed when verifying identities and stability of the lines. Complications arise in cell lines such as instability (microsatellite instability), amplification of a region, deletion, duplication, cell line mixing and non-human cell line contamination. While STR technology generates genomic profiles that offers a coded pattern, the interpretation of the coded pattern often result in more than one possible inference, thus, leaving scientists with great amount of time to consume retrieving and analyzing these data.
Overall, it is a good practice to observe the long-term stability of a line. Thus far, cell line plays as the only available modeling tool to monitor the efficacy of anti-cancer agents. However, it is expected that such model may behave slightly different than it would behave inside its natural host "the human body" for diseased and/or free of disease host as opposed to a line that is cloned and preserved in a media. For example, it is expected that the loss of
default characteristics of a model line may be carried faster, as opposed to the loss of default characteristics of a similar characterized cell line living in its natural host. Overall, the only means to enable new innovative experimental studies to test such hypothesis is to have a cost effective system in place that delivers fast, effective, and reliable mechanisms to decimate barriers that may otherwise obscure and/or overburden the translation of such large high throughput genomic data such as the chaos of the sample and/or line miss-identified because of redundant and repetitive manual manipulation and the lack of standards being followed.
Cultured cell lines are one of the most used prototypes in pathway discovery, and drug pre-screening research. It is estimated that over fifty percent of the cell lines used do not correctly confer the characteristics assumed for these lines that were published and/or used in research [Rees at al. 1981, Drexler et al. 1999, Drexler et al. 2003]. The direct and/or indirect clinical impact of such consequences is misleading and can be devastating.
Using incorrect cell lines not only affects researchers who have had to retract papers but also has implications on data generated in the past [Lucey et al. 2006]. Incorrect cell line use has resulted in erroneous conclusions and has set back research in Mesenchymal Stem Cell (MSC) transplantation, thyroid cancer, and esophageal cancers [Boonstra et al. 2010, Vogel et al. 2010]. Incorrect cell lines have also had an effect on clinical practice; data from cell lines of the wrong tumor type have been used to justify clinical trials. Further, lists of misidentified cell cultures can be found on web sites from company that characterize and sell cultures such as German Collection of Microorganisms and Cell Cultures (DMSZ), Japanese Collection of Research Bio-resources (JCRB), and American Type Culture Collection (ATCC). Because of these high profile failures, many journals now require that a cell line be validated prior to publication of any data [Lorenzi et al. 2009].
The problem is not limited to lines that are misidentified only, but rather included lines that are cross contaminated [Bubenik et al. 2000, Drexler et al. 2003, Lacroix et al. 2008, Cape Davis, Macleod et al. 1999]. Cross-contaminated lines may account for as low as fourteen percent to thirty percent of the cell used in a cohort [Masters et al. 2002, Roland M. Nardone et al. 2007]. Another significant characteristic is the long-term stability of the lines as it should be observed while pre-screening cell culture. Without these three critical elements (correct identification, stability, and non-contamination) any diagnostic value acquired and/or gained from using cell culture, as a prototype should be considered as false claim as discussed in many retrospective studies [Harris et al 1981, Nelson Rees et al. 1981,
Buehring et al. 2004, Melcher 2005, Masters et al. 2005, Liscovitch et al. 2006, Chatterjee et al. 2007].
There exist several laboratory methods for the detection of inter-species and intra- species identification [Masters 2001, Drexler et al. 1999, Melcher et al. 2005]. Different DNA analysis can be used in detecting intra-species cross-contamination [Gignac et al. 1993, Roland M. Nardone et al. 2007-2008, Lacroix et al. 2008]. In recent years, short tandem repeat (STR) is been increasingly used by many researchers to replace some of these DNA analyses for the detection of the intra-species cross-contamination [Masters et al. 2001, Muller et al. 2004, Parson et al. 2005]. Although several studies evaluated the statistical extent of the STR profile in use to determine paternity, in the identification of individual in the US, and in identifying potential suspects in a crime scene, and showed that STR profiles have a positive impact. With the new Power Plex 16 core CODIS system the statistical power of using the STR profile to identify an individual is equal to the odd chances of one match in trillion individuals [FBI 2008]. Thus, Short Tandem Repeat (STR) technology has enabled fast DNA informative data for further analysis. These data however, need to be carefully interpreted and analyzed when verifying identities and stability of the lines in culture.
Similarity techniques, also known as matching techniques, are measures of identifications that are heavily used in artificial intelligence such as classification networks, and fuzzy logic [Weaver et al. 1949, Miller 1955, Watson et al. 1980, Gibbs et al. 2002]. The goal is to acquire a similarity score between a target feature vector and a source or reference feature vector. Apart from the string matching techniques that may involve ontology driven techniques, most of the matching techniques adopted on genomic data is based on the concept of the probability of shared information to the overall information. There are a wide variety of similarity coefficients that are frequently used with genomic data these are based on either the use of Euclidean, Jaccard, Cosine, or Pearson Correlation among others. Correlation Techniques often assume linear relations among variables. Jaccard and Cosine coefficients for example are best to use with asymmetric data versus Euclidean coefficient that is rather good with symmetric data [Steele et al. 1986]. A recent study compared the late methods among other similarity coefficients are in the work of [Meyer et al. 2004].
The Entropy methods of matching are of broader scenario from the mutual information (MI) theory that also accounts as one of the basic techniques adopted in many researches for making a decision of identification off of two feature vectors that are extracted to measure the inter- independency between two random variables [Thomas et al. 1991, Nimenman et al. 2002]. The entropy represents measure of uncertainty for the amount of information a random variable contains, which is based on the probability mass function. Therefore, in most cases entropy can perform less efficiently than MI that takes into account all information dependencies in all clusters i.e. with low or high score [Meyer et al. 2004].
High throughput data such as the genomic data, however, may often inherit different conceptual information and a single method can hardly capture the broad complex definitions that often exist in genomic coded data. Therefore, the resultant feature vector from one application to another are often different, thus, requires different adaptations to an existing method or new method. Thus far, none of these methods take into considerations the conceptual mechanism that often exists in such genomic coding, like the coded values that exist in the STR profile data i.e. while accounting for all the complex interpretations that these coded alleles are bound to.
Cell line misidentification is been the talk of the hour for few years. The subject is very important to the scientific community. American Type Culture Collection Standards Development (ATCC) Organization Workgroup ASN-0002, and other companies such as Genetica, Human Cell line Authentication Laboratory offer services to authenticate Cell Lines using the STR DNA Technology, which is the method mostly used to certify the identity of cell lines. However, this authentication does not include evaluating the stability and cross-contamination of cell lines. The process of obtaining a DNA profile using highly polymorphic STR repeats includes the isolation of DNA from a cell line specimen, the amplification of targeted areas of that DNA, and analysis of the amplified product using a genetic analyzer.
Thus far, there exists no automated protocol in place for the translation of such complex STR profiles for matching cell lines. An automated system to detect miss-identified lines and to provide measures for the stability of the cell lines is a critical step to bring any further in-vitro investigation to a solid ground. Embodiments of the invention include one or more of the following: 1 -Integrate the short tandem repeat profile of one or more target lines into a DBMS architecture designed to host the high throughput STR profile data generated. 2-
Match each target STR profile to other Reference STR profiles, and disseminate if the line is an alias to one or more pre-recognized lines from the public and/or private repositories that are housed in a DBMS. 3-The system uses the high throughput STR profile data to measure the stability of the line and inspect the potential of any cross-contamination. If the target line is deemed unique and stable a reference profile of the line is then created, thus, allowing a dynamic platform for future matching. An Enterprise In-House Laboratory Information Management System (LIMS) is designed and implemented to allow scientists to periodically observe the long-term stability of the line and to detect miss-identified lines prior to using the line to selectively experiment any targeted molecule interactions. Thus, an embodiment of the invention is an automated logic that detects identity, instability, and potential of cross- contamination in a cell line.
An automated system to detect the miss-identified lines and provide decision support for the stability of the line is a critical step to bring any further in-vitro investigation to a solid ground. The system described here is deemed successful to deliver solutions to the cell culture dilemma. Embodiments of the system detects if a target line is mis identified, detects if the target line is stable or unstable, and checks if there is probable cross-contamination from an extracted pattern. This pattern represents the proportional changes in the complex tetra- nucleotides alleles. If the target line is deemed unique and stable a reference profile of the line is then created thus allowing dynamic platform for future matching. The system offers a dynamic, robust, and efficient detection of misidentified, and/or unstable cell lines, and saves a great deal of scientists' time and efforts as it ensures the delivery of eradicated and disseminated information to the scientific community.
BRIEF SUMMARY OF THE INVENTION
Embodiments of the invention relate to a system, method and a computer readable medium embodied by a computer readable medium, which detects misidentified, unstable and cross-contaminated cell lines using STR profiles.
An embodiment of the invention is a method comprising a) obtaining a sample comprising a target cell line; b) testing the sample to determine an STR profile of the target cell line; c) receiving information corresponding to a STR profile of a reference cell line; d) comparing the alleles in the STR profile of the target to the alleles in the reference STR profile and calculating a probable identity match; e) calculating probable instability of the
target cell line; calculating probable cross-contamination of the target cell line; f) and generating an output comprising the probable identify match of the target cell line with the cross reference cell line, g) the probable instability of the target line, and the probable cross- contamination of the target line. The method may additionally comprise identifying the target cell line as the reference cell line or a different cell line from the probable identity match. The method may also additionally comprise receiving information corresponding to the instability of the reference cell line and the potential of cross-contamination of the reference cell line. In a specific embodiment of the invention calculating the probable instability of the target cell line comprises equation 3. In another specific embodiment of the invention calculating the probable cross-contamination of the target line comprises equations 3 and 4. In a further embodiment of the invention, comparing the alleles in the STR profile of the target to the alleles in the reference STR profile and calculating the probable identify match comprises equation 1. Additionally, the information corresponding to the STR profile of the target cell may comprise a different number of STR loci than the corresponding number of STR loci that are available in the STR profile of the reference cell line. In another embodiment of the invention, the information corresponding to a STR profile of a reference cell line is received from an electronic database management system of STR profiles or other data set obtained from public STR profiles. Further, the output may be sent to an electronic database management system for storage. In an embodiment of the invention, the method is performed for each reference cell line in the database.
Another embodiment of the invention is a system for cell line discernment for use by a software application, the system comprising a processor in communication with a memory, where: the memory stores processor-executable program code; and the processor is configured to be operative in conjunction with the processor-executable program code to: a) receive information corresponding to an STR profile of a target cell line; b) receive information corresponding to an STR profile of a reference cell line; c) compare the alleles in the STR profile of the target to the alleles in the reference STR profile and to calculate a probable identity match; d) calculate probable instability of the target cell line; e) calculate probable cross-contamination of the target cell line; and f) generate an output comprising the probable identify match of the target cell line with the cross reference cell line, the probable instability of the target line, and the probable cross-contamination of the target line. The system may also include a server comprising an electronic database system that is interfaced with the processor. The embodiment may additionally comprise identifying the target cell
line as the reference cell line or a different cell line from the probable identity match. Further, additional information may be received corresponding to the instability of the reference cell line and the potential cross-contamination of the reference cell line. In a specific embodiment of the invention the calculating instability comprises calculating the probable instability of the target cell line with equation 3. In another specific embodiment of the invention, calculating the probable cross-contamination of the target line comprises equations 3 and 4. In another embodiment of the invention, the comparing the alleles comprise calculating the probable identity match with equation 1. Further, comparing the information corresponding to the STR profile of the target cell may comprise a different number of STR loci than the corresponding number of STR loci that are available in the STR profile of the reference cell line. In a further embodiment, the additional information is received from an electronic database management system of STR profiles or other data set obtained from public STR profiles. Additionally, the output generated may sent to an electronic database management system. The embodiment may additionally comprise repeating the executable embodiments of the invention for each reference cell line in the electronic database management system.
Another embodiment of the invention is a computer readable medium comprising computer-usable program code executable to perform operations comprising: a) receiving information corresponding to a STR profile of a target cell line; b) receiving information corresponding to a STR profile of a reference cell line; c) comparing the alleles in the STR profile of the target to the alleles in the reference STR profile and to calculate a probable identity match; d) calculating probable instability of the target cell line; e) calculating probable cross-contamination of the target cell line; and f) sending the probable identify match of the target cell line with the cross reference cell line, the probable instability of the target line, and the probable cross-contamination of the target line. The computer readable medium may additionally comprise identifying the target cell line as the reference cell line or a different cell line from the probable identity match. The computer readable medium may also additionally comprise receiving information corresponding to the instability of the reference cell line and the potential cross-contamination of the reference cell line. In a specific embodiment of the invention calculating the probable instability of the target cell line comprises equation 3. In another specific embodiment of the invention, calculating the probable cross-contamination of the target line comprises equations 3 and 4. In a further specific embodiment of the invention comparing the alleles comprises equation 1. The
information corresponding to the STR profile of the target cell may comprise a different number of STR loci than the corresponding number of STR loci that are available in the STR profile of the reference cell line. The computer readable medium may additional comprise interfacing with an electronic database system. In another embodiment of the invention the receiving information corresponding to a STR profile of a reference cell line receives information from an electronic database management system of STR profiles or other data set obtained from public STR profiles. In another embodiment of the invention, the data is sent to an electronic database management system. The computer readable medium may further comprise executing the method for each reference cell line in the database. The foregoing has outlined rather broadly the features and technical advantages of the present invention in order that the detailed description of the invention that follows may be better understood. Additional features and advantages of the invention will be described hereinafter, which form the subject of the claims of the invention. It should be appreciated by those skilled in the art that the conception and specific embodiment disclosed might be readily utilized as a basis for modifying or designing other structures for carrying out the same purposes of the present invention. It should also be realized by those skilled in the art that such equivalent constructions do not depart from the spirit and scope of the invention as set forth in the appended claims. The novel features which are believed to be characteristic of the invention, both as to its organization and method of operation, together with further objects and advantages will be better understood from the following description when considered in connection with the accompanying figures. It is to be expressly understood, however, that each of the figures is provided for the purpose of illustration and description only and is not intended as a definition of the limits of the present invention.
BRIEF DESCRIPTION OF THE DRAWINGS For a more complete understanding of the present invention, reference is now made to the following descriptions taken in conjunction with the accompanying drawing, in which:
Fig. 1 shows the association of global hit score, instability score, and potential of cross-contamination score including proper dependencies affecting each score and associated rules.
Fig. 2 illustrates a 16 locus analysis using hit score, instability score and potential of cross-contamination score showing the alignment STR channels for target line A and an instant reference line B.
Fig. 3 illustrates the reference and target locus and the hit score, instability score and potential of cross-contamination score, and associated embodiments of mathematical equations.
Fig. 4 is a schematic flow which illustrates the steps of an embodiment of the algorithm.
Fig. 5 illustrates STR object models and their relations as used in an embodiment of the system.
Fig. 6 illustrates a summary of the top-down work flow and example input files.
Fig. 7 is an embodiment of the system architecture illustrating different architecture layers and elements.
Fig. 8 is schematic of the layers within an example system architecture showing a) a cross section of the system architecture and b) a layout of an embodiment of the network.
Fig. 9 is a screenshot of a) a log in screen and b) a STR enterprise web interface.
Fig. 10 is a screenshot of a STR LIMS submission panel.
Fig. 1 lis a screenshot of a STR LIMS self-hit results panel.
Fig. 12 is a screenshot from the STR LIMS of a target/reference preview panel. Fig. 13 is a screenshot of a sample STR LIMS instant STR profile search.
Fig. 14 is a screen shot of the sample queue enterprise LIMS showing sample submission.
Fig. 15 is the linguistic cluster results of the target line C8161-C9.
Fig. 16 is a table of selected target lines that are highly matched to other references in the DBMS.
Fig. 17 is a table of selected target lines highly matched to other references in the DBMS and were marked as unstable.
Fig. 18 is a table of selected target line (HCT1 16-B) marked as unstable and probable mixed with reference lines SKOV3 and ARTFLOX. Fig. 19 is a table of the sensitivity (se) and specificity (sp) analysis for the instability and probable mix detection evaluation that is prepared using pre-selected 77 cultures.
Fig. 20 is a linguistic cluster of a selected target line AU565.
Fig. 21 is a table of selected target lines detected as alias to other references and marked as stable. Fig. 22 is a table of four selected target lines that were marked as unstable.
Fig. 23 is a table of selected target lines matched with other reference lines and a) marked as alias and b) partial match, unstable and probable mix.
Fig. 24 is graph analysis of the HCT116-A and HCT116-B family of cultured cells.
Fig. 25 is graph analysis for line AU565 hits with other cultured lines. Fig. 26 is graph analysis of the CAKI1 family with other cultured lines.
Fig. 27 is graph analysis of the IGROV1 family with other cultured lines.
Fig. 28 is a table of the collection of thirty seven cultures that were mixed and not mixed and used for probable mix testing.
Fig. 29 is a table of a linguistic cluster of a selected target lines CAKI1, HCT116, and TK-10 with selected reference lines.
DETAILED DESCRIPTION OF THE INVENTION
Various features and advantageous details are explained more fully with reference to the non-limiting embodiments that are illustrated in the accompanying drawings and detailed in the following description. Descriptions of well known starting materials, processing techniques, components, and equipment are omitted so as not to unnecessarily obscure the
invention in detail. It should be understood, however, that the detailed description and the specific examples, while indicating embodiments of the invention, are given by way of illustration only, and not by way of limitation. Various substitutions, modifications, additions, and/or rearrangements within the spirit and/or scope of the underlying inventive concept will become apparent to those skilled in the art from this disclosure.
In keeping with long-standing patent law convention, the words "a" and "an" when used in the present specification in concert with the word comprising, including the claims, denote "one or more." Some embodiments of the invention may consist of or consist essentially of one or more elements, method steps, and/or methods of the invention. It is contemplated that any method or composition described herein can be implemented with respect to any other method or composition described herein.
System Flow
The genotype STR profile collected is very useful to inspect identity, monitor long term stability of a line, and can be used to infer mutation or mixed lines being observed either across different cultures and/or across large family of passages in an assorted culture. The following section is used to outline an embodiment of the method and system flow.
A 16 -core CODIS system may be used. Table 1 illustrates a 16 CODIS generated (i.e. STR profile). Other systems may be used such as FFFL, GammaSTR, and AmpF/STR systems.
As an example, the two STR profiles to be matched are produced by the same platform to exclude any potential shifts per locus and that may exist with some loci due to different platforms such as the difference between Power Plex 1.1, and Power Plex 16 for some of the loci, and that all other precautions in handling low DNA are preserved and taken according to all guidelines. In another embodiment of the invention, different systems may be matched. Reference lines may come from any STR databases such as ATCC, DSMZ, RIKEN, JCRB and NIST, in addition to institutional databases.
Target line STR profiles are examined against reference STR line profiles stored as references in the DBMS creating a node in the target linguistic cluster being inspected. A bank of the 16 channels STR profiles is prepared for both target and reference line being tested for a matching hit; the system first inspects the influence Ia, and K for target line and reference line respectively. Therefore, off the 16 channels, each channel is first examined locally producing, overall, three scores per local channels that is total of 16x3 = 48 local scores. The first 16 scores represents the local hit scores between the target line (A) and the reference line (B). The second 32 scores represent the target line, where 16 scores indicates stability of 16 local locus followed by 16 score inspecting cross-contamination for the target line (A). The first 16 local hit scores are then weighted to calculate the global weighted hit scores for the target (A). The mathematical notations are described in the following section. Matching Algorithm μ, τ, and v serve as three indicators for the global hit score, stability score, and potential of cross-contamination score respectively for a line. Fig. 1 shows the association of the three measures used to describe a line during an instant hit between a target line and a reference line. The figure simplifies the inclusions of each score with regard to the other two. (I) Misidentification
In an embodiment using a 16-core locus, a single amplification routine generates coded repeats pattern, which are used as inputs to the Algorithm. The Algorithm then uses these sixteen channels, aiming at three objectives namely; 1- Calculate global weighted score between a target line A and a reference line B mathematically, noted as μ. 2- State if a line is stable or unstable based on a weighted stability score noted as τ. 3- Inspect if a line expresses
potential of cross-contamination based on weighted score noted as v, as illustrated in Fig. 2. It is important to note that even though μ is considered as a property of the matching procedure between two STR profiles that belongs to two instant lines, two and three are properties for each individual line therefore, τ and v are independent from μ but vise versa is not true, that is, μ is affected by any instability mark the line may exhibit. Equation 1 and 2 are given below.
After aligning the 16 channels between a target STR (A) and a reference STR (B) for every locus, every allele in the target domain is inspected against every allele in the corresponding reference allelic domain to calculate μ; where, aj represents all target instant allele being visited and, thus, domain of all instances of circumstances are accounted
for. The notation used to extract locus's differentiated alleles μ; per locus is presented in equation 1 and 2. The value m represents all distinct alleles appeared in both target and corresponding reference allelic domain that are instantly being matched. Equation 1 depicts both absolute and conditional partial match while providing local score between finite sets A (target) and B (reference). The aj is the new pattern between instant pair a, and b allele as it defines the concise match without creating anarchy. The later is important to account for small aberrations between the default pair of alleles being matched and ensures they are accounted for and not ignored. For example target allele 9 matches reference allele 9.3 with 0.7 if reference lack instance of 9. Note that the notations used denote the cardinality of A and B sets.
A Global Hit score μ is then calculated and weighted with a mutual influence notated as ω = Ia / lb where, Ia<Ib is satisfied. It is important to note that the values Ia, and ¾, represent the inferred influence of the relevant profiles known as target and reference, respectively (i.e. number of available informative locus (i.e. fully matched, and partially matched). In addition,
the value Ia represents the inferred influence rather than the abstract influence across all loci with available information. Thus making ω defines the ratio between the inferred to the abstract influences.
To further illustrate the notion of mutual influence, assume we have matched 2 sets of profiles together and they resolved in set of local hit scores as { 1, 0.5, 0, 0} where there are available 3 loci with information for the target {(3), (1,3), (5), (0)}, and 3 loci for the reference {(3), (3), (0), (5)} . The min mutual Influence will be
and the maximum mutual influence will be ¾ = 0.75, whereas abstract influence is
Overall, the possible hit scores have to confer 2 ranges between minimum and maximum global scores. After calculating local scores for the available markers a minimum and maximum global hits are calculated and used to categorize the hits. In the case illustrated above, are
calculated as they define minimum global hit and maximum global hit.
(II) Instability Measures:
The clinical significance of agents that are pre-screened from stable cultured lines is thought of to be remarkable. It is expected that agents prescreened with stable cultures will introduce tremendous difference compared to agents that were prescreened from unstable or mixed cultures. Thus, introducing more reliable prompts to researchers. Therefore, a computerized decision support for culture instability is critical. Further, monitoring culture during modified protocol by itself is important. Cell cultures often exhibit instability, as characteristics are lost due to many reasons.
Although STR allelic profiles may exhibit compound repeats, every individuals is either homozygous or heterozygous i.e. showing inherited alleles from one parent or both parents. Despite the fact that these locus are chosen to represent alleles with high variations in a population, it does not change the reality that one individual cannot possibly represent all incidence of populated events. Therefore, a template filter is designed to extract how variant subsequent repeats are presented per locus STR. The filter is sensitive to the smallest variations and preserves heterozygosis events as the default. Equation 3 gives power not only to the variation in the allelic repeats but also to their spatial order from the first 2 alleles. The notation τ; depict the per-locus differentiated repeats as a measure of variation to highlight any instability per locus notated as τ; a mathematical description of τ; is in equation 3. As a result for every locus, ¾ : used to evaluate the pattern extracted to detect instability as notated
in equation (3). Stability is then preserved if
for all loci thus preserving heterozygosis events, and
maintains the order of the variant alleles where, the index j is the repeats order as its variation is being evaluated.
In order for a line to indicate potential of cross-contamination, then this line must have expressed instability already However, vice versa is not true by means that
stable line cannot be x-contaminated or even showing any sign of mixed (cross- contamination) phase i.e. loses its characteristics. Even though STR allelic profiles may exhibit compound repeats, every individual is either homozygous or heterozygous i.e. showing inherited alleles from one parent or both parents. Despite the fact that these locus are chosen to represent alleles with high variations in a population, one individual cannot possibly represent all incidence of populated events. Therefore, in one embodiment, a filter template is used to read how variant subsequent alleles are presented in every local STR given power not only to the variation in the allelic repeats but also to their spatial order from the first 2 alleles. The variation to highlight any instability per locus is notated as
Equation 3 is given below.
The later preserves heterozygosis events, and maintains the order of the variant alleles where, the index j is the allele order as its variation is being evaluated. For every locus used to evaluate the pattern extracted to detect instability as notated in
equation (3). These measures are again calculated for each locus separately. However,
is monitored based on incremental events represented as τ. While τ being incrementally evaluated until the first accumulation that brings the condition that marks
instability for the instant culture.
The threshold τ° is examined for each locus denoting an unstable line if τ° is reached from one or more locus. Thus far, τ° is defined empirically using several sets of data. The final decision for instability is then presented to the scientists as one of two categories and marked as stable versus unstable.
(III) Probable mixes
A cross-contaminated line can be caused from a variety of reasons. One example is if regulations in handling low DNA are not met. This can cause the scientist DNA to show in the STR profile leaving a cross-contamination call to that STR. Further, intra-species cross- contamination is expected particularly in handling cell culture.
Further, tetra-nucleotide alleles are associated with certain incomplete polymorphisms that appear as x.2, x.3 in the pattern. Some STR alleles are compound tetra-nucleotides. If a line expressed frequent fragments of incomplete polymorphisms that are forming series of compound tetra-nucleotides per locus that raises the possibility that it captures sign of cross- contaminations thus loosing pattern of solid and discrete repeats.
However, even with this consideration one filter to capture variations in compound repeats per locus is not enough as not all locus are designed to capture compound repeats. Therefore, equation 3 and equation 4 are examined for any indication of probable mix allowing sensing all variations that may be presented with simple and/or compound repeats in a locus. Thus, probability of having potential mixed lines will increase if a target will show sign of instability from equation 3 where
alone with or without
Mathematically, represents a filter that collects the proportional variation in patterns that is associated with the complex repeats rather than the simple repeats in the set of local alleles. Where yj denotes the fraction in a compound repeats whereas xj is the whole compound repeats appears as x.2. The notation
further indicates potential of cross- contamination (see equation 4).
Mathematically, it can be proved that cross-contamination score can be derived from a stability score i.e.
(i.e. v entails τ). A filter to collect the ratio between variations in the irrational domain of an allele to the variation in its rational domain was designed to indicate cross-contamination score as follows in Equation 4:
Therefore, if one or more locus were detected unstable as τ; extracted from equation 3, are closely re-examined (re-visited for all locus) with equation 4 to evaluate probable cross-contamination. Equation 4 extract patterns are associated with the complex alleles rather than the simple alleles in the set of local alleles. Furthermore this pattern examines the proportional change in consecutive complex alleles.
To explain one case of such scenario, if an instant target set of alleles are found as { 11, 12, 13, 19, 23.2, 23.3} then will be reached specifying the line as unstable but
moreover, will be equal to (l)2. In another example, with allelic set as { 1 1, 12, 13, 19, 9.3, 23.2} , is equal to (0.0071942446)2 where, The lower found the higher
probability that the tetra nucleotide fragments presented comes from another line raising the probability of a mixture by one at this locus. Exactly as with τ, however, τ depicts the consecutive variations for every 2 consecutive compound alleles in the set that is denoting either simple or complex repeats. Furthermore, in one embodiment τ assumes heterozygosis as a default. The symbol yj denotes the fraction in the complex allele if there exists complex alleles while xj represents the complex allele.
Overall, and due to the fact that not all loci are designed for compound repeats but rather simple repeats, the final decision is inferred from both equation 3 and 4 and the number of loci expressing both phenomena. Given the assumption that for a line to indicate potential of probable mix then this line must of expressed instability already
however, vise versa is not true by means that stable line will not express such potential as previously, discussed in the hypothesis section.
Studies have shown that if an individual is exposed to genetic mutations, the homeostatic mechanisms that maintain regularity of paired structures tend to break down. That is, the tendency to maintain relatively stable equilibrium between interdependent elements will break down. Thus, cell lines will represent at time point certain state. Differentiation of this state is more likely to have occurred on a faster pace compared to other lines that were subjected to modification or compared to same cell type and origin that is still harbored within its natural host. However, it takes a long observation and recurrent experiments to determine this hypothesis. Each STR profile is then considered exhibiting instability if τ; > τ° for one or two locus. As far as a mixed line, instability will be disrupting
the majority of the 16 loci thus disrupting the normal state of homozygosis and/or heterozygosis structure indicating potential of mixed lines being involved in the STR pattern. On the other hand, contamination may imply that germ and/or viral line is involved which is far from being detectable by the STR patterns and other experiments are designed to validate typical type of contamination. The relationship between these equations is illustrated in Fig. 3.
Algorithm
An embodiment of the algorithm or procedure is illustrated in Fig. 4 and given below. After initializing the pivot of elements for an instant target line A, which in part, either determines if the current detectable target is informative enough to include or exclude and, thus, flag it as empty vector of target loci. STR object model has been designed to reproduce all elements needed to describe an instant profile with 16 pre-known loci in a map where each node in the map has a set of alleles (see Fig. 5). After initializing a linguistic cluster list mapped by the target, and initializing STR target object and for every pulled instant reference in the list also initialize STR reference object. An instant pivot to match the target to all references in a collected list of STR references that are called back from the DBMS works as follows: Line one initiates and starts locus procedure. Line 2-5 initializes parameters needed for every locus at a time for local target hit for both objects instantly, for the local instability, and local cross-contamination where, all being examined at a time. Line seven states that for each allele being observed in the current locus between target and reference while the stated conditions are met, line 8-1 1 marks the equations and their conditions for the STR Hit Match.
Line 13-14 lists the instant calculation to scale the severity of cross-contamination if conditions stated are met for the calculation. According to line 16 if there exists more than a pair alleles in the instant locus being examined then line 17 evaluates the differential change that is also being accumulated. However, another parameter τ° is incremented as τ; > τ°. Thus, the target is marked as unstable if 1 or more locus is detected meeting τ° value. Line 18-23 is calculated after all alleles in the Loci domain are examined. Specifically, line 22-23 calculates Local hit and updates the global hit μ. The above steps are repeated until all 16 loci are exhausted. Line 25-27 equations establish the weight. Line 28 computes weighted global hit. A distance value between target and instant reference is then calculated in line 29. Finally, line 30-31 pushes the final linguistic cluster between the target A and the instant reference STR object Bk while updating the linguistic cluster map for this Target A as the
routine is being exhausted over all references collected that ends with a return of the linguistic cluster for a target line A as LCA.
The system allows public STR profiles generated by other facilities to be matched, which may enable smart logic while mining lines from different laboratories. In order to reason the behavior of the cell line model that was subjected to different modification, and/or may been altered due to different handling technique from different laboratories and over the long run (certain period of time) of characterizing a line, a panel was designed to allow scientists to retrieved matched lines that were subjected to varieties of modification and/or panel of lines collected from across culture i.e. different passages. System framework
Certain units described in this specification have been labeled as modules, in order to more particularly emphasize their implementation independence. A module is "[a] self- contained hardware or software component that interacts with a larger system." Alan Freedman, "The Computer Glossary" 268 (8th ed. 1998). A module comprises a machine or machines executable instructions. For example, a module may be implemented as a hardware circuit comprising custom VLSI circuits or gate arrays, off-the-shelf semiconductors such as logic chips, transistors, or other discrete components. A module may also be implemented in programmable hardware devices such as field programmable gate arrays, programmable array logic, programmable logic devices or the like. Modules may also include software-defined units or instructions, that when executed by a processing machine or device, transform data stored on a data storage device from a first state to a second state. An identified module of executable code may, for instance, comprise one or more physical or logical blocks of computer instructions, which may be organized as
an object, procedure, or function. Nevertheless, the executables of an identified module need not be physically located together, but may comprise disparate instructions stored in different locations which, when joined logically together, comprise the module, and when executed by the processor, achieve the stated data transformation. Indeed, a module of executable code may be a single instruction, or many instructions, and may even be distributed over several different code segments, among different programs, and across several memory devices. Similarly, operational data may be identified and illustrated herein within modules, and may be embodied in any suitable form and organized within any suitable type of data structure. The operational data may be collected as a single data set, or may be distributed over different locations including over different storage devices.
In the following description, numerous specific details are provided, such as examples of programming, software modules, user selections, network transactions, database queries, database structures, hardware modules, hardware circuits, hardware chips, etc., to provide a thorough understanding of the present embodiments. One skilled in the relevant art will recognize, however, that the invention may be practiced without one or more of the specific details, or with other methods, components, materials, and so forth. In other instances, well- known structures, materials, or operations are not shown or described in detail to avoid obscuring aspects of the invention. The system for cell line discernment using STR profiles may include a server, a data storage device, a network, and a user interface device. In a further embodiment, the system may include a storage controller, or storage server configured to manage data communications between the data storage device, and the server or other components in communication with the network. In an alternative embodiment, the storage controller may be coupled to the network.
In one embodiment, the user interface device is referred to broadly and is intended to encompass a suitable processor-based device such as a desktop computer, a laptop computer, a Personal Digital Assistant (PDA), a mobile communication device or organizer device having access to a network. In a further embodiment, the user interface device may access the Internet to access a web application or web services hosted by the server and provide a
user interface for enabling a user to enter or receive information. For example, the user may enter text files such as the one shown in Fig. 6 into the system.
The network may facilitate communications of data between the server and the user interface device, the server and a database, and/or a user and a database. The network may include any type of communications network including, but not limited to, a direct PC to PC connection, a local area network (LAN), a wide area network (WAN), a modem to modem connection, the Internet, a combination of the above, or any other communications network now known or later developed within the networking arts which permits two or more computers to communicate, one with another. STR profiles against reference cell lines STR profiles to determine the identity, stability, and probable of cross-contamination of the target cell lines. Additionally, the server may access data stored in a data storage device via a Storage Area Network (SAN) connection, a LAN, a data bus, or the like. In another embodiment, a desktop computer is configured to analyze target cell lines STR profiles against reference cell lines STR profiles to determine the identity, stability and probable of cross-contamination of the target cell lines.
The data storage device may include a hard disk, including hard disks arranged in an Redundant Array of Independent Disks (RAID) array, a tape storage drive comprising a magnetic tape data storage device, an optical storage device, or the like. In one embodiment, the data storage device may store reference STR profiles and/or target cell line STR profiles. The data may be arranged in an electronic database and accessible through Structured Query Language (SQL) queries, or other database query languages or operations. A data storage device may be located on a desktop type computer, on a server, and/or on a dedicated database device.
The server may host a software application configured for cell line discernment using STR profiles. In another embodiment, a desktop type computer may host a software application configured for cell line discernment using STR profiles. The software application may further include modules for interfacing with the data storage devices, interfacing with a network, interfacing with a user, and the like. In a further embodiment, the server or desktop type computer may host an engine, application plug-in, or application -programming interface (API). In another embodiment, the server may host a web service or web accessible
software application. The software application may be embodied on a computer readable medium.
An example of the system architecture is shown in Fig. 7. A user can interact with the system through multiple web-based interfaces. STR profiles of target cell lines are input into the system through an input file and/or directly from a STR profile-generating machine. The STR target cell line profiles are compared to reference cell line STR profiles that are located in one or more database. The process of inputting files and comparing the STR profiles may be automated. The automation may occur at the system that hosts the database, in the users system, or in a system that bridges the user and the database. The input files may include one or more STR profile.
Another example of the system architecture is shown in Fig. 8. Fig. 8a illustrates the layers of the architecture. The exterior of the circle represents user interaction with the system such as users from lab groups or institutions who interacts with the system through a web browser and/or client APIs, for example. The interactions with the system may be through an enterprise Laboratory Information Management System (LIMS) and/or web portal. Scheduled processes layer which interface with one or more database handles intermediate transactions that do not need an interface with users but rather needs layer of automated scheduled processes. Web services layer handles all users and scheduled processes requests in connections with DBMS layers or AI layers. Fig. 8b illustrates an example of a layout of a cancer network. Cell profiles may come from any type of cell line in the cancer network.
Fig. 9 is an example of a system graphical user interface. Once a user logs into the system, the user can perform various STR related activities such as submitting an STR profile, matching an STR profile or searching for an STR profile. In this example, the user submits a STR profile and identifies where the profile comes from, the species the cell line belongs to, the technology used to generate the STR profile, the lab ID, the set ID and a description. After submitting the STR profile, the user can preview self -hit results and preview the STR client sample details. Multiple target reference cell STR profiles can be submitted in batch. This is shown in more detail in Fig. 10-12. Fig. 10 illustrates submitting STR information for a target cell, as described above. Fig. 11 shows a view of the preview of STR client self hit results with target cell lines matched against multiple reference cell lines. The self hit results show the identity of the target and reference cell lines, the min hit %, the
average distance, the hit score, and the probable stability. Fig. 12 shows a screen shot of the target/reference preview window and/or panel. This window illustrates the STR profiles of the input target cells STR profiles, in addition to the hit results. Fig. 13 is a screenshot showing an instant search of STR samples. The example search is based on any STR with completed analysis. Fig. 14 is an illustration showing a multiple sample submission. The sample names, cell lines, tissue type, state, species and STR profiles are shown. Fig. 6, as discussed above, is an example STR profile input file generated from GeneMapper and another format of input file that is used or been transformed to, and the system top-down flow to output the match, stability and potential cross-contamination. EXAMPLE 1
Exemplary Materials and Data Set
Cell lines From the National Cancer Institute (NCI) were used. A similar panel of the sixty cell lines' STR profiles was described in Lorenzi et al. 2009. In addition, forty-seven cultured cell lines, some of which are well known to be unstable, were used to test the instability results. Further, a set of forty-nine cell lines were mixed with other cell cultures making different mixture for the purpose of evaluating the detection system in detecting probable cross-contamination. List of used cultures that were pre-prepared by an expert biologist for the purpose of testing the mixed lines are shown in Fig. 28.
A 16 core CODIS system is used. Table 1 represents the 16 CODIS generated (i.e. STR profile) for two lines. DNA was extracted from cell lines using QiaAMP maxi preps (Qiagen) and DNA was quantified by Nanodrop spectrometry. Cell line profiles obtained by STR DNA fingerprinting using the AmpFfSTR identifier kit according to manufacturer instructions (Applied Biosystems cat 4322288).
American Type Culture Collection Standards Development (ATCC) Organization Workgroup ASN-0002, and other companies offer services to authenticate Cell Lines using the STR DNA Technology which is the method mostly used to certify the identity of cell lines. The STR profiles were compared to known ATCC fingerprints, to the Cell Line Integrated Molecular Authentication (CLIMA) database version 0.1.200808, (Nucleic Acids Research 37:D925-D932 PMCID: PMC2686526), and to the complete 16 loci and to the MD Anderson fingerprint database where the STR profiles either matched known DNA fingerprints or were unique.
Further, the lines HCT-1 16 (NCI-60), and IGROV-1, X, Y were obtained as part of the NCI-60 cell line collection from the National Cancer Institute. Furthermore, OVCAR429, OVCAR432, OVACR433, and HCT-1 16 cell lines were provided by several labs. NCI-60 cell lines were grown in RPMI [Lorenzi et al. 2009, Stults et al. 2011], 10% FBS in a Co2 incubator. Same line is introduced in Savas et al. 2011 as it was obtained from the ATCC Organization Workgroup.
EXAMPLE 2
Exemplary results
The true positive fraction of the system was tested based on expert biologists and their definition of instability. Further, set of prepared mixtures were designed by expert biologist for the purpose of evaluating the instability and probable mix decisions placed, by the detection schema. A total of 77 cultures were used in this study. Amongst those a panel was prepared with 49 mixture lines which, were used to evaluate the se and sp of cross contamination detection. The culture profiles were presented to two expert biologists for the determination of instability those were compared to the system detection result. Overall, total of 77 lines were used for the se and sp analysis. Fig. 19 shows the se and sp analysis for the instability detection and probable cross contamination.
The detection schema offers robust solution to the misidentification, instability, and probable cross contamination. Overall, the true positive fraction (TPF) for the instability detection was 0.98 at a false positive fraction (FPF) of 2/55. The positive predictive value (ppv) for the instability detection is 0.962 and the negative predictive value (npv) for the instability detection is 0.667. Since the instability result is evaluated across the 16 loci with at least one locus is found unstable to determine if the line is unstable, the instability detection is influenced by subjective assessment. To assess the performance of the system in detecting the pre-prepared probable mixed lines, the probable mix detection performance was analyzed. Thus, for the 49 mixed lines pre-prepared by biologists the TPF was 0.957 at FPF 3/49. The positive predictive value (ppv) for the probable mix detection is 0.937.
Fig. 15 shows an example of matching results showing selected Target Line (C8161- C9) as it hits several reference lines. Fig. 16 shows matched cell lines where the target was misidentified. Fig. 17 shows examples of target lines that were highly matched an instant reference and both found unstable in a, and b examples. Fig. 18 shows a case of cross
contamination. The HCT1 16-A and HCT116-B lines were also identified as unstable lines. However, HCT1 16-B was marked as probable cross-contamination in the system while proposing 2 other lines (SKOV3 and ARTFLOX) as potential of being cross contaminated with HCT1 16-B. Unlike HCT116-A the system detected the culture as unstable. Fig. 29 shows an example of matching results showing selected target line (CAKI1) as it hits several reference lines.
EXAMPLE 3
Exemplary Discussion
A matching technique is presented based on pattern extracted from the STR profiles. The system offers a dynamic, and efficient detection of misidentified, and unstable cell lines. The system saves a great deal of scientists' time and efforts thus ensuring the welfare of the research as it ensures the delivery of disseminated information to the scientific community. The robust performance of the detection techniques suggests that the system provides an efficient and more advanced starting point for the automated cell cultures validation and authentication efforts that may or may not currently, exist for the detection of misidentification.
Thus far, cell line plays as the only available modeling tool in monitoring the efficacy of anti-cancer agents. For this purpose, good practice is to observe the long-term stability of a line. Overall, it is expected that such model may behave slightly different than it would behave inside its natural host "the human body" for diseased and/or free of disease host as opposed to a line that are cloned and preserved in a media. For example, it is expected that the lost of default characteristics of a model line may be carried faster opposed to the lost of default characteristics of a similar characterized cell line living in its natural host. In conclusion, the only means to enable further studies to test such hypothesis is to have such a cost effective system in place that deliver fast, effective, and reliable mechanisms to decimate barriers that may otherwise obscure and/or overburden the translation of such large high throughput genomic data.
EXAMPLE 4
Various STR profiles of cell lines from a variety of sources were analyzed using the methods described above. Fig. 20 shows the linguistic cluster of the target cell line AU565
matching results. The matching results show the reference and target cell lines identification, as well as the min hit, min distance, weighted hit and the stability of the reference cell line.
Fig. 21 are selected target lines detected as alias to reference lines and marked as stable. The matches are shown to be exact, exact, and exact with difference for the three cultures displayed. Fig. 22 shows the results of four target lines that were marked as unstable and Fig. 23 shows target lines that were matched with reference lines and marked as alias (a) and (b) partial match, unstable, and probably mix. Figs. 24-27 are graph analysis of families of cultured lines namely: HCT116, AU565, CAKI1, and IGROV.
Although the present invention and its advantages have been described in detail, it should be understood that various changes, substitutions and alterations can be made herein without departing from the spirit and scope of the invention as defined by the appended claims. Moreover, the scope of the present application is not intended to be limited to the particular embodiments of the process, machine, manufacture, composition of matter, means, methods and steps described in the specification. As one of ordinary skill in the art will readily appreciate from the disclosure of the present invention, processes, machines, manufacture, compositions of matter, means, methods, or steps, presently existing or later to be developed that perform substantially the same function or achieve substantially the same result as the corresponding embodiments described herein may be utilized according to the present invention. Accordingly, the appended claims are intended to include within their scope such processes, machines, manufacture, compositions of matter, means, methods, or steps.
REFERENCES
All patents and publications cited herein are hereby incorporated by reference in their entirety herein.
American Type Culture Collection (ATCC) Standards Development Organization Workgroup ASN-0002. Nature Reviews 10(6): 441-8, 2010.
An Open Letter, July. The Discovery Center for Cell and Molecular Biology, The Catholic University of America, 2007.
Boonstra, J.J., et al. J. Natl Cancer Inst., Vol. 102(4):271-274, 2010.
Bubenik, J., et al, Br. J. Cancer 44:692, 1981.
Bubenik, J., et al, Intl. J. Cancer, 1 1, 765; 1973.
Butler, J.M. National Institute of Standards and Technology, 2005.
Butler, J.M. Academic Press, London. 2001.
Butler, J.M. Elsevier Academic Press, Lib. Of Cong. Cat. No. 2004114214; 2005.
Butler, J.M. Elsevier Sciences, Feb. 2005.
Butler, J.M., et al. J. Forensic Sci. 48: 1054-1064. 2003.
Butler, John et al. NIST Biotechnology Division, 1997 '.
Capes-Davis, A. et al. Intl J. of Cancer. Jul, 127(1): 1-8, 2010.
Chatterjee, R. Science 315:928-931, 2007.
Clayton, T.M. et al. Forensic Science International 76, 7-15, 1995.
Cover, Thomas M. et al. John Wiley & Sons, Inc., 1991.
Defendi, V., R.E. et al. J. Natl. Cancer Inst. 25:359-385, 1960.
Dirks, W.G., et al. ALTEX. 22(2): 103-9, 2005.
Dirks, W.G., et al, Methods Mol Med. 2004;88:43-55, 2004.
Drexler, H.G., et al. Journal oj 'the Leukemia Socety of 'America. Feb. 17(2):416-26, 2003
Drexler, H.G., et al. Journal of the Leukemia Socety of America, Oct. 13(10): 1601-7, 1999.
Evans, W.E., et al. Science. 286:487, 1999.
Faivre, L., et al. Orphanet Encyclopedia, April 2005.
Federal Bureau of Investigation, "CODIS-NDIS Statistics.", 2008.
Gartler, S. M. NCI Monogr. 26: 167-195, 1967.
Gartler, S.M. Nature (Lond.), 217:750-751, 1968.
Gibbs, A., et al. Inter- national Statistical Review, 70, 419-436. 2002.
Gignac, S.M. et al. Leukemia Lymphoma, Jul, 10(4-5):359-68, 1993.
Gilbert, D.A. et al. American Journal of Human Genetics. 47:499-514, 1990.
Hagelberg et al. Nature 352, 427-429, 1991.
Harris, E.D., et al. (1981). Arch. Biochem. Biophys., 206, 105; 1981.
Kasai et al. Journal of Forensic Sciences, 35, 1196-1200. 1990.
Lacroix, M. . J. Cancer. 122: 1-4, 2008.
Lorenzi, Philip L., et al. Mol. Cancer Ther. 2009; 8(4). April 2009.
Lucey, D.J. et al. Laryngoscope. 116: 161-162. 2006.
Lucy, B. P. et al. Archives of Pathology and Laboratory Medicine. 133 : 1463-1467, 2009.
Masood, E. Nature 398:545^16, 1999.
Masters, J.R. Nat. Rev. Cancer 2002. 2:3 15-319, 2002.
Matsuo, Y. et al. Journal of Human Cell Research Society, 1999, Sep., 12(3): 149-54, 1999. Meyer, Andreia de Silva et al. Genet. Mol. Biol, vol.27 no.1, San Pao, 2004
Miller, G. , Information Theory in Psychology II-B (pp. 95-100), 1955
Muller, P., et al. The American Society of Biochemistry and Molecular Biology Inc. 2004.
Nardone, R.M. Department of Biology , Catholic University of America, Washington, DC.
Nardone, R.M. BioTechniques 45:221-227, 2008.
Nardone, R.M. Cell Biology Toxicol (23 :367-372). 2007.
Nelson, M.S. et al. Journal of Forensic Sciences, 47:345-349, 2002.
Nelson-Rees, W.A. et al. Science \U: 1093-1096, 1974.
Nelson-Rees, W.A., et al. Science. 209:719-720, 1980.
Nelson-Rees, W.A., et al, Science. 212:446-452, 1981.
Nelson-Rees, W.A., et al. Science. 195: 1343-1344, 1977.
Nemenman, I., et al. Advances in Neural Information Processing, MIT Press, Vol. 14, 2002.
NIH Pub. No. 91-3260, March 1999.
Puers, C. et al, American Journals of Human Genetics, 53, 953-958; 1993.
Ringel, M. D. J Clin Endocrinol Metab 93(1 1): 4226-7, 2008.
Sarkar, N., et al. Forensic Science Intl. 128: 196-201, 2002.
Savas, S., et al. (201 1) NCI60 PLoS ONE 6(4): el 8306. doi: 10.137 l/journal.pone.OO 18306. Schweppe, R. E., et al. J Clin Endocrinol Metab 93(1 1): 4331-41, 2008.
Steele, J. Annals of Statistics, 14, 753-758, 1986.
Stults Dawn M., et al. BMC Molecular Biolog)- 201 1, 12-23, 201 1.
Tanabe H., et al. Tissue Culture Res. Commun., 18:329-338, 1999.
Vogel C, et al. Mol Syst Biol. 6:400, 2010.
Watson, G. New York: Wiley, 1980.
Weaver, W., et al. Urbana: University of Illinois Press, 1949.
Wickelgren, I. Science. 285:998, 1999.
Claims
CLAIMS 1. An in vitro method comprising:
obtaining a sample comprising a target cell line;
testing the sample to determine an STR profile of the target cell Line; receiving information corresponding to a STR profile of a reference cell line; comparing the alleles in the STR profile of the target to the alleles in the
reference STR profile and calculating a probable identity match; calculating probable instability of the target cell line; calculating probable cross-contamination of the target cell line; and generating an output comprising the probable identify match of the target cell line with the cross reference cell line, the probable instability of the target line, and the probable cross-contamination of the target line.
2. The method of claim 1 , additionally comprising identifying the target cell line as the reference cell line or as an alias to a different cell line from the probable identity match.
3. The method of claim 1, additionally comprising receiving information corresponding to the instability of the reference cell line and the potential of cross-contamination of the reference cell line.
4. The method of claim 1, wherein calculating the probable instability of the target cell line comprises equation 3.
5. The method of claim 1, wherein calculating the probable cross-contamination of the target line comprises equations 3 and 4.
6. The method of claim 1 , wherein comparing the alleles in the STR profile of the target to the alleles in the reference STR profile and calculating the probable identify match comprises equation 1.
7. The method of claim 1, wherein the information corresponding to the STR profile of the target cell comprises a different number of STR loci than the corresponding number of STR loci that are available in the STR profile of the reference cell line.
8. The method of claim 1, wherein the information corresponding to a STR profile of a reference cell line is received from an electronic database management system of STR profiles or other data set obtained from public STR profiles.
9. The method of claim 1, wherein the output is sent to an electronic database management system for storage.
10. The method of claim 8, wherein the method is performed for each reference cell line in the database.
1 1. A system for cell line discernment for use by a software application, the system comprising a processor in communication with a memory, where: the memory stores processor-executable program code; and the processor is configured to be operative in conjunction with the processor-executable program code to:
a) receive information corresponding to a STR profile of a target cell line;
b) receive information corresponding to a STR profile of a reference cell line; c) compare the alleles in the STR profile of the target to the alleles in the reference STR profile and to calculate a probable identity match;
d) calculate probable instability of the target cell line;
e) calculate probable cross-contamination of the target cell line; and
f) generate an output comprising the probable identify match of the target cell line with the cross reference cell line, the probable instability of the target line, and the probable cross-contamination of the target line.
12. The system of claim 11, additionally comprising identifying the target cell line as the reference cell line or a different cell line from the probable identity match.
13. The system of claim 1 1, wherein the additionally comprising receiving and/or calculating additional information corresponding to the instability of the reference cell line and the potential cross-contamination of the reference cell line.
14. The system of claim 11, additionally comprising calculating the probable instability of the target cell line with equation 3.
15. The system of claim 11, wherein calculating the cross-contamination comprises calculates the probable cross-contamination of the target line with equations 3 and 4.
16. The system of claim 1 1, comparing the alleles comprises calculating the probable identify match with equation 1.
17. The system of claim 1 1, wherein the information corresponding to the STR profile of the target cell comprises a different number of STR loci than the corresponding number of STR loci that are available in the STR profile of the reference cell line.
18. The system of claim 11, additionally comprising a server connected to the processor comprising an electronic database system wherein the processor receives reference cell line information from the database system.
19. The system of claim 1 1, wherein the output module is configured to send data to an electronic database management system.
20. The system of claim 18, additionally comprising executing the system of claim 1 for each reference cell line in the electronic database management system.
21. A computer readable medium comprising computer-usable program code
executable to perform operations comprising:
a) receiving information corresponding to a STR profile of a target cell line; b) receiving information corresponding to a STR profile of a reference cell line;
c) comparing the alleles in the STR profile of the target to the alleles in the reference STR profile and to calculate a probable identity match;
d) calculating probable instability of the target cell line;
e) calculating probable cross-contamination of the target cell line; and f) sending the probable identify match of the target cell line with the cross reference cell line, the probable instability of the target line, and the probable cross-contamination of the target line.
22. The computer readable medium of claim 21, additionally comprising identifying the target cell line as the reference cell line or a different cell line from the probable identity match.
23. The computer readable medium of claim 21, additionally comprising receiving information corresponding to the instability of the reference cell line and the potential cross-contamination of the reference cell line.
24. The computer readable medium of claim 21, wherein calculating the probable instability of the target cell line comprises equation 3.
25. The computer readable medium of claim 21, wherein calculating the probable cross-contamination of the target line comprises equations 3 and 4.
26. The computer readable medium of claim 21, wherein the comparing the alleles comprises equation 1.
27. The computer readable medium of claim 21, wherein the information corresponding to the STR profile of the target cell comprises a different number of STR loci than the corresponding number of STR loci that are available in the STR profile of the reference cell line.
28. The computer readable medium of claim 21, wherein receiving information corresponding to a STR profile of a reference cell line is configured to receive information from an electronic database management system of STR profiles or other data set obtained from public STR profiles.
29. The computer readable medium of claim 21, wherein outputing sends data to an electronic database management system.
30. The computer readable medium of claim 28, further comprising executing the method for each reference cell line in the database.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201161525793P | 2011-08-21 | 2011-08-21 | |
US61/525,793 | 2011-08-21 |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2013028699A2 true WO2013028699A2 (en) | 2013-02-28 |
WO2013028699A3 WO2013028699A3 (en) | 2013-05-02 |
Family
ID=47747069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2012/051746 WO2013028699A2 (en) | 2011-08-21 | 2012-08-21 | Cell line discernment using short tandem repeat |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2013028699A2 (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014074611A1 (en) * | 2012-11-07 | 2014-05-15 | Good Start Genetics, Inc. | Methods and systems for identifying contamination in samples |
WO2018150378A1 (en) * | 2017-02-17 | 2018-08-23 | Grail, Inc. | Detecting cross-contamination in sequencing data using regression techniques |
WO2019005877A1 (en) * | 2017-06-27 | 2019-01-03 | Grail, Inc. | Detecting cross-contamination in sequencing data |
WO2019241913A1 (en) * | 2018-06-19 | 2019-12-26 | 深圳华大基因科技有限公司 | Digital identification generating method, device and system and storage medium |
US12006533B2 (en) | 2017-02-17 | 2024-06-11 | Grail, Llc | Detecting cross-contamination in sequencing data using regression techniques |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090068646A1 (en) * | 2004-10-22 | 2009-03-12 | Promega Corporation | Methods and kits for detecting mutations |
WO2010129793A1 (en) * | 2009-05-06 | 2010-11-11 | Ibis Biosciences, Inc. | Methods for rapid forensic dna analysis |
-
2012
- 2012-08-21 WO PCT/US2012/051746 patent/WO2013028699A2/en active Application Filing
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090068646A1 (en) * | 2004-10-22 | 2009-03-12 | Promega Corporation | Methods and kits for detecting mutations |
WO2010129793A1 (en) * | 2009-05-06 | 2010-11-11 | Ibis Biosciences, Inc. | Methods for rapid forensic dna analysis |
Non-Patent Citations (3)
Title |
---|
BARALLON, R. ET AL.: 'Recommendation of short tandem repeat profiling for authenticating human cell lines, stem cells, and tissues' IN VITRO CELL. DEV. BIOL. ANIM. vol. 46, no. 9, 08 July 2010, pages 727 - 732, XP055066088 * |
ELTONSY, N. ET AL.: 'Detection algorithm for the validation of human cell lines' INT. J. CANCER. vol. 131, no. 6, 12 April 2012, pages E1024 - E1030, XP055066086 * |
NIMS, R. W. ET AL.: 'Short tandem repeat profiling: part of an overall strategy for reducing the frequency of cell misidentification' IN VITRO CELL. DEV. BIOL. ANIM. vol. 46, no. 10, 07 October 2010, pages 811 - 819, XP055066087 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014074611A1 (en) * | 2012-11-07 | 2014-05-15 | Good Start Genetics, Inc. | Methods and systems for identifying contamination in samples |
WO2018150378A1 (en) * | 2017-02-17 | 2018-08-23 | Grail, Inc. | Detecting cross-contamination in sequencing data using regression techniques |
US12006533B2 (en) | 2017-02-17 | 2024-06-11 | Grail, Llc | Detecting cross-contamination in sequencing data using regression techniques |
WO2019005877A1 (en) * | 2017-06-27 | 2019-01-03 | Grail, Inc. | Detecting cross-contamination in sequencing data |
WO2019241913A1 (en) * | 2018-06-19 | 2019-12-26 | 深圳华大基因科技有限公司 | Digital identification generating method, device and system and storage medium |
US11822629B2 (en) | 2018-06-19 | 2023-11-21 | Bgi Shenzhen Co., Limited | Method and apparatus for generating digital identity and storage medium |
Also Published As
Publication number | Publication date |
---|---|
WO2013028699A3 (en) | 2013-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dorrity et al. | The regulatory landscape of Arabidopsis thaliana roots at single-cell resolution | |
Zoffmann et al. | Machine learning-powered antibiotics phenotypic drug discovery | |
Kotopka et al. | Model-driven generation of artificial yeast promoters | |
Bloom et al. | Finding the sources of missing heritability in a yeast cross | |
Auer et al. | Statistical design and analysis of RNA sequencing data | |
Faure et al. | DiMSum: an error model and pipeline for analyzing deep mutational scanning data and diagnosing common experimental pathologies | |
Vincent et al. | Next-generation sequencing (NGS) in the microbiological world: How to make the most of your money | |
Lohmueller et al. | Proportionally more deleterious genetic variation in European than in African populations | |
Brettschneider et al. | Quality assessment for short oligonucleotide microarray data | |
Selmecki et al. | Polyploidy can drive rapid adaptation in yeast | |
Navarro et al. | Chromosomal speciation and molecular divergence--accelerated evolution in rearranged chromosomes | |
Conant et al. | Asymmetric sequence divergence of duplicate genes | |
Noor et al. | Consequences of recombination rate variation on quantitative trait locus mapping studies: simulations based on the Drosophila melanogaster genome | |
Mather et al. | CADD score has limited clinical validity for the identification of pathogenic variants in noncoding regions in a hereditary cancer panel | |
Heyse et al. | Coculturing bacteria leads to reduced phenotypic heterogeneities | |
Martiniuk et al. | Impact of commercial strain use on Saccharomyces cerevisiae population structure and dynamics in Pinot Noir vineyards and spontaneous fermentations of a Canadian winery | |
JP2005512175A (en) | A method for identifying genetic features of complex genetic classifiers | |
WO2013028699A2 (en) | Cell line discernment using short tandem repeat | |
US20230141128A1 (en) | Molecular technology for predicting a phenotypic trait of a bacterium from its genome | |
Baker | The changes that count | |
Yin et al. | Validation of preimplantation genetic tests for aneuploidy (PGT-A) with DNA from spent culture media (SCM): concordance assessment and implication | |
Logsdon et al. | A novel variational Bayes multiple locus Z-statistic for genome-wide association studies with Bayesian model averaging | |
Eyre et al. | Prediction of minimum inhibitory concentrations of antimicrobials for Neisseria gonorrhoeae using whole-genome sequencing | |
Daw Elbait et al. | A population-specific major allele reference genome from the United Arab Emirates population | |
Rapti et al. | CoverageMaster: comprehensive CNV detection and visualization from NGS short reads for genetic medicine applications |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 12825048 Country of ref document: EP Kind code of ref document: A2 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 12825048 Country of ref document: EP Kind code of ref document: A2 |