Search NASA⌕ Search

SEARCH · Search NASA

Results for “single cell genomics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Single-cell RNA sequencing of batch Chlamydomonas cultures reveals heterogeneity in their diurnal cycle phase

The photosynthetic unicellular alga Chlamydomonas (Chlamydomonas reinhardtii) is a versatile reference for algal biology because of its ease of culture in the laboratory. Genomic and systems biology approaches have previously described transcriptome responses to environmental changes using bulk data, thus representing the average behavior from pools of cells. Here, we apply single-cell RNA sequencing (scRNA-seq) to probe the heterogeneity of Chlamydomonas cell populations under three environments and in two genotypes differing by the presence of a cell wall. First, we determined that RNA can be extracted from single algal cells with or without a cell wall, offering the possibility to sample natural algal communities. Second, scRNA-seq successfully separated single cells into nonoverlapping cell clusters according to their growth conditions. Cells exposed to iron or nitrogen deficiency were easily distinguished despite a shared tendency to arrest photosynthesis and cell division to economize resources. Notably, these groups of cells not only recapitulated known patterns observed with bulk RNA-seq but also revealed their inherent heterogeneity. A substantial source of variation between cells originated from their endogenous diurnal phase, although cultures were grown in constant light. We exploited this result to show that circadian iron responses may be conserved from algae to land plants. We document experimentally that bulk RNA-seq data represent an average of typically hidden heterogeneity in the population.

59 BASIC BIOLOGICAL SCIENCES↗

Unraveling Vulnerabilities in Endocrine Therapy-Resistant HER2+/ER+ Breast Cancer

Abstract Breast tumors overexpressing human epidermal growth factor receptor (HER2) confer intrinsic resistance to endocrine therapy (ET), and patients with HER2/estrogen receptor–positive (HER2+/ER+) breast cancer (BCa) are less responsive to ET than HER2–/ER+. However, real-world evidence reveals that a large subset of patients with HER2+/ER+ receive ET as monotherapy, positioning this treatment pattern as a clinical challenge. In the present study, we developed and characterized 2 in vitro models of ET-resistant (ETR) HER2+/ER+ BCa to identify possible therapeutic vulnerabilities. To mimic ETR to aromatase inhibitors (AIs), we developed 2 long-term estrogen deprivation (LTED) cell lines from BT-474 (BT474) and MDA-MB-361 (MM361). Growth assays, PAM50 subtyping, and genomic and transcriptomic analyses, followed by validation and functional studies, were used to identify targetable differences between ET-responsive parental and ETR-LTED HER2+/ER+ cells. Compared to their parental cells, MM361 LTEDs grew faster, lost ER, and increased HER2 expression, whereas BT474 LTEDs grew slower and maintained ER and HER2 expression. Both LTED variants had reduced responsiveness to fulvestrant. Whole-genome sequencing of aggressive MM361 LTEDs identified mutations in genes encoding transcription factors and chromatin modifiers. Single-cell RNA sequencing demonstrated a shift towards non-luminal phenotypes, and revealed metabolic remodeling of MM361 LTEDs, with upregulated lipid metabolism and ferroptosis-associated antioxidant genes, including GPX4. Combining a GPX4 inhibitor with anti-HER2 agents induced significant cell death in both MM361 and BT474 LTEDs. The BT474 and MM361 AI-resistant models capture distinct phenotypes of HER2+/ER+ BCa and identify altered lipid metabolism and ferroptosis remodeling as vulnerabilities of this type of ETR BCa.

60 APPLIED LIFE SCIENCES↗

High-throughput chemical imaging for optimizing biofuel synthesis using synthetic biology (Final Technical Report)

