Estimation of the incidence of a rare genetic disease through a two-tier mutation survey.

Recent attempts to detect mutations involving single base changes or small deletions that are specific to genetic diseases provide an opportunity to develop a two-tier mutation-screening program through which incidence of rare genetic disorders and gene carriers may be precisely estimated. A two-tier survey consists of mutation screening in a sample of patients with specific genetic disorders and in a second sample of newborns from the same population in which mutation frequency is evaluated. We provide the statistical basis for evaluating the incidence of affected and gene carriers in such two-tier mutation-screening surveys, from which the precision of the estimates is derived. Sample-size requirements of such two-tier mutation-screening surveys are evaluated. Considering examples of cystic fibrosis (CF) and medium-chain acyl-CoA dehydrogenase deficiency (MCAD), the two most frequent autosomal recessive disease in Caucasian populations and the two most frequent mutations (delta F508 and G985) that occur on these disease allele-bearing chromosomes, we show that, with 50-100 patients and a 20-fold larger sample of newborns screened for these mutations, the incidence of such diseases and their gene carriers in a population may be quite reliably estimated. The theory developed here is also applicable to rare autosomal dominant diseases for which disease-specific mutations are found.

Mutation and association analysis of the PVR and PVRL2 genes in patients with non-syndromic cleft lip and palate.

Orofacial clefts (OFC; MIM 119530) are among the most common major birth defects. Here, we carried out mutation screening of the PVR and PVRL2 genes, which are both located at an OFC linkage region at 19q13 (OFC3) and are closely related to PVRL1, which has been associated with both syndromic and non-syndromic cleft lip and palate (nsCLP). We screened a total of 73 nsCLP patients and 105 non-cleft controls from the USA for variants in PVR and PVRL2, including all exons and encompassing all isoforms. We identified four variants in PVR and five in PVRL2. One non-synonymous PVR variant, A67T, was more frequent among nsCLP patients than among normal controls, but this difference did not achieve statistical significance.

Mutation analysis of the PVRL1 gene in caucasians with nonsyndromic cleft lip/palate.

Nonsyndromic cleft lip with or without cleft palate (nsCL/P, MIM 119530) is perhaps the most common major birth defect. Homozygous PVRL1 loss-of-function mutations result in an autosomal recessive CL/P syndrome, CLPED1, and a PVRL1 nonsense mutation is associated with sporadic nsCL/P in Northern Venezuela. To address the more general role of PVRL1 variation in risk of nsCL/P, we carried out mutation analysis of PVRL1 in North American and Australian nsCL/P cases and population-matched controls. We identified a total of 15 variants, 5 of which were seen in both populations and 1 of which, an in-frame insertion at Glu442, was more frequent in patients than in controls in both populations, though the difference was not statistically significant. Another variant, which is specific to the PVRL1 beta (HIgR) isoform, S447L, was marginally associated with nsCL/P in North American Caucasian patients, but not in Australian patients, and overall variants that affect the beta-isoform were significantly more frequent among North American patients. One Australian patient had a splice junction mutation of PVRL1. Our results suggest that PVRL1 may play a minor role in susceptibility to the occurrence of nsCL/P in some Caucasian populations, and that variation involving the beta (HIgR) isoform might have particular importance for risk of orofacial clefts. Nevertheless, these results underscore the need for studies that involve very large numbers when assessing the possible role of rare variants in risk of complex traits such as nsCL/P.

Severe aortic and arterial aneurysms associated with a TGFBR2 mutation.

BACKGROUND: A 24-year-old man presented with previously diagnosed Marfan's syndrome. Since the age of 9 years, he had undergone eight cardiovascular procedures to treat rapidly progressive aneurysms, dissection and tortuous vascular disease involving the aortic root and arch, the thoracoabdominal aorta, and brachiocephalic, vertebral, internal thoracic and superior mesenteric arteries. Throughout this extensive series of cardiovascular surgical repairs, he recovered without stroke, paraplegia or renal impairment. INVESTIGATIONS: CT scans, arteriogram, genetic mutation screening of transforming growth factor beta receptors 1 and 2. DIAGNOSIS: Diffuse and rapidly progressing vascular disease in a patient who met the diagnostic criteria for Marfan's syndrome, but was later rediagnosed with Loeys-Dietz syndrome. Genetic testing also revealed a de novo mutation in transforming growth factor beta receptor 2. MANAGEMENT: Regular cardiovascular surveillance for aneurysms and dissections, and aggressive surgical treatment of vascular disease.

