Thomas J. White

A Large-Scale Genetic Association Study Confirms IL12B and Leads to the Identification of IL23R as Psoriasis-Risk Genes

We performed a multitiered, case-control association study of psoriasis in three independent sample sets of white North American individuals (1,446 cases and 1,432 controls) with 25,215 genecentric single-nucleotide polymorphisms (SNPs) and found a highly significant association with an IL12B 3-untranslated-region SNP (rs3212227), confirming the results of a small Japanese study. This SNP was significant in all three sample sets (odds ratio [OR](common) 0.64, combined P [Pcomb]=7.85x10(-10)). A Monte Carlo simulation to address multiple testing suggests that this association is not a type I error. The coding regions of IL12B were resequenced in 96 individuals with psoriasis, and 30 additional IL12B-region SNPs were genotyped. Haplotypes were estimated, and genotype-conditioned analyses identified a second risk allele (rs6887695) located approximately 60 kb upstream of the IL12B coding region that exhibited association with psoriasis after adjustment for rs3212227. Together, these two SNPs mark a common IL12B risk haplotype (OR(common) 1.40, Pcomb=8.11x10(-9)) and a less frequent protective haplotype (OR(common) 0.58, Pcomb=5.65x10(-12)), which were statistically significant in all three studies. Since IL12B encodes the common IL-12p40 subunit of IL-12 and IL-23, we individually genotyped 17 SNPs in the genes encoding the other chains of these cytokines (IL12A and IL23A) and their receptors (IL12RB1, IL12RB2, and IL23R). Haplotype analyses identified two IL23R missense SNPs that together mark a common psoriasis-associated haplotype in all three studies (OR(common) 1.44, Pcomb=3.13x10(-6)). Individuals homozygous for both the IL12B and the IL23R predisposing haplotypes have an increased risk of disease (OR(common) 1.66, Pcomb=1.33x10(-8)). These data, and the previous observation that administration of an antibody specific for the IL-12p40 subunit to patients with psoriasis is highly efficacious, suggest that these genes play a fundamental role in psoriasis pathogenesis.

A Scan for Positively Selected Genes in the Genomes of Humans and Chimpanzees

Since the divergence of humans and chimpanzees about 5 million years ago, these species have undergone a remarkable evolution with drastic divergence in anatomy and cognitive abilities. At the molecular level, despite the small overall magnitude of DNA sequence divergence, we might expect such evolutionary changes to leave a noticeable signature throughout the genome. We here compare 13,731 annotated genes from humans to their chimpanzee orthologs to identify genes that show evidence of positive selection. Many of the genes that present a signature of positive selection tend to be involved in sensory perception or immune defenses. However, the group of genes that show the strongest evidence for positive selection also includes a surprising number of genes involved in tumor suppression and apoptosis, and of genes involved in spermatogenesis. We hypothesize that positive selection in some of these genes may be driven by genomic conflict due to apoptosis during spermatogenesis. Genes with maximal expression in the brain show little or no evidence for positive selection, while genes with maximal expression in the testis tend to be enriched with positively selected genes. Genes on the X chromosome also tend to show an elevated tendency for positive selection. We also present polymorphism data from 20 Caucasian Americans and 19 African Americans for the 50 annotated genes showing the strongest evidence for positive selection. The polymorphism analysis further supports the presence of positive selection in these genes by showing an excess of high-frequency derived nonsynonymous mutations.

Natural selection on protein-coding genes in the human genome