Fatty acids can be produced biosynthetically in microbes and these compounds can serve as precursors to biodiesels and other high value oleochemicals. However, progress on engineering fatty acid biosynthesis, and biofuel synthesis more generally, has been hindered by current quantification methods that are either indirect or not amenable to high-throughput or single-cell resolution screening. In this project, we assembled an interdisciplinary team with complimentary expertise in synthetic biology and microscopy, metabolic engineering, and chemical imaging to address these challenges. We used chemical imaging to directly measure lipid biosynthesis in Escherichia coli engineered to produce fatty acids, obtaining detailed single-cell resolution measurements. We deployed stimulated Raman scattering (SRS) microscopy in concert with multiplexed genome engineering and gene circuit design strategies from synthetic biology to optimize production of fatty acids. These results provided novel insight into cell-to-cell heterogeneity present in biofuel production strains. In addition, we introduced new chemical imaging methods which are label-free and do not require fluorescent reporters. These efforts were complemented by other studied developing foundational tools for regulation and control, which offer excellent potential for advancing researchers’ ability to rapidly design, build, and test strains for enhanced biofuel synthesis.

09 BIOMASS FUELS↗

Dynamics of phage-host interactions in Bacteroides fragilis resolved by single-cell transcriptomics

The interactions between lytic phages and their hosts are typically studied in bulk culture, which obscures cell-cell differences in infection susceptibility or expression of protective factors. Here, we use bacterial single-cell RNA sequencing to profile the transcriptomes of ~50,000 cells from cultures of a human pathobiont, Bacteroides fragilis, infected with a lytic bacteriophage. From a single sampling, we quantified the asynchronous progression of phage infection in individual bacterial cells and reconstructed the infection timeline, characterizing both host and phage transcriptomic changes as infection unfolded. Further, we discovered phenotypic subpopulations of bacteria that remained uninfected. Each cell’s vulnerability to phage infection was influenced by expression of multiple genetic loci, most prominently phase-variable capsular polysaccharide (CPS) biosynthesis pathways and an operon predicted to encode fimbrial genes. These findings uncovered genome-wide phase variation and stochasticity that enable bacterial survival and re-growth without acquiring additional mutations. Overall, we establish bacterial single-cell RNA sequencing as a powerful platform for investigating the dynamics of host-phage interactions and revealing the roles of phase variation and stochasticity in bacterial defenses.

Bacteria↗

Spatial co-transcriptomics reveals discrete stages of the arbuscular mycorrhizal symbiosis

Abstract The symbiotic interaction of plants with arbuscular mycorrhizal (AM) fungi is ancient and widespread. Plants provide AM fungi with carbon in exchange for nutrients and water, making this interaction a prime target for crop improvement. However, plant–fungal interactions are restricted to a small subset of root cells, precluding the application of most conventional functional genomic techniques to study the molecular bases of these interactions. Here we used single-nucleus and spatial RNA sequencing to explore both Medicago truncatula and Rhizophagus irregularis transcriptomes in AM symbiosis at cellular and spatial resolution. Integrated, spatially registered single-cell maps revealed infected and uninfected plant root cell types. We observed that cortex cells exhibit distinct transcriptome profiles during different stages of colonization by AM fungi, indicating dynamic interplay between both organisms during establishment of the cellular interface enabling successful symbiosis. Our study provides insight into a symbiotic relationship of major agricultural and environmental importance and demonstrates a paradigm combining single-cell and spatial transcriptomics for the analysis of complex organismal interactions.

59 BASIC BIOLOGICAL SCIENCES↗

Prediction of condition-specific regulatory genes using machine learning

Recent advances in genomic technologies have generated data on large-scale protein–DNA interactions and open chromatin regions for many eukaryotic species. How to identify condition-specific functions of transcription factors using these data has become a major challenge in genomic research. To solve this problem, we have developed a method called ConSReg, which provides a novel approach to integrate regulatory genomic data into predictive machine learning models of key regulatory genes. Using Arabidopsis as a model system, we tested our approach to identify regulatory genes in data sets from single cell gene expression and from abiotic stress treatments. Our results showed that ConSReg accurately predicted transcription factors that regulate differentially expressed genes with an average auROC of 0.84, which is 23.5–25% better than enrichment-based approaches. To further validate the performance of ConSReg, we analyzed an independent data set related to plant nitrogen responses. ConSReg provided better rankings of the correct transcription factors in 61.7% of cases, which is three times better than other plant tools. We applied ConSReg to Arabidopsis single cell RNA-seq data, successfully identifying candidate regulatory genes that control cell wall formation. Our methods provide a new approach to define candidate regulatory genes using integrated genomic data in plants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Perspectives for self-driving labs in synthetic biology

