TERRA ONTseq: a long-read-based sequencing pipeline to study the human telomeric transcriptome.

The long noncoding RNA TERRA is transcribed from telomeres in virtually all eukaryotes with linear chromosomes. In humans, TERRA transcription is driven in part by promoters comprising CpG dinucleotide-rich repeats of 29 bp repeats, believed to be present in half of the subtelomeres. Thus far, TERRA expression has been analyzed mainly using molecular biology-based approaches that only generate partial and somehow biased results. Here, we present a novel experimental pipeline to study human TERRA based on long-read sequencing (TERRA ONTseq). By applying TERRA ONTseq to different cell lines, we show that the vast majority of human telomeres produce TERRA and that the cellular levels of TERRA transcripts vary according to their chromosomes of origin. Using TERRA ONTseq, we also identified regions containing TERRA transcription start sites (TSSs) in more than half of human subtelomeres. TERRA TSS regions are generally found immediately downstream from 29 bp repeat-related sequences, which appear to be more widespread than previously estimated. Finally, we isolated a novel TERRA promoter from the highly expressed subtelomere of the long arm of Chromosome 7. With the development of TERRA ONTseq, we provide a refined picture of human TERRA biogenesis and expression and we equip the scientific community with an invaluable tool for future studies.

Dinucleotide biases in RNA viruses that infect vertebrates or invertebrates.

IMPORTANCE: Akin to a molecular signature, dinucleotide composition can be exploited by the zinc-finger antiviral protein (ZAP) to restrict CpG-rich (and UpA-rich) RNA viruses. ZAP evolved in tetrapods, and it is not encoded by invertebrates and fish. Because a systematic analysis is missing, we analyzed the genomes of RNA viruses that infect vertebrates or invertebrates. We show that vertebrate single-stranded (ss) RNA(+) viruses and, to a lesser extent, double-stranded RNA viruses tend to have stronger CpG bias than invertebrate viruses. Conversely, ssRNA(-) viruses have similar dinucleotide composition whether they infect vertebrates or invertebrates. Analysis of ssRNA(+) viruses that infect mammals, reptiles, and fish indicated that ZAP is unlikely to be a major driver of CpG depletion. We also show that, compared to other coronaviruses, the genome of SARS-CoV-2 is not homogeneously CpG-depleted. Our study provides new insights into virus evolution and strategies for recoding RNA virus genomes.

Strong Selectional Forces Fine-Tune CpG Content in Genes Involved in Neurological Disorders as Revealed by Codon Usage Patterns.

Neurodegenerative disorders cause irreversible damage to the neurons and adversely affect the quality of life. Protein misfolding and their aggregation in specific parts of the brain, mitochondrial dysfunction, calcium load, proteolytic stress, and oxidative stress are among the causes of neurodegenerative disorders. In addition, altered metabolism has been associated with neurodegeneration as evidenced by reductions in glutamine and alanine in transient global amnesia patients, higher homocysteine-cysteine disulfide, and lower methionine decline in serum urea have been observed in Alzheimer's disease patients. Neurodegeneration thus appears to be a culmination of altered metabolism. The study's objective is to analyze various attributes like composition, physical properties of the protein, and factors like selectional and mutational forces, influencing codon usage preferences in a panel of genes involved directly or indirectly in metabolism and contributing to neurodegeneration. Various parameters, including gene composition, dinucleotide analysis, Relative synonymous codon usage (RSCU), Codon adaptation index (CAI), neutrality and parity plots, and different protein indices, were computed and analyzed to determine the codon usage pattern and factors affecting it. The correlation of intrinsic protein properties such as the grand average of hydropathicity index (GRAVY), isoelectric point, hydrophobicity, and acidic, basic, and neutral amino acid content has been found to influence codon usage. In genes up to 800 amino acids long, the GC3 content was highly variable, while GC12 content was relatively constant. An optimum CpG content is present in genes to maintain a high expression level as required for genes involved in metabolism. Also observed was a low codon usage bias with a higher protein expression level. Compositional parameters and nucleotides at the second position of codons played essential roles in explaining the extent of bias. Overall analysis indicated that the dominance of selection pressure and compositional constraints and mutational forces shape codon usage.

