The human thymine-DNA glycosylase (TDG) initiates the base excision repair (BER) pathway to remove spontaneous and induced DNA base damage. It was first biochemically characterized for its ability to remove T mispaired with G in CpG context. TDG is involved in the epigenetic regulation of gene expressions by protecting CpG-rich promoters from de novo DNA methylation. Here we demonstrate that TDG initiates aberrant repair by excising T when it is paired with a damaged adenine residue in DNA duplex. TDG targets the non-damaged DNA strand and efficiently excises T opposite of hypoxanthine (Hx), 1,N(6)-ethenoadenine, 7,8-dihydro-8-oxoadenine and abasic site in TpG/CpX context, where X is a modified residue. In vitro reconstitution of BER with duplex DNA containing Hx•T pair and TDG results in incorporation of cytosine across Hx. Furthermore, analysis of the mutation spectra inferred from single nucleotide polymorphisms in human population revealed a highly biased mutation pattern within CpG islands (CGIs), with enhanced mutation rate at CpA and TpG sites. These findings demonstrate that under experimental conditions used TDG catalyzes sequence context-dependent aberrant removal of thymine, which results in TpG, CpA→CpG mutations, thus providing a plausible mechanism for the putative evolutionary origin of the CGIs in mammalian genomes.

Molecular barcoding technologies that uniquely identify single cells are hampered by limitations in barcode measurement. Readout by sequencing does not preserve the spatial organization of cells in tissues, whereas imaging methods preserve spatial structure but are less sensitive to barcode sequence. Here we introduce a system for image-based readout of short (20-base-pair) DNA barcodes. In this system, called Zombie, phage RNA polymerases transcribe engineered barcodes in fixed cells. The resulting RNA is subsequently detected by fluorescent in situ hybridization. Using competing match and mismatch probes, Zombie can accurately discriminate single-nucleotide differences in the barcodes. This method allows in situ readout of dense combinatorial barcode libraries and single-base mutations produced by CRISPR base editors without requiring barcode expression in live cells. Zombie functions across diverse contexts, including cell culture, chick embryos and adult mouse brain tissue. The ability to sensitively read out compact and diverse DNA barcodes by imaging will facilitate a broad range of barcoding and genomic recording strategies.

Base-pairing interactions mediate many intermolecular target recognition events. Even a single base-pair mismatch can cause a substantial difference in activity but how such changes influence the target search kinetics in vivo is unknown. Here, we use high-throughput sequencing and quantitative super-resolution imaging to probe the mutants of bacterial small RNA, SgrS, and their regulation of ptsG mRNA target. Mutations that disrupt binding of a chaperone protein, Hfq, and are distal to the mRNA annealing region still decrease the rate of target association, kon, and increase the dissociation rate, koff, showing that Hfq directly facilitates sRNA-mRNA annealing in vivo. Single base-pair mismatches in the annealing region reduce kon by 24-31% and increase koff by 14-25%, extending the time it takes to find and destroy the target by about a third. The effects of disrupting contiguous base-pairing are much more modest than that expected from thermodynamics, suggesting that Hfq buffers base-pair disruptions.

Lambda exonuclease (λ exo) plays an important role in the resection of DNA ends for DNA repair. Currently, it is also a widely used enzymatic tool in genetic engineering, DNA-binding protein mapping, nanopore sequencing and biosensing. Herein, we disclose two noncanonical properties of this enzyme and suggest a previously undescribed hydrophobic interaction model between λ exo and DNA substrates. We demonstrate that the length of the free portion of the substrate strand in the dsDNA plays an essential role in the initiation of digestion reactions by λ exo. A dsDNA with a 5' non-phosphorylated, two-nucleotide-protruding end can be digested by λ exo with very high efficiency. Moreover, we show that when a conjugated structure is covalently attached to an internal base of the dsDNA, the presence of a single mismatched base pair at the 5' side of the modified base may significantly accelerate the process of digestion by λ exo. A detailed comparison study revealed additional π-π stacking interactions between the attached label and the amino acid residues of the enzyme. These new findings not only broaden our knowledge of the enzyme but will also be very useful for research on DNA repair and in vitro processing of nucleic acids.

