THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

Decomposition of $\rm R\times C$ contingency table chi -square: Application to binned DNA fragment size data and population structure analysis

The purpose of this research is to develop a new statistical method to determine the minimum set of rows (R) in a R x C contingency table of discrete data that explains the dependence of observations. The statistical power of the method will be empirically determined by computer simulation to judge its efficiency over the presently existing methods. The method will be applied to data on DNA fragment length variation at six VNTR loci in over 72 populations from five major racial groups of human (total sample size is over 15,000 individuals; each sample having at least 50 individuals). DNA fragment lengths grouped in bins will form the basis of studying inter-population DNA variation within the racial groups are significant, will provide a rigorous re-binning procedure for forensic computation of DNA profile frequencies that takes into account intra-racial DNA variation among populations. ^

SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.

A bacterial cytotoxin identifies the RhoA exchange factor Net1 as a key effector in the response to DNA damage.

BACKGROUND: Exposure of adherent cells to DNA damaging agents, such as the bacterial cytolethal distending toxin (CDT) or ionizing radiations (IR), activates the small GTPase RhoA, which promotes the formation of actin stress fibers and delays cell death. The signalling intermediates that regulate RhoA activation and promote cell survival are unknown. PRINCIPAL FINDINGS: We demonstrate that the nuclear RhoA-specific Guanine nucleotide Exchange Factor (GEF) Net1 becomes dephosphorylated at a critical inhibitory site in cells exposed to CDT or IR. Expression of a dominant negative Net1 or Net1 knock down by iRNA prevented RhoA activation, inhibited the formation of stress fibers, and enhanced cell death, indicating that Net1 activation is required for this RhoA-mediated responses to genotoxic stress. The Net1 and RhoA-dependent signals involved activation of the Mitogen-Activated Protein Kinase p38 and its downstream target MAPK-activated protein kinase 2. SIGNIFICANCE: Our data highlight the importance of Net1 in controlling RhoA and p38 MAPK mediated cell survival in cells exposed to DNA damaging agents and illustrate a molecular pathway whereby chronic exposure to a bacterial toxin may promote genomic instability.

E2F1 AND TUMOR SUPPRESSION: THE ROLE OF p21, miRNAS, AND THE DNA DAMAGE RESPONSE

E2F1 is a multi-faceted protein that has roles in a number of important cellular processes including cell cycle regulation, apoptosis, proliferation, and the DNA damage response (DDR). Moreover, E2F1 has opposing roles in tumor development, acting as either a tumor suppressor or an oncogene depending on the context. In human cancer, E2F1 is often deregulated through aberrations in the Rb-p16INK4a-cyclin D1 pathway. In these studies we examined three mechanisms by which E2F1 might mediate its tumor suppressive properties: p21-induced senescence, miRNAs, and the DNA damage response. We found that E2F1 acts as a tumor suppressor in response to ras activation through a non-apoptotic mechanism requiring ARF and p53, but not p21. However, p21-loss inhibited two-stage chemical carcinogenesis in FVB mice. In response to E2F1 overexpression, we found that 22 miRNAs are differentially regulated in mouse epidermis, including let-7a, let-7c, and miR-301. Additionally, regulation of miR-301 involves binding of E2F1 to its promoter. Finally, our data indicate a role for E2F1 at sites of DNA damage requiring E2F1’s phosphorylation at serine 31 which may involve DNA repair. Further, this role in the DDR may affect tumor aggressiveness and multiplicity. In all, we have explored three mechanisms for E2F1-induced tumor suppression and identified E2F1’s role in the DNA damage response as a likely contributor to this phenomenon.

POLG DNA testing as an emerging standard of care before instituting valproic acid therapy for pediatric seizure disorders.

PURPOSE: To review our clinical experience and determine if there are appropriate signs and symptoms to consider POLG sequencing prior to valproic acid (VPA) dosing in patients with seizures. METHODS: Four patients who developed VPA-induced hepatotoxicity were examined for POLG sequence variations. A subsequent chart review was used to describe clinical course prior to and after VPA dosing. RESULTS: Four patients of multiple different ethnicities, age 3-18 years, developed VPA-induced hepatotoxicity. All were given VPA due to intractable partial seizures. Three of the patients had developed epilepsia partialis continua. The time from VPA exposure to liver failure was between 2 and 3 months. Liver failure was reversible in one patient. Molecular studies revealed homozygous p.R597W or p.A467T mutations in two patients. The other two patients showed compound heterozygous mutations, p.A467T/p.Q68X and p.L83P/p.G888S. Clinical findings and POLG mutations were diagnostic of Alpers-Huttenlocher syndrome. CONCLUSION: Our cases underscore several important findings: POLG mutations have been observed in every ethnic group studied to date; early predominance of epileptiform discharges over the occipital region is common in POLG-induced epilepsy; the EEG and MRI findings varying between patients and stages of the disease; and VPA dosing at any stage of Alpers-Huttenlocher syndrome can precipitate liver failure. Our data support an emerging proposal that POLG gene testing should be considered in any child or adolescent who presents or develops intractable seizures with or without status epilepticus or epilepsia partialis continua, particularly when there is a history of psychomotor regression.

