Estimation of male-to-female ratios of mutation rate in primates and rodents

Haldane (1935) developed a method for estimating the male-to-female ratio of mutation rate ($\alpha$) by using sex-linked recessive genetic disease, but in six different studies using hemophilia A data the estimates of $\alpha$ varied from 1.2 to 29.3. Direct genomic sequencing is a better approach, but it is laborious and not readily applicable to non-human organisms. To study the sex ratios of mutation rate in various mammals, I used an indirect method proposed by Miyata et al. (1987). This method takes advantage of the fact that different chromosomes segregate differently between males and females, and uses the ratios of mutation rate in sequences on different chromosomes to estimate the male-to-female ratio of mutation rate. I sequenced the last intron of ZFX and ZFY genes in 6 species of primates and 2 species of rodents; I also sequenced the partial genomic sequence of the Ube1x and Ube1y genes of mice and rats. The purposes of my study in addition to estimation of $\alpha$'s in different mammalian species, are to test the hypothesis that most mutations are replication dependent and to examine the generation-time effect on $\alpha$. The $\alpha$ value estimated from the ZFX and ZFY introns of the six primate specise is ${\sim}$6. This estimate is the same as an earlier estimate using only 4 species of primates, but the 95% confidence interval has been reduced from (2, 84) to (2, 33). The estimate of $\alpha$ in the rodents obtained from Zfx and Zfy introns is ${\sim}$1.9, and that deriving from Ube1x and Ube1y introns is ${\sim}$2. Both estimates have a 95% confidence interval from 1 to 3. These two estimates are very close to each other, but are only one-third of that of the primates, suggesting a generation-time effect on $\alpha$. An $\alpha$ of 6 in primates and 2 in rodents are close to the estimates of the male-to-female ratio of the number of germ-cell divisions per generation in humans and mice, which are 6 and 2, respectively, assuming the generation time in humans is 20 years and that in mice is 5 months. These findings suggest that errors during germ-cell DNA replication are the primary source of mutation and that $\alpha$ decreases with decreasing length of generation time. ^

Structure-function analyses of {\it PAX6\/} and the molecular bases of aniridia

PAX6 is a transcription activator that regulates eye development in animals ranging from Drosophila to human. The C-terminal region of PAX6 is proline/serine/threonine-rich (PST) and functions as a potent transactivation domain when attached to a heterologous DNA-binding domain of the yeast transcription factor, GAL4. The PST region comprises 152 amino acids encoded by four exons. The transactivation function of the PST region has not been defined and characterized in detail by in vitro mutagenesis. I dissected the PST domain in two independent systems, a heterologous system using a GAL4 DNA-binding site and the native system of PAX6. In both systems, the results show consistently that all four constituent exons of the PST domain are responsible for the transactivation function. The four exon fragments act cooperatively to stimulate transcription, although none of them can function individually as an independent transactivation domain. Combinations of two or more exon fragments can reconstitute substantial transactivation activity when fused to the DNA-binding domain of GAL4, but they surprisingly do not produce much activity in the context of native PAX6 even though the mutant PAX6 proteins are stable and their DNA-binding function remains unaffected. I conclude that the PAX6 protein contains an unusually large transactivation domain that is evolutionarily conserved to a high degree, and that its full transactivation activity relies on the cooperative action of the four exon fragments.^ Most PAX6 mutations detected in patients with aniridia result in truncations of the protein. Some of the truncation mutations occur in the PST region of PAX6, resulting in mutant proteins that retain their DNA-binding ability but have no significant transactivation activity. It is not clear whether such mutants are true loss-of-function or dominant-negative mutants. I show that these mutants are dominant-negative if they are coexpressed with wild-type PAX6 in cultured cells and that the dominant-negative effects result from enhanced DNA-binding ability of these mutants due to removal of the PST domain. These mutants are able to repress the wild-type PAX6 activity not only at target genes with paired domain binding sites but also at target genes with homeodomain binding sites.^ Mutations in the human PAX6 gene produce various phenotypes, including aniridia, Peters' anomaly, autosomal dominant keratitis, and familial foveal dysplasia. The various phenotypes may arise from different mutations in the same gene. To test this theory, I performed a functional analysis of two missense mutations in the paired domain: the R26G mutation reported in a case of Peters' anomaly, and the I87R mutation identified in a patient with aniridia. While both the R26 and the I87 positions are conserved in the paired boxes of all known PAX genes, X-ray crystallography has shown that only R26 makes contact with DNA. I found that the R26G mutant failed to bind a subset of paired domain binding sites but, surprisingly, bound other sites and successfully transactivated promoters containing those sites. In contrast, the I87R mutant had lost the ability to bind DNA at all tested sites and failed to transactivate promoters. My data support the haploinsufficiency hypothesis of aniridia, and the hypothesis that R26G is a hypomorphic allele. ^

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