Structural characterization of dicyanopyridine containing DNMT1-selective, non-nucleoside inhibitors.

DNMT1 maintains the parental DNA methylation pattern on newly replicated hemimethylated DNA. The failure of this maintenance process causes aberrant DNA methylation that affects transcription and contributes to the development and progression of cancers such as acute myeloid leukemia. Here, we structurally characterized a set of newly discovered DNMT1-selective, reversible, non-nucleoside inhibitors that bear a core 3,5-dicyanopyridine moiety, as exemplified by GSK3735967, to better understand their mechanism of inhibition. All of the dicyanopydridine-containing inhibitors examined intercalate into the hemimethylated DNA between two CpG base pairs through the DNA minor groove, resulting in conformational movement of the DNMT1 active-site loop. In addition, GSK3735967 introduces two new binding sites, where it interacts with and stabilizes the displaced DNMT1 active-site loop and it occupies an open aromatic cage in which trimethylated histone H4 lysine 20 is expected to bind. Our work represents a substantial step in generating potent, selective, and non-nucleoside inhibitors of DNMT1.

Differential ETS1 binding to T:G mismatches within a CpG dinucleotide contributes to C-to-T somatic mutation rate of the IDH2 hotspot at codon Arg140.

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.

Msuite2: All-in-one DNA methylation data analysis toolkit with enhanced usability and performance.

DNA methylation is an important epigenetic regulator that plays crucial roles in various biological processes. Recent developments in experimental approaches and dramatic expansion of sequencing capacities have imposed new challenges in the analysis of large-scale, cross-species DNA methylation data. Hence, user-friendly toolkits with high usability and performance are in urgent need. In this work, we present Msuite2, an easy-to-use, all-in-one, and universal toolkit for DNA methylation data analysis and visualization with high flexibility, usability, and performance. Msuite2 is among the fastest tools in read alignment (in particular, it runs as much as 5x faster than its predecessor, Msuite1) with low computing resource usage. In addition, Msuite2 shows both balanced and high performance in terms of mapping efficiency and accuracy, demonstrating high potential to facilitate the investigation and application of large-scale DNA methylation analysis in various biomedical studies. Msuite2 is freely available at https://github.com/hellosunking/Msuite2/.

Generation of CpG-Recoded Zika Virus Vaccine Candidates.

Experimental increase of cytosine-phosphate-guanine (CpG) dinucleotides in an RNA virus genome impairs infection. Beneficially, this weak infection may lead to robust antiviral host immunity providing a cutting-edge approach for vaccines. For example, we have recently demonstrated that recoded Zika virus variants with the increased CpG content showed considerable attenuated infection phenotypes and protection against lethal challenge in mice. Here, we describe the workflow for the design and generation of CpG-recoded Zika virus vaccine candidates. The workflow can be adapted for other viruses.

Evolutionary Signatures Governing the Codon Usage Bias in Coronaviruses and Their Implications for Viruses Infecting Various Bat Species.

Many viruses that cause serious diseases in humans and animals, including the betacoronaviruses (beta-CoVs), such as SARS-CoV, MERS-CoV, and the recently identified SARS-CoV-2, have natural reservoirs in bats. Because these viruses rely entirely on the host cellular machinery for survival, their evolution is likely to be guided by the link between the codon usage of the virus and that of its host. As a result, specific cellular microenvironments of the diverse hosts and/or host tissues imprint peculiar molecular signatures in virus genomes. Our study is aimed at deciphering some of these signatures. Using a variety of genetic methods we demonstrated that trends in codon usage across chiroptera-hosted CoVs are collaboratively driven by geographically different host-species and temporal-spatial distribution. We not only found that chiroptera-hosted CoVs are the ancestors of SARS-CoV-2, but we also revealed that SARS-CoV-2 has the codon usage characteristics similar to those seen in CoVs infecting the Rhinolophus sp. Surprisingly, the envelope gene of beta-CoVs infecting Rhinolophus sp., including SARS-CoV-2, had extremely high CpG levels, which appears to be an evolutionarily conserved trait. The dissection of the furin cleavage site of various CoVs infecting hosts revealed host-specific preferences for arginine codons; however, arginine is encoded by a wider variety of synonymous codons in the murine CoV (MHV-A59) furin cleavage site. Our findings also highlight the latent diversity of CoVs in mammals that has yet to be fully explored.

