Methods and Systems for Genomic Analysis Using Ancestral Data

Information

  • Patent Application
  • 20090099789
  • Publication Number
    20090099789
  • Date Filed
    September 26, 2008
    15 years ago
  • Date Published
    April 16, 2009
    15 years ago
Abstract
The present disclosure provides methods and systems for assessing an individual's genotype correlations to a phenotype by analyzing the individual's genomic profile and using ancestral data to determine the correlations between genotypes and phenotypes.
Description
BACKGROUND

Sequencing of the human genome and other recent developments in human genomics has revealed that the genomic makeup between any two humans can be over 99.9% similarity. The relatively small number of variations in DNA between individuals gives rise to differences in phenotypic traits, and is related to many human diseases, susceptibility to various diseases, and response to treatment of disease. Variations in DNA between individuals occur in both coding and non-coding regions, and include changes in bases at a particular locus in genomic DNA sequences, as well as insertions and deletions of DNA. Changes that occur at single base positions in the genome are referred to as single nucleotide polymorphisms, or “SNPs.”


While SNPs are relatively rare in the human genome, they account for a majority of DNA sequence variations between individuals, occurring approximately once every 1,200 base pairs in the human genome (see International HapMap Project, www.hapmap.org). As more human genetic information becomes available, the complexity of SNPs is beginning to be understood. In turn, the occurrences of SNPs in the genome are becoming correlated to the presence of and/or susceptibility to various diseases and conditions.


As these correlations and other advances in human genetics are being made, medicine and personal health in general are moving toward a customized approach in which a patient will make appropriate medical and other choices in consideration of his or her genomic information, among other factors. An important factor that may affect considerations is an individual's ancestral data (ancestry) or ethnicity. For example, different populations may have different linkage disequilibrium patterns due to various possible reasons such as variation in recombination rates, selection pressure, or population bottleneck. Thus, if a study has been done on population A, yielding a specific odds ratio in that population for a genetic variation correlated with a phenotype, the same odds ratio cannot be assumed in population B. Thus, there is a need to provide individuals and their care-givers with information specific to the individual's personal genome, incorporating ancestral data, toward providing personalized medical and other decisions.


SUMMARY

The present disclosure provides a method of assessing genotype correlations to a phenotype of an individual comprising: a) obtaining a genetic sample of the individual, b) generating a genomic profile for the individual, c) determining the individual's genotype correlations with phenotypes by comparing the individual's genomic profile to a current database of human genotype correlations with phenotypes, d) reporting the results from step c) to the individual or a health care manager of the individual, e) updating the database of human genotype correlations with an additional human genotype correlation as the additional human genotype correlation becomes known, f) updating the individual's genotype correlations by comparing the individual's genomic profile from step c) or a portion thereof to the additional human genotype correlation and determining an additional genotype correlation of the individual, and g) reporting the results from step f) to the individual or the health care manager of the individual.


The present disclosure further provides a business method of assessing genotype correlations of an individual comprising: a) obtaining a genetic sample of the individual; b) generating a genomic profile for the individual; c) determining the individual's genotype correlations by comparing the individual's genomic profile to a database of human genotype correlations; d) providing results of the determining of the individual's genotype correlations to the individual in a secure manner; e) updating the database of human genotype correlations with an additional human genotype correlation as the additional human genotype correlation becomes known; f) updating the individual's genotype correlations by comparing the individual's genomic profile or a portion thereof to the additional human genotype correlation and determining an additional genotype correlation of the individual; and g) providing results of the updating of the individual's genotype correlations to the individual of the health care manager of the individual.


Another aspect of the present disclosure is a method generating a phenotype profile for an individual comprising: a) providing a rule set comprising rules, each rule indicating a correlation between at least one genotype and at least one phenotype, b) providing a data set comprising genomic profiles of each of a plurality of individuals, wherein each genomic profile comprises a plurality of genotypes; c) periodically updating the rule set with at least one new rule, wherein the at least one new rule indicates a correlation between a genotype and a phenotype not previously correlated with each other in the rule set; d) applying each new rule to the genomic profile of at least one of the individuals, thereby correlating at least one genotype with at least one phenotype for the individual, and optionally, e) generating a report comprising the phenotype profile of the individual.


The present disclosure also provides a system comprising a) a rule set comprising rules, each rule indicating a correlation between at least one genotype and at least one phenotype; b) code that periodically updates the rule set with at least one new rule, wherein the at least one new rule indicates a correlation between a genotype and a phenotype not previously correlated with each other in the rule set; c) a database comprising genomic profiles of a plurality of individuals; d) code that applies the rule set to the genomic profiles of individuals to determine phenotype profiles for the individuals; and e) code that generates reports for each individual.


The present disclosure further provides a method of assessing genotype correlations of an individual comprising: (a) comparing (i) a first linkage disequilibrium (LD) pattern comprising a genetic variation correlated with a phenotype, wherein the first LD pattern is of a first population of individuals; and, (ii) a second LD pattern comprising the genetic variation, wherein the second LD pattern is of a second population of individuals; (b) determining a probability of the genetic variation being correlated with the phenotype in said second population from said comparing in (a); (c) assessing a genotype correlation of said phenotype from a genomic profile of the individual comprising using the probability of step (b); and, (d) reporting results comprising the genotype correlation from step c) to the individual or a health care manager of the individual. In some embodiments, the methods further comprise (e) updating said results with additional genetic variations.


The probability can be an odds ratio (OR), wherein the OR can be derived from a known OR. For example, the known OR can be for the genetic variation correlated with the phenotype for the first population, such as an OR published for a genetic variation, such as a SNP, in a scientific journal. In some embodiments, the first population and the second population have similar LD patterns. Also provided herein is a method of assessing genotype correlations of an individual comprising: (a) determining a causal genetic variation probability for each of a plurality of genetic variations in a first population of individuals; (b) identifying each of said probability in step (a) as a probability for each of said plurality of genetic variations in a second population of individuals; (c) assessing a genotype correlation from a genomic profile of the individual comprising using the probability of step (b); and, (d) reporting results comprising the genotype correlation from step (c) to the individual or a health care manager of the individual. In some embodiments, the methods further comprise (e) updating said results with additional genetic variations.


The known genetic variation, such as a SNP, can be a genetic variation with an OR published in a scientific journal. The probability can be an odds ratio (OR) and each of the genetic variations of step (a) can be proximal to a known genetic variation correlated to a phenotype in the first population. For example, each of the genetic variations can be in linkage disequilibrium to the known genetic variation.


In some embodiments of the methods and systems disclosed herein, the genotype correlation is reported as a GCI score. The second population is typically of an ancestry different from the first population, and the individual is of an ancestry of the second population. In some embodiments, the causal genetic variation is unknown. The genetic variation can be a single nucleotide polymorphism (SNP).


Another aspect of the present disclosure is transmission over a network, in a secure or non-secure manner, the methods and systems described above. The reporting can be through an on-line portal, by paper or by e-mail. The genomic profile used can be generated and from a genetic sample. A third party can generate the genomic profile, obtain the genetic sample, or both obtain the sample and generate the genomic profile. The genetic sample can be DNA or RNA and obtained from a biological sample selected from the group consisting of: blood, hair, skin, saliva, semen, urine, fecal material, sweat, and buccal sample. The genomic profile can be deposited into a secure database or vault. Furthermore the genomic profile can be a single nucleotide polymorphism profile, and in some embodiments, the genomic profile can comprise truncations, insertions, deletions, or repeats. The genomic profile can be generated by using a high density DNA microarray, RT-PCR, DNA sequencing, or a combination of techniques.


The method of the invention also includes the populations comprising any of the HapMap populations (YRI,CEU,CHB,JPT,ASW,CHD,GIH,LWK,MEX,MKK,TSI), or to any other population such as, but not limited to African American, Caucasian, Ashkenazi Jewish, Sepharadic Jewish, Indian, Pacific islanders, middle eastern, Druze, Bedouins, south Europeans, Scandinavians, eastern Europeans, North Africans, Basques, West Africans, or East Africans.


INCORPORATION BY REFERENCE

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.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 is a flow chart illustrating aspects of the method herein.



FIG. 2 is an example of a genomic DNA quality control measure.



FIG. 3 is an example of a hybridization quality control measure.



FIG. 4 are tables of representative genotype correlations from published literature with test SNPs and effect estimates. A-I) represents single locus genotype correlations; J) represents a two locus genotype correlation; K) represents a three locus genotype correlation; L) is an index of the ethnicity and country abbreviations used in A-K; M) is an index of the abbreviations of the Short Phenotype Names in A-K, the heritability, and the references for the heritability.



FIG. 5A-J are tables of representative genotype correlations with effect estimates.



FIG. 6A-F are tables of representative genotype correlations and estimated relative risks.



FIG. 7 is a sample report.



FIG. 8 is a schematic of a system for the analysis and transmission of genomic and phenotype profiles over a network.



FIG. 9 is a flow chart illustrating aspects of the business method herein.



