Epigenetic patterns within the haplotype phased fig (Ficus carica L.) genome

Due to DNA heterozygosity and repeat content, assembly of non‐model plant genomes is challenging. Herein, we report a high‐quality genome reference of one of the oldest known domesticated species, fig (Ficus carica L.), using Pacific Biosciences single‐molecule, real‐time sequencing. The fig genome is ~333 Mbp in size, of which 80% has been anchored to 13 chromosomes. Genome‐wide analysis of N⁶‐methyladenine and N⁴‐methylcytosine revealed high methylation levels in both genes and transposable elements, and a prevalence of methylated over non‐methylated genes. Furthermore, the characterization of N⁶‐methyladenine sites led to the identification of ANHGA, a species‐specific motif, which is prevalent for both genes and transposable elements. Finally, exploiting the contiguity of the 13 pseudomolecules, we identified 13 putative centromeric regions. The high‐quality reference genome and the characterization of methylation profiles, provides an important resource for both fig breeding and for fundamental research into the relationship between epigenetic changes and phenotype, using fig as a model species.     

Interactions of Quinoxaline Antibiotic and DNA.: The Molecular Structure of a Triostin A—d(GCGTACGC) Complex

Abstract

The crystal structure of a DNA. octamer d(GCGTA.CGC) complexed to an antitumor antibiotic, triostin A, has been solved and refined to 2.2 Å resolution by x-ray diffraction analysis. The antibiotic molecule acts as a true bis intercalator surrouding the d(CpG) sequence at either end of the unwound right-handed DNA. double helix. A.s previously observed in the structure of triostin A.—d(CGTA.CG) complex (A.H.-J. Wang, et. al., Science, 225,1115–1121 (1984)), the alanine amino acid residues of the drug molecule form sequence-specific hydrogen bonds to guanines in the minor groove. The two central A · T base pairs are in Hoogsteen configuration with adenine in the syn conformation. In addition, the two terminal G · C base pairs flanking the quinoxaline rings are also held together by Hoogsteen base pairing. This is the first observation in an oligonucleotide of. Hoogsteen G · C base pairs where the cytosine is protonated. The principal functional components of a bis-intercalative compound are discussed.     

New information content in RNA base pairing deduced from quantitative analysis of high-resolution structures

Non-canonical base pairs play important roles in organizing the complex three-dimensional folding of RNA. Here, we outline methodology developed both to analyze the spatial patterns of interacting base pairs in known RNA structures and to reconstruct models from the collective experimental information. We focus attention on the structural context and deformability of the seven pairing patterns found in greatest abundance in the helical segments in a set of well-resolved crystal structures, including (i–ii) the canonical A·U and G·C Watson–Crick base pairs, (iii) the G·U wobble pair, (iv) the sheared G·A pair, (v) the A·U Hoogsteen pair, (vi) the U·U wobble pair, and (vii) the G·A Watson–Crick-like pair. The non-canonical pairs stand out from the canonical associations in terms of apparent deformability, spanning a broader range of conformational states as measured by the six rigid-body parameters used to describe the spatial arrangements of the interacting bases, the root-mean-square deviations of the base-pair atoms, and the fluctuations in hydrogen-bonding geometry. The deformabilties, the modes of base-pair deformation, and the preferred sites of occurrence depend on sequence. We also characterize the positioning and overlap of the base pairs with respect to the base pairs that stack immediately above and below them in double-helical fragments. We incorporate the observed positions of the bases, base pairs, and intervening phosphorus atoms in models to predict the effects of the non-canonical interactions on overall helical structure.     

Metastable secondary structures in ribosomal RNA: a new method for analyzing the titration behavior of rRNA