Evolutionary perspectives and adaptation dynamics of human seasonal influenza viruses from 2009 to 2019: An insight from codon usage.

The annually recurrent seasonal influenza viruses, namely, influenza A viruses (H1N1/pdm2009 and H3N2) and influenza B viruses, contribute substantially to human disease burden. Elucidation of host adaptation, population dynamics and evolutionary patterns of these viruses contribute to better control of current epidemic situation and bolster efforts towards pandemic preparedness. Present study has been addressed at unraveling the signatures of codon usage and dinucleotide distribution of these seasonal influenza viruses associating with their fitness and ongoing adaptive evolution in human population. Thorough analysis of codon usage adaptation revealed that H3N2 has been exhibited best adapted to human cellular system, which correlate with its highest epidemic intensity as compared with the other seasonal influenza viruses. CpG dinucleotide was found to be strongly avoided among the seasonal influenza viruses with more restraint among influenza B viruses than influenza A viruses, and might be accounted to the strategy of the viral pathogens in evading human immune signals. Dynamic scenes of ongoing evolution in codon usage and elimination of CpG motif among the viruses, which correlate with their distinct host adaption state, signifying the marked impact of selective force operational on the viral genomes, aimed at proficient circulation, enhanced fitness and successful infective manifestations in humans.

Targeted Restriction of Viral Gene Expression and Replication by the ZAP Antiviral System.

The zinc finger antiviral protein (ZAP) restricts the replication of a broad range of RNA and DNA viruses. ZAP directly binds viral RNA, targeting it for degradation and inhibiting its translation. While the full scope of RNA determinants involved in mediating selective ZAP activity is unclear, ZAP binds CpG dinucleotides, dictating at least part of its target specificity. ZAP interacts with many cellular proteins, although only a few have been demonstrated to be essential for its antiviral activity, including the 3'-5' exoribonuclease exosome complex, TRIM25, and KHNYN. In addition to inhibiting viral gene expression, ZAP also directly and indirectly targets a subset of cellular messenger RNAs to regulate the innate immune response. Overall, ZAP protects a cell from viral infection by restricting viral replication and regulating cellular gene expression. Further understanding of the ZAP antiviral system may allow for novel viral vaccine and anticancer therapy development.

CpG Methylation Profiles of HIV-1 Pro-Viral DNA in Individuals on ART.

The latent HIV-1 reservoir is comprised of stably integrated and intact proviruses with limited to no viral transcription. It has been proposed that latent infection may be maintained by methylation of pro-viral DNA. Here, for the first time, we investigate the cytosine methylation of a replication competent provirus (AMBI-1) found in a T cell clone in a donor on antiretroviral therapy (ART). Methylation profiles of the AMBI-1 provirus were compared to other proviruses in the same donor and in samples from three other individuals on ART, including proviruses isolated from lymph node mononuclear cells (LNMCs) and peripheral blood mononuclear cells (PBMCs). We also evaluated the apparent methylation of cytosines outside of CpG (i.e., CpH) motifs. We found no evidence for methylation in AMBI-1 or any other provirus tested within the 5' LTR promoter. In contrast, CpG methylation was observed in the env-tat-rev overlapping reading frame. In addition, we found evidence for differential provirus methylation in cells isolated from LNMCs vs. PBMCs in some individuals, possibly from the expansion of infected cell clones. Finally, we determined that apparent low-level methylation of CpH cytosines is consistent with occasional bisulfite reaction failures. In conclusion, our data do not support the proposition that latent HIV infection is associated with methylation of the HIV 5' LTR promoter.