Self-driving labs (SDLs) combine fully automated experiments and data collection with artificial intelligence (AI) and control algorithms that decide not only the set of parameters for the next experiment, but also potentially which scientific hypotheses to test. Taken to their ultimate expression, SDLs could usher a new paradigm of scientific research, where the world is probed, interpreted, and explained by machines for human benefit. Whereas there are functioning SDLs in the fields of chemistry and materials science, we contend that synthetic biology provides a unique opportunity since the genome provides a single, easily accessible, target for affecting the incredibly wide repertoire of biological cell behavior. Since they can provide large amounts of high-quality data, SDLs can be a platform for AI to develop approaches to systematically convert data into scientific knowledge systems. These knowledge systems can be used both to understand the biological world and to design bioengineered systems to fit a desired specification (inverse design). However, the level of investment required for the creation of biological SDLs is only warranted if directed towards solving difficult and enabling biological questions. Here, we discuss challenges and opportunities in creating SDLs for synthetic biology.

59 BASIC BIOLOGICAL SCIENCES↗

CRISPR–Cas9-mediated nuclear transport and genomic integration of nanostructured genes in human primary cells

DNA nanostructures are a promising tool to deliver molecular payloads to cells. DNA origami structures, where long single-stranded DNA is folded into a compact nanostructure, present an attractive approach to package genes; however, effective delivery of genetic material into cell nuclei has remained a critical challenge. Here, we describe the use of DNA nanostructures encoding an intact human gene and a fluorescent protein encoding gene as compact templates for gene integration by CRISPR-mediated homology-directed repair (HDR). Our design includes CRISPR–Cas9 ribonucleoprotein binding sites on DNA nanostructures to increase shuttling into the nucleus. We demonstrate efficient shuttling and genomic integration of DNA nanostructures using transfection and electroporation. These nanostructured templates display lower toxicity and higher insertion efficiency compared to unstructured double-stranded DNA templates in human primary cells. Furthermore, our study validates virus-like particles as an efficient method of DNA nanostructure delivery, opening the possibility of delivering nanostructures in vivo to specific cell types. Together, these results provide new approaches to gene delivery with DNA nanostructures and establish their use as HDR templates, exploiting both their design features and their ability to encode genetic information. This work also opens a door to translate other DNA nanodevice functions, such as biosensing, into cell nuclei.

59 BASIC BIOLOGICAL SCIENCES↗

Novel candidate genes for lignin structure identified through genome-wide association study of naturally varying Populus trichocarpa

Populus is a promising lignocellulosic feedstock for biofuels and bioproducts. However, the cell wall biopolymer lignin is a major barrier in conversion of biomass to biofuels. To investigate the variability and underlying genetic basis of the complex structure of lignin, a population of 409 three-year-old, naturally varying Populus trichocarpa genotypes were characterized by heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR). A subsequent genome-wide association study (GWAS) was conducted using approximately 8.3 million single nucleotide polymorphisms (SNPs), which identified 756 genes that were significantly associated (−log 10 ( p -value)>6) with at least one lignin phenotype. Several promising candidate genes were identified, many of which have not previously been reported to be associated with lignin or cell wall biosynthesis. These results provide a resource for gaining insights into the molecular mechanisms of lignin biosynthesis and new targets for future genetic improvement in poplar.

59 BASIC BIOLOGICAL SCIENCES↗

Comparison of repair of DNA double-strand breaks in identical sequences in primary human fibroblast and immortal hamster-human hybrid cells harboring a single copy of human chromosome 11