DYSREGULATION OF MEOX2 FOLLOWING WT1 MUTATION IN KIDNEY DEVELOPMENT AND WILMS TUMORIGENESIS

Wilms tumor (WT) is a childhood tumor of the kidney and a productive model for understanding the role of genetic alteration and interactions in tumorigenesis. The Wilms tumor gene 1 (WT1) is a transcriptional factor and one of the few genes known to have genetic alterations in WT and has been shown be inactivated in 20% of WTs. However, the mechanisms of how WT1 mutations lead to Wilms tumorigenesis and its influence on downstream genes are unknown. Since it has been established that WT1 is a transcriptional regulator, it has been hypothesized that the loss of WT1 leads to the dysregulation of downstream genes, in turn result in the formation of WTs. To identify the dysregulated downstream genes following WT1 mutations, an Affymetrix GeneChip Human Genome Array was previously conducted to assess the differentially expressed genes in the WT1-wildtype human and WT1-mutant human WTs. Approximately 700 genes were identified as being significantly dysregulated. These genes were further prioritized based on their statistical significance, fold change, chromosomal region, spatial pattern of gene expression and known or putative cellular functions. Mesenchyme homeobox 2 (MEOX2) was one of the most significantly upregulated genes in WT1-mutant WT. MEOX2 is known to play a role in cell proliferation, apoptosis, and differentiation. In addition to its biological roles, it is expressed during early kidney development in the condensed mesenchyme similar to WT1. Furthermore, the use of the Match® web-based tool from the BIOBASE Biological Data base identified a significant predicted WT1 binding site within the first intron of MEOX2. The similarity in spatial gene expression in the developing kidney and the significant predicted WT1 binding site found in the first intron of MEOX2 lead to the development of my hypothesis that MEOX2 is upregulated via a WT1-dependent manner. Here as a part of my master’s work, I have validated the Affymetrix GeneChip Human Genome Array data using an independent set of Wilms tumors. MEOX2 remained upregulated in the mutant WT1 Wilms tumor by 41-fold. Wt1 and Meox2 gene expression were assessed in murine newborn kidney; both Wt1 and Meox2 were expressed in the condensed, undifferentiated metanephric mesenchyme. I have shown that the in vivo ablation of Wt1 during embryonic development at embryonic day (E) 13.5 resulted in the slight increase of Meox2 gene expression by two fold. In order to functionally demonstrate the effect of the loss of Wt1 on Meox2 gene expression in undifferentiated metanephric mesenchyme, I have generated a kidney mesenchymal cell line to genetically ablate Wt1 in vitro by adenoviral infection. The ablation of Wt1 in the kidney mesenchymal cell line resulted in the upregulation of Meox2 by 61-fold. Moreover, the upregulation of Meox2 resulted in the significant induction of p21 and Itgb5. In addition to the dysregulation of these genes the ablation of Wt1 in the kidney mesenchymal cells resulted in decrease in cell growth and loss of cellular adherence. However, it is uncertain whether the upregulation of Meox2 caused this particular cellular phenotype. Overall, I have demonstrated that the upregulation of Meox2 is Wt1-dependent during early kidney development.

Detection and sequence analysis of a nickel-induced mutation in a retroviral model system

We have developed a novel way to assess the mutagenicity of environmentally important metal carcinogens, such as nickel, by creating a positive selection system based upon the conditional expression of a retroviral transforming gene. The target gene is the v-mos gene in MuSVts110, a murine retrovirus possessing a growth temperature dependent defect in expression of the transforming gene due to viral RNA splicing. In normal rat kidney cells infected with MuSVts110 (6m2 cells), splicing of the MuSVts110 RNA to form the mRNA from which the transforming protein, p85$\sp{\rm gag-mos}$, is translated is growth-temperature dependent, occurring at 33 C and below but not at 39 C and above. This splicing "defect" is mediated by cis-acting viral sequences. Nickel chloride treatment of 6m2 cells followed by growth at 39 C, allowed the selection of "revertant" cells which constitutively express p85$\sp{\rm gag-mos}$ due to stable changes in the viral RNA splicing phenotype, suggesting that nickel, a carcinogen whose mutagenicity has not been well established, could induce mutations in mammalian genes. We also show by direct sequencing of PCR-amplified integrated MuSVts110 DNA from a 6m2 nickel-revertant cell line that the nickel-induced mutation affecting the splicing phenotype is a cis-acting 70-base duplication of a region of the viral DNA surrounding the 3$\sp\prime$ splice site. These findings provide the first example of the molecular basis for a nickel-induced DNA lesion and establish the mutagenicity of this potent carcinogen. ^

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

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

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