FIG. 10 is a schematic of a published SNP in CEU (Caucasian ancestry/ethnicity) with a specific odds ratio cannot be assumed to be the same in a different population of a different ancestral background, YRI (Yoruban ancestry/ethnicity see HapMap project (http://hapmap.org/hapmappopulations.html.en)).





DETAILED DESCRIPTION

The present disclosure provides methods and systems for generating phenotype profiles based on a stored genomic profile of an individual or group of individuals, and for readily generating original and updated phenotype profiles based on the stored genomic profiles. Genomic profiles are generated by determining genotypes from biological samples obtained from individuals. Biological samples obtained from individuals may be any sample from which a genetic sample may be derived. Samples may be from buccal swabs, saliva, blood, hair, or any other type of tissue sample. Genotypes may then be determined from the biological samples. Genotypes may be any genetic variant or biological marker, for example, single nucleotide polymorphisms (SNPs), haplotypes, or sequences of the genome. The genotype may be the entire genomic sequence of an individual. The genotypes may result from high-throughput analysis that generates thousands or millions of data points, for example, microarray analysis for most or all of the known SNPs. In other embodiments, genotypes may also be determined by high throughput sequencing.


The genotypes form a genomic profile for an individual. The genomic profile is stored digitally and is readily accessed at any point of time to generate phenotype profiles. Phenotype profiles are generated by applying rules that correlate or associate genotypes with phenotypes. Rules can be made based on scientific research that demonstrates a correlation between a genotype and a phenotype. The correlations may be curated or validated by a committee of one or more experts. By applying the rules to a genomic profile of an individual, the association between an individual's genotype and a phenotype may be determined. The phenotype profile for an individual will have this determination. The determination may be a positive association between an individual's genotype and a given phenotype, such that the individual has the given phenotype, or will develop the phenotype. Alternatively, it may be determined that the individual does not have, or will not develop, a given phenotype. In other embodiments, the determination may be a risk factor, estimate, or a probability that an individual has, or will develop a phenotype.


The determinations may be made based on a number of rules, for example, a plurality of rules may be applied to a genomic profile to determine the association of an individual's genotype with a specific phenotype. The determinations may also incorporate factors that are specific to an individual, such as ethnicity, gender, lifestyle, age, environment, family medical history, personal medical history, and other known phenotypes. The incorporation of the specific factors may be by modifying existing rules to encompass these factors. Alternatively, separate rules may be generated by these factors and applied to a phenotype determination for an individual after an existing rule has been applied.


Phenotypes may include any measurable trait or characteristic, such as susceptibility to a certain disease or response to a drug treatment. Other phenotypes that may be included are physical and mental traits, such as height, weight, hair color, eye color, sunburn susceptibility, size, memory, intelligence, level of optimism, and general disposition. Phenotypes may also include genetic comparisons to other individuals or organisms. For example, an individual may be interested in the similarity between their genomic profile and that of a celebrity. They may also have their genomic profile compared to other organisms such as bacteria, plants, or other animals.


In another aspect of the disclosure information about the association of multiple genetic markers with one or more diseases or conditions is combined and analyzed to produce a Genetic Composite Index (GCI) score (such as described in PCT Publication No. WO2008/067551, which is herein incorporated by reference). This score incorporates known risk factors, as well as other information and assumptions such as the allele frequencies and the prevalence of a disease. The GCI can be used to qualitatively estimate the association of a disease or a condition with the combined effect of a set of Genetic markers. The GCI score can be used to provide people not trained in genetics with a reliable (i.e., robust), understandable, and/or intuitive sense of what their individual risk of a disease is compared to a relevant population based on current scientific research. The GCI score may be used to generate GCI Plus scores, as described in PCT Publication No. WO2008/067551. The GCI Plus score may contain all the GCI assumptions, including risk (such as lifetime risk), age-defined prevalence, and/or age-defined incidence of the condition. The lifetime risk for the individual may then be calculated as a GCI Plus score which is proportional to the individual's GCI score divided by the average GCI score. The average GCI score may be determined from a group of individuals of similar ancestral background, for example a group of Caucasians, Asians, East Indians, or other group with a common ancestral background. Groups may comprise of at least 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, or 60 individuals. In some embodiments, the average may be determined from at least 75, 80, 95, or 100 individuals. The GCI Plus score may be determined by determining the GCI score for an individual, dividing the GCI score by the average relative risk and multiplying by the lifetime risk for a condition or phenotype. For example, using data from PCT Publication No. WO2008/067551, such as FIG. 22 and/or FIG. 25 with information in FIG. 24 to calculate GCI Plus scores such as in FIG. 19.


The present disclosure encompasses using the GCI score as described herein, and one of ordinary skill in the art will readily recognize the use of GCI Plus scores or variations thereof, in place of GCI scores as described herein. In one embodiment a GCI score is generated for each disease or condition of interest. These GCI scores may be collected to form a risk profile for an individual. The GCI scores may be stored digitally so that they are readily accessible at any point of time to generate risk profiles. Risk profiles may be broken down by broad disease classes, such as cancer, heart disease, metabolic disorders, psychiatric disorders, bone disease, or age on-set disorders. Broad disease classes may be further broken down into subcategories. For example for a broad class such as a cancer, sub-categories of cancer may be listed such as by type (sarcoma, carcinoma or leukemia, etc.) or by tissue specificity (neural, breast, ovaries, testes, prostate, bone, lymph nodes, pancreas, esophagus, stomach, liver, brain, lung, kidneys, etc.).


In another embodiment a GCI score is generated for an individual, which provides them with easily comprehended information about the individual's risk of acquiring or susceptibility to at least one disease or condition. In one embodiment multiple GCI scores are generated for different diseases or conditions. In another embodiment at least one GCI score is accessible by an on-line portal. Alternatively, at least one GCI score may be provided in paper form, with subsequent updates also provided in paper form. In one embodiment access to at least one GCI score is provided to a subscriber, which is an individual who subscribes to the service. In an alternative embodiment access is provided to non-subscribers, wherein they may have limited access to at least one of their GCI scores, or they may have an initial report on at least one of their GCI scores generated, but updated reports will be generated only with purchase of a subscription. In another embodiment health care managers and providers, such as caregivers, physicians, and genetic counselors may also have access to at least one of an individual's GCI scores.


Together, the collection of correlated phenotypes determined for an individual comprises the phenotype profile for the individual. The phenotype profile may be accessible by an on-line portal. Alternatively, the phenotype profile as it exists at a certain time may be provided in paper form, with subsequent updates also provided in paper form. The phenotype profile may also be provided by an on-line portal. The on-line portal may optionally be a secure on-line portal. Access to the phenotype profile may be provided to a subscriber, which is an individual who subscribes to the service that generates rules on correlations between phenotypes and genotypes, determines the genomic profile of an individual, applies the rules to the genomic profile, and generates a phenotype profile of the individual. Access may also be provided to non-subscribers, wherein they may have limited access to their phenotype profile and/or reports, or may have an initial report or phenotype profile generated, but updated reports will be generated only with purchase of a subscription. Health care managers and providers, such as caregivers, physicians, and genetic counselors may also have access to the phenotype profile.


In another aspect of the disclosure the genomic profile may be generated for subscribers and non-subscribers and stored digitally but access to the phenotype profile and reports may be limited to subscribers. In another variation, both subscribers and non-subscribers may access their genotype and phenotype profiles, but have limited access, or have a limited report generated for non-subscribers, whereas subscribers have full access and may have a full report generated. In another embodiment, both subscribers and non-subscribers may have full access initially, or full initial reports, but only subscribers may access updated reports based on their stored genomic profile.


There may also be a basic subscription model. A basic subscription may provide a phenotype profile where the subscriber may choose to apply all existing rules to their genomic profile, or a subset of the existing rules, to their genomic profile. For example, they may choose to apply only the rules for disease phenotypes that are actionable. The basic subscription may have different levels within the subscription class. For example, different levels may be dependent on the number of phenotypes a subscriber wants correlated to their genomic profile, or the number of people that may access their phenotype profile. Another level of basic subscription may be to incorporate factors specific to an individual, such as already known phenotypes such as age, gender, or medical history, to their phenotype profile.


Still another level of the basic subscription may allow an individual to generate at least one GCI score for a disease or condition. A variation of this level may further allow an individual to specify for an automatic update of at least one GCI score for a disease or condition to be generated if their is any change in at least one GCI score due to changes in the analysis used to generate at least one GCI score. In some embodiments the individual may be notified of the automatic update by email, voice message, text message, mail delivery, or fax.


Subscribers may also generate reports that have their phenotype profile as well as information about the phenotypes, such as genetic and medical information about the phenotype. For example, the prevalence of the phenotype in the population, the genetic variant that was used for the correlation, the molecular mechanism that causes the phenotype, therapies for the phenotype, treatment options for the phenotype, and preventative actions, may be included in the report. In other embodiments, the reports may also include information such as the similarity between an individual's genotype and that of other individuals, such as celebrities or other famous people. The information on similarity may be, but are not limited to, percentage homology, number of identical variants, and phenotypes that may be similar. These reports may further contain at least one GCI score.


The report may also provide links to other sites with further information on the phenotypes, links to on-line support groups and message boards of people with the same phenotype or one or more similar phenotypes, links to an on-line genetic counselor or physician, or links to schedule telephonic or in-person appointments with a genetic counselor or physician, if the report is accessed on-line. If the report is in paper form, the information may be the website location of the aforementioned links, or the telephone number and address of the genetic counselor or physician. The subscriber may also choose which phenotypes to include in their phenotype profile and what information to include in their report. The phenotype profile and reports may also be accessible by an individual's health care manager or provider, such as a caregiver, physician, psychiatrist, psychologist, therapist, or genetic counselor. The subscriber may be able to choose whether the phenotype profile and reports, or portions thereof, are accessible by such individual's health care manager or provider.


The present disclosure may also include a premium level of subscription. The premium level of subscription maintains their genomic profile digitally after generation of an initial phenotype profile and report, and provides subscribers the opportunity to generate phenotype profiles and reports with updated correlations from the latest research. In another embodiment, subscribers have the opportunity to generate risk profile and reports with updated correlations from the latest research. As research reveals new correlations between genotypes and phenotypes, disease or conditions, new rules will be developed based on these new correlations and can be applied to the genomic profile that is already stored and being maintained. The new rules may correlate genotypes not previously correlated with any phenotype, correlate genotypes with new phenotypes, or modify existing correlations, or provide the basis for adjustment of a GCI score based on a newly discovered association between a genotype and disease or condition. Subscribers may be informed of new correlations via e-mail or other electronic means, and if the phenotype is of interest, they may choose to update their phenotype profile with the new correlation. Subscribers may choose a subscription where they pay for each update, or for a number of updates or an unlimited number of updates for a designated time period (e.g. three months, six months, or one year). Another subscription level may be where a subscriber has their phenotype profile or risk profile automatically updated, instead of where the individual chooses when to update their phenotype profile or risk profile, whenever a new rule is generated based on a new correlation.


In another aspect of the subscription, subscribers may refer non-subscribers to the service that generates rules on correlations between phenotypes and genotypes, determines the genomic profile of an individual, applies the rules to the genomic profile, and generates a phenotype profile of the individual. Referral by a subscriber may give the subscriber a reduced price on subscription to the service, or upgrades to their existing subscriptions. Referred individuals may have free access for a limited time or have a discounted subscription price.


Phenotype profiles and reports as well as risk profiles and reports may be generated for individuals that are human and non-human. For example, individuals may include other mammals, such as bovines, equines, ovines, canines, or felines. Subscribers, as used herein, are human individuals who subscribe to a service by purchase or payment for one or more services. Services may include, but are not limited to, one or more of the following: having their or another individual's, such as the subscriber's child or pet, genomic profile determined, obtaining a phenotype profile, having the phenotype profile updated, and obtaining reports based on their genomic and phenotype profile.


In another aspect of the disclosure, “field-deployed” mechanisms may be gathered from individuals to generate phenotype profiles for individuals. In preferred embodiments, an individual may have an initial phenotype profile generated based on genetic information. For example, an initial phenotype profile is generated that includes risk factors for different phenotypes as well as suggested treatments or preventative measures. For example, the profile may include information on available medication for a certain condition, and/or suggestions on dietary changes or exercise regimens. The individual may choose to see, or contact via a web portal or phone call, a physician or genetic counselor, to discuss their phenotype profile. The individual may decide to take a certain course of action, for example, take specific medications, change their diet, etc.


The individual may then subsequently submit biological samples to assess changes in their physical condition and possible change in risk factors. Individuals may have the changes determined by directly submitting biological samples to the facility (or associated facility, such as a facility contracted by the entity generating the genetic profiles and phenotype profiles us) that generates the genomic profiles and phenotype profiles. Alternatively, the individuals may use a “field-deployed” mechanism, wherein the individual may submit their saliva, blood, or other biological sample into a detection device at their home, analyzed by a third party, and the data transmitted to be incorporated into another phenotype profile. For example, an individual may have received an initial phenotype report based on their genetic data reporting the individual having an increased lifetime risk of myocardial infarction (MI). The report may also have suggestions on preventative measures to reduce the risk of MI, such as cholesterol lowering drugs and change in diet. The individual may choose to contact a genetic counselor or physician to discuss the report and the preventative measures and decides to change their diet. After a period of being on the new diet, the individual may see their personal physician to have their cholesterol level measured. The new information (cholesterol level) may be transmitted (for example, via the Internet) to the entity with the genomic information, and the new information used to generate a new phenotype profile for the individual, with a new risk factor for myocardial infarction, and/or other conditions.


The individual may also use a “field-deployed” mechanism, or direct mechanism, to determine their individual response to specific medications. For example, an individual may have their response to a drug measured, and the information may be used to determine more effective treatments. Measurable information include, but are not limited to, metabolite levels, glucose levels, ion levels (for example, calcium, sodium, potassium, iron), vitamins, blood cell counts, body mass index (BMI), protein levels, transcript levels, heart rate, etc., can be determined by methods readily available and can be factored into an algorithm to combine with initial genomic profiles to determine a modified overall risk estimate score.


The term “biological sample” refers to any biological sample from which a genetic sample of an individual can be isolated.


As used herein, a “genetic sample” refers to DNA and/or RNA obtained or derived from an individual.


As used herein, the term “genome” is intended to mean the full complement of chromosomal DNA found within the nucleus of a human cell. The term “genomic DNA” refers to one or more chromosomal DNA molecules occurring naturally in the nucleus of a human cell, or a portion of the chromosomal DNA molecules.


The term “genomic profile” refers to a set of information about an individual's genes, such as the presence or absence of specific SNPs or mutations. Genomic profiles include the genotypes of individuals. Genomic profiles may also be substantially the complete genomic sequence of an individual. In some embodiments, the genomic profile may be at least 60%, 80%, or 95% of the complete genomic sequence of an individual. The genomic profile may be approximately 100% of the complete genomic sequence of an individual. In reference to a genomic profile, “a portion thereof” refers to the genomic profile of a subset of the genomic profile of an entire genome.


The term “genotype” refers to the specific genetic makeup of an individual's DNA. The genotype may include the genetic variants and markers of an individual. Genetic markers and variants may include nucleotide repeats, nucleotide insertions, nucleotide deletions, chromosomal translocations, chromosomal duplications, or copy number variations. Copy number variation may include microsatellite repeats, nucleotide repeats, centromeric repeats, or telomeric repeats. The genotypes may also be SNPs, haplotypes, or diplotypes. A haplotype may refer to a locus or an allele. A haplotype is also referred to as a set of single nucleotide polymorphisms (SNPs) on a single chromatid that are statistically associated. A diplotype is a set of haplotypes. The term single nucleotide polymorphism or “SNP” refers to a particular locus on a chromosome which exhibits variability such as at least one percent (1%) with respect to the identity of the nitrogenous base present at such locus within the human population For example, where one individual might have adenosine (A) at a particular nucleotide position of a given gene, another might have cytosine (C), guanine (G), or thymine (T) at this position, such that there is a SNP at that particular position.


As used herein, the terminology “SNP genomic profile” refers to the base content of a given individual's DNA at SNP sites throughout the individual's entire genomic DNA sequence. A “SNP profile” can refer to an entire genomic profile, or may refer to a portion thereof, such as a more localized SNP profile which can be associated with a particular gene or set of genes.


The term “phenotype” is used to describe a quantitative trait or characteristic of an individual. Phenotypes include, but are not limited to, medical and non-medical conditions. Medical conditions include diseases and disorders. Phenotypes may also include physical traits, such as hair color, physiological traits, such as lung capacity, mental traits, such as memory retention, emotional traits, such as ability to control anger, ethnicity, such as ethnic background, ancestry, such as an individual's place of origin, and age, such as age expectancy or age of onset of different phenotypes. Phenotypes may also be monogenic, wherein it is thought that one gene may be correlated with a phenotype, or multigenic, wherein more than one gene is correlated with a phenotype.


A “rule” is used to define the correlation between a genotype and a phenotype. The rules may define the correlations by a numerical value, for example by a percentage, risk factor, or confidence score. A rule may incorporate the correlations of a plurality of genotypes with a phenotype. A “rule set” comprises more than one rule. A “new rule” may be a rule that indicates a correlation between a genotype and a phenotype for which a rule does not currently exist. A new rule may correlate an uncorrelated genotype with a phenotype. A new rule may also correlate a genotype that is already correlated with a phenotype to a phenotype it had not been previously correlated to. A “new rule” may also be an existing rule that is modified by other factors, including another rule. An existing rule may be modified due to an individual's known characteristics, such as ethnicity, ancestry, geography, gender, age, family history, or other previously determined phenotypes.


Use of “genotype correlation” herein refers to the statistical correlation between an individual's genotype, such as presence of a certain mutation or mutations, and the likelihood of being predisposed to a phenotype, such as a particular disease, condition, physical state, and/or mental state. The frequency with which a certain phenotype is observed in the presence of a specific genotype determines the degree of genotype correlation or likelihood of a particular phenotype. For example, as detailed herein, SNPs giving rise to the apolipoprotein E4 isoform are correlated with being predisposed to early onset Alzheimer's disease. Genotype correlations may also refer to correlations wherein there is not a predisposition to a phenotype, or a negative correlation. The genotype correlations may also represent an estimate of an individual to have a phenotype or be predisposed to have a phenotype. The genotype correlation may be indicated by a numerical value, such as a percentage, a relative risk factor, an effects estimate, or confidence score.


The term “phenotype profile” refers to a collection of a plurality of phenotypes correlated with a genotype or genotypes of an individual. Phenotype profiles may include information generated by applying one or more rules to a genomic profile, or information about genotype correlations that are applied to a genomic profile. Phenotype profiles may be generated by applying rules that correlate a plurality of genotypes with a phenotype. The probability or estimate may be expressed as a numerical value, such as a percentage, a numerical risk factor or a numerical confidence interval. The probability may also be expressed as high, moderate, or low. The phenotype profiles may also indicate the presence or absence of a phenotype or the risk of developing a phenotype. For example, a phenotype profile may indicate the presence of blue eyes, or a high risk of developing diabetes. The phenotype profiles may also indicate a predicted prognosis, effectiveness of a treatment, or response to a treatment of a medical condition.


The term risk profile refers to a collection of GCI scores for more than one disease or condition. GCI scores are based on analysis of the association between an individual's genotype with one or more diseases or conditions. Risk profiles may display GCI scores grouped into categories of disease. Further the Risk profiles may display information on how the GCI scores are predicted to change as the individual ages or various risk factors are adjusted. For example, the GCI scores for particular diseases may take into account the effect of changes in diet or preventative measures taken (smoking cessation, drug intake, double radical mastectomies, hysterectomies). The GCI scores may be displayed as a numerical measure, a graphical display, auditory feedback or any combination of the preceding.


As used herein, the term “on-line portal” refers to a source of information which can be readily accessed by an individual through use of a computer and internet website, telephone, or other means that allow similar access to information. The on-line portal may be a secure website. The website may provide links to other secure and non-secure websites, for example links to a secure website with the individual's phenotype profile, or to non-secure websites such as a message board for individuals sharing a specific phenotype.


The practice of the present disclosure may employ, unless otherwise indicated, conventional techniques and descriptions of molecular biology, cell biology, biochemistry, and immunology, which are within the skill of the art. Such conventional techniques include nucleic acid isolation, polymer array synthesis, hybridization, ligation, and detection of hybridization using a label. Specific illustrations of suitable techniques are exemplified and referenced herein. However, other equivalent conventional procedures can also be used. Other conventional techniques and descriptions can be found in standard laboratory manuals and texts such as Genome Analysis: A Laboratory Manual Series (Vols. I-IV), PCR Primer: A Laboratory Manual, Molecular Cloning: A Laboratory Manual (all from Cold Spring Harbor Laboratory Press); Stryer, L. (1995) Biochemistry (4th Ed.) Freeman, N.Y.; Gait, “Oligonucleotide Synthesis: A Practical Approach” 1984, IRL Press, London, Nelson and Cox (2000); Lehninger, Principles of Biochemistry 3rd Ed., W.H. Freeman Pub., New York, N.Y.; and Berg et al. (2002) Biochemistry, 5th Ed., W.H. Freeman Pub., New York, N.Y., all of which are herein incorporated in their entirety by reference for all purposes.


The methods of the present disclosure involve analysis of an individual's genomic profile to provide the individual with molecular information relating to a phenotype. As detailed herein, the individual provides a genetic sample, from which a personal genomic profile is generated. The data of the individual's genomic profile is queried for genotype correlations by comparing the profile against a database of established and validated human genotype correlations. The database of established and validated genotype correlations may be from peer-reviewed literature and further judged by a committee of one or more experts in the field, such as geneticists, epidemiologists, or statisticians, and curated. In preferred embodiments, rules are made based on curated genotype correlations and are applied to an individual's genomic profile to generate a phenotype profile. Results of the analysis of the individual's genomic profile, phenotype profile, along with interpretation and supportive information, are provided to the individual of the individual's health care manager, to empower personalized choices for the individual's health care.


The method of the disclosure is detailed as in FIG. 1, where an individual's genomic profile is first generated. An individual's genomic profile will contain information about an individual's genes based on genetic variations or markers. Genetic variations are genotypes, which make up genomic profiles. Such genetic variations or markers include, but are not limited to, single nucleotide polymorphisms, single and/or multiple nucleotide repeats, single and/or multiple nucleotide deletions, microsatellite repeats (small numbers of nucleotide repeats with a typical 5-1,000 repeat units), di-nucleotide repeats, tri-nucleotide repeats, sequence rearrangements (including translocation and duplication), copy number variations (both loss and gains at specific loci), and the like. Other genetic variations include chromosomal duplications and translocations as well as centromeric and telomeric repeats.


Genotypes may also include haplotypes and diplotypes. In some embodiments, genomic profiles may have at least 100,000, 300,000, 500,000, or 1,000,000 genotypes. In some embodiments, the genomic profile may be substantially the complete genomic sequence of an individual. In other embodiments, the genomic profile is at least 60%, 80%, or 95% of the complete genomic sequence of an individual. The genomic profile may be approximately 100% of the complete genomic sequence of an individual. Genetic samples that contain the targets include, but are not limited to, unamplified genomic DNA or RNA samples or amplified DNA (or cDNA). The targets may be particular regions of genomic DNA that contain genetic markers of particular interest.


In step 102 of FIG. 1, a genetic sample of an individual is isolated from a biological sample of an individual. Such biological samples include, but are not limited to, blood, hair, skin, saliva, semen, urine, fecal material, sweat, buccal, and various bodily tissues. In some embodiments, tissues samples may be directly collected by the individual, for example, a buccal sample may be obtained by the individual taking a swab against the inside of their cheek. Other samples such as saliva, semen, urine, fecal material, or sweat, may also be supplied by the individual themselves. Other biological samples may be taken by a health care specialist, such as a phlebotomist, nurse or physician. For example, blood samples may be withdrawn from an individual by a nurse. Tissue biopsies may be performed by a health care specialist, and kits are also available to health care specialists to efficiently obtain samples. A small cylinder of skin may be removed or a needle may be used to remove a small sample of tissue or fluids.


In some embodiments, kits are provided to individuals with sample collection containers for the individual's biological sample. The kit may also provide instructions for an individual to directly collect their own sample, such as how much hair, urine, sweat, or saliva to provide. The kit may also contain instructions for an individual to request tissue samples to be taken by a health care specialist. The kit may include locations where samples may be taken by a third party, for example kits may be provided to health care facilities who in turn collect samples from individuals. The kit may also provide return packaging for the sample to be sent to a sample processing facility, where genetic material is isolated from the biological sample in step 104.


A genetic sample of DNA or RNA may be isolated from a biological sample according to any of several well-known biochemical and molecular biological methods, see, e.g., Sambrook, et al., Molecular Cloning: A Laboratory Manual (Cold Spring Harbor Laboratory, New York) (1989). There are also several commercially available kits and reagents for isolating DNA or RNA from biological samples, such as those available from DNA Genotek, Gentra Systems, Qiagen, Ambion, and other suppliers. Buccal sample kits are readily available commercially, such as the MasterAmp™ Buccal Swab DNA extraction kit from Epicentre Biotechnologies, as are kits for DNA extraction from blood samples such as Extract-N-Amp™ from Sigma Aldrich. DNA from other tissues may be obtained by digesting the tissue with proteases and heat, centrifuging the sample, and using phenol-chloroform to extract the unwanted materials, leaving the DNA in the aqueous phase. The DNA can then be further isolated by ethanol precipitation.


In a preferred embodiment, genomic DNA is isolated from saliva. For example, using DNA self collection kit technology available from DNA Genotek, an individual collects a specimen of saliva for clinical processing. The sample conveniently can be stored and shipped at room temperature. After delivery of the sample to an appropriate laboratory for processing, DNA is isolated by heat denaturing and protease digesting the sample, typically using reagents supplied by the collection kit supplier at 50° C. for at least one hour. The sample is next centrifuged, and the supernatant is ethanol precipitated. The DNA pellet is suspended in a buffer appropriate for subsequent analysis.


In another embodiment, RNA may be used as the genetic sample. In particular, genetic variations that are expressed can be identified from mRNA. The term “messenger RNA” or “mRNA” includes, but is not limited to pre-mRNA transcript(s), transcript processing intermediates, mature mRNA(s) ready for translation and transcripts of the gene or genes, or nucleic acids derived from the mRNA transcript(s). Transcript processing may include splicing, editing and degradation. As used herein, a nucleic acid derived from an mRNA transcript refers to a nucleic acid for whose synthesis the mRNA transcript or a subsequence thereof has ultimately served as a template. Thus, a cDNA reverse transcribed from an mRNA, a DNA amplified from the cDNA, an RNA transcribed from the amplified DNA, etc., are all derived from the mRNA transcript. RNA can be isolated from any of several bodily tissues using methods known in the art, such as isolation of RNA from unfractionated whole blood using the PAXgene™ Blood RNA System available from PreAnalytiX. Typically, mRNA will be used to reverse transcribe cDNA, which will then be used or amplified for gene variation analysis.


Prior to genomic profile analysis, a genetic sample will typically be amplified, either from DNA or cDNA reverse transcribed from RNA. DNA can be amplified by a number of methods, many of which employ PCR. See, for example, PCR Technology: Principles and Applications for DNA Amplification (Ed. H. A. Erlich, Freeman Press, NY, N.Y., 1992); PCR Protocols: A Guide to Methods and Applications (Eds. Innis, et al., Academic Press, San Diego, Calif., 1990); Mattila et al., Nucleic Acids Res. 19, 4967 (1991); Eckert et al., PCR Methods and Applications 1, 17 (1991); PCR (Eds. McPherson et al., IRL Press, Oxford); and U.S. Pat. Nos. 4,683,202, 4,683,195, 4,800,159 4,965,188, and 5,333,675, and each of which is incorporated herein by reference in their entireties for all purposes.


Other suitable amplification methods include the ligase chain reaction (LCR) (for example, Wu and Wallace, Genomics 4, 560 (1989), Landegren et al., Science 241, 1077 (1988) and Barringer et al. Gene 89:117 (1990)), transcription amplification (Kwoh et al., Proc. Natl. Acad. Sci. USA 86:1173-1177 (1989) and WO88/10315), self-sustained sequence replication (Guatelli et al., Proc. Nat. Acad. Sci. USA, 87:1874-1878 (1990) and WO90/06995), selective amplification of target polynucleotide sequences (U.S. Pat. No. 6,410,276), consensus sequence primed polymerase chain reaction (CP-PCR) (U.S. Pat. No. 4,437,975), arbitrarily primed polymerase chain reaction (AP-PCR) (U.S. Pat. Nos. 5,413,909, 5,861,245) nucleic acid based sequence amplification (NABSA), rolling circle amplification (RCA), multiple displacement amplification (MDA) (U.S. Pat. Nos. 6,124,120 and 6,323,009) and circle-to-circle amplification (C2CA) (Dahl et al. Proc. Natl. Acad. Sci 101:4548-4553 (2004)). (See, U.S. Pat. Nos. 5,409,818, 5,554,517, and 6,063,603, each of which is incorporated herein by reference). Other amplification methods that may be used are described in, U.S. Pat. Nos. 5,242,794, 5,494,810, 5,409,818, 4,988,617, 6,063,603 and 5,554,517 and in U.S. Ser. No. 09/854,317, each of which is incorporated herein by reference.


Generation of a genomic profile in step 106 is performed using any of several methods. Generation of a genomic profile can be performed using any of several methods. Several methods are known in the art to identify genetic variations, and include, but are not limited to, DNA sequencing by any of several methodologies, PCR based methods, fragment length polymorphism assays (restriction fragment length polymorphism (RFLP), cleavage fragment length polymorphism (CFLP)) hybridization methods using an allele-specific oligonucleotide as a template (e.g., TaqMan assays and microarrays, further described herein), methods using a primer extension reaction, mass spectrometry (such as, MALDI-TOF/MS method), and the like, such as described in Kwok, Pharmocogenomics 1:95-100 (2000). Other methods include invader methods, such as monoplex and biplex invader assays (e.g. available from Third Wave Technologies, Madison, Wis. and described in Olivier et al., Nucl. Acids Res. 30:e53 (2002)).


In one embodiment, a high density DNA array is used for SNP identification and profile generation. Such arrays are commercially available from Affymetrix and Illumina (see Affymetrix GeneChip® 500K Assay Manual, Affymetrix, Santa Clara, Calif. (incorporated by reference); Sentrix® humanHap650Y genotyping beadchip, Illumina, San Diego, Calif.).


For example, a SNP profile can be generated by genotyping more than 900,000 SNPs using the Affymetrix Genome Wide Human SNP Array 6.0. Alternatively, more than 500,000 SNPs through whole-genome sampling analysis may be determined by using the Affymetrix GeneChip Human Mapping 500K Array Set. In these assays, a subset of the human genome is amplified through a single primer amplification reaction using restriction enzyme digested, adaptor-ligated human genomic DNA. As shown in FIG. 2, the concentration of the ligated DNA may then be determined. The amplified DNA is then fragmented and the quality of the sample determined prior to continuing with step 106. If the samples meet the PCR and fragmentation standards, the sample is denatured, labeled, and then hybridized to a microarray consisting of small DNA probes at specific locations on a coated quartz surface. The amount of label that hybridizes to each probe as a function of the amplified DNA sequence is monitored, thereby yielding sequence information and resultant SNP genotyping.


Use of the Affymetrix GeneChip 500K Assay is carried out according to the manufacturer's directions. Briefly, isolated genomic DNA is first digested with either a NspI or StyI restriction endonuclease. The digested DNA is then ligated with a NspI or StyI adaptor oligonucleotide that respectively anneals to either the NspI or StyI restricted DNA. The adaptor-containing DNA following ligation is then amplified by PCR to yield amplified DNA fragments between about 200 and 1100 base pairs, as confirmed by gel electrophoresis. PCR products that meet the amplification standard are purified and quantified for fragmentation. The PCR products are fragmented with DNase I for optimal DNA chip hybridization. Following fragmentation, DNA fragments should be less than 250 base pairs, and on average, about 180 base pairs, as confirmed by gel electrophoresis. Samples that meet the fragmentation standard are then labeled with a biotin compound using terminal deoxynucleotidyl transferase. The labeled fragments are next denatured and then hybridized into a GeneChip 250K array. Following hybridization, the array is stained prior to scanning in a three step process consisting of a streptavidin phycoerythin (SAPE) stain, followed by an antibody amplification step with a biotinylated, anti-streptavidin antibody (goat), and final stain with streptavidin phycoerythin (SAPE). After labeling, the array is covered with an array holding buffer and then scanned with a scanner such as the Affymetrix GeneChip Scanner 3000.


Analysis of data following scanning of an Affymetrix GeneChip Human Mapping 500K Array Set is performed according to the manufacturer's guidelines, as shown in FIG. 3. Briefly, acquisition of raw data using GeneChip Operating Software (GCOS) occurs. Data may also be acquired using Affymetrix GeneChip Command Console™. The acquisition of raw data is followed by analysis with GeneChip Genotyping Analysis Software (GTYPE). For purposes of the present disclosure, samples with a GTYPE call rate of less than 80% are excluded. Samples are then examined with BRLMM and/or SNiPer algorithm analyses. Samples with a BRLMM call rate of less than 95% or a SNiPer call rate of less than 98% are excluded. Finally, an association analysis is performed, and samples with a SNiPer quality index of less than 0.45 and/or a Hardy-Weinberg p-value of less than 0.00001 are excluded.


As an alternative to or in addition to DNA microarray analysis, genetic variations such as SNPs and mutations can be detected by other hybridization based methods, such as the use of TaqMan methods and variations thereof. TaqMan PCR, iterative TaqMan, and other variations of real time PCR (RT-PCR), such as those described in Livak et al., Nature Genet., 9, 341-32 (1995) and Ranade et al. Genome Res., 11, 1262-1268 (2001) can be used in the methods disclosed herein. In some embodiments, probes for specific genetic variations, such as SNPs, are labeled to form TaqMan probes. The probes are typically approximately at least 12, 15, 18 or 20 base pairs in length. They may be between approximately 10 and 70, 15 and 60, 20 and 60, or 18 and 22 base pairs in length. The probe is labeled with a reporter label, such as a fluorophore, at the 5′ end and a quencher of the label at the 3′ end: The reporter label may be any fluorescent molecule that has its fluorescence inhibited or quenched when in close proximity, such as the length of the probe, to the quencher. For example, the reporter label can be a fluorophore such as 6-carboxyfluorescein (FAM), tetracholorfluorescin (TET), or derivatives thereof, and the quencher tetramethylrhodamine (TAMRA), dihydrocyclopyrroloindole tripeptide (MGB), or derivatives thereof.


As the reporter fluorophore and quencher are in close proximity, separated by the length of the probe, the fluorescence is quenched. When the probe anneals to a target sequence, such as a sequence comprising a SNP in a sample, DNA polymerase with 5′ to 3′ exonuclease activity, such as Taq polymerase, can extend the primer and the exonuclease activity cleaves the probe, separating the reporter from the quencher, and thus the reporter can fluoresce. The process can be repeated, such as in RT-PCR. The TaqMan probe is typically complementary to a target sequence that is located between two primers that are designed to amplify a sequence. Thus, the accumulation of PCR product can be correlated to the accumulation of released fluorophore, as each probe can hybridize to newly generated PCR product. The released fluorophore can be measured and the amount of target sequence present can be determined. RT-PCR methods for high throughput genotyping, can be employed.


Genetic variations can also be identified by DNA sequencing. DNA sequencing may be used to sequence a substantial portion, or the entire, genomic sequence of an individual. Traditionally, common DNA sequencing has been based on polyacrylamide gel fractionation to resolve a population of chain-terminated fragments (Sanger et al., Proc. Natl. Acad. Sci. USA 74:5463-5467 (1977)). Alternative methods have been and continue to be developed to increase the speed and ease of DNA sequencing. For example, high throughput and single molecule sequencing platforms are commercially available or under development from 454 Life Sciences (Branford, Conn.) (Margulies et al., Nature 437:376-380 (2005)); Solexa (Hayward, Calif.); Helicos BioSciences Corporation (Cambridge, Mass.) (U.S. application Ser. No. 11/167,046, filed Jun. 23, 2005), and Li-Cor Biosciences (Lincoln, Nebr.) (U.S. application Ser. No. 11/118,031, filed Apr. 29, 2005).


After an individual's genomic profile is generated in step 106, the profile is stored digitally in step 108, such profile may be stored digitally in a secure manner. The genomic profile is encoded in a computer readable format to be stored as part of a data set and may be stored as a database, where the genomic profile may be “banked”, and can be accessed again later. The data set comprises a plurality of data points, wherein each data point relates to an individual. Each data point may have a plurality of data elements. One data element is the unique identifier, used to identify the individual's genomic profile. It may be a bar code. Another data element is genotype information, such as the SNPs or nucleotide sequence of the individual's genome. Data elements corresponding to the genotype information may also be included in the data point. For example, if the genotype information includes SNPs identified by microarray analysis, other data elements may include the microarray SNP identification number, the SNP rs number, and the polymorphic nucleotide. Other data elements may be chromosome position of the genotype information, quality metrics of the data, raw data files, images of the data, and extracted intensity scores.


The individual's specific factors such as physical data, medical data, ethnicity, ancestry, geography, gender, age, family history, known phenotypes, demographic data, exposure data, lifestyle data, behavior data, and other known phenotypes may also be incorporated as data elements. For example, factors may include, but are not limited to, individual's: birthplace, parents and/or grandparents, relatives' ancestry, location of residence, ancestors' location of residence, environmental conditions, known health conditions, known drug interactions, family health conditions, lifestyle conditions, diet, exercise habits, marital status, and physical measurements, such as weight, height, cholesterol level, heart rate, blood pressure, glucose level and other measurements known in the art The above mentioned factors for an individual's relatives or ancestors, such as parents and grandparents, may also be incorporated as data elements and used to determine an individual's risk for a phenotype or condition.


The specific factors may be obtained from a questionnaire or from a health care manager of the individual. Information from the “banked” profile can then be accessed and utilized as desired. For example, in the initial assessment of an individual's genotype correlations, the individual's entire information (typically SNPs or other genomic sequences across, or taken from an entire genome) will be analyzed for genotype correlations. In subsequent analyses, either the entire information can be accessed, or a portion thereof, from the stored, or banked genomic profile, as desired or appropriate.


Comparison of Genomic Profile with Database of Genotype Correlations.


In step 110, genotype correlations are obtained from scientific literature. Genotype correlations for genetic variations are determined from analysis of a population of individuals who have been tested for the presence or absence of one or more phenotypic traits of interest and for genotype profile. The alleles of each genetic variation or polymorphism in the profile are then reviewed to determine whether the presence or absence of a particular allele is associated with a trait of interest. Correlation can be performed by standard statistical methods and statistically significant correlations between genetic variations and phenotypic characteristics are noted. For example, it may be determined that the presence of allele A1 at polymorphism A correlates with heart disease. As a further example, it might be found that the combined presence of allele A1 at polymorphism A and allele B1 at polymorphism B correlates with increased risk of cancer. The results of the analyses may be published in peer-reviewed literature, validated by other research groups, and/or analyzed by a committee of experts, such as geneticists, statisticians, epidemiologists, and physicians, and may also be curated.


In FIGS. 4, 5, and 6 are examples of correlations between genotypes and phenotypes from which rules to be applied to genomic profiles may be based. For example, in FIGS. 4A and B, each row corresponds to a phenotype/locus/ethnicity, wherein FIGS. 4C through I contains further information about the correlations for each of these rows. As an example, in FIG. 4A, the “Short Phenotype Name” of BC, as noted in FIG. 4M, an index for the names of the short phenotypes, is an abbreviation for breast cancer. In row BC4, which is the generic name for the locus, the gene LSP1 is correlated to breast cancer. The published or functional SNP identified with this correlation is rs3817198, as shown in FIG. 4C, with the published risk allele being C, the nonrisk allele being T. The published SNP and alleles are identified through publications such as seminal publications as in FIGS. 4E-G. In the example of LSP1 in FIG. 4E, the seminal publication is Easton et al., Nature 447:713-720 (2007).


Alternatively, the correlations may be generated from the stored genomic profiles. For example, individuals with stored genomic profiles may also have known phenotype information stored as well. Analysis of the stored genomic profiles and known phenotypes may generate a genotype correlation. As an example, 250 individuals with stored genomic profiles also have stored information that they have previously been diagnosed with diabetes. Analysis of their genomic profiles is performed and compared to a control group of individuals without diabetes. It is then determined that the individuals previously diagnosed with diabetes have a higher rate of having a particular genetic variant compared to the control group, and a genotype correlation may be made between that particular genetic variant and diabetes.


In step 112, rules are made based on the validated correlations of genetic variants to particular phenotypes. Rules may be generated based on the genotypes and phenotypes correlated as listed in Table 1, for example. Rules based on correlations may incorporate other factors such as gender (e.g. FIG. 4) or ethnicity (FIGS. 4 and 5), to generate effects estimates, such as those in FIGS. 4 and 5. Other measures resulting from rules may be estimated relative risk increase such as in FIG. 6. The effects estimates and estimated relative risk increase may be from the published literature, or calculated from the published literature. Alternatively, the rules may be based on correlations generated from stored genomic profiles and previously known phenotypes.


In a preferred embodiment, the genetic variants will be SNPs. While SNPs occur at a single site, individuals who carry a particular SNP allele at one site often predictably carry specific SNP alleles at other sites. A correlation of SNPs and an allele predisposing an individual to disease or condition occurs through linkage disequilibrium, in which the non-random association of alleles at two or more loci occur more or less frequently in a population than would be expected from random formation through recombination.


Other genetic markers or variants, such as nucleotide repeats or insertions, may also be in linkage disequilibrium with genetic markers that have been shown to be associated with specific phenotypes. For example, a nucleotide insertion is correlated with a phenotype and a SNP is in linkage disequilibrium with the nucleotide insertion. A rule is made based on the correlation between the SNP and the phenotype. A rule based on the correlation between the nucleotide insertion and the phenotype may also be made. Either rules or both rules may be applied to a genomic profile, as the presence of one SNP may give a certain risk factor, the other may give another risk factor, and when combined may increase the risk.


Through linkage disequilibrium, a disease predisposing allele cosegregates with a particular allele of a SNP or a combination of particular alleles of SNPs. A particular combination of SNP alleles along a chromosome is termed a haplotype, and the DNA region in which they occur in combination can be referred to as a haplotype block. While a haplotype block can consist of one SNP, typically a haplotype block represents a contiguous series of 2 or more SNPs exhibiting low haplotype diversity across individuals and with generally low recombination frequencies. An identification of a haplotype can be made by identification of one or more SNPs that lie in a haplotype block. Thus, a SNP profile typically can be used to identify haplotype blocks without necessarily requiring identification of all SNPs in a given haplotype block.


Genotype correlations between SNP haplotype patterns and diseases, conditions or physical states are increasingly becoming known. For a given disease, the haplotype patterns of a group of people known to have the disease are compared to a group of people without the disease. By analyzing many individuals, frequencies of polymorphisms in a population can be determined, and in turn these frequencies or genotypes can be associated with a particular phenotype, such as a disease or a condition. Examples of known SNP-disease correlations include polymorphisms in Complement Factor H in age-related macular degeneration (Klein et al., Science: 308:385-389, (2005)) and a variant near the INSIG2 gene associated with obesity (Herbert et al., Science: 312:279-283 (2006)). Other known SNP correlations include polymorphisms in the 9p21 region that includes CDKN2A and B, such as) such as rs10757274, rs2383206, rs13333040, rs2383207, and rs10116277 correlated to myocardial infarction (Helgadottir et al., Science 316:1491-1493 (2007); McPherson et al., Science 316:1488-1491 (2007))


The SNPs may be functional or non-functional. For example, a functional SNP has an effect on a cellular function, thereby resulting in a phenotype, whereas a non-functional SNP is silent in function, but may be in linkage disequilibrium with a functional SNP. The SNPs may also be synonymous or non-synonymous. SNPs that are synonymous are SNPs in which the different forms lead to the same polypeptide sequence, and are non-functional SNPs. If the SNPs lead to different polypeptides, the SNP is non-synonymous and may or may not be functional. SNPs, or other genetic markers, used to identify haplotypes in a diplotype, which is 2 or more haplotypes, may also be used to correlate phenotypes associated with a diplotype. Information about an individual's haplotypes, diplotypes, and SNP profiles may be in the genomic profile of the individual.


In preferred embodiments, for a rule to be generated based on a genetic marker in linkage disequilibrium with another genetic marker that is correlated with a phenotype, the genetic marker may have a r2 or D′ score, scores commonly used in the art to determine linkage disequilibrium, of greater than 0.5. In preferred embodiments, the score is greater than 0.6, 0.7, 0.8, 0.90, 0.95 or 0.99. As a result, in the present disclosure, the genetic marker used to correlate a phenotype to an individual's genomic profile may be the same as the functional or published SNP correlated to a phenotype, or different. For example, using BC4, the test SNP and published SNP are the same, as are the test risk and nonrisk alleles are the same as the published risk and nonrisk alleles (FIGS. 4A and C). However, for BC5, CASP8 and its correlation to breast cancer, the test SNP is different from its functional or published SNP, as are the test risk and nonrisk alleles to the published risk and nonrisk alleles. The test and published alleles are oriented relative to the plus strand of the genome, and from these columns, it can be inferred the homozygous risk or nonrisk genotype, which may generate a rule to be applied to the genomic profile of individuals such as subscribers.


The test SNPs may be “DIRECT” or “TAG” SNPs (FIGS. 4E-G, FIG. 5). Direct SNPs are the test SNPs that are the same as the published or functional SNP, such as for BC4. Direct SNPs may also be used for FGFR2 correlation with breast cancer, using the SNP rs1073640 in Europeans and Asians, where the minor allele is A and the other allele is G (Easton et al., Nature 447:1087-1093 (2007)). Another published or functional SNP for FGFR2 correlation to breast cancer is rs1219648, also in Europeans and Asians (Hunter et al., Nat. Genet. 39:870-874 (2007)). Tag SNPs are where the test SNP is different from that of the functional or published SNP, as in for BC5. Tag SNPs may also be used for other genetic variants such as SNPs for CAMTA1 (rs4908449), 9p21 (rs10757274, rs2383206, rs13333040, rs2383207, rs10116277), COL1A1 (rs1800012), FVL (rs6025), HLA-DQA1 (rs4988889, rs2588331), eNOS (rs1799983), MTHFR (rs1801133), and APC (rs28933380).


Databases of SNPs are publicly available from, for example, the International HapMap Project (see www.hapmap.org, The International HapMap Consortium, Nature 426:789-796 (2003), and The International HapMap Consortium, Nature 437:1299-1320 (2005)), the Human Gene Mutation Database (HGMD) public database (see www.hgmd.org), and the Single Nucleotide Polymorphism database (dbSNP) (see www.ncbi.nlm.nih.gov/SNP/). These databases provide SNP haplotypes, or enable the determination of SNP haplotype patterns. Accordingly, these SNP databases enable examination of the genetic risk factors underlying a wide range of diseases and conditions, such as cancer, inflammatory diseases, cardiovascular diseases, neurodegenerative diseases, and infectious diseases. The diseases or conditions may be actionable, in which treatments and therapies currently exist. Treatments may include prophylactic treatments as well as treatments that ameliorate symptoms and conditions, including lifestyle changes.


Many other phenotypes such as physical traits, physiological traits, mental traits, emotional traits, ethnicity, ancestry, and age may also be examined. Physical traits may include height, hair color, eye color, body, or traits such as stamina, endurance, and agility. Mental traits may include intelligence, memory performance, or learning performance. Ethnicity and ancestry may include identification of ancestors or ethnicity, or where an individual's ancestors originated from. The age may be a determination of an individual's real age, or the age in which an individual's genetics places them in relation to the general population. For example, an individual's real age is 38 years of age, however their genetics may determine their memory capacity or physical well-being may be of the average 28 year old. Another age trait may be a projected longevity for an individual.


Other phenotypes may also include non-medical conditions, such as “fun” phenotypes. These phenotypes may include comparisons to well known individuals, such as foreign dignitaries, politicians, celebrities, inventors, athletes, musicians, artists, business people, and infamous individuals, such as convicts. Other “fun” phenotypes may include comparisons to other organisms, such as bacteria, insects, plants, or non-human animals. For example, an individual may be interested to see how their genomic profile compares to that of their pet dog, or to a former president.


At step 114, the rules are applied to the stored genomic profile to generate a phenotype profile of step 116. For example, information in FIG. 4, 5, or 6 may form the basis of rules, or tests, to apply to an individual's genomic profile. The rules may encompass the information on test SNP and alleles, and the effect estimates of FIG. 4, where the UNITS for effect estimate is the units of the effect estimate, such as OR, or odds-ratio (95% confidence interval) or mean. The effects estimate may be a genotypic risk (FIGS. 4C-G) in preferred embodiments, such as the risk for homozygotes (homoz or RR), risk heterozygotes (heteroz or RN), and nonrisk homozygotes (homoz or NN). In other embodiments, the effect estimate may be carrier risk, which is RR or RN vs NN. In yet other embodiments, the effect estimate may be based on the allele, an allelic risk such as R vs. N. There may also be two locus (FIG. 4J) or three locus (FIG. 4K) genotypic effect estimate (e.g. RRRR, RRNN, etc for the 9 possible genotype combinations for a two locus effect estimate). The test SNP frequency in the public HapMap is also noted in FIGS. 4H and I.


The estimated risk for a condition may be based on the SNPs as listed in US Patent Application Publication No. 20080131887 and PCT Publication No. WO2008/067551. In some embodiments, the risk for a condition may be based on at least one SNP. For example, assessment of an individual's risk for Alzheimers (AD), colorectal cancer (CRC), osteoarthritis (OA) or exfoliation glaucoma (XFG), may be based on 1 SNP (for example, rs4420638 for AD, rs6983267 for CRC, rs4911178 for OA and rs2165241 for XFG). For other conditions, such as obesity (BMIOB), Graves' disease (GD), or hemochromatosis (HEM), an individual's estimated risk may be based on at least 1 or 2 SNPs (for example, rs9939609 and/or rs9291171 for BMIOB; DRB1*0301 DQA1*0501 and/or rs3087243 for GD; rs1800562 and/or rs129128 for HEM). For conditions such as, but not limited to, myocardial infarction (MI), multiple sclerosis (MS), or psoriasis (PS), 1, 2, or 3 SNPs may be used to assess an individual's risk for the condition (for example, rs1866389, rs1333049, and/or rs6922269 for MI; rs6897932, rs12722489, and/or DRB1*1501 for MS; rs6859018, rs11209026, and/or HLAC*0602 for PS). For estimating an individual's risk of restless legs syndrome (RLS) or celiac disease (Ce1D), 1, 2, 3, or 4 SNPs (for example, rs6904723, rs2300478, rs1026732, and/or rs9296249 for RLS; rs6840978, rs11571315, rs2187668, and/or DQA1*0301 DQB1*0302 for Ce1D). For prostate cancer (PC) or lupus (SLE), 1, 2, 3, 4, or 5 SNPs may be used to estimate an individual's risk for PC or SLE (for example, rs4242384, rs6983267, rs16901979, rs17765344, and/or rs4430796 for PC; rs12531711, rs10954213, rs2004640, DRB1*0301, and/or DRB1*1501 for SLE). For estimating an individual's lifetime risk of macular degeneration (AMD) or rheumatoid arthritis (RA), 1, 2, 3, 4, 5, or 6 SNPs, may be used (for example, rs10737680, rs10490924, rs541862, rs2230199, rs1061170, and/or rs9332739 for AMD; rs6679677, rs11203367, rs6457617, DRB*0101, DRB1*0401, and/or DRB1*0404 for RA). For estimating an individual's lifetime risk of breast cancer (BC), 1, 2, 3, 4, 5, 6 or 7 SNPs may be used (for example, rs3803662, rs2981582, rs4700485, rs3817198, rs17468277, rs6721996, and/or rs3803662). For estimating an individual's lifetime risk of Crohn's disease (CD) or Type 2 diabetes (T2D), 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, or 11 SNPs may be used (for example, rs2066845, rs5743293, rs10883365, rs17234657, rs10210302, rs9858542, rs11805303, rs1000113, rs17221417, rs2542151, and/or rs10761659 for CD; rs13266634, rs4506565, rs10012946, rs7756992, rs10811661, rs12288738, rs8050136, rs111875, rs4402960, rs5215, and/or rs1801282 for T2D). In some embodiments, the SNPs used as a basis for determining risk may be in linkage disequilibrium with the SNPs as mentioned above, or other SNPs, such as in US Patent Publication No. 20080131887 and PCT Publication No. WO2008/067551.


The phenotype profile of an individual may comprise a number of phenotypes. In particular, the assessment of a patient's risk of disease or other conditions such as likely drug response including metabolism, efficacy and/or safety, by the methods of the present disclosure allows for prognostic or diagnostic analysis of susceptibility to multiple, unrelated diseases and conditions, whether in symptomatic, presymptomatic or asymptomatic individuals, including carriers of one or more disease/condition predisposing alleles. Accordingly, these methods provide for general assessment of an individual's susceptibility to disease or condition without any preconceived notion of testing for a specific disease or condition. For example, the methods of the present disclosure allow for assessment of an individual's susceptibility to any of the several conditions listed in Tables 1, FIG. 4, 5, or 6, based on the individual's genomic profile. The assessment preferably provides information for 2 or more of these conditions, and more preferably, 3, 4, 5, 10, 20, 50, 100 or even more of these conditions. In preferred embodiments, the phenotype profile results from the application of at least 20 rules to the genomic profile of an individual. In other embodiments, at least 50 rules are applied to the genomic profile of an individual. A single rule for a phenotype may be applied for monogenic phenotypes. More than one rule may also be applied for a single phenotype, such as a multigenic phenotype or a monogenic phenotype wherein multiple genetic variants within a single gene affects the probability of having the phenotype.


Following an initial screening of an individual patient's genomic profile, updates of an individual's genotype correlations are made (or are available) through comparisons to additional nucleotide variants, such as SNPs, when such additional nucleotide variants become known. For example, step 110 may be performed periodically, for example, daily, weekly, or monthly by one or more people of ordinary skill in the field of genetics, who scan scientific literature for new genotype correlations. The new genotype correlations may then be further validated by a committee of one or more experts in the field. Step 112 may then also be periodically updated with new rules based on the new validated correlations.


The new rule may encompass a genotype or phenotype without an existing rule. For example, a genotype not correlated with any phenotype is discovered to correlate with a new or existing phenotype. A new rule may also be for a correlation between a phenotype for which no genotype has previously been correlated to. New rules may also be determined for genotypes and phenotypes that have existing rules. For example, a rule based on the correlation between genotype A and phenotype A exists. New research reveals genotype B correlates with phenotype A, and a new rule based on this correlation is made. Another example is phenotype B is discovered to be associated with genotype A, and thus a new rule may be made.


Rules may also be made on discoveries based on known correlations but not initially identified in published scientific literature. For example, it may be reported genotype C is correlated with phenotype C. Another publication reports genotype D is correlated with phenotype D. Phenotype C and D are related symptoms, for example phenotype C may be shortness of breath, and phenotype D is small lung capacity. A correlation between genotype C and phenotype D, or genotype D with phenotype C, may be discovered and validated through statistical means with existing stored genomic profiles of individuals with genotypes C and D, and phenotypes C and D, or by further research. A new rule may then be generated based on the newly discovered and validated correlation. In another embodiment, stored genomic profiles of a number of individuals with a specific or related phenotype may be studied to determine a genotype common to the individuals, and a correlation may be determined. A new rule may be generated based on this correlation.


Rules may also be made to modify existing rules. For example, correlations between genotypes and phenotypes may be partly determined by a known individual characteristic, such as ethnicity, ancestry, geography, gender, age, family history, or any other known phenotypes of the individual. Rules based on these known individual characteristics may be made and incorporated into an existing rule, to provide a modified rule. The choice of modified rule to be applied will be dependent on the specific individual factor of an individual. For example, a rule may be based on the probability an individual who has phenotype E is 35% when the individual has genotype E. However, if an individual is of a particular ethnicity, the probability is 5%. A new rule may be generated based on this result and applied to individuals with that particular ethnicity. Alternatively, the existing rule with a determination of 35% may be applied, and then another rule based on ethnicity for that phenotype is applied. The rules based on known individual characteristics may be determined from scientific literature or determined based on studies of stored genomic profiles. New rules may be added and applied to genomic profiles in step 114, as the new rules are developed, or they may be applied periodically, such as at least once a year.


Information of an individual's risk of disease can also be expanded as technology advances allow for finer resolution SNP genomic profiles. As indicated above, an initial SNP genomic profile readily can be generated using microarray technology for scanning of 500,000 SNPs. Given the nature of haplotype blocks, this number allows for a representative profile of all SNPs in an individual's genome. Nonetheless, there are approximately 10 million SNPs estimated to occur commonly in the human genome (the International HapMap Project; www.hapmap.org). As technological advances allow for practical, cost-efficient resolution of SNPs at a finer level of detail, such as microarrays of 1,000,000, 1,500,000, 2,000,000, 3,000,000, or more SNPs, or whole genomic sequencing, more detailed SNP genomic profiles can be generated. Likewise, cost-efficient analysis of finer SNP genomic profiles and updates to the master database of SNP-disease correlations will be enabled by advances in computational analytical methodology.


After generation of phenotype profile at step 116, a subscriber or their health care manager may access their genomic or phenotype profiles via an on-line portal or website as in step 118. Reports containing phenotype profiles and other information related to the phenotype and genomic profiles may also be provided to the subscriber or their health care manager, as in steps 120 and 122. The reports may be printed, saved on the subscriber's computer, or viewed on-line.


A sample on-line report is shown in FIG. 7. The subscriber may choose to display a single phenotype, or more than one phenotype. The subscriber may also have different viewing options, for example, as shown in FIG. 7, a “Quick View” option. The phenotype may be a medical condition and different treatments and symptoms in the quick report may link to other web pages that contain further information about the treatment. For example, by clicking on a drug, it will lead to website that contains information about dosages, costs, side effects, and effectiveness. It may also compare the drug to other treatments. The website may also contain a link leading to the drug manufacturer's website. Another link may provide an option for the subscriber to have a pharmacogenomic profile generated, which would include information such as their likely response to the drug based on their genomic profile. Links to alternatives to the drug may also be provided, such as preventative action such as fitness and weight loss, and links to diet supplements, diet plans, and to nearby health clubs, health clinics, health and wellness providers, day spas and the like may also be provided. Educational and informational videos, summaries of available treatments, possible remedies, and general recommendations may also be provided.


The on-line report may also provide links to schedule in-person physician or genetic counseling appointments or to access an on-line genetic counselor or physician, providing the opportunity for a subscriber to ask for more information regarding their phenotype profile. Links to on-line genetic counseling and physician questions may also be provided on the on-line report.


Reports may also be viewed in other formats such as a comprehensive view for a single phenotype, wherein more detail for each category is provided. For example, there may be more detailed statistics about the likelihood of the subscriber developing the phenotype, more information about the typical symptoms or phenotypes, such as sample symptoms for a medical condition, or the range of a physical non-medical condition such as height, or more information about the gene and genetic variant, such as the population incidence, for example in the world, or in different countries, or in different age ranges or genders. In another embodiment, the report may be of a “fun” phenotype, such as the similarity of an individual's genomic profile to that of a famous individual, such as Albert Einstein. The report may display a percentage similarity between the individual's genomic profile to that of Einstein's, and may further display a predicted IQ of Einstein and that of the individual's. Further information may include how the genomic profile of the general population and their IQ compares to that of the individual's and Einstein's.


In another embodiment, the report may display all phenotypes that have been correlated to the subscriber's genomic profile. In other embodiments, the report may display only the phenotypes that are positively correlated with an individual's genomic profile. In other formats, the individual may choose to display certain subgroups of phenotypes, such as only medical phenotypes, or only actionable medical phenotypes. For example, actionable phenotypes and their correlated genotypes, may include Crohn's disease (correlated with IL23R and CARD 15), Type 1 diabetes (correlated with HLA-DR/DQ), lupus (correlated HLA-DRB1), psoriasis (HLA-C), multiple sclerosis (HLA-DQA1), Graves disease (HLA-DRB1), rheumatoid arthritis (HLA-DRB1), Type 2 diabetes (TCF7L2), breast cancer (BRCA2), colon cancer (APC), episodic memory (KIBRA), and osteoporosis (COL1A1). The individual may also choose to display subcategories of phenotypes in their report, such as only inflammatory diseases for medical conditions, or only physical traits for non-medical conditions.


Information submitted by and conveyed to an individual may be secure and confidential, and access to such information may be controlled by the individual. Information derived from the complex genomic profile may be supplied to the individual as regulatory agency approved, understandable, medically relevant and/or high impact data. Information may also be of general interest, and not medically relevant. Information can be securely conveyed to the individual by several means including, but not restricted to, a portal interface and/or mailing. More preferably, information is securely (if so elected by the individual) provided to the individual by a portal interface, to which the individual has secure and confidential access. Such an interface is preferably provided by on-line, internet website access, or in the alternative, telephone or other means that allow private, secure, and readily available access. The genomic profiles, phenotype profiles, and reports are provided to an individual or their health care manager by transmission of the data over a network.


Accordingly, FIG. 8 is a block diagram showing a representative example logic device through which a phenotype profile and report may be generated. FIG. 8 shows a computer system (or digital device) 800 to receive and store genomic profiles, analyze genotype correlations, generate rules based on the analysis of genotype correlations, apply the rules to the genomic profiles, and produce a phenotype profile and report. The computer system 800 may be understood as a logical apparatus that can read instructions from media 811 and/or network port 805, which can optionally be connected to server 809 having fixed media 812. The system shown in FIG. 8 includes CPU 801, disk drives 803, optional input devices such as keyboard 815 and/or mouse 816 and optional monitor 807. Data communication can be achieved through the indicated communication medium to a server 809 at a local or a remote location. The communication medium can include any means of transmitting and/or receiving data. For example, the communication medium can be a network connection, a wireless connection or an internet connection. Such a connection can provide for communication over the World Wide Web. It is envisioned that data relating to the present disclosure can be transmitted over such networks or connections for reception and/or review by a party 822. The receiving party 822 can be but is not limited to an individual, a subscriber, a health care provider or a health care manager. In one embodiment, a computer-readable medium includes a medium suitable for transmission of a result of an analysis of a biological sample or a genotype correlation. The medium can include a result regarding a phenotype profile of an individual subject, wherein such a result is derived using the methods described herein.


A personal portal will preferably serve as the primary interface with an individual for receiving and evaluating genomic data. A portal will enable individuals to track the progress of their sample from collection through testing and results. Through portal access, individuals are introduced to relative risks for common genetic disorders based on their genomic profile. The subscriber may choose which rules to apply to their genomic profile through the portal.


In one embodiment, one or more web pages will have a list of phenotypes and next to each phenotype a box in which a subscriber may select to include in their phenotype profile. The phenotypes may be linked to information on the phenotype, to help the subscriber make an informed choice about the phenotype they want included in their phenotype profile. The webpage may also have phenotypes organized by disease groups, for example as actionable diseases or not. For example, a subscriber may choose actionable phenotypes only, such as HLA-DQA1 and celiac disease. The subscriber may also choose to display pre or post symptomatic treatments for the phenotypes. For example, the individual may choose actionable phenotypes with pre-symptomatic treatments (outside of increased screening), for celiac disease, a pre-symptomatic treatment of gluten free diet. Another example may be Alzheimer's, the pre-symptomatic treatment of statins, exercise, vitamins, and mental activity. Thrombosis is another example, with a pre-symptomatic treatment of avoid oral contraceptives and avoid sitting still for long periods of time. An example of a phenotype with an approved post symptomatic treatment is wet AMD, correlated with CFH, wherein individuals may obtain laser treatment for their condition.


The phenotypes may also be organized by type or class of disease or conditions, for example neurological, cardiovascular, endocrine, immunological, and so forth. Phenotypes may also be grouped as medical and non-medical phenotypes. Other groupings of phenotypes on the webpage may be by physical traits, physiological traits, mental traits, or emotional traits. The webpage may further provide a section in which a group of phenotypes are chosen by selection of one box. For example, a selection for all phenotypes, only medically relevant phenotypes, only non-medically relevant phenotypes, only actionable phenotypes, only non-actionable phenotypes, different disease group, or “fun” phenotypes. “Fun” phenotypes may include comparisons to celebrities or other famous individuals, or to other animals or even other organisms. A list of genomic profiles available for comparison may also be provided on the webpage for selection by the subscriber to compare to the subscriber's genomic profile.


The on-line portal may also provide a search engine, to help the subscriber navigate the portal, search for a specific phenotype, or search for specific terms or information revealed by their phenotype profile or report. Links to access partner services and product offerings may also be provided by the portal. Additional links to support groups, message boards, and chat rooms for individuals with a common or similar phenotype may also be provided. The on-line portal may also provide links to other sites with more information on the phenotypes in a subscriber's phenotype profile. The on-line portal may also provide a service to allow subscribers to share their phenotype profile and reports with friends, families, or health care managers. Subscribers may choose which phenotypes to show in the phenotype profile they want shared with their friends, families, or health care managers.


The phenotype profiles and reports provide a personalized genotype correlation to an individual. The genotype correlations provided to an individual can be used in determining personal health care and lifestyle choices. If a strong correlation is found between a genetic variant and a disease for which treatment is available, detection of the genetic variant may assist in deciding to begin treatment of the disease and/or monitoring of the individual. In the case where a statistically significant correlation exists but is not regarded as a strong correlation, an individual can review the information with a personal physician and decide an appropriate, beneficial course of action. Potential courses of action that could be beneficial to an individual in view of a particular genotype correlation include administration of therapeutic treatment, monitoring for potential need of treatment or effects of treatment, or making life-style changes in diet, exercise, and other personal habits/activities. For example, an actionable phenotype such as celiac disease may have a pre-symptomatic treatment of a gluten-free diet. Likewise, genotype correlation information could be applied through pharmacogenomics to predict the likely response an individual would have to treatment with a particular drug or regimen of drugs, such as the likely efficacy or safety of a particular drug treatment.


Subscribers may choose to provide the genomic and phenotype profiles to their health care managers, such as a physician or genetic counselor. The genomic and phenotype profiles may be directly accessed by the healthcare manager, by the subscriber printing out a copy to be given to the healthcare manager, or have it directly sent to the healthcare manager through the on-line portal, such as through a link on the on-line report.


Delivery of this pertinent information will empower patients to act in concert with their physician. In particular, discussions between patients and their physicians can be empowered through an individual's portal and links to medical information, and the ability to tie patient's genomic information into their medical records. Medical information may include prevention and wellness information. The information provided to the individual patient by the present disclosure will enable patients to make informed choices for their health care. In this manner, patients will be able to make choices that may help them avoid and/or delay diseases that their individual genomic profile (inherited DNA) makes more likely. In addition, patients will be able to employ a treatment regime that personally fits their specific medical needs. Individuals also will have the ability to access their genotype data should they develop an illness and need this information to help their physician form a therapeutic strategy.


Genotype correlation information could also be used in cooperation with genetic counseling to advise couples considering reproduction, and potential genetic concerns to the mother, father and/or child. Genetic counselors may provide information and support to subscribers with phenotype profiles that display an increased risk for specific conditions or diseases. They may interpret information about the disorder, analyze inheritance patterns and risks of recurrence, and review available options with the subscriber. Genetic counselors may also provide supportive counseling refer subscribers to community or state support services. Genetic counseling may be included with specific subscription plans. In some embodiments, genetic counseling may be scheduled within 24 hours of request and available during of hours such as evenings, Saturdays, Sundays, and/or holidays.


An individual's portal will also facilitate delivery of additional information beyond an initial screening. Individuals will be informed about new scientific discoveries that relate to their personal genetic profile, such as information on new treatments or prevention strategies for their current or potential conditions. The new discoveries may also be delivered to their healthcare managers. In preferred embodiments, the subscribers, or their healthcare providers are informed of new genotype correlations and new research about the phenotypes in the subscriber's phenotype profiles, by e-mail. In other embodiments, e-mails of “fun” phenotypes are sent to subscribers, for example, an e-mail may inform them that their genomic profile is 77% identical to that of Abraham Lincoln and that further information is available via an on-line portal.


The present disclosure also provides a system of computer code for generating new rules, modifying rules, combining rules, periodically updating the rule set with new rules, maintaining a database of genomic profile securely, applying the rules to the genomic profiles to determine phenotype profiles, and for generating reports. Computer code for notifying subscribers of new or revised correlations new or revised rules, and new or revised reports, for example with new prevention and wellness information, information about new therapies in development, or new treatments available.


Business Method

The present disclosure provides a business method of assessing an individual's genotype correlations based on comparison of the patient's genome profile against a clinically-derived database of established, medically relevant nucleotide variants. The present disclosure further provides a business method for using the stored genomic profile of the individual for assessing new correlations that were not initially known, to generate updated phenotype profiles for an individual, without the requirement of the individual submitting another biological sample. A flow chart illustrating the business method is in FIG. 9.


A revenue stream for the subject business method is generated in part at step 101, when an individual initially requests and purchases a personalized genomic profile for genotype correlations for a multitude of common human diseases, conditions, and physical states. A request and purchase can be made through any number of sources, including but not limited to, an on-line web portal, an on-line health service, and an individual's personal physician or similar source of personal medical attention. In an alternative embodiment, the genomic profile may be provided free, and the revenue stream is generated at a later step, such as step 103.


A subscriber, or customer, makes a request for purchase of a phenotype profile. In response to a request and purchase, a customer is provided a collection kit for a biological sample used for genetic sample isolation at step 103. When a request is made on-line, by telephone, or other source in which a collection kit is not readily physically available to the customer, a collection kit is provided by expedited delivery, such as courier service that provides same-day or overnight delivery. Included in the collection kit is a container for a sample, as well as packaging materials for expedited delivery of the sample to a laboratory for genomic profile generation. The kit may also include instructions for sending the sample to the sample processing facility, or laboratory, and instructions for accessing their genomic profile and phenotype profile, which may occur through an on-line portal.


As detailed above, genomic DNA can be obtained from any of a number of types of biological samples. Preferably, genomic DNA is isolated from saliva, using a commercially available collection kit such as that available from DNA Genotek. Use of saliva and such a kit allows for a non-invasive sample collection, as the customer conveniently provides a saliva sample in a container from a collection kit and then seals the container. In addition, a saliva sample can be stored and shipped at room temperature.


After depositing a biological sample into a collection or specimen container, a customer will deliver the sample to a laboratory for processing at step 105. Typically, the customer may use packaging materials provided in the collection kit to deliver/send the sample to a laboratory by expedited delivery, such as same-day or overnight courier service.


The laboratory that processes the sample and generates the genomic profile may adhere to appropriate governmental agency guidelines and requirements. For example, in the United States, a processing laboratory may be regulated by one or more federal agencies such as the Food and Drug Administration (FDA) or the Centers for Medicare and Medicaid Services (CMS), and/or one or more state agencies. In the United States, a clinical laboratory may be accredited or approved under the Clinical Laboratory Improvement Amendments of 1988 (CLIA).


At step 107, the laboratory processes the sample as previously described to isolate the genetic sample of DNA or RNA. Analysis of the isolated genetic sample and generation of a genomic profile is then performed at step 109. Preferably, a genomic SNP profile is generated. As described above, several methodologies may be used to generate a SNP profile. Preferably, a high density array, such as the commercially available platforms from Affymetrix or Illumina, is used for SNP identification and profile generation. For example, a SNP profile may be generated using an Affymetrix GeneChip assay, as described above in more detail. As technology evolves, there may be other technology vendors who can generate high density SNP profiles. In another embodiment, a genomic profile for a subscriber will be the genomic sequence of the subscriber.


Following generation of an individual's genomic profile, the genotype data is preferably encrypted, imported at step 111, and deposited into a secure database or vault at step 113, where the information is stored for future reference. The genomic profile and related information may be confidential, with access to this proprietary information and the genomic profile limited as directed by the individual and/or his or her personal physician. Others, such as family and the genetic counselor of the individual may also be permitted access by the subscriber.


The database or vault may be located on-site with the processing laboratory. Alternatively, the database may be located at a separate location. In this scenario, the genomic profile data generated by the processing lab can be imported at step 111 to a separate facility that contains the database.


After an individual's genomic profile is generated, the individual's genetic variations are then compared against a clinically-derived database of established, medically relevant genetic variants in step 115. Alternatively, the genotype correlations may not be medically relevant but still incorporated into the database of genotype correlations, for example, physical traits such as eye color, or “fun” phenotypes such as genomic profile similarity to a celebrity.


The medically relevant SNPs may have been established through the scientific literature and related sources. The non-SNP genetic variants may also be established to be correlated with phenotypes. Generally, the correlation of SNPs to a given disease is established by comparing the haplotype patterns of a group of people known to have the disease to a group of people without the disease. By analyzing many individuals, frequencies of polymorphisms in a population can be determined, and in turn these genotype frequencies can be associated with a particular phenotype, such as a disease or a condition. Alternatively, the phenotype may be a non-medical condition.


The relevant SNPs and non-SNP genetic variants may also be determined through analysis of the stored genomic profiles of individuals rather than determined by available published literature. Individuals with stored genomic profiles may disclose phenotypes that have previously been determined. Analysis of the genotypes and disclosed phenotypes of the individuals may be compared to those without the phenotypes to determine a correlation that may then be applied to other genomic profiles. Individuals that have their genomic profiles determined may fill out questionnaires about phenotypes that have previously been determined. Questionnaires may contain questions about medical and non-medical conditions, such as diseases previously diagnosed, family history of medical conditions, lifestyle, physical traits, mental traits, age, social life, environment and the like.


In one embodiment, an individual may have their genomic profile determined free of charge if they fill out a questionnaire. In some embodiments, the questionnaires are to be filled out periodically by the individuals in order to have free access to their phenotype profile and reports. In other embodiments, the individuals that fill out the questionnaires may be entitled to a subscription upgrade, such that they have more access than their previous subscription level, or they may purchase or renew a subscription at a reduced cost.


All information deposited in the database of medically relevant genetic variants at step 121 is first approved by a research/clinical advisory board for scientific accuracy and importance, coupled with review and oversight by an appropriate governmental agency if warranted at step 119. For example, in the United States, the FDA may provide oversight through approval of algorithms used for validation of genetic variant (typically SNP, transcript level, or mutation) correlative data. At step 123, scientific literature and other relevant sources are monitored for additional genetic variant-disease or condition correlations, and following validation of their accuracy and importance, along with governmental agency review and approval, these additional genotype correlations are added to the master database at step 125.


The database of approved, validated medically-relevant genetic variants, coupled with a genome-wide individual profile, will advantageously allow genetic risk-assessment to be performed for a large number of diseases or conditions. Following compilation of an individual's genomic profile, individual genotype correlations can be determined through comparison of the individual's nucleotide (genetic) variants or markers with a database of human nucleotide variants that have been correlated to a particular phenotype, such as a disease, condition, or physical state. Through comparison of an individual's genomic profile to the master database of genotype correlations, the individual can be informed whether they are found to be positive or negative for a genetic risk factor, and to what degree. An individual will receive relative risk and/or predisposition data on a wide range of scientifically validated disease states (e.g., Alzheimer's, cardiovascular disease, blood clotting). For example, genotype correlations in Table 1 may be included. In addition, SNP disease correlations in the database may include, but are not limited to, those correlations shown in FIG. 4. Other correlations from FIGS. 5 and 6 may also be included. The subject business method therefore provides analysis of risk to a multitude of diseases and conditions without any preconceived notion of what those diseases and conditions might entail.