Association of Zinc Finger Antiviral Protein Binding to Viral Genomic RNA with Attenuation of Replication of Echovirus 7.

Previous studies have implicated both zinc finger antiviral protein (ZAP) and oligoadenylate synthetase 3 (OAS3)/RNase L in the attenuation of RNA viruses with elevated CpG and UpA dinucleotides. Mechanisms and interrelationships between these two pathways were investigated using an echovirus 7 (E7) replicon with compositionally modified sequences inserted into the 3' untranslated region. ZAP and OAS3 immunoprecipitation (IP) assays provided complementary data on dinucleotide composition effects on binding. Elevated frequencies of alternative pyrimidine/purine (CpA and UpG) and reversed (GpC and ApU) dinucleotides showed no attenuating effect on replication or specific binding to ZAP by IP. However, the bases 3' and 5' of CpG motifs influenced replication and ZAP binding; UCGU enhanced CpG-mediated attenuation and ZAP binding, while A residues shielded CpGs from ZAP recognition. Attenuating effects of elevated frequencies of UpA on replication occurred independently of CpG dinucleotides and bound noncompetitively with CpG-enriched RNA, consistent with a separate recognition site from CpG. Remarkably, immunoprecipitation with OAS3 antibody reproduced the specific binding to CpG- and UpA-enriched RNA sequences. However, OAS3 and ZAP were coimmunoprecipitated in both ZAP and OAS3 IP and colocalized with E7 and stress granules (SGs) by confocal microscopy analysis of infected cells. ZAP's association with larger cellular complexes may mediate the recruitment of OAS3/RNase L, KHNYN, and other RNA degradation pathways.IMPORTANCE We recently discovered that the OAS3/RNase L antiviral pathway is essential for restriction of CpG- and UpA-enriched viruses, in addition to the requirement for zinc finger antiviral protein (ZAP). The current study provides evidence for the specific dinucleotide and wider recognition contexts associated with virus recognition and attenuation. It further documents the association of ZAP and OAS3 and association with stress granules and a wider protein interactome that may mediate antiviral effects in different cellular compartments. The study provides a striking reconceptualization of the pathways associated with this aspect of antiviral defense.

The Long-Term Evolutionary History of Gradual Reduction of CpG Dinucleotides in the SARS-CoV-2 Lineage.

Recent studies suggested that the fraction of CG dinucleotides (CpG) is severely reduced in the genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The CpG deficiency was predicted to be the adaptive response of the virus to evade degradation of the viral RNA by the antiviral zinc finger protein that specifically binds to CpG nucleotides. By comparing all representative genomes belonging to the genus Betacoronavirus, this study examined the potential time of origin of CpG depletion. The results of this investigation revealed a highly significant correlation between the proportions of CpG nucleotide (CpG content) of the betacoronavirus species and their times of divergence from SARS-CoV-2. Species that are distantly related to SARS-CoV-2 had much higher CpG contents than that of SARS-CoV-2. Conversely, closely related species had low CpG contents that are similar to or slightly higher than that of SARS-CoV-2. These results suggest a systematic and continuous reduction in the CpG content in the SARS-CoV-2 lineage that might have started since the Sarbecovirus + Hibecovirus clade separated from Nobecovirus, which was estimated to be 1213 years ago. This depletion was not found to be mediated by the GC contents of the genomes. Our results also showed that the depletion of CpG occurred at neutral positions of the genome as well as those under selection. The latter is evident from the progressive reduction in the proportion of arginine amino acid (coded by CpG dinucleotides) in the SARS-CoV-2 lineage over time. The results of this study suggest that shedding CpG nucleotides from their genome is a continuing process in this viral lineage, potentially to escape from their host defense mechanisms.

