Parameters for reliable results in genetic association studies in common disease

Genetic Variants in IL-12B and IL-27 in the Polish Patients with Systemic Lupus Erythematosus

To investigate the potential association between IL-12B and IL-27 gene polymorphisms and systemic lupus erythematosus (SLE), we performed a case-control study based on the Polish population. Patients with SLE and healthy individuals were examined for -6415 CTCTAA/GC (rs17860508) and +1188A/C (rs3212227) in IL-12B and -924A/G (rs153109) and 4730T/C (rs181206) in IL-27 gene polymorphisms using the high-resolution melting method, PCR-RFLP method and TaqMan SNP genotyping assay, respectively. An increased frequency of GC/GC genotype as well as GC allele of the IL-12B rs17860508 was found in patients with SLE, as compared with healthy subjects (P < 0.001). We did not find differences in genotype and allele frequencies of the IL-12B rs3212227 and IL-27 rs153109 and rs181206 variants between patients with SLE and controls. IL-27 haplotype rs181206C/rs153109G indicated higher risk for SLE (P = 0.002), whereas haplotype rs181206T/rs153109G indicated reduced risk for SLE (P = 0.005). The IL-12B rs3212227 A/C polymorphism was associated with the mean value of the platelets (PLT), urea and complement C3 level. Furthermore, IL-12B rs17860508 genetic variant showed correlation with PLT, prothrombin time, international normalized ratio and alkaline phosphatase. Our results revealed that IL-12B rs17860508 and IL-27 haplotype CG are genetic risk factors for SLE and that both IL-12B rs17860508 and rs3212227 predict disease phenotype.

Apolipoprotein E ε4 Allele was Associated With Nonlesional Mesial Temporal Lobe Epilepsy in Han Chinese Population

Apolipoprotein E (APOE) gene has been implicated as one of the genes susceptible to temporal lobe epilepsy (TLE), but the association is inconsistent. We carried out a study to investigate the association of APOEε4 allele with a subtype of TLE-nonlesional mesial temporal lobe epilepsy (NLMTLE) in Han Chinese people.T he study consisted of total 308 NLMTLE patients and 302 controls in Han Chinese. The APOE polymorphisms were genotyped using polymerase chain reaction (PCR) DNA sequencing. We compared the frequency of APOEε4 allele and carrying status between NLMTLE patients and control subjects to test for the association of APOEε4 allele with NLMTLE clinical status. Carrying status of APOEε4 allele was significantly associated with the risk of NLMTLE. No effect of APOEε4 allele was found on the age of onset, duration of epilepsy, or frequency of seizure. Moreover, there was no association between APOEε4 allele and hippocampal sclerosis (HS) or febrile convulsion (FC) history.O ur study provided an evidence that APOEε4 allele was a possible risk factor for NLMTLE, and further study with a larger sample is needed to warrant this finding.

A distinct and replicable variant of the squamous cell carcinoma gene inositol polyphosphate-5-phosphatase modifies the susceptibility of arsenic-associated skin lesions in Bangladesh

BACKGROUND

Single-nucleotide polymorphisms (SNPs) in inflammation, one-carbon metabolism, and skin cancer genes might influence susceptibility to arsenic-induced skin lesions.

METHODS

A case-control study was conducted in Pabna, Bangladesh (2001-2003), and the drinking-water arsenic concentration was measured for each participant. A panel of 25 candidate SNPs was analyzed in 540 cases and 400 controls. Logistic regression was used to estimate the association between each SNP and the potential for gene-environment interactions in the skin lesion risk, with adjustments for relevant covariates. Replication testing was conducted in an independent Bangladesh population with 488 cases and 2,794 controls.

RESULTS

In the discovery population, genetic variants in the one-carbon metabolism genes phosphatidylethanolamine N-methyltransferase (rs2278952, P for interaction  = .004; rs897453, P for interaction = .05) and dihydrofolate reductase (rs1650697, P for interaction = .02), the inflammation gene interleukin 10 (rs3024496, P for interaction =.04), and the skin cancer genes inositol polyphosphate-5-phosphatase (INPP5A; rs1133400, P for interaction = .03) and xeroderma pigmentosum complementation group C (rs2228000, P for interaction = .01) significantly modified the association between arsenic and skin lesions after adjustments for multiple comparisons. The significant gene-environment interaction between a SNP in the INPP5A gene (rs1133400) and water arsenic with respect to the skin lesion risk was successfully replicated in an independent population (P for interaction = .03).

CONCLUSIONS

Minor allele carriers of the skin cancer gene INPP5A modified the odds of arsenic-induced skin lesions in both main and replicative populations. Genetic variation in INPP5A appears to have a role in susceptibility to arsenic toxicity.

Integration of genomic information in the clinical management of HCC

Molecular profiling of hepatocellular carcinoma (HCC) is enabling the advancement of novel approaches to disease diagnosis and management. Accurate prognosis prediction in HCC is specially critical. Clinical staging systems for HCC support clinical decision-making (e.g., BCLC algorithm) might be complemented by molecular-based information in the near future. Molecular signatures derived from tumour and non-tumour samples are associated with patient recurrence an outcome. Single nucleotide polymorphisms have been linked with HCC development. Next generation sequencing studies have brought to light the genomic diversity of this disease. Gens recurrently altered in HCC and susceptible to be targeted belong to signalling pathways including telomere maintenance, cell cycle, chromatin remodelling, Wnt/beta-catenin, RAS/RAF/MAPK and PI3K/AKT/mTOR pathways. Oncogenic loops are unknown but might include some of the already discovered aberrations. Despite the intratumoral heterogeneity observed in HCC tumours, studies including large number of samples can identify key genetic drivers and contribute to the development of novel treatments and a personalized medicine.

Detecting signals in pharmacogenomic genome-wide association studies

In one common pharmacogenomic scenario, outcome measures are compared for treated and untreated subjects across genotype defined subgroups. The key question is whether treatment benefit (or harm) is particularly strong in certain subgroups, and therefore statistical analysis focuses on the interaction between treatment and genotype. However, genome-wide analysis in such scenarios requires careful statistical thought since, in addition to the usual problems of multiple testing, the marker-defined sample sizes, and therefore power, vary across the individual genotypes being evaluated. The variability in power means the usual practice of using a common p-value threshold across tests has difficulties. The reason is that the use of a fixed threshold, with variable power, implies that the costs of type I and type II errors are varying across tests in a manner which is implicit rather than dictated by the analyst. In this paper we discuss this problem and describe an easily implementable solution based on Bayes factors. We pay particular attention to the specification of priors, which is not a straightforward task. The methods are illustrated using data from a randomized controlled clinical trial in which homocysteine levels are compared in individuals receiving low and high doses of folate supplements and across marker subgroups. The method we describe is implemented in the R computing environment with code available from http://faculty.washington.edu/jonno/cv.html.