We have optimized a pulsed-field gel electrophoresis assay that measures induction and repair of double-strand breaks (DSBs) in specific regions of the genome (Lobrich et al., Proc. Natl. Acad. Sci. USA 92, 12050-12054, 1995). The increased sensitivity resulting from these improvements makes it possible to analyze the size distribution of broken DNA molecules immediately after the introduction of DSBs and after repair incubation. This analysis shows that the distribution of broken DNA pieces after exposure to sparsely ionizing radiation is consistent with the distribution expected from randomly induced DSBs. It is apparent from the distribution of rejoined DNA pieces after repair incubation that DNA ends continue to rejoin between 3 and 24 h postirradiation and that some of these rejoining events are in fact misrejoining events, since novel restriction fragments both larger and smaller than the original fragment are generated after repair. This improved assay was also used to study the kinetics of DSB rejoining and the extent of misrejoining in identical DNA sequences in human GM38 cells and human-hamster hybrid A(L) cells containing a single human chromosome 11. Despite the numerous differences between these cells, which include species and tissue of origin, levels of TP53, expression of telomerase, and the presence or absence of a homologous chromosome for the restriction fragments examined, the kinetics of rejoining of radiation-induced DSBs and the extent of misrejoining were similar in the two cell lines when studied in the G(1) phase of the cell cycle. Furthermore, DSBs were removed from the single-copy human chromosome in the hamster A(L) cells with similar kinetics and misrejoining frequency as at a locus on this hybrid's CHO chromosomes.

NASA Program Biomedical Research and Countermeasur↗

Immunogenicity and therapeutic targeting of a public neoantigen derived from mutated PIK3CA

Public neoantigens (NeoAgs) represent an elite class of shared cancer-specific epitopes derived from recurrently mutated driver genes. Here we describe a high-throughput platform combining single-cell transcriptomic and T cell receptor (TCR) sequencing to establish whether mutant PIK3CA, among the most frequently genomically altered driver oncogenes, generates an immunogenic public NeoAg. Using this strategy, we developed a panel of TCRs that recognize an endogenously processed neopeptide encompassing a common PIK3CA hotspot mutation restricted by the prevalent human leukocyte antigen (HLA)-A*03:01 allele. Mechanistically, immunogenicity to this public NeoAg arises from enhanced neopeptide/HLA complex stability caused by a preferred HLA anchor substitution. Structural studies indicated that the HLA-bound neopeptide presents a comparatively ‘featureless’ surface dominated by the peptide’s backbone. To bind this epitope with high specificity and affinity, we discovered that a lead TCR clinical candidate engages the neopeptide through an extended interface facilitated by an unusually long CDR3β loop. In patients with diverse malignancies, we observed NeoAg clonal conservation and spontaneous immunogenicity to the neoepitope. Finally, adoptive transfer of TCR-engineered T cells led to tumor regression in vivo in mice bearing PIK3CA-mutant tumors but not wild-type PIK3CA tumors. Together, these findings establish the immunogenicity and therapeutic potential of a mutant PIK3CA-derived public NeoAg.

60 APPLIED LIFE SCIENCES↗

Clustered DNA damages induced in isolated DNA and in human cells by low doses of ionizing radiation

Clustered DNA damages-two or more closely spaced damages (strand breaks, abasic sites, or oxidized bases) on opposing strands-are suspects as critical lesions producing lethal and mutagenic effects of ionizing radiation. However, as a result of the lack of methods for measuring damage clusters induced by ionizing radiation in genomic DNA, neither the frequencies of their production by physiological doses of radiation, nor their repairability, nor their biological effects are known. On the basis of methods that we developed for quantitating damages in large DNAs, we have devised and validated a way of measuring ionizing radiation-induced clustered lesions in genomic DNA, including DNA from human cells. DNA is treated with an endonuclease that induces a single-strand cleavage at an oxidized base or abasic site. If there are two closely spaced damages on opposing strands, such cleavage will reduce the size of the DNA on a nondenaturing gel. We show that ionizing radiation does induce clustered DNA damages containing abasic sites, oxidized purines, or oxidized pyrimidines. Further, the frequency of each of these cluster classes is comparable to that of frank double-strand breaks; among all complex damages induced by ionizing radiation, double-strand breaks are only about 20%, with other clustered damage constituting some 80%. We also show that even low doses (0.1-1 Gy) of high linear energy transfer ionizing radiation induce clustered damages in human cells.

