Abstract: Systems, methods, software and computer-usable media for annotating biomolecule-related data are disclosed. In certain exemplified embodiments, the biomolecules can be nucleic acids and the data can be sequence-related data. In various embodiments, systems can include one or more public or private biological attributes (e.g., annotation information databases, data storage devices and systems, etc.) sources, one or more genomic features data sources (e.g., genomic variant tools, genomic variant databases, genomic variant data storage devices and systems, etc.), a computing device (e.g., workstation, server, personal computer, mobile device, etc.) hosting an annotations module and/or a reporting module, and a client terminal.

Abstract: A targeted panel with low sample input requirements from a tumor only sample may be processed to estimate mutation load in a tumor sample. The method may include detecting variants in nucleic acid sequence reads corresponding to targeted locations in the tumor sample genome; annotating detected variants with an annotation information from a population database; filtering the detected variants, wherein the filtering rule set retains the somatic variants and removes germ-line variants; counting the identified somatic variants to give a number of somatic variants; determining a number of bases in covered regions of the targeted locations in the tumor sample genome; and calculating a number of somatic variants per megabase, provides an estimate of the mutation load per megabase in the tumor sample genome.

Abstract: A method for detecting a gene fusion includes amplifying a nucleic acid sample in the presence of primer pool to produce a plurality of amplicons. The primer pool includes primers targeting a plurality of exon-exon junctions of a driver gene. The amplicons correspond to the exon-exon junctions. The amplicons are sequenced and aligned to a reference sequence. The number of reads corresponding to each amplicon is normalized to give a normalized read count. A baseline correction is applied to the normalized read counts for the amplicons to form corrected read counts. A binary segmentation score is calculated for each corrected read count. A predicted breakpoint for the gene fusion is determined based on the amplicon index corresponding to the maximum absolute binary segmentation score. Gene fusion events may be detected in a partner agnostic manner, i.e. without prior knowledge of the specific fusion partner genes or specific breakpoint information.

Abstract: Systems and method for identifying gene fusions can obtain sequencing information for a plurality of amplicons from a nucleic acid sample. The sequencing information can include a plurality of reads that are initially partially mapped to a reference sequence. Fragments may be generated by splitting the partially mapped reads into mapped and unmapped fragments, and the fragments may be remapped to the reference sequence. Gene fusions can be identified based on reads where the first fragment maps to a first gene and the second fragment maps to a second gene.

Abstract: Systems and method for determining variants can receive mapped reads, align flow space information to a flow space representation of a corresponding portion of the reference. Reads spanning a position with a potential variant can be evaluated in a context specific manner. A list of probable variants can be provided.

Abstract: The method includes compressing numbers of reads data for targeted genes of a gene expression assay performed on a test sample. The targeted genes are organized into categories. Each category represents a functional context associated with the targeted genes in that category. The numbers of reads corresponding to targeted genes each category is compressed to form a compressed value for the category. The compressed value is compared to a baseline value for the category to determine an enrichment or a loss of a signature corresponding to the functional context of the category. The method may include analyzing information from multiple assays performed on the test sample, assigning a score value to each assay result and predicting a response to immune-oncology treatment based on the assigned scores.

Abstract: A method for detecting a gene fusion includes amplifying a nucleic acid sample in the presence of primer pool to produce a plurality of amplicons. The primer pool includes primers targeting a plurality of exon-exon junctions of a driver gene. The amplicons correspond to the exon-exon junctions. The amplicons are sequenced and aligned to a reference sequence. The number of reads corresponding to each amplicon is normalized to give a normalized read count. A baseline correction is applied to the normalized read counts for the amplicons to form corrected read counts. A binary segmentation score is calculated for each corrected read count. A predicted breakpoint for the gene fusion is determined based on the amplicon index corresponding to the maximum absolute binary segmentation score. Gene fusion events may be detected in a partner agnostic manner, i.e. without prior knowledge of the specific fusion partner genes or specific breakpoint information.

Abstract: The method includes compressing numbers of reads data for targeted genes of a gene expression assay performed on a test sample. The targeted genes are organized into categories. Each category represents a functional context associated with the targeted genes in that category. The numbers of reads corresponding to targeted genes each category is compressed to form a compressed value for the category. The compressed value is compared to a baseline value for the category to determine an enrichment or a loss of a signature corresponding to the functional context of the category. The method may include analyzing information from multiple assays performed on the test sample, assigning a score value to each assay result and predicting a response to immune-oncology treatment based on the assigned scores.

Abstract: A method for detecting large rearrangements in BRCA1 and BRCA2 genes includes amplifying a nucleic acid sample in the presence of a primer pool to produce amplicons, where the primer pool includes target specific primers targeting regions of exons of the BRCA1 and BRCA2 genes. The method further includes sequencing the amplicons to generate a plurality of reads, mapping the reads to a reference sequence, determining a number of reads per amplicon for the amplicons associated with the exons of the BRCA and the BRCA2 genes, determining exon copy numbers for the exons of the BRCA1 and BRCA2 genes based on the number of reads per amplicon, detecting an exon deletion or duplication based on the exon copy numbers, and detecting a whole gene deletion of the BRCA1 or BRCA2 gene based on the number of reads per amplicon associated with the exons of the BRCA1 and BRCA2 genes.