DNA damage response (DDR) genes have low mutation rates, which may restrict their clinical applications in predicting the outcomes of immune checkpoint inhibitor (ICI) treatment. Thus, a systemic analysis of multiple DDR genes is needed to identify potential biomarkers of ICI efficacy.

The uracil DNA glycosylase superfamily consists of several distinct families. Family 2 mismatch-specific uracil DNA glycosylase (MUG) from Escherichia coli is known to exhibit glycosylase activity on three mismatched base pairs, T/U, G/U and C/U. Family 1 uracil N-glycosylase (UNG) from E. coli is an extremely efficient enzyme that can remove uracil from any uracil-containing base pairs including the A/U base pair. Here, we report the identification of an important structural determinant that underlies the functional difference between MUG and UNG. Substitution of a Lys residue at position 68 with Asn in MUG not only accelerates the removal of uracil from mismatched base pairs but also enables the enzyme to gain catalytic activity on A/U base pairs. Binding and kinetic analysis demonstrate that the MUG-K68N substitution results in enhanced ground state binding and transition state interactions. Molecular modeling reveals that MUG-K68N, UNG-N123 and family 5 Thermus thermophiles UDGb-A111N can form bidentate hydrogen bonds with the N3 and O4 moieties of the uracil base. Genetic analysis indicates the gain of function for A/U base pairs allows the MUG-K68N mutant to remove uracil incorporated into the genome during DNA replication. The implications of this study in the origin of life are discussed.

We report here the genome-wide changes resulting from low N (N-W+), low water (N+W-)) and dual stresses (N-W-) in root and shoot tissues of two rice genotypes, namely, IR 64 (IR64) and Nagina 22 (N22), and their association with the QTLs for nitrogen use efficiency. For all the root parameters, except for root length under N-W+, N22 performed better than IR64. Chlorophyll a, b and carotenoid content were higher in IR64 under N+W+ treatment and N-W+ and N+W- stresses; however, under dual stress, N22 had higher chlorophyll b content. While nitrite reductase, glutamate synthase (GS) and citrate synthase assays showed better specific activity in IR64, glutamate dehydrogenase showed better specific activity in N22 under dual stress (N-W-); the other N and C assimilating enzymes showed similar but low specific activities in both the genotypes. A total of 8926 differentially expressed genes (DEGs) were identified compared to optimal (N+W+) condition from across all treatments. While 1174, 698 and 903 DEGs in IR64 roots and 1197, 187 and 781 in N22 roots were identified, nearly double the number of DEGs were found in the shoot tissues; 3357, 1006 and 4005 in IR64 and 4004, 990 and 2143 in N22, under N-W+, N+W- and N-W- treatments, respectively. IR64 and N22 showed differential expression in 15 and 11 N-transporter genes respectively, under one or more stress treatments, out of which four showed differential expression also in N+W- condition. The negative regulators of N- stress, e.g., NIGT1, OsACTPK1 and OsBT were downregulated in IR64 while in N22, OsBT was not downregulated. Overall, N22 performed better under dual stress conditions owing to its better root architecture, chlorophyll and porphyrin synthesis and oxidative stress management. We identified 12 QTLs for seed and straw N content using 253 recombinant inbred lines derived from IR64 and N22 and a 5K SNP array. The QTL hotspot region on chromosome 6 comprised of 61 genes, of which, five were DEGs encoding for UDP-glucuronosyltransferase, serine threonine kinase, anthocyanidin 3-O-glucosyltransferase, and nitrate induced proteins. The DEGs, QTLs and candidate genes reported in this study can serve as a major resource for both rice improvement and functional biology.