IL-18 Gene Promoter Region 607C/A Polymorphism in HIV-1 Infected North Indian Population

Several host genetic factors play an important role in susceptibility to human immunodeficiency virus type 1 (HIV-1) infection and in its progression to acquired immune deficiency syndrome (AIDS). The interleukin-18 (IL-18) is a multifunctional proinflammatory cytokine that regulates immune responses and plays a pathogenic role in HIV-1 infection by enhancing viral replication. Single nucleotide polymorphisms (SNPs) in the IL-18 gene promoter region may lead to altered transcriptional activity and IL-18 production, and may account for variation in the risk of HIV-1 infection. We have investigated the association between IL-18 promoter polymorphism −607C>A and HIV-1 infection through a case-control study of 500 patients with HIV-1/AIDS and an equal number of age and sex matched controls in a north Indian population. Genotyping using sequence specific primer-polymerase chain reaction (SSP-PCR) showed a statistically significant reduced risk of HIV-1 infection for the A>A genotype [odds ratio (OR) = 0.57, 95% confidence interval (95% CI) = 0.33–0.98, p = 0.040], but not for the C>A genotype (OR = 0.87, 95% CI = 0.66–1.14, p = 0.321). We concluded that the −607A allele of the IL-18 gene promoter polymorphism may play a protective role against the progression of HIV-1 infection in this population.

Differential susceptibility to prevention: GABAergic, dopaminergic, and multilocus effects

BACKGROUND

Randomized prevention trials provide a unique opportunity to test hypotheses about the interaction of genetic predispositions with contextual processes to create variations in phenotypes over time.

METHODS

Using two longitudinal, randomized prevention trials, molecular genetic and alcohol use outcome data were gathered from more than 900 youths to determine whether prevention program participation would, across 2 years, moderate genetic risk for increased alcohol use conferred by the dopaminergic and GABAergic systems.

RESULTS

We found that (a) variance in dopaminergic (DRD2, DRD4, ANKK1) and GABAergic (GABRG1, GABRA2) genes forecast increases in alcohol use across 2 years, and (b) youths at genetic risk who were assigned to the control condition displayed greater increases in alcohol use across 2 years than did youths at genetic risk who were assigned to the prevention condition or youths without genetic risk who were assigned to either condition.

CONCLUSIONS

This study is unique in combining data from two large prevention trials to test hypotheses regarding genetic main effects and gene × prevention interactions. Focusing on gene systems purported to confer risk for alcohol use and abuse, the study demonstrated that participation in efficacious prevention programs can moderate genetic risk. The results also support the differential susceptibility hypothesis that some youths, for genetic reasons, are more susceptible than others to both positive and negative contextual influences.

Diabetes-associated common genetic variation and its association with GLP-1 concentrations and response to exogenous GLP-1

The mechanisms by which common genetic variation predisposes to type 2 diabetes remain unclear. The disease-associated variants in TCF7L2 (rs7903146) and WFS1 (rs10010131) have been shown to affect response to exogenous glucagon-like peptide 1 (GLP-1), while variants in KCNQ1 (rs151290, rs2237892, and rs2237895) alter endogenous GLP-1 secretion. We set out to validate these observations using a model of GLP-1-induced insulin secretion. We studied healthy individuals using a hyperglycemic clamp and GLP-1 infusion. In addition, we measured active and total GLP-1 in response to an oral challenge in nondiabetic subjects. After genotyping the relevant single nucleotide polymorphisms, generalized linear regression models and repeated-measures ANCOVA models incorporating potential confounders, such as age and BMI, were used to assess the associations, if any, of response with genotype. These variants did not alter GLP-1 concentrations in response to oral intake. No effects on β-cell responsiveness to hyperglycemia and GLP-1 infusion were apparent. Diabetes-associated variation (T allele at rs7903146) in TCF7L2 may impair the ability of hyperglycemia to suppress glucagon (45 ± 2 vs. 47 ± 2 vs. 60 ± 5 ng/L for CC, CT, and TT, respectively, P = 0.02). In nondiabetic subjects, diabetes-associated genetic variation does not alter GLP-1 concentrations after an oral challenge or its effect on insulin secretion.

Allele Frequencies and ther2Measure of Linkage Disequilibrium: Impact on Design and Interpretation of Association Studies

The design and interpretation of genetic association studies depends on the relationship between the genotyped variants and the underlying functional variant, often parameterized as the squared correlation orr2measure of linkage disequilibrium between two loci. While it has long been recognized that placing a constraint on ther2between two loci also places a constraint on the difference in frequencies between the coupled alleles, this constraint has not been quantified. Here, quantification of this severe constraint is presented. For example, forr2≥ .8, the maximum difference in allele frequency is ± .06 which occurs when one locus has allele frequency .5. Forr2≥ .8 and allele frequency at one locus of .1, the maximum difference in allele frequency at the second locus is only ± .02. The impact on the design and interpretation of association studies is discussed.

Multilocus genetic profile for dopamine signaling predicts ventral striatum reactivity

Research integrating neuroimaging and molecular genetics has yielded important insights into how variability in brain chemistry predicts individual differences in brain function, behavior and related risk for psychopathology. However, existing studies have been limited by their focus on the independent effects of single polymorphisms with modest impact on brain chemistry. Here, we explored the effects of five functional polymorphisms affecting dopamine (DA) signaling on reward-related ventral striatum (VS) reactivity, measured with BOLD fMRI, in a sample of 69 Caucasians. We also compiled individual multilocus genetic profile scores reflecting the additive effects of alleles conferring relatively increased DA signaling across the five polymorphic loci: DAT1 9-repeat, DRD4 7-repeat, DRD2 -141C Del, DRD2 Taq1A C (A2), and COMT (158)Met. These multilocus DA profile scores accounted for 10.9% of the inter-individual variability in reward-related VS reactivity. In contrast, none of the individual polymorphisms accounted for significant variability. Our results show that biologically informed multilocus genetic profiles have unique promise as indices of variability in brain chemistry that may yield advances in mapping individual differences in behaviorally relevant brain function. In turn, such genetic profiles may fuel gene-environment interactions research establishing trajectories of risk for psychopathology.

A functional polymorphism in the SCN1A gene does not influence antiepileptic drug responsiveness in Italian patients with focal epilepsy