The pH titration behavior of E. coli rRNA in the acid range has been analyzed by combining spectrophotometric and potentiometric titration data. The “simplest” model for the system, which considers as possible reactions the protonation of adenine (A), cytosine (C), and guanine (G) residues along with the opening of A·U and G·C base pairs, does not adequately account for the titration properties. It is postulated that extra reactions may occur in addition to those in the “simplest” model, and a new analytical method was developed to deal with this situation. Our approach yields the ultraviolet spectral changes which accompany the extra reactions, from which the nature of these reactions can in principle be deduced. The calculations also give, at each pH, the extents of the extra reactions as well as the extents of those reactions which comprise the “simplest” model. We infer that in acidic RNA solutions of 0.1M ionic strength there occur at least two extra reactions, each of which involves G residues. We propose that in the pH range 6.0 ≥ pH ≥ 3.8 triple-stranded helical sequences, presumably protonated G·C·G, are formed. These regions are replaced at lower pH by acid-stable structures involving G·G and A·A base pairs. In solutions of lower ionic strength (I = 0.01M) no triple strands are formed, but G·G and A·A regions seem to develop even at pH values as high as 6.0. At I = 0.1M, an acid–base titration cycle between pH 7 and 2.8 is not reversible; rRNA shows true hysteresis behavior. We conclude that in ribosomal RNA's, which are generally G-rich, guanine residues may participate in hitherto unpredicted conformations, some of which may be metastable while others are equilibrium structures.     

Effects of N2,N2-dimethylguanosine on RNA structure and stability: crystal structure of an RNA duplex with tandem m2 2G:A pairs

Methylation of the exocyclic amino group of guanine is a relatively common modification in rRNA and tRNA. Single methylation (N²-methylguanosine, m²G) is the second most frequently encountered nucleoside analog in Escherichia coli rRNAs. The most prominent case of dual methylation (N²,N²-dimethylguanosine, m²₂G) is found in the majority of eukaryotic tRNAs at base pair m²₂G26:A44. The latter modification eliminates the ability of the N² function to donate in hydrogen bonds and alters its pairing behavior, notably vis-à-vis C. Perhaps a less obvious consequence of the N²,N²-dimethyl modification is its role in controlling the pairing modes between G and A. We have determined the crystal structure of a 13-mer RNA duplex with central tandem m²₂G:A pairs. In the structure both pairs adopt an imino-hydrogen bonded, pseudo-Watson–Crick conformation. Thus, the sheared conformation frequently seen in tandem G:A pairs is avoided due to a potential steric clash between an N²-methyl group and the major groove edge of A. Additionally, for a series of G:A containing self-complementary RNAs we investigated how methylation affects competitive hairpin versus duplex formation based on UV melting profile analysis.     

Structure determination of DNA methylation lesions N1-meA and N3-meC in duplex DNA using a cross-linked protein–DNA system

N¹-meA and N³-meC are cytotoxic DNA base methylation lesions that can accumulate in the genomes of various organisms in the presence of S_N2 type methylating agents. We report here the structural characterization of these base lesions in duplex DNA using a cross-linked protein–DNA crystallization system. The crystal structure of N¹-meA:T pair shows an unambiguous Hoogsteen base pair with a syn conformation adopted by N¹-meA, which exhibits significant changes in the opening, roll and twist angles as compared to the normal A:T base pair. Unlike N¹-meA, N³-meC does not establish any interaction with the opposite G, but remains partially intrahelical. Also, structurally characterized is the N⁶-meA base modification that forms a normal base pair with the opposite T in duplex DNA. Structural characterization of these base methylation modifications provides molecular level information on how they affect the overall structure of duplex DNA. In addition, the base pairs containing N¹-meA or N³-meC do not share any specific characteristic properties except that both lesions create thermodynamically unstable regions in a duplex DNA, a property that may be explored by the repair proteins to locate these lesions.     

Human METTL16 is a N6‐methyladenosine (m6A) methyltransferase that targets pre‐mRNAs and various non‐coding RNAs

N⁶‐methyladenosine (m⁶A) is a highly dynamic RNA modification that has recently emerged as a key regulator of gene expression. While many m⁶A modifications are installed by the METTL3–METTL14 complex, others appear to be introduced independently, implying that additional human m⁶A methyltransferases remain to be identified. Using crosslinking and analysis of cDNA (CRAC), we reveal that the putative human m⁶A “writer” protein METTL16 binds to the U6 snRNA and other ncRNAs as well as numerous lncRNAs and pre‐mRNAs. We demonstrate that METTL16 is responsible for N⁶‐methylation of A43 of the U6 snRNA and identify the early U6 biogenesis factors La, LARP7 and the methylphosphate capping enzyme MEPCE as METTL16 interaction partners. Interestingly, A43 lies within an essential ACAGAGA box of U6 that base pairs with 5′ splice sites of pre‐mRNAs during splicing, suggesting that METTL16‐mediated modification of this site plays an important role in splicing regulation. The identification of METTL16 as an active m⁶A methyltransferase in human cells expands our understanding of the mechanisms by which the m⁶A landscape is installed on cellular RNAs.     