NH groups in proteins or nucleic acids are the most challenging target for chemical shift prediction. Here we show that the RNA base pair triplet motif dictates imino chemical shifts in its central base pair. A lookup table is established that links each type of base pair triplet to experimental chemical shifts of the central base pair, and can be used to predict imino chemical shifts of RNAs to remarkable accuracy. Strikingly, the semiempirical method can well interpret the variations of chemical shifts for different base pair triplets, and is even applicable to non-canonical motifs. This finding opens an avenue for predicting chemical shifts of more complicated RNA motifs. Furthermore, we combine the imino chemical shift prediction with NMR relaxation dispersion experiments targeting both 15N and 1HN of the imino group, and verify a previously characterized excited state of P5abc subdomain including an earlier speculated non-native G•G mismatch.

Oxidative damage to 5-methylcytosine in DNA, followed by deamination, yields thymine glycol (Tg), 5,6-dihydroxy-5,6-dihydrothymine, mispaired with deoxyguanosine. The structure of the 5R Tg.G mismatch pair has been refined using a combination of simulated annealing and isothermal molecular dynamics calculations restrained by NMR-derived distance restraints and torsion angle restraints in 5'-d(G(1)T(2)G(3)C(4)G(5)Tg(6)G(7)T(8)T(9)T(10)G(11)T(12))-3'.5'-d(A(13)C(14)A(15)A(16)A(17)C(18)G(19)C(20)G(21)C(22)A(23)C(24))-3'; Tg = 5R Tg. In this duplex the cis-5R,6S:trans-5R,6R equilibrium favors the cis-5R,6S epimer [Brown, K. L., Adams, T., Jasti, V. P., Basu, A. K., and Stone, M. P. (2008) J. Am. Chem. Soc. 130, 11701-11710]. The cis-5R,6S Tg lesion is in the wobble orientation such that Tg(6) O(2) is proximate to G(19) N1H and Tg(6) N3H is proximate to G(19) O(6). Both Tg(6) and the mismatched nucleotide G(19) remain stacked in the helix. The Tg(6) nucleotide shifts toward the major groove and stacks below the 5'-neighbor base G(5), while its complement G(19) stacks below the 5'-neighbor C(20). In the 3'-direction, stacking between Tg(6) and the G(7).C(18) base pair is disrupted. The solvent-accessible surface area of the Tg nucleotide increases as compared to the native Watson-Crick hydrogen-bonded T.A base pair. An increase in T(2) relaxation rates for the Tg(6) base protons is attributed to puckering of the Tg base, accompanied by increased disorder at the Tg.G mismatch pair. The axial vs equatorial conformation of the Tg(6) CH(3) group cannot be determined with certainty from the NMR data. The rMD trajectories suggest that in either the axial or equatorial conformations the cis-5R,6S Tg lesion does not form strong intrastrand hydrogen bonds with the imidazole N7 atom of the 3'-neighbor purine G(7). The wobble pairing and disorder of the Tg.G mismatch correlate with the reduced thermodynamic stability of the mismatch and likely modulate its recognition by DNA base excision repair systems.

Phytochromes are the main plant photoreceptors regulating multiple developmental processes. However, the regulatory network of phytochrome-mediated plant reproduction has remained largely unexplored. There are three phytochromes in rice, phyA, phyB and phyC. No changes in fertility are observed in the single mutants, whereas the seed-setting rate of the phyA phyB double mutant is significantly reduced. Histological and cytological analyses showed that the reduced fertility of the phyA phyB mutant was due to defects in both anther and pollen development. The four anther lobes in the phyA phyB mutant were developed at different stages with fewer pollen grains, most of which were aborted. At the mature stage, more than one lobe in the double mutant was just consisted of several cell layers. To identify genes involved in phytochrome-mediated anther development, anther transcriptomes of phyA, phyB and phyA phyB mutants were compared to that of wild-type rice respectively. Analysis of 2,241 double-mutant-specific differentially expressed transcripts revealed that the metabolic profiles, especially carbohydrate metabolism, were altered greatly, and heat-shock responses were activated in the double mutant. This study firstly provides valuable insight into the complex regulatory networks underlying phytochrome-mediated anther and pollen development in plants, and offers novel clues for hybrid rice breeding.