In other embodiments, the genotype correlations that are coupled to the genome wide individual profile are non-medically relevant phenotypes, such as “fun” phenotypes or physical traits such as hair color. In preferred embodiments, a rule or rule set is applied to the genomic profile or SNP profile of an individual, as described above. Application of the rules to a genomic profile generates a phenotype profile for the individual.


Accordingly, the master database of human genotype correlations is expanded with additional genotype correlations as new correlations become discovered and validated. An update can be made by accessing pertinent information from the individual's genomic profile stored in a database as desired or appropriate. For example, a new genotype correlation that becomes known may be based on a particular gene variant. Determination of whether an individual may be susceptible to that new genotype correlation can then be made by retrieving and comparing just that gene portion of the individual's entire genomic profile.


The results of the genomic query preferably are analyzed and interpreted so as to be presented to the individual in an understandable format. At step 117, the results of an initial screening are then provided to the patient in a secure, confidential form, either by mailing or through an on-line portal interface, as detailed above.


The report may contain the phenotype profile as well as genomic information about the phenotypes in the phenotype profile, for example basic genetics about the genes involved or the statistics of the genetic variants in different populations. Other information based on the phenotype profile that may be included in the report are prevention strategies, wellness information, therapies, symptom awareness, early detection schemes, intervention schemes, and refined identification and sub-classification of the phenotypes. Following an initial screening of an individual's genomic profile, controlled, moderated updates are or can be made.