Comparisons of DNA polymorphism within species to divergence between species enables the discovery of molecular adaptation in evolutionarily constrained genes as well as the differentiation of weak from strong purifying selection. The extent to which weak negative and positive darwinian selection have driven the molecular evolution of different species varies greatly, with some species, such as Drosophila melanogaster, showing strong evidence of pervasive positive selection, and others, such as the selfing weed Arabidopsis thaliana, showing an excess of deleterious variation within local populations. Here we contrast patterns of coding sequence polymorphism identified by direct sequencing of 39 humans for over 11,000 genes to divergence between humans and chimpanzees, and find strong evidence that natural selection has shaped the recent molecular evolution of our species. Our analysis discovered 304 (9.0%) out of 3,377 potentially informative loci showing evidence of rapid amino acid evolution. Furthermore, 813 (13.5%) out of 6,033 potentially informative loci show a paucity of amino acid differences between humans and chimpanzees, indicating weak negative selection and/or balancing selection operating on mutations at these loci. We find that the distribution of negatively and positively selected genes varies greatly among biological processes and molecular functions, and that some classes, such as transcription factors, show an excess of rapidly evolving genes, whereas others, such as cytoskeletal proteins, show an excess of genes with extensive amino acid polymorphism within humans and yet little amino acid divergence between humans and chimpanzees.

Assessing the evolutionary impact of amino acid mutations in the human genome.

Quantifying the distribution of fitness effects among newly arising mutations in the human genome is key to resolving important debates in medical and evolutionary genetics. Here, we present a method for inferring this distribution using Single Nucleotide Polymorphism (SNP) data from a population with non-stationary demographic history (such as that of modern humans). Application of our method to 47,576 coding SNPs found by direct resequencing of 11,404 protein coding-genes in 35 individuals (20 European Americans and 15 African Americans) allows us to assess the relative contribution of demographic and selective effects to patterning amino acid variation in the human genome. We find evidence of an ancient population expansion in the sample with African ancestry and a relatively recent bottleneck in the sample with European ancestry. After accounting for these demographic effects, we find strong evidence for great variability in the selective effects of new amino acid replacing mutations. In both populations, the patterns of variation are consistent with a leptokurtic distribution of selection coefficients (e.g., gamma or log-normal) peaked near neutrality. Specifically, we predict 27–29% of amino acid changing (nonsynonymous) mutations are neutral or nearly neutral (|s|<0.01%), 30–42% are moderately deleterious (0.01%<|s|<1%), and nearly all the remainder are highly deleterious or lethal (|s|>1%). Our results are consistent with 10–20% of amino acid differences between humans and chimpanzees having been fixed by positive selection with the remainder of differences being neutral or nearly neutral. Our analysis also predicts that many of the alleles identified via whole-genome association mapping may be selectively neutral or (formerly) positively selected, implying that deleterious genetic variation affecting disease phenotype may be missed by this widely used approach for mapping genes underlying complex traits.

Proportionally more deleterious genetic variation in European than in African populations

Quantifying the number of deleterious mutations per diploid human genome is of crucial concern to both evolutionary and medical geneticists. Here we combine genome-wide polymorphism data from PCR-based exon resequencing, comparative genomic data across mammalian species, and protein structure predictions to estimate the number of functionally consequential single-nucleotide polymorphisms (SNPs) carried by each of 15 African American (AA) and 20 European American (EA) individuals. We find that AAs show significantly higher levels of nucleotide heterozygosity than do EAs for all categories of functional SNPs considered, including synonymous, non-synonymous, predicted ‘benign’, predicted ‘possibly damaging’ and predicted ‘probably damaging’ SNPs. This result is wholly consistent with previous work showing higher overall levels of nucleotide variation in African populations than in Europeans. EA individuals, in contrast, have significantly more genotypes homozygous for the derived allele at synonymous and non-synonymous SNPs and for the damaging allele at ‘probably damaging’ SNPs than AAs do. For SNPs segregating only in one population or the other, the proportion of non-synonymous SNPs is significantly higher in the EA sample (55.4%) than in the AA sample (47.0%; Pu2009<u20092.3u2009×u200910-37). We observe a similar proportional excess of SNPs that are inferred to be ‘probably damaging’ (15.9% in EA; 12.1% in AA; Pu2009<u20093.3u2009×u200910-11). Using extensive simulations, we show that this excess proportion of segregating damaging alleles in Europeans is probably a consequence of a bottleneck that Europeans experienced at about the time of the migration out of Africa.