DNA polymerase III mis-insertion may, where not corrected by its 3'→ 5' exonuclease or the mismatch repair (MMR) function, result in all possible non-cognate base pairs in DNA generating base substitutions. The most thermodynamically unstable base pair, the cytosine (C)⋅C mismatch, destabilizes adjacent base pairs, is resistant to correction by MMR in Escherichia coli, and its repair mechanism remains elusive. We present here in vitro evidence that C⋅C mismatch can be processed by base excision repair initiated by the E. coli formamidopyrimidine-DNA glycosylase (Fpg) protein. The k cat for C⋅C is, however, 2.5 to 10 times lower than for its primary substrate 8-oxoguanine (oxo8G)⋅C, but approaches those for 5,6-dihydrothymine (dHT)⋅C and thymine glycol (Tg)⋅C. The K M values are all in the same range, which indicates efficient recognition of C⋅C mismatches in DNA. Fpg activity was also exhibited for the thymine (T)⋅T mismatch and for N 4- and/or 5-methylated C opposite C or T, Fpg activity being enabled on a broad spectrum of DNA lesions and mismatches by the flexibility of the active site loop. We hypothesize that Fpg plays a role in resolving C⋅C in particular, but also other pyrimidine⋅pyrimidine mismatches, which increases survival at the cost of some mutagenesis.

Base pair mismatches in genomic DNA can result in mutagenesis, and consequently in tumorigenesis. To investigate how mismatch repair deficiency increases mutagenicity under oxidative stress, we examined the type and frequency of mutations arising in the mucosa of the small intestine of mice carrying a reporter gene encoding guanine phosphoribosyltransferase (gpt) and in which the Msh2 gene, which encodes a component of the mismatch repair system, was either intact (Msh2+/+::gpt/0; Msh2-bearing) or homozygously knockout (KO) (Msh2-/-::gpt/0; Msh2-KO).

Smad-interacting protein-1 (SIP1) has been implicated in the development of Mowat-Wilson syndrome whose patients exhibit Hirschsprung disease, an aganglionosis of the large intestine, as well as other phenotypes. We have identified and cloned two sip1 orthologues in zebrafish. Both sip1 orthologues are expressed maternally and have dynamic zygotic expression patterns that are temporally and spatially distinct. We have investigated the function of both orthologues using translation and splice-blocking morpholino antisense oligonucleotides. Knockdown of the orthologues causes axial and neural patterning defects consistent with the previously described function of SIP1 as an inhibitor of BMP signaling. In addition, knockdown of both genes leads to a significant reduction/loss of the post-otic cranial neural crest. This results in a subsequent absence of neural crest precursors in the posterior pharyngeal arches and a loss of enteric precursors in the intestine.

Artificially Expanded Genetic Information Systems (AEGIS) add independently replicable unnatural nucleotide pairs to the natural G:C and A:T/U pairs found in native DNA, joining the unnatural pairs through alternative modes of hydrogen bonding. Whether and how AEGIS pairs are recognized and processed by multi-subunit cellular RNA polymerases (RNAPs) remains unknown. Here, we show that E. coli RNAP selectively recognizes unnatural nucleobases in a six-letter expanded genetic system. High-resolution cryo-EM structures of three RNAP elongation complexes containing template-substrate UBPs reveal the shared principles behind the recognition of AEGIS and natural base pairs. In these structures, RNAPs are captured in an active state, poised to perform the chemistry step. At this point, the unnatural base pair adopts a Watson-Crick geometry, and the trigger loop is folded into an active conformation, indicating that the mechanistic principles underlying recognition and incorporation of natural base pairs also apply to AEGIS unnatural base pairs. These data validate the design philosophy of AEGIS unnatural basepairs. Further, we provide structural evidence supporting a long-standing hypothesis that pair mismatch during transcription occurs via tautomerization. Together, our work highlights the importance of Watson-Crick complementarity underlying the design principles of AEGIS base pair recognition.

