Patent Grant
8700337
Technological advance has greatly reduced the cost of genetic information to the point where it is possible to contemplate individual genome sequencing. Coupled with the enormous increase in available genetic data, there has been an ongoing effort into the assessment of associations between human genetic variation, physiology, and disease risk. As an increasing number of individual genome sequences are made available, there can be an evaluation of the biological, population, medical, and physical data for association statistics. This pool of combined data may provide for further information on the kinds and levels of variation that exist generally throughout and between individual genomes.
Individual genome information has the potential for personal health benefits by providing the information content of a large number of individual genetic tests, which may predict risk for serious disease. The availability of such knowledge can be utilized in further testing, medical intervention, and long term surveillance for signs of disease development.
However, the explanatory power and path to clinical translation of risk estimates for common variants reported in genome-wide association studies remain unclear. Much of the reason lies in the presence of rare and structural genetic variation. With the availability of rapid an inexpensive sequencing of complete genomes, comprehensive genetic risk assessment and individualization of treatment might be possible. However, present analytical methods are insufficient to make genetic data accessible in a clinical context, and the clinical usefulness of these data for individual patients has not been formally assessed. The present invention addresses the integrated analysis of a complete human genome in a clinical context.
Methods and systems are provided for the computation and display of an individual's personalized health risk based on the individual genome sequence, known etiological interactions between diseases for which the individual is determined to have genetic risk factors, and environmental etiological factors associated with the same diseases that represent potentially modifiable disease risk modifiers. The analysis is performed with novel computational methods, and is optionally displayed with a visualization approach useful in the clinical interpretation of risk factors personalized for the individual.
In a first component of the invention, an etiology integration engine takes as input the genetic data from the genome of an individual, utilizing per-disease combined genetic risks. Using a knowledge base system, the etiology integration engine (i) makes etiological connections between diseases for which the individual has been determined to have a genetic risk and (ii) makes etiological connections between these diseases and environmental factors that are disease risk modifiers.
In a second component, using network visualization methods, relationships between components, e.g. the disease and environmental conditions involved in integrated genomic and etiological risks for an individual, are graphically displayed, to aid in the identification of such relationships.
The risk assessment of the invention is useful in guiding treatment and preventive care for a patient. In some embodiments of the invention an individual is provided with an output that comprises a risk assessment diagram. In some embodiments an individual is further provided with a personalized guide for prevention and/or treatment of the risks identified by the methods of the invention.
In other embodiments, kits are provided for assessing integrated genetic and etiological risk in a human individual, the kit comprising reagents for data input and analysis for the methods described herein.
The above summary is not intended to include all features and aspects of the present invention nor does it imply that the invention must include all features and aspects discussed in this summary.
All publications and patent applications mentioned in this specification are herein incorporated by reference to the same extent as if each individual publication or patent application was specifically and individually indicated to be incorporated by reference.
Genome Sequence Information.
As known in the art, the complete genome of an individual comprises the total chromosomal genetic sequence information. For the purposes of the present invention, the genome sequence information comprises all or a portion of the total genome sequence from the individual. In some embodiments the genome sequence information comprises the total non-redundant sequence of an individual genome. In other embodiments the genome sequence information comprises a dataset of single nucleotide polymorphisms (SNP) data from an individual obtained by sequencing or hybridization or a combination thereof, preferably by sequencing, e.g. at least about 105 single nucleotide polymorphisms, usually at least about 5×105 SNPs, more usually at least about 106 SNPs, and may be 2×106 SNPs or more. The genome sequence information may further comprise a dataset for copy number variations in the individual genome, e.g. at least about 100 copy number variations, at least about 200 copy number variations, at least about 500 copy number variations, or more.
Sequence analysis can also be used to detect specific polymorphisms in a nucleic acid. A test sample of DNA or RNA is obtained from the test individual. PCR or other appropriate methods can be used to amplify the gene or nucleic acid, and/or its flanking sequences, if desired. The sequence of a nucleic acid, or a fragment of the nucleic acid, or cDNA, or fragment of the cDNA, or mRNA, or fragment of the mRNA, is determined, using standard methods. The sequence of the nucleic acid, nucleic acid fragment, cDNA, cDNA fragment, mRNA, or mRNA fragment is compared with the known nucleic acid sequence of the gene or cDNA or mRNA, as appropriate.
Sequencing platforms that can be used in the present disclosure include but are not limited to: pyrosequencing, sequencing-by-synthesis, single-molecule sequencing, second-generation sequencing, nanopore sequencing, sequencing by ligation, or sequencing by hybridization. Preferred sequencing platforms are those commercially available from Illumina (RNA-Seq) and Helicos (Digital Gene Expression or “DGE”). “Next generation” sequencing methods include, but are not limited to those commercialized by: 1) 454/Roche Lifesciences including but not limited to the methods and apparatus described in Margulies et al., Nature (2005) 437:376-380 (2005); and U.S. Pat. Nos. 7,244,559; 7,335,762; 7,211,390; 7,244,567; 7,264,929; 7,323,305; 2) Helicos BioSciences Corporation (Cambridge, Mass.) as described in U.S. application Ser. No. 11/167,046, and U.S. Pat. Nos. 7,501,245; 7,491,498; 7,276,720; and in U.S. Patent Application Publication Nos. US20090061439; US20080087826; US20060286566; US20060024711; US20060024678; US20080213770; and US20080103058; 3) Applied Biosystems (e.g. SOLiD sequencing); 4) Dover Systems (e.g., Polonator G.007 sequencing); 5) IIlumina as described U.S. Pat. Nos. 5,750,341; 6,306,597; and 5,969,119; and 6) Pacific Biosciences as described in U.S. Pat. Nos. 7,462,452; 7,476,504; 7,405,281; 7,170,050; 7,462,468; 7,476,503; 7,315,019; 7,302,146; 7,313,308; and US Application Publication Nos. US20090029385; US20090068655; US20090024331; and US20080206764. All references are herein incorporated by reference. Such methods and apparatuses are provided here by way of example and are not intended to be limiting.
Allele-specific oligonucleotides can also be used to detect the presence of a polymorphism in a nucleic acid, through the use of dot-blot hybridization of amplified oligonucleotides with allele-specific oligonucleotide (ASO) probes (see, for example, Saiki, R. et al., Nature 324:163-166 (1986)). An “allele-specific oligonucleotide” (also referred to herein as an “allele-specific oligonucleotide probe”) is an oligonucleotide of approximately 10-50 base pairs that specifically hybridizes to a nucleic acid, and that contains a polymorphism. An allele-specific oligonucleotide probe that is specific for particular polymorphisms in a nucleic acid can be prepared, using standard methods. With the addition of such analogs as locked nucleic acids (LNAs), the size of primers and probes can be reduced to as few as 8 bases.
In another aspect, arrays of oligonucleotide probes that are complementary to target nucleic acid sequence segments from an individual can be used to identify the presence of polymorphic alleles in a nucleic acid. For example, in one aspect, an oligonucleotide array can be used. Oligonucleotide arrays typically comprise a plurality of different oligonucleotide probes that are coupled to a surface of a substrate in different known locations. These oligonucleotide arrays, also described as “Genechips™,” have been generally described in the art, for example, U.S. Pat. No. 5,143,854 and PCT patent publication Nos. WO 90/15070 and 92/10092. These arrays can generally be produced using mechanical synthesis methods or light directed synthesis methods that incorporate a combination of photolithographic methods and solid phase oligonucleotide synthesis methods. See Fodor et al., Science 251:767-777 (1991), Pirrung et al., U.S. Pat. No. 5,143,854 (see also PCT Application No. WO 90/15070) and Fodor et al., PCT Publication No. WO 92/10092 and U.S. Pat. No. 5,424,186, the entire teachings are incorporated by reference herein. Techniques for the synthesis of these arrays using mechanical synthesis methods are described in, e.g., U.S. Pat. No. 5,384,261; the entire teachings are incorporated by reference herein. In another example, linear arrays can be utilized.
Once an oligonucleotide array is prepared, a nucleic acid of interest is hybridized with the array and scanned for polymorphisms. Hybridization and scanning are generally carried out by methods described herein and also in, e.g., published PCT Application Nos. WO 92/10092 and WO 95/11995, and U.S. Pat. No. 5,424,186, the entire teachings are incorporated by reference herein. In brief, a target nucleic acid sequence that includes one or more previously identified polymorphic markers is amplified by well-known amplification techniques, e.g., PCR. Typically, this involves the use of primer sequences that are complementary to the two strands of the target sequence both upstream and downstream from the polymorphism. Asymmetric PCR techniques may also be used. Amplified target, generally incorporating a label, is then hybridized with the array under appropriate conditions. Upon completion of hybridization and washing of the array, the array is scanned to determine the position on the array to which the target sequence hybridizes. The hybridization data obtained from the scan is typically in the form of fluorescence intensities as a function of location on the array.
Polymorphism, as used herein refers to variants in the gene sequence. Such variants may include single nucleotide polymorphisms, splice variants, insertions, deletions and transpositions. The polymorphisms can be those variations (DNA sequence differences) that are generally found between individuals or different ethnic groups and geographic locations which, while having a different sequence, produce functionally equivalent gene products. Polymorphisms also encompass variations which can be classified as alleles and/or mutations which can produce gene products which may have an altered function, i.e. variants in the sequence which can lead to gene products that are not functionally equivalent. Polymorphisms also encompass variations which can be classified as alleles and/or mutations which either produce no gene product, an inactive gene product or increased gene product. Further, the term is also used interchangeably with allele as appropriate.
Where a polymorphic site is a single nucleotide in length, the site is referred to as a single nucleotide polymorphism (“SNP”). SNP nomenclature used herein refers to the official Reference SNP (rs) ID identification tag as assigned to each unique SNP by the National Center for Biotechnological Information (NCBI), although other references may find use.
Risk Factor.
As used herein, a risk factor is a genetic or environmental exposure or characteristic, which, on the basis of epidemiologic evidence, is known to be associated with a health-related condition considered important to prevent. The term “risk factor”, as described herein, means primarily increased susceptibility to a human disease, e.g. diabetes, cancer, cardiovascular disease, and the like. Thus, particular genetic variants or environmental exposures may be characteristic of increased susceptibility of disease, as characterized by a relative risk of greater than one.
Likelihood Ratio.
For genetic variants, e.g. SNP variants, copy number variants, etc., a likelihood ratio (LR) of disease risk may be calculated published information reporting association with disease, which association is optionally stratified by sex, ethnicity, age, etc. For each variant the LR may be calculated as the probability of the genotype appearing in a disease, or case population/probability of the genotype appearing in control population.
To calculate the likelihood ratio (LR), a genetic association database is utilized that provides case/control associations with diseases in a population of interest. For every disease SNP, we calculated the LR for each genotype using the following equation:
Odds Ratio.
The odds ratio is a measure of effect size, describing the strength of association or non-independence between two binary data values. It is used as a descriptive statistic, and plays an important role in logistic regression. Unlike other measures of association for paired binary data such as the relative risk, the odds ratio treats the two variables being compared symmetrically, and can be estimated using some types of non-random samples.
The odds ratio is the ratio of the odds of an event occurring in one group to the odds of it occurring in another group. An odds ratio of 1 indicates that the condition or event under study is equally likely to occur in both groups. An odds ratio greater than 1 indicates that the condition or event is more likely to occur in the first group. And an odds ratio less than 1 indicates that the condition or event is less likely to occur in the first group. The odds ratio must be greater than or equal to zero if it is defined. In clinical studies the parameter of greatest interest is often the relative risk rather than the odds ratio. The relative risk is best estimated using a population sample, but if the rare disease assumption holds, the odds ratio is a good approximation to the relative risk—the odds is p/(1−p), so when p moves towards zero, 1−p moves towards 1, meaning that the odds approaches the risk, and the odds ratio approaches the relative risk.
A feature of the invention is the generation of a database of risk associations for a variety of genetic variants and environmental factors. Such a database will typically comprise associations as described above, for a number of variants and factors, which may be selected and arranged according to various criteria. The databases may be provided in a variety of media to facilitate their use. “Media” refers to a manufacture that contains the datasets of the present invention. The datasets can be recorded on computer readable media, e.g. any medium that can be read and accessed directly by a computer. Such media include, but are not limited to: magnetic storage media, such as floppy discs, hard disc storage medium, and magnetic tape; optical storage media such as CD-ROM; electrical storage media such as RAM and ROM; and hybrids of these categories such as magnetic/optical storage media. One of skill in the art can readily appreciate how any of the presently known computer readable mediums can be used to create a manufacture comprising a recording of the present database information. “Recorded” refers to a process for storing information on computer readable medium, using any such methods as known in the art. Any convenient data storage structure may be chosen, based on the means used to access the stored information. A variety of data processor programs and formats can be used for storage, e.g. word processing text file, database format, etc.
As used herein, “a computer-based system” refers to the hardware means, software means, and data storage means used to analyze the information of the present invention. The minimum hardware of the computer-based systems of the present invention comprises a central processing unit (CPU), input means, output means, and data storage means. A skilled artisan can readily appreciate that any one of the currently available computer-based system are suitable for use in the present invention. The data storage means may comprise any manufacture comprising a recording of the present information as described above, or a memory access means that can access such a manufacture.
As used herein, the term “nucleic acid probe” refers to a molecule capable of sequence specific hybridization to a nucleic acid, and includes analogs of nucleic acids, as are known in the art, e.g. DNA, RNA, peptide nucleic acids, and the like, and may be double-stranded or single-stranded. Also included are synthetic molecules that mimic nucleic acid molecules in the ability to bind to a designated sequence via hydrogen bonding and other chemical interactions. Such molecules include, for example, those in which peptide linkages substitute for phosphate linkages in the backbone of the molecule.
“Specific hybridization,” as used herein, refers to the ability of a first nucleic acid to hybridize to a second nucleic acid in a manner such that the first nucleic acid does not hybridize to any nucleic acid other than to the second nucleic acid. “Stringency conditions” for hybridization is a term of art which refers to the incubation and wash conditions, e.g., conditions of temperature and buffer concentration, which permit hybridization of a particular nucleic acid to a second nucleic acid; the first nucleic acid may be perfectly (i.e., 100%) complementary to the second, or the first and second may share some degree of complementarity which is less than perfect (e.g., 70%, 75%, 85%, 90%, 95%). The percent homology or identity of two nucleotide or amino acid sequences can be determined by aligning the sequences for optimal comparison purposes (e.g., gaps can be introduced in the sequence of a first sequence for optimal alignment). The nucleotides or amino acids at corresponding positions are then compared, and the percent identity between the two sequences is a function of the number of identical positions shared by the sequences (i.e., % identity=# of identical positions/total # of positions×100). When a position in one sequence is occupied by the same nucleotide or amino acid residue as the corresponding position in the other sequence, then the molecules are homologous at that position. As used herein, nucleic acid or amino acid “homology” is equivalent to nucleic acid or amino acid “identity”. A preferred, non-limiting example of such a mathematical algorithm is described in Karlin et al., Proc. Natl. Acad. Sci. USA 90:5873-5877 (1993). Such an algorithm is incorporated into the NBLAST and XBLAST programs (version 2.0) as described in Altschul et al., Nucleic Acids Res. 25:389-3402 (1997). When utilizing BLAST and Gapped BLAST programs, the default parameters of the respective programs (e.g., NBLAST) can be used. In one aspect, parameters for sequence comparison can be set at score=100, wordlength=12, or can be varied (e.g., W=5 or W=20).
A “marker”, as described herein, refers to a genomic sequence characteristic of a particular allele at a polymorphic site.
A “haplotype,” as described herein, refers to a segment of a genomic DNA strand that is characterized by a specific combination of genetic markers (“alleles”) arranged along the segment. In a certain embodiment, the haplotype can comprise one or more alleles, two or more alleles, three or more alleles, four or more alleles, or five or more alleles.
“Comparable cell” shall mean a cell whose type is identical to that of another cell to which it is compared. Examples of comparable cells are cells from the same cell line.
“Inhibiting” the onset of a disorder shall mean either lessening the likelihood of the disorder's onset, or preventing the onset of the disorder entirely. In the preferred embodiment, inhibiting the onset of a disorder means preventing its onset entirely.
“Treating” a disorder shall mean slowing, stopping or reversing the disorder's progression. In the preferred embodiment, treating a disorder means reversing the disorder's progression, ideally to the point of eliminating the disorder itself. As used herein, ameliorating a disorder and treating a disorder are equivalent.
“Inhibiting” the expression of a gene in a cell shall mean either lessening the degree to which the gene is expressed, or preventing such expression entirely. “Specifically inhibit” the expression of a protein shall mean to inhibit that protein's expression (a) more than the expression of any other protein, or (b) more than the expression of all but 10 or fewer other proteins.
“Subject” or “patient” shall mean any animal, such as a human, non-human primate, mouse, rat, guinea pig or rabbit.
“Suitable conditions” shall have a meaning dependent on the context in which this term is used. That is, when used in connection with an antibody, the term shall mean conditions that permit an antibody to bind to its corresponding antigen. When this term is used in connection with nucleic acid hybridization, the term shall mean conditions that permit a nucleic acid of at least 15 nucleotides in length to hybridize to a nucleic acid having a sequence complementary thereto. When used in connection with contacting an agent to a cell, this term shall mean conditions that permit an agent capable of doing so to enter a cell and perform its intended function. In one embodiment, the term “suitable conditions” as used herein means physiological conditions.
Unless otherwise apparent from the context, all elements, steps or features of the invention can be used in any combination with other elements, steps or features.
General methods in molecular and cellular biochemistry can be found in such standard textbooks as Molecular Cloning: A Laboratory Manual, 3rd Ed. (Sambrook et al., Harbor Laboratory Press 2001); Short Protocols in Molecular Biology, 4th Ed. (Ausubel et al. eds., John Wiley & Sons 1999); Protein Methods (Bollag et al., John Wiley & Sons 1996); Nonviral Vectors for Gene Therapy (Wagner et al. eds., Academic Press 1999); Viral Vectors (Kaplift & Loewy eds., Academic Press 1995); Immunology Methods Manual (I. Lefkovits ed., Academic Press 1997); and Cell and Tissue Culture: Laboratory Procedures in Biotechnology (Doyle & Griffiths, John Wiley & Sons 1998). Reagents, cloning vectors, and kits for genetic manipulation referred to in this disclosure are available from commercial vendors such as BioRad, Stratagene, Invitrogen, Sigma-Aldrich, and ClonTech.
The present invention has been described in terms of particular embodiments found or proposed by the present inventor to comprise preferred modes for the practice of the invention. It will be appreciated by those of skill in the art that, in light of the present disclosure, numerous modifications and changes can be made in the particular embodiments exemplified without departing from the intended scope of the invention. For example, due to codon redundancy, changes can be made in the underlying DNA sequence without affecting the protein sequence. Moreover, due to biological functional equivalency considerations, changes can be made in protein structure without affecting the biological action in kind or amount. All such modifications are intended to be included within the scope of the appended claims.
The following examples are put forth so as to provide those of ordinary skill in the art with a complete disclosure and description of how to make and use the present invention, and are not intended to limit the scope of what the inventors regard as their invention nor are they intended to represent that the experiments below are all or the only experiments performed. Efforts have been made to ensure accuracy with respect to numbers used (e.g. amounts, temperature, etc.) but some experimental errors and deviations should be accounted for. Unless indicated otherwise, parts are parts by weight, molecular weight is weight average molecular weight, temperature is in degrees Centigrade, and pressure is at or near atmospheric.
All publications and patent applications cited in this specification are herein incorporated by reference as if each individual publication or patent application were specifically and individually indicated to be incorporated by reference.
The present application may make reference to information provided in Ashley et al. (2010) Lancet 375:1525-35, including supplemental materials provided therein, which is herein specifically incorporated by reference in its entirety.
Methods of integrating genetic and environmental risk factors to provide a personalized overall risk assessment for an individual are provided herein. By evaluating an individual genotype for the presence of polymorphisms in relevant genes, integrated with environmental factors influencing the individual, the susceptibility to development of disease can be predicted. The knowledge about integrated risks allows for the ability to maintain specific testing and thus to diagnose the disease at an early stage, to provide information to the clinician about prognosis for disease in order to be able to apply the most appropriate treatment, to provide assessment for lifestyle changes commensurate with risk factors, and the like. In some embodiments of the invention an individual is provided with an output that comprises a risk assessment diagram. In some embodiments an individual is further provided with a personalized guide for prevention and/or treatment of the risks identified by the methods of the invention.
In some embodiments, the methods of the invention comprise inputting genome sequence information from an individual into an etiology integration engine to make connections between diseases for which the individual has been determined to have a genetic risk, and between these diseases and environmental factors that are disease risk modifiers. The genome sequence information comprises all or a portion of the total genome sequence from the individual. The genome sequence information may comprise the total sequence of an individual genome. Alternatively the input genome sequence information may comprise a set of SNP data, e.g. at least about 105 SNPs, usually at least about 5×105 SNPs, and may be 106 SNPs, 2×106 SNPs or more.
The genome sequence information is initially analyzed for variants in coding and non-coding regions, utilizing public databases of human genetic sequence information. Analysis is performed on four areas: (i) variants associated with genes for mendelian disease; (ii) novel mutations; (iii) variants known to modulate response to pharmacotherapy; and (iv) single nucleotide polymorphisms previously associated with complex disease.
The analysis includes a comparison with variants having validated reference SNPs, and an analysis for novel variants. Public databases of sequence information include, without limitation, dbSNP from the National Institutes of Health, the SNP consortium, European SNP database, and the human genome variation database. A comparison is made with the validated SNPs and the match data for the genome being analyzed is input to the etiology integration engine.
In some embodiments of the invention a high-quality disease-associated SNP database is utilized, where the database is built from the publicly available SNP information at dbSNP from those SNPs relevant to human disease. The database entries for each SNP may include, without limitation, a disease name, specific phenotype, study population, case and control population, genotyping technology, major/minor/risk alleles, likelihood ratio, likelihood ratio, odds ratio, 95% confidence interval of the odds ratio, published p-value, and genetic model for each included, statistically significant genotype comparison. Preferred information includes a disease name and odds ratio for the SNP. To enable the integration of multiple studies on similar diseases and phenotypes, the disease/phenotype names in our association database are mapped to the Unified Medical Language System (UMLS) Concept Unique Identifiers (CUIs). The SNP data may further be subjected to an algorithm to correctly identify the strand direction and suitably annotated.
The analysis of novel variants may utilize one or more independent analysis tools, usually two or more, three, four or more analysis tools. Tools include the (a) “Sorting Intolerant from Tolerant” (SIFT) algorithm, which predicts the effects of non-synonymous polymorphisms on protein function based on homology, conservation, and physical properties of amino acid substitutions; (b) the “Polymorphism Phenotyping” (PolyPhen) tool, which predicts the impact of amino acid changes on protein function using an algorithm that incorporates information on site of substitution, i.e., whether an amino acid change occurs in one of several sites of functional importance such as binding sites or trans-membrane regions, sequence alignment, and known protein structural changes; (c) query for rare coding variants using the Universal Protein Resource (UniProt) database consisting of annotated protein sequence variation data with experimentally proven or computer-predicted data to support phenotypic association; and analysis for coding region variants that produce premature stop codons or read-throughs in existing stop codons, e.g. using the PolyDoms database, which incorporates SNP-phenotype associations from several different sources, including SIFT, PolyPhen, Uniprot, Database of Secondary Structure in Proteins, Protein Databank, and Database of Secondary Structure in Proteins.
The genome sequence information is also analyzed for variants in on-coding regions, e.g. the SNP databases referenced above, and other public databases, which include without limitation, Online Mendelian Inheritance of Man; Human Gene Mutation Database; Pubmed, and the like.
The analysis of genetic information thus obtained provides an initial estimate of disease association between the genome information of the individual and known disease correlations. A likelihood ratio (LR) of disease risk is calculated for each SNP having a disease association. The pre-test probability for disease is calculated from the LR using criteria specific for the individual, which criteria may include, without limitation, age, sex, weight, ethnicity, and the like.
A post-test calculation of developing disease is then made. The post-test probability is calculated as follows. For SNPs with multiple LR from multiple studies, the mean LR is calculated, and weighted by the square root of sample sizes. The human genome is partitioned into Haplotype blocks, and for each haplotype block, the highest LR SNP is used. LR from all SNPs are multiplied to report the cumulative LR for the individual, and the cumulative LR calculated for all relevant diseases. The pre-test probability is translated into pre-test odds, multiplied by the cumulative LR to get post-test odds, and then converted to the post-test probability using the following equations:
pre-test odds=pre-test probability/1−pre-test probability
post-test odds=pre-test odds×LR
pest-test probability=post-test odds/1+post-test odds
As a second component of the risk integration analysis, gene-environment interaction and conditionally dependent risk are analyzed. A database of established links between diseases and known etiological factors is obtained by mapping Medical Subject Heading (MeSH) annotations in the MEDLINE database to the Unified Medical Language System (UMLS) semantic network, keeping only those publications with annotations representing diseases and environmental factors. A compilation of drug-related genotype-phenotype associations is then drawn from a public database, utilizing, for example, the Pharmacogenomics Knowledgebase.
An output of the combined risk factor thus obtained may be presented to the individual, including the post-test calculations of disease and environmental interactions. In a preferred embodiment, the information thus obtained is presented in a conditional dependency diagram for diseases represented in the patient's genetic-risk profile. Such a diagram may represent diseases having a significant post-test risk probability, e.g. 5%, 10%, 15%, etc., with connections between diseases if one disease predisposes a patient to the other. Environmental factors, e.g. those that are potentially modifiable may be presented, including showing an association between an environmental factor and a disease if there is a reported association.
In addition to distance and association visualization, the display of information may include other classification schemes to aid in analysis. Each point, which represents a disease in the analysis matrix, may be arbitrarily assigned features, such as color, size, shape, etc. where the assignment provides information about the condition. For example, the size of the point may represent the risk of disease; or may convey the seriousness of the disease. Colors and shapes may be used in various ways, e.g. to represent classes of diseases, compounds or environmental factors, such as steroids, lipids, polypeptides, polynucleotides, and the like; species of origin or gene families; signaling pathways; and the like.
Such additional information may also be conveyed by the use of multiple visualization windows. In addition to the graphic display of clustering information, the windows may contain text annotation of the profile; different spatial views of the matrix, different features, selected regions, and the like.
The risk analysis may be implemented in hardware or software, or a combination of both. In one embodiment of the invention, a machine-readable storage medium is provided, the medium comprising a data storage material encoded with machine readable data which, when using a machine programmed with instructions for using said data, is capable of displaying a any of the datasets and data comparisons of this invention. Such data may be used for a variety of purposes. Preferably, the invention is implemented in computer programs executing on programmable computers, comprising a processor, a data storage system (including volatile and non-volatile memory and/or storage elements), at least one input device, and at least one output device. Program code is applied to input data to perform the functions described above and generate output information. The output information is applied to one or more output devices, in known fashion. The computer may be, for example, a personal computer, microcomputer, or workstation of conventional design.
Each program is preferably implemented in a high level procedural or object oriented programming language to communicate with a computer system. However, the programs can be implemented in assembly or machine language, if desired. In any case, the language may be a compiled or interpreted language. Each such computer program is preferably stored on a storage media or device (e.g., ROM or magnetic diskette) readable by a general or special purpose programmable computer, for configuring and operating the computer when the storage media or device is read by the computer to perform the procedures described herein. The system may also be considered to be implemented as a computer-readable storage medium, configured with a computer program, where the storage medium so configured causes a computer to operate in a specific and predefined manner to perform the functions described herein.
A variety of structural formats for the input and output means can be used to input and output the information in the computer-based systems of the present invention. One format for an output means test datasets possessing varying degrees of similarity to a trusted profile. Such presentation provides a skilled artisan with a ranking of similarities and identifies the degree of similarity contained in the test pattern.
The following methods were used in the examples that are described further below to illustrate embodiments of the present invention.
Before the subject invention is further described, it is to be understood that the invention is not limited to the particular embodiments of the invention described below, as variations of the particular embodiments may be made and still fall within the scope of the appended claims. It is also to be understood that the terminology employed is for the purpose of describing particular embodiments, and is not intended to be limiting. Instead, the scope of the present invention will be established by the appended claims.
In this specification and the appended claims, the singular forms “a,” “an,” and “the” include plural reference unless the context clearly dictates otherwise. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood to one of ordinary skill in the art to which this invention belongs.
All publications and patent applications cited in this specification are herein incorporated by reference as if each individual publication or patent application were specifically and individually indicated to be incorporated by reference. The citation of any publication is for its disclosure prior to the filing date and should not be construed as an admission that the present invention is not entitled to antedate such publication by virtue of prior invention.
Although the foregoing invention has been described in some detail by way of illustration and example for purposes of clarity of understanding, it is readily apparent to those of ordinary skill in the art in light of the teachings of this invention that certain changes and modifications may be made thereto without departing from the spirit or scope of the appended claims.
The above-described analytical methods may be embodied as a program of instructions executable by computer to perform the different aspects of the invention. Any of the techniques described above may be performed by means of software components loaded into a computer or other information appliance or digital device. When so enabled, the computer, appliance or device may then perform the above-described techniques to assist the analysis of sets of values associated with a plurality of genes in the manner described above, or for comparing such associated values. The software component may be loaded from a fixed media or accessed through a communication medium such as the internet or other type of computer network. The above features are embodied in one or more computer programs may be performed by one or more computers running such programs.
The following examples are offered by way of illustration and not by way of limitation.
The cost of genomic information has fallen steeply, but the clinical translation of genetic risk estimates remains unclear. We aimed to undertake an integrated analysis of a complete human genome in a clinical context. We assessed a patient with a family history of vascular disease and early sudden death. Clinical assessment included analysis of this patient's full genome sequence, risk prediction for coronary artery disease, screening for causes of sudden cardiac death, and genetic counseling. Genetic analysis included the development of novel methods for the integration of whole genome and clinical risk. Disease and risk analysis focused on prediction of genetic risk of variants associated with mendelian disease, recognized drug responses, and pathogenicity for novel variants. We queried disease-specific mutation databases and pharmacogenomics databases to identify genes and mutations with known associations with disease and drug response. We estimated post-test probabilities of disease by applying likelihood ratios derived from integration of multiple common variants to age-appropriate and sex-appropriate pretest probabilities. We also accounted for gene-environment interactions and conditionally dependent risks.
Analysis of 2.6 million single nucleotide polymorphisms and 752 copy number variations showed increased genetic risk for myocardial infarction, type 2 diabetes, and some cancers. We discovered rare variants in three genes that are clinically associated with sudden cardiac death—TMEM43, DSP, and MYBPC3. A variant in LPA was consistent with a family history of coronary artery disease. The patient had a heterozygous null mutation in CYP2C19 suggesting probable clopidogrel resistance, several variants associated with a positive response to lipid-lowering therapy, and variants in CYP4F2 and VKORC1 that suggest he might have a low initial dosing requirement for warfarin. Many variants of uncertain importance were reported.
Our results show that whole-genome sequencing can yield useful and clinically relevant information for individual patients.
Methods
Patient.
A patient with a family history of vascular disease and early sudden death was assessed at Stanford's Center for Inherited Cardiovascular Disease by a cardiologist and a board-certified genetic counselor. We took the patient's medical history and he was clinically assessed. A four-generation pedigree was drawn. In view of his family history, he underwent electrocardiography, an echocardiogram, and a cardiopulmonary exercise test.
Genome Analysis.
Technical details of genome sequencing for this patient have been described previously. In brief, genomic DNA was purified from 2 mL of whole blood and sequenced with a Heliscope (Helicos BioSciences, Cambridge, Mass., USA) genome sequencer. We mapped sequence data to the National Center for Biotechnology Information reference human genome build 36 using the open-source aligner IndexDP (Helicos BioSciences, Cambridge, Mass., USA). Base calling was done with the UMKA algorithm. A subset of single nucleotide polymorphism calls were independently validated with the Illumina BeadArray (San Diego, Calif., USA) and all variants reported here and discussed with the patient were validated with Sanger sequencing. A subset of copy number variation calls were independently validated with digital PCR.
Disease and Risk Analysis.
Analysis focused on four areas: (i) variants associated with genes for mendelian disease; (ii) novel mutations; (iii) variants known to modulate response to pharmacotherapy; and (iv) single nucleotide polymorphisms previously associated with complex disease. Database queries, biophysical prediction algorithms, and analyses of non-coding regions were used to screen for rare and novel variants in the genome. We examined disease specific mutation databases, the human genome mutation database, and Online Mendelian Inheritance in Man to identify genes and mutations with known associations to monogenic diseases. We applied prediction algorithms to weight the likelihood of variant pathogenicity on the basis of allele frequency, conservation, and protein domain disruption. Additionally, we developed algorithms to index variants affecting or creating start sites, stop sites, splice sites, and microRNAs (
To integrate common variant genetic risk across a range of human disease, we built a manually curated disease and single-nucleotide-polymorphism database. Diseases and phenotypes were mapped to Unified Medical Language System Concept Unique Identifiers. Since strand direction was variably reported between studies, we identified strand direction by comparing with major or minor alleles in the appropriate HapMap population. Odds ratios were available for allele comparisons in most cases; however, to generate a medically relevant post-test probability of disease from integrated environmental and genetic risk, we calculated likelihood ratios (LRs) for the most important single nucleotide polymorphism from every haplotype block. Pre-test probability was derived from published sources and the LR was applied to the pre-test odds of disease, which were calculated from age appropriate and sex-appropriate population prevalence. Investigators did not always provide frequency data for genotype that allowed calculation of the LR. The study was reviewed by the institutional review board of Stanford University and the patient gave written consent. The patient received education and counseling before signing the consent form and throughout testing and follow-up.
Results
The patient was a 40-year-old man who presented with a family history of coronary artery disease and sudden death. His medical history was not clinically significant and the patient exercised regularly without symptoms. He was taking no prescribed medications and appeared well. Clinical characteristics were within normal limits (table 1). Electrocardiography showed sinus rhythm, normal axis, and high praecordial voltage with early repolarisation. An echocardiogram revealed normal right and left ventricular size, systolic, diastolic, and valvular function. There were no wall motion abnormalities during maximum exercise and 1-5 mm ST depression was upsloping. Maximum oxygen uptake was 50 mL/kg per min. A four-generation family pedigree (
Sequencing of genomic DNA resulted in an output of 148 GB of raw sequence, with an average read length of 33 bases. Base calling detected 2.6 million single nucleotide polymorphisms and 752 copy number variations. An important benefit of sequencing compared with DNA chip-based methods of genotyping is the identification of rare or novel variants. We searched for evidence of rare or novel variants that would predispose the patient or his family to disease. Specific to cardiovascular disease, we discovered rare variants in three genes that are clinically associated with sudden cardiac death—TMEM43, DSP, and MYBPC3. The MYBPC3 variant, encoding an arginine-to-glutamine change at position 326 of the cardiac myosin-binding protein C, was originally associated with late-onset hypertrophic cardiomyopathy. Subsequently, this variant has also been reported in several independent control populations without known hypertrophic cardiomyopathy, suggesting that it might be benign. Mutations in TMEM4330 or DSP31 have been associated with familial arrhythmogenic right-ventricular dysplasia or cardiomyopathy. Review of previous clinical assessment of extended family members revealed minor criteria for this disease in one first cousin, whose son died suddenly in his teens. By contrast with the findings for the identified rare MYBPC3 variant, the TMEM43 variant, encoding a methionine-to-valine change at position 41 of transmembrane protein 43, has not been previously published, but was seen in one of 150 probands with arrhythmogenic right-ventricular dysplasia or cardiomyopathy. The identified DSP variant, encoding an arginine-to-histidine change to aminoacid 1838 of the desmoplakin protein, is entirely novel. Control populations from clinical testing laboratories (more than 1000 total chromosomes) did not contain either the DSP or TMEM43 variants.
Analysis of the patient's genome revealed three novel and potentially damaging variants in two related genes that were previously associated with development of hemochromatosis. Subsequent to these findings, detailed review of personal and family history did not identify a history of hemochromatosis in the patient or family members. Echocardiogram results and liver function tests did not show evidence of the disease. Justification for continued surveillance and testing with serum iron studies was explored with the patient.
Additionally, the patient had a novel stop mutation in a gene implicated in hyperparathyroidism and parathyroid tumors. This variant might increase probability of future development of hyperparathyroidism or parathyroid tumors through a loss-of-heterozygosity mechanism. Consistent with a variant in a gene previously associated with osteoarthritis, there was a family history of osteoarthritis and the patient reported chronic knee pain without a formal diagnosis.
We noted 63 clinically relevant previously described pharmacogenomic variants and six novel, non-conservative, amino acid changing single nucleotide polymorphisms in genes that are important for drug response. There was a heterozygous null mutation in CYP2C19, the gene product of which is important for metabolism of many drugs, including proton-pump inhibitors, antiepileptic drugs, and the antiplatelet agent clopidogrel. Notably, the rate of cardiovascular events is raised in patients with CYP2C19 loss-of-function mutations who take clopidogrel. Additionally, the patient had two types of distinct genetic variations related to decreased maintenance dosing of warfarin. The patient had the single most important variant in VKORC1 associated with a low maintenance dose, and was homozygous for a CYP4F2 single nucleotide polymorphism that is associated with reduced dosing. Thus, if prescription of warfarin became necessary, loading could be individually tailored for this patient, with lowered expected doses. The patient had several variants that are associated with good response to statins (including reduced risk of myopathy) and one variant suggesting that he might need a raised dose to achieve a good response. Finally, the patient was wild type (with no copy number variations) for genes for important drug-metabolizing enzymes (CYP2D6, CYP2C9, and CYP3A4) affecting hundreds of drug responses.
Although genome-wide association studies have provided strong association of many common variants with disease, integration of these small odds ratios in the context of the individual patient remains challenging. In particular, additive or multiplicative models of even strongly associated single nucleotide polymorphisms can add little to the classified status of the patient. Furthermore, these approaches take no account of previous probability of disease. To counter some of these concerns, we adopted established methods from within evidence-based medicine that have rarely been applied to clinical genetics. We estimated pre-test probabilities from referenced sources for 121 diseases. Of the 55 diseases for which we could estimate a post-test probability, genetic risk was consistently increased (LR>2) for eight diseases and decreased (<0-5) for seven diseases (
Raised genetic risk did not always translate into high post-test probability. Post-test probabilities that were an order of magnitude higher or lower than pre-test probabilities were rare. Any decision towards acting on these predictions will necessarily be a function of the post-test probability threshold for action (e.g., the post-test probability of type 2 diabetes), the consequences of action (e.g., regular testing for fasting blood sugar), and the usefulness and effectiveness of action.
Increased genetic risk for myocardial infarction took the form of five single nucleotide polymorphisms associated with susceptibility to myocardial infarction and two protective polymorphisms (
We placed disease-associated genetic risk into the context of environmental and behavioral modifiers, as well as predisposing disorders (
During genetic counseling, we discussed the possibility that clinical assessment incorporating a personal genome might uncover high risk of a serious disease, including some for which there is no treatment. Additionally, we described the reproductive implications of heterozygous status for autosomal recessive diseases such as cystic fibrosis, which might not be predictable from family history (table 2,
We provide an approach to comprehensive analysis of a human genome in a defined clinical context. We assessed whole-genome genetic risk, focusing on variants in genes that are associated with mendelian disease, novel and rare variants across the genome, and variants of pharmacogenomic importance. Additionally, we developed an approach to the integration of disease risk across several common polymorphisms. The results provide proof that clinically meaningful information can be derived about disease risk and response to drugs in patients with whole genome sequence data.
Prominent aspects of the patient's family history (
Our patient reported a prominent family history of vascular disease including aortic aneurysm and coronary artery disease (
In view of a predisposition to coronary artery disease and other diseases on which risk is conditionally dependent (
By contrast, our patient did not report a family history of hemochromatosis or parathyroid tumors, yet has some genetic risk for these disorders. In consideration of future screening studies, integrated clinical and genetic risks were assessed.
Important limitations remain in our ability to comprehensively integrate genetic information into clinical care. For example, a comprehensive database of rare mutations is needed. Since risk estimates change as studies are completed, a continually updated pipeline is necessary. There are imperfections in all human genomes published to date—false positive and false negative SNP calls, incomplete measurement of structural variation, and little direct haplotype data. Finally, gene-environment interactions are challenging to quantify and have been little studied.
As whole-genome sequencing becomes increasingly widespread, availability of genomic information will no longer be the limiting factor in application of genetics to clinical medicine. Development of methods integrating genetic and clinical data will assist clinical decision making and represent a large step towards individualized medicine. The transition to a new era of genome-informed medical care will need a team approach incorporating medical and genetics professionals, ethicists, and health-care delivery organizations. Further evidence is provided by Dewey et al. (2011) Phased Whole-Genome Genetic Risk in a Family Quartet Using a Major Allele Reference Sequence. PLoS Genet 7(9):e1002280.
Methods
Rare variants. We searched the set of 2.6 million SNPs from the patient for rare variants in coding and important non-coding regions. The heuristic approach to this strategy is presented in
Manual curation was performed for each novel non-synonymous coding variant associated with a known or suspected disease gene, as found in Online Mendelian Inheritance in Man (OMIM). For inherited rare cardiovascular diseases, novel variants in genes currently analyzed by CLIA approved laboratories for clinical genetic testing were evaluated. Manual search of individual international locus specific mutation databases, curated dbSNP and OMIM entries, the public version of the Human Gene Mutation Database, and the published medical literature (PubMed) was performed to determine if prior information was available for putatively novel variants.
SIFT and PolyPhen were used to evaluate the potential effect of novel changes on predicted protein structure and/or the degree of evolutionary conservation of the altered base pair. Putative mRNA splice sites based on consensus coding determination sequences and the sequence of known mature microRNA coding regions were also searched for rare variants by queries matching splice site and miRNA databases to patient variants by chromosome location. This search yielded 394 splice site variants and one rare miRNA variant, which was not found in the putative miRNA seed sequence.
Pharmacogenomics. The compilation of drug-related genotype-phenotype associations was drawn from the Pharmacogenomics Knowledgebase (PharmGKB.org) whose mission is to collect, encode, and disseminate knowledge about the impact of human genetic variations on drug response. From the larger pool of gene-drug-disease literature annotations in PharmGKB, over 2500 manually curated genotypic associations were reduced to approximately 650 genotypic associations with drugs. The PharmGKB curators reviewed the drug-related variant annotations for their clinical relevance to this patient (e.g., patient's genotype, gender, age, family and medical history). A level of evidence was assigned to each variant annotation based on a clinician's appraisal of the impact of the variant based upon the strength of the results in the literature, the effect size of the phenotype, and the availability of drugs to use as alternative treatment. Results include literature evidence for the likelihood of good drug response, inefficacy, side effects, or dosing recommendations based upon the patient's genotype (Table 3).
Disease risk. Disease-associated SNP database. To analyze the disease risk across the spectrum of human disease, a high-quality disease-associated SNP database was built. Starting with a list of all SNPs in dbSNP that were measured in the HapMap 2 and 3 projects, relevant papers were found by searching for rsIDs from within the abstracts of all papers in MEDLINE, eliminating papers covering non-human genetics. A total of 2,671 papers were manually curated and a database generated, including from each full-text paper the disease name (e.g. coronary artery disease), specific phenotype (e.g. acute coronary syndrome in coronary artery disease), study population (e.g. Finnish), case and control population (e.g. 2508 patients with angiographically proven coronary artery disease), genotyping technology, major/minor/risk alleles, odds ratio, 95% confidence interval of the odds ratio, published p-value, and genetic model for each included, statistically significant genotype comparison, including those involving any non-HapMap SNPs within these papers. If not present from the search-based sample, disease associations were recorded from the full text of all papers referenced in the professional version of the Human Gene Mutation Database (HGMD) and the NHGRI GWAS catalog that included dbSNP IDs and included in the association database.
Studies were then separated into two categories: case/control studies (e.g. coronary artery disease vs. control) and cohort studies (e.g. Increase in P-selectin glycoprotein ligand (SELPLG) levels in coronary artery disease vs. control). For each study, the frequency of each genotype and allele in the case and control populations was recorded. As of this writing, 55,258 associations, including 37,409 genotype comparisons, 15,333 case/control genotype frequencies and 7,648 cohort genotype frequencies have been recorded. In the genotype comparison data, 15,333 out of 37,409 total comparisons (40.99%) contained the odds ratio for a specific genotype comparison.
Categorizing studies on similar diseases and phenotypes. To enable the integration of multiple studies on similar diseases and phenotypes, the disease/phenotype names in our association database were mapped to the Unified Medical Language System (UMLS) Concept Unique Identifiers (CUIs). To improve the quality of matching disease names to UMLS CUI, the best descriptive disease names from the Medical Subject Heading (MeSH) terms were manually selected and associated with each paper during curation, and the matching. UMLS CUIs were then manually examined.
The 55,258 identified associations were categorized into 813 different diseases and phenotypes for 9,649 specific dbSNP IDs. As a comparison, the professional version of HGMD contains 3,725 disease-associated dbSNP IDs, and does not include any genotype comparison information, such as the odds ratio. The NHGRI GWAS catalog includes 2062 records from 379 papers, relating 1683 dbSNP IDs to 251 diseases and phenotypes, including the odds ratio, but does not include corresponding genotype comparisons.
Identifying strand direction in association studies. It was noted that genotypes at the same SNP were reported in both the positive and negative strand in different studies, especially for genes transcribed in the negative strand. For example, the alleles for rs1004819 in IL23R was reported as UT for Crohn's disease and G/A for Ankylosing spondylitis. Strand information was rarely found to be reported in publications, so an algorithm to correctly identify the strand direction automatically was developed to compare study reported major/minor alleles with the major/minor alleles as found in the most similar population to the patient available in HapMap data. A total of 6,196 records were identified that reported genotypes in the negative strand.
Calculating likelihood ratio (LR) of disease risk for SNP genotypes. To calculate the likelihood ratio (LR), the genetic association database was limited to the set of case/control studies reporting associations with diseases in populations that included Caucasian males. We retrieved 4,137 LR for 735 SNPs on 141 diseases, curated from a total of 480 publications. For every disease SNP, we calculated the LR for each genotype using the following equation.
Most publications and related, publicly available datasets report the odds ratios for allele comparisons without exact genotype frequency data, and thus preclude calculation of a likelihood ratio from available data. As a result, many potentially relevant publications were excluded from the provided LR estimates. For this work, likelihood ratios were calculated using only the genotype frequencies reported in the literature.
Pre-test probability for diseases. Pre-test probabilities of lifetime risk of disease were calculated for a wide range of conditions for a person matching the patient's characteristics (age, sex, and ethnicity) using a combination of sources
For some conditions, it was possible to match the attributes of the patient very closely, while for many conditions, prevalence and lifetime risk in the general population or in a non-identical (but similar) age or demographic group were used to infer a pre-test probability. The majority of genetic association studies from which LR were derived were case-control, rather than cohort studies which might allow better risk estimates of future disease. Pre and posttest probabilities were used, with recognition of the limitations of this methodology, to provide a more useful construct to display and determine future risk than a simple display of odds ratios alone.
Calculating post-test probability of disease risk for the patient. The patient's full genome sequence was mapped to dbSNP build 129 from the UCSC genome browser to identify dbSNP IDs. Unmapped positions were mapped to dbSNP build 130 to identify additional dbSNP IDs. The LR for the patient's genotype was then determined at each SNP for each disease. For the LR reported in the negative strands, the genotypes were converted into positive strand genotypes before retrieval. We retrieved 1,141 LR on 121 diseases for 524 SNPs.
For each of 121 diseases, the post-test probability of developing disease for the patient was calculated as follows. First, for SNPs with multiple LR from multiple studies, the mean LR was calculated, weighted by the square root of sample sizes. Second, the human genome was partitioned into Haplotype blocks using HapBlock software; for each haplotype block, the highest LR SNP was used. Third, LR from all SNPs were multiplied to report the cumulative LR for the patient. The cumulative LR was calculated for 121 diseases. Finally, the pre-test probability was translated into pre-test odds, multiplied by the cumulative LR to get post-test odds, and then converted to the post-test probability using the following equation. Post-test probabilities were calculable for 55 diseases, while post-test probabilities were calculated using more than one SNP for 32 diseases (
Calculating the odds ratios of disease risk. Calculation of the LR of disease risk requires the frequency of three genotypes in the case and control populations. This was available for 1016 out of 2671 papers. In order to include most published genetic associations, the disease risks of the patient were calculated using the odds ratio as extracted from the literature. The allele and genotype frequencies in the case and control population were manually checked to identify the actual comparisons. Study authors were contacted when reported associations and genotype frequencies were discordant from population frequencies. Genetic associations were then retrieved from the disease-associated SNP database for each SNP. Studies reporting associations with diseases (as opposed to other human traits and drug effects) in Caucasian male population, and to those where odds ratios for single SNPs were used for further analysis. A total of 6,958 genetic associations for 2,325 dbSNPs on 235 diseases, curated from 900 publications were used. For each disease/SNP/study combination, the best matching record yielding the highest odds ratio for the patient as compared to the reference sequence was used. In many disease-associated SNPs, the patient has the same genotype as the reference sequence, giving an odds ratio of 1 for the patient versus the reference.
Gene-environment interaction and conditionally dependent risk. A database of established links between diseases and known aetiological factors was obtained as previously described. In brief, Medical Subject Heading (MeSH) annotations in the MEDLINE database were mapped to the Unified Medical Language System (UMLS) semantic network, keeping only those publications with annotations representing diseases and environmental factors. 18 Publications that contained a major MeSH term and sub-heading
(indicating a primary focus of the paper) on the “etiology” or “chemical induction” of a disease were selected. These publications also had to contain a major MeSH term and sub-heading on the “adverse effect” or “complication” or “toxicity” or “poisoning” related to an environmental factor. The frequency of finding each disease-factor association across the entire MEDLINE corpus was noted. For the purpose of constructing the aetiological “circuit diagram”, this database was filtered to aetiological links that were curated from 10 or more independent publications. For example, if 27 publications in MEDLINE contained the MeSH annotations “Renal Dialysis/adverse effects” and “Atherosclerosis/etiology”, and if only 25 publications with “Dietary Fats/adverse effects” and “Atherosclerosis/etiology”, then dialysis was considered as an environmental factor slightly more strongly associated with atherosclerosis than dietary fats.
| Entry |
|---|
| Ashley et al. Clinical assessment incorporating a personal genome The Lancet vol. 375 pp. 1525-1535 and Supplementary webappendix pp. 1-20 (2010). |
| Ashley; et al., “Clinical assessment incorporating a personal genome”, Lancet (2010), 375(9725):1525-1535. |
| Capriotti; et al., “Predicting the insurgence of human genetic diseases associated to single point protein mutations with support vector machines and evolutionary information”, Bioinformatics (2006), 22(22):2729-2734. |
| Jakobsdottir; et al., “Interpretation of Genetic Association Studies: Markers with Replicated Highly Significant Odds Ratios May Be Poor Classifiers”, PLOS Genetics (2009), 5(2):e1000337, 8 pages. |
| Jegga; et al., “PolyDoms: a whole genome database for the identification of non-synonymous coding SNPs with the potential to impact disease”, Nucleic, Acids Research (2007), 35:D700-D706. |
| Pushkarev; et al., “Single-molecule sequencing on an individual human genome”, Nature Biotechnology (2009), 27 (9):847-852. |
| Klein; et al., “Integrating genotype and phenotype information: an overview of the PharmGKB project”, The Pharmacogenomics Journal (2001), 1:167-170. |
| Tucker; et al., “Massively Parallel Sequencing: The Next Big Thing in Genetic Medicine”, The American Journal of Human Genetics (2009), 85:142-154. |
| UKPMC Funders Group, “Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls”, Nature (2007), 447(7145):661-678. |
| Number | Date | Country | |
|---|---|---|---|
| 20120101736 A1 | Apr 2012 | US |
| Number | Date | Country | |
|---|---|---|---|
| 61406447 | Oct 2010 | US |