A 7 gene signature identifies the risk of developing cirrhosis in patients with chronic hepatitis C

Clinical factors such as age, gender, alcohol use, and age‐at‐infection influence the progression to cirrhosis but cannot accurately predict the risk of developing cirrhosis in patients with chronic hepatitis C (CHC). The aim of this study was to develop a predictive signature for cirrhosis in Caucasian patients. All patients had well‐characterized liver histology and clinical factors; DNA was extracted from whole blood for genotyping. We validated all significant markers from a genome scan in the training cohort, and selected 361 markers for the signature building. Using a “machine learning” approach, a signature consisting of markers most predictive for cirrhosis risk in Caucasian patients was developed in the training set (N = 420). The Cirrhosis Risk Score (CRS) was calculated to estimate the risk of developing cirrhosis for each patient. The CRS performance was then tested in an independently enrolled validation cohort of 154 Caucasian patients. A CRS signature consisting of 7 markers was developed for Caucasian patients. The area‐under‐the‐ROC curves (AUC) of the CRS was 0.75 in the training cohort. In the validation cohort, AUC was only 0.53 for clinical factors, increased to 0.73 for CRS, and 0.76 when CRS and clinical factors were combined. A low CRS cutoff of <0.50 to identify low‐risk patients would misclassify only 10.3% of high‐risk patients, while a high cutoff of >0.70 to identify high‐risk patients would misclassify 22.3% of low‐risk patients. Conclusion: CRS is a better predictor than clinical factors in differentiating high‐risk versus low‐risk for cirrhosis in Caucasian CHC patients. Prospective studies should be conducted to further validate these findings. (HEPATOLOGY 2007.)

Phylogeography of the fungal pathogen Histoplasma capsulatum.

Until recently, Histoplasma capsulatum was believed to harbour three varieties, var. capsulatum (chiefly a New World human pathogen), var. duboisii (an African human pathogen) and var. farciminosum (an Old World horse pathogen), which varied in clinical manifestations and geographical distribution. We analysed the phylogenetic relationships of 137 individuals representing the three varieties from six continents using DNA sequence variation in four independent protein‐coding genes. At least eight clades were idengified: (i) North American class 1 clade; (ii) North American class 2 clade; (iii) Latin American group A clade; (iv) Latin American group B clade; (v) Australian clade; (vi) Netherlands (Indonesian?) clade; (vii) Eurasian clade and (viii) African clade. Seven of eight clades represented genetically isolated groups that may be recognized as phylogenetic species. The sole exception was the Eurasian clade which originated from within the Latin American group A clade. The phylogenetic relationships among the clades made a star phylogeny. Histoplasma capsulatum var. capsulatum individuals were found in all eight clades. The African clade included all of the H. capsulatum var. duboisii individuals as well as individuals of the other two varieties. The 13 individuals of var. farciminosum were distributed among three phylogenetic species. These findings suggest that the three varieties of Histoplasma are phylogenetically meaningless. Instead we have to recognize the existence of genetically distinct geographical populations or phylogenetic species. Combining DNA substitution rates of protein‐coding genes with the phylogeny suggests that the radiation of Histoplasma started between 3 and 13 million years ago in Latin America.

Targets of Balancing Selection in the Human Genome

Balancing selection is potentially an important biological force for maintaining advantageous genetic diversity in populations, including variation that is responsible for long-term adaptation to the environment. By serving as a means to maintain genetic variation, it may be particularly relevant to maintaining phenotypic variation in natural populations. Nevertheless, its prevalence and specific targets in the human genome remain largely unknown. We have analyzed the patterns of diversity and divergence of 13,400 genes in two human populations using an unbiased single-nucleotide polymorphism data set, a genome-wide approach, and a method that incorporates demography in neutrality tests. We identified an unbiased catalog of genes with signatures of long-term balancing selection, which includes immunity genes as well as genes encoding keratins and membrane channels; the catalog also shows enrichment in functional categories involved in cellular structure. Patterns are mostly concordant in the two populations, with a small fraction of genes showing population-specific signatures of selection. Power considerations indicate that our findings represent a subset of all targets in the genome, suggesting that although balancing selection may not have an obvious impact on a large proportion of human genes, it is a key force affecting the evolution of a number of genes in humans.