Genomic stability depends in part on an efficient DNA lesion recognition and correction by the DNA mismatch repair (MMR) system. We investigated mutations arising spontaneously in rice OsMsh6 mutants by specific-locus amplified fragment sequencing. Totally 994 single-nucleotide mutations were identified in three mutants and on average the mutation density is about 1/136.72 Kb per mutant line. These mutations were relatively randomly distributed in genome and might be accumulated in generation-dependent manner. All possible base transitions and base transversions could be seen and the ratio of transitions to transversions was about 3.12. We also observed the nearest-neighbor bias around the mutated base. Our data suggests that OsMsh6 (LOC_Os09g24220) is important in ensuring genome stability by recognizing mismatches that arise spontaneously and provides useful information for investigating the function of the OsMsh6 gene in DNA repair and exploiting MMR mutants in rice induced mutation breeding.

Cytosine to thymine (C>T) somatic mutation is highly enriched in certain types of cancer, and most commonly occurs via deamination of a 5-methylcytosine (5mC) to thymine, in the context of a CpG dinucleotide. In theory, deamination should occur at equal rates to both 5mC nucleotides on opposite strands. In most cases, the resulting T:G or G:T mismatch can be repaired by thymine DNA glycosylase activities. However, while some hotspot-associated CpG mutations have approximately equal numbers of mutations that resulted either from C>T or G>A in a CpG dinucleotide, many showed strand bias, being skewed toward C>T of the first base pair or G>A of the second base pair. Using the IDH2 Arg140 codon as a case study, we show that the two possible T:G mismatches at the codon-specific CpG site have differing effects on transcription factor ETS1 binding affinity, differentially affecting access of a repair enzyme (MBD4) to the deamination-caused T:G mismatch. Our study thus provides a plausible mechanism for exclusion of repair enzymes by the differential binding of transcription factors affecting the rate at which the antecedent opposite-strand mutations occur.

We previously developed reporter-peptide nucleic acid (PNA)-peptides for sequence-specific radioimaging and fluorescence imaging of particular mRNAs in cells and tumors. However, a direct test for PNA-peptide hybridization with RNA in the cytoplasm would be desirable. Thiazole orange (TO) dye at the 5' end of a hybridization agent shows a strong increase in fluorescence quantum yield when stacked upon a 5' terminal base pair, in solution and in cells. We hypothesized that hybridization agents with an internal TO could distinguish a single base mutation in RNA. Thus, we designed KRAS2 PNA-IGF1 tetrapeptide agents with an internal TO adjacent to the middle base of the 12th codon, a frequent site of cancer-initiating mutations. Our molecular dynamics calculations predicted a disordered bulge with weaker hybridization resulting from a single RNA mismatch. We observed that single-stranded PNA-IGF1 tetrapeptide agents with an internal TO showed low fluorescence, but fluorescence escalated 5-6-fold upon hybridization with KRAS2 RNA. Circular dichroism melting curves showed ∼10 °C higher Tm for fully complementary vs single base mismatch TO-PNA-peptide agent duplexes with KRAS2 RNA. Fluorescence measurements of treated human lung cancer cells similarly showed elevated cytoplasmic fluorescence intensity with fully complementary vs single base mismatch agents. Sequence-specific elevation of internal TO fluorescence is consistent with our hypothesis of detecting cytoplasmic PNA-peptide:RNA hybridization if a mutant agent encounters the corresponding mutant mRNA.