Sequence specificity of the P1 modification methylase (M.Eco P1) and the DNA methylase (M.Eco dam) controlled by the Escherichia coli dam gene.

Labeled oligonucleotides have been fractionated from pancreatic DNase digests of DNA that had been methylated in vitro with the P1 modification enzyme (M·Eco P1) or with the DNA-adenine methylase (M·Eco dam) controlled by the Escherichia coli dam gene. The sequences of methylated oligonucleotides were established for M·Eco dam modification of calf thymus DNA. The results show that M·Eco dam inethylates adenine residues contained in the twofold symmetrical sequence, 5′ … G-A-T-C … 3′. The sequence for the site methylated by M·Eco P1 has also been deduced; we propose that M·Eco P1 modification produces the following methylated pentameric sequence: 5′ … A-G-A¹-C-Py … 3′ (where $\text{A}{}^1\text{= N}{}^6$ methyladenine and Py is C or T).     

Structure-activity relationship of cytokinins: crystal structure and conformation of 6-furfurylaminopurine (kinetin).

The crystal structure of 6-furfurylaminopurine (kinetin) was determined from three-dimensional x-ray diffraction data. The N⁶-substituent is distal to the imidazole ring. The molecules are linked across centers of inversion by pairs of N(6)-H···N(7) and N(9)-H···N(3) hydrogen bonds, utilizing the Hoogsteen sites for base pairing. From an analysis of the conformational preferences of cytokinins, an “active conformation” favourable for eliciting maximal cytokinin activity is proposed. The loss of cytokinin activity due to various chemical modifications such as the substitution at N¹, alteration of the methylene bridge, is correlated to the inability of cytokinins to achieve the active conformation.     

Sequence and conformational effects on imino proton exchange in A.T- and A.U-containing DNA and RNA duplexes

¹H-nmr relaxation has been used to study the effect of sequence and conformation on imino proton exchange in adenine–thymine (A · T) and adenine–uracil (A · U) containing DNA and RNA duplexes. At low temperature, relaxation is caused by dipolar interactions between the imino and the adenine amino and AH2 protons, and at higher temperature, by exchange with the solvent protons. Although room temperature exchange rates vary between 3 and 12s^?1, the exchange activation energies (E_α) are insensitive to changes in the duplex sequence (alternating vs homopolymer duplexes), the conformation (B-form DNA vs A-form RNA), and the identity of the pyrimidine base (thymine vs uracil). The average value of the activation energy for the five duplexes studied, poly[d(A-T)], poly[d(A) · d(T)], poly[d(A-U)], Poly[d(A) · d(U)], and poly[r(A) · r(U)], was 16.8 ± 1.3 kcal/mol. In addition, we find that the average E_α for the A.T base pairs in a 43-base-pair restriction fragment is 16.4 ± 1.0 kcal/mol. This result is to be contrasted with the observation that the E_α of cytosine-containing duplexes depends on the sequence, conformation, and substituent groups on the purine and pyrimidine bases. Taken together, the data indicate that there is a common low-energy pathway for the escape of the thymine (uracil) imino protons from the double helix. The absolute values of the exchange rates in the simple sequence polymers are typically 3–10 times faster than in DNAs containing both A · T and G · C base pairs.     

DNA methylation in vegetative and conjugating cells of a protozoan ciliate: Blepharisma japonicum

