Compositional differences within and between eukaryotic genomes

Eukaryotic genome similarity relationships are inferred using sequence information derived from large aggregates of genomic sequences. Comparisons within and between species sample sequences are based on the profile of dinucleotide relative abundance values (The profile is ρ*XY = f*XY/f*Xf*Y for all XY, where f*X denotes the frequency of the nucleotide X and f*XY denotes the frequency of the dinucleotide XY, both computed from the sequence concatenated with its inverted complement). Previous studies with respect to prokaryotes and this study document that profiles of different DNA sequence samples (sample size ≥50 kb) from the same organism are generally much more similar to each other than they are to profiles from other organisms, and that closely related organisms generally have more similar profiles than do distantly related organisms. On this basis we refer to the collection {ρ*XY} as the genome signature. This paper identifies ρ*XY extremes and compares genome signature differences for a diverse range of eukaryotic species. Interpretations on the mechanisms maintaining these profile differences center on genome-wide replication, repair, DNA structures, and context-dependent mutational biases. It is also observed that mitochondrial genome signature differences between species parallel the corresponding nuclear genome signature differences despite large differences between corresponding mitochondrial and nuclear signatures. The genome signature differences also have implications for contrasts between rodents and other mammals, and between monocot and dicot plants, as well as providing evidence for similarities among fungi and the diversity of protists.

Phylogenetic assessment of length variation at a microsatellite locus

Sixty-six haplotypes at a locus containing a simple dinucleotide (CA)n microsatellite repeat were isolated by PCR–single-strand conformational polymorphism from populations of the horseshoe crab Limulus polyphemus. These haplotypes were sequenced to assess nucleotide variation directly. Thirty-four distinct sequences (alleles) were identified in a region 570 bp long that included the microsatellite motif. In the repeat region itself, CA-number varied in integer values from 5 to 11 across alleles, except that a (CA)8 class was not observed. Differences among alleles were due also to polymorphisms at 22 sites in regions immediately flanking the microsatellite repeats. Nucleotide substitutions in these regions were used to estimate phylogenetic relationships among alleles, and the gene phylogeny was used to trace the evolution of length variation and CA repeat numbers. A low correlation between size variation and genealogical relationships among alleles suggests that absolute fragment size (as normally scored in microsatellite assays) is an unreliable indicator of historical affinities among alleles. This finding on the molecular fine structure of microsatellite variation suggests the need for caution in the use of repeat counts at microsatellite loci as secure indicators of allelic relationships.

Alternative genetic pathways in colorectal carcinogenesis

The comparative typing of matched tumor and blood DNAs at dinucleotide repeat (microsatellite) loci has revealed in tumor DNA the presence of alleles that are not observed in normal DNA. The occurrence of these additional alleles is possibly due to replication errors (RERs). Although this observation has led to the recognition of a subtype of colorectal cancer with a high incidence of RERs (caused by a deficiency in DNA mismatch repair), a thorough analysis of the RER frequency in a consecutive series of colorectal cancers had not been reported. It is shown here that the extensive typing of 88 colorectal tumors reveals a bimodal distribution for the frequency of RER at microsatellite loci. Within the major mode (75 tumors, RER− subtype), the probability that a locus exhibited instability did not differ significantly among loci and tumors, being 0.02. The subsequent development of a statistical test for an operational discrimination between the RER− and RER+ subtypes indicated that the probability of misclassification did not exceed 0.001 in this series. The frequency of K-ras mutation was found to be equivalent in the two subtypes. However, in the RER+ tumors, the p53 gene mutation was less frequently detected, the adenomatous polyposis coli (APC) mutation was rare, and the biallelic inactivation of either of these genes was not observed. Furthermore, the concomitant occurrence of APC and tumor growth factor β receptor type II gene alterations was found only once. These data suggest that the repertoires of genes that are frequently altered in RER+ and RER− tumors may be more different than previously thought.

Gene cloning, sequence analysis, and expression of 2-methyl-3-hydroxypyridine-5-carboxylic acid oxygenase

The gene encoding 2-methyl-3-hydroxypyridine-5-carboxylic acid oxygenase (MHPCO; EC 1.14.12.4) was cloned by using an oligonucleotide probe corresponding to the N terminus of the enzyme to screen a DNA library of Pseudomonas sp. MA-1. The gene encodes for a protein of 379 amino acid residues corresponding to a molecular mass of 41.7 kDa, the same as that previously estimated for MHPCO. MHPCO was expressed in Escherichia coli and found to have the same properties as the native enzyme from Pseudomonas sp. MA-1. This study shows that MHPCO is a homotetrameric protein with one flavin adenine dinucleotide bound per subunit. Sequence comparison of the enzyme with other hydroxylases reveals regions that are conserved among aromatic flavoprotein hydroxylases.