Regulation of RNA splicing by nonsense mutations

Translation termination as a result of premature nonsense codon-incorporation in a RNA transcript can lead to the production of aberrant proteins with gain-of-function or dominant negative properties that could have deletrious effects on the cell. T-cell Receptor (TCR) genes acquire premature termination codons two-thirds of the time as a result of the error-prone programmed rearrangement events that normally occur during T-cell development. My studies have focused on the fate of TCR precursor mRNAs in response to in-frame nonsense mutations. ^ Previous published studies from our laboratory have shown that TCR precursor mRNAs are subject to nonsense mediated upregulation of pre-mRNA (NMUP). In this dissertation, I performed substitution and deletion analysis to characterize specific regions of TCR which are required to elicit NMUP. I performed frame- and factor-dependence studies to determine its relationship with other nonsense codon induced responses using several approaches including (i) translation dependence studies (ii) deletion and mutational analysis, as well as (iii) siRNA mediated knockdown of proteins involved. I also addressed the underlying molecular mechanism for this pre-mRNA upregulation by (i) RNA half-life studies using a c-fos inducible promoter, and (ii) a variety of assays to determine pre-mRNA splicing efficiency. ^ Using these approaches, I have identified a region of TCR that is both necessary and sufficient to elicit (NMUP). I have also found that neither cytoplasmic translation machinery nor the protein UPF1 are involved in eliciting this nuclear event. I have shown that the NMUP can be induced not only by nonsense and frameshift mutations, but also missense mutations that disrupt a cis splicing element in the exon that contains the mutation. However, the effect of nonsense mutations on pre-mRNA is unique and distinguishable from that of missense mutations in that nonsense mutations can upregulate pre-mRNA in a frame-dependent manner. Lastly, I provide evidence that NMUP occurs by a mechanism in which nonsense mutations inhibit the splicing of introns. In summary, I have found that TCR precursor mRNAs are subject to multiple forces involving both RNA splicing and translation that can either increase or decrease the levels of these precursor mRNAs. ^

A coalescent analysis for modeling the mutation process in colorectal cancer

Colorectal cancer is the forth most common diagnosed cancer in the United States. Every year about a hundred forty-seven thousand people will be diagnosed with colorectal cancer and fifty-six thousand people lose their lives due to this disease. Most of the hereditary nonpolyposis colorectal cancer (HNPCC) and 12% of the sporadic colorectal cancer show microsatellite instability. Colorectal cancer is a multistep progressive disease. It starts from a mutation in a normal colorectal cell and grows into a clone of cells that further accumulates mutations and finally develops into a malignant tumor. In terms of molecular evolution, the process of colorectal tumor progression represents the acquisition of sequential mutations. ^ Clinical studies use biomarkers such as microsatellite or single nucleotide polymorphisms (SNPs) to study mutation frequencies in colorectal cancer. Microsatellite data obtained from single genome equivalent PCR or small pool PCR can be used to infer tumor progression. Since tumor progression is similar to population evolution, we used an approach known as coalescent, which is well established in population genetics, to analyze this type of data. Coalescent theory has been known to infer the sample's evolutionary path through the analysis of microsatellite data. ^ The simulation results indicate that the constant population size pattern and the rapid tumor growth pattern have different genetic polymorphic patterns. The simulation results were compared with experimental data collected from HNPCC patients. The preliminary result shows the mutation rate in 6 HNPCC patients range from 0.001 to 0.01. The patients' polymorphic patterns are similar to the constant population size pattern which implies the tumor progression is through multilineage persistence instead of clonal sequential evolution. The results should be further verified using a larger dataset. ^