Pirahy virus: Identification of a new and potential emerging arbovirus in South Brazil.

Genomic and epidemiological surveillance are paramount for the discovery of new viruses with the potential to cross species barriers. Here, we present a new member of the genus Alphavirus found in Trichoprosopon and Wyeomia mosquitoes, tentatively named Pirahy virus (PIRAV). PIRAV was isolated from mosquito pools collected in a rural area of Piraí do Sul, South Brazil. In vitro assays revealed that PIRAV replicates and causes cytopathic effects in vertebrate cell lines such as Vero E6, SH-SY5Y, BHK-21 and UMNSAH/DF-1. Genomic signature analysis supports these results showing a dinucleotide and codon usage balance compatible with several hosts. Phylogenetic analyses placed PIRAV basal to the Venezuelan equine encephalitis complex. Genome analyses, electron microscopy, and biological characterization show findings that may alert for the emergence of a new arbovirus in South America.

Msuite: A High-Performance and Versatile DNA Methylation Data-Analysis Toolkit.

DNA methylation is a pervasive and important epigenetic regulator in mammalian genome. For DNA methylome profiling, emerging bisulfite-free methods have demonstrated desirable superiority over the conventional bisulfite-treatment-based approaches, although current analysis software could not make full use of their advantages. In this work, we present Msuite, an easy-to-use, all-in-one data-analysis toolkit. Msuite implements a unique 4-letter analysis mode specifically optimized for emerging protocols; it also integrates quality controls, methylation call, and data visualizations. Msuite demonstrates substantial performance improvements over current state-of-the-art tools as well as fruitful functionalities, thus holding the potential to serve as an optimal toolkit to facilitate DNA methylome studies. Source codes and testing datasets for Msuite are freely available at https://github.com/hellosunking/Msuite/.

Use of a small DNA virus model to investigate mechanisms of CpG dinucleotide-induced attenuation of virus replication.

Suppression of the CpG dinucleotide is widespread in RNA viruses infecting vertebrates and plants, and in the genomes of retroviruses and small mammalian DNA viruses. The functional basis for CpG suppression in the latter was investigated through the construction of mutants of the parvovirus, minute virus of mice (MVM) with increased CpG or TpA dinucleotides in the VP gene. CpG-high mutants displayed extraordinary attenuation in A9 cells compared to wild-type MVM (>six logs), while TpA elevation showed no replication effect. Attenuation was independent of Toll-like receptor 9 and STING-mediated DNA recognition pathways and unrelated to effects on translation efficiency. While translation from codon-optimized VP RNA was enhanced in a cell-free assay, MVM containing this sequence was highly attenuated. Further mutational analysis indicated that this arose through its increased numbers of CpG dinucleotides (7→70) and separately from its increased G+C content (42.3→57.4 %), which independently attenuated replication. CpG-high viruses showed impaired NS mRNA expression by qPCR and reduced NS and particularly VP protein expression detected by immunofluorescence and replication in A549 cells, effects reversed in zinc antiviral protein (ZAP) knockout cells, even though nuclear relocalization of VP remained defective. The demonstrated functional basis for CpG suppression in MVM and potentially other small DNA viruses and the observed intolerance of CpGs in coding sequences, even after codon optimization, has implications for the use of small DNA virus vectors in gene therapy and immunization.

Genome-wide evolution analysis reveals low CpG contents of fast-evolving genes and identifies antiviral microRNAs.

Noncoding RNAs (ncRNAs) play important roles in many biological processes and provide materials for evolutionary adaptations beyond protein-coding genes, such as in the arms race between the host and pathogen. However, currently, a comprehensive high-resolution analysis of primate genomes that includes the latest annotated ncRNAs is not available. Here, we developed a computational pipeline to estimate the selections that act on noncoding regions based on comparisons with a large number of reference sequences in introns adjacent to the interested regions. Our method yields result comparable with those of the established codon-based method and phyloP method for coding genes; thus, it provides a holistic framework for estimating the selection on the entire genome. We further showed that fast-evolving protein-coding genes and their corresponding 5' UTRs have a significantly lower frequency of the CpG dinucleotides than those evolving at an average pace, and these fast-evolving genes are enriched in the process of immunity and host defense. We also identified fast-evolving miRNAs with antiviral functions in cells. Our results provide a resource for high-resolution evolution analysis of the primate genomes.