Association of the Trp719Arg Polymorphism in Kinesin-Like Protein 6 With Myocardial Infarction and Coronary Heart Disease in 2 Prospective Trials The CARE and WOSCOPS Trials

OBJECTIVESnWe asked whether 35 genetic polymorphisms, previously found to be associated with cardiovascular disease, were associated with myocardial infarction (MI) in the CARE (Cholesterol and Recurrent Events) trial and with coronary heart disease (CHD) in the WOSCOPS (West of Scotland Coronary Prevention Study) trial and whether the risk associated with these polymorphisms could be reduced by pravastatin treatment.nnnBACKGROUNDnIdentification of genetic polymorphisms associated with CHD may improve assessment of CHD risk and understanding of disease pathophysiology.nnnMETHODSnWe tested the association between genotype and recurrent MI in the CARE study and between genotype and primary CHD in the WOSCOPS trial using regression models that adjusted for conventional risk factors: Cox proportional hazards models for the CARE study and conditional logistic regression models for a nested case-control study of the WOSCOPS trial.nnnRESULTSnWe found that Trp719Arg (rs20455) in KIF6 was associated with coronary events. KIF6 encodes kinesin-like protein 6, a member of the molecular motor superfamily. In placebo-treated patients, carriers of the KIF6 719Arg allele (59.4% of the CARE trial cohort) had a hazard ratio of 1.50 (95% confidence interval [CI] 1.05 to 2.15) in the CARE trial and an odds ratio of 1.55 (95% CI 1.14 to 2.09) in the WOSCOPS trial. Among carriers, the absolute risk reduction by pravastatin was 4.89% (95% CI 1.81% to 7.97%) in the CARE trial and 5.49% (95% CI 3.52% to 7.46%) in the WOSCOPS trial.nnnCONCLUSIONSnIn both the CARE and the WOSCOPS trials, carriers of the KIF6 719Arg allele had an increased risk of coronary events, and pravastatin treatment substantially reduced that risk.

Estrogen Receptor Genotypes and Haplotypes Associated with Breast Cancer Risk

Nearly one in eight US women will develop breast cancer in their lifetime. Most breast cancer is not associated with a hereditary syndrome, occurs in postmenopausal women, and is estrogen and progesterone receptor-positive. Estrogen exposure is an epidemiologic risk factor for breast cancer and estrogen is a potent mammary mitogen. We studied single nucleotide polymorphisms (SNPs) in estrogen receptors in 615 healthy subjects and 1011 individuals with histologically confirmed breast cancer, all from New York City. We analyzed 13 SNPs in the progesterone receptor gene (PGR), 17 SNPs in estrogen receptor 1 gene (ESR1), and 8 SNPs in the estrogen receptor 2 gene (ESR2). We observed three common haplotypes in ESR1 that were associated with a decreased risk for breast cancer [odds ratio (OR), ∼ O.4; 95% confidence interval (CI), 0.2–0.8; P < 0.01]. Another haplotype was associated with an increased risk of breast cancer (OR, 2.1; 95% CI, 1.2–3.8; P < 0.05). A unique risk haplotype was present in ∼7% of older Ashkenazi Jewish study subjects (OR, 1.7; 95% CI, 1.2–2.4; P < 0.003). We narrowed the ESR1 risk haplotypes to the promoter region and first exon. We define several other haplotypes in Ashkenazi Jews in both ESR1 and ESR2 that may elevate susceptibility to breast cancer. In contrast, we found no association between any PGR variant or haplotype and breast cancer. Genetic epidemiology study replication and functional assays of the haplotypes should permit a better understanding of the role of steroid receptor genetic variants and breast cancer risk.