NRAS and BRAF mutations in primary cutaneous melanoma: A comparison of mutation rates between radial and vertical growth phases in individual tumors

Primary cutaneous melanoma is a cancer arising from melanocytes in the skin. In recent decades the incidence of this malignancy has increased significantly. Mortality rates are high for patients with tumors measuring over a few millimeters in thickness. Response rates to conventional radiation and chemotherapy are very low in patients with metastatic melanoma. New therapies targeting melanoma’s aberrant cell signaling pathways such as the MAP Kinase pathway are being developed. Mutations of NRAS and BRAF genes are quite common in cutaneous melanoma and lead to constitutive activation of the MAP Kinase pathway. This study tests the hypothesis that NRAS and BRAF mutations increase as a tumor progresses from the noninvasive radial growth phase (RGP) to the invasive vertical growth phase (VGP). Laser capture microdissection was used to obtain separate, pure tumor DNA samples from the RGP and VGP of thirty primary cutaneous melanomas. PCR was used to amplify NRAS exon 2 and BRAF exon 15 tumor DNA. The amplified DNA was sequenced and analyzed for mutations. An overall mutation rate of 74% was obtained for the twenty-three melanomas in which there were complete sequence results. With the exception of one melanoma NRAS and BRAF mutations were mutually exclusive. All seven NRAS exon 2 mutations involved codon 61. Three of these melanomas had mutations in both the RGP and VGP. The remaining four tumors were wild type for NRAS exon 2 in the RGP but mutated in the VGP. Of the fifteen BRAF exon 15 mutated melanomas all but one involved codon 600. Twelve of the fifteen BRAF exon 15 mutations were the T1799A type. Nine of the fifteen BRAF mutated tumors had the same mutation in both the RGP and VGP. Five of fifteen melanomas had wild type RGP DNA and BRAF exon 15 mutated VGP DNA. A single melanoma had BRAF exon 15 mutated DNA in the RGP and wild type DNA in the VGP. Overall, these results suggest a trend toward the acquisition of NRAS and BRAF mutations as cutaneous melanomas change from a noninvasive to an invasive, potentially deadly cancer.^

Therapeutic potential of p53 restoration in spontaneous tumors

Over 80% of p53 mutations found in human cancers are p53 missense mutations. Recent studies have shown that p53 restoration leads to tumor regression in mice with p53 deletions, but the therapeutic efficacy of p53 restoration in tumors containing p53 missense mutations has not been evaluated. Since p53 mutant such as p53R172H has gain-of-function activities and dominant-negative effect that repress wild type p53, the activity of restored wild-type p53 might be compromised by the mutant p53 in tumors. We hypothesized that p53 restoration in tumors with the p53R172H mutation may be less therapeutically effective as p53 restoration in tumors null for p53. I tested this hypothesis by comparison of the therapeutic outcomes of p53 restoration in mice with spontaneous tumors that either lacked p53 or contained the p53R172H mutation. While p53 restoration causes tumor regression in mice lacking p53, the same p53 restoration halts tumor progression in mice with the p53R172H mutation. This phenotypic difference suggests a dominant-negative activity of the mutant p53. Moreover, I showed that the mutant p53 only inhibits part of the activity of the restored wild-type p53 and that the remaining wild-type activity still causes a delay in tumor progression. We conclude that p53 restoration has therapeutic potential in p53R172H tumors via suppression of tumor progression. This knowledge is of critical importance for p53 targeted cancer therapy because many patients with cancers harbor p53 missense mutations rather p53-null mutations. Since p53R172H mutation represents one of the most frequent and potent p53 missense mutations observed in human cancers, the current findings implicates that p53 restoration may be therapeutically important not only in human cancers characterized by loss of p53 alleles but also in those in which p53 missense mutations play an important pathogenetic role. ^

Functional analysis of p53 phosphorylation and a p53 hot spot mutation

