Global patterns of DNA polymorphism in noncoding regions in human populations

DNA sequence variation is currently a major source of data for studying human origins, evolution, and demographic history, and for detecting linkage association of complex diseases. In this dissertation, I investigated DNA variation in worldwide populations from two ∼10 kb autosomal regions on 22q11.2 (noncoding) and 1q24 (introns). A total of 75 variant sites were found among 128 human sequences in the 22q11.2 region, yielding an estimate of 0.088% for nucleotide diversity (π), and a total of 52 variant sites were found among 122 human sequences in the 1q24 region with an estimated π value of 0.057%. The data from these two regions and a 10 kb noncoding region on Xq13.3 all show a strong excess of low-frequency variants in comparison to that expected from an equilibrium population, indicating a relatively recent population expansion. The effective population sizes estimated from the three regions were 11,000, 12,700, and 8,600, respectively, which are close to the commonly used value of 10,000. In each of the two autosomal regions, the age of the most recent common ancestor (MRCA) was estimated to be older than 1 million years among all the sequences and ∼600,000 years among non-African sequences, providing first evidence from autosomal noncoding or intronic regions for a genetic history of humans much more ancient than the emergence of modern humans. The ancient genetic history of humans indicates no severe bottleneck during the evolution of humans in the last half million years; otherwise, much of the ancient genetic history would have been lost during a severe bottleneck. This study strongly suggests that both the “out of Africa” and the multiregional models are too simple for explaining the evolution of modern humans. A compilation of genome-wide data revealed that nucleotide diversity is highest in autosomal regions, intermediate in X-linked regions, and lowest in Y-linked regions. The data suggest the existence of background selection or selective sweep on Y-linked loci. In general, the nucleotide diversity in humans is low compared to that in chimpanzee and Drosophila populations. ^

Single nucleotide polymorphism (SNP) selection for genotype-phenotype association studies

With hundreds of single nucleotide polymorphisms (SNPs) in a candidate gene and millions of SNPs across the genome, selecting an informative subset of SNPs to maximize the ability to detect genotype-phenotype association is of great interest and importance. In addition, with a large number of SNPs, analytic methods are needed that allow investigators to control the false positive rate resulting from large numbers of SNP genotype-phenotype analyses. This dissertation uses simulated data to explore methods for selecting SNPs for genotype-phenotype association studies. I examined the pattern of linkage disequilibrium (LD) across a candidate gene region and used this pattern to aid in localizing a disease-influencing mutation. The results indicate that the r2 measure of linkage disequilibrium is preferred over the common D′ measure for use in genotype-phenotype association studies. Using step-wise linear regression, the best predictor of the quantitative trait was not usually the single functional mutation. Rather it was a SNP that was in high linkage disequilibrium with the functional mutation. Next, I compared three strategies for selecting SNPs for application to phenotype association studies: based on measures of linkage disequilibrium, based on a measure of haplotype diversity, and random selection. The results demonstrate that SNPs selected based on maximum haplotype diversity are more informative and yield higher power than randomly selected SNPs or SNPs selected based on low pair-wise LD. The data also indicate that for genes with small contribution to the phenotype, it is more prudent for investigators to increase their sample size than to continuously increase the number of SNPs in order to improve statistical power. When typing large numbers of SNPs, researchers are faced with the challenge of utilizing an appropriate statistical method that controls the type I error rate while maintaining adequate power. We show that an empirical genotype based multi-locus global test that uses permutation testing to investigate the null distribution of the maximum test statistic maintains a desired overall type I error rate while not overly sacrificing statistical power. The results also show that when the penetrance model is simple the multi-locus global test does as well or better than the haplotype analysis. However, for more complex models, haplotype analyses offer advantages. The results of this dissertation will be of utility to human geneticists designing large-scale multi-locus genotype-phenotype association studies. ^

An Sp1/Sp3 site polymorphism associated with hypermethylation of the candidate tumor suppressor gene RIL in cancer