KHNYN is essential for the zinc finger antiviral protein (ZAP) to restrict HIV-1 containing clustered CpG dinucleotides.

CpG dinucleotides are suppressed in most vertebrate RNA viruses, including HIV-1, and introducing CpGs into RNA virus genomes inhibits their replication. The zinc finger antiviral protein (ZAP) binds regions of viral RNA containing CpGs and targets them for degradation. ZAP does not have enzymatic activity and recruits other cellular proteins to inhibit viral replication. We found that KHNYN, a protein with no previously known function, interacts with ZAP. KHNYN overexpression selectively inhibits HIV-1 containing clustered CpG dinucleotides and this requires ZAP and its cofactor TRIM25. KHNYN requires both its KH-like domain and NYN endonuclease domain for antiviral activity. Crucially, depletion of KHNYN eliminated the deleterious effect of CpG dinucleotides on HIV-1 RNA abundance and infectious virus production and also enhanced the production of murine leukemia virus. Overall, we have identified KHNYN as a novel cofactor for ZAP to target CpG-containing retroviral RNA for degradation.

DamMet: ancient methylome mapping accounting for errors, true variants, and post-mortem DNA damage.

BACKGROUND: Recent computational advances in ancient DNA research have opened access to the detection of ancient DNA methylation footprints at the genome-wide scale. The most commonly used approach infers the methylation state of a given genomic region on the basis of the amount of nucleotide mis-incorporations observed at CpG dinucleotide sites. However, this approach overlooks a number of confounding factors, including the presence of sequencing errors and true variants. The scale and distribution of the inferred methylation measurements are also variable across samples, precluding direct comparisons. FINDINGS: Here, we present DamMet, an open-source software program retrieving maximum likelihood estimates of regional CpG methylation levels from ancient DNA sequencing data. It builds on a novel statistical model of post-mortem DNA damage for dinucleotides, accounting for sequencing errors, genotypes, and differential post-mortem cytosine deamination rates at both methylated and unmethylated sites. To validate DamMet, we extended gargammel, a sequence simulator for ancient DNA data, by introducing methylation-dependent features of post-mortem DNA decay. This new simulator provides direct validation of DamMet predictions. Additionally, the methylation levels inferred by DamMet were found to be correlated to those inferred by epiPALEOMIX and both on par and directly comparable to those measured from whole-genome bisulphite sequencing experiments of fresh tissues. CONCLUSIONS: DamMet provides genuine estimates for local DNA methylation levels in ancient individual genomes. The returned estimates are directly cross-sample comparable, and the software is available as an open-source C++ program hosted at https://gitlab.com/KHanghoj/DamMet along with a manual and tutorial.

The Effect of Single CpG Demethylation on the Pattern of DNA-Protein Binding.

Epidermal differentiation is a complex process and its regulation may involve epigenetic factors. Analysis of DNA methylation in 20 selected regions within the epidermal differentiation complex (EDC) gene cluster by targeted next-generation sequencing (NGS) detected no or only minor changes in methylation, mostly slight demethylation, occurring during the course of keratinocyte differentiation. However, a single CpG pair within the exon of the PGLYRP3 gene underwent a pronounced demethylation concomitant with an increase in PGLYRP3 expression. We have employed a DNA-affinity precipitation assay (DAPA) and mass spectrometry to examine changes in the composition of proteins that bind to DNA containing either methylated or unmethylated CpG. We found that the unmethylated probe attracted mostly RNA binding proteins, including splicing factors, which suggests that demethylation of this particular CpG may facilitate PGLYRP3 transcription and/or pre-mRNA splicing.