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. ^

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.

Resting state functional connectivity patterns associated with ApoE genotype. A MEG study in Mild Cognitive Impairment (MCI)

Alteration of brain communication due to abnormal patterns of synchronization is nowadays one of the most suitable mechanisms for having a better understanding of brain pathologies. Very recently, it has been proved that abnormal changes in both local and long range functional interactions underlie the cognitive deficits associated with different brain disorders. Mild cognitive impairment (MCI) is a state characterized for cognitive dysfunction, such as the memory. The study of the spatial and dynamic alterations in MCI subjects' functional networks could provide important evidences of the brain mechanisms responsible for such impairment.

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.

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.

Role of lipid polymorphism in G protein-membrane interactions: Nonlamellar-prone phospholipids and peripheral protein binding to membranes

Heterotrimeric G proteins (peripheral proteins) conduct signals from membrane receptors (integral proteins) to regulatory proteins localized to various cellular compartments. They are in excess over any G protein-coupled receptor type on the cell membrane, which is necessary for signal amplification. These facts account for the large number of G protein molecules bound to membrane lipids. Thus, the protein-lipid interactions are crucial for their cellular localization, and consequently for signal transduction. In this work, the binding of G protein subunits to model membranes (liposomes), formed with defined membrane lipids, has been studied. It is shown that although G protein α-subunits were able to bind to lipid bilayers, the presence of nonlamellar-prone phospholipids (phosphatidylethanolamines) enhanced their binding to model membranes. This mechanism also appears to be used by other (structurally and functionally unrelated) peripheral proteins, such as protein kinase C and the insect protein apolipophorin III, indicating that it could constitute a general mode of protein-lipid interactions, relevant in the activity and translocation of some peripheral (amphitropic) proteins from soluble to particulate compartments. Other factors, such as the presence of cholesterol or the vesicle surface charge, also modulated the binding of the G protein subunits to lipid bilayers. Conversely, the binding of G protein-coupled receptor kinase 2 and the G protein β-subunit to liposomes was not increased by hexagonally prone lipids. Their distinct interactions with membrane lipids may, in part, explain the different cellular localizations of all of these proteins during the signaling process.

A CDKN2-like polymorphism in Xiphophorus LG V is associated with UV-B-induced melanoma formation in platyfish–swordtail hybrids

The genetic basis of spontaneous melanoma formation in spotted dorsal (Sd) Xiphophorus platyfish–swordtail hybrids has been studied for decades, and is adequately explained by a two-gene inheritance model involving a sex-linked oncogene, Xmrk, and an autosomal tumor suppressor, DIFF. The Xmrk oncogene encodes a receptor tyrosine kinase related to EGFR; the nature of the DIFF tumor suppressor gene is unknown. We analyzed the genetic basis of UV-B-induced melanoma formation in closely related, spotted side platyfish–swordtail hybrids, which carry a different sex-linked pigment pattern locus, Sp. We UV-irradiated spotted side Xiphophorus platyfish–swordtail backcross hybrids to induce melanomas at frequencies 6-fold higher than occur spontaneously in unirradiated control animals. To identify genetic determinants of melanoma susceptibility in this UV-inducible Xiphophorus model, we genotyped individual animals from control and UV-irradiated experimental regimes using allozyme and DNA restriction fragment length polymorphisms and tested for joint segregation of genetic markers with pigmentation phenotype and UV-induced melanoma formation. Joint segregation results show linkage of a CDKN2-like DNA polymorphism with UV-B-induced melanoma formation in these hybrids. The CDKN2-like polymorphism maps to Xiphophorus linkage group V and exhibits recombination fractions with ES1 and MDH2 allozyme markers consistent with previous localization of the DIFF tumor suppressor locus. Our results indicate that the CDKN2-like sequence we have cloned and mapped is a candidate for the DIFF tumor suppressor gene.

RNA-controlled polymorphism in the in vivo assembly of 180-subunit and 120-subunit virions from a single capsid protein

Repeated, specific interactions between capsid protein (CP) subunits direct virus capsid assembly and exemplify regulated protein–protein interactions. The results presented here reveal a striking in vivo switch in CP assembly. Using cryoelectron microscopy, three-dimensional image reconstruction, and molecular modeling, we show that brome mosaic virus (BMV) CP can assemble in vivo two remarkably distinct capsids that selectively package BMV-derived RNAs in the absence of BMV RNA replication: a 180-subunit capsid indistinguishable from virions produced in natural infections and a previously unobserved BMV capsid type with 120 subunits arranged as 60 CP dimers. Each such dimer contains two CPs in distinct, nonequivalent environments, in contrast to the quasi-equivalent CP environments throughout the 180-subunit capsid. This 120-subunit capsid utilizes most of the CP interactions of the 180-subunit capsid plus nonequivalent CP–CP interactions. Thus, the CP of BMV, and perhaps other viruses, can encode CP–CP interactions that are not apparent from mature virions and may function in assembly or disassembly. Shared structural features suggest that the 120- and 180-subunit capsids share assembly steps and that a common pentamer of CP dimers may be an important assembly intermediate. The ability of a single CP to switch between distinct capsids by means of alternate interactions also implies reduced evolutionary barriers between different capsid structures. The in vivo switch between alternate BMV capsids is controlled by the RNA packaged: a natural BMV genomic RNA was packaged in 180-subunit capsids, whereas an engineered mRNA containing only the BMV CP gene was packaged in 120-subunit capsids. RNA features can thus direct the assembly of a ribonucleoprotein complex between alternate structural pathways.