A splice site variation (c.603-91G>A or rs3812718) in the SCN1A gene has been claimed to influence efficacy and dose requirements of carbamazepine and phenytoin. We investigated the relationship between c.603-91G>A polymorphism and response to antiepileptic drugs (AEDs) in 482 patients with drug-resistant and 401 patients with drug-responsive focal epilepsy. Most commonly used AEDs were carbamazepine and oxcarbazepine. The distribution of c.603-91G>A genotypes was similar among drug-resistant and drug-responsive subjects, both in the entire population and in the groups treated with carbamazepine or oxcarbazepine. There was no association between the c.603-91G>A genotype and dosages of carbamazepine or oxcarbazepine. These findings rule out a major role of the SCN1A polymorphism as a determinant of AED response.

ASIC1a polymorphism is associated with temporal lobe epilepsy

Recent in vitro and in vivo data show that acid-sensing ion channel 1a (ASIC1a) activation enhances neuronal excitability in the hippocampus and neocortex, indicating that ASIC1a might play a role in the generation and maintenance of epileptic seizures. The aim of this study was to investigate association of the ASIC1a gene with temporal lobe epilepsy (TLE) for the first time. Six tag single-nucleotide polymorphisms (SNPs) of the ASIC1a gene were selected and genotyped using polymerase chain reaction-restriction fragment length polymorphism in 560 TLE patients and 401 healthy controls. There was a significant allelic and genotypic association between rs844347:A>C and TLE compared with controls. The rs844347-A allele frequency was 88.1% in the patients and 83.0% in control subjects (OR=1.516, 95% CI 1.142-2.013, p=0.004). Furthermore, the haplotype analysis revealed a significant association with TLE. The results of this study demonstrate for the first time an association between an ASC1a variant allele and TLE in a Han Chinese population.

Hepatocellular carcinoma: novel molecular approaches for diagnosis, prognosis, and therapy

The genomic era is changing the understanding of cancer, although translation of the vast amount of data available into decision-making algorithms is far from reality. Molecular profiling of hepatocellular carcinoma (HCC), the most common cause of death among cirrhotic patients and a fast-growing malignancy in Western countries, is enabling the advancement of novel approaches to disease diagnosis and management. Most HCCs arise on a cirrhotic liver, and predictably, an accurate genomic characterization will allow the identification of procarcinogenic signals amenable to selective targeting by chemopreventive strategies. Molecular diagnosis is currently feasible for small tumors, but it has not yet been formalized by scientific guidelines. Molecular treatment is a reality: Sorafenib confers unprecedented survival benefits in patients at advanced stages. Genomic information from tumor and nontumoral tissue will aid prognosis prediction and facilitate the identification of oncogene addiction loops, providing the opportunity for more personalized medicine.

Graves' disease and gene polymorphism of TNF-α, IL-2, IL-6, IL-12, and IFN-γ

The role of genetic factors in the pathogenesis of Graves' disease (GD) is not clear. The purpose of this study was to investigate the association between single nucleotide polymorphisms in pro-inflammatory cytokine genes and GD in Iranian patients. A case-control hospital-based study was carried out on 107 GD patients and 140 healthy controls. Cytokine typing was performed by polymerase chain reaction with sequence-specific primers (PCR-SSP) assay. The allele and genotype frequencies of the following cytokine genes were determined: TNF-α (-308A/G, -238A/G), IL-2 (-330T/G, +166G/T), IL-6 (-174C/G, A/G nt565), IL-12 (-1188A/C), and IFN-γ (UTR 5644A/T). The following alleles and genotypes were significantly overrepresented in patients: TNF-α -308A allele (P < 0.01) and AA genotype (P < 0.05), IL-2 -330G allele (P < 0.01) and GG genotype (P < 0.01), IL-6 -174C allele (P < 0.01) and CC genotype (P < 0.01), IL-12 -1188C allele (P < 0.01) and CC genotype (P < 0.01), IFN-γ UTR5644T allele (P < 0.01) and TT genotype (P < 0.01). In conclusion, this is the first study to show a significant association between GD and IL-2 -330G, IL-12 -1188C, and IFN-γ UTR 5644T alleles. Our results support the hypothesis that polymorphism in pro-inflammatory cytokines might be involved in predisposition to GD.

Analysis of 19 genes for association with type I diabetes in the Type I Diabetes Genetics Consortium families

In recent years the pace of discovery of genetic associations with type I diabetes (T1D) has accelerated, with the total number of confirmed loci, including the major histocompatibility complex (MHC) region, reaching 43. However, much of the deciphering of the associations at these, and the established T1D loci, has yet to be performed in sufficient numbers of samples or with sufficient markers. Here, 257 single-nucleotide polymorphisms (SNPs) have been genotyped in 19 candidate genes (INS, PTPN22, IL2RA, CTLA4, IFIH1, SUMO4, VDR, PAX4, OAS1, IRS1, IL4, IL4R, IL13, IL12B, CEACAM21, CAPSL, Q7Z4c4(5Q), FOXP3, EFHB) in 2300 affected sib-pair families and tested for association with T1D as part of the Type I Diabetes Genetics Consortium's candidate gene study. The study had approximately 80% power at alpha=0.002 and a minor allele frequency of 0.2 to detect an effect with a relative risk (RR) of 1.20, which drops to just 40% power for a RR of 1.15. At the INS gene, rs689 (-23 HphI) was the most associated SNP (P=3.8 x 10(-31)), with the estimated RR=0.57 (95% confidence interval, 0.52-0.63). In addition, rs689 was associated with age-at-diagnosis of T1D (P=0.001), with homozygosity for the T1D protective T allele, delaying the onset of T1D by approximately 2 years in these families. At PTPN22, rs2476601 (R620W), in agreement with previous reports, was the most significantly associated SNP (P=6.9 x 10(-17)), with RR=1.55 (1.40-1.72). Evidence for association with T1D was observed for the IFIH1 SNP, rs1990760 (P=7.0 x 10(-4)), with RR=0.88 (0.82-0.95) and the CTLA4 SNP rs1427676 (P=0.0005), with RR=1.14 (1.06-1.23). In contrast, no convincing evidence of association was obtained for SUMO4, VDR, PAX4, OAS1, IRS1, IL4, IL4R, IL13, IL12B, CEACAM21 or CAPSL gene regions (http://www.T1DBase.org).

Bayes factors for genome-wide association studies: comparison with P-values

The Bayes factor is a summary measure that provides an alternative to the P-value for the ranking of associations, or the flagging of associations as "significant". We describe an approximate Bayes factor that is straightforward to use and is appropriate when sample sizes are large. We consider various choices of the prior on the effect size, including those that allow effect size to vary with the minor allele frequency (MAF) of the marker. An important contribution is the description of a specific prior that gives identical rankings between Bayes factors and P-values, providing a link between the two approaches, and allowing the implications of the use of P-values to be more easily understood. As a summary measure of noteworthiness P-values are difficult to calibrate since their interpretation depends on MAF and, crucially, on sample size. A consequence is that a consistent decision-making procedure using P-values requires a threshold for significance that reduces with sample size, contrary to common practice.