Updates of an individual's genomic profile are made or are available in conjunction with updates to the master database as new genotype correlations emerge and are both validated and approved. New rules based on the new genotype correlations may be applied to the initial genomic profile to provide updated phenotype profiles. An updated genotype correlation profile can be generated by comparing the relevant portion of the individual's genomic profile to a new genotype correlation at step 127. For example, if a new genotype correlation is found based on variation in a particular gene, then that gene portion of the individual's genomic profile can be analyzed for the new genotype correlation. In such a case, one or more new rules may be applied to generate an updated phenotype profile, rather than an entire rule set with rules that had already been applied. The results of the individual's updated genotype correlations are provided in a secure manner at step 129.


Initial and updated phenotype profiles may be a service provided to subscribers or customers. Varying levels of subscriptions to genomic profile analysis and combinations thereof can be provided. Likewise, subscription levels can vary to provide individuals choices of the amount of service they wish to receive with their genotype correlations. Thus, the level of service provided would vary with the level of service subscription purchased by the individual.


An entry level subscription for a subscriber may include a genomic profile and an initial phenotype profile. This may be a basic subscription level. Within the basic subscription level may be varying levels of service. For example, a particular subscription level could provide references for genetic counseling, physicians with particular expertise in treating or preventing a particular disease, and other service options. Genetic counseling may be obtained on-line or by telephone. In another embodiment, the price of the subscription may depend on the number of phenotypes an individual chooses for their phenotype profile. Another option may be whether the subscriber chooses to access on-line genetic counseling.