Gene silencing due to promoter methylation is an alternative to mutations and deletions, which inactivate tumor suppressor genes (TSG) in cancer. We identified RIL by Methylated CpG Island Amplification technique as a novel aberrantly methylated gene. RIL is expressed in normal tissues and maps to the 5q31 region, frequently deleted in leukemias. We found methylation of RIL in 55/80 (69%) cancer cell lines, with highest methylation in leukemia and colon. We also observed methylation in 46/80 (58%) primary tumors, whereas normal tissues showed substantially lower degrees of methylation. RIL expression was lost in 13/16 cancer cell lines and was restored by demethylating agent. Screening of 38 cell lines and 13 primary cancers by SSCP revealed no mutations in RIL, suggesting that methylation and LOH are the primary inactivation mechanisms. Stable transfection of RIL into colorectal cancer cells resulted in reduction in cell growth, clonogenicity, and increased apoptosis upon UVC treatment, suggesting that RIL is a good candidate TSG. ^ In searching for a cause of RIL hypermethylation, we identified a 12-bp polymorphic sequence around the transcription start site of the gene that creates a long allele containing 3CTC repeat. Evolutionary studies suggested that the long allele appeared late in evolution due to insertion. Using bisulfite sequencing, in cancers heterozygous for RIL, we found that the short allele is 4.4-fold more methylated than the long allele (P = 0.003). EMSA results suggested binding of factor(s) to the inserted region of the long allele, but not to the short. EMSA mutagenesis and competition studies, as well as supershifts using nuclear extracts or recombinant Sp1 strongly indicated that those DNA binding proteins are Sp1 and Sp3. Transient transfections of RIL allele-specific expression constructs showed less than 2-fold differences in luciferase activity, suggesting no major effects of the additional Sp1 site on transcription. However, stable transfection resulted in 3-fold lower levels of transcription from the short allele 60 days post-transfection, consistent with the concept that the polymorphic Sp1 site protects against time-dependent silencing. Thus, an insertional polymorphism in the RIL promoter creates an additional Sp1/Sp3 site, which appears to protect it from silencing and methylation in cancer. ^

G870A polymorphism of cyclin D1 gene influences cancer risk - From epidemiological study to mechanism analysis

Alternate splicing of the cyclin D1 gene gives rise to transcript a and b which encode two protein isoforms cyclin D1a and cyclin D1b. Through testing transcript a and transcript b in a series of human samples, we found that cyclin D1 transcript b is ubiquitously expressed as transcript a but in the lower abundance compared to transcript a. Epidemiological studies have reported that the cyclin D1 gene (CCND1) G870A polymorphism influences the risk for a variety of cancer. In this investigation, we examined the cyclin D1b levels in tumor samples with different genotypes and found that higher levels of cyclin D1b are expressed from the A allele than the G allele. Cyclin D1 is known as a cell cycle regulator facilitating the progression of the cell cycle from G1 to S phase in response to the mitogenic signals. It also interacts with several transcription factors and transcriptional coregulators to modulate their activities. It has been reported that cyclin D1a can substitute for estrogen to activate estrogen receptor α (ERα) mediated transcription and can induce the proliferation of estrogen responsive tissues. However the biological role of cyclin D1b in ERα transcriptional regulation has not been previously explored. In this study, we determined that cyclin D1b antagonizes the action of cyclin D1a on ERα mediated transcription. Cell proliferation assays provided the evidence that cyclin D1b negatively regulates estrogen responsive breast cancer cell growth. Taken together, our findings show that the CCND1 G870A polymorphism is correlated with increased levels of cyclin D1b and that cyclin D1b antagonizes the action of cyclin D1a on ERα mediated transcription providing evidence for the mechanism by which the CCND1 G870A polymorphism may be protective in certain types of breast cancer. ^

THE CYTOPLASMIC TAIL OF MHC CLASS I MOLECULES PLAYS A CRITICAL ROLE IN DENDRITIC CELL-INDUCED T CELL IMMUNITY