The medical and economic roles of pipeline pharmacogenetics: Alzheimer's disease as a model of efficacy and HLA-B(*)5701 as a model of safety

Pharmacogenetics (PGX) is the study of drug response as a function of an individual's DNA. PGX is often viewed as an extension of disease association genetics, and although this information may be related, it is not the study of drug response. Although medicines are used to treat diseases, the value of strategies that identify and incorporate DNA biomarkers associated with clinical efficacy, or DNA biomarkers for untoward clinical responses, can be applied directly to pharmaceutical pipelines. The growth of adverse event PGX studies involving marketed medicines generally uses relatively large numbers of affected patients, but has been productive. However, the two critical strategies for pipeline genetics must make use of fewer patients: (1) the early identification of efficacy signals so that they can be applied early in development for targeted therapies and (2) identification of safety signals that can subsequently be validated prospectively during development using the least number of patients with adverse responses. Assumptions are often made that large numbers of patients are necessary to recognize PGX hypotheses and to validate DNA biomarkers. In some ways, pipeline pharmacogenetics may be viewed as the opposite of current genome-wide scanning designs. The goal is to obtain PGX signals in as few patients as possible, and then validate PGX hypotheses for specificity and sensitivity as development trials go forward--not using hundreds of thousand of markers to detect strong linkage disequilibrium signals in thousands of patients and their controls. Drug development takes 5-7 years for a drug candidate to traverse to registration--and this is similar to the timeframe for validating genetic biomarkers using sequential clinical trials. Two important examples are discussed, the association of APOE genotypes to the demonstration of actionable efficacy signals for the use of rosiglitazone for Alzheimer's disease; and the identification of HLA-B(*)5701 as a highly sensitive and specific predictive marker for abacavir treated patients who will develop hypersensitivity syndrome (HSS). The rosiglitazone study prevented pipeline attrition by changing the interpretation of a critical Phase IIB proof of concept study (2005) from a failed study, to a positive efficacy response in a genetically predictable proportion of patients. Now, three years later, a Phase III program of clinical trials using pharmacogenetic designs is months away from completion (late08). If successfully registered (early09), millions of patients could benefit, and efficacy PGX would have achieved its first prospective block-buster. The use of safety candidate gene association genetics in patients who received abacavir therapy and developed HSS starting in 1998 culminated in a double blind clinical trial that determined sensitivity > 97% and specificity >99% in 2007. Clinical consensus panels rapidly recommended abacavir as the preferred therapy along with HLA-B(*)5701 pre-testing, immediately increasing the market share of abacavir with respect to other reverse transcriptases that are associated with there own adverse events. Targeting of medicines during drug development is now possible, practical, and profitable.

Correlation of population parameters leading to power differences in association studies with population stratification

The power of statistical tests to measure effect sizes in the presence of population stratification is an important issue for the design and analysis of population-based association studies. Comparisons of statistical tests have shown that the power of different statistical approaches varies in different genetic scenarios. However, the impact of stratified population parameters on statistical power is not yet understood in a general statistical framework, particularly the impact of correlated population parameters. To investigate such impact in detail, we implemented a genetic model for population-based association studies with stratified samples and evaluated the impact on power with different genetic scenarios. The investigation shows that correlation between disease prevalence and risk allele frequency among subpopulations impacts statistical power. In a model with five subpopulations and moderate population divergence (Fst= 0.01), the correlation accounts for more than 85% of power difference. Our results also show that the estimation of genetic effect for candidate loci is biased by population divergence. Beneficial alleles could be wrongly characterized as risk alleles when prevalence differences and divergences of risk loci are large among subpopulations.

False-positive results in cancer epidemiology: a plea for epistemological modesty

False-positive results are inherent in the scientific process of testing hypotheses concerning the determinants of cancer and other human illnesses. Although much of what is known about the etiology of human cancers has arisen from well-conducted epidemiological studies, epidemiology has been increasingly criticized for producing findings that are often sensationalized in the media and fail to be upheld in subsequent studies. Herein we describe examples from cancer epidemiology of likely false-positive findings and discuss conditions under which such results may occur. We suggest general guidelines or principles, including the endorsement of editorial policies requiring the prominent listing of study caveats, which may help reduce the reporting of misleading results. Increased epistemological humility regarding findings in epidemiology would go a long way to diminishing the detrimental effects of false-positive results on the allocation of limited research resources, on the advancement of knowledge of the causes and prevention of cancer, and on the scientific reputation of epidemiology and would help to prevent oversimplified interpretations of results by the media and the public.

Controlling false discoveries in genetic studies

A false discovery occurs when a researcher concludes that a marker is involved in the etiology of the disease whereas in reality it is not. In genetic studies the risk of false discoveries is very high because only few among the many markers that can be tested will have an effect on the disease. In this article, we argue that it may be best to use methods for controlling false discoveries that would introduce the same ratio of false discoveries divided by all rejected tests into the literature regardless of systematic differences between studies. After a brief discussion of traditional "multiple testing" methods, we show that methods that control the false discovery rate (FDR) may be more suitable to achieve this goal. These FDR methods are therefore discussed in more detail. Instead of merely testing for main effects, it may be important to search for gene-environment/covariate interactions, gene-gene interactions or genetic variants affecting disease subtypes. In the second section, we point out the challenges involved in controlling false discoveries in such searches. The final section discusses the role of replication studies for eliminating false discoveries and the complexities associated with the definition of what constitutes a replication and the design of these studies.

Density of common complex ocular traits in the aging eye: analysis of secondary traits in genome-wide association studies

Genetic association studies are identifying genetic risks for common complex ocular traits such as age-related macular degeneration (AMD). The subjects used for discovery of these loci have been largely from clinic-based, case-control studies. Typically, only the primary phenotype (e.g., AMD) being studied is systematically documented and other complex traits (e.g., affecting the eye) are largely ignored. The purpose of this study was to characterize these other or secondary complex ocular traits present in the cases and controls of clinic-based studies being used for genetic study of AMD. The records of 100 consecutive new patients (of any diagnosis) age 60 or older for which all traits affecting the eye had been recorded systematically were reviewed. The average patient had 3.5 distinct diagnoses. A subset of 10 complex traits was selected for further study because they were common and could be reliably diagnosed. The density of these 10 complex ocular traits increased by 0.017 log-traits/year (P = 0.03), ranging from a predicted 2.74 at age 60 to 4.45 at age 90. Trait-trait association was observed only between AMD and primary vitreomacular traction (P = 0.0009). Only 1% of subjects age 60 or older had no common complex traits affecting the eye. Extrapolations suggested that a study of 2000 similar subjects would have sufficient power to detect genetic association with an odds ratio of 2.0 or less for 4 of these 10 traits. In conclusion, the high prevalence of complex traits affecting the aging eye and the inherent biases in referral patterns leads to the potential for confounding by undocumented secondary traits within case-control studies. In addition to the primary trait, other common ocular phenotypes should be systematically documented in genetic association studies so that adjustments for potential trait-trait associations and other bias can be made and genetic risk variants identified in secondary analyses.

