GENOMIC SERVICES PLATFORM SUPPORTING MULTIPLE APPLICATION PROVIDERS
Systems, platforms, methods and media for providing genomic services axe disclosed. In one example, a system for providing genomic services comprises genomic sequencing equipment configured to generate sequence reads based upon a biological sample obtained from a user, store the sequence reads in a FASTQ storage file, and communicate the FASTQ file electronically to a recipient. The system also includes a genomic services platform which includes a network interface through which the sequence reads are received, and a bioinformatics processing pipeline. The bioinformatics processing pipeline includes a read alignment module configured to generate observed sequence data, and a variant calling module is operative to identify observed variants in the observed sequence data. A variant storage module is disposed to receive a query from network infrastructure of a partner application provider and to provide selected ones of refined variants in response to a query.
1 . A system for providing genomic services, the system comprising:
genomic sequencing equipment wherein the genomic sequencing equipment is configured to:
generate sequence reads based upon a biological sample obtained from a user,
store the sequence reads in a FASTQ storage file; and
communicate the FASTQ file electronically to a recipient;
a genomic services platform, the platform including:
a network interlace through which are received the sequence reads,
a bioinformatics processing pipeline, the bioinformatics processing pipeline including:
a read alignment module configured to:
receive the FASTQ file from the genomic sequencing equipment and use data contained therein to generate observed sequence data by aligning the sequence reads relative to a reference sequence; and
store the alignments in a file in a compressed format,
a variant calling module operative to:
identify observed variants in the observed sequence data, the variants including one or more of a single nucleotide polymorphism, an insertion, and a deletion; and
store the observed variants in a variant calling file (VCF) format, or a genomic variant call file (GVCF) format; and a variant refinement module for producing a set of refined variants associated with the user,
genomic data storage containing the set of refined variants; and
a variant storage module disposed to receive a query from network infrastructure of a partner application provider and to provide selected ones of the refined variants in response to the query.
2 . The system of claim 1 , wherein read alignment module utilizes a mapping algorithm to compare the sequence of a given read to that of the reference sequence to locate a potentially unique location in the reference sequence that matches the read.
3 . The system of claim 1 , wherein the read alignment module is further configured to store the alignments in a compressed format in a Binary Alignment Map (BAM) file, and index the BAM file relative to the reference sequence.
4 . The system of claim wherein the variant calling module is further operative to process the BAM file to identify the existence of the observed variants in the observed sequence data.
5 . The system of claim 1 , wherein the variant calling module is further operative to produce at least two variant data files, the at least two variant data files including:
a first genomic variant data file (GVCF) providing genomic information about all sites in the sequence reads obtained from the user genome, the sites including both sites with variants and reference call sites without variants, and
a second variant data file (VCF) providing genomic information about sites in the sequence reads obtained from the user genome, the genomic information including posterior genotype likelihoods for sites with variants, but not including information about reference call sites without variants.
6 . The system of claim 5 , wherein the variant refinement module is operative to merge the first GVCF and second VCF files to produce a single variant call file, the single variant call file including reference calls from the first GVCF file and variants calls with posterior genotype likelihoods from the second VCF file.
7 . The system of claim 1 wherein the genomic services platform further includes file storage including one or more Browser Extensible Data (BED) files defining one or more genomic windows associated with the partner application provider or a partner application, and wherein the selected ones of the refined variants are associated with reference locations within one of the one or more genomic windows.
8 . A method for providing genomic services, the method comprising:
generating or receiving sequence reads based upon a biological sample obtained from a user;
storing the sequence reads in a FASTQ storage tile and communicating the FASTQ file electronically to a recipient,
using data stored in the FASTQ file to generate observed sequence data by aligning the sequence reads relative to a reference sequence;
storing the alignments in a rile in a compressed format;
identifying observed variants in the observed sequence data, the variants including one or more of a single nucleotide polymorphism, an insertion, and a deletion;
storing the observed variants in a variant calling file (VCF) format, or a genomic variant call tile (GVCF) format;
producing and storing a set of refined variants associated with the user; and receiving a query from network infrastructure of a partner application provider and providing selected ones of the refined variants in response to the query.
9 . The method of claim 8 , further comprising utilizing a mapping algorithm to compare the sequence of a given read to that of the reference sequence to locate a potentially unique location in the reference sequence that matches the read.
10 . The method of claim 8 , further comprising storing the alignments in a compressed format in a Binary Alignment Map (BAM) file, and indexing the BAM file relative to the reference sequence.
11 The method of claim 10 , further comprising processing the BAM file to identify the existence of the observed variants in the observed sequence data.
12 . The method of claim 8 , further comprising producing at least two variant data files, the at least two variant data files including:
a first genomic variant data file (GVCF) providing genomic information about all sites in the sequence reads obtained from the user genome, the sites including both sites with variants and reference call sites without variants; and
a second variant data file (VCF) providing genomic information about sites in the sequence reads obtained from the user genome, the genomic information including posterior genotype likelihoods for sites with variants, but not including information about reference call sites without variants.
13 . The method, of claim 12 , further comprising merging the first GVCF and second VCF files to produce a single variant call file, the single variant call file including reference calls from the first GVCF file and variants calls with posterior genotype likelihoods from the second VCF file.
14 . The method of claim 7 , further comprising using or storing one or more Browser Extensible Data (BED) files defining one or more genomic windows associated with the partner application provider or a partner application, and wherein the selected ones of the refined variants are associated with reference locations within one of the one or more genomic windows.
15 . A non-transitory machine-readable medium comprising instructions for implementing a method for providing genomic services, the instructions, when read by a machine, causing the machine to perform operations comprising, at least.
generating or receiving sequence reads based upon a biological sample obtained from a user;
storing the sequence reads in a FASTQ storage file and communicating the FASTQ file electronically to a recipient;
using data stored in the FASTQ file to generate observed sequence data by aligning the sequence reads relative to a reference sequence;
storing the alignments in a file in a compressed format;
identifying observed variants in the observed sequence data, the variants including one or more of a single nucleotide polymorphism, an insertion, and a deletion;
storing the observed variants in a variant calling file (VCF) format, or a genomic variant call file (GVCF) format;
producing and storing a set of refined variants associated with the user; and
receiving a query from network infrastructure of a partner application provider and providing selected ones of the refined variants in response to the query
16 . The medium of claim 15 , wherein the operations further comprise utilizing a mapping algorithm to compare the sequence of a given read to that of the reference sequence to locate a potentially unique location in the reference sequence that matches the read.
17 . The medium of claim 15 , wherein the operations further comprise storing the alignments in a compressed format in a Binary Alignment Map (BAM) file, and indexing the BAM file relative to the reference sequence
18 . The medium of claim 17 , wherein the operations further comprise processing the BAM file to identify the existence of the observed variants in the observed sequence data.
19 . The medium of claim 15 , wherein the operations further comprise producing at least two variant data files, the two variant data files including:
a first genomic variant data file (GVCF) providing genomic information about all sites in the sequence reads obtained from the user genome, the sites including both sites with variants and reference call sites without variants, and
a second variant data file (VCF) providing genomic information about sites in the sequence reads obtained from the user genome, the genomic information including posterior genotype likelihoods for sites with variants, but not including information about reference call sites without variants.
20 . The medium of claim 19 , wherein the operations further comprise merging the first GVCF and second VCF files to produce a single variant call file, the single variant call file including reference calls from the first GVCF file and variants calls with posterior genotype likelihoods from the second VCF file.
21 . The medium of claim 15 , wherein the operations further comprise using or storing one or more Browser Extensible Data (BED) files defining one or more genomic windows associated with the partner application provider or a partner application, and wherein the selected ones of the refined variants are associated with reference locations within one of the one or more genomic windows.