Poly(ADP-ribose) polymerase is a mediator of necrotic cell death by ATP depletion

Apoptotic and necrotic cell death are well characterized and are influenced by intracellular ATP levels. Poly(ADP-ribose) polymerase (PARP), a nuclear enzyme activated by DNA strand breaks, physiologically participates in DNA repair. Overactivation of PARP after cellular insults can lead to cell death caused by depletion of the enzyme’s substrate β-nicotinamide adenine dinucleotide and of ATP. In this study, we have differentially elicited apoptosis or necrosis in mouse fibroblasts. Fibroblasts from PARP-deficient (PARP−/−) mice are protected from necrotic cell death and ATP depletion but not from apoptotic death. These findings, together with cell death patterns in PARP−/− animals receiving other types of insults, indicate that PARP activation is an active trigger of necrosis, whereas other mechanisms mediate apoptosis.

Enhancement of DNA repair in human skin cells by thymidine dinucleotides: Evidence for a p53-mediated mammalian SOS response

Thymidine dinucleotide (pTpT) stimulates melanogenesis in mammalian pigment cells and intact skin, mimicking the effects of UV irradiation and UV-mimetic DNA damage. Here it is shown that, in addition to tanning, pTpT induces a second photoprotective response, enhanced repair of UV-induced DNA damage. This enhanced repair results in a 2-fold increase in expression of a UV-damaged chloramphenicol acetyltransferase expression vector transfected into pTpT-treated skin fibroblasts and keratinocytes, compared with diluent-treated cells. Direct measurement of thymine dimers and (6–4) photoproducts by immunoassay demonstrates faster repair of both of these UV-induced photoproducts in pTpT-treated fibroblasts. This enhanced repair capacity also improves cell survival and colony-forming ability after irradiation. These effects of pTpT are accomplished, at least in part, by the up-regulation of a set of genes involved in DNA repair (ERCC3 and GADD45) and cell cycle inhibition (SDI1). At least two of these genes (GADD45 and SDI1) are known to be transcriptionally regulated by the p53 tumor suppressor protein. Here we show that pTpT activates p53, leading to nuclear accumulation of this protein, and also increases the specific binding of this transcription factor to its DNA consensus sequence.

Three-dimensional structure of human electron transfer flavoprotein to 2.1-Å resolution

Mammalian electron transfer flavoproteins (ETF) are heterodimers containing a single equivalent of flavin adenine dinucleotide (FAD). They function as electron shuttles between primary flavoprotein dehydrogenases involved in mitochondrial fatty acid and amino acid catabolism and the membrane-bound electron transfer flavoprotein ubiquinone oxidoreductase. The structure of human ETF solved to 2.1-Å resolution reveals that the ETF molecule is comprised of three distinct domains: two domains are contributed by the α subunit and the third domain is made up entirely by the β subunit. The N-terminal portion of the α subunit and the majority of the β subunit have identical polypeptide folds, in the absence of any sequence homology. FAD lies in a cleft between the two subunits, with most of the FAD molecule residing in the C-terminal portion of the α subunit. Alignment of all the known sequences for the ETF α subunits together with the putative FixB gene product shows that the residues directly involved in FAD binding are conserved. A hydrogen bond is formed between the N5 of the FAD isoalloxazine ring and the hydroxyl side chain of αT266, suggesting why the pathogenic mutation, αT266M, affects ETF activity in patients with glutaric acidemia type II. Hydrogen bonds between the 4′-hydroxyl of the ribityl chain of FAD and N1 of the isoalloxazine ring, and between αH286 and the C2-carbonyl oxygen of the isoalloxazine ring, may play a role in the stabilization of the anionic semiquinone. With the known structure of medium chain acyl-CoA dehydrogenase, we hypothesize a possible structure for docking the two proteins.

Mutation rate varies among alleles at a microsatellite locus: Phylogenetic evidence

The understanding of the mutational mechanism that generates high levels of variation at microsatellite loci lags far behind the application of these genetic markers. A phylogenetic approach was developed to study the pattern and rate of mutations at a dinucleotide microsatellite locus tightly linked to HLA-DQB1 (DQCAR). A random Japanese population (n = 129) and a collection of multiethnic samples (n = 941) were typed at the DQB1 and DQCAR loci. The phylogeny of DQB1 alleles was then reconstructed and DQCAR alleles were superimposed onto the phylogeny. This approach allowed us to group DQCAR alleles that share a common ancestor. The results indicated that the DQCAR mutation rate varies drastically among alleles within this single microsatellite locus. Some DQCAR alleles never mutated during a long period of evolutionary time. Sequencing of representative DQCAR alleles showed that these alleles lost their ability to mutate because of nucleotide substitutions that shorten the length of uninterrupted CA repeat arrays; in contrast, all mutating alleles had relatively longer perfect CA repeat sequences.