Consortium analysis of 7 candidate SNPs for ovarian cancer

The Ovarian Cancer Association Consortium selected 7 candidate single nucleotide polymorphisms (SNPs), for which there is evidence from previous studies of an association with variation in ovarian cancer or breast cancer risks. The SNPs selected for analysis were F31I (rs2273535) in AURKA, N372H (rs144848) in BRCA2, rs2854344 in intron 17 of RB1, rs2811712 5' flanking CDKN2A, rs523349 in the 3' UTR of SRD5A2, D302H (rs1045485) in CASP8 and L10P (rs1982073) in TGFB1. Fourteen studies genotyped 4,624 invasive epithelial ovarian cancer cases and 8,113 controls of white non-Hispanic origin. A marginally significant association was found for RB1 when all studies were included [ordinal odds ratio (OR) 0.88 (95% confidence interval (CI) 0.79-1.00) p = 0.041 and dominant OR 0.87 (95% CI 0.76-0.98) p = 0.025]; when the studies that originally suggested an association were excluded, the result was suggestive although no longer statistically significant (ordinal OR 0.92, 95% CI 0.79-1.06). This SNP has also been shown to have an association with decreased risk in breast cancer. There was a suggestion of an association for AURKA, when one study that caused significant study heterogeneity was excluded [ordinal OR 1.10 (95% CI 1.01-1.20) p = 0.027; dominant OR 1.12 (95% CI 1.01-1.24) p = 0.03]. The other 5 SNPs in BRCA2, CDKN2A, SRD5A2, CASP8 and TGFB1 showed no association with ovarian cancer risk; given the large sample size, these results can also be considered to be informative. These null results for SNPs identified from relatively large initial studies shows the importance of replicating associations by a consortium approach.

Racial disparity in maternal-fetal genetic epistasis in spontaneous preterm birth

OBJECTIVE

To understand the differences in genetic interactions among tumor necrosis factor-alpha, interleukin-6 and their receptor gene variants between black and white patients in spontaneous preterm birth.

STUDY DESIGN

Maternal and fetal DNA (n = 1195) were collected from cases (preterm birth < 36 weeks' gestation; n = 448), controls (> 37 weeks' gestation; n = 747), and genotyped for single nucleotide polymorphisms in tumor necrosis factor-alpha, tumor necrosis factor receptor 1, and tumor necrosis factor receptor 2, interleukin-6, and interleukin-6 receptor loci. Multifactor dimensionality reduction analysis was used to test all single and multilocus combinations for the ability to predict pregnancy outcome.

RESULTS

In white patients, multilocus interactions in maternal DNA between single nucleotide polymorphisms at -7227 (interleukin-6), 22,215 (interleuki-6 receptor) and -3448 (tumor necrosis factor-alpha) was predictive of approximately 59.1% (P < .02; odds ratio, 2.3 [95% confidence interval = 1.6-3.4]) of pregnancy outcome. In white fetal DNA and black maternal DNA, no significant interactive models were observed. In black patients, the best epistatic model was in fetal DNA between single nucleotide polymorphisms at 17,691 (tumor necrosis factor-receptor 1) and at -3448 (tumor necrosis factor-alpha) and was predictive of pregnancy outcome 68.3% of the time (P < .01; odds ratio, 5.0 [95% confidence interval = 2.6-9.6]).

CONCLUSION

Analyses of multilocus interactions found/associated different models in black and white patients in both maternal and fetal DNA with preterm birth as outcome. Significant maternal-fetal interactions were not detected in either race.

Genetic developments in autoimmune thyroid disease: an evolutionary process

The identification of genes placing individuals at an increased risk for the development of autoimmune thyroid disease (AITD) has been a slow process. However, over the last 20 years or so real progress has been made with the mapping of novel loci, via a number of different approaches. First, through the use of traditional immunological methods, Human Leucocyte Antigen (HLA)/Major Histocompatibility Complex (MHC) was the first gene region to be associated with AITD and consistent replications have been reported. Second, the CTLA-4 gene region on 2q33 was the first non-MHC replicated locus to be primarily identified using the candidate gene method. Third, family-based linkage studies led to the mapping of a new type 1 diabetes locus, the PTPN22 gene, which has subsequently been independently replicated as a susceptibility gene for Graves' disease (GD). Fourth, despite many unsuccessful attempts at implicating the TSHR gene as a susceptibility locus for GD, a recent approach of 'tagging' all the common variation within the gene has led to its identification as the first GD specific locus. Moreover, the use of tag single nucleotide polymorphisms (SNPs) has also been used to implicate the recently identified type 1 diabetes locus, CD25 as a susceptibility gene for GD. Finally, large scale, ongoing genome-wide association studies in multiple autoimmune diseases (AID) states, including AITD seem likely to lead to the identification of additional MHC and non-MHC susceptibility loci.

Methodological Issues in Genetic Association Studies of Inherited Thrombophilia: Original Report of Recent Practice

The aims of this article are to evaluate the methodological quality of genetic association studies on the inherited thrombophilia published during 2003 to 2005, to identify the most common mistakes made by authors of those studies, and to examine if overall quality of the article correlates with the quality of the journal. Articles were evaluated by 2 independent reviewers using the checklist of 16 items. A total of 58 eligible studies were identified. Average total score was 7.59 ± 1.96. Total article score did not correlate with the journal impact factor (r = 0.3971; 95% confidence interval [CI], 0.1547-0.5944, P = .002). Total score did not differ across years ( P = .624). Finally, it is concluded that methodological quality of genetic association studies is not optimal, and it does not depend on the quality of the journal. Journals should adopt methodological criteria for reporting the genetic association studies, and editors should encourage authors to strictly adhere to those criteria.

The candidate genes TAF5L, TCF7, PDCD1, IL6 and ICAM1 cannot be excluded from having effects in type 1 diabetes

BACKGROUND

As genes associated with immune-mediated diseases have an increased prior probability of being associated with other immune-mediated diseases, we tested three such genes, IL23R, IRF5 and CD40, for an association with type 1 diabetes. In addition, we tested seven genes, TAF5L, PDCD1, TCF7, IL12B, IL6, ICAM1 and TBX21, with published marginal or inconsistent evidence of an association with type 1 diabetes.

METHODS