The presentation of MHC class I (MHC-I)/peptide complexes by dendritic cells (DCs) is critical for the maintenance of central tolerance to self and for the regulation of cytotoxic T lymphocytes (CTL)-mediated adaptive immune responses against pathogens and cancer cells. Interestingly, several findings have suggested that the cytoplasmic tail of MHC class I plays a functional role in the regulation of CTL immune responses. For example, our previous studies demonstrated that exon 7-deleted MHC-I molecules not only showed extended DC cell surface half-lives but also induced significantly increased CTL responses to viral challange invivo. Although exon 7-deleted variant of MHC-I does not occur naturally in humans, the animal studies prompted us to examine whether exon 7-deleted MHC-I molecules could generate augmented CTL responses in a therapeutic DC-based vaccine setting. To examine the stimulatory capacity of exon 7-deleted MHC-I molecules, we generated a lentivirus-mediated gene transfer system to induce the expression of different MHC-I cytoplasmic tail isoforms in both mouse and human DCs. These DCs were then used as vaccines in a melanoma mouse tumor model and in a human invitro co-culture system. In this thesis, we show that DCs expressing exon 7-deleted MHC-I molecules, stimulated remarkably higher levels of T-cell cytokine production and significantly increased the proliferation of meanoma-specific (Pmel-1) T cells compared with DCs expressing wild type MHC-I. We also demonstrate that, in combination with adoptive transfer of Pmel-1 T-cell, DCs expressing exon 7-deleted Db molecules induced greater anti-tumor responses against established B16 melanoma tumors, significantly extending mouse survival as compared to DCs expressing wild-type Db molecules. Moreover, we also observed that human DCs expressing exon 7-deleted HLA-A2 molecules showed similarly augmented CTL stimulatory ability. Mechanistic studies suggest that exon 7-deleted MHC-I molecules showed impaired lateral membrane movement and extended cell surface half-lives within the DC/T-cell interface, leading to increased spatial availability of MHC-I/peptide complexes for recognition by CD8+ T cells. Collectively, these results suggesr that targeting exon 7 within the cytoplasmic tail of MHC-I molecules in DC vaccines has the potential to enhance CD8+ T cell stimulatory capacity and improve clinical outcomes in patients with cancer or viral infections.

HKT2;2/1, a K+-permeable transporter identified in a salt tolerant rice cultivar through surveys of natural genetic polymorphism

We have investigated OsHKT2;1 natural variation in a collection of 49 cultivars with different levels of salt tolerance and geographical origins. The effect of identified polymorphism on OsHKT2;1 activity was analysed through heterologous expression of variants in Xenopus oocytes. OsHKT2;1 appeared to be a highly conserved protein with only five possible amino acid substitutions that have no substantial effect on functional properties. Our study, however, also identified a new HKT isoform, No-OsHKT2;2/1 in Nona Bokra, a highly salt-tolerant cultivar. No-OsHKT2;2/1 probably originated from a deletion in chromosome 6, producing a chimeric gene. Its 5¢ region corresponds to that of OsHKT2;2, whose full-length sequence is not present in Nipponbare but has been identified in Pokkali, a salt-tolerant rice cultivar. Its 3¢ region corresponds to that of OsHKT2;1. No-OsHKT2;2/1 is essentially expressed in roots and displays a significant level of expression at high Na+ concentrations, in contrast to OsHKT2;1. Expressed in Xenopus oocytes or in Saccharomyces cerevisiae, No-OsHKT2;2/1 exhibited a strong permeability to Na+ and K+, even at high external Na+ concentrations, like OsHKT2;2, and in contrast to OsHKT2;1. Our results suggest that No-OsHKT2;2/1 can contribute to Nona Bokra salt tolerance by enabling root K+ uptake under saline conditions.

Dopamine receptor 4 promoter polymorphism modulates memory and neuronal responses to salience