We report here the presence of N⁶-methyladenine (MeAde) in the macronuclear DNA (maDNA) of Blepharisma japonicum vegetative cells. We have further investigated the relationship between DNA methylation and cell union in cells activated for conjugation. Such activation was induced by treating cells of mating type I with complementary gamone 2. We found a reduction of about 24% of MeAde content in gamone-treated cells ready for cell union. First indications of the presence and reduction of MeAde content came from electrophoresis of maDNA digested by appropriate restriction endonucleases. Chromatographic determination of the amount of methylated base by HPLC substantiated these observations. In vegetative cells, 1.576 ± 0.02% of total adenine was found to be methylated as opposed to 1.193 ± 0.04% in activated cells. The HPLC analysis of maDNA also revealed a peak with a retention time corresponding to that of 5-hydroxymethyluracil, already found in some species of dinoflagellates. In that gamone treatment is correlated with a differential gene expression (indicated by a differential RNA and protein synthesis), our results suggest that there is a relationship between macronuclear genome activation and demethylation of maDNA. This is the first report of a correlation between gene activation and adenine demethylation in a eukaryotic organism.     

Parallel nucleic acid helices with hoogsteen base pairing: Symmetry and structure

Molecular structures for parallel DNA and RNA double helices with Hoogsteen pairing are proposed for the first time. The DNA helices have sugars in the C2′-endo region and the phosphodiester conformations are (trans, gauche^?), and the RNA helices have sugars in the C3′-endo region and the phosphodiester conformations are (gauche^?, gauche^?). A pseudorotational symmetry relates the two parallel strands of DNA helices and a screw symmetry relates the two strands of RNA helices, which have an associated tilt of the The conformational space of parallel helices with Hoogsteen base pairing, unlike the Watson-Crick duplex, is highly restricted due to the unique positioning of the symmetry axis in the former case. The features of the parallel double helix with Hoogsteen pairing are compared with the Watson-Crick duplex and the corresponding triple helix. © 1994 John Wiley & Sons, Inc.     

Reactions of plant copper/topaquinone amine oxidases with N6-aminoalkyl derivatives of adenine

Plant copper/topaquinone-containing amine oxidases (CAOs, EC 1.4.3.6) are enzymes oxidising various amines. Here we report a study on the reactions of CAOs from grass pea (Lathyrus sativus), lentil (Lens esculenta) and Euphorbia characias, a Mediterranean shrub, with N⁶-aminoalkyl adenines representing combined analogues of cytokinins and polyamines. The following compounds were synthesised: N⁶-(3-aminopropyl)adenine, N⁶-(4-aminobutyl)adenine, N⁶-(4-amino-trans-but-2-enyl)adenine, N⁶-(4-amino-cis-but-2-enyl)adenine and N⁶-(4-aminobut-2-ynyl)adenine. From these, N⁶-(4-aminobutyl)adenine and N⁶-(4-amino-trans-but-2-enyl)adenine were found to be substrates for all three enzymes (K_m～10^?4?M). Absorption spectroscopy demonstrated such an interaction with the cofactor topaquinone, which is typical for common diamine substrates. However, only the former compound provided a regular reaction stoichiometry. Anaerobic absorption spectra of N⁶-(3-aminopropyl)adenine, N⁶-(4-amino-cis-but-2-enyl)adenine and N⁶-(4-aminobut-2-ynyl)adenine reactions revealed a similar kind of initial interaction, although the compounds finally inhibited the enzymes. Kinetic measurements allowed the determination of both inhibition type and strength; N⁶-(3-aminopropyl)adenine and N⁶-(4-amino-cis-but-2-enyl)adenine produced reversible inhibition (K_i～10^?5–10^?4?M) whereas, N⁶-(4-aminobut-2-ynyl)adenine could be considered a powerful inactivator.     

Methylation of DNA in various organs of c57b1 mice by a carcinogenic dose of N-methyl-N-nitrosourea and stability of some methylation products up to 18 hours

Female 6–8-week-old C57B1 mice were injected i.p. with N-methyl-N-nitrosourea (MNUA) (¹⁴C or ³H-methyl-labelled) in saline (80 mg/kg) and DNA was isolated from bone marrow, small bowel, kidneys, liver, lungs, spleen and thymus at various times thereafter up to 18 h. Methylation of DNA was found in all organs examined, and by analyses using column or paper chromatography of DNA hydrolysates, the extent of methylation of DNA purines was determined.Methylated guanine residues (at N-3, N-7 and O-6 positions) were stable in DNA up to 18 h, but methylated adenines (at N-3 or N-7) were removed from DNA of all organs examined; the overall half-life of methyladenines was about 3 h, but removal appeared to occur in a biphasic manner, with a proportion of methyladenine remaining relatively stable. This relative stability was somewhat more marked in bone marrow than in other organs.     