We genotyped reported polymorphisms of the ten genes, nonsynonymous SNPs (nsSNPs) and, for the IL12B and IL6 regions, tag SNPs in up to 7,888 case, 8,858 control and 3,142 parent-child trio samples. In addition, we analysed data from the Wellcome Trust Case Control Consortium genome-wide association study to determine whether there was any further evidence of an association in each gene region.

RESULTS

We found some evidence of associations between type 1 diabetes and TAF5L, PDCD1, TCF7 and IL6 (ORs = 1.05 - 1.13; P = 0.0291 - 4.16 x 10-4). No evidence of an association was obtained for IL12B, IRF5, IL23R, ICAM1, TBX21 and CD40, although there was some evidence of an association (OR = 1.10; P = 0.0257) from the genome-wide association study for the ICAM1 region.

CONCLUSION

We failed to exclude the possibility of some effect in type 1 diabetes for TAF5L, PDCD1, TCF7, IL6 and ICAM1. Additional studies, of these and other candidate genes, employing much larger sample sizes and analysis of additional polymorphisms in each gene and its flanking region will be required to ascertain their contributions to type 1 diabetes susceptibility.

Mechanisms of disease: The genetic basis of coronary heart disease

Since completion of the human genome sequence, considerable progress has been made in determining the genetic basis of human diseases. Understanding the genetic basis of coronary heart disease (CHD), the leading cause of mortality in developed countries, is a priority. Here we provide an update on the genetic basis of CHD, focusing mainly on the clinical manifestations rather than the risk factors, most of which are heritable and also influenced by genetic factors. The challenges faced when identifying clinically relevant genetic determinants of CHD include phenotypic and genetic heterogeneity, and gene-gene and gene-environment interactions. In addition, the etiologic spectrum includes common genetic variants with small effects, as well as rare genetic variants with large effects. Advances such as the cataloging of human genetic variation, new statistical approaches for analyzing massive amounts of genetic data, and the development of high-throughput single-nucleotide polymorphism genotyping platforms, will increase the likelihood of success in the search for genetic determinants of CHD. Such knowledge could refine cardiovascular risk stratification and facilitate the development of new therapies.

Non-replication of association studies: "pseudo-failures" to replicate?

Recently, serious doubts have been cast on the usefulness of association studies as a means to genetically dissect complex diseases because most initial findings fail to replicate in subsequent studies. The reasons usually invoked are population stratification, genetic heterogeneity, and inflated Type I errors. In this article, we argue that, even when these problems are addressed, the scientific community usually has unreasonably high expectations on replication success, based on initial low P values, a phenomenon known as the replication fallacy. We present a modified formula that gives the replication power of a second association study based on the P value of an initial study. When both studies have similar sample sizes, this formula shows that: (1) a P value only slightly lower than the nominal alpha results in only approximately 50% replication power; (2) very low P values are required to achieve a replication power of at least 80% (e.g., at alpha = 0.05, a P value of <0.005 is required). Because many initially significant findings result in low replication power, replication failure should not be surprising or be interpreted as necessarily refuting the initial findings. We refer to replication failures for which the replication power is low as "pseudo-failures."

Study on the association of BoLA-DRB3.2 alleles with clinical mastitis in Norwegian Red cows

Genotyping of bovine leucocyte antigen DRB3.2 (BoLA-DRB3.2) in a total of 523 Norwegian Red (NR) cows from two groups selected for high protein yield and low clinical mastitis, respectively, identified 27 previously reported BoLA-DRB3.2 alleles across the groups. Significant differences in BoLA-DRB3.2 allele frequencies were found between the selection groups. Alleles *13, *18, *22 and *27 had a significantly higher frequency in cows selected for low clinical mastitis, while alleles *3, *9, *11 and *26 had a higher frequency in cows selected for high protein yield. Associations between BoLA-DRB3.2 alleles and clinical mastitis were analysed based on mastitis data from 741,072 first-lactation NR cows, of which 452 were genotyped. Alleles *22 and *26 were found to be associated with increased clinical mastitis, while alleles *7, *11, *18 and *24 had a favourable effect on mastitis resistance. Contradictory results from different studies investigating associations between BoLA-DRB3.2 alleles and mastitis indicate that future studies should focus on associations of mastitis with BoLA haplotypes rather than with single BoLA genes.

Susceptibility genes in Graves' ophthalmopathy: searching for a needle in a haystack?

The variety of clinical presentations of eye changes in patients with Graves' disease suggests that complex interactions between genetic, environmental, endogenous and local factors influence the development/severity of Graves' ophthalmopathy (GO). At present, the role of genetic factors in the development of GO remains unknown. Based on small case-control association studies with candidate genes, several susceptibility loci in GO have been proposed. These are human leucocyte antigen (HLA, 6p21.3), cytotoxic T-lymphocyte antigen-4 (CTLA-4, 2q33), tumour necrosis factor (TNF, 6p21.3), interferon-gamma (IFN-gamma, 12q14), intercellular adhesion molecule-1 (ICAM-1, 19p13), and thyroid stimulating hormone receptor gene (TSH-R, 14q31). Unfortunately, these results were either not confirmed or require replication in larger studies. There are many reasons for the lack of reproducibility of association studies in GO, including poor characterization of the studied groups and small sample sizes, which may result in both false positive and negative results. Thus, the genetic background of GO remains to be elucidated in future research. However, the possibility that GO may be a genetically heterogeneous disorder, or that the development of GO may be predominantly influenced by environmental factors such as cigarette smoking, can not be disregarded.

Isolation of two regions on rat chromosomes 5 and 18 affecting mammary cancer susceptibility

We previously mapped several quantitative trait loci (QTLs) controlling DMBA-induced mammary tumor development in female rats derived from a SPRD-Cu3 (susceptible strain) x WKY (resistant strain) cross. Two of these QTLs were assigned to chromosomes 5 and 18. In the present study, we generated and characterized congenic strains in which a segment of WKY chromosomes 5 or 18 was introduced in the SPRD-Cu3 genetic background, thereby physically demonstrating that each of these two chromosomes controls mammary tumor multiplicity. The chromosome 5 QTL (Mcstm1) accounts for 7 tumors per animal (versus a total of 11 tumors per SPRD-Cu3 rat). The chromosome 18 QTL (Mcstm2) accounts for 3 tumors per animal and is the first chemically-induced mammary cancer susceptibility locus assigned to this chromosome. In addition, the Mcstm1 region was shown to also controls tumor latency. These loci thus play a major role in chemically-induced mammary tumor development. QTLs controlling chemically-induced or estrogen-induced mammary tumor development have independently been identified on chromosomes 5 and 18, using susceptible strains others than SPRD-Cu3. Therefore the haplotype structure of the relevant chromosome regions was analyzed in the different strains. Some chromosome regions were found to be highly mosaic (haplotype blocks < 1 Mb), while one region showed an apparently conserved haplotype block of 7.5 Mb. This analysis points to limited regions that could harbor the causative genes and also indicates that at least Mcstm2 is a novel QTL.