Animal models and human functional imaging data implicate the dopamine system in mediating enhanced encoding of novel stimuli into human memory. A separate line of investigation suggests an association between a functional polymorphism in the promoter region for the human dopamine 4 receptor gene (DRD4) and sensitivity to novelty. We demonstrate, in two independent samples, that the -521Cmayor queT DRD4 promoter polymorphism determines the magnitude of human memory enhancement for contextually novel, perceptual oddball stimuli in an allele dose-dependent manner. The genotype-dependent memory enhancement conferred by the C allele is associated with increased neuronal responses during successful encoding of perceptual oddballs in the ventral striatum, an effect which is again allele dose-dependent. Furthermore, with repeated presentations of oddball stimuli, this memory advantage decreases, an effect mirrored by adaptation of activation in the hippocampus and substantia nigra/ventral tegmental area in C carriers only. Thus, a dynamic modulation of human memory enhancement for perceptually salient stimuli is associated with activation of a dopaminergic-hippocampal system, which is critically dependent on a functional polymorphism in the DRD4 promoter region.

MCT1 T1470A polymorphism influences adherence to strength training in overweight and obese men following a weight loss program

The monocarboxylate transporter (MCT) family member MCT1 transports lactate into and out of myocytes. Oxidative cells import lactate through MCT1 as a substrate, being the role of MCT1 in glycolysis-derived lactate efflux less clear. MCT1 T1470A polymorphism (rs1049434), which has been related with lactate metabolism and sports specialty 1, 2, could be an influencing factor for exercise adherence. Therefore the aim of this study was to relate the adherence to different training modalities with the T1470A MCT1 polymorphism in overweight and obese men following a weight loss program (WLP).

Computerized polymorphic marker identification: Experimental validation and a predicted human polymorphism catalog

A computational system for the prediction of polymorphic loci directly and efficiently from human genomic sequence was developed and verified. A suite of programs, collectively called pompous (polymorphic marker prediction of ubiquitous simple sequences) detects tandem repeats ranging from dinucleotides up to 250 mers, scores them according to predicted level of polymorphism, and designs appropriate flanking primers for PCR amplification. This approach was validated on an approximately 750-kilobase region of human chromosome 3p21.3, involved in lung and breast carcinoma homozygous deletions. Target DNA from 36 paired B lymphoblastoid and lung cancer lines was amplified and allelotyped for 33 loci predicted by pompous to be variable in repeat size. We found that among those 36 predominately Caucasian individuals 22 of the 33 (67%) predicted loci were polymorphic with an average heterozygosity of 0.42. Allele loss in this region was found in 27/36 (75%) of the tumor lines using these markers. pompous provides the genetic researcher with an additional tool for the rapid and efficient identification of polymorphic markers, and through a World Wide Web site, investigators can use pompous to identify polymorphic markers for their research. A catalog of 13,261 potential polymorphic markers and associated primer sets has been created from the analysis of 141,779,504 base pairs of human genomic sequence in GenBank. This data is available on our Web site (pompous.swmed.edu) and will be updated periodically as GenBank is expanded and algorithm accuracy is improved.

Definition of MHC and T cell receptor contacts in the HLA-DR4restricted immunodominant epitope in type II collagen and characterization of collagen-induced arthritis in HLA-DR4 and human CD4 transgenic mice

Rheumatoid arthritis (RA) is an autoimmune disease associated with the HLA-DR4 and DR1 alleles. The target autoantigen(s) in RA is unknown, but type II collagen (CII) is a candidate, and the DR4- and DR1-restricted immunodominant T cell epitope in this protein corresponds to amino acids 261–273 (CII 261–273). We have defined MHC and T cell receptor contacts in CII 261–273 and provide strong evidence that this peptide corresponds to the peptide binding specificity previously found for RA-associated DR molecules. Moreover, we demonstrate that HLA-DR4 and human CD4 transgenic mice homozygous for the I-Abβ0 mutation are highly susceptible to collagen-induced arthritis and describe the clinical course and histopathological changes in the affected joints.