EFFECTS OF THE ACTA2 R258C MUTATION ON VASCULAR SMOOTH MUSCLE CELL PHENOTYPE AND PROPERTIES

Thoracic Aortic Aneurysms and Dissections (TAAD) are the fifteenth leading cause of death in the United States. About 15% of TAAD patients have family history of the disease. The most commonly mutated gene in these families is ACTA2, encoding smooth muscle-specific α-actin. ACTA2 missense mutations predispose individuals both to TAAD and to vascular occlusive disease of small, muscular arteries. Mice carrying an Acta2 R258C mutant transgene with a wildtype Acta2 promoter were generated and bred with Acta2-/- mice to decrease the wildtype: mutant Acta2 ratio. Acta2+/+ R258C TGmice have decreased aortic contractility without aortic disease. Acta2+/- R258C TG mice, however, have significant aortic dilatations by 12 weeks of age and a hyperproliferative response to injury. We characterized smooth muscle cells (SMCs) from bothmouse models under the hypothesis that mutant α-actin has a dominant negative effect, leading to impaired contractile filament formation/stability, improper focal adhesion maturation and increased proliferation. Explanted aortic SMCs from Acta2+/+ R258C TG mice are differentiated - they form intact filaments, express higher levels of contractile markers compared to wildtype SMCs and have predominantly nuclear Myocardin-Related Transcription Factor A (MRTF-A) localization. However, ultracentrifugation assays showed large unpolymerized actin fractions, suggesting that the filaments are brittle. In contrast, Acta2+/- R258C TG SMCs are less well-differentiated, with pools of unpolymerized actin, more cytoplasmic MRTF-A and decreased contractile protein expression compared to wildtype cells. Ultracentrifugation assays after treating Acta2+/- R258C TGSMCs with phalloidin showed actin filament fractions, indicating that mutant α-actin can polymerize into filaments. Both Acta2+/+ R258C TGand Acta2+/- R258C TGSMCs have larger and more peripheral focal adhesions compared to wildtype SMCs. Rac1 was more activated in Acta2+/+ R258C TGSMCs; both Rac1 and RhoA were less activated in Acta2+/- R258C TG SMCs, and FAK was more activated in both transgenic SMC lines compared to wildtype. Proliferation in both cell lines was significantly increased compared to wildtype cells and could be partially attenuated by inhibition of FAK or PDGFRβ. These data support a dominant negative effect of the Acta2 R258C mutation on the SMC phenotype, with increasing phenotypic severity when wildtype: mutant α-actin levels are decreased.

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

Gain-of-function of the p53R172H mutation in a mouse model of Li -Fraumeni syndrome

Missense mutations in the p53 tumor-suppressor gene are the most common alterations of p53 in somatic tumors and in patients with Li-Fraumeni syndrome. p53 missense mutations occur in the DNA binding region and disrupt the ability of p53 to activate transcription. In vitro studies have shown that some p53 missense mutants have a gain-of-function or dominant-negative activity. ^ The p53 175 Arg-to-His (p53 R175H) mutation in humans has been shown to have dominant-negative and gain-of-function properties in vitro. This mutation is observed in the germline of individuals with Li-Fraumeni syndrome. To accurately model Li-Fraumeni syndrome and to examine the mechanistic nature of a gain-of-function missense mutation on in vivo tumorigenesis, we generated and characterized a mouse with the corresponding mutation, p53 R172H. p53R172H homozygous and heterozygous mice developed similar tumor spectra and survival curves as p53 −/− and p53+/− mice, respectively. However, tumors in p53+/R172H mice metastasized to various organs with high frequency, suggesting a gain-of-function phenotype by p53R172H in vivo. Mouse embryonic fibroblasts (MEFs) from p53R172H mice also showed gain-of-function phenotypes in cell proliferation, DNA synthesis, and transformation potential, while cells from p53+/− and p53−/− mice did not. ^ To mechanistically characterize the gain-of-function phenotype of the p53R172H mutant, the role of p53 family members, p63 and p73, was analyzed. Disruption of p63 and p73 by siRNAs in p53 −/− MEFs increased transformation potential and reinitiated DNA synthesis to levels observed in p53R172H/R172H cells. Additionally, p63 and p73 were bound and functionally inactivated by p53R172H in metastatic p53 R172H tumor-derived cell lines, indicating a role for the p53 family members in the gain-of-function phenotype. This study provides in vivo evidence for the gain-of-function effect of p53 missense mutations and more accurately models the Li-Fraumeni syndrome. ^