The glu298asp polymorphism in the nitric oxide synthase 3 gene is associated with the risk of ischemic stroke in two large independent case–control studies

The search for genes involved in the pathogenesis of stroke has been highlighted as a field of needs. We followed the concept, that stroke represents a complex genetic disorder, and analyzed the contribution of 106 informative single nucleotide polymorphisms (SNPs) from 63 candidate genes for cardiovascular diseases for the risk of stroke. We conducted two independent case-control studies in two different German regions and recruited a total of 1,901 hospitalized stroke cases and 1,747 regional population controls. The smaller of both studies was used as the replication study. Multiplex PCR in combination with allele-specific hybridization was used for genotype determination. Descriptive statistics, permutations and multivariable logistic regression were used in the analyses. After permutation testing 5 SNPs, located in the nitric oxide synthase 3, the alpha 2 integrin, the interleukin 13, the selectin P and the chemokine receptor 2 genes, had a significant allele difference between cases and controls in the larger study. For one of these SNPs, the glu298asp polymorphism in the nitric oxide synthase 3 gene, an association with ischemic stroke was replicated in the second study and also in a combined analysis of both studies. This association was independent of age, gender, hypertension, diabetes and hypercholesterolemia in both studies. Using large sample sizes and a replication study approach, we found evidence for a role of a polymorphism in the nitric oxide synthase 3 gene in stroke onset.

Interleukin-12B gene polymorphism does not confer susceptibility to Graves' ophthalmopathy in Japanese population

Graves' disease (GD) is an autoimmune disorder with genetic predisposition and frequently associated with Graves' ophthalmopathy (GO). Interleukin 12 (IL-12) is an important mediator of inflammatory immune responses and is expressed in the thyroid and orbit. IL-12B gene, which encodes the p40 subunit of IL-12, is located at chromosome 5q31-33. The aim of the present study was to investigate whether IL-12B gene polymorphism is associated with the development of GD or GO. IL-12B gene polymorphism was studied in Japanese GD patients (n = 329) and healthy control subjects without anti-thyroid autoantibodies or a family history of autoimmune disorders (n = 226). The A/C polymorphism at position 1188 of the 3' untranslated region (3'UTR) of the IL-12B gene was analyzed using the polymerase chain reaction--restriction fragment length polymorphism method. There was no difference in allele or genotype frequency of the IL-12B gene polymorphism (1188A/C) between GD patients and control subjects. There was no association of the IL-12B gene polymorphism with ophthalmopathy, severity of hyperthyroidism or serum IgE levels. There was no association of the IL-12B gene polymorphism with serum IL-12 levels, which were significantly elevated in hyperthyroid phase of GD. In conclusion, IL-12B gene 1188A/C polymorphism is not associated with GD or GO susceptibility in Japanese.

Commonly studied single-nucleotide polymorphisms and breast cancer: results from the Breast Cancer Association Consortium

BACKGROUND

The Breast Cancer Association Consortium (BCAC) is an international collaboration that was established to provide large sample sizes for examining genetic associations. We conducted combined analyses on all single-nucleotide polymorphisms (SNPs) whose associations with breast cancer have been investigated by at least three participating groups.

METHODS