Abstract: Systems and method for determining variants can receive mapped reads, and call variants. In embodiments, flow space information for the reads can be aligned to a flow space representation of a corresponding portion of the reference. Reads spanning a position with a potential variant can be grouped and a score can be calculated for the variant. Based on the scores, a list of probable variants can be provided. In various embodiments, low frequency variants can be identified where multiple potential variants are present at a position.

Abstract: A computer-implemented method for classifying alignments of paired nucleic acid sequence reads is disclosed. A plurality of paired nucleic acid sequence reads is received, wherein each read is comprised of a first tag and a second tag separated by an insert region. Potential alignments for the first and second tags of each read to a reference sequence is determined, wherein the potential alignments satisfies a minimum threshold mismatch constraint. Potential paired alignments of the first and second tags of each read are identified, wherein a distance between the first and second tags of each potential paired alignment is within an estimated insert size range. An alignment score is calculated for each potential paired alignment based on a distance between the first and second tags and a total number of mismatches for each tag.

Abstract: Systems and method for identifying variants associated with a genetic disease can include obtaining sequencing reads for a plurality of individuals for a list of variant positions. The reads can be compared to identify variants that are found in affected individuals and absent in non-affected individuals. Such variants can include loss of heterozygosity, trans-phased compound heterozygotes, increased frequency mitochondrial variants, homozygous recessive variants, de novo variants, sex-linked variants, and combinations thereof.

Abstract: Systems and methods for analyzing overlapping sequence information can obtain first and second overlapping sequence information for a polynucleotide, align the first and second sequence information, determine a degree of agreement between the first and second sequence information for a location along the polynucleotide, and determine a base call and a quality value for the location.

Abstract: The method includes compressing numbers of reads data for targeted genes of a gene expression assay performed on a test sample. The targeted genes are organized into categories. Each category represents a functional context associated with the targeted genes in that category. The numbers of reads corresponding to targeted genes each category is compressed to form a compressed value for the category. The compressed value is compared to a baseline value for the category to determine an enrichment or a loss of a signature corresponding to the functional context of the category. The method may include analyzing information from multiple assays performed on the test sample, assigning a score value to each assay result and predicting a response to immune-oncology treatment based on the assigned scores.

Abstract: Systems and method for identifying variants associated with a genetic disease can include obtaining calls for a plurality of individuals for a list of variant positions. The calls can be compared to identify variants that are found in affected individuals and absent in non-affected individuals. Such variants can include loss of heterozygosity, trans-phased compound heterozygotes, increased frequency mitochondrial variants, homozygous recessive variants, de novo variants, sex-linked variants, and combinations thereof.

Abstract: The method includes compressing numbers of reads data for targeted genes of a gene expression assay performed on a test sample. The targeted genes are organized into categories. Each category represents a functional context associated with the targeted genes in that category. The numbers of reads corresponding to targeted genes each category is compressed to form a compressed value for the category. The compressed value is compared to a baseline value for the category to determine an enrichment or a loss of a signature corresponding to the functional context of the category. The method may include analyzing information from multiple assays performed on the test sample, assigning a score value to each assay result and predicting a response to immune-oncology treatment based on the assigned scores.

Abstract: Systems and method for identifying somatic mutations can receive first and second sequence information, determine if a variant present in the first sequencing information is also present in the second sequence information, and identify variants present in the first sequence information are somatic mutations when the variant is either not present in the second sequence information or the presence of the variant in the second sequence information is likely due to a sequencing error.

Abstract: Systems and methods for analyzing overlapping sequence information can obtain first and second overlapping sequence information for a polynucleotide, align the first and second sequence information, determine a degree of agreement between the first and second sequence information for a location along the polynucleotide, and determine a base call and a quality value for the location.

Abstract: Systems and method for identifying somatic mutations can receive first and second sequence information, determine if a variant present in the first sequencing information is also present in the second sequence information, and identify variants present in the first sequence information are somatic mutations when the variant is either not present in the second sequence information or the presence of the variant in the second sequence information is likely due to a sequencing error.

Abstract: Identification of exon junctions includes obtaining a first read sequence based on a detected plurality of signals of a first sequence. A list of exon prefix and suffix sequences are generated by identifying exons of the human genome with a prefix sequence mapping to a suffix sequence of the first read sequence and by identifying exons with a suffix sequence mapping to a prefix sequence of the first read sequence. A pair of exon sequences is selected, with a first exon sequence being one of the exon suffix sequences and a second exon sequence being one of the exon prefix sequences. Summing a number of sequence elements of the first exon sequence that overlap the prefix of the first read sequence, a number of sequence elements of the second exon sequence that overlap the suffix of the first read sequence, and a constant is used to identify a fusion junction.