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

Elucidation of the role of 53BP1 in DNA double strand break (DSB) repair and tumor suppression

The protein p53 binding protein one (53BP1) was discovered in a yeast two-hybrid screen that used the DNA binding domain of p53 as bait. Cloning of full-length 53BP1 showed that this protein contains several protein domains which help make up the protein, which include two tandem BRCT domains and a amino-terminal serine/glutamine cluster domain (SCD). These are two protein domains are often seen in factors that are involved in the cellular response to DNA damage and control of cell cycle checkpoints and we hypothesize that 53BP1 is involved in the cellular response to DNA damage. In support of this hypothesis we observe that 53BP1 is phosphorylated and undergoes a dramatic nuclear re-localization in response to DNA damaging agents. 53BP1 also interacts with several factors that are important in the cellular response to DNA damage, such as the BRCA1 tumor suppressor, ATM and Rad3 related (ATR), and the phosphorylated version of the histone variant H2AX. Mice deficient in 53BP1 display increased sensitivity ionizing radiation (IR), a DNA damaging agent that introduces DNA double strand breaks (DSBs). In addition, 53BP1-deficient mice do not properly undergo the process of class switch recombination (CSR). We also observe that when a defect in 53BP1 is combined with a defect in p53; the resulting mice have an increased rate of formation of spontaneous tumors, notably the formation of B and T lineage lymphomas. The T lineage tumors arise by two distinct mechanisms: one driven by defects in cell cycle regulation and a second driven by defects in the ability to repair DNA DSBs. The B lineage tumors arise by the inability to repair DNA damage and over-expression of the oncogene c-myc. ^ With these observations, we conclude that not only does 53BP1 function in the cellular response to DNA damage, but it also works in concert with p53 to suppress tumor formation. ^

Emerging roles for the double strand break repair protein 53BP1 in transcriptional regulation

Disruption of the mechanisms that regulate cell-cycle checkpoints, DNA repair, and apoptosis results in genomic instability and often leads to the development of cancer. In response to double stranded breaks (DSBs) as induced by ionizing radiation (IR), generated during DNA replication, or through immunoglobulin heavy chain (IgH) rearrangements in T and B cells of lymphoid origin, the protein kinases ATM and ATR are central players that activate signaling pathways leading to DSB repair. p53 binding protein 1 (53BP1) participates in the repair of DNA double stranded breaks (DSBs) where it is recruited to or near sites of DNA damage. In addition to its well established role in DSB repair, multiple lines of evidence implicate 53BP1 in transcription which stem from its initial discovery as a p53 binding protein in a yeast two-hybrid screen. However, the mechanisms behind the role of 53BP1 in these processes are not well understood. ^ 53BP1 possesses several motifs that are likely important for its role in DSB repair including two BRCA1 C-terminal repeats, tandem Tudor domains, and a variety of phosphorylation sites. In addition to these motifs, we identified a glycine and arginine rich region (GAR) upstream of the Tudor domains, a sequence that is oftentimes serves as a site for protein arginine methylation. The focus of this project was to characterize the methylation of 53BP1 and to evaluate how methylation influenced the role of 53BP1 as a tumor suppressor. ^ Using a variety of biochemical techniques, we demonstrated that 53BP1 is methylated by the PRMT1 methyltransferase in vivo. Moreover, GAR methylation occurs on arginine residues in an asymmetric manner. We further show that sequences upstream of the Tudor domains that do not include the GAR stretch are sufficient for 53BP1 oligomerization in vivo. While investigating the role of arginine methylation in 53BP1 function, we discovered that 53BP1 associates with proteins of the general transcription apparatus as well as to other factors implicated in coordinating transcription with chromatin function. Collectively, these data support a role for 53BP1 in regulating transcription and provide insight into the possible mechanisms by which this occurs. ^

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