In another scenario, a subscription could provide for an initial genome-wide, genotype correlation, with maintenance of the individual's genomic profile in a database; such database may be secure if so elected by the individual. Following this initial analysis, subsequent analyses and additional results could be made upon request and additional payment by the individual. This may be a premium level of subscription.


In one embodiment of the subject business method, updates of an individual's risks are performed and corresponding information made available to individuals on a subscription basis. The updates may be available to subscribers who purchase the premium level of subscription. Subscription to genotype correlation analysis can provide updates with a particular category or subset of new genotype correlations according to an individual's preferences. For example, an individual might only wish to learn of genotype correlations for which there is a known course of treatment or prevention. To aid an individual in deciding whether to have an additional analysis performed, the individual can be provided with information regarding additional genotype correlations that have become available. Such information can be conveniently mailed or e-mailed to a subscriber.


Within the premium subscription, there may be further levels of service, such as those mentioned in the basic subscription. Other subscription models may be provided within the premium level. For example, the highest level may provide a subscriber to unlimited updates and reports. The subscriber's profile may be updated as new correlations and rules are determined. At this level, subscribers may also permit access to unlimited number of individuals, such as family members and health care managers. The subscribers may also have unlimited access to on-line genetic counselors and physicians.


The next level of subscription within the premium level may provide more limited aspects, for example a limited number of updates. The subscriber may have a limited number of updates for their genomic profile within a subscription period, for example, 4 times a year. In another subscription level, the subscriber may have their stored genomic profile updated once a week, once a month, or once a year. In another embodiment, the subscriber may only have a limited number of phenotypes they may choose to update their genomic profile against.


A personal portal will also conveniently allow an individual to maintain a subscription to risk or correlation updates and information updates or alternatively, make requests for updated risk assessment and information. As described above, varying subscription levels could be provided to allow individuals choices of various levels of genotype correlation results and updates and may different subscription levels may be chosen by the subscriber via their personal portal.


Any of these subscription options will contribute to the revenue stream for the subject business method. The revenue stream for the subject business method will also be added by the addition of new customers and subscribers, wherein the new genomic profiles are added to the database.









TABLE 1







Representative genes having genetic variants


correlated with a phenotype.








Gene
Phenotype





A2M
Alzheimer's Disease


ABCA1
cholesterol, HDL


ABCB1
HIV


ABCB1
epilepsy


ABCB1
kidney transplant complications


ABCB1
digoxin, serum concentration


ABCB1
Crohn's disease; ulcerative colitis


ABCB1
Parkinson's disease


ABCC8
Type 2 diabetes


ABCC8
diabetes, type 2


ABO
myocardial infarct


ACADM
medium-chain acyl-CoA dehydrogenase deficiency


ACDC
Type 2 diabetes


ACE
Type 2 diabetes


ACE
hypertension


ACE
Alzheimer's Disease


ACE
myocardial infarction


ACE
cardiovascular


ACE
left ventricular hypertrophy


ACE
coronary artery disease


ACE
atherosclerosis, coronary


ACE
retinopathy, diabetic


ACE
systemic lupus erythematosus


ACE
blood pressure, arterial


ACE
erectile dysfunction


ACE
Lupus


ACE
polycystic kidney disease


ACE
stroke


ACP1
diabetes, type 1


ACSM1 (LIP)c
cholesterol levels


ADAM33
asthma


ADD1
hypertension


ADD1
blood pressure, arterial


ADH1B
alcohol abuse


ADH1C
alcohol abuse


ADIPOQ
diabetes, type 2


ADIPOQ
obesity


ADORA2A
panic disorder


ADRB1
hypertension


ADRB1
heart failure


ADRB2
asthma


ADRB2
hypertension


ADRB2
obesity


ADRB2
blood pressure, arterial


ADRB2
Type 2 Diabetes


ADRB3
obesity


ADRB3
Type 2 Diabetes


ADRB3
hypertension


AGT
hypertension


AGT
Type 2 diabetes


AGT
Essential Hypertension


AGT
myocardial infarction


AGTR1
hypertension


AGTR2
hypertension


AHR
breast cancer


ALAD
lead toxicity


ALDH2
alcoholism


ALDH2
alcohol abuse


ALDH2
colorectal cancer


ALDRL2
Type 2 diabetes


ALOX5
asthma


ALOX5AP
asthma


APBB1
Alzheimer's Disease


APC
colorectal cancer


APEX1
lung cancer


APOA1
atherosclerosis, coronary


APOA1
cholesterol, HDL


APOA1
coronary artery disease


APOA1
Type 2 diabetes


APOA4
Type 2 diabetes


APOA5
triglycerides


APOA5
atherosclerosis, coronary


APOB
hypercholesterolemia


APOB
obesity


APOB
cardiovascular


APOB
coronary artery disease


APOB
coronary heart disease


APOB
Type 2 diabetes


APOC1
Alzheimer's Disease


APOC3
triglycerides


APOC3
Type 2 Diabetes


APOE
Alzheimer's Disease


APOE
Type 2 diabetes


APOE
multiple sclerosis


APOE
atherosclerosis, coronary


APOE
Parkinson's disease


APOE
coronary heart disease


APOE
myocardial infarction


APOE
stroke


APOE
Alzheimer's disease


APOE
coronary artery disease


APP
Alzheimer's Disease


AR
prostate cancer


AR
breast cancer


ATM
breast cancer


ATP7B
Wilson disease


ATXN8OS
spinocerebellar ataxia


BACE1
Alzheimer's Disease


BCHE
Alzheimer's Disease


BDKRB2
hypertension


BDNF
Alzheimer's Disease


BDNF
bipolar disorder


BDNF
Parkinson's disease


BDNF
schizophrenia


BDNF
memory


BGLAP
bone density


BRAF
thyroid cancer


BRCA1
breast cancer


BRCA1
breast cancer; ovarian cancer


BRCA1
ovarian cancer


BRCA2
breast cancer


BRCA2
breast cancer; ovarian cancer


BRCA2
ovarian cancer


BRIP1
breast cancer


C4A
systemic lupus erythematosus


CALCR
bone density


CAMTA1
episodic memory


CAPN10
diabetes, type 2


CAPN10
Type 2 diabetes


CAPN3
muscular dystrophy


CARD15
Crohn's disease


CARD15
Crohn's disease; ulcerative colitis


CARD15
Inflammatory Bowel Disease


CART
obesity


CASR
bone density


CCKAR
schizophrenia


CCL2
systemic lupus erythematosus


CCL5
HIV


CCL5
asthma


CCND1
colorectal cancer


CCR2
HIV


CCR2
HIV infection


CCR2
hepatitis C


CCR2
myocardial infarct


CCR3
Asthma


CCR5
HIV


CCR5
HIV infection


CCR5
hepatitis C


CCR5
asthma


CCR5
multiple sclerosis


CD14
atopy


CD14
asthma


CD14
Crohn's disease


CD14
Crohn's disease; ulcerative colitis


CD14
periodontitis


CD14
Total IgE


CDH1
prostate cancer


CDH1
colorectal cancer


CDKN2A
melanoma


CDSN
psoriasis


CEBPA
leukemia, myeloid


CETP
atherosclerosis, coronary


CETP
coronary heart disease


CETP
hypercholesterolemia


CFH
macular degeneration


CFTR
cystic fibrosis


CFTR
pancreatitis


CFTR
Cystic Fibrosis


CHAT
Alzheimer's Disease


CHEK2
breast cancer


CHRNA7
schizophrenia


CMA1
atopic dermatitis


CNR1
schizophrenia


COL1A1
bone density


COL1A1
osteoporosis


COL1A2
bone density


COL2A1
Osteoarthritis


COMT
schizophrenia


COMT
breast cancer


COMT
Parkinson's disease


COMT
bipolar disorder


COMT
obsessive compulsive disorder


COMT
alcoholism


CR1
systemic lupus erythematosus


CRP
C-reactive protein


CST3
Alzheimer's Disease


CTLA4
Type 1 diabetes


CTLA4
Graves' disease


CTLA4
multiple sclerosis


CTLA4
rheumatoid arthritis


CTLA4
systemic lupus erythematosus


CTLA4
lupus erythematosus


CTLA4
celiac disease


CTSD
Alzheimer's Disease


CX3CR1
HIV


CXCL12
HIV


CXCL12
HIV infection


CYBA
atherosclerosis, coronary


CYBA
hypertension


CYP11B2
hypertension


CYP11B2
left ventricular hypertrophy


CYP17A1
breast cancer


CYP17A1
prostate cancer


CYP17A1
endometriosis


CYP17A1
endometrial cancer


CYP19A1
breast cancer


CYP19A1
prostate cancer


CYP19A1
endometriosis


CYP1A1
lung cancer


CYP1A1
breast cancer


CYP1A1
Colorectal Cancer


CYP1A1
prostate cancer


CYP1A1
esophageal cancer


CYP1A1
endometriosis


CYP1A1
cytogenetic studies


CYP1A2
schizophrenia


CYP1A2
colorectal cancer


CYP1B1
breast cancer


CYP1B1
glaucoma


CYP1B1
prostate cancer


CYP21A2
21-hydroxylase deficiency


CYP21A2
congenital adrenal hyperplasia


CYP21A2
adrenal hyperplasia, congenital


CYP2A6
smoking behavior


CYP2A6
nicotine


CYP2A6
lung cancer


CYP2C19

H. pylori infection



CYP2C19
phenytoin


CYP2C19
gastric disease


CYP2C8
malaria, plasmodium falciparum


CYP2C9
anticoagulant complications


CYP2C9
warfarin sensitivity


CYP2C9
warfarin therapy, response to


CYP2C9
colorectal cancer


CYP2C9
phenytoin


CYP2C9
acenocoumarol response


CYP2C9
coagulation disorder


CYP2C9
hypertension


CYP2D6
colorectal cancer


CYP2D6
Parkinson's disease


CYP2D6
CYP2D6 poor metabolizer phenotype


CYP2E1
lung cancer


CYP2E1
colorectal cancer


CYP3A4
prostate cancer


CYP3A5
prostate cancer


CYP3A5
esophageal cancer


CYP46A1
Alzheimer's Disease


DBH
schizophrenia


DHCR7
Smith-Lemli-Opitz syndrome


DISC1
schizophrenia


DLST
Alzheimer's Disease


DMD
muscular dystrophy


DRD2
alcoholism


DRD2
schizophrenia


DRD2
smoking behavior


DRD2
Parkinson's disease


DRD2
tardive dyskinesia


DRD3
schizophrenia


DRD3
tardive dyskinesia


DRD3
bipolar disorder


DRD4
attention deficit hyperactivity disorder


DRD4
schizophrenia


DRD4
novelty seeking


DRD4
ADHD


DRD4
personality traits


DRD4
heroin abuse


DRD4
alcohol abuse


DRD4
alcoholism


DRD4
personality disorders


DTNBP1
schizophrenia


EDN1
hypertension


EGFR
lung cancer


ELAC2
prostate cancer


ENPP1
Type 2 diabetes


EPHB2
prostate cancer


EPHX1
lung cancer


EPHX1
colorectal cancer


EPHX1
cytogenetic studies


EPHX1
chronic obstructive pulmonary disease/COPD


ERBB2
breast cancer


ERCC1
lung cancer


ERCC1
colorectal cancer


ERCC2
lung cancer


ERCC2
cytogenetic studies


ERCC2
bladder cancer


ERCC2
colorectal cancer


ESR1
bone density


ESR1
bone mineral density


ESR1
breast cancer


ESR1
endometriosis


ESR1
osteoporosis


ESR2
bone density


ESR2
breast cancer


estrogen receptor
bone mineral density


F2
coronary heart disease


F2
stroke


F2
thromboembolism, venous


F2
preeclampsia


F2
thrombosis


F5
thromboembolism, venous


F5
preeclampsia


F5
myocardial infarct


F5
stroke


F5
stroke, ischemic


F7
atherosclerosis, coronary


F7
myocardial infarct


F8
hemophilia


F9
hemophilia


FABP2
Type 2 diabetes


FAS
Alzheimer's Disease


FASLG
multiple sclerosis


FCGR2A
systemic lupus erythematosus


FCGR2A
lupus erythematosus


FCGR2A
periodontitis


FCGR2A
rheumatoid arthritis


FCGR2B
lupus erythematosus


FCGR2B
systemic lupus erythematosus


FCGR3A
systemic lupus erythematosus


FCGR3A
lupus erythematosus


FCGR3A
periodontitis


FCGR3A
arthritis


FCGR3A
rheumatoid arthritis


FCGR3B
periodontitis


FCGR3B
periodontal disease


FCGR3B
lupus erythematosus


FGB
fibrinogen


FGB
myocardial infarction


FGB
coronary heart disease


FLT3
leukemia, myeloid


FLT3
leukemia


FMR1
Fragile X syndrome


FRAXA
Fragile X Syndrome


FUT2

H. pylori infection



FVL
Factor V Leiden


G6PD
G6PD deficiency


G6PD
hyperbilirubinemia


GABRA5
bipolar disorder


GBA
Gaucher disease


GBA
Parkinson's disease


GCGR (FAAH,
body mass/obesity


ML4R, UCP2)


GCK
Type 2 diabetes


GCLM (F12,
atherosclerosis, myocardial infarction


TLR4)


GDNF
schizophrenia


GHRL
obesity


GJB1
Charcot-Marie-Tooth disease


GJB2
deafness


GJB2
hearing loss, sensorineural nonsyndromic


GJB2
hearing loss, sensorineural


GJB2
hearing loss/deafness


GJB6
hearing loss, sensorineural nonsyndromic


GJB6
hearing loss/deafness


GNAS
hypertension


GNB3
hypertension


GPX1
lung cancer


GRIN1
schizophrenia


GRIN2B
schizophrenia


GSK3B
bipolar disorder


GSTM1
lung cancer


GSTM1
colorectal cancer


GSTM1
breast cancer


GSTM1
prostate cancer


GSTM1
cytogenetic studies


GSTM1
bladder cancer


GSTM1
esophageal cancer


GSTM1
head and neck cancer


GSTM1
leukemia


GSTM1
Parkinson's disease


GSTM1
stomach cancer


GSTP1
Lung cancer


GSTP1
colorectal cancer


GSTP1
breast cancer


GSTP1
cytogenetic studies


GSTP1
prostate cancer


GSTT1
lung cancer


GSTT1
colorectal cancer


GSTT1
breast cancer


GSTT1
prostate cancer


GSTT1
Bladder Cancer


GSTT1
cytogenetic studies


GSTT1
asthma


GSTT1
benzene toxicity


GSTT1
esophageal cancer


GSTT1
head and neck cancer


GYS1
Type 2 diabetes


HBB
thalassemia


HBB
thalassemia, beta


HD
Huntington's disease


HFE
Hemochromatosis


HFE
iron levels


HFE
colorectal cancer


HK2
Type 2 diabetes


HLA
rheumatoid arthritis


HLA
Type1 diabetes


HLA
Behcet's Disease


HLA
celiac disease


HLA
psoriasis


HLA
Graves disease


HLA
multiple sclerosis


HLA
schizophrenia


HLA
asthma


HLA
diabetes mellitus


HLA
Lupus


HLA-A
leukemia


HLA-A
HIV


HLA-A
diabetes, type 1


HLA-A
graft-versus-host disease


HLA-A
multiple sclerosis


HLA-B
leukemia


HLA-B
Behcet's Disease


HLA-B
celiac disease


HLA-B
diabetes, type 1


HLA-B
graft-versus-host disease


HLA-B
sarcoidosis


HLA-C
psoriasis


HLA-DPA1
measles


HLA-DPB1
diabetes, type 1


HLA-DPB1
Asthma


HLA-DQA1
diabetes, type 1


HLA-DQA1
celiac disease


HLA-DQA1
cervical cancer


HLA-DQA1
asthma


HLA-DQA1
multiple sclerosis


HLA-DQA1
diabetes, type 2; diabetes, type 1


HLA-DQA1
lupus erythematosus


HLA-DQA1
pregnancy loss, recurrent


HLA-DQA1
psoriasis


HLA-DQB1
diabetes, type 1


HLA-DQB1
celiac disease


HLA-DQB1
multiple sclerosis


HLA-DQB1
cervical cancer


HLA-DQB1
lupus erythematosus


HLA-DQB1
pregnancy loss, recurrent


HLA-DQB1
arthritis


HLA-DQB1
asthma


HLA-DQB1
HIV


HLA-DQB1
lymphoma


HLA-DQB1
tuberculosis


HLA-DQB1
rheumatoid arthritis