Molecular dynamics simulations and coupled nucleotide substitution experiments indicate the nature of A·A base pairing and a putative structure of the coralyne-induced homo-adenine duplex

Coralyne is an alkaloid drug that binds homo-adenine DNA (and RNA) oligonucleotides more tightly than it does Watson–Crick DNA. Hud’s laboratory has shown that poly(dA) in the presence of coralyne forms an anti-parallel duplex, however attempts to determine the structure by NMR spectroscopy and X-ray crystallography have been unsuccessful. Assuming adenine–adenine hydrogen bonding between the two poly(dA) strands, we constructed 40 hypothetical homo-(dA) anti-parallel duplexes and docked coralyne into the six most favorable duplex structures. The two most stable structures had trans glycosidic bonds, but distinct pairing geometries, i.e. either Watson–Crick Hoogsteen (transWH) or Watson–Crick Watson–Crick (transWW) with stability of transWH > transWW. To narrow down the possibilities, 7-deaza adenine base substitutions (dA→7) were engineered into homo-(dA) sequences. These substitutions significantly reduced the thermal stability of the coralyne-induced homo-(dA) structure. These experiments strongly suggest the involvement of N7 in the coralyne-induced A·A base pairs. Moreover, due to the differential effect on melting as a function of the location of the dA→7 mutations, these results are consistent with the N1–N7 base pairing of the transWH pairs. Together, the simulation and base substitution experiments predict that the coralyne-induced homo-(dA) duplex structure adopts the transWH geometry.     

Adenine Methylation in Eukaryotic DNA

N⁶-Methyladenine (m⁶A) has been found in DNAs of various eukaryotes (algae, fungi, protozoa, and higher plants). Like bacterial DNA, DNAs of these organisms are subject to enzymatic modification (methylation) not only at cytosine, but also at adenine bases. There is indirect evidence that adenine methylation of the genome occurs in animals as well. In plants, m⁶A was detected in total, mitochondrial, and nuclear DNAs. It was observed that both adenines and cytosines can be methylated in one gene (DRM2). Open reading frames coding for homologs of bacterial adenine DNA methyltransferases were revealed in protozoan, yeast, higher plant, insect, nematode, and vertebrate genomes, suggesting the presence of adenine DNA methyltransferases in evolutionarily distant eukaryotes. The first higher-eukaryotic adenine DNA N⁶-methyltransferase (wad-mtase) was isolated from vacuolar vesicles of wheat coleoptiles. The enzyme depends on Mg²⁺ or Ca²⁺ and, in the presence of S-adenosyl-L-methionine, methylates de novo the first adenine of the sequence TGATCA in single- and double-stranded DNAs, preferring the former. Adenine methylation of eukaryotic DNA is probably involved in regulating gene expression and replication, including that of mitochondrial DNA; plays a role in controlling the persistence of foreign DNA in the cell; and acts as a component of a plant restriction— modification system. Thus, the eukaryotic cell has at least two different systems for enzymatic methylation of DNA (at adenines and at cytosines) and a special mechanism regulating the functions of genes via a combinatorial hierarchy of these interdependent modifications of the genome.__________Translated from Molekulyarnaya Biologiya, Vol. 39, No. 4, 2005, pp. 557–566.Original Russian Text Copyright © 2005 by Vanyushin.To the memory of my teacher, Academician Andrei Nikolaevich Belozersky     

Methylated nucleotide content of mitochondrial ribosomal RNA from hamster cells

We have examined the state of methylation of mitochondrial ribosomal RNA from cultured hamster (BHK-21) cells. Ethidium-sensitive, and hence mitochondrion-specific, methylation levels were determined using multiple isotope techniques and improved purification procedures. The larger mitochondrial rRNA species, 17 S RNA, was found to contain 0.13 methyl group per 100 nucleotides and the smaller, 13 S RNA, 0.37. Methylated nucleotide and base analysis indicated that 17 S RNA contained one ribose-methylated residue (UmUp) per molecule and one unidentified residue; and 13 S RNA contained one methylated cytosine residue, one N⁶-dimethyladenine residue and one thymine residue per molecule. Possible evolutionary implications of these findings have been discussed.