Targeted Hybridization Capture of SARS-CoV-2 and Metagenomics Enables Genetic Variant Discovery and Nasal Microbiome Insights.

The emergence of novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genetic variants that may alter viral fitness highlights the urgency of widespread next-generation sequencing (NGS) surveillance. To profile genetic variants of the entire SARS-CoV-2 genome, we developed and clinically validated a hybridization capture SARS-CoV-2 NGS assay, integrating novel methods for panel design using double-stranded DNA (dsDNA) biotin-labeled probes, and built accompanying software. This test is the first hybrid capture-based NGS assay given Food and Drug Administration (FDA) emergency use authorization for detection of the SARS-CoV-2 virus. The positive and negative percent agreement (PPA and NPA, respectively) were defined in comparison to the results for an orthogonal real-time reverse transcription polymerase chain reaction (RT-PCR) assay (PPA and NPA, 96.7 and 100%, respectively). The limit of detection was established to be 800 copies/ml with an average fold enrichment of 46,791. Furthermore, utilizing the research-use-only analysis to profile the variants, we identified 55 novel mutations, including 11 in the functionally important spike protein. Finally, we profiled the full nasopharyngeal microbiome using metagenomics and found overrepresentation of 7 taxa and evidence of macrolide resistance in SARS-CoV-2-positive patients. This hybrid capture NGS assay, coupled with optimized software, is a powerful approach to detect and comprehensively map SARS-CoV-2 genetic variants for tracking viral evolution and guiding vaccine updates. IMPORTANCE This is the first FDA emergency-use-authorized hybridization capture-based next-generation sequencing (NGS) assay to detect the SARS-CoV-2 genome. Viral metagenomics and the novel hybrid capture NGS-based assay, along with its research-use-only analysis, can provide important genetic insights into SARS-CoV-2 and other emerging pathogens and improve surveillance and early detection, potentially preventing or mitigating new outbreaks. Better understanding of the continuously evolving SARS-CoV-2 viral genome and the impact of genetic variants may provide individual risk stratification, precision therapeutic options, improved molecular diagnostics, and population-based therapeutic solutions.

Repair of thymine glycol by hNth1 and hNeil1 is modulated by base pairing and cis-trans epimerization.

Oxidation of thymine yields 5,6-dihydroxy-5,6-dihydrothymine (thymine glycol. Tg) which, as cis 5S,6R and 5R,6S 2'-deoxyribonucleoside diastereoisomers (dTg1, dTg2), are in equilibrium with their trans 5S,6S and 5R,6R epimers. The stereoselective excision of Tg from DNA by the mammalian orthologs of E. coli DNA N-glycosylase/AP lyases Nth and Nei was reported using substrates in which Tg opposed adenine. Since we showed that Tg is the major product of oxidation of 5-methylcytosine, we asked if the opposing purine influenced stereospecific enzymatic excision. The human ortholog hNth1 released Tg2 much more rapidly than Tg1 regardless of the opposing purine. In contrast, hNeil1 released Tg non-stereoselectively, but the rate of excision was much greater when Tg opposed guanine. Remarkably, the kinetics of excision of Tg by hNth1 and hNeil1 were biphasic, describing a double exponential curve which yielded two rate constants. We suggest that the greater rate constant describes the rate of enzymatic excision of Tg. The smaller rate constant represents the equilibrium constant for the cis and trans epimerization of dTg1 and dTg2 in high molecular weight DNA. Thus, only one of the epimers of dTg1 and dTg2 are enzymatically processed but it is not yet known whether it is cis or trans. Thus, base excision repair of Tg in mammals is mediated by at least two DNA N-glycosylase/AP lyases which are affected by the nature of the diastereoisomer of dTg, the rate of cis-trans epimerization of each diastereoisomer, and the nature of the opposing purine.

Targeted deletion of mNth1 reveals a novel DNA repair enzyme activity.