Major histocompatibility complex (MHC) class I KbDb −/− deficient mice possess functional CD8+ T cells and natural killer cells

We obtained mice deficient for major histocompatibility complex (MHC) molecules encoded by the H-2K and H-2D genes. H-2 KbDb −/− mice express no detectable classical MHC class I-region associated (Ia) heavy chains, although β2-microglobulin and the nonclassical class Ib proteins examined are expressed normally. KbDb −/− mice have greatly reduced numbers of mature CD8+ T cells, indicating that selection of the vast majority (>90%) of CD8+ T cells cannot be compensated for by β2-microglobulin-associated molecules other than classical H-2K and D locus products. In accord with the greatly reduced number of CD8+ T cells, spleen cells from KbDb −/− mice do not generate cytotoxic responses in primary mixed-lymphocyte cultures against MHC-disparate (allogeneic) cells. However, in vivo priming of KbDb −/− mice with allogeneic cells resulted in strong CD8+ MHC class Ia-specific allogeneic responses. Thus, a minor population of functionally competent peripheral CD8+ T cells capable of strong cytotoxic activity arises in the complete absence of classical MHC class Ia molecules. KbDb −/− animals also have natural killer cells that retain their cytotoxic potential.

A natural polymorphism in β-lactamase is a global suppressor

A M182T substitution was discovered as a second-site suppressor of a missense mutation in TEM-1 β-lactamase. The combination of the M182T substitution with other substitutions in the enzyme indicates the M182T substitution is a global suppressor of missense mutations in β-lactamase. The M182T substitution also is found in natural variants of TEM-1 β-lactamase with altered substrate specificity that have evolved in response to antibiotic therapy. The M182T substitution may have been selected in natural isolates as a suppressor of folding or stability defects resulting from mutations associated with drug resistance. This pathway of protein evolution may occur in other targets of antimicrobial drugs such as the HIV protease.

Molecular keys to speciation: DNA polymorphism and the control of genetic exchange in enterobacteria

Speciation involves the establishment of genetic barriers between closely related organisms. The extent of genetic recombination is a key determinant and a measure of genetic isolation. The results reported here reveal that genetic barriers can be established, eliminated, or modified by manipulating two systems which control genetic recombination, SOS and mismatch repair. The extent of genetic isolation between enterobacteria is a simple mathematical function of DNA sequence divergence. The function does not depend on hybrid DNA stability, but rather on the number of blocks of sequences identical in the two mating partners and sufficiently large to allow the initiation of recombination. Further, there is no obvious discontinuity in the function that could be used to define a level of divergence for distinguishing species.

Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination

One-third of humans are infected with Mycobacterium tuberculosis, the causative agent of tuberculosis. Sequence analysis of two megabases in 26 structural genes or loci in strains recovered globally discovered a striking reduction of silent nucleotide substitutions compared with other human bacterial pathogens. The lack of neutral mutations in structural genes indicates that M. tuberculosis is evolutionarily young and has recently spread globally. Species diversity is largely caused by rapidly evolving insertion sequences, which means that mobile element movement is a fundamental process generating genomic variation in this pathogen. Three genetic groups of M. tuberculosis were identified based on two polymorphisms that occur at high frequency in the genes encoding catalase-peroxidase and the A subunit of gyrase. Group 1 organisms are evolutionarily old and allied with M. bovis, the cause of bovine tuberculosis. A subset of several distinct insertion sequence IS6110 subtypes of this genetic group have IS6110 integrated at the identical chromosomal insertion site, located between dnaA and dnaN in the region containing the origin of replication. Remarkably, study of ≈6,000 isolates from patients in Houston and the New York City area discovered that 47 of 48 relatively large case clusters were caused by genotypic group 1 and 2 but not group 3 organisms. The observation that the newly emergent group 3 organisms are associated with sporadic rather than clustered cases suggests that the pathogen is evolving toward a state of reduced transmissability or virulence.