Specific interactions between T7 RNA polymerase and its promoter have been studied by a simple steady-state kinetic assay using synthetic oligonucleotide promoters that produce a short five-base message. A series of promoters with upstream lengths extending to promoter positions -19, -17, -14, and -12 show that promoters extending to -19 and -17 produce very specific transcripts with initiation rate constant Kcat = 50 min-1 and a Michaelis constant Km = 0.02 microM, indicating that the consensus sequence to position -17 is sufficient for maximum promoter usage. Shortening the upstream region of the promoter to -14 substantially increases Km (0.3 microM) but does not significantly reduce the maximum velocity (kcat = 30 min-1). Finally, truncation of the promoter at position -12 results in extremely low levels of specific transcription. The coding and noncoding strands appear to make different contributions to promoter recognition. Although the double-stranded promoter of upstream length -12 is very poor as a transcription template, extension of only the noncoding strand to -17 very significantly improves both Kcat and Km. In contrast, extension of only the coding strand results in no significant improvement. Substitution of an AT base pair at position -10 by CG (as found in T3 RNA polymerase promoters) produces a 10-fold increase in Km, with little effect on Kcat. Comparison of two promoters containing a base pair mismatch at this site (AG or CT) demonstrates that promoter recognition is very sensitive to the nature of the base on the noncoding strand and is only slightly affected by the presence of a mismatch created by a wrong base in the coding strands.(ABSTRACT TRUNCATED AT 250 WORDS)

RNA editing of the hepatitis delta virus (HDV) is essential for generating the large delta antigen, which is crucial for virion assembly. In HDV genotype 1 (HDV-1), editing occurs within the context of the unbranched rod-like structure characteristic of HDV RNA, while RNA editing in HDV-3 requires a branched double-hairpin structure. The regulation of RNA editing in HDV-2 and HDV-4 remains uncertain. Based on predictions of the unbranched rod-like RNA structures of HDV-2 and HDV-4, the editing site occurs as an A.C mismatch pair, surrounded by four base pairs upstream and two base pairs downstream of the editing site, respectively. To investigate HDV-2 and HDV-4 RNA editing, cultured cells were transfected with non-replicating editing reporters carrying wild-type sequences or specific mutations. The results revealed that the editing rates observed for wild-type HDV-2 and HDV-4 were fairly similar, albeit lower than that of HDV-1. Like HDV-1, both HDV-2 and HDV-4 showed a reduction in editing rate when the A.C mismatch pair and the immediately upstream base-paired region were disturbed. Notably, extending the downstream base-paired region from two to three or four (forming a structure identical to that of HDV-1) base pairs increased editing rate. Furthermore, we presented novel evidence that indicates the importance of the first bulge's size, located upstream of the editing site, and the base-pairing length within 7-13 and 28-39 nucleotides downstream of the editing site in influencing the HDV-4 editing rate. To summarize, our analyses suggest that the unbranched rod-like structures surrounding the editing site of HDV-2 and HDV-4 play a crucial role in regulating their RNA editing rates.

The specificity of small interfering RNA (siRNA)-mediated gene silencing is a critical consideration for the application of RNA interference (RNAi). While the discovery of potential off-target effects by siRNAs is of concern, no systematic analysis has been conducted to explore the specificity of RNAi. Here, we present a study where a functionally validated siRNA (siCD46) was examined for silencing specificity on all possible 57 permutated target sites, each carrying a single-nucleotide mutation that would generate a mismatch when paired with siRNA antisense strand. We found that it was not only the position of the mismatched base pair, but also the identity of the nucleotides forming the mismatch that influenced silencing. Surprisingly, mismatches formed between adenine (A) and cytosine (C), in addition to the G:U wobble base pair, were well tolerated and target sites containing such mismatches were silenced almost as efficiently as its fully matched counterpart by siCD46. Northern blots showed that the silencing of fusion genes harboring the mutated target sites involved target mRNA degradation. This study provides direct evidence that the target recognition of siRNA is far more degenerative than previously considered. This finding is instrumental in the understanding of RNAi specificity and may aid the computational prediction of RNA secondary structure.