{\it In vivo\/} induction of DNA changes in cervicovaginal epithelium by perinatal estrogen exposure

Epidemiological studies have associated estrogens with human neoplasm such as the endometrium, cervix, vagina, breast, and liver. Perinatal exposure to natural (17$\beta$-estradiol (17$\beta$-E$\sb2)\rbrack$ and synthetic (diethylstilbestrol (DES)) estrogens induces neoplastic changes in humans and rodents. Previous studies demonstrated that neonatal 17$\beta$-E$\sb2$ treatment increased the nuclear DNA content of mouse cervicovaginal epithelium that preceded histologically evident neoplasia. In order to determine whether this effect was specific to 17$\beta$-E$\sb2,$ associated with chromosomal changes, and relevant to the human, female BALB/c mice were treated neonatally with either 17$\alpha$-estradiol (17$\alpha$-E$\sb2)$ and 5$\beta$-dihydrotestosterone ($5\beta$-DHT), both inactive steroids in adult reproductive tissue, or 17$\beta$-E$\sb2.$ Ten-day-old mice received pellet implants of 17$\beta$-E$\sb2,$ 17$\alpha$-E$\sb2,$ $5\beta$-DHT, or cholesterol. Seventy-day-old cervicovaginal tracts were examined histologically and flow cytometrically. 17$\beta$-E$\sb2$-treated animals were evaluated by fluorescent in situ hybridization (FISH) using a probe specific for chromosome 1. Trisomy of chromosomes 1, 7, 11, and 17 was evaluated by FISH in cervicovaginal material from 19 DES-exposed and 19 control patients.^ $17\beta$-E$\sb2, 17\alpha$-E$\sb2$, and $5\beta$-DHT-induced dramatic developmental and histological changes in the cervicovaginal tract, including hypospadia, hyperplasia, and persistent cornification. The changes induced by 17$\alpha$-E$\sb2$ were equivalent to 17$\beta$-E$\sb2.$ Neonatal 17$\alpha$-E$\sb2$-induced adenosquamous cervicovaginal tumors at 24 months. 17$\alpha$-E$\sb2$ and $5\beta$-DHT significantly increased the nuclear DNA content over control animals, but at significantly lower levels than 17$\beta$-E$\sb2.$ DNA ploidy changes were highest (80%) in animals treated neonatally and secondarily with 17$\beta$-E$\sb2.$ Secondary 17$\alpha$-E$\sb2$ and $5\beta$-DHT administration, unlike 17$\beta$-E$\sb2,$ didn't significantly increase DNA content. Chromosome 1 trisomy incidence was 66% in neonatal 17$\beta$-E$\sb2$-treated animals. Trisomy was evident in 4 DES-exposed patients: one patient with trisomy of chromosomes 1, 7, and 11; one patient with chromosome 7 trisomy; and two patients with chromosome 1 trisomy. These data demonstrated the biological effects of 17$\alpha$-E$\sb2$ and $5\beta$-DHT were age-dependent, 17$\alpha$-E$\sb2$ was equivalent to 17$\beta$-E$\sb2$ and tumorigenic when administered neonatally, and histological changes were not steroid specific. Chromosomal changes were associated with increased nuclear DNA content and chromosomal changes may be an early event in the development of tumors in human DES-exposed tissues. ^

Functional studies on DNA double-strand break repair proteins (MmRad51, Ku80) in mice

RecA in Escherichia coli and it's homologue, ScRad51 in Saccharomyces cerevisiae, play important roles in recombinational repair. ScRad51 homologues have been discovered in a wide range of organisms including Schizosaccharomyces pombe, lily, chicken, mouse and human. To date there is no direct evidence to describe that mouse Rad51(MmRad51) is involved in DNA double-strand break repair. In order to elucidate the role of MmRad51 in vivo, it was mutated by the embryonic stem (ES) cell/gene targeting technology in mice. The mutant embryos arrested in development shortly after implantation. There was a decrease in cell proliferation followed by programmed cell death, and trophectoderm-derived cells were sensitive to $\gamma$-radiation. Severe chromosome loss was observed in most mitotically dividing cells. The mutant embryos lived longer and developed further in a p53 mutant background; however, double-mutant embryonic fibroblasts failed to proliferate in tissue culture, reflecting the embryos limited life span. Based on these data, MmRad51 repairs DNA damage induced by $\gamma$-radiation, is needed to maintain euplody, and plays an important role in proliferating cells.^ Ku is a heterodimer of 70 and 80 kDs subunit, which binds to DNA ends and other altered DNA structures such as hairpins, nicks, and gaps. In addition, Ku is required for DNA-PK activity through a direct association. Although the biochemical properties of Ku and DNA-PKcs have been characterized in cells, their physiological functions are not clear. In order to understand the function of Ku in vivo, we generated mice homozygous for a mutation of the Ku80 gene. Ku80-deficient mice, like scid mice, showed severe immunodeficiency due to a impairment of V(D)J recombination. Mutant mice were semiviable and runted, cells derived from mutant embryos displayed hypersensitivity to $\gamma$-radiation, a decreased growth rate, a slow entry into S phase, altered colony size distributions, and a short life span. Based on these results, mutant cells and mice appeared to prematurely age. ^

Statistical characterization and development of summary measures of non-normal distributions of nuclear morphometry data from prostate cancer cells for use as prognostic indicators

Nuclear morphometry (NM) uses image analysis to measure features of the cell nucleus which are classified as: bulk properties, shape or form, and DNA distribution. Studies have used these measurements as diagnostic and prognostic indicators of disease with inconclusive results. The distributional properties of these variables have not been systematically investigated although much of the medical data exhibit nonnormal distributions. Measurements are done on several hundred cells per patient so summary measurements reflecting the underlying distribution are needed.^ Distributional characteristics of 34 NM variables from prostate cancer cells were investigated using graphical and analytical techniques. Cells per sample ranged from 52 to 458. A small sample of patients with benign prostatic hyperplasia (BPH), representing non-cancer cells, was used for general comparison with the cancer cells.^ Data transformations such as log, square root and 1/x did not yield normality as measured by the Shapiro-Wilks test for normality. A modulus transformation, used for distributions having abnormal kurtosis values, also did not produce normality.^ Kernel density histograms of the 34 variables exhibited non-normality and 18 variables also exhibited bimodality. A bimodality coefficient was calculated and 3 variables: DNA concentration, shape and elongation, showed the strongest evidence of bimodality and were studied further.^ Two analytical approaches were used to obtain a summary measure for each variable for each patient: cluster analysis to determine significant clusters and a mixture model analysis using a two component model having a Gaussian distribution with equal variances. The mixture component parameters were used to bootstrap the log likelihood ratio to determine the significant number of components, 1 or 2. These summary measures were used as predictors of disease severity in several proportional odds logistic regression models. The disease severity scale had 5 levels and was constructed of 3 components: extracapsulary penetration (ECP), lymph node involvement (LN+) and seminal vesicle involvement (SV+) which represent surrogate measures of prognosis. The summary measures were not strong predictors of disease severity. There was some indication from the mixture model results that there were changes in mean levels and proportions of the components in the lower severity levels. ^

Delineation of a novel genetic region at 5q11 harboring tumor suppressor gene(s) and identification of the telomeric DNA as a marker for human prostate cancer

Prostate cancer is the most commonly diagnosed cancer and the second leading cause of cancer mortality in American men. The distinction between those cases of prostate cancer destined to progress rapidly to lethal metastatic disease and those with little likelihood of causing morbidity and mortality is a major goal of current research. Some type of diagnostic method is urgently needed to identify which histological prostate cancers have completed the progression to a stage that will produce a life-threatening disease, thus requiring immediate therapeutic intervention. The objectives of this dissertation are to delineate a novel genetic region harboring tumor suppressor gene(s) and to identify a marker for prostate tumorigenesis. I first established an in vitro cell model system from a human prostate epithelial cells derived from tissue fragments surrounding a prostate tumor in a patient with prostatic adenocarcinoma. Since chromosome 5 abnormality was present in early, middle and late passages of this cell model system, I examined long-term established prostate cancer cell lines for this chromosome abnormality. The results implicated the region surrounding marker D5S2068 as the locus of interest for further experimentation and location of a tumor suppressor gene in human prostate cancer. ^ Cancer is a group of complex genetic diseases with uncontrolled cell; division and prostate cancer is no exception. I determined if telomeric DNA, and telomerase activity, alone or together, could serve as biomarkers of prostate tumorigenesis. I studied three newly established human prostate cancer cell lines and three fibroblast cell cultures derived from prostate tissues. In conclusion, my data reveal that in the presence of telomerase activity, telomeric repeats are maintained at a certain optimal length, and analysis of telomeric DNA variations might serve as early diagnostic and prognostic biomarkers for prostate cancer. (Abstract shortened by UMI.)^

CHROMOSOMAL MAPPING, LINKAGE RELATIONSHIPS, AND EVOLUTIONARY STUDIES OF THE HUMAN ANONYMOUS DNA CLONE D1S1 (IN SITU HYBRIDIZATION, PRIMATES, GENE MAPPING, AUTOSOMAL DOMINANT RETINITIS PIGMENTOSA, GENE DUPLICATION)

D1S1, an anonymous human DNA clone originally called (lamda)Ch4-H3 or (lamda)H3, was the first single copy mapped to a human chromosome (1p36) by in situ hybridization. The chromosomal assignment has been confirmed in other laboratories by repeating the in situ hybridization but not by another method. In the present study, hybridization to a panel of hamster-human somatic cell hybrids revealed copies of D1S1 on both chromosomes 1 and 3. Subcloning D1S1 showed that the D1S1 clone itself is from chromosome 3, and the sequence detected by in situ hybridization is at least two copies of part of the chromosome 3 copy. This finding demonstrates the importance of verifying gene mapping with two methods and questions the accuracy of in situ hybridization mapping.^ Non-human mammals have only one copy of D1S1, and the non-human primate D1S1 map closely resembles the human chromosome 3 copy. Thus, the human chromosome 1 copies appear to be part of a very recent duplication that occurred after the divergence between humans and the other great apes.^ A moderately informative HindIII D1S1 RFLP was mapped to chromosome 3. This marker and 12 protein markers were applied to a linkage study of autosomal dominant retinitis pigmentosa (ADRP). None of the markers proved linkage, but adding the three families examined to previously published data raises the ADRP:Rh lod score to 1.92 at (THETA) = 0.30. ^

CHROMOSOMAL ORIGIN OF DM-DNA

Double minutes (dm) are small chromatin particles of 0.3 microns diameter found only in the metaphase cells of human and murine tumors. Dm are unique cytogenetic structures since their numbers per cell show wide variation. At cell division, dm are retained despite the lack of centromeres. In squash preparations, dm show clustering often in association with chromosomes. Human carcinoma cell line SW613-S18 was found to have large numbers of dm and biological characteristics favorable for mitotic synchronization and chromosome isolation experiments.^ S18 cells were synchronized to mitosis with metabolic and mitotic blocking compounds. Mitotic cells were lysed to release chromosomes and dm from the mitotic spindle and the resulting suspensions were fractionated to enrich for dm. The DNA in enriched fractions was characterized. The reassociation kinetics of dm-DNA driven with placental human DNA was similar to the reassociation curve of labeled placental DNA under similar conditions. In situ hybridization of dm-DNA to tumor and normal metaphase cells showed grain localization over the entire karyotype. Dm-DNA was shown by pulse chase DNA replication experiments to replicate during early and mid S-phase of the cell cycle, but not in late S-phase. In addition, BrdUrd incorporation studies showed that dm-DNA replicates only once during the S-phase. Premature chromosome condensation studies suggest the basis of numerical heterogeneity of dm is nondisjunction, not anomalous or unscheduled DNA replication.^ These data and previous cytochemical banding studies of dm in SW613-S18 indicate that dm-DNA is chromosomal in origin. No evidence of gene amplification was found in the DNA reassociation data. It is likely that dm-DNA represents the pale-staining G-band regions of the human karyotype in this cell line. ^

CYTOGENETIC AND DNA CONTENT ANALYSIS IN MULTIPLE MYELOMA

In this dissertation, the cytogenetic characteristics of bone marrow cells from 41 multiple myeloma patients were investigated. These cytogenetic data were correlated with the total DNA content as measured by flow cytometry. Both the cytogenetic information and DNA content were then correlated with clinical data to determine if diagnosis and prognosis of multiple myeloma could be improved.^ One hundred percent of the patients demonstrated abnormal chromosome numbers per metaphase. The average chromosome number per metaphase ranged from 42 to 49.9, with a mean of 44.99. The percent hypodiploidy ranged from 0-100% and the percent hyperdiploidy from 0-53%. Detailed cytogenetic analyses were very difficult to perform because of the paucity of mitotic figures and the poor chromosome morphology. Thus, detailed chromosome banding analysis on these patients was impossible.^ Thirty seven percent of the patients had normal total DNA content, whereas 63% had abnormal amounts of DNA (one patient with less than normal amounts and 25 patients with greater than normal amounts of DNA).^ Several clinical parameters were used in the statistical analyses: tumor burden, patient status at biopsy, patient response status, past therapy, type of treatment and percent plasma cells. Only among these clinical parameters were any statistically significant correlations found: pretreatment tumor burden versus patient response, patient biopsy status versus patient response and past therapy versus patient response.^ No correlations were found between percent hypodiploid, diploid, hyperdiploid or DNA content, and the patient response status, nor were any found between those patients with: (a) normal plasma cells, low pretreatment tumor mass burden and more than 50% of the analyzed metaphases with 46 chromosomes; (b) normal amounts of DNA, low pretreatment tumor mass burden and more than 50% of the metaphases with 46 chromosomes; (c) normal amounts of DNA and normal quantities of plasma cells; (d) abnormal amounts of DNA, abnormal amounts of plasma cells, high pretreatment tumor mass burden and less than 50% of the metaphases with 46 chromosomes.^ Technical drawbacks of both cytogenetic and DNA content analysis in these multiple myeloma patients are discussed along with the lack of correlations between DNA content and chromosome number. Refined chromosome banding analysis awaits technical improvements before we can understand which chromosome material (if any) makes up the "extra" amounts of DNA in these patients. None of the correlations tested can be used as diagnostic or prognostic aids for multiple myeloma. ^

STUDIES OF THE MOLECULAR MECHANISMS OF CHROMOSOME ABERRATION FORMATION (DNA REPAIR, CROSSLINKS, MITOMYCIN C)

The objective of this research has been to study the molecular basis for chromosome aberration formation. Predicated on a variety of data, Mitomycin C (MMC)-induced DNA damage has been postulated to cause the formation of chromatid breaks (and gaps) by preventing the replication of regions of the genome prior to mitosis. The basic protocol for these experiments involved treating synchronized Hela cells in G(,1)-phase with a 1 (mu)g/ml dose of MMC for one hour. After removing the drug, cells were then allowed to progress to mitosis and were harvested for analysis by selective detachment. Utilizing the alkaline elution assay for DNA damage, evidence was obtained to support the conclusion that Hela cells can progress through S-phase into mitosis with intact DNA-DNA interstrand crosslinks. A higher level of crosslinking was observed in those cells remaining in interphase compared to those able to reach mitosis at the time of analysis. Dual radioisotope labeling experiments revealed that, at this dose, these crosslinks were associated to the same extent with both parental and newly replicated DNA. This finding was shown not to be the result of a two-step crosslink formation mechanism in which crosslink levels increase with time after drug treatment. It was also shown not to be an artefact of the double-labeling protocol. Using neutral CsCl density gradient ultracentrifugation of mitotic cells containing BrdU-labeled newly replicated DNA, control cells exhibited one major peak at a heavy/light density. However, MMC-treated cells had this same major peak at the heavy/light density, in addition to another minor peak at a density characteristic for light/light DNA. This was interpreted as indicating either: (1) that some parental DNA had not been replicated in the MMC treated sample or; (2) that a recombination repair mechanism was operational. To distinguish between these two possibilities, flow cytometric DNA fluorescence (i.e., DNA content) measurements of MMC-treated and control cells were made. These studies revealed that the mitotic cells that had been treated with MMC while in G(,1)-phase displayed a 10-20% lower DNA content than untreated control cells when measured under conditions that neutralize chromosome condensation effects (i.e., hypotonic treatment). These measurements were made under conditions in which the binding of the drug, MMC, was shown not to interfere with the stoichiometry of the ethidium bromide-mithramycin stain. At the chromosome level, differential staining techniques were used in an attempt to visualize unreplicated regions of the genome, but staining indicative of large unreplicated regions was not observed. These results are best explained by a recombinogenic mechanism. A model consistent with these results has been proposed.^