HLA-DQB1
diabetes, type 2


HLA-DQB1
graft-versus-host disease


HLA-DQB1
narcolepsy


HLA-DQB1
arthritis, rheumatoid


HLA-DQB1
cholangitis, sclerosing


HLA-DQB1
diabetes, type 2; diabetes, type 1


HLA-DQB1
Graves' disease


HLA-DQB1
hepatitis C


HLA-DQB1
hepatitis C, chronic


HLA-DQB1
malaria


HLA-DQB1
malaria, plasmodium falciparum


HLA-DQB1
melanoma


HLA-DQB1
psoriasis


HLA-DQB1
Sjogren's syndrome


HLA-DQB1
systemic lupus erythematosus


HLA-DRB1
diabetes, type 1


HLA-DRB1
multiple sclerosis


HLA-DRB1
systemic lupus erythematosus


HLA-DRB1
rheumatoid arthritis


HLA-DRB1
cervical cancer


HLA-DRB1
arthritis


HLA-DRB1
celiac disease


HLA-DRB1
lupus erythematosus


HLA-DRB1
sarcoidosis


HLA-DRB1
HIV


HLA-DRB1
tuberculosis


HLA-DRB1
Graves' disease


HLA-DRB1
lymphoma


HLA-DRB1
psoriasis


HLA-DRB1
asthma


HLA-DRB1
Crohn's disease


HLA-DRB1
graft-versus-host disease


HLA-DRB1
hepatitis C, chronic


HLA-DRB1
narcolepsy


HLA-DRB1
sclerosis, systemic


HLA-DRB1
Sjogren's syndrome


HLA-DRB1
Type 1 diabetes


HLA-DRB1
arthritis, rheumatoid


HLA-DRB1
cholangitis, sclerosing


HLA-DRB1
diabetes, type 2; diabetes, type 1


HLA-DRB1

H. pylori infection



HLA-DRB1
hepatitis C


HLA-DRB1
juvenile arthritis


HLA-DRB1
leukemia


HLA-DRB1
malaria


HLA-DRB1
melanoma


HLA-DRB1
pregnancy loss, recurrent


HLA-DRB3
psoriasis


HLA-G
pregnancy loss, recurrent


HMOX1
atherosclerosis, coronary


HNF4A
diabetes, type 2


HNF4A
type 2 diabetes


HSD11B2
hypertension


HSD17B1
breast cancer


HTR1A
depressive disorder, major


HTR1B
alcohol dependence


HTR1B
alcoholism


HTR2A
memory


HTR2A
schizophrenia


HTR2A
bipolar disorder


HTR2A
depression


HTR2A
depressive disorder, major


HTR2A
suicide


HTR2A
Alzheimer's Disease


HTR2A
anorexia nervosa


HTR2A
hypertension


HTR2A
obsessive compulsive disorder


HTR2C
schizophrenia


HTR6
Alzheimer's Disease


HTR6
schizophrenia


HTRA1
wet age-related macular degeneration


IAPP
Type 2 Diabetes


IDE
Alzheimer's Disease


IFNG
tuberculosis


IFNG
Type 1 diabetes


IFNG
graft-versus-host disease


IFNG
hepatitis B


IFNG
multiple sclerosis


IFNG
asthma


IFNG
breast cancer


IFNG
kidney transplant


IFNG
kidney transplant complications


IFNG
longevity


IFNG
pregnancy loss, recurrent


IGFBP3
breast cancer


IGFBP3
prostate cancer


IL10
systemic lupus erythematosus


IL10
asthma


IL10
graft-versus-host disease


IL10
HIV


IL10
kidney transplant


IL10
kidney transplant complications


IL10
hepatitis B


IL10
juvenile arthritis


IL10
longevity


IL10
multiple sclerosis


IL10
pregnancy loss, recurrent


IL10
rheumatoid arthritis


IL10
tuberculosis


IL12B
Type 1 diabetes


IL12B
asthma


IL13
asthma


IL13
atopy


IL13
chronic obstructive pulmonary disease/COPD


IL13
Graves' disease


IL1A
periodontitis


IL1A
Alzheimer's Disease


IL1B
periodontitis


IL1B
Alzheimer's Disease


IL1B
stomach cancer


IL1R1
Type 1 diabetes


IL1RN
stomach cancer


IL2
asthma; eczema; allergic disease


IL4
Asthma


IL4
atopy


IL4
HIV


IL4R
asthma


IL4R
atopy


IL4R
Total serum IgE


IL6
Bone Mineralization


IL6
kidney transplant


IL6
kidney transplant complications


IL6
Longevity


IL6
multiple sclerosis


IL6
bone density


IL6
bone mineral density


IL6
Colorectal Cancer


IL6
juvenile arthritis


IL6
rheumatoid arthritis


IL9
asthma


INHA
premature ovarian failure


INS
Type 1 diabetes


INS
Type 2 diabetes


INS
diabetes, type 1


INS
obesity


INS
prostate cancer


INSIG2
obesity


INSR
Type 2 diabetes


INSR
hypertension


INSR
polycystic ovary syndrome


IPF1
diabetes, type 2


IRS1
Type 2 diabetes


IRS1
diabetes, type 2


IRS2
diabetes, type 2


ITGB3
myocardial infarction


ITGB3
atherosclerosis, coronary


ITGB3
coronary heart disease


ITGB3
myocardial infarct


KCNE1
EKG, abnormal


KCNE2
EKG, abnormal


KCNH2
EKG, abnormal


KCNH2
long QT syndrome


KCNJ11
diabetes, type 2


KCNJ11
Type 2 Diabetes


KCNN3
schizophrenia


KCNQ1
EKG, abnormal


KCNQ1
long QT syndrome


KIBRA
episodic memory


KLK1
hypertension


KLK3
prostate cancer


KRAS
colorectal cancer


LDLR
hypercholesterolemia


LDLR
hypertension


LEP
obesity


LEPR
obesity


LIG4
breast cancer


LIPC
atherosclerosis, coronary


LPL
Coronary Artery Disease


LPL
hyperlipidemia


LPL
triglycerides


LRP1
Alzheimer's Disease


LRP5
bone density


LRRK2
Parkinson's disease


LRRK2
Parkinsons disease


LTA
type 1 diabetes


LTA
Asthma


LTA
systemic lupus erythematosus


LTA
sepsis


LTC4S
Asthma


MAOA
alcoholism


MAOA
schizophrenia


MAOA
bipolar disorder


MAOA
smoking behavior


MAOA
personality disorders


MAOB
Parkinson's disease


MAOB
smoking behavior


MAPT
Parkinson's disease


MAPT
Alzheimer's Disease


MAPT
dementia


MAPT
Frontotemporal dementia


MAPT
progressive supranuclear palsy


MC1R
melanoma


MC3R
obesity


MC4R
obesity


MECP2
Rett syndrome


MEFV
Familial Mediterranean Fever


MEFV
amyloidosis


MICA
Type 1 diabetes


MICA
Behcet's Disease


MICA
celiac disease


MICA
rheumatoid arthritis


MICA
systemic lupus erythematosus


MLH1
colorectal cancer


MME
Alzheimer's Disease


MMP1
Lung Cancer


MMP1
ovarian cancer


MMP1
periodontitis


MMP3
myocardial infarct


MMP3
ovarian cancer


MMP3
rheumatoid arthritis


MPO
lung cancer


MPO
Alzheimer's Disease


MPO
breast cancer


MPZ
Charcot-Marie-Tooth disease


MS4A2
asthma


MS4A2
atopy


MSH2
colorectal cancer


MSH6
colorectal cancer


MSR1
prostate cancer


MTHFR
colorectal cancer


MTHFR
Type 2 diabetes


MTHFR
neural tube defects


MTHFR
homocysteine


MTHFR
thromboembolism, venous


MTHFR
atherosclerosis, coronary


MTHFR
Alzheimer's Disease


MTHFR
esophageal cancer


MTHFR
preeclampsia


MTHFR
pregnancy loss, recurrent


MTHFR
stroke


MTHFR
thrombosis, deep vein


MT-ND1
diabetes, type 2


MTR
colorectal cancer


MT-RNR1
hearing loss, sensorineural nonsyndromic


MTRR
neural tube defects


MTRR
homocysteine


MT-TL1
diabetes, type 2


MUTYH
colorectal cancer


MYBPC3
cardiomyopathy


MYH7
cardiomyopathy


MYOC
glaucoma, primary open-angle


MYOC
glaucoma


NAT1
colorectal cancer


NAT1
Breast Cancer


NAT1
bladder cancer


NAT2
colorectal cancer


NAT2
bladder cancer


NAT2
breast cancer


NAT2
Lung Cancer


NBN
breast cancer


NCOA3
breast cancer


NCSTN
Alzheimer's Disease


NEUROD1
Type 1 diabetes


NF1
neurofibromatosis1


NOS1
Asthma


NOS2A
multiple sclerosis


NOS3
hypertension


NOS3
coronary heart disease


NOS3
atherosclerosis, coronary


NOS3
coronary artery disease


NOS3
myocardial infarction


NOS3
acute coronary syndrome


NOS3
blood pressure, arterial


NOS3
preeclampsia


NOS3
nitric oxide


NOS3
Alzheimer's Disease


NOS3
asthma


NOS3
Type 2 diabetes


NOS3
cardiovascular disease


NOS3
Behcet's Disease


NOS3
erectile dysfunction


NOS3
kidney failure, chronic


NOS3
lead toxicity


NOS3
left ventricular hypertrophy


NOS3
pregnancy loss, recurrent


NOS3
retinopathy, diabetic


NOS3
stroke


NOTCH4
schizophrenia


NPY
alcohol abuse


NQO1
lung cancer


NQO1
colorectal cancer


NQO1
benzene toxicity


NQO1
bladder cancer


NQO1
Parkinson's Disease


NR3C2
hypertension


NR4A2
Parkinson's disease


NRG1
schizophrenia


NTF3
schizophrenia


OGG1
lung cancer


OGG1
colorectal cancer


OLR1
Alzheimer's Disease


OPA1
glaucoma


OPRM1
alcohol abuse


OPRM1
substance dependence


OPTN
glaucoma, primary open-angle


P450
drug metabolism


PADI4
rheumatoid arthritis


PAH
phenylketonuria/PKU


PAI1
coronary heart disease


PAI1
asthma


PALB2
breast cancer


PARK2
Parkinson's disease


PARK7
Parkinson's disease


PDCD1
lupus erythematosus


PINK1
Parkinson's disease


PKA
memory


PKC
memory


PLA2G4A
schizophrenia


PNOC
schizophrenia


POMC
obesity


PON1
atherosclerosis, coronary


PON1
Parkinson's disease


PON1
Type 2 Diabetes


PON1
atherosclerosis


PON1
coronary artery disease


PON1
coronary heart disease


PON1
Alzheimer's Disease


PON1
longevity


PON2
atherosclerosis, coronary


PON2
preterm delivery


PPARG
Type 2 Diabetes


PPARG
obesity


PPARG
diabetes, type 2


PPARG
Colorectal Cancer


PPARG
hypertension


PPARGC1A
diabetes, type 2


PRKCZ
Type 2 diabetes


PRL
systemic lupus erythematosus


PRNP
Alzheimer's Disease


PRNP
Creutzfeldt-Jakob disease


PRNP
Jakob-Creutzfeldt disease


PRODH
schizophrenia


PRSS1
pancreatitis


PSEN1
Alzheimer's Disease


PSEN2
Alzheimer's Disease


PSMB8
Type 1 diabetes


PSMB9
Type 1 diabetes


PTCH
skin cancer, non-melanoma


PTGIS
hypertension


PTGS2
colorectal cancer


PTH
bone density


PTPN11
Noonan syndrome


PTPN22
rheumatoid arthritis


PTPRC
multiple sclerosis


PVT1
end stage renal disease


RAD51
breast cancer


RAGE
retinopathy, diabetic


RB1
retinoblastoma


RELN
schizophrenia


REN
hypertension


RET
thyroid cancer


RET
Hirschsprung's disease


RFC1
neural tube defects


RGS4
schizophrenia


RHO
retinitis pigmentosa


RNASEL
prostate cancer


RYR1
malignant hyperthermia


SAA1
amyloidosis


SCG2
hypertension


SCG3
obesity


SCGB1A1
asthma


SCN5A
Brugada syndrome


SCN5A
EKG, abnormal


SCN5A
long QT syndrome


SCNN1B
hypertension


SCNN1G
hypertension


SERPINA1
COPD


SERPINA3
Alzheimer's Disease


SERPINA3
COPD


SERPINA3
Parkinson's disease


SERPINE1
myocardial infarct


SERPINE1
Type 2 Diabetes


SERPINE1
atherosclerosis, coronary


SERPINE1
obesity


SERPINE1
preeclampsia


SERPINE1
stroke


SERPINE1
hypertension


SERPINE1
pregnancy loss, recurrent


SERPINE1
thromboembolism, venous


SLC11A1
tuberculosis


SLC22A4
Crohn's disease; ulcerative colitis


SLC22A5
Crohn's disease; ulcerative colitis


SLC2A1
Type 2 diabetes


SLC2A2
Type 2 diabetes


SLC2A4
Type 2 diabetes


SLC3A1
cystinuria


SLC6A3
attention deficit hyperactivity disorder


SLC6A3
Parkinson's disease


SLC6A3
smoking behavior


SLC6A3
alcoholism


SLC6A3
schizophrenia


SLC6A4
depression


SLC6A4
depressive disorder, major


SLC6A4
schizophrenia


SLC6A4
suicide


SLC6A4
alcoholism


SLC6A4
bipolar disorder


SLC6A4
personality traits


SLC6A4
attention deficit hyperactivity disorder


SLC6A4
Alzheimer's Disease


SLC6A4
personality disorders


SLC6A4
panic disorder


SLC6A4
alcohol abuse


SLC6A4
affective disorder


SLC6A4
anxiety disorder


SLC6A4
smoking behavior


SLC6A4
depressive disorder, major; bipolar disorder


SLC6A4
heroin abuse


SLC6A4
irritable bowel syndrome


SLC6A4
migraine


SLC6A4
obsessive compulsive disorder


SLC6A4
suicidal behavior


SLC7A9
cystinuria


SNAP25
ADHD


SNCA
Parkinson's disease


SOD1
ALS/amyotrophic lateral sclerosis


SOD2
breast cancer


SOD2
lung cancer


SOD2
prostate cancer


SPINK1
pancreatitis


SPP1
multiple sclerosis


SRD5A2
prostate cancer


STAT6
asthma


STAT6
Total IgE


SULT1A1
breast cancer


SULT1A1
colorectal cancer


TAP1
Type 1 diabetes


TAP1
lupus erythematosus


TAP2
Type 1 diabetes


TAP2
diabetes, type 1


TBX21
asthma


TBXA2R
asthma


TCF1
diabetes, type 2


TCF1
Type 2 diabetes


TF
Alzheimer's Disease


TGFB1
breast cancer


TGFB1
kidney transplant


TGFB1
kidney transplant complications


TH
schizophrenia


THBD
myocardial infarction


TLR4
asthma


TLR4
Crohn's disease; ulcerative colitis


TLR4
sepsis


TNF
asthma


TNFA
cerebrovascular disease


TNF
Type 1 diabetes


TNF
rheumatoid arthritis


TNF
systemic lupus erythematosus


TNF
kidney transplant


TNF
psoriasis


TNF
sepsis


TNF
Type 2 Diabetes


TNF
Alzheimer's Disease


TNF
Crohn's disease


TNF
diabetes, type 1


TNF
hepatitis B


TNF
kidney transplant complications


TNF
multiple sclerosis


TNF
schizophrenia


TNF
celiac disease


TNF
obesity


TNF
pregnancy loss, recurrent


TNFRSF11B
bone density


TNFRSF1A
rheumatoid arthritis


TNFRSF1B
Rheumatoid Arthritis


TNFRSF1B
systemic lupus erythematosus


TNFRSF1B
arthritis


TNNT2
cardiomyopathy


TP53
lung cancer


TP53
breast cancer


TP53
Colorectal Cancer


TP53
prostate cancer


TP53
cervical cancer


TP53
ovarian cancer


TP53
smoking


TP53
esophageal cancer


TP73
lung cancer


TPH1
suicide


TPH1
depressive disorder, major


TPH1
suicidal behavior


TPH1
schizophrenia


TPMT
thiopurine methyltransferase activity


TPMT
leukemia


TPMT
inflammatory bowel disease


TPMT
thiopurine S-methyltransferase phenotype


TSC1
tuberous sclerosis


TSC2
tuberous sclerosis


TSHR
Graves' disease


TYMS
colorectal cancer


TYMS
stomach cancer


TYMS
esophageal cancer


UCHL1
Parkinson's disease


UCP1
obesity


UCP2
obesity


UCP3
obesity


UGT1A1
hyperbilirubinemia


UGT1A1
Gilbert syndrome


UGT1A6
colorectal cancer


UGT1A7
colorectal cancer


UTS2
diabetes, type 2


VDR
bone density


VDR
prostate cancer


VDR
bone mineral density


VDR
Type 1 diabetes


VDR
osteoporosis


VDR
bone mass


VDR
breast cancer


VDR
lead toxicity


VDR
tuberculosis


VDR
Type 2 diabetes


VEGF
breast cancer


vit D rec
idiopathic short stature


VKORC1
warfarin therapy, response to


WNK4
hypertension


XPA
lung cancer


XPC
lung cancer


XPC
cytogenetic studies


XRCC1
lung cancer


XRCC1
cytogenetic studies


XRCC1
breast cancer


XRCC1
bladder cancer


XRCC2
breast cancer


XRCC3
breast cancer


XRCC3
cytogenetic studies


XRCC3
lung cancer


XRCC3
bladder cancer


ZDHHC8
schizophrenia









Incorporating Ancestral Data

The present disclosure also provides methods and systems, such as described herein, that correlated phenotypes using genomic profiles by incorporating ancestral data. Thus, assessing an individual's genotype correlation may be expressed or reported as a GCI score, and incorporate ancestral data in generating the GCI score. For example, OR used in determining GCI scores may be modified based on an individual's ancestry or ethnicity.


The risk of an individual to develop a certain condition is typically based on the individual's genetics and environment. When trying to estimate the risk based on genetics, current studies can be limited by the fact that only a subset of all genetics markers or variations, such as SNPs, can be measured. Particularly, for complex diseases, the complex interaction of many genetic and environmental factors can lead to the development of a condition, and therefore there can be many genetic variations, such as SNPs, that marginally contribute to the risk. Current Whole-Genome-Association (WGA) studies normally consider each region in the genome in isolation, and try to answer the question as to what is the effect of a mutation in a specific SNP in that region on the risk for the condition, when keeping all other genetic factors and environmental factors as unknown. Mathematically, these studies essentially estimate the marginal distribution of the risk probability as a function of a SNP (these distributions as referred herein as the effect of a SNP).


The risk for developing the condition can be affected not by one genetic variation or SNP, but by many SNPs or other genetic variations, and environmental factors. Therefore, if two populations differ in their allelic distribution across the genome, and in the environmental factors affecting them, there may be a potential difference in the effect of a specific genetic variation, such as a SNP, in each of the populations. This is particularly the case when there is a gene-gene or gene-environment interaction between this SNP and another SNP, other genetic variants, or environmental factor. However, even in cases where there is no interaction, a different ‘background distribution’ of the other genetic and environmental factors can affect the effect of a genetic variation, such as a SNP. Thus, without being bound by theory, different populations can have different effect sizes for the same genetic variant, such as a SNP. In practice, however, almost all known conditions in which there is a SNP whose effect size was measured in more than one population, the effects measured were either very close to each other, or at least within the 95% CI of each other. As a result, in some embodiments, a simplifying assumption that can be used herein, is the effect size of a genetic variation, such as a causal SNP, is in fact the same across all populations.


Unfortunately, even with the assumption that the effect size is the same across populations, a limitation is the fact that the causal genetic variation may be unknown, for example, the causal SNP can not or has not been genotyped. Fortunately, SNPs or other genetic variations in close proximity on the genome can be correlated, such as in LD, and therefore even if the causal SNP is not measured, a tag SNP can be used as a proxy to the causal (see FIG. 10 for example). However, different populations can have different linkage disequilibrium patterns due to various possible reasons such as variation in recombination rates, selection pressure, or population bottleneck. Thus, in some embodiments, if a study has been done on population A, yielding a specific odds ratio in that population, the same odds ratio cannot be assumed in population B. This can be illustrated by the following example (see FIG. 10). For example, study has been performed on a Caucasian (CEU) population, and a large effect size has been reported for one of the SNP (the ‘published SNP’). In the example, the published SNP belongs to an LD block which is shared with the causal SNP, so the r̂2 (the square of the correlation coefficient) between the causal and the published SNP is 1; put differently, the published SNP and the causal SNPs are perfectly correlated in the CEU population. However, it may be the case that in another population (in this case, YRI, Yoruban), the published SNP and the causal SNP are in different LD blocks. In the extreme case, they can have r̂2=0, in which case they are independent of each other in that population. Under such a scenario, if the same study had been done on the YRI population, no effect would have been detected for the published SNP. It would therefore be wrong to estimate the risk of an YRI individual by ignoring the LD patterns in the underlying population of the customer and of the originally studied population. The example in FIG. 10 is an extreme case, but in reality similar patterns can happen with less extreme consequences.