The tumor suppressor p53 is a phosphoprotein which functions as a transcriptional activator. By monitoring the transcriptional activity, we studied how p53 functions is regulated in relation to cell growth and contact inhibition. When cells were arrested at G1 phase of the cell cycle by contact inhibition, we found that p53 transactivation function was suppressed. When contact inhibition was overridden by cyclin E overexpression which stimulates cell cycle progression, p53 function was restored. This observation led to the development of a cell density assay to study the regulation of p53 function during cell cycle for the functional significance of p53 phosphorylation. The murine p53 is phosphorylated at serines 7, 9, 12, 18, 37, 312 and 389. To understand the role of p53 phosphorylation, we generated p53 constructs encoding serine-to-alanine or serine-to-glutamate mutations at these codons. The transcriptional activity were measured in cells capable of contact inhibition. In low-density cycling cells, no difference in transcriptional activity was found between wild type p53 and any of the mutants. In contact-inhibited cells, however, only mutations of p53 at serine 389 resulted in altered responses to cell cycle arrest and to cyclin E overexpression. The mutant with serine-to-glutamate substitution at codon 389 retained its function in contact inhibited cells. Cyclin E overexpression in these cells induced p53 phosphorylation at serine 389. Furthermore, we showed that phosphorylation at serine 389 regulates p53 DNA binding activity. Our findings implicate that phosphorylation is an important mechanism for p53 activation.^ p53 is the most frequently mutated gene in human tumors. To study the mechanism of p53 inactivation by mutations, we carried out detailed analysis of a murine p53 mutation with an arginine-to-tryptophane substitution at codon 245. The corresponding human p53 mutation at amino acid 248 is the most frequently mutated codon in tumors. We showed that this mutant is inactive in suppressing focus formation, binding to DNA and transactivation. Structural analysis revealed that this mutant assumes the wild type protein conformation. These findings define a novel class of p53 mutations and help to understand structure-function relationship of p53. ^

Altered tumorigenic phenotype in mice with a hypomorphic p53 mutant allele

Mutations in the p53 tumor suppressor gene are found in over 50% of human tumors and in the germline of Li-Fraumeni syndrome families. About 80% of these mutations are missense in nature. In order to study how p53 missense mutations affect tumorigenesis in vivo, we focused on the murine p53 arg-to-his mutation at amino acid 172, which corresponds to the human hot spot mutation at amino acid 175. The double replacement procedure was employed to introduce the p53 R172H mutation into the p53 locus of ES cells and mice were generated. An additional 1bp deletion in the intron 2 splice acceptor site was detected in the same allele in mice. We named this allele p53R172HΔg. This allele makes a small amount of full length p53 mutant protein. ^ Spontaneous tumor formation and survival were studied in these mice. Mice heterozygous for the p53R172HΔg allele showed 50% survival at 17 months of age, similar to the p53+/− mice. Moreover, the p53R172HΔg/+ mice showed a distinct tumor spectrum: 55% sarcomas, including osteosarcoms, fibrosarcomas and angiosarcomas; 27% carcinomas, including lung adenocarcinomas, squamous cell carcinomas, hepatocellular carcinomas and islet cell carcinomas; and 18% lymphomas. Compared to the p53+/− mice, there was a clear increase in the frequency of carcinoma development and a decrease in lymphoma incidence. Among the sarcomas that developed, fibrosarcomas in the skin were also more frequently observed. More importantly, osteosarcomas and carinomas that developed in the p53R172HΔg/+ mice metastasized at very high frequency (64% and 67%, respectively) compared with less than 10% in the p53+/− mice. The metastatic lesions were usually found in lung and liver, and less frequently in other tissues. The altered tumor spectrum in the mice and increased metastatic potential of the tumors suggested that the p53R172H mutation represents a gain-of-function. ^ Mouse embryonic fibroblasts (MEFs) from the mice homozygous and heterozygous for the p53R172HΔg allele were studied for growth characteristics, immortalization potential and genomic instability. All of the p53R172HΔg /+ MEF lines are immortalized under a 3T3 protocol while under the same protocol p53+/− MEFs are not immortalized. Karyotype analysis showed a persistent appearance of chromosome end-to-end fusion in the MEFs both homozygous and heterozygous for the p53R172HΔg allele. These observations suggest that increased genomic instability in the cells may cause the altered tumor phenotypes. ^