NASA Discipline Radiation Health↗

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Abstract Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN)cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genesMAFB, PHOX2A, CHN1, andEBF3– as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Science & Technology - Other Topics↗

Genomic instability, bystander effect, cytoplasmic irradiation and other phenomena that may achieve fame without fortune

The possible risk of induced malignancies in astronauts, as a consequence of the radiation environment in space, is a factor of concern for long term missions. Cancer risk estimates for high doses of low LET radiation are available from the epidemiological studies of the A-bomb survivors. Cancer risks at lower doses cannot be detected in epidemiological studies and must be inferred by extrapolation from the high dose risks. The standard setting bodies, such as the ICRP recommend a linear, no-threshold extrapolation of risks from high to low doses, but this is controversial. A study of mechanisms of carcinogenesis may shed some light on the validity of a linear extrapolation. The multi-step nature of carcinogenesis suggests that the role of radiation may be to induce a mutation leading to a mutator phenotype. High energy Fe ions, such as those encountered in space are highly effective in inducing genomic instability. Experiments involving the single particle microbeam have demonstrated a "bystander effect", ie a biological effect in cells not themselves hit, but in close proximity to those that are, as well as the induction of mutations in cells where only the cytoplasm, and not the nucleus, have been traversed by a charged particle. These recent experiments cast doubt on the validity of a simple linear extrapolation, but the data are so far fragmentary and conflicting. More studies are necessary. While mechanistic studies cannot replace epidemiology as a source of quantitative risk estimates, they may shed some light on the shape of the dose response relationship and therefore on the limitations of a linear extrapolation to low doses.

NASA Discipline Radiation Health↗

Selective Whole-Genome Amplification as a Tool to Enrich Specimens with Low Treponema pallidum Genomic DNA Copies for Whole-Genome Sequencing

Downstream next-generation sequencing (NGS) of the syphilis spirochete Treponema pallidum subspecies pallidum (T. pallidum) is hindered by low bacterial loads and the overwhelming presence of background metagenomic DNA in clinical specimens. In this study, we investigated selective whole-genome amplification (SWGA) utilizing multiple displacement amplification (MDA) in conjunction with custom oligonucleotides with an increased specificity for the T. pallidum genome and the capture and removal of 5'-C-phosphate-G-3' (CpG) methylated host DNA using the NEBNext Microbiome DNA enrichment kit followed by MDA with the REPLI-g single cell kit as enrichment methods to improve the yields of T. pallidum DNA in isolates and lesion specimens from syphilis patients. Sequencing was performed using the Illumina MiSeq v2 500 cycle or NovaSeq 6000 SP platform. These two enrichment methods led to 93 to 98% genome coverage at 5 reads/site in 5 clinical specimens from the United States and rabbit-propagated isolates, containing >14 T. pallidum genomic copies/μL of sample for SWGA and >129 genomic copies/μL for CpG methylation capture with MDA. Variant analysis using sequencing data derived from SWGA-enriched specimens showed that all 5 clinical strains had the A2058G mutation associated with azithromycin resistance. SWGA is a robust method that allows direct whole-genome sequencing (WGS) of specimens containing very low numbers of T. pallidum, which has been challenging until now.

59 BASIC BIOLOGICAL SCIENCES↗

Marine DNA methylation patterns are associated with microbial community composition and inform virus-host dynamics