Role of adenine functional groups in the recognition of the 3′-splice-site AG during the second step of pre-mRNA splicing

The AG dinucleotide at the 3′ splice sites of metazoan nuclear pre-mRNAs plays a critical role in catalytic step II of the splicing reaction. Previous studies have shown that replacement of the guanine by adenine in the AG (AG → GG) inhibits this step. We find that the second step was even more severely inhibited by cytosine (AG → CG) or uracil (AG → UG) substitutions at this position. By contrast, a relatively moderate inhibition was observed with a hypoxanthine substitution (AG → HG). When adenine was replaced by a purine base (AG → PG) or by 7-deazaadenine (AG → c7AG), little effect on the second step was observed, suggesting that the 6-NH2 and N7 groups do not play a critical role in adenine recognition. Finally, replacement of adenine by 2-aminopurine (AG → 2-APG) had no effect on the second step. Taken together, our results suggest that the N1 group of adenine functions as an essential determinant in adenine recognition during the second step of pre-mRNA splicing.

DNA methylation is a reversible biological signal

The pattern of DNA methylation plays an important role in regulating different genome functions. To test the hypothesis that DNA methylation is a reversible biochemical process, we purified a DNA demethylase from human cells that catalyzes the cleavage of a methyl residue from 5-methyl cytosine and its release as methanol. We show that similar to DNA methyltransferase, DNA demethylase shows CpG dinucleotide specificity, can demethylate mdCpdG sites in different sequence contexts, and demethylates both fully methylated and hemimethylated DNA. Thus, contrary to the commonly accepted model, DNA methylation is a reversible signal, similar to other physiological biochemical modifications.

Crystal structures of bovine milk xanthine dehydrogenase and xanthine oxidase: Structure-based mechanism of conversion

Mammalian xanthine oxidoreductases, which catalyze the last two steps in the formation of urate, are synthesized as the dehydrogenase form xanthine dehydrogenase (XDH) but can be readily converted to the oxidase form xanthine oxidase (XO) by oxidation of sulfhydryl residues or by proteolysis. Here, we present the crystal structure of the dimeric (Mr, 290,000) bovine milk XDH at 2.1-Å resolution and XO at 2.5-Å resolution and describe the major changes that occur on the proteolytic transformation of XDH to the XO form. Each molecule is composed of an N-terminal 20-kDa domain containing two iron sulfur centers, a central 40-kDa flavin adenine dinucleotide domain, and a C-terminal 85-kDa molybdopterin-binding domain with the four redox centers aligned in an almost linear fashion. Cleavage of surface-exposed loops of XDH causes major structural rearrangement of another loop close to the flavin ring (Gln 423—Lys 433). This movement partially blocks access of the NAD substrate to the flavin adenine dinucleotide cofactor and changes the electrostatic environment of the active site, reflecting the switch of substrate specificity observed for the two forms of this enzyme.

Structural code for DNA recognition revealed in crystal structures of papillomavirus E2-DNA targets

Transcriptional regulation in papillomaviruses depends on sequence-specific binding of the regulatory protein E2 to several sites in the viral genome. Crystal structures of bovine papillomavirus E2 DNA targets reveal a conformational variant of B-DNA characterized by a roll-induced writhe and helical repeat of 10.5 bp per turn. A comparison between the free and the protein-bound DNA demonstrates that the intrinsic structure of the DNA regions contacted directly by the protein and the deformability of the DNA region that is not contacted by the protein are critical for sequence-specific protein/DNA recognition and hence for gene-regulatory signals in the viral system. We show that the selection of dinucleotide or longer segments with appropriate conformational characteristics, when positioned at correct intervals along the DNA helix, can constitute a structural code for DNA recognition by regulatory proteins. This structural code facilitates the formation of a complementary protein–DNA interface that can be further specified by hydrogen bonds and nonpolar interactions between the protein amino acids and the DNA bases.

A terbenzimidazole that preferentially binds and conformationally alters structurally distinct DNA duplex domains: A potential mechanism for topoisomerase I poisoning

