Classifications

• • 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
  • G16B50/00—ICT programming tools or database systems specially adapted for bioinformatics
• • 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
  • G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
• • 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
  • G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
  • G16B30/10—Sequence alignment; Homology search
• • 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
  • G16B50/00—ICT programming tools or database systems specially adapted for bioinformatics
  • G16B50/10—Ontologies; Annotations

Definitions

• the present disclosure relates in part to the computerized analysis of sequencing data. More particularly, the present disclosure relates in part to the computerized process of identifying and analyzing genome modifications such as transgene insertion sites.
• transgene flanking sequences may be needed for the commercialization and registration of products that contain transgene sequences.
• the identification and characterization of transgene flanking sequences may also be important for other types of activities, like characterization of events generated by EXZACTTM Precision Technology brand genome modification technology. For example, EXZACTTM Precision
• ZFNs zinc finger nucleases
• EXZACTTM brand technologies can be used to generate ZFN-promoted double strand breaks within the genome of an organism, thereby resulting in the targeted insertion of transgenes at a specific loci of interest in a DNA sequence.
• the transgene flanking sequence consists of a chromosomal flanking region of the genomic integration site and the integrated transgene.
• the transgene flanking sequences may contain deletions, inversions, or insertions which result from the integration of the transgene into a specific location of the chromosome. Regions of nucleic acid similarity may exist between the transgene DNA, the cloning vector used in sequencing, primers and/or adapters used to isolate the transgene flanking region sequence, the chromosomal sequence in which the transgene has integrated, and other unrelated DNA fragments which have been inserted into the genome via unexpected rearrangements.
• Various methods can be used to isolate a transgene flanking region sequence. This transgene flanking region sequence can then be sequenced using traditional dideoxy sequencing methods, chain termination sequencing methods, or via Next Generation Sequencing methods.
• DNA sequence analysis can be used to determine the nucleotide sequence of the isolated and amplified fragment.
• the amplified fragments can be isolated and sub-cloned into a vector and sequenced using chain-terminator method (also referred to as Sanger sequencing) or Dye-terminator sequencing.
• the amplicon can be sequenced with Next Generation Sequencing.
• NGS technologies do not require the sub-cloning step, and multiple sequencing reads can be completed in a single reaction.
• Three NGS platforms are commercially available, the Genome Sequencer FLX from 454 Life Sciences/Roche, the Illumina Genome Analyser from Solexa and Applied Biosystems' SOLiD (acronym for:
• tSMS Single Molecule Sequencing
• SMRT Single Molecule Real Time se menting
• the Genome Sequencer FLX which is marketed by 454 Life Sciences/Roche is a long read NGS, which uses emulsion PCR and pyrosequencing to generate sequencing reads. DNA fragments of 300 - 800 bp or libraries containing fragments of 3 -20 kbp can be used. The reactions can produce over a million reads of about 250 to 400 bases per run for a total yield of 250 to 400 megabases. This technology produces the longest reads but the total sequence output per run is low compared to other NGS technologies.
• the Illumina Genome Analyser which is marketed by Solexa is a short read NGS which uses sequencing by synthesis approach with fluorescent dye-labeled reversible terminator nucleotides and is based on solid-phase bridge PCR. Construction of paired end sequencing libraries containing DNA fragments of up to lOkb can be used. The reactions produce over 100 million short reads that are 35 - 76 bases in length. This data can produce from 3 - 6 gigabases per run.
• the Sequencing by Oligo Ligation and Detection (SOLiD) system marketed by Applied Biosystems is a short read technology.
• This NGS technology uses fragmented double stranded DNA that are up to 10 kbp in length.
• the system uses sequencing by ligation of dye-labeled oligonucleotide primers and emulsion PCR to generate one billion short reads that result in a total sequence output of up to 30 gigabases per run.
• tSMS of Helicos Bioscience and SMRT of Pacific Biosciences apply a different approach which uses single DNA molecules for the sequence reactions.
• the tSMS Helicos system produces up to 800 million short reads that result in 21 gigabases per run. These reactions are completed using fluorescent dye-labeled virtual terminator nucleotides that is described as a 'sequencing by synthesis' approach.
• the SMRT Next Generation Sequencing system marketed by Pacific Biosciences uses a real time sequencing by synthesis. This technology can produce reads of up to 1000 bp in length as a result of not being limited by reversible terminators. Raw read throughput that is equivalent to one-fold coverage of a diploid human genome can be produced per day using this technology.
• a high-throughput method is needed to confirm that a transgene is integrated into the genome, and for identifying the specific chromosomal location of a transgene, if inserted through random integration or targeted to a site specific locus via homologous recombination.
• a flexible, high-throughput transgene flanking sequence analysis system is provided to analyze sequence data and define transgene insertion sites within the genome of an organism.
• the method includes steps to identify and annotate the transgene and the transgene flanking sequence, including the chromosomal flanking sequence, within a contiguous DNA fragment of, for example and without limitation, a complete genome.
• the analysis system contains, in an embodiment, a graphical user interface, an analysis pipeline, and a summary display for input sequences.
• the present disclosure includes a method for analysis.
• the method comprises: electronically receiving sequence data, electronically receiving one or more reference data sequences related to at least an expression vector, associating the sequence data with at least one of the reference data sequences to identify a transgene flanking sequence, searching a genome for one or more insertion sites of the transgene flanking sequence, and annotating the genome and the one or more insertion sites within the genome when one or more insertion sites are found.
• the reference data is further related to at least one primer. In a further embodiment of any of the above embodiments, the reference data is further related to at least one adapter. In a further embodiment of any of the above embodiments, the reference data is related to at least a primer and an adapter. In a further embodiment of any of the above embodiments, the reference data is further related to at least one cloning vector. In a further embodiment of any of the above embodiments, the reference data is further related to a right cloning vector and a left cloning vector.
• the reference data is further related to at least one of a left cloning vector, a primer, an adapter, a right cloning vector, and a transgene expression vector sequence.
• the reference data is further related to a cloning vector, a primer, and an adapter. In another further embodiment of any of the above embodiments, the reference data is further related to a left cloning vector, a right cloning vector, a primer, and an adapter.
• the method further includes searching the sequence data for a first reference data sequence; and searching the sequence data for a second reference data sequence when said first reference data sequence is located.
• the first reference data sequence is selected from the group consisting of: an expression vector, an adapter, a primer, and a cloning vector sequence.
• the second reference data sequence is selected from the group consisting of: an expression vector, an adapter, a primer, and a cloning vector, sequence, the second reference data sequence being selected independently of the first reference data sequence.
• the first reference data sequence is an expression vector and the second reference data sequence is an adapter.
• the first and second reference data sequences are independently selected from the group consisting of: a primer and an adapter.
• associating the sequence data with the reference data sequence includes finding the exact sequence of the reference data sequence. In another further embodiment of any of the above embodiments, associating the sequence data with the reference data sequence includes finding the sequence within a margin of error of five percent of the base pairs in the reference data sequence.
• the present disclosure includes a system for analysis.
• the system includes a module for receiving sequence data, a module for receiving one or more reference sequences related to at least an expression vector, and a calculation module operable to associate the sequence data with at least one of the reference data sequences to identify a transgene flanking sequence, search a genome for one or more insertion sites of the transgene flanking sequence, and annotate the genome and the one or more insertion sites within the genome when the one or more insertion sites are found.
• the reference sequences are further related to at least one primer.
• the reference sequences are further related to at least one adapter. In a further embodiment of any of the above embodiments, the reference sequences are related to at least a primer and an adapter. In a further embodiment of any of the above embodiments, the reference sequences are further related to at least one expression vector sequence. In a further embodiment of any of the above embodiments, the reference sequences are further related to at least one cloning vector. In a further embodiment of any of the above embodiments, the reference sequences are further related to a right cloning vector and a left cloning vector. [0022] In a further embodiment of any of the above embodiments, the reference sequences are further related to at least one of a left cloning vector, a primer, an adapter, a right cloning vector, and an expression vector sequence.
• the reference sequences are further related to at least a cloning vector, a primer, and an adapter. In another further embodiment of any of the above embodiments, the reference sequences are further related to at least a right cloning vector, a left cloning vector, a primer, and an adapter.
• the computation module is further operable to search the sequence data for a first reference data sequence; and search the sequence data for a second reference data sequence when said first reference data sequence is located.
• the first reference data sequence is selected from the group consisting of: an expression vector, an adapter, a primer, and a cloning vector sequence.
• the second reference data sequence is selected from the group consisting of: an expression vector, an adapter, a primer, and a cloning vector sequence, the second reference data sequence being selected independently of the first reference data sequence.
• the first reference data sequence is an expression vector and the second reference data sequence is an adapter.
• the first and second reference data sequences are independently selected from the group consisting of: a primer and an adapter.
• associating the sequence data with the reference data sequence includes finding the exact sequence of the reference data sequence.
• associating the sequence data with the reference data sequence includes finding the sequence within a margin of error of five percent of the base pairs in the reference data sequence.
• Figure 1 A is an exemplary diagram showing a typical sequence which is produced, comprising a left cloning vector, a primer, a expression vector , a transgene flanking region sequence, an adapter, and a right cloning vector according to an embodiment of the present disclosure.
• Figure IB is an exemplary diagram showing a transgene insertion within the genome comprising an expression vector, a primer sequence and a transgene flanking region sequence inserted between sections of genome sequence according to an embodiment of the present disclosure.
• Figure 2A shows the flow of data and samples from sample input to the analysis system according to an embodiment of the present disclosure.
• Figure 2B shows a flow chart showing a method of data analysis according to an embodiment of the present disclosure.
• Figure 3 is a system diagram of a data analyzer according to an embodiment of the present disclosure.
• Figure 4 is a flow chart showing a method of data analysis according to an embodiment of the present disclosure.
• Figure 5A is a flow chart showing a flanking sequence identification processing sequence or method according to the flow chart of Figure 4.
• Figure 5B is a flow chart showing a method of identifying and marking a transgene flanking sequence.
• Figure 5C is a flow chart showing another embodiment of a method of identifying a transgene flanking sequence according to the flow chart of Figure 5 A.
• Figure 6 is an exemplary sequence according to an embodiment of the present disclosure.
• Figure 7 is an exemplary input screen of an identification system according to an embodiment of the present disclosure.
• Figure 8 is an exemplary output from the analysis system according to an embodiment of the present disclosure.
• Figure 9A is an exemplary screen showing the position of an expression vector, adapter, primer, and transgene flanking sequence.
• Figure 9B is an input sequence graphically identified in Figure 9 A.
• Figure 9C is a transgene expression vector 103 sequence graphically identified in Figure 9A.
• Figure 9D is an adapter sequence graphically identified in Figure 9A.
• Figure 9E is a primer sequence graphically identified in Figure 9A.
• Figure 9F is the genomic sequence flanking the transgene identified from the input sequence of Figure 9B.
• Figure 10 is an exemplary screen showing a transgene flanking sequence with a primer, but no right cloning vector.
• Figure 11 is an exemplary screen shot showing a transgene flanking sequence with an expression vector sequence, but no cloning vectors.
• Corresponding reference characters indicate corresponding parts throughout the several views. The exemplifications set out herein illustrate exemplary embodiments of the disclosure and such exemplifications are not to be construed as limiting the scope of the disclosure in any manner.
• transgene flanking sequences are discussed, the teachings herein may be applied to the analysis of other sequences.
• the systems and methods described may be applicable to output from any molecular method for identifying and characterizing transgene flanking sequences, and the systems and methods provide an automated way of locating the transgene insertion site or sites within a genome.
• the methods and systems also provide neighboring sequences and a local environment surrounding the insertion site, to determine if there are rearrangements in the local environment at or near the insertion site.
• An ideal isolated insertion sequence includes a left cloning vector 101, a primer 105, transgene flanking region sequence 107 transgene expression vector sequence 103, an adapter 109, and a right cloning vector 111.
• the left cloning vector 101 and right cloning vector 111 are parts of a cloning vector, which is a first sequence of DNA that a second sequence of DNA may be inserted into.
• the insertion of the second sequence of DNA divides the cloning vector into a right (3' portion) cloning vector 111 and a left (5' portion) cloning vector 101.
• the digestion of a cloning vector is completed by a restriction enzyme or via another method known in the art, thereby resulting in a cleaved DNA fragment.
• the digestion of the cloning vector at a single specific site generally yields a known left cloning vector 101 and right cloning vector 111 sequence.
• the insertion sequence inserted into a genome sequence is shown with respect to
• the expression vector 103 is a sequence that is used to introduce a gene into a target cell.
• a primer 105 is a short DNA sequence used to begin the process of DNA synthesis.
• the expression vector 103 is generally a sequence used for integration of a transgene into a genome.
• the transgene flanking region sequence 107 is the genomic sequence immediately upstream or downstream of the transgene insertion site; in the embodiment this sequence may either be known or unknown.
• An adapter 109 is a short oligonucleotide sequence which is ligated or annealed to the end of the transgene flanking sequence 107.
• the sequence of the adapter 109 is known, and is used to mark the end of the sequence and can also be used to amplify or sequence the unknown transgene flanking sequence 107.
• the transgene flanking sequence 107 consists of a chromosomal flanking region of the genomic integration site flanking the integrated transgene.
• the transgene flanking sequence may contain deletions, inversions, or insertions which result from the integration of the transgene into a specific location of the chromosome.
• the isolated sequence is ordered as a left cloning vector 101, a primer 105, an expression vector sequence 103, a transgene flanking region sequence 107, an adapter 109, and a right cloning vector 111, as illustrated in Figure 1A, however, the order of the sequence is not limited to those illustrated in Figures 1 A and IB.
• primer 105, expression vector 103, transgene flanking region sequence 107 are inserted into a genome sequence, and appear within the genome sequence.
• the adapter sequence is incorporated later as part of a method used to isolate the transgene flanking sequence.
• the resulting transgene flanking sequence as depicted in Figure 1 A is then subsequently analyzed using data analysis methods shown below.
• the sequences of the left cloning vector 101, the expression vector 103, the primer 105, the adapter 109, and the right cloning vector 111 are all known. In practice, one or more of the sections of the ideal sequence may be missing or may contain alterations.
• Figure 2A shows the flow of data and samples from sample input to the analysis system 207.
• Figure 2B shows a flow chart 220 showing a method of data analysis according to an embodiment of the present disclosure.
• input samples 201 are prepared with, for example and without limitation, a ZFN-initiated transgene insertion protocol.
• a ZFN-initiated transgene insertion protocol In the protocol, one or more portions of known sequences, such as a primer 105 or adapter 109, are added to a target genome whose sequence is also known.
• the samples may also be prepared by other methods of transgene insertion.
• the transgene insertion process creates modified sequences, with insertions at one or more sites in the genome.
• An exemplary modified sequence is provided in Figure IB.
• one or more sequencers 205 generate sequence data from one or more input samples 201.
• the sequencers 205 determine the transgene flanking region sequence which is used to identify the location of the insertion in the genome, and confirm the specific sequence of the transgene insertion.
• the sample data in the embodiment, is in the form of one or more text files including sequence data.
• the input samples 201 are loaded into a sequencer 205 according to a protocol or operating instructions of the sequencer 205. For example, a Solexa ILLUMINA brand
• the sequencer 205 generates data related to the sequences 201.
• the data may include, but is not limited to, one or more text files, Standard Flowgram Format ("SFF") or similar files, images files, or other data files containing information related to the sequences of the DNA strands in the input samples 201.
• the sequence information also includes confidence data, so that each base in a sequence may have a confidence interval associated with it, or each sequence has a confidence interval associated with it.
• the confidence interval is a mathematical calculation calculated by the sequencer, and may include the strength of the read of the particular base by the sequencer 205. In one illustrative example, the confidence interval is an integer from one to nine.
• a confidence interval of one indicates that the sequencer 205 has relatively low confidence that the base reported was the base in the DNA strand.
• a confidence interval of nine indicates that the sequencer 205 has relatively high confidence that the base reported was the base in the DNA strand.
• the sequencer 205 also reports other information in addition to the confidence interval. For example, the sequencer 205 may report when a base could not be read.
• the data from the sequencer 205 is provided to the analysis system 207.
• the data is provided by a network or a dedicated connection between the sequencer and the analysis system 207, or by a removable storage from the sequencer to the analysis system 207.
• the sequencer prints the data to a screen or to a printer, and the data is input into the analysis system 207 from, for example and without limitation, a keyboard or a scanner.
• the analysis system 207 is a part of the sequencer.
• the reference sample information 203 is transmitted to the analysis system 207.
• the reference sample information 203 may include, but is not limited to, the sequences of the left and right cloning vectors, which may be provided as a single sequence, the expression vector 103, the primer 105, and the adapter 109.
• the sequence information in an embodiment, is transferred to the analysis system 207 via a network.
• the reference sample information 203 is transmitted to the analysis system 207 with the sequence information from the sequencers 205.
• the analysis system 207 receives the sequence data from the one or more sequencers 205, and analyzes the sequence data, as described more fully below.
• the analysis system 207 also takes reference sample data 203 as an input.
• the reference sample data 203 may include, for example and without limitation, sequence information of the adapter 109, the primer 105, the left 101 and/or right cloning vectors 111, the expression vector 103, or the target genome sequence information.
• the entire target genome sequence data is provided to the analysis system 207.
• a subset of the entire target genome sequence is provided to the analysis system 207.
• the analysis system 207 sends a request for all or a portion of the target genome sequence to another system.
• the matched sequence data and other data produced by the analysis system 207 undergoes additional processing. Additional processing may include, but is not limited to, visualization, quantification, aggregation with data from other samples or other trials, or comparisons to a target genome sequence.
• the additional processing in an embodiment, is carried out by another system.
• the analysis system 207 carries out all or a portion of the additional processing. Additional processing is described below.
• Figure 3 shows a component view of the analysis system 207 according to an
• the analysis system 207 may include an input module 303, a calculation module 305, an output module 307, and a visualization module 311, which, in an embodiment, reside in memory 315 of the analysis system 207.
• the modules may be executed by a controller 325 of analysis system 207.
• the controller 325 is one or more processors, and the controller 325 includes operating system software to control access to the controller 325 and the memory 315.
• the memory 315 includes computer readable media.
• Computer-readable media may be any available media that may be accessed by one or more processors of the analysis system 207 and includes both volatile and non- volatile media. Further, computer readable-media may be one or both of removable and non-removable media.
• computer-readable media may include, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, Digital Versatile Disk (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which may be used to store the desired information and which may be accessed by analysis system 207.
• the analysis system 207 may be a single system, or may be two or more systems in communication with each other.
• the analysis system 207 includes one or more input devices, one or more output devices, one or more processors, and memory associated with the one or more processors.
• the memory associated with the one or more processors may include, but is not limited to, memory associated with the execution of the modules, and memory associated with the storage of data.
• the analysis system 207 is associated with one or more networks, and communicates with one or more additional systems via the one or more networks.
• the modules may be implemented in hardware or software, or a combination of hardware and software.
• the analysis system 207 also includes additional hardware and/or software to allow the analysis system 207 to access the input devices, the output devices, the processors, the memory, and the modules.
• the modules, or a combination of the modules may be associated with a different processor and/or memory, for example on distinct systems, and the systems may be located separately from one another.
• the modules are executed on the same system as one or more processes or services.
• the modules are operable to communicate with one another and to share information.
• the functions of two or more modules may instead be executed in the same process, or in the same system.
• the input module 303 receives data from an input device 301.
• the input module 303 may also receive data over a network from another system.
• the input module 303 receives one or more signals from a computer over one or more networks.
• the input module 303 receives data from the input device 301, and may rearrange or reprocess the data into a format recognizable by the calculation module 305, so that the data may be interpreted by the calculation module 305.
• the input device 301 may, in an embodiment, be a client 304, which a user interacts with to send signals to and receive signals from the analysis system 207.
• the client 304 may communicate with the analysis system 207 via one or more networks 302.
• the network 302 may include one or more of: a local area network, a wide area network, a radio network such as a radio network using an IEEE 802.1 lx communications protocol, a cable network, a fiber network or other optical network, a token ring network, or any other kind of packet-switched network may be used.
• the network 302 may include the Internet, or may include any other type of public or private network.
• the use of the term "network" does not limit the network to a single style or type of network, or imply that one network is used.
• a combination of networks of any communications protocol or type may be used. For example, two or more packet-switched networks may be used, or a packet-switched network may be in communication with a radio network.
• the input device 301 may communicate with the input module 303 via a dedicated connection or any other type of connection.
• the input device 301 may be in communication with the input module 303 via a Universal Serial Bus
• USB universal serial or parallel connection
• the transmission may also occur via one or more physical objects.
• the sequencer generates one or more files, and the sequencer or a user copies the one or more files to a removable storage device, such as a USB storage device or a hard drive, and a user may remove the removable storage device from the sequencer and attach it to the input module 303 of the analysis system 207.
• Any communications protocol may be used to communicate between the input device 301 and the input module 303.
• a USB protocol or a Bluetooth protocol may be used.
• the input device 301 is a sequencer.
• the sequencer analyzes one or more samples and generates sequence data regarding the one or more samples.
• the sequencer may communicate the sequence data to the input module 303 over a wireless or wired connection.
• the data is in the form of one or more files, or the sequencer may print the data to a screen or a printer, and the data is input into the analysis system 207 by, for example and without limitation, a keyboard, mouse, or scanner.
• the sequencer also includes additional data describing the samples.
• the calculation module 305 receives inputs from the input module 303, and executes one or more processing sequences based on the inputs. For example, and without limitation, the calculation module 305 receives sequence information and reference sample information for the sequences.
• Sample data includes the sequence information, for example and without limitation, the primer 105, the left and/or right cloning vectors 111, the expression vector 103, and/or the target genome.
• the sample data may be provided to the analysis system 207 by the user, by the sequencer, by a third party system, by another system associated with the analysis system 207, by a combination of two or more of these inputs or other suitable sources.
• the sample data may be provided to the analysis system 207 as a text file in a standard format.
• the text file may be formatted in the FASTA format.
• the sample data information may be input into the analysis system 207 by typing or pasting information into one or more text entry fields.
• the information may be formatted in the FASTA format, or another standardized format.
• other formats may be used.
• the Genbank® format may be used, or another format.
• the analysis system 207 may receive the sample data in a particular format, and may reformat the data to be further analyzed by the analysis system 207.
• the calculation module 305 applies one or more algorithms to identify the vector and/or adapter 109 within the input sequence, identify the orientation of the input sequence, locate the transgene flanking sequence within the input sequence, based on the vector and/or adapter 109 within the input sequence, if possible, receives the genome information related to the input sequence, and attempts to map the flanking sequence to the genome.
• the algorithms generate additional quantitative and qualitative data related to the input sequences. Additionally, in an embodiment, the input sequences are annotated and analyzed and/or visualized. The algorithms and processes used to identify and annotate input sequences are described with respect to the flow charts shown in Figures 4, 5A, 5B, and 5C.
• the calculation module 305 provides as an output, for example, data regarding the sequences and their position in a genome, and/or additional data to be used by a visualization module to visualize one or more of the sequences.
• the visualization module 311 receives data as input regarding the input sequences and the annotations from the calculation module 305.
• the visualization module 311 allows a user to visualize and/or manipulate the sequences and/or annotations.
• the visualization module 311 allows a user to visualize and/or manipulate the sequences and/or annotations.
• visualization module 311 may use Gbrowse, or a modified version of Gbrowse. Other sequence visualization software programs may be used in additional embodiments.
• a user may have the ability to manipulate a visual representation of the target sequences, or the target sequences and the genome.
• the visualization module allows the user to view the location of the target sequences in the genome, or the location of other sequences of interest within the genome.
• the visualization step allows a user to locate the target sequence within the genome and the location or changes to other sequences of the genome. This visualization may be helpful for providing an analysis of the transgene flanking sequence.
• the output module 307 receives an input, and transmits the input to an output device 309.
• the output module 307 receives the input from the calculation module 305, the visualization device 311 , or both the calculation module 305 and the visualization device 311.
• the received data may be in the form of alphanumeric data, and reformats the data to a format understandable to the output device 309, and transmits the data to the output device 309.
• the output module 307 and the output device 309 are in communication with one another.
• the output module 307 and the output device 309 is in communication via a network, or is in communication via a dedicated connection, such as a cable or radio link.
• the output module 307 may also reformat the data received from the calculation module 305 into a format usable by the output device 309.
• the output module 307 may create one or more files that may be read by the output device 309.
• the output device 309 is, in an embodiment, a visualization system, another data analysis system 207, or a data storage system.
• the output module 307 communicates with the output device 309 by transmitting one or more electronic files to the output device 309.
• the transmission may occur over a dedicated link, for example a USB connection or a serial connection, or may occur over one or more network connections.
• the transmission may also occur via one or more physical objects.
• the output module 307 may generate one or more files, and may copy the one or more files to a removable storage device, such as a USB storage device or a hard drive, and a user may remove the removable storage device from the analysis system 207 and attach it to the visualization system, another data analysis system 207, or the data storage system.
• Figure 4 shows a flow chart showing a method of data analysis according to an embodiment of the present disclosure.
• the samples are prepared according to one or more preparation protocols, and unknown samples are created with transgene insertions.
• the unknown samples are sequenced. Sequencing may occur according to a protocol or operating instructions of the sequencer. For example, a Solexa ILLUMINA brand sequencing machine or a Roche 454 brand sequencing machine may be used.
• the sequencer generates data related to the sequences.
• the data may include, but is not limited to, one or more text files or other data files containing information related to the sequences of the DNA strands in the samples.
• the sequence information also includes confidence data, so that each base in a sequence may have a confidence interval associated with it, or each sequence has a confidence interval associated with it.
• the confidence interval is a mathematical calculation calculated by the sequencer, and may include the strength of the read of the particular base by the sequencer.
• the confidence interval is an integer from one to nine.
• a confidence interval of one indicates that the sequencer has relatively low confidence that the base reported was the base in the DNA strand.
• a confidence interval of nine indicates that the sequencer has relatively high confidence that the base reported was the base in the DNA strand.
• the sequencer also reports other information in addition to the confidence interval. For example, the sequencer may report when a base could not be read.
• the data from the sequencer is input into the analysis system 207, and the system locates and identifies the flanking sequences in each of the sequenced input sequences. Flanking sequences may not be present in each of the input sequences, or the system may not be able to identify the location of a flanking sequence in an input sequence. Sequences where the flanking sequence is located and identified are noted by the system, and sequences where the flanking sequence is not located, or is located but not identified, are also noted by the system. The system generates output data based on the sequence data and the analysis conducted by the system. Exemplary analysis of sequence data is also described below with reference to Figures 5A-5C.
• FIG. 407 the system performs post-processing analysis on the sequence data and the flanking sequence location information as determined by the system.
• the sequence data, the target genome, and/or the flanking sequence location information may be visualized, qualitative measurements may be made with the data, and/or quantitative measurements may be made with the data.
• Figure 5A is a flow chart showing an exemplary method executed by analysis system 207 for flanking sequence identification.
• the expression vector 103 that is used as a part of the protocol to generate the input sequences is input into the system.
• one or more of the sequences for the right and left cloning vectors, the primer 105, and/or the adapter 109 are also provided.
• each of the sequences for the right and left cloning vectors, the primer 105, , and the adapter 109 are also provided.
• the sequences for the cloning vectors, the expression vector 103, the primer 105, and the adapter 109 are typically known, so that they can be identified and located within the genome.
• the information for the known sequences is input into the system to allow for identification of the sequences when compared to the input sequences.
• the input sequences are received from the sequencers or from one or more files.
• the one or more files may be transmitted to the system via, for example, a network, or may be provided to the system in another way.
• sequence information is received from the sequencers, it may be transmitted to the system via, for example, a network.
• the sequence information is in an electronic form that can be transmitted to the system and read by the system.
• the sequence information may, in an embodiment, include verification data or other additional data to ensure that the sequence information has not been corrupted or altered during transmission.
• the sequence information is stored in one or more databases, and the sequence information is transmitted from the one or more databases to the system via, for example, a network. Additionally, the genome information may be received from another database across a network.
• the genome information may be stored in a publicly accessible database, or a privately accessible database, and the genome information may be requested by the system, and the entire genome or a requested portion of the genome may be transmitted to the system based at least in part on the request.
• the analysis system 207 searches the input sequence for similarities with the known sequences including expression vector 103. If provided in step 501, the analysis system 207 may further search similarties with the cloning vectors, primer 105, and/or adapter 109 sequences. If one or more of these sequences is not provided in step 501, the analysis system 207 treats the sequence as not found. The analysis system 207 may use different search parameters to search for different sequences.
• the analysis system 207 may use a more stringent set of search parameters to identify the primer 105 and adapter 109, as they are shorter sequences and less likely to have been modified.
• the analysis system 207 may use comparatively less stringent search parameters to search for the other sequences in the input sequence, as they are longer and/or more likely to have been altered during the integration of the transgene into the genome.
• the analysis system 207 must find the exact sequence to identify the expression vector 103.
• the analysis system 207 identifies the expression vector 103 if the sequence for the expression vector 103 is found to within a margin of error. For example, the margin of error may be five percent of the base pairs in the expression vector 103 sequence.
• the analysis system 207 uses the LASTZ alignment program and algorithms to search for sequence similarity between the input sequence and the known sequences consisting of the cloning vector, transgene expression vector 103, primer 105, and/or adapter 109 sequences.
• the LASTZ program is described in Harris, R.S. (2007) Improved pairwise alignment of genomic DNA. Ph.D. Thesis, The Pennsylvania State University, the disclosure of which is hereby incorporated by reference in its entirety.
• the LASTZ program performs two kinds of sequence similarity searches.
• the first kind of sequence similarity search is an "exact search" which is a specific parameter setting of the LASTZ program.
• An “exact search” requires 95% identity, no gaps in the sequence, and at least 15 perfect character matches within the sequence.
• a scoring matrix is used to determine a "score" for the sequence, with the matrix including 1 for a match with the target sequence and -10 for mismatch with the target sequence.
• This search is used to identify the primer 105 and the adapter 109 within the input sequence if provided, since the primer 105 and adapter 109 in the input sequence are expected to be exactly the same as the primer 105 and adapter 109 sample sequences, as the primer 105 and adapter 109 sequences are short and therefore unlikely to have been modified during the experiment.
• the second kind of sequence similarity search is a "loose search.”
• the "loose search” does not have the same stringent requirements as the "exact search.”
• This search uses the default parameters for LASTZ, and is deployed for finding the transgene expression vector 103 and cloning vector sequence similarities in the input sequence.
• a “loose search” is used for the transgene expression vector 103 and cloning vector sequences, as they are longer and therefore more likely to have been modified during the experiment.
• Subsequences, within the input sequence, which share sequence similarity with a reference data sequence are labeled as a "type.”
• transgene expression vector 103 highly similar sequences between the input sequence and any of the selected primer 105 sequences are labeled or associated as the "primer 105 type.”
• the user selects 15 transgene expression vector 103 sequences to be included in the analysis and each has 30 homologies to subsequences within the input sequence, all 450 sequences will be associated with the type "transgene expression vector 103.”
• sequences that align with the highest levels of sequence similarity and alignment length to primer 105 sequences are classified as "primer 105 type.”
• sequences that align with highest levels of sequence similarity and alignment length to adapter 109 sequences are classified as "adapter 109 type.”
• sequence "type" is chosen arbitrarily from all of the tied sequences.
• the analysis system 207 searches the input sequence for the transgene expression vector 103 which shares the most sequence similarity. This search is conducted in one of two different ways, depending on whether or not a sequence similar to the primer 105 was identified. If a primer 105 sequence was identified in the input sequence, the best match containing the primer 105 is identified. In one embodiment, if the primer 105 was not provided in step 501 or identified in step 507, or none of the transgene expression vector 103 sequences contain a sequence which shares similarity with the "primer 105 type," the best overall match is considered and the transgene expression vector 103 with the highest sequence similarity is chosen. "Best overall match” in this context means choosing the match with the highest levels of sequence similarity and alignment lengths.
• transgene expression vector 103 Once the transgene expression vector 103 is located and identified, location and identification of the cloning vector sequence via sequence similarity alignments to known cloning vectors is attempted. Once a putative transgene expression vector 103 sequence is identified, the sequences upstream and downstream of this sequence are further characterized. The upstream cloning vector sequence is queried to identify cloning vectors which share sequence similarity at the start and end coordinates. The previously annotated sequences (transgene expression vector 103, primer 105, and adapter 109) are not queried. As such, the analysis system 207 searches all possible cloning vectors for sequence similarity with the region upstream from the previously identified feature. Then the analysis system 207 searches identified cloning vector sequence information for sequence similarity with the region downstream from the previously identified feature cloning vector in a similar manner. The vectors are identified by choosing the match with the highest levels of sequence similarity and alignment lengths.
• the orientation of the input sequence is identified, if possible.
• the analysis system 207 attempts to order input sequences in a left hand to right hand orientation; that is, with the 5' end of the sequence on the left side and the 3' end of the sequence on the right side.
• the sequencer may have sequenced the antisense strand of the DNA, in which case the sequence has to be reverse complemented.
• the system uses this information to identify and/or orient the input sequence. Orientation is determined by the location of the primer 105 and adapter 109 sequences. A forward orientation, wherein the primer 105 is located before the adapter 109 is preferred because of ease of visualization.
• FIG. 6 An example of an input sequence from the antisense strand is shown in Figure 6.
• the sequence of the primer 105 is known to the analysis system 207 as "TAAACA.”
• the analysis system 207 may initially not find either the primer 603 sequence in the input sequence 605.
• the analysis system 207 reverse complements the input sequence 605 to resolve a reverse
• the analysis system 207 system finds an exact match of the primer 603 to subsequences within the reverse complemented sequence 607.
• the analysis system 207 isolates the sequence 609 from the known primer 603, and proceeds with analysis of the reverse complemented sequence 607.
• the analysis system 207 instead compares reverse complemented sequences for the known primer 603 to the sequence 605, and, having identified the reverse complemented primer sequence 603, may reverse complement the entire sequence to yield a reverse complemented sequence 607, and may proceed with processing with the reverse complemented sequence 607.
• the transgene flanking sequence is located within the input sequence or the reverse complemented sequence, if the sequence was reverse complemented in the previous step. Exemplary location methods are described more fully with respect to Figures 5B and 5C.
• the transgene flanking sequence if found in the previous step, is located within the genome.
• the transgene flanking sequence is located in an integration site within the genome and is upstream or downstream of the transgene insertion site and contiguous with the expression vector sequence.
• the integration site is determined using a matching algorithm.
• BLAST Basic Local Alignment Search Tool
• the BLAST algorithm is described in Altschul S.F, et al., "Basic local alignment search tool.” J Mol Biol. 1990 Oct 5;215(3):403-10, the disclosure of which is hereby incorporated by reference in its entirety.
• the inputs for the BLAST search are the transgene flanking sequence and the genome.
• the BLAST search locates, if possible, the site or sites of integration of the transgene flanking sequence into the genome.
• the output of the BLAST search is a list of possible integration sites and a score for the fit. All masking and low complexity filtering is disabled for this homology search, to identify as many integration sites as possible.
• the output is parsed to find the top hit, which has the highest score for the fit. Once a top hit is identified, this region is considered to be the putative integration site of the transgene.
• linked endogenous upstream and downstream genes which are annotated in the genome are identified using a computer script.
• the input file of genome annotations is parsed, and the genes are indexed by chromosome and sorted by start coordinate.
• the system identifies the appropriate list of gene coordinates and performs a binary search to identify the correct insertion point for the integration site.
• the sorted list of coordinates for the transgene integration site will appear.
• the list is searched forward until a sequence greater than 10 kilobase pairs from the integration site is located. Then the list is searched backward until a sequence greater than 10 kilobase (kb) pairs from the integration site is located.
• the distance parameter can be varied, for example and without limitation, to >10 kb or ⁇ 10 kb of the integration site. Other ranges from the integration site may also be used.
• the analysis system 207 calculates the amount of overlap that exists between the chromosomal flanking sequence and any other sequence "types" used in any of the previously mentioned processes. This measure is calculated as the ratio of the number of bases in the input sequence similarity that are unique and not overlapped by any other sequence similarity ( unique _ bases ) and the total number of bases in the input sequence similarity ( total _ bases ). unique _ bases
• This ratio gives a quantitative value to the integration site.
• FIG. 9 A The annotated data from the previous boxes in Figure 5A may, in an embodiment, be presented for visual inspection in box 517. Examples of visualization are shown in Figures 9 A and 10. Additionally, the input sequence, the transgene flanking sequence, and/or additional information regarding the cloning vectors, the expression vector 103, the primer 105, the adapter 109, or the input sequence, is presented for visualization. Data regarding the transgene flanking sequence, the cloning vectors, the expression vector 103, the primer 105, the adapter 109, or the input sequence is also saved to one or more electronic files. [0089] Figure 5B is a flow chart showing a generalized method of marking a transgene flanking sequence 850.
• the expression vector 103 that is used as a part of the protocol to generate the input sequences is input into the system.
• one or more of the sequences for the right and left cloning vectors, the primer 105, the transgene expression vector sequence 103, and the adapter 109 are also provided.
• each of the sequences for the right and left cloning vectors, the primer 105, the transgene expression vector sequence 103, and the adapter 109 are also provided.
• the sequences for the cloning vectors, the expression vector 103, the primer 105, and the adapter 109 are typically known, so that they can be identified and located within the input unknown sequence.
• the information for the known sequences is input into the system to allow for identification of the sequences when compared to the input sequences.
• the input sequences are received from the sequencers or from one or more files.
• the one or more files may be transmitted to the system via, for example, a network, or may be provided to the system in another way. If sequence information is received from the sequencers, it may be transmitted to the system via, for example, a network.
• the sequence information is in an electronic form that can be transmitted to the system and read by the system.
• the sequence information may, in an embodiment, include verification data or other additional data to ensure that the sequence information has not been corrupted or altered during transmission.
• the sequence information is stored in one or more databases, and the sequence information is transmitted from the one or more databases to the system via, for example, a network.
• the genome information may be received from another database across a network.
• the genome information may be stored in a publicly accessible database, or a privately accessible database, and the genome information may be requested by the system, and the entire genome or a requested portion of the genome may be transmitted to the system based at least in part on the request.
• the analysis system 207 searches the input sequence for similarities with the known sequences including a first reference sequence, illustratively expression vector 103. If the expression vector 103 is not found in box 858, the method proceeds to box 860. The lack of expression vector 103 may indicate an error in the creation or the processing of the input sequence.
• the input sequence is marked as failing and is not matched against the genome. In an embodiment, the sequence is marked as red when the sequences are visualized.
• the method 850 proceeds to box 862.
• the analysis system 207 must find the exact sequence of expression vector 103 to proceed to box 862.
• the analysis system 207 may proceed to box 862 if the sequences for the expression vector 103 is found to within a margin of error.
• the margin of error may be five percent of the base pairs in the expression vector 103 sequence. In another embodiment, the margin of error is greater or smaller than five percent.
• the analysis system 207 searches the input sequence for similarities with the known sequences including a second reference sequence, illustratively adapter sequence 109.
• the method proceeds to box 866. If the adapter sequence 109 is found, in box 864 the method proceeds to box 880.
• the analysis system 207 must find the exact sequence of adapter sequence 109 to proceed to box 866. In another embodiment, the analysis system 207 may proceed to box 866 if the sequence for the adapter sequence 109 is found to within a margin of error. For example, the margin of error may be five percent of the base pairs in the adapter sequence 109. In another embodiment, the margin of error is greater or smaller than five percent.
• the method 550 proceeds to box 866.
• analysis system 207 attempts to identify the unknown sequence input in box 854.
• the known adapter is removed from the unknown sequence prior to further processing. In another embodiment, the known adapter is not removed from the unknown sequence prior to further processing. If the unknown sequence is identified, the method proceeds to box 870. If the unknown sequence is not identified, the method proceeds to box 878. The failure to identify the unknown sequence may indicate an error in the creation or the processing of the sequence.
• the input sequence is marked as failing processing. In an embodiment, the sequence is marked as red when the sequences are visualized.
• the input sequence is searched against the genome.
• the BLAST search algorithm is used to attempt to match the reduced input sequence to the genome.
• the method proceeds to box 874. If the reduced input sequence is not matched to any position in the genome, then the method proceeds to box 876.
• the input sequence matches against a portion of the genome.
• the analysis system 207 notes the location of the input sequence in the genome, and also notes the regions of interest in neighboring regions of the location. In an embodiment, the analysis system 207 notes regions of interest within 200 kilobase pairs of the location.
• the analysis system 207 notes regions of interest within a larger or smaller amount of base pairs.
• the user is able to specify the size of the neighboring region that the analysis system 207 notes around the location.
• the sequence is marked as green when the sequences are visualized.
• the input sequence is marked as failing to match against the genome.
• the reduced input sequence may have been damaged during sequencing, or may have been
• sequenced incorrectly In an embodiment, the sequence is marked as orange when the sequences are visualized.
• the method 850 proceeds to box 880.
• analysis system 207 attempts to identify the unknown sequence input in box 854. If the unknown sequence is identified in box 882, the method proceeds to box 886. If the unknown sequence is not identified, the method proceeds to box 884. The failure to identify the unknown sequence may indicate an error in the creation or the processing of the sequence.
• the input sequence is marked as failing processing. In an embodiment, the sequence is marked as red when the sequences are visualized.
• the input sequence is searched against the genome. In one embodiment, the BLAST search algorithm is used to attempt to match the reduced input sequence to the genome.
• the method proceeds to box 890. If the reduced input sequence is not matched to any position in the genome, then the method proceeds to box 892. [00100]
• box 890 the input sequence matches against a portion of the genome.
• the analysis system 207 notes the location of the input sequence in the genome, and also notes the regions of interest in neighboring regions of the location. In an embodiment, the analysis system 207 notes regions of interest within 200 kilobase pairs of the location. In other embodiments, the analysis system 207 notes regions of interest within a larger or smaller amount of base pairs. In an embodiment, the user is able to specify the size of the neighboring region that the analysis system 207 notes around the location. In an embodiment, the sequence is marked as green when the sequences are visualized. [00101] In box 892, the input sequence is marked as failing to match against the genome.
• the reduced input sequence may have been damaged during sequencing, or may have been sequenced incorrectly.
• the sequence is marked as orange when the sequences are visualized.
• Figure 5C is a flow chart showing another method of marking a transgene flanking sequence 507 according to the flow chart of Figure 5 A in which the known sequence for the primer 105, adapter 109, or both are provided in step 501.
• the analysis system 207 searches for the sequences identified as the primer 105 and the adapter 109 in the input sequence.
• the analysis system 207 searches for the adapter 109 and the primer 105 within the input sequence. If both the adapter 109 and the primer 105 sequences were provided in step 501 and are found within the input sequence, the method proceeds to box 559. If either the adapter 109 or the primer 105 sequences are not found within the input sequence, or if either the adapter 109 or the primer 105 sequences are not provided in step 501, the method proceeds to box 555. In an embodiment, the analysis system 207 must find the exact sequence of both the adapter 109 and the primer 105 to proceed to box 559.
• the analysis system 207 may proceed to box 559 if the sequences for the adapter 109 and the primer 105 are found to within a margin of error.
• the margin of error may be five percent of the base pairs in the adapter 109 or the primer 105 sequences. In another embodiment, the margin of error is greater or smaller than five percent. In another embodiment, the margin of error for the primer 105 and the margin of error for the adapter 109 are different.
• the known sequences for the adapter 109 and the primer 105 are removed from the input sequence, so that the input sequence is reduced to the sequence between the adapter 109 and the primer 105.
• the reduced input sequence is searched against the genome.
• the BLAST search algorithm is used to attempt to match the reduced input sequence to the genome.
• box 563 if the reduced input sequence is matched against the genome, the method proceeds to box 571. If the reduced input sequence is not matched to any position in the genome, then the method proceeds to box 565, and the input sequence is marked as failing to match against the genome.
• the reduced input sequence may have been damaged during sequencing, or may have been sequenced incorrectly, or the adapter 109 and the primer 105 may have abutted one another in the sequence, leaving no reduced input sequence. In an embodiment, the sequence is marked as orange when the sequences are visualized.
• box 571 the reduced input sequence matches against a portion of the genome.
• the analysis system 207 notes the location of the input sequence in the genome, and also notes the regions of interest in neighboring regions of the location. In an embodiment, the analysis system 207 notes regions of interest within 200 kilobase pairs of the location. In other embodiments, the analysis system 207 notes regions of interest within a larger or smaller amount of base pairs. In an embodiment, the user is able to specify the size of the neighboring region that the analysis system 207 notes around the location. In an embodiment, the sequence is marked as green when the sequences are visualized.
• the method proceeds from box 553 to box 555.
• the analysis system 207 determines if either of the adapter 109 or the primer 105 sequences are found in the input sequence. If either of the adapter 109 or the primer 105 sequences are found in the input sequence, the method proceeds to box 561. If both of the adapter 109 and the primer 105 sequences are not found in the input sequence, the method proceeds to box 557. [00108] In box 557, neither the adapter 109 nor the primer 105 were found within the input sequence. The lack of primer 105 and adapter 109 may indicate an error in the creation or the processing of the input sequence. The input sequence is marked as failing, and is not matched against the genome. In an embodiment, the sequence is marked as red when the sequences are visualized.
• either the adapter 109 or the primer 105 sequences are found within the input sequence.
• the adapter 109 or the primer 105 sequences are found within the input sequence to within a margin of error.
• the missing adapter 109 or primer 105 sequences indicates that the input sequence of the input sequence extends to either the 5 ' or the 3' end of the input sequence, and so the input sequence may not have captured the entire sequence of the input sequence.
• the known adapter 109 or the known primer 105 whichever is present in the input sequence, is removed from the input sequence so that the input sequence is reduced to the sequence between the adapter 109 and the primer 105.
• the reduced input sequence is searched against the genome, shown in box 567. In one embodiment, a BLAST search algorithm is used to attempt to match the reduced input sequence to the genome.
• box 567 if the reduced input sequence is matched against the genome, the method proceeds to box 573. If the reduced input sequence is not matched to any position in the genome, then the method proceeds to box 569, and the input sequence is marked as failing to match against the genome.
• the reduced input sequence may have been damaged during sequencing, or may have been sequenced incorrectly, or the adapter 109 and the primer 105 may have abutted one another in the sequence, leaving no reduced input sequence. In an embodiment, the sequence is marked as orange when the sequences are visualized.
• box 573 the reduced input sequence matches against a portion of the genome.
• the analysis system 207 notes the location of the input sequence in the genome, and also notes the regions of interest in neighboring regions of the location. In an embodiment, the analysis system 207 notes regions of interest within 200 kilobase pairs of the location. In other embodiments, the analysis system 207 notes regions of interest within a larger or smaller amount of base pairs. In an embodiment, the user is able to specify the size of the neighboring region that the analysis system 207 notes around the location. Regions of interest may include sequences encoding genes or other genomic information. Regions of interest may be received from a third party system, for example the system from which the analysis system 207 received the genome sequence information. In an embodiment, the sequence is marked as yellow when the sequences are visualized.
• Figure 7 shows a sample input screen for the analysis system 207.
• the user may select a series of input sequences in box 701.
• the input sequences may be in a standard form for providing sequence information, or may be a form that the analysis system 207 can parse and identify.
• the user may also select an organism's genome to map the input sequences against.
• the genome may be provided by the analysis system 207, so that the user identifies one or more genomes available to the analysis system 207, or the user may provide a path to an electronic file that contains sequence information for the organism's genome.
• the genome may be complete or partial.
• the user in box 705, selects one or more expression vectors 103 used in the experiment and which should be present in the input sequences.
• the user in boxes 707, 709, and 711, selects the vector sequences, the primer 105 sequences, and the adapter 109 sequences, respectively, that were used in the experiment and which should be present in the input sequences. The user then presses the "Submit” button to begin the data importation process and the analysis.
• Figure 8 shows an exemplary output of the analysis system 207 according to an embodiment of the present disclosure.
• the rows of the table labeled T indicate input sequences in which a chromosomal flanking sequence was identified correctly by the analysis system 207. These rows may be color coded, for example color coded green, for differentiation from the other rows.
• the rows of the table labeled '2' indicate input sequences in which a chromosomal flanking sequence was identified, but the analysis contains anomalies because all known sequences searched could not be identified so that, for example, the adapter 109 could not be located within the input sequence.
• rows may be coded as a different color than the rows of the table labeled ⁇ .
• the rows of the table labeled '3' indicate input sequences in which a chromosomal flanking sequence could not be identified. These rows are color coded as red.
• the Neighbors column indicate genes from a genomic sequence which proximal to the integration site.
• Figure 9A shows a summary display of the analysis system 207 which provides a graphical display of the integration site analysis for a particular input sequence from exemplary Soybean Event 416. At the top of the image, the coordinates of the input sequence are displayed. The remaining sequences that are shown within this summary display are annotated relative to these coordinates.
• the input reference sequence, in the exemplary screen, are oriented so that the primer 105 and transgene expression vector 103 appear on the left hand side of the screen, and the genomic flanking sequence and adapter 109 appear on the right hand side of the screen.
• the graphic display shows the input sequence for Event 416 (SEQ ID NO: l) (shown as Figure 9B) that has been annotated to identify the transgene expression vector 103 ("pDAB4468"; SEQ ID NO:2) (shown as Figure 9C), adapter 109 ("Soybe-”; SEQ ID NO:3 ) (shown as Figure 9D) and primer 105 ("soybean_primer "; SEQ ID NO:4) (shown as Figure 9E) sequences within it.
• the identified chromosomal flanking sequence is annotated as a solid line (SEQ ID NO:5)
• the analysis system 207 has aligned the chromosomal flanking sequence with the Glycine max genome.
• the chromosomal flanking sequence aligns to region 46003248, 46004030 of chromosome 4 with a sequence similarity score of 780; region 11825430, 11825559 of chromosome 6 with a sequence similarity score of 96; region 24517407, 24517435 of chromosome 15 with a sequence similarity score of 29; and region 37323425, 37323452 of chromosome 5 with a sequence similarity score of 28.
• the input sequence, the transgene expression vector 103, the adapter 109, and the primer 105 are graphically represented in the figure.
• Figure 10 shows the application of the analysis system 207 for use in Arabidopsis thaliana. Illustrated is the summary display of the analysis system 207 which provides an intuitive graphical display of the integration site analysis for an input sequence. At the top of the image, the coordinates of the input sequence are displayed. The remaining sequences that are shown within this summary display are annotated relative to these coordinates.
• the graphic display shows the input sequence for the event that has been annotated to identify the cloning vector ("pCR2.1-TOP") and adapter 109 (" ImAdp-Pri"). The identified chromosomal flanking sequence is annotated as a solid line.
• the analysis system 207 has aligned the chromosomal flanking sequence with the Arabidopsis genome sequence.
• FIG. 10 shows a transgene flanking sequence with a primer 105, but no right cloning vector 111.
• Figure 11 shows the application of the analysis system 207 for use in maize.
• the graphic display shows the input sequence for the event that has been annotated to identify the expression vector 103 ("pEPS1027").
• the identified chromosomal flanking sequence is annotated as a solid line.
• the analysis system 207 has aligned the chromosomal flanking sequence with the maize genome sequence. The chromosomal flanking sequence is aligned to a specific region of the Zea genomic sequence identifier 5337731, 5338124 and a sequence similarity score of 728 is reported.
• Figure 11 shows a transgene flanking sequence with an expression vector 103, but no right or left cloning vector s 101, 111.

Landscapes

• Physics & Mathematics (AREA)
• Life Sciences & Earth Sciences (AREA)
• Health & Medical Sciences (AREA)
• Engineering & Computer Science (AREA)
• Evolutionary Biology (AREA)
• General Health & Medical Sciences (AREA)
• Theoretical Computer Science (AREA)
• Bioinformatics & Cheminformatics (AREA)
• Bioinformatics & Computational Biology (AREA)
• Biotechnology (AREA)
• Spectroscopy & Molecular Physics (AREA)
• Biophysics (AREA)
• Medical Informatics (AREA)
• Bioethics (AREA)
• Databases & Information Systems (AREA)
• Chemical & Material Sciences (AREA)
• Analytical Chemistry (AREA)
• Proteomics, Peptides & Aminoacids (AREA)
• Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
• Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Priority Applications (10)

Applications Claiming Priority (4)

Publications (1)

Family

ID=48946332

Family Applications (1)

Country Status (14)

Cited By (15)

Families Citing this family (11)

Citations (3)

Family Cites Families (5)

• 2013
  • 2013-02-07 AR ARP130100389A patent/AR089934A1/es not_active Application Discontinuation
  • 2013-02-07 CA CA2863524A patent/CA2863524A1/en not_active Abandoned
  • 2013-02-07 AU AU2013217079A patent/AU2013217079B2/en not_active Ceased
  • 2013-02-07 WO PCT/US2013/025087 patent/WO2013119770A1/en active Application Filing
  • 2013-02-07 TW TW102104862A patent/TWI596493B/zh not_active IP Right Cessation
  • 2013-02-07 KR KR1020147021853A patent/KR20140119723A/ko not_active Application Discontinuation
  • 2013-02-07 EP EP13746881.5A patent/EP2812831A4/en not_active Withdrawn
  • 2013-02-07 BR BR112014019047A patent/BR112014019047A2/pt not_active Application Discontinuation
  • 2013-02-07 IN IN5963DEN2014 patent/IN2014DN05963A/en unknown
  • 2013-02-07 JP JP2014556652A patent/JP6314091B2/ja active Active
  • 2013-02-07 CN CN201380008411.9A patent/CN104272311B/zh not_active Expired - Fee Related
  • 2013-02-07 US US13/761,711 patent/US20130211729A1/en not_active Abandoned
• 2014
  • 2014-07-27 IL IL233819A patent/IL233819A0/en unknown
• 2015
  • 2015-02-09 HK HK15101413.0A patent/HK1201951A1/xx not_active IP Right Cessation

Patent Citations (3)

Non-Patent Citations (1)

Also Published As

Similar Documents

Legal Events