Background: DNA methylation in prokaryotes is involved in many different cellular processes including cell cycle regulation and defense against viruses. To date, most prokaryotic methylation systems have been studied in culturable microorganisms, resulting in a limited understanding of DNA methylation from a microbial ecology perspective. Here, we analyze the distribution patterns of several microbial epigenetics marks in the ocean microbiome through genome-centric metagenomics across all domains of life. Results: We reconstructed 15,056 viral, 252 prokaryotic, 56 giant viral, and 6 eukaryotic metagenome-assembled genomes from northwest Pacific Ocean seawater samples using short- and long-read sequencing approaches. These metagenome-derived genomes mostly represented novel taxa, and recruited a majority of reads. Thanks to single-molecule real-time (SMRT) sequencing technology, base modification could also be detected for these genomes. This showed that DNA methylation can readily be detected across dominant oceanic bacterial, archaeal, and viral populations, and microbial epigenetic changes correlate with population differentiation. Furthermore, our genome-wide epigenetic analysis of Pelagibacter suggests that GANTC, a DNA methyltransferase target motif, is related to the cell cycle and is affected by environmental conditions. Yet, the presence of this motif also partitions the phylogeny of the Pelagibacter phages, possibly hinting at a competitive co-evolutionary history and multiple effects of a single methylation mark. Conclusions: Overall, this study elucidates that DNA methylation patterns are associated with ecological changes and virus-host dynamics in the ocean microbiome.

59 BASIC BIOLOGICAL SCIENCES↗

Mathematical modeling of hepatitis C RNA replication, exosome secretion and virus release

Hepatitis C virus (HCV) causes acute hepatitis C and can lead to life-threatening complications if it becomes chronic. The HCV genome is a single plus strand of RNA. Its intracellular replication is a spatiotemporally coordinated process of RNA translation upon cell infection, RNA synthesis within a replication compartment, and virus particle production. While HCV is mainly transmitted via mature infectious virus particles, it has also been suggested that HCV-infected cells can secrete HCV RNA carrying exosomes that can infect cells in a receptor independent manner. In order to gain insight into these two routes of transmission, we developed a series of intracellular HCV replication models that include HCV RNA secretion and/or virus assembly and release. Fitting our models to in vitro data, in which cells were infected with HCV, suggests that initially most secreted HCV RNA derives from intracellular cytosolic plus-strand RNA, but subsequently secreted HCV RNA derives equally from the cytoplasm and the replication compartments. Furthermore, our model fits to the data suggest that the rate of virus assembly and release is limited by host cell resources. Including the effects of direct acting antivirals in our models, we found that in spite of decreasing intracellular HCV RNA and extracellular virus concentration, low level HCV RNA secretion may continue as long as intracellular RNA is available. This may possibly explain the presence of detectable levels of plasma HCV RNA at the end of treatment even in patients that ultimately attain a sustained virologic response.

59 BASIC BIOLOGICAL SCIENCES↗

A minimal complex of KHNYN and zinc-finger antiviral protein binds and degrades single-stranded RNA

Detecting viral infection is a key role of the innate immune system. The genomes of some RNA viruses have a high CpG dinucleotide content relative to most vertebrate cell RNAs, making CpGs a molecular marker of infection. The human zinc-finger antiviral protein (ZAP) recognizes CpG, mediates clearance of the foreign CpG-rich RNA, and causes attenuation of CpG-rich RNA viruses. While ZAP binds RNA, it lacks enzymatic activity that might be responsible for RNA degradation and thus requires interacting cofactors for its function. One of these cofactors, KHNYN, has a predicted nuclease domain. Using biochemical approaches, we found that the KHNYN NYN domain is a single-stranded RNA ribonuclease that does not have sequence specificity and digests RNA with or without CpG dinucleotides equivalently in vitro. We show that unlike most KH domains, the KHNYN KH domain does not bind RNA. Indeed, a crystal structure of the KH region revealed a double-KH domain with a negatively charged surface that accounts for the lack of RNA binding. Rather, the KHNYN C-terminal domain (CTD) interacts with the ZAP RNA-binding domain (RBD) to provide target RNA specificity. We define a minimal complex composed of the ZAP RBD and the KHNYN NYN-CTD and use a fluorescence polarization assay to propose a model for how this complex interacts with a CpG dinucleotide-containing RNA. In the context of the cell, this module would represent the minimum ZAP and KHNYN domains required for CpG-recognition and ribonuclease activity essential for attenuation of viruses with clusters of CpG dinucleotides.

Yeoh, Zoe C. (ORCID:0000000226949068)↗