Thus, in some embodiments, the present disclosure provides a method of assessing genotype correlations of an individual comprising comparing loci between populations of different ancestry. For example, odds ratios taken for a first population may be applied, or varied, to a second population, depending on such factors as LD patters. For example, for AS (Asians), the odds ratios used may be that of studies of AS, YRI (Yoruban), CEU (Caucasian/European) ancestry/ethnicity, in this order, since YRI has a lower LD than CEU. In some embodiments, locus-specific ancestry may be used for admixed populations.


In some embodiments, the populations of the first and second populations could comprise, but not be limited to any other population such as African American, Caucasian, Ashkenazi Jewish, Sepharadic Jewish, Indian, Pacific islanders, middle eastern, Druze, Bedouins, south Europeans, Scandinavians, eastern Europeans, North Africans, Basques, West Africans, East Africans. Otherwise stated, the populations of the first and second populations could comprise, but not limited to any of the HapMap populations (YRI,CEU,CHB,JPT, ASW, CHD, GIH, LWK, MEX, MKK,TSI). The description of the HapMap populations can be found in http://hapmap.org/hapmappopulations.html.en and in enclosed document.
















Number




of


Label
Population Sample
Samples

















ASW
African ancestry in Southwest USA
90


CEU
Utah residents with Northern and Western European
180



ancestry from the CEPH collection


CHB
Han Chinese in Beijing, China
90


CHD
Chinese in Metropolitan Denver, Colorado
100


GIH
Gujarati Indians in Houston, Texas
100


JPT
Japanese in Tokyo, Japan
91


LWK
Luhya in Webuye, Kenya
100


MEX
Mexican ancestry in Los Angeles, California
90


MKK
Maasai in Kinyawa, Kenya
180


TSI
Toscans in Italy
100


YRI
Yoruba in Ibadan, Nigeria
180









In some embodiments, methods for assessing an individual's genotype correlations to a phenotype may comprise comparing a first linkage disequilibrium (LD) pattern comprising a genetic variation, such as a SNP, correlated with a phenotype, wherein the first LD pattern is of a first population of individuals; and, a second LD pattern comprising the genetic variation (such as the SNP), wherein the second LD pattern is of a second population of individuals; determining a probability of the genetic variation being correlated with the phenotype in the second population from the comparing; and assessing a genotype correlation of the phenotype from a genomic profile of the individual comprising using the probability; and, reporting results comprising said genotype correlation from to said individual or a health care manager of said individual.


For example, assuming that a published SNP P has been reported for a first population, A, with odds ratios OR[P,A], and that the causal SNP C is unknown. Also in this example, is the assumption that for a second population, B, the odds ratios of C in A and B are the same, that is OR[C,A]=OR[C,B], if the sample size is large enough. Thus, if the location of C is known, and OR[C,A], the LD patterns in B can be used to estimate the best tag SNP to capture C in the population. However, in some embodiments, the location of C is unknown as is the odds ratio of C. However, for every SNP S and value X, the probability can be computed: Prob[S=C, OR[C,A]=X]|r̂2(P,S) in A, P, OR[P,A]], i.e., the probability that S is the causal SNP, with odds ratio of X (assuming an infinite sample size), given the correlation coefficient between S and P in population A, and given that P is the published SNP with odds ratio OR[P,A]. In order to calculate this probability, the fact that in the actual study the odds ratios of S is lower than the odds ratio of C is used and the question of what is the probability for this to happen given that OR[S,A] should approach X for a large enough sample size can be answered. Given the distribution of causal SNPs and their effect sizes, the expected effect size of a tag SNP can be determined by computing the expectation (the weighted average) of the effect sizes resulting from the different SNPs being causal.


In the example given in FIG. 10, since the LD block is of a perfect LD, all SNPs in the CEU block have the same probability of being causal with the same distribution of effect size (i.e., the log odds ratio is Normally distributed, where the confidence intervals determine its standard deviation). However, when the published SNP in YRI that aims at tagging the causal, the expected odds ratio of this SNP will be the weighted average between the published odds ratio and 1, where the weight corresponds to the length of the LD blocks involved.


Thus, in some embodiments, modifications of ORs may be determined by methods including, but not limited to, determining a causal genetic variation probability, such as an OR, for each of a plurality of genetic variations in a first population of individuals, or reference population, such as CEU as described in the above example. The OR may be then be used in assessing a genotype correlation from a genomic profile of an individual of a another population of individuals or reference group, such as YRI, reporting results comprising said genotype correlation from step (c) to said individual or a health care manager of said individual. Thus, the each of the genetic variations used in calculating their probability of being the causal genetic variant (such as a causal SNP), is typically proximal to a known genetic variation correlated to a phenotype in the first population, such as the published genetic variation, such as a published SNP. In some embodiments, each of each of the genetic variations used in calculating their probability of being the causal genetic variant (such as a causal SNP), is in linkage disequilibrium to the known or published genetic variation.


For example, again assuming that a published SNP P has been reported for a first population, A, with odds ratios OR[P,A], and that the causal SNP C is unknown, and that for a second population, B, the odds ratios of C in A and B are the same, that is OR[C,A]=OR[C,B], if the sample size is large enough. Another assumption in the example, the LD patterns are known for the studied population and for an individual's population. For example, it is assumed that the study has been done on the CEU population (first population), and that the individual is of the YRI population (second population), although the example can be extended to other populations. In every position, it is assumed that there is a risk allele R, and a non-risk allele N. The three possible genotypes in a given SNP are RR, RN, and NN. For a given SNP S, a genotype G (which is either RR, RN, or NN), and a group of individuals I, is denoted by F(S,I,G) the number of individuals in I with genotype G at SNP S. Thus, the odds ratios measured on the CEU population in the published SNP P is given by








OR


(

P
,
CEU
,
G

)


=



F


(

P
,
CA
,
G

)




F


(

P
,
CT
,
NN

)





F


(

P
,
CA
,
NN

)




F


(

P
,
CT
,
G

)





,




where CA and CT represent the case and control populations. Similarly, it is denoted by f(S,I,G) the frequency of the genotype G in population I. For a pair of SNPs S1 and S2, it is denoted by PCEU(S1,G1|S2, G2) the probability that an individual has genotype G1 at SNP S1, given that the individual has G2 at S2 (in CEU). A similar notation is used for YRI, the second population.


In some embodiments, an algorithm is used to determine an OR for the second population, and thus use in assessing an individual's genotype correlation to a phenotype, such as through the use of a GCI score. For example, the input and output of the algorithm disclosed herein may have the following information provided:


1) A list of SNPs (such as those disclosed in the HapMap), and another special SNP P, which is the published SNP.


2) The list of SNPs from the above SNPs that are measured in the study.


3) For the published SNP P, one of the following is assumed to be known:

    • (a) The genotype counts from the study for the cases and the controls, that is, the values of F(P,CA,G) and F(P,CT,G) are known for every genotype G.
    • (b) Alternatively, it is assumed that the genotypic odds ratios are known at SNP P, their confidence intervals, and the total number of cases and controls.


4) For every pair of SNPs S1,S2, and every pair of genotypes G1,G2, the algorithm will be provided with PCEU(S1,G1|S2, G2) and with PYRI(S1,G1|S2, G2)—this information can be found from the HapMap or other reference dataset.


The algorithm can then output for every SNP S in the proximity of P, the expected odds ratio of the SNP under the assumption that the number of individuals in the study is very large (approaching infinity). The algorithm will make the assumption that the odds ratios of the causal C for CEU (i.e. first population) and YRI (i.e. second population) approach the same number when the sample size approaches infinity.


Thus, the algorithm disclosed herein can include the following major steps (see also Example 5):


1) Find the LD probabilities based on the reference dataset.


2) Find the counts F(P,CA,G), F(P,CT,G) if they are not given as an input (alternative 1-b).


3) Sample n (n very large, e.g. >>1,000,000) instances for the genotype frequencies of the cases and controls at P, in CEU. The sampling is based on the posterior distribution of f(P,CA,G), f(P,CT,G), given the counts.


4) For each instance of the frequencies, and for each SNP S:

    • (a) Calculate f(S,CA,G) and f(S,CT,G) based on the frequencies in P and on PCEU.
    • (b) Generate an instance of F(S,CA,G), F(S,CT,G) based on the sampled allele frequencies in S.
    • (c) Calculate the p-value for S based on F(S), and based on f(S) (the latter is the p-value in the asymptotic sense).
    • (d) Find the min p-value based on F(S) across all measured SNPs S. If this is not P, then this instance is rejected.
    • (e) If the instance is not rejected, the instance is kept, together with the min p-value based on f(S); this will be the causal SNP of that instance.


5) The previous phase results in a set of causal SNPsC1, . . . , Cn, and their corresponding odds ratios. For each such causal, it is assumed that the same odds ratio holds for the YRI population.

    • (a) This information, together with the genotype frequencies in YRI at Ci is used to estimate the frequencies fYRI(Ci,CA,G) for every genotype G.
    • (b) The LD information is used to estimate fYRI(S,CA,G),fYRI(S,CT,G) for every SNP S, and calculate the asymptotic odds ratios based on these frequencies.
    • (c) For each SNP S, the asymptotic odds ratios are averaged across all instances, resulting in the expected asymptotic odds ratios.


In some embodiments, ancestral data may be used to assess an individual for their sub-group, for example, the present disclosure provides a method of assessing a reference sub-group of an individual comprising: obtaining a genetic sample of the individual; generating a genomic profile for the individual; determining the individual's one or more reference sub-groups by comparing the individual's genomic profile to a current database of human genotype correlations with ethnicity, geographic origin, or ancestry; and, reporting the results from step c) to the individual or a health care manager of the individual.


In one aspect a reference data set comprising multiple sets of genotyping data from individuals, wherein substantially the entire genome is used in the present disclosure. In one embodiment, the reference data contains genotyping data from substantially the entire genome of multiple individuals. Wherein in one embodiment, substantially the entire genome means that genetic markers are detected that cover at least 80% of an individuals genome, including but not limited to at least 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99% or 100% genomic coverage. In another embodiment at least 75% of the sets of genotyping data from the individuals included in the reference data include information from genetic markers that cover at least 80% of each individual's genome. In a further embodiment greater than 75% (including but not limited to greater than 76%, 77%, 78%, 79%, 80%, 81%, 82% 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99% or 100%) of the sets of genotyping data from the individuals included in the reference data include information from genetic markers that cover at least 80% of each individual's genome.


In one embodiment, the reference data set includes information on multiple genetic markers including but not limited to nucleotide repeats, nucleotide insertions, nucleotide deletions, chromosomal translocations, chromosomal duplications, copy number variations, microsatellite repeats, nucleotide repeats, centromeric repeats, or telomeric repeats or SNPs. In another embodiment the reference data set includes information which is substantially limited to a single genetic marker, such as SNPS or microsatellites. Wherein at least 80% of the genetic markers included in a reference set are of the same type, including but not limited to at least 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99% or 100% of the genetic markers.


In another embodiment, the reference data set consists essentially of whole genome SNP genotyping data. In some embodiments the SNP data is derived from analyses of indviduals' genomes using a high density DNA array for SNP identification and profile generation. Such arrays include but are not limited to those commercially available from Affymetrix and Illumina (see Affymetrix GeneChip® 500K Assay Manual, Affymetrix, Santa Clara, Calif. (incorporated by reference); Sentrix® humanHap650Y genotyping beadchip, Illumina, San Diego, Calif.). In some embodiments a reference set consists essentially of SNP data generated by genotyping more than 900,000 SNPs using the Affymetrix Genome Wide Human SNP Array 6.0. In alternative embodiment, more than 500,000 SNPs through whole-genome sampling analysis may be determined by using the Affymetrix GeneChip Human Mapping 500K Array Set.