The terbenzimidazoles are a class of synthetic ligands that poison the human topoisomerase I (TOP1) enzyme and promote cancer cell death. It has been proposed that drugs of this class act as TOP1 poisons by binding to the minor groove of the DNA substrate of TOP1 and altering its structure in a manner that results in enzyme-mediated DNA cleavage. To test this hypothesis, we characterize and compare the binding properties of a 5-phenylterbenzimidazole derivative (5PTB) to the d(GA4T4C)2 and d(GT4A4C)2 duplexes. The d(GA4T4C)2 duplex contains an uninterrupted 8-bp A⋅T domain, which, on the basis of x-ray crystallographic data, should induce a highly hydrated “A-tract” conformation. This duplex also exhibits anomalously slow migration in a polyacrylamide gel, a feature characteristic of a noncanonical global conformational state frequently described as “bent.” By contrast, the d(GT4A4C)2 duplex contains two 4-bp A⋅T tracts separated by a TpA dinucleotide step, which should induce a less hydrated “B-like” conformation. This duplex also migrates normally in a polyacrylamide gel, a feature further characteristic of a global, canonical B-form duplex. Our data reveal that, at 20°C, 5PTB exhibits an ≈2.3 kcal/mol greater affinity for the d(GA4T4C)2 duplex than for the d(GT4A4C)2 duplex. Significantly, we find this sequence/conformational binding specificity of 5PTB to be entropic in origin, an observation consistent with a greater degree of drug binding-induced dehydration of the more solvated d(GA4T4C)2 duplex. By contrast with the differential duplex affinity exhibited by 5PTB, netropsin and 4′,6-diamidino-2-phenylindole (DAPI), two AT-specific minor groove binding ligands that are inactive as human TOP1 poisons, bind to both duplexes with similar affinities. The electrophoretic behaviors of the ligand-free and ligand-bound duplexes are consistent with 5PTB-induced bending and/or unwinding of both duplexes, which, for the d(GA4T4C)2 duplex, is synergistic with the endogenous sequence-directed electrophoretic properties of the ligand-free duplex state. By contrast, the binding to either duplex of netropsin or DAPI induces little or no change in the electrophoretic mobilities of the duplexes. Our results demonstrate that the TOP1 poison 5PTB binds differentially to and alters the structures of the two duplexes, in contrast to netropsin and DAPI, which bind with similar affinities to the two duplexes and do not significantly alter their structures. These results are consistent with a mechanism for TOP1 poisoning in which drugs such as 5PTB differentially target conformationally distinct DNA sites and induce structural changes that promote enzyme-mediated DNA cleavage.

Methylation of the minimal promoter of an embryonic globin gene silences transcription in primary erythroid cells

Methylation of cytosines in the dinucleotide CpG has been shown to suppress transcription of a number of tissue-specific genes, yet the precise mechanism is not fully understood. The vertebrate globin genes were among the first examples in which an inverse correlation was shown between CpG methylation and transcription. We studied the methylation pattern of the 235-bp ρ-globin gene promoter in genomic DNA from primary chicken erythroid cells using the sodium bisulfite conversion technique and found all CpGs in the promoter to be methylated in erythroid cells from adult chickens in which the ρ-globin gene is silent but unmethylated in 5-day (primitive) embryonic red cells in which the gene is transcribed. To elucidate further the mechanism of methylation-induced silencing, an expression construct consisting of 235 bp of 5′ promoter sequence of the ρ-globin gene along with a strong 5′ erythroid enhancer driving a chloramphenicol acetyltransferase reporter gene, ρ-CAT, was transfected into primary avian erythroid cells derived from 5-day embryos. Methylation of just the 235-bp ρ-globin gene promoter fragment at every CpG resulted in a 20- to 30-fold inhibition of transcription, and this effect was not overridden by the presence of potent erythroid-specific enhancers. The ability of the 235-bp ρ-globin gene promoter to bind to a DNA Methyl Cytosine binding Protein Complex (MeCPC) was tested in electrophoretic mobility shift assays utilizing primary avian erythroid cell nuclear extract. The results were that fully methylated but not unmethylated 235-bp ρ-globin gene promoter fragment could compete efficiently for MeCPC binding. These results are a direct demonstration that site-specific methylation of a globin gene promoter at the exact CpGs that are methylated in vivo can silence transcription in homologous primary erythroid cells. Further, these data implicate binding of MeCPC to the promoter in the mechanism of silencing.

Replacement of Fhit in cancer cells suppresses tumorigenicity

The candidate tumor suppressor gene, FHIT, encompasses the common human chromosomal fragile site at 3p14.2, the hereditary renal cancer translocation breakpoint, and cancer cell homozygous deletions. Fhit hydrolyzes dinucleotide 5′,5‴-P1,P3-triphosphate in vitro and mutation of a central histidine abolishes hydrolase activity. To study Fhit function, wild-type and mutant FHIT genes were transfected into cancer cell lines that lacked endogenous Fhit. No consistent effect of exogenous Fhit on growth in culture was observed, but Fhit and hydrolase “dead” Fhit mutant proteins suppressed tumorigenicity in nude mice, indicating that 5′,5‴-P1,P3-triphosphate hydrolysis is not required for tumor suppression.