Data from up to 12 studies were pooled for each SNP (ADH1C I350V, AURKA F31I, BRCA2 N372H, CASP8 D302H, ERCC2 D312N, IGFBP3 -202 c>a, LIG4 D501D, PGR V660L, SOD2 V16A, TGFB1 L10P, TP53 R72P, XRCC1 R399Q, XRCC2 R188H, XRCC3 T241M, XRCC3 5' UTR, and XRCC3 IVS7-14). Genotype frequencies in case and control subjects were compared, and genotype-specific odds ratios for the risk of breast cancer in heterozygotes and homozygotes for the rare allele compared with homozygotes for the common allele were estimated with logistic regression. Statistical tests were two-sided.

RESULTS

The total number of subjects for analysis of each SNP ranged from 12,013 to 31,595. For five SNPs--CASP8 D302H, IGFBP3 -202 c>a, PGR V660L, SOD2 V16A, and TGFB1 L10P--the associations with breast cancer were of borderline statistical significance (P = .016, .060, .047, .056, and .0088 respectively). The remaining 11 SNPs were not associated with breast cancer risk; genotype-specific odds ratios were close to unity. There was some evidence for between-study heterogeneity (P<.05) for four of the 11 SNPs (ADH1C I350V, ERCC2 D312N, XRCC1 R399Q, and XRCC3 IVS5-14).

CONCLUSION

Pooling data within a large consortium has helped to clarify associations of SNPs with breast cancer. In the future, consortia such as the BCAC will be important in the analysis of rare polymorphisms and gene x gene or gene x environment interactions, for which individual studies have low power to identify associations, and in the validation of associations identified from genome-wide association studies.

Association of the PTPN22/LYP gene with type 1 diabetes

OBJECTIVES

The goal of this study was to verify the association between type 1 diabetes (T1D) and the protein tyrosine phosphatase, non-receptor type 22 (PTPN22) gene in non-Hispanic whites (NHWs) and Hispanics from Colorado.

SUBJECTS AND METHODS

The C1858T single-nucleotide polymorphism within the PTPN22 gene was genotyped in 753 patients with T1D ascertained from the diabetes clinic at the Barbara Davis Center in Denver and 662 control population.

RESULTS

Both the PTPN22 CT genotype [odds ratio (OR) = 1.96; p < 0.0001] and TT genotype (OR = 4.41; p = 0.02) were significantly associated with T1D in the NHW population. While the association was stronger in subjects with non-HLA-DR3/4 genotypes than in those with the HLA-DR3/4 genotype, regression analyses did not reveal significant interaction between PTPN22 genotypes and HLA-DR3/4. The strength of the association was similar in males and females, patients diagnosed before and after age 10 yr, and in Hispanics and NHWs.

CONCLUSION

In this study, we confirm that PTPN22 is associated with T1D in the Colorado population.

Case-control studies in the genomic era: a clinician's guide

The goal of case-control association studies is to find genetic variants in the human genome that influence common traits. The Human Genome and HapMap projects have added fresh impetus to this goal by cataloguing the raw genetic data behind human DNA variation. Studies that associate these genetic variants with phenotype improve both molecular diagnostics and drug discovery and offer clinicians important opportunities to improve care of patients. In this review I focus on case-control studies, which are the most widely used design and expected to be the most powerful. I also address the problem of case-control non-replication, which is widespread despite enormous effort and use of resources. Important causes of non-replication include inadequate statistical power to detect small and moderate effects, phenotype heterogeneity, population stratification, publication bias, and multiple comparison testing.

Estimation and testing of genotype and haplotype effects in case-control studies: comparison of weighted regression and multiple imputation procedures

A popular approach for testing and estimating genotype and haplotype effects associated with a disease outcome is to conduct a population-based case/control study, in which haplotypes are not directly observed but may be inferred probabilistically from unphased genotype data. A variety of methods exist to analyse the resulting data while accounting for the uncertainty in haplotype assignment, but most focus on the issue of testing the global null hypothesis that no genotype or haplotype effects exist. A more interesting question, once a region of disease association has been identified, is to estimate the relevant genotypic or haplotypic effects and to perform tests of complex null hypotheses such as the hypothesis that some loci, but not others, are associated with disease. Here I examine the assumptions behind, and the performance of, two classes of methods for addressing this question. The first is a weighted regression approach in which posterior probabilities of haplotype assignments are used as weights in a logistic regression analysis, generating a test based on either a weighted pseudo-likelihood, or a weighted log-likelihood. The second is a multiple imputation approach using either an improper procedure in which the posterior probabilities are used to generate replicate imputed data sets, or a proper data augmentation procedure. I compare these approaches to a simple expectation substitution (haplotype trend regression) approach. In simulations, all methods gave unbiased parameter estimation but the weighted pseudo-likelihood, expectation substitution and multiple imputation methods had superior confidence interval coverage. For the weighted pseudo-likelihood and expectation substitution methods it was necessary to estimate posterior haplotype assignment probabilities using the combined case/control data, whereas for the multiple imputation approaches it was necessary to estimate these probabilities in the case and control groups separately. Overall, multiple imputation was easiest approach to implement in standard statistical software and to extend to more complex models such as those that include gene-gene or gene-environment interactions.

Interleukin (IL)-12, IL-2, and IL-6 gene polymorphisms in Takayasu's arteritis from Turkey

Takayasu's arteritis (TA) is a chronic arterial inflammation of unknown etiology involving mainly the aorta and its major branches. Genetic polymorphisms of cytokines are screened as susceptibility factors for TA in Turkey. A total of 94 patients with TA were investigated for the genetic polymorphisms of the interleukin genes IL12, IL2,and IL6 and were compared with 108 healthy control subjects using polymerase chain reaction-sequence-specific primer method. The frequencies of IL12B 1188 C allele (p = 0.03, OR = 1.7) and CC genotype (p = 0.007, OR = 3.7) were both higher in TA patients than in control subjects. TT genotype at IL2-330 (p = 0.006, OR = 2.4) and GG genotype at IL6-174 (p = 0.04, OR = 1.9) were more frequent in TA patients. Lower prevalence of GT genotype at IL2-330 (p = 0.005, OR = 0.4), CG genotype at IL6-174 (p = 0.001, OR = 0.4), and AG genotypes at IL6-598 (p = 0.01, OR = 0.4) were also detected. The polymorphism of IL-12 as well as IL-6 and IL-2 genes may contribute to susceptibility and pathogenesis of TA by altering cytokine production and inducing inflammation.

Statistical false positive or true disease pathway?

Three very recent reports provide convincing statistical evidence (P < 10(-8)), at a genome-wide level, of the association of common polymorphisms with three different common diseases: systemic lupus erythematosus (IRF5), prostate cancer and type 1 diabetes (IFIH1 region). This adds to the trickle--soon to be a flood--of disease association results that are highly unlikely to be false positives. There are other convincing examples in the last 12 months: age-related macular degeneration (CFH), type 1 diabetes (IL2RA, also known as CD25) and type 2 diabetes (TCF7L2). Given 20 years of a literature full of irreproducible results, what has changed?

A systematic analysis of disease-associated variants in the 3' regulatory regions of human protein-coding genes II: the importance of mRNA secondary structure in assessing the functionality of 3' UTR variants

In an attempt both to catalogue 3' regulatory region (3' RR)-mediated disease and to improve our understanding of the structure and function of the 3' RR, we have performed a systematic analysis of disease-associated variants in the 3' RRs of human protein-coding genes. We have previously analysed the variants that have occurred in two specific domains/motifs of the 3' untranslated region (3' UTR) as well as in the 3' flanking region. Here we have focused upon 83 known variants within the upstream sequence (USS; between the translational termination codon and the upstream core polyadenylation signal sequence) of the 3' UTR. To place these variants in their proper context, we first performed a comprehensive survey of known cis-regulatory elements within the USS and the mechanisms by which they effect post-transcriptional gene regulation. Although this survey supports the view that RNA regulatory elements function within the context of specific secondary structures, there are no general rules governing how secondary structure might exert its influence. We have therefore addressed this question by systematically evaluating both functional and non-functional (based upon in vitro reporter gene and/or electrophoretic mobility shift assay data) USS variant-containing sequences against known cis-regulatory motifs within the context of predicted RNA secondary structures. This has allowed us not only to establish a reliable and objective means to perform secondary structure prediction but also to identify consistent patterns of secondary structural change that could potentiate the discrimination of functional USS variants from their non-functional counterparts. The resulting rules were then used to infer potential functionality in the case of some of the remaining functionally uncharacterized USS variants, from their predicted secondary structures. This not only led us to identify further patterns of secondary structural change but also several potential novel cis-regulatory motifs within the 3' UTRs studied.

Toward the 1,000 dollars human genome

Revolutionary new technologies, capable of transforming the economics of sequencing, are providing an unparalleled opportunity to analyze human genetic variation comprehensively at the whole-genome level within a realistic timeframe and at affordable costs. Current estimates suggest that it would cost somewhere in the region of 30 million US dollars to sequence an entire human genome using Sanger-based sequencing, and on one machine it would take about 60 years. Solexa is widely regarded as a company with the necessary disruptive technology to be the first to achieve the ultimate goal of the so-called 1,000 dollars human genome - the conceptual cost-point needed for routine analysis of individual genomes. Solexa's technology is based on completely novel sequencing chemistry capable of sequencing billions of individual DNA molecules simultaneously, a base at a time, to enable highly accurate, low cost analysis of an entire human genome in a single experiment. When applied over a large enough genomic region, these new approaches to resequencing will enable the simultaneous detection and typing of known, as well as unknown, polymorphisms, and will also offer information about patterns of linkage disequilibrium in the population being studied. Technological progress, leading to the advent of single-molecule-based approaches, is beginning to dramatically drive down costs and increase throughput to unprecedented levels, each being several orders of magnitude better than that which is currently available. A new sequencing paradigm based on single molecules will be faster, cheaper and more sensitive, and will permit routine analysis at the whole-genome level.