In another embodiment, the reference data set contains information about the ethnicity, geographic origin and/or ancestry of each individual whose genotype data is included. In one embodiment said information is present in a reference data set, such as the HapMap or the Genographic Project (https://www3.nationalgeographic.com/genographic/). In another embodiment said information is self-reported, such as by subscribers or non-subscribers. In another embodiment subscribers may receive an incentive to self-report information about their ethnicity, geographic origin and/or ancestry. In another embodiment, subscribers may receive an incentive to self-report information about their disease status (such as information about any diseases or conditions they may display symptoms of or have a hereditary pre-disposition for). In another embodiment the individual receives an incentive to allow the use of this information and the individual's genotype in at least one reference data set. In some embodiments the incentive may be a financial incentive a discount on services offered, an offer of free services, an offer of a service upgrade (such as an increase in subscriber status, from basic to a premium membership category), an offer of free or discounted services for a relative, or an offer of discounted, free or credited services with a 3rd party vendor (such as Amazon, Starbucks, WebMD). In a related embodiment subscribers or non-subscribers who disclose information related to their ethnicity, geographic origin and/or ancestry, or disease status may be advised about the possible uses of such disclosed information and given the opportunity to supply or to withhold their informed consent.


In one embodiment, the reference data set contains information from multiple individuals with different ethnicities, geographic origins and/or ancestries. In another embodiment, the reference data set contains more than one individual from each class of ethnicity, geographic origin and/or ancestry represented in said reference data set. In another embodiment the reference data set contains more than five individual from each class of ethnicity, geographic origin and/or ancestry represented in said reference data set. In another embodiment the reference data set contains more than ten individuals from each class of ethnicity, geographic origin and/or ancestry represented in said reference data set. In another embodiment the reference data set contains more than twenty individuals from each class of ethnicity, geographic origin and/or ancestry represented in said reference data set.


In another embodiment, the assembled data in a reference set is analyzed to correlate ethnicity, geographic origin and/or ancestry, with at least one disease or condition and genetic marker associations. In another embodiment self-reported ethnicity, geographic origin and/or ancestry may be used to flag specific diseases or conditions for risk analysis. In another embodiment, an individual's ethnicity, geographic origin and/or ancestry is correlated with their genotype for further analysis (such as in silico population genetics studies) of associations between genetic markers and a disease or condition within a sub-grouping of individuals with a similar or shared ethnicity, geographic origin and/or ancestry. For example, it is known certain groups of individuals with a shared ethnicity, geographic origin and/or ancestry, such as the Ashkenazi Jews have a much higher likelihood of having children with diseases such as Tay Sachs. The analysis of an individual who self identifies as an Ashkenazi Jew could be modified to take this information into account when analyzing the individual's genetic markers.


In another embodiment, the data in reference data set can stratified into reference data sub-groups. A population, when considered as a whole, may contain multiple sub-groups, which may have different allele frequencies. The presence of multiple subgroups with different allele frequencies within a population can make association studies less informative. The different underlying allele frequencies in sampled subgroups may be independent of a disease or condition within each group, and they can lead to erroneous conclusions of linkage disequilibrium or disease relevance. Comparison of an individual's genotype to a reference data sub-group rather than to the entire reference data set can reduce the likelihood of errors created by spurious allelic associations. The data in each reference data sub-group may be organized by at least one shared feature, such as shared ethnicity, geographic origin and/or ancestry. The genotypes of individuals whose data is comprised within each sub-group can be further analyzed to identify common genetic markers that are of indicative of a specific ethnicity, geographic origin and/or ancestry. In an alternative embodiment assembled data in a reference set can be used to genetic markers which are associated with at least one disease or condition, and that are also associated with at least one ethnicity, geographic origin and/or ancestry.


In one embodiment, an individual's at least one self reported ethnic, geographic origin and/or ancestral trait is used to modify the analysis of the individual's genotype. A modified analysis may focus on genetic markers that are associated with a disease or condition, which are also common to at least one self identified ethnic, geographic origin and/or ancestral sub-group. In an alternative embodiment information related to an individual's ethnicity, geographic origin and/or ancestry is determined based on the individual's genotype. For example, an individual's genotype is compared to at least one reference data set and used to determine information about the individual's ethnicity, geographic origin and/or ancestry. This information is then incorporated into the analysis of the individual's genotype for association with at least one disease or condition. The analysis may focus on genetic markers associated with at least one disease or condition, which may also be common to at least one ethnicity, geographic origin and/or ancestry.


In another embodiment both information about an individual's ethnicity, geographic origin and/or ancestry, and information derived from analysis of the individual's genetic markers is used to determine the likelihood that the individual shares a specific ethnicity, geographic origin and/or ancestry. By combining both types of information the information obtained from the genotype analysis can be used to verify the individual's self-reported ethnicity, geographic origin and/or ancestry and to correct for any inaccuracies. In one embodiment information about an individual's ethnicity, geographic origin and/or ancestry is self-reported. In an alternative embodiment information about an individual's ethnicity, geographic origin and/or ancestry is estimated. Estimating an individual's ethnicity, geographic origin and/or ancestry can provide a continuous measure to assess population structure in the study of complex diseases or conditions. There can be a fair amount of heterogeneity in ethnic, geographic origin and/or ancestral groupings based on individuals' self reported information. For example, individual ethnic, geographic origin and/or ancestral proportions (such as European, North African, Aboriginal, etc.) can be estimated based on published allele frequencies. The estimated example individual ethnic, geographic origin and/or ancestral proportion can be used as a surrogate for self-reported information to investigate an association between at least one genetic marker and at least one disease or condition. Genetic risk models can then be used to determine if adjusting for an estimated individual ethnic, geographic origin and/or ancestral proportion provides a better fit to the data compared to a model with no adjustment for ethnicity, geographic origin and/or ancestry or one based on self-reported information. The model that provides the best fit can then be used to determine an individual's risk of acquiring at least one disease or condition.


In another embodiment, ethnicity, geographic origin and/or ancestry information from an individual that is based on genotype and/or self-reported data may be used to mathematically determine the closest reference sub-group or sub-groups to the individual, in terms of contribution to the individual's global genome. For example, if it may be determined that an individual's genotype suggests that he/she shares genetic markers indicative of more than one ethnicity, geographic origin and/or ancestry. This determination may include likelihoods, and optionally confidence intervals (such as there is an X %±Y), that at least one of an individual's relatives was from a specific ethnic, geographic and/or ancestral origin. This determination than can be used to inform an individual of the genetic markers typically associated with at least one disease or condition in individual's who share a similar ethnic, geographic and/or ancestral origin and their risk of acquiring said at least one disease or condition. In another embodiment a report may be generated which includes information on the contribution to an individual's entire genome from various ethnic sources, geographic origins and/or ancestral sources. For example a report may describe aggregate ancestral origins over an individual's entire genome in percentages, such as 20% from Africa, 30% from Asia, 50% from Europe. In a further embodiment such a report may optionally include confidence intervals (such as 20%±3 from Africa, 30%±5 from Asia, 50%±2 from Europe).


In another embodiment, an individual's determined ethnicity, geographic origin and/or ancestry may be used to determine an individual's risk of acquiring at least one disease or condition based on analysis of specific loci. In a related embodiment, a report may generated for at least one locus that characterizes the likelihood that an individual inherited said locus from a relative with a specific ethnicity, geographic origin and/or ancestry and the association of an allele at said locus with at least one disease or condition. In another embodiment at least two locus specific association results may be aggregated to determine an individual's combined risk of acquiring at least one disease or condition.


In another embodiment the risk of acquiring at least one disease or condition may be determined for an individual who has an ethnicity, geographic origin and/or ancestry that differs from those of individuals previously reported in association studies. In another embodiment the risk of acquiring at least one disease or condition may be determined for an individual who has a unique or rare ethnicity, geographic origin and/or ancestry that makes it difficult or impossible to find a reference data sub-group to compare the individual's genotype to. For example an individual may want to know his/her risk of acquiring an inherited disease which may directly related to his ethnicity, geographic origin and/or ancestry. Some exceedingly rare diseases, such as oculopharyngeal muscular dystrophy, are found only within small localized groups in a population. Often diseases of this nature can be traced back to a single founder or to a limited number of past disease carriers. For diseases of this nature it is often possible to exclude an individual from an at-risk group if it can be determined that the individual is not related to the original founder or disease carriers. In one embodiment it may beneficial to conduct one or more association studies of other individuals with a shared genetic background or shared ethnicity, geographic origin and/or ancestry. Wherein the individuals' ethnicity, geographic origin and/or ancestry is determined by estimation or by self-reported information. These studies can combine information on individual's genotype, ethnicity, geographic origin and/or ancestry, and status of at least one disease or condition. Results obtained from at least two studies can be compared to determine if a similar association between an allele of a genetic marker and at least one disease or condition is observed. Results may depend on the correlation structure and allele frequencies in each of the populations studied and the relationship between them. Further, said studies can be used to identify genetic markers that are associated with susceptibility to said at least one disease or condition. In one embodiment the absence of at least one allele for at least one genetic marker is used to exclude an individual from being at risk for at least one disease or condition. In an alternative embodiment the presence of at least one allele for at least one genetic marker is used to categorize an individual as being at risk for at least one disease or condition.


The following examples illustrate and explain the disclosure. The scope of the disclosure is not limited by these examples.


EXAMPLES
Example 1
Generation and Analysis of SNP Profile

The individual is provided a sample tube in the kit, such as that available from DNA Genotek, into which the individual deposits a sample of saliva (approximately 4 mls) from which genomic DNA will be extracted. The saliva sample is sent to a CLIA certified laboratory for processing and analysis. The sample is typically sent to the facility by overnight mail in a shipping container that is conveniently provided to the individual in the collection kit.


In a preferred embodiment, genomic DNA is isolated from saliva. For example, using DNA self collection kit technology available from DNA Genotek, an individual collects a specimen of about 4 ml saliva for clinical processing. After delivery of the sample to an appropriate laboratory for processing, DNA is isolated by heat denaturing and protease digesting the sample, typically using reagents supplied by the collection kit supplier at 50° C. for at least one hour. The sample is next centrifuged, and the supernatant is ethanol precipitated. The DNA pellet is suspended in a buffer appropriate for subsequent analysis.


The individual's genomic DNA is isolated from the saliva sample, according to well known procedures and/or those provided by the manufacturer of a collection kit. Generally, the sample is first heat denatured and protease digested. Next, the sample is centrifuged, and the supernatant is retained. The supernatant is then ethanol precipitated to yield a pellet containing approximately 5-16 ug of genomic DNA. The DNA pellet is suspended in 10 mM Tris pH 7.6, 1 mM EDTA (TE). A SNP profile is generated by hybridizing the genomic DNA to a commercially available high density SNP array, such as those available from Affymetrix or Illumina, using instrumentation and instructions provided by the array manufacturer. The individual's SNP profile is deposited into a secure database or vault.


The patient's data structure is queried for risk-imparting SNPs by comparison to a clinically-derived database of established, medically relevant SNPs whose presence in a genome correlates to a given disease or condition. The database contains information of the statistical correlation of particular SNPs and SNP haplotypes to particular diseases or conditions. For example, as shown in Example III, polymorphisms in the apolipoprotein E gene give rise to differing isoforms of the protein, which in turn correlate with a statistical likelihood of developing Alzheimer's Disease. As another example, individuals possessing a variant of the blood clotting protein Factor V known as Factor V Leiden have an increased tendency to clot. A number of genes in which SNPs have been associated to a disease or condition phenotype are shown in Table 1. The information in the database is approved by a research/clinical advisory board for its scientific accuracy and importance, and may be reviewed with governmental agency oversight. The database is continually updated as more SNP-disease correlations emerge from the scientific community.


The results of the analysis of an individual's SNP profile is securely provided to patient by an on-line portal or mailings. The patient is provided interpretation and supportive information, such as the information shown for Factor V Leiden in Example IV. Secure access to the individual's SNP profile information, such as through an on-line portal, will facilitate discussions with the patient's physician and empower individual choices for personalized medicine.


Example 2
Update of Genotype Correlations

In response to a request for an initial determination of an individual's genotype correlations, a genomic profile is generated, genotype correlations are made, and the results are provided to the individual as described in Example I. Following an initial determination of an individual's genotype correlations, subsequent, updated correlations are or can be determined as additional genotype correlations become known. The subscriber has a premium level subscription and their genotype profile and is maintained in a secure database. The updated correlations are performed on the stored genotype profile.


For example, an initial genotype correlation, such as described above in Example I, could have determined that a particular individual does not have ApoE4 and thus is not predisposed to early-onset Alzheimer's Disease, and that this individual does not have Factor V Leiden. Subsequent to this initial determination, a new correlation could become known and validated, such that polymorphisms in a given gene, hypothetically gene XYZ, are correlated to a given condition, hypothetically condition 321. This new genotype correlation is added to the master database of human genotype correlations. An update is then provided to the particular individual by first retrieving the relevant gene XYZ data from the particular individual's genomic profile stored in a secure database. The particular individual's relevant gene XYZ data is compared to the updated master database information for gene XYZ. The particular individual's susceptibility or genetic predisposition to condition 321 is determined from this comparison. The results of this determination are added to the particular individual's genotype correlations. The updated results of whether or not the particular individual is susceptible or genetically predisposed to condition 321 is provided to the particular individual, along with interpretative and supportive information.


Example 3
Correlation of ApoE4 Locus and Alzheimer's Disease

The risk of Alzheimer's disease (AD) has been shown to correlate with polymorphisms in the apolipoprotein E (APOE) gene, which gives rise to three isoforms of APOE referred to as ApoE2, ApoE3, and ApoE4. The isoforms vary from one another by one or two amino acids at residues 112 and 158 in the APOE protein. ApoE2 contains 112/158 cys/cys; ApoE3 contains 112/158 cys/arg; and ApoE4 contains 112/158 arg/arg. As shown in Table 2, the risk of Alzeimer's disease onset at an earlier age increases with the number of APOE ε4 gene copies. Likewise, as shown in Table 3, the relative risk of AD increases with number of APOE ε4 gene copies.









TABLE 2







Prevalence of AD Risk Alleles (Corder et al., Science: 261: 921-3, 1993)














Alzheimer's




APOE ε4 Copies
Prevalence
Risk
Onset Age







0
73%
20%
84



1
24%
47%
75



2
 3%
91%
68

















TABLE 3







Relative Risk of AD with ApoE4


(Farrer et al., JAMA: 278: 1349-56, 1997)










APOE Genotype
Odds Ratio














ε2ε2
0.6



ε2ε3
0.6



ε3ε3
1.0



ε2ε4
2.6



ε3ε4
3.2



ε4ε4
14.9










Example 4
Information for Factor V Leiden Positive Patient

The following information is exemplary of information that could be supplied to an individual having a genomic SNP profile that shows the presence of the gene for Factor V Leiden. The individual may have a basic subscription in which the information may be supplied in an initial report.


What is Factor V Leiden?

Factor V Leiden is not a disease, it is the presence of a particular gene that is passed on from one's parents. Factor V Leiden is a variant of the protein Factor V (5) which is needed for blood clotting. People who have a Factor V deficiency are more likely to bleed badly while people with Factor V Leiden have blood that has an increased tendency to clot.


People carrying the Factor V Leiden gene have a five times greater risk of developing a blood clot (thrombosis) than the rest of the population. However, many people with the gene will never suffer from blood clots. In Britain and the United States, 5 percent of the population carry one or more genes for Factor V Leiden, which is far more than the number of people who will actually suffer from thrombosis.


How do You Get Factor V Leiden?

The genes for the Factor V are passed on from one's parents. As with all inherited characteristics, one gene is inherited from the mother and one from the father. So, it is possible to inherit:—two normal genes or one Factor V Leiden gene and one normal gene—or two Factor V Leiden genes. Having one Factor V Leiden gene will result in a slightly higher risk of developing a thrombosis, but having two genes makes the risk much greater.


What are the Symptoms of Factor V Leiden?

There are no signs, unless you have a blood clot (thrombosis).


What are the Danger Signals?

The most common problem is a blood clot in the leg. This problem is indicated by the leg becoming swollen, painful and red. In rarer cases a blood clot in the lungs (pulmonary thrombosis) may develop, making it hard to breathe. Depending on the size of the blood clot this can range from being barely noticeable to the patient experiencing severe respiratory difficulty. In even rarer cases the clot might occur in an arm or another part of the body. Since these clots formed in the veins that take blood to the heart and not in the arteries (which take blood from the heart), Factor V Leiden does not increase the risk of coronary thrombosis.


What can be Done to Avoid Blood Clots?

Factor V Leiden only slightly increases the risk of getting a blood clot and many people with this condition will never experience thrombosis. There are many things one can do to avoid getting blood clots. Avoid standing or sitting in the same position for long periods of time. When traveling long distances, it is important to exercise regularly—the blood must not ‘stand still’. Being overweight or smoking will greatly increase the risk of blood clots. Women carrying the Factor V Leiden gene should not take the contraceptive pill as this will significantly increase the chance of getting thrombosis. Women carrying the Factor V Leiden gene should also consult their doctor before becoming pregnant as this can also increase the risk of thrombosis.


How does a Doctor Find Out if You have Factor V Leiden?


The gene for Factor V Leiden can be found in a blood sample.


A blood clot in the leg or the arm can usually be detected by an ultrasound examination.


Clots can also be detected by X-ray after injecting a substance into the blood to make the clot stand out. A blood clot in the lung is harder to find, but normally a doctor will use a radioactive substance to test the distribution of blood flow in the lung, and the distribution of air to the lungs. The two patterns should match—a mismatch indicates the presence of a clot.


How is Factor V Leiden Treated?

People with Factor V Leiden do not need treatment unless their blood starts to clot, in which case a doctor will prescribe blood-thinning (anticoagulant) medicines such as warfarin (e.g. Marevan) or heparin to prevent further clots. Treatment will usually last for three to six months, but if there are several clots it could take longer. In severe cases the course of drug treatment may be continued indefinitely; in very rare cases the blood clots may need to be surgically removed.


How is Factor V Leiden Treated During Pregnancy?

Women carrying two genes for Factor V Leiden will need to receive treatment with a heparin coagulant medicine during pregnancy. The same applies to women carrying just one gene for Factor V Leiden who have previously had a blood clot themselves or who have a family history of blood clots.


All women carrying a gene for Factor V Leiden may need to wear special stockings to prevent clots during the last half of pregnancy. After the birth of the child they may be prescribed the anticoagulant drug heparin.


Prognosis

The risk of developing a clot increases with age, but in a survey of people over the age of 100 who carry the gene, it was found that only a few had ever suffered from thrombosis. The National Society for Genetic Counselors (NSGC) can provide a list of genetic counselors in your area, as well as information about creating a family history. Search their on-line database at www.nsgc.org/consumer.


Example 5
Generating Odds Ratio for an Individual of a Different Ancestry

1. Find the LD probabilities based on the reference dataset.


The number of HapMap individuals with genotype pair (G1,G2) at SNPs S1,S2 is counted to generate the joint distribution of the two SNPs. The marginal distributions of each of the SNPs is combined, using Bayes law, to estimate PCEU(S1,G1|S2, G2) (CEU is the published, or first population) and with PYRI(S1,G1|S2, G2) (YRI is the second population, ancestry of the individual)


2. Find the counts F(P, CA, G), F(P, CT, G) if they are not given as an input (alternative 1-b).


If the counts are not given as an input, the following set of equations is used to find them:








F


(

P
,
CA
,
NN

)


+

F


(

P
,
CA
,
NR

)


+

F


(

P
,
CA
,
RR

)



=
N








F


(

P
,
CT
,
NN

)


+

F


(

P
,
CT
,
NR

)


+

F


(

P
,
CT
,
RR

)



=
M







OR


(

P
,
CEU
,
RR

)


=





F


(

P
,
CA
,
RR

)




F


(

P
,
CT
,
NN

)





F


(

P
,
CA
,
NN

)




F


(

P
,
CT
,
RR

)




.





OR


(

P
,
CEU
,
RN

)



=



F


(

P
,
CA
,
RN

)




F


(

P
,
CT
,
NN

)





F


(

P
,
CA
,
NN

)




F


(

P
,
CT
,
RN

)













1

F


(

P
,
CA
,
NN

)



+

1

F


(

P
,
CA
,
RR

)



+

1

F


(

P
,
CT
,
NN

)



+

1

F


(

P
,
CT
,
RR

)




=


(


log


(


UB


(

P
,
CEU
,
RR

)



OR


(

P
,

CEU
.
RR


)



)


1.96

)

2









1

F


(

P
,
CA
,
NN

)



+

1

F


(

P
,
CA
,
RN

)



+

1

F


(

P
,
CT
,
NN

)



+

1

F


(

P
,
CT
,
RN

)




=


(


log


(


UB


(

P
,
CEU
,
RN

)



OR


(

P
,

CEU
.
RN


)



)


1.96

)

2





In the above equations, UB(P,CEU,G) is the upper bound on the confidence interval of the odds ratios for genotype G at the published SNP P. M and N are the number of controls and cases in the study respectively.


These are six equations with six variables. Enumeration over all values of F(P,CA,NN) and F(P,CA,RN) is performed. For each such pair of values, 2-4 equations determining the rest of the variables is present by solving a set of linear equations, and the last two equations are used for validation. The running time is bounded by N2.


3. Sample n (n very large, e.g. >>1,000,000) instances for the genotype frequencies of the cases and controls at P, in CEU. The sampling is based on the posterior distribution of f(P, CA, G), f(P, CT, G), given the counts.


Given f(P), the likelihood of seeing F(P) can be calculated under the assumption of a Multinomial distribution. By assuming a uniform prior on the possible values of f(P), it is known that the probability Prob(f(P)|F(P)) αProb(F(P)|f(P)). An MCMC approach is used to sample from this distribution using a Gibbs Sampler.


4. For each instance of the frequencies, and for each SNP S:


a) Calculate f(S, CA, G) and f(S, CT, G) based on the frequencies in P and on PCEU. The formula







f


(

S
,
CA
,
G

)


=




G






f


(

P
,
CA
,

G



)





P
CEU



(

S
,

G

P

,

G



)








is used to estimate the frequencies at S in the cases. A similar formula can be used for the controls.


b) Generate an instance of F(S,CA,G), F(S,CT,G) based on the sampled allele frequencies in S. This is done by assuming a multinomial random variable that represents the genotype at S.


c) Calculate the p-value for S based on F(S), and based on f(S) (the latter is the p-value in the asymptotic sense).


The Armitage-Trend test is used to calculate the p-value based on F(S). In order to calculate the asymptotic p-value, it is assumed a sample size of N cases and N controls, with counts that match the expectation, e.g., F(S,CA,G) will be assumed to be Nf(S,CA,G).


d) Find the min p-value based on F(S) across all measured SNPs S. If this is not P, then this instance is rejected.


e) If the instance is not rejected, the instance is kept, together with the min p-value based on f(S); this is the causal SNP of that instance.


5. The previous phase results in a set of causal SNPsC1, . . . , Cn, and their corresponding odds ratios. For each such causal, it is assumed that the same odds ratio holds for the YRI population.


a) This information is used, together with the genotype frequencies in YRI at Ci to estimate the frequencies fYRI(Ci,CA,G) for every genotype G.


To do so, the following equation is solved:









f
YRI



(

S
,
CA
,
NN

)


+


f

YR





1




(

S
,
CA
,
RN

)


+


f
YRI



(

S
,
CA
,
RR

)



=
1







OR


(
RR
)


=




f
YRI



(

S
,
CA
,
RR

)





f
YRI



(

S
,
CT
,
NN

)






f
YRI



(

S
,
CA
,
NN

)





f
YRI



(

S
,
CT
,
TT

)











OR


(
RN
)


=




f
YRI



(

S
,
CA
,
RN

)





f
YRI



(

S
,
CT
,
NN

)






f
YRI



(

S
,
CA
,
NN

)





f
YRI



(

S
,
CT
,
RN

)








There are three missing variables in this equation, since fYRI(S,CT,G) are assumed to be known from the reference population (HapMap). The above set of equations is therefore a set of linear equations and can be solved efficiently.


b) The LD information is used to estimate fYRI(S,CA,G), fYRI(S,CT,G) for every SNP S, and calculate the asymptotic odds ratios based on these frequencies. This is done in a similar manner to step 4(a).


c) For each SNP S, the asymptotic odds ratios is averaged across all instances, resulting in the expected asymptotic odds ratios.


The odds ratios can then be used in determining an individual's genotype correlation.


While preferred embodiments of the present disclosure have been shown and described herein, it will be obvious to those skilled in the art that such embodiments are provided by way of example only. Numerous variations, changes, and substitutions will now occur to those skilled in the art without departing from the disclosure. It should be understood that various alternatives to the embodiments of the disclosure described herein may be employed in practicing the disclosure. It is intended that the following claims define the scope of the disclosure and that methods and structures within the scope of these claims and their equivalents be covered thereby.

Claims
  • 1. A method of assessing genotype correlations of an individual to a phenotype comprising: (a) comparing: (i) a first linkage disequilibrium (LD) pattern comprising a genetic variation correlated with a phenotype, wherein said first LD pattern is of a first population of individuals; and,(ii) a second LD pattern comprising said genetic variation, wherein said second LD pattern is of a second population of individuals;(b) determining a probability of said genetic variation being correlated with said phenotype in said second population from said comparing in (a);(c) assessing a genotype correlation of said phenotype from a genomic profile of said individual comprising using said probability of step (b); and,(d) reporting results comprising said genotype correlation from step c) to said individual or a health care manager of said individual.
  • 2. The method of claim 1, wherein in step (b), said probability is either an allelic or genotypic odds ratio (OR).
  • 3. The method of claim 2, wherein said OR is derived from a known OR, wherein said known OR is for said genetic variation correlated with said phenotype for said first population.
  • 4. The method of claim 2, wherein said first population and said second population have similar LD patterns.
  • 5. A method of assessing genotype correlations of an individual to a phenotype comprising: (a) determining a causal genetic variation probability for each of a plurality of genetic variations in a first population of individuals;(b) identifying each of said probability in step (a) as a probability for each of said plurality of genetic variations in a second population of individuals;(c) assessing a genotype correlation to a phenotype from a genomic profile of said individual comprising using said probability of step (b); and,(d) reporting results comprising said genotype correlation to a phenotype from step (c) to said individual or a health care manager of said individual.
  • 6. The method of claim 5, wherein said probability in step (a) is an OR.
  • 7. The method of claim 5, wherein each of said genetic variations of step (a) is proximal to a known genetic variation correlated to a phenotype in said first population.
  • 8. The method of claim 7, wherein each of said genetic variations of step (a) is in linkage disequilibrium to said known genetic variation.
  • 9. The method of claim 1 or 5, wherein said genotype correlation to a phenotype is reported as a GCI score.
  • 10. The method of claim 1 or 5, wherein said second population is of an ancestry different from said first population.
  • 11. The method of claim 1 or 5, wherein said individual is of an ancestry of said second population.
  • 12. The method of claim 1 or 5, wherein a causal genetic variation is unknown.
  • 13. The method of claim 1 or 5, wherein said genetic variation is single nucleotide polymorphism (SNP).
  • 14. The method of claim 1 or 5, wherein said reporting comprises transmission of said results over a network.
  • 15. The method of claim 1 or 5, wherein said reporting is through an on-line portal.
  • 16. The method of claim 1 or 5, wherein said reporting is by paper or by e-mail.
  • 17. The method of claim 1 or 5, wherein said reporting comprises reporting in a secure manner.
  • 18. The method of claim 1 or 5, wherein said reporting comprises reporting in a non-secure manner.
  • 19. The method of claim 1 or 5, wherein generating said genomic profile is by a third party.
  • 20. The method of claim 1 or 5, wherein said genomic profile is generated from a genetic sample.
  • 21. The method of claim 20, wherein a third party obtains said genetic sample.
  • 22. The method of claim 20, wherein said genetic sample is DNA.
  • 23. The method of claim 20, wherein said genetic sample is RNA.
  • 24. The method of claim 20, wherein said genetic sample is from a biological sample selected from the group consisting of: blood, hair, skin, saliva, semen, urine, fecal material, sweat, and buccal sample.
  • 25. The method of claim 1 or 5, wherein said genomic profile is deposited into a secure database or vault.
  • 26. The method of claim 1 or 5, wherein said genomic profile is a single nucleotide polymorphism profile.
  • 27. The method of claim 1 or 5, wherein said genomic profile comprises truncations, insertions, deletions, or repeats.
  • 28. The method of claim 1 or 5, wherein said genomic profile is generated using a high density DNA microarray.
  • 29. The method of claim 1 or 5, wherein said genomic profile is generated using RT-PCR.
  • 30. The method of claim 1 or 5, wherein said genomic profile is generated using DNA sequencing.
  • 31. The method of claim 1 or 5, further comprising (e) updating said results with additional genetic variations.
  • 32. The method of claim 1 or 5 wherein the populations of claims 1 or 2 comprise any of the HapMap populations (YRI,CEU,CHB,JPT,ASW,CHD,GIH,LWK,MEX,MKK,TSI), or to any other population such as, but not limited to African American, Caucasian, Ashkenazi Jewish, Sepharadic Jewish, Indian, Pacific islanders, middle eastern, Druze, Bedouins, south Europeans, Scandinavians, eastern Europeans, North Africans, Basques, West Africans, or East Africans.
CROSS-REFERENCE

This application claims the benefit of U.S. provisional application 60/975,495 filed Sep. 26, 2007, which is herein incorporated by reference in its entirety.

Provisional Applications (1)
Number Date Country
60975495 Sep 2007 US