DNA N-glycosylase/AP (apurinic/apyrimidinic) lyase enzymes of the endonuclease III family (nth in Escherichia coli and Nth1 in mammalian organisms) initiate DNA base excision repair of oxidized ring saturated pyrimidine residues. We generated a null mouse (mNth1(-/-)) by gene targeting. After almost 2 years, such mice exhibited no overt abnormalities. Tissues of mNth1(-/-) mice contained an enzymatic activity which cleaved DNA at sites of oxidized thymine residues (thymine glycol [Tg]). The activity was greater when Tg was paired with G than with A. This is in contrast to Nth1, which is more active against Tg:A pairs than Tg:G pairs. We suggest that there is a back-up mammalian repair activity which attacks Tg:G pairs with much greater efficiency than Tg:A pairs. The significance of this activity may relate to repair of oxidized 5-methyl cytosine residues (5meCyt). It was shown previously (S. Zuo, R. J. Boorstein, and G. W. Teebor, Nucleic Acids Res. 23:3239-3243, 1995) that both ionizing radiation and chemical oxidation yielded Tg from 5meCyt residues in DNA. Thus, this previously undescribed, and hence novel, back-up enzyme activity may function to repair oxidized 5meCyt residues in DNA while also being sufficient to compensate for the loss of Nth1 in the mutant mice, thereby explaining the noninformative phenotype.

Expression of the human DNA glycosylase hSMUG1 in Trypanosoma brucei causes DNA damage and interferes with J biosynthesis.

In kinetoplastid flagellates such as Trypanosoma brucei, a small percentage of the thymine residues in the nuclear DNA is replaced by the modified base beta-D-glucosyl-hydroxymethyluracil (J), mostly in repetitive sequences like the telomeric GGGTTA repeats. In addition, traces of 5-hydroxymethyluracil (HOMeUra) are present. Previous work has suggested that J is synthesised in two steps via HOMedU as an intermediate, but as J synthesising enzymes have not yet been identified, the biosynthetic pathway remains unclear. To test a model in which HOMeUra functions as a precursor of J, we introduced an inducible gene for the human DNA glycosylase hSMUG1 into bloodstream form T.brucei. In higher eukaryotes SMUG1 excises HOMeUra as part of the base excision repair system. We show that expression of the gene in T.brucei leads to massive DNA damage in J-modified sequences and results in cell cycle arrest and, eventually, death. hSMUG1 also reduces the J content of the trypanosome DNA. This work supports the idea that HOMeUra is a precursor of J, freely accessible to a DNA glycosylase.

Substrate specificity of human endonuclease III (hNTH1). Effect of human APE1 on hNTH1 activity.

Base excision repair of oxidized pyrimidines in human DNA is initiated by the DNA N-glycosylase/apurinic/apyrimidinic (AP) lyase, human NTH1 (hNTH1), the homolog of Escherichia coli endonuclease III (Nth). In contrast to Nth, the DNA N-glycosylase activity of hNTH1 is 7-fold greater than its AP lyase activity when the DNA substrate contains a thymine glycol (Tg) opposite adenine (Tg:A) (Marenstein, D. R., Ocampo, M. T. A., Chan, M. K., Altamirano, A., Basu, A. K., Boorstein, R. J., Cunningham, R. P., and Teebor, G. W. (2001) J. Biol. Chem. 276, 21242-21249). When Tg is opposite guanine (Tg:G), the two activities are of the same specific activity as the AP lyase activity of hNTH1 against Tg:A (Ocampo, M. T. A., Chaung, W., Marenstein, D. R., Chan, M. K., Altamirano, A., Basu, A. K., Boorstein, R. J., Cunningham, R. P., and Teebor, G. W. (2002) Mol. Cell. Biol. 22, 6111-6121). We demonstrate here that hNTH1 was inhibited by the product of its DNA N-glycosylase activity directed against Tg:G, the AP:G site. In contrast, hNTH1 was not as inhibited by the AP:A site arising from release of Tg from Tg:A. Addition of human APE1 (AP endonuclease-1) increased dissociation of hNTH1 from the DNA N-glycosylase-generated AP:A site, resulting in abrogation of AP lyase activity and an increase in turnover of the DNA N-glycosylase activity of hNTH1. Addition of APE1 did not abrogate hNTH1 AP lyase activity against Tg:G. The stimulatory protein YB-1 (Marenstein et al.), added to APE1, resulted in an additive increase in both activities of hNTH1 regardless of base pairing. Tg:A is formed by oxidative attack on thymine opposite adenine. Tg:G is formed by oxidative attack on 5-methylcytosine opposite guanine (Zuo, S., Boorstein, R. J., and Teebor, G. W. (1995) Nucleic Acids Res. 23, 3239-3243). It is possible that the in vitro substrate selectivity of mammalian NTH1 and the concomitant selective stimulation of activity by APE1 are indicative of selective repair of oxidative damage in different regions of the genome.