Search NASA⌕ Search

SEARCH · Search NASA

Results for “genomic analysis”

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 451 records · Page 25

A Rhodopseudomonas strain with a substantially smaller genome retains the core metabolic versatility of its genus

ABSTRACT Rhodopseudomonas are a group of phototrophic microbes with a marked metabolic versatility and flexibility that underpins their potential use in the production of value-added products, bioremediation, and plant growth promotion. Members of this group have an average genome size of about 5.5 Mb, but two closely related strains have genome sizes of about 4.0 Mb. To identify the types of genes missing in a reduced genome strain, we compared strain DSM127 with other Rhodopseudomonas isolates at the genomic and phenotypic levels. We found that DSM127 can grow as well as other members of the Rhodopseudomonas genus and retains most of their metabolic versatility, but it has many fewer genes associated with high-affinity transport of nutrients, iron uptake, nitrogen metabolism, and biodegradation of aromatic compounds. This analysis indicates genes that can be deleted in genome reduction campaigns and suggests that DSM127 could be a favorable choice for biotechnology applications using Rhodopseudomonas or as a strain that can be engineered further to reside in a specialized natural environment. IMPORTANCE Rhodopseudomonas are a cohort of phototrophic bacteria with broad metabolic versatility. Members of this group are present in diverse soil and water environments, and some strains are found associated with plants and have plant growth-promoting activity. Motivated by the idea that it may be possible to design bacteria with reduced genomes that can survive well only in a specific environment or that may be more metabolically efficient, we compared Rhodopseudomonas strains with typical genome sizes of about 5.5 Mb to a strain with a reduced genome size of 4.0 Mb. From this, we concluded that metabolic versatility is part of the identity of the Rhodopseudomonas group, but high-affinity transport genes and genes of apparent redundant function can be dispensed with.

59 BASIC BIOLOGICAL SCIENCES↗

From 2D to 4D: a containerized workflow and browser to explore dynamic chromatin architecture

Background Characterizing the physical organization of the genome is essential for understanding long-range gene regulation, chromatin compartmentalization, and epigenetic accessibility. Hi-C experiments generate two-dimensional (2D) genome-wide contact maps of chromatin interactions by capturing the spatial proximity between genomic loci, which reveal interaction frequencies but lack the spatial resolution needed to interpret the three-dimensional (3D) genome structure(s). Emerging evidence suggests that epigenetic regulation is closely linked to 3D genome architecture, and that structural changes over time (4D) drive key biological processes in development, disease, and environmental response. Thus, integrating 3D structure with functional data is critical for a more complete understanding of genome regulation. Previous work, most notably the 4DHiC chromosome modeling framework, has shown that physical multi-dimensional modeling approaches rooted in polymer physics and molecular dynamics can resolve these structures at biologically meaningful resolutions by integrating temporal Hi-C data with physical constraints to uncover dynamic chromosome reorganization. Thus, molecular dynamics simulations, constrained by Hi-C contact matrices, can resolve fine-scale structural changes and reveal functionally significant transitions in chromatin conformation. Results Herein, we present the 4D Genome Browser Workflow (4DGBWorkflow) and the 4D Genome Browser (4DGB). The algorithm is based on the 4DHiC method, and the containerized tool is an end-to-end workflow that can transform, filter, and view 4D epigenomics and chromatin datasets, allowing non-specialists to apply three-dimensional modeling principles to diverse datasets and experimental conditions. The software executes on a laptop running macOS, Linux or Windows. From input Hi-C files (.hic), the 4DGBWorkflow produces 3D reconstructions of chromosomes, integrates the reconstruction with track data (e.g., epigenetic marks, transcriptome profiles), and provides comparative visualization of the results in a single workflow. Conclusions The 4DGBWorkflow and 4D Genome Browser are open-source tools for comparative analysis and visualization of 4D chromosome datasets, including chromatin architecture and epigenomic signals. Automatic integration of Hi-C data with molecular dynamics democratizes the construction of time resolved 3D genome structures, simplifying complex simulations and data integration schemes.

3D Genome Browser↗

Inferring microbial interactions with their environment from genomic and metagenomic data

Microbial communities assemble through a complex set of interactions between microbes and their environment, and the resulting metabolic impact on the host ecosystem can be profound. Microbial activity is known to impact human health, plant growth, water quality, and soil carbon storage which has lead to the development of many approaches and products meant to manipulate the microbiome. In order to understand, predict, and improve microbial community engineering, genome-scale modeling techniques have been developed to translate genomic data into inferred microbial dynamics. However, these techniques rely heavily on simulation to draw conclusions which may vary with unknown parameters or initial conditions, rather than more robust qualitative analysis. To better understand microbial community dynamics using genome-scale modeling, we provide a tool to investigate the network of interactions between microbes and environmental metabolites over time. Using our previously developed algorithm for simulating microbial communities from genome-scale metabolic models (GSMs), we infer the set of microbe-metabolite interactions within a microbial community in a particular environment. Because these interactions depend on the available environmental metabolites, we refer to the networks that we infer as metabolically contextualized , and so name our tool MetConSIN: Met abolically Con textualized S pecies I nteraction N etworks.

59 BASIC BIOLOGICAL SCIENCES↗

Pseudodesulfovibrio mercurii sp. nov., a mercury-methylating bacterium isolated from sediment

The sulfate-reducing, mercury-methylating strain ND132 T was isolated from the brackish anaerobic bottom sediments of Chesapeake Bay, USA. Capable of high levels of mercury (Hg) methylation, ND132 T has been widely used as a model strain to study the process and to determine the genetic basis of Hg methylation. Originally called Desulfovibrio desulfuricans ND132 T on the basis of an early partial 16S rRNA sequence, the strain has never been formally described. Phylogenetic and physiological traits place this strain within the genus Pseudodesulfovibrio, in the recently reclassified phylum Desulfobacterota (formerly Deltaproteobacteria ). ND132 T is most closely related to Pseudodesulfovibrio hydrargyri BerOc1 T and Pseudodesulfovibrio indicus J2T. Analysis of average nucleotide identity (ANI) of whole-genome sequences showed roughly 88 % ANI between P. hydrargyri BerOc1T and ND132T, and 84 % similarity between ND132T and P. indicus J2 T . These cut-off scores <95 %, along with a multi-gene phylogenetic analysis of members of the family Desulfovibrionacea, and differences in physiology indicate that all three strains represent separate species. The Gram-stain-negative cells are vibrio-shaped, motile and not sporulated. ND132 T is a salt-tolerant mesophile with optimal growth in the laboratory at 32 °C, 2 % salinity, and pH 7.8. The DNA G+C content of the genomic DNA is 65.2 %. It is an incomplete oxidizer of short chain fatty acids, using lactate, pyruvate and fumarate with sulfate or sulfite as the terminal electron acceptors. ND132 T can respire fumarate using pyruvate as an electron donor. The major fatty acids are iso-C 15 :0 , anteiso-C 15 : 0 , iso-C 17 : 0 , iso-C 17 : 1 ω9c and anteiso-C 17 : 0 . We propose the classification of strain ND132 T (DSM 110689, ATCC TSD-224) as the type strain Pseudodesulfovibrio mercurii sp. nov.

59 BASIC BIOLOGICAL SCIENCES↗

ldrd_virus_work

This is a Python code base that takes openly-available genetic information on known viruses and performs supervised machine learning and feature importance analysis on the relationship of the viral genomes to the competence to infect humans or bind to a specific host cell receptor.

Reddy, Tyler [LANL]↗

Oncogenic Ras and ΔNp63α cooperate to recruit immunosuppressive polymorphonuclear myeloid-derived suppressor cells in a mouse model of squamous cancer pathogenesis

Introduction Amplification of human chromosome 3q26-29, which encodes oncoprotein ΔNp63 among other isoforms of the p63 family, is a feature common to squamous cell carcinomas (SCCs) of multiple tissue origins. Along with overexpression of ΔNp63, activation of the protooncogene, RAS , whether by overexpression or oncogenic mutation, is frequently observed in many cancers. In this study, analysis of transcriptome data from The Cancer Genome Atlas (TCGA) demonstrated that expression of TP63 mRNA , particularly ΔNp63 isoforms, and HRAS are significantly elevated in advanced squamous cell carcinomas of the head and neck (HNSCCs), suggesting pathological significance. However, how co-overexpressed ΔNp63 and HRAS affect the immunosuppressive tumor microenvironment (TME) is incompletely understood. Methods Here, we established and characterized an immune competent mouse model using primary keratinocytes with retroviral-mediated overexpression of ΔNp63α and constitutively activated HRAS (v-ras Ha G12R) to evaluate the role of these oncogenes in the immune TME. Results In this model, orthotopic grafting of wildtype syngeneic keratinocytes expressing both v-ras Ha and elevated levels of ΔNp63α consistently yield carcinomas in syngeneic hosts, while cells expressing v-ras Ha alone yield predominantly papillomas. We found that polymorphonuclear (PMN) myeloid cells, experimentally validated to be immunosuppressive and thus representing myeloid-derived suppressor cells (PMN-MDSCs), were significantly recruited into the TME of carcinomas arising early following orthotopic grafting of ΔNp63α/v-ras Ha -expressing keratinocytes. ΔNp63α/v-ras Ha -driven carcinomas expressed higher levels of chemokines implicated in recruitment of MDSCs compared to v-ras Ha -initiated tumors, providing a heretofore undescribed link between ΔNp63α/HRAS-driven carcinomas and the development of an immunosuppressive TME. Conclusion These results support the utilization of a genetic carcinogenesis model harboring specific genomic drivers of malignancy to study mechanisms underlying the development of local immunosuppression.

60 APPLIED LIFE SCIENCES↗

Genome-wide association study identifies four pan-ancestry loci for suicidal ideation in the Million Veteran Program

Suicidal ideation (SI) often precedes and predicts suicide attempt and death, is the most common suicidal phenotype and is over-represented in veterans. The genetic architecture of SI in the absence of suicide attempt (SA) is unknown, yet believed to have distinct and overlapping risk with other suicidal behaviors. We performed the first GWAS of SI without SA in the Million Veteran Program (MVP), identifying 99,814 SI cases from electronic health records without a history of SA or suicide death (SD) and 512,567 controls without SI, SA or SD. GWAS was performed separately in the four largest ancestry groups, controlling for sex, age and genetic substructure. Ancestry-specific results were combined via meta-analysis to identify pan-ancestry loci. Four genome-wide significant (GWS) loci were identified in the pan-ancestry meta-analysis with loci on chromosomes 6 and 9 associated with suicide attempt in an independent sample. Pan-ancestry gene-based analysis identified GWS associations with DRD2, DCC, FBXL19, BCL7C, CTF1, ANNK1, and EXD3. Gene-set analysis implicated synaptic and startle response pathways (q’s<0.05). European ancestry (EA) analysis identified GWS loci on chromosomes 6 and 9, as well as GWS gene associations in EXD3, DRD2, and DCC. No other ancestry-specific GWS results were identified, underscoring the need to increase representation of diverse individuals. The genetic correlation of SI and SA within MVP was high (r G = 0.87; p = 1.09e-50), as well as with post-traumatic stress disorder (PTSD; r G = 0.78; p = 1.98e-95) and major depressive disorder (MDD; r G = 0.78; p = 8.33e-83). Conditional analysis on PTSD and MDD attenuated most pan-ancestry and EA GWS signals for SI without SA to nominal significance, with the exception of EXD3 which remained GWS. Our novel findings support a polygenic and complex architecture for SI without SA which is largely shared with SA and overlaps with psychiatric conditions frequently comorbid with suicidal behaviors.

59 BASIC BIOLOGICAL SCIENCES↗

Population genomics provides insights into the genetic basis of adaptive evolution in the mushroom-forming fungus Lentinula edodes

Introduction: Mushroom-forming fungi comprise diverse species that develop complex multicellular structures. In cultivated species, both ecological adaptation and artificial selection have driven genome evolution. However, little is known about the connections among genotype, phenotype and adaptation in mushroom-forming fungi. Objectives: This study aimed to (1) uncover the population structure and demographic history of Lentinula edodes, (2) dissect the genetic basis of adaptive evolution in L. edodes, and (3) determine if genes related to fruiting body development are involved in adaptive evolution. Methods: We analyzed genomes and fruiting body-related traits (FBRTs) in 133 L. edodes strains and conducted RNA-seq analysis of fruiting body development in the YS69 strain. Combined methods of genomic scan for divergence, genome-wide association studies (GWAS), and RNA-seq were used to dissect the genetic basis of adaptive evolution. Results: We detected three distinct subgroups of L. edodes via single nucleotide polymorphisms, which showed robust phenotypic and temperature response differentiation and correlation with geographical distribution. Demographic history inference suggests that the subgroups diverged 36,871 generations ago. Moreover, L. edodes cultivars in China may have originated from the vicinity of Northeast China. A total of 942 genes were found to be related to genetic divergence by genomic scan, and 719 genes were identified to be candidates underlying FBRTs by GWAS. Integrating results of genomic scan and GWAS, 80 genes were detected to be related to phenotypic differentiation. A total of 364 genes related to fruiting body development were involved in genetic divergence and phenotypic differentiation. Conclusion: Adaptation to the local environment, especially temperature, triggered genetic divergence and phenotypic differentiation of L. edodes. A general model for genetic divergence and phenotypic differentiation during adaptive evolution in L. edodes, which involves in signal perception and transduction, transcriptional regulation, and fruiting body morphogenesis, was also integrated here.

59 BASIC BIOLOGICAL SCIENCES↗

Quantitative analysis of bristle number in Drosophila mutants identifies genes involved in neural development

BACKGROUND: The identification of the function of all genes that contribute to specific biological processes and complex traits is one of the major challenges in the postgenomic era. One approach is to employ forward genetic screens in genetically tractable model organisms. In Drosophila melanogaster, P element-mediated insertional mutagenesis is a versatile tool for the dissection of molecular pathways, and there is an ongoing effort to tag every gene with a P element insertion. However, the vast majority of P element insertion lines are viable and fertile as homozygotes and do not exhibit obvious phenotypic defects, perhaps because of the tendency for P elements to insert 5' of transcription units. Quantitative genetic analysis of subtle effects of P element mutations that have been induced in an isogenic background may be a highly efficient method for functional genome annotation. RESULTS: Here, we have tested the efficacy of this strategy by assessing the extent to which screening for quantitative effects of P elements on sensory bristle number can identify genes affecting neural development. We find that such quantitative screens uncover an unusually large number of genes that are known to function in neural development, as well as genes with yet uncharacterized effects on neural development, and novel loci. CONCLUSIONS: Our findings establish the use of quantitative trait analysis for functional genome annotation through forward genetics. Similar analyses of quantitative effects of P element insertions will facilitate our understanding of the genes affecting many other complex traits in Drosophila.

Non-NASA Center↗

Decomposing a San Francisco estuary microbiome using long-read metagenomics reveals species- and strain-level dominance from picoeukaryotes to viruses

ABSTRACT Although long-read sequencing has enabled obtaining high-quality and complete genomes from metagenomes, many challenges still remain to completely decompose a metagenome into its constituent prokaryotic and viral genomes. This study focuses on decomposing an estuarine metagenome to obtain a more accurate estimate of microbial diversity. To achieve this, we developed a new bead-based DNA extraction method, a novel bin refinement method, and obtained 150 Gbp of Nanopore sequencing. We estimate that there are ~500 bacterial and archaeal species in our sample and obtained 68 high-quality bins (>90% complete, <5% contamination, ≤5 contigs, contig length of >100 kbp, and all ribosomal and tRNA genes). We also obtained many contigs of picoeukaryotes, environmental DNA of larger eukaryotes such as mammals, and complete mitochondrial and chloroplast genomes and detected ~40,000 viral populations. Our analysis indicates that there are only a few strains that comprise most of the species abundances. IMPORTANCE Ocean and estuarine microbiomes play critical roles in global element cycling and ecosystem function. Despite the importance of these microbial communities, many species still have not been cultured in the lab. Environmental sequencing is the primary way the function and population dynamics of these communities can be studied. Long-read sequencing provides an avenue to overcome limitations of short-read technologies to obtain complete microbial genomes but comes with its own technical challenges, such as needed sequencing depth and obtaining high-quality DNA. We present here new sampling and bioinformatics methods to attempt decomposing an estuarine microbiome into its constituent genomes. Our results suggest there are only a few strains that comprise most of the species abundances from viruses to picoeukaryotes, and to fully decompose a metagenome of this diversity requires 1 Tbp of long-read sequencing. We anticipate that as long-read sequencing technologies continue to improve, less sequencing will be needed.

Lui, Lauren M.↗

Multiscale Simulation of Deployable Composite Structures

In this paper, a multiscale simulation method for analyzing deployable composite structures is presented. Effective shell properties of the composites are obtained based on Mechanics of Structure Genome (MSG) homogenization, and then implemented into a user-subroutine UGENS for structural simulation with shell elements in Abaqus. The column bending test(CBT) of a flat thin flexure and lenticular composite boom in a simplified deployer structure are studied for demonstration. A viscoelastic material model with direct integration is adopted in this paper. The CBT simulation shows good agreement with experiments during relaxation, while errors are observed when comparing residual deformation. It is shown that this CBT model can be calibrated to CBT test results. For the lenticular boom analysis, the complete process of flattening, coiling, stowage, deployment and recovery is simulated with the viscoelastic shell model.

Finite element analysis↗

Deuterostome Evolution: Large Data Set Analysis

This award allowed us to develop novel hardware for phylogenetics, collect genomic data and produce several phylogenies of deuterostome organisms, communicate the results publicly, release software into the public domain, publish textbooks and papers, and prepare for the next research projects. There are no resulting subject inventions to report. We review these activities in three sections: 1) Hardware and software and development; 2) Evolutionary biology research; 3) Our proposed future direction, predictive analysis of pathogens in support of the NASA mission.

Janies, Daniel↗

Rational Optimization of Microbial Processing for High Yield CO 2 -to-Isopropanol Conversion: Cooperative Research and Development Final Report, CRADA Number CRD-20-17114

This project focuses on the production of the fuel blendstock isopropanol using a CO 2 -fixing Clostridium by metabolic engineering and process optimizations. The project will initiate from a baseline isopropanol producer and pursue isopropanol production at high carbon-conversion efficiency. We will lead engineering work by in-depth pathway analyses including thermodynamics optimization, enzyme expense analysis, metabolic robustness analysis, and -omics analysis. The isopropanol production will be optimized via genome editing followed by fermentation optimizations. This project will deliver a novel microbial process that efficiently converts waste CO 2 to isopropanol at ~g/L titer level within 18-months. This project will layout a solid knowledge basis and technology platform for renewable CO 2 valorization to bio-blendstock that help achieve Co-Optima and Shell's goals.

09 BIOMASS FUELS↗

Threat Reduction Research Networks: Fostering Sustainable Collaborations Through Trainings for Genomics for Biosurveillance

Scientific research communities can be represented as heterogeneous or multidimensional networks encompassing multiple types of entities and relationships. These networks might include researchers, institutions, meetings, and publications, connected by relationships like authorship, employment, and attendance. We describe a method for efficiently and flexibly capturing, storing, and extracting information from multidimensional scientific networks using a graph database. The database structure is based on an ontology that captures allowable types of entities and relationships. This allows us to construct a variety of projections of the underlying multidimensional graph through database queries to answer specific research questions. We demonstrate this process through a study of the U.S. Biological Threat Reduction Program (BTRP), which seeks to develop Threat Reduction Networks to build and strengthen a sustainable international community of biosecurity, biosafety, and biosurveillance experts to address shared biological threat reduction challenges. Networks like these create connectional intelligence among researchers and institutions around the world, and are central to the concept of cooperative threat reduction. Our analysis focuses on a series of seven BTRP genome sequencing training workshops, showing how they created a growing network of participants and countries over time, which is also reflected in coauthorship relationships among attendees. By capturing concept and relationship hierarchies, our ontology-based approach allows us to pose general or specific questions about networks within the same framework. This approach can be applied to other research communities or multidimensional social networks to capture, analyze, and visualize different types of interactions and how they change over time.

59 BASIC BIOLOGICAL SCIENCES↗

Comparative Phenotyping of Two Commonly Used Chlamydomonas reinhardtii Background Strains: CC-1690 (21gr) and CC-5325 (The CLiP Mutant Library Background)

The unicellular green alga Chlamydomonas reinhardtii is an excellent model organism to investigate many essential cellular processes in photosynthetic eukaryotes. Two commonly used background strains of Chlamydomonas are CC-1690 and CC-5325. CC-1690, also called 21gr, has been used for the Chlamydomonas genome project and several transcriptome analyses. CC-5325 is the background strain for the Chlamydomonas Library Project (CLiP). Photosynthetic performance in CC-5325 has not been evaluated in comparison with CC-1690. Additionally, CC-5325 is often considered to be cell-wall deficient, although detailed analysis is missing. The circadian rhythms in CC-5325 are also unclear. To fill these knowledge gaps and facilitate the use of the CLiP mutant library for various screens, we performed phenotypic comparisons between CC-1690 and CC-5325. Our results showed that CC-5325 grew faster heterotrophically in dark and equally well in mixotrophic liquid medium as compared to CC-1690. CC-5325 had lower photosynthetic efficiency and was more heat-sensitive than CC-1690. Furthermore, CC-5325 had an intact cell wall which had comparable integrity to that in CC-1690 but appeared to have reduced thickness. Additionally, CC-5325 could perform phototaxis, but could not maintain a sustained circadian rhythm of phototaxis as CC1690 did. Finally, in comparison to CC-1690, CC-5325 had longer cilia in the medium with acetate but slower swimming speed in the medium without nitrogen and acetate. Our results will be useful for researchers in the Chlamydomonas community to choose suitable background strains for mutant analysis and employ the CLiP mutant library for genome-wide mutant screens under appropriate conditions, especially in the areas of photosynthesis, thermotolerance, cell wall, and circadian rhythms.

59 BASIC BIOLOGICAL SCIENCES↗

Impact of genotype-calling methodologies on genome-wide association and genomic prediction in polyploids

This dataset contains all data used in the paper "Impact of genotype-calling methodologies on genome-wide association and genomic prediction in polyploids". The dataset includes genotypes and phenotypic data from two autotetraploid species Miscanthus sacchariflorus and Vaccinium corymbosum that was used used for genome wide association studies and genomic prediction and the scripts used in the analysis. In this V2, 2 files have the raw data are added: "Miscanthus_sacchariflorus_RADSeq.vcf" is the VCF file with the raw SNP calls of the Miscanthus sacchariflorus data used for genotype calling using the 6 genotype calling methods. "Blueberry_data_read_depths.RData" is the a RData file with the read depth data that was used for genotype calling in the Blueberry dataset.

allelic dosage↗

Biomass yields, reproductive fertility, compositional analysis, and genetic diversity of newly developed triploid giant miscanthus hybrids

Abstract Miscanthus × giganteus (giant miscanthus), first found as a naturally occurring hybrid, has shown promise as a bioenergy/biomass crop throughout much of the temperate world. This allotriploid (2 n = 3 x = 57) hybrid resulted from a cross between tetraploid Miscanthus sacchariflorus (2 n = 4 x = 76) and diploid Miscanthus sinensis (2 n = 2 x = 38) and is particularly desirable due to its low fertility that minimizes reseeding and potential invasiveness. However, there is limited genetic diversity in commonly grown cultivars of triploid M. × giganteus and breeding and development efforts to improve and domesticate this crop have been minimal. Here, we report on newly developed M. × giganteus hybrids compared with the industry standard M. × giganteus '1993‐1780'. Dry biomass yields of new hybrids ranged from 19.5 to 32.4 Mg/ha/year for the fourth growing season, compared with 21.0 Mg/ha/year for M. × giganteus '1993‐1780'. Plant reproductive fertility remained low for all accessions with overall fertility [(seed set × seed germination)/100] ranging from 0.3% to 4.5% for new hybrids compared to 0.4% for M. × giganteus '1993‐1780'. Culm density and height varied among accessions and were positively correlated with increased biomass. Based on compositional analyses, theoretical ethanol yields ranged from 9, 740 to 16,278 L/ha/year for new hybrids compared to 10,406 L/ha/year for M. × giganteus '1993‐1780'. Relative feed value indices were low overall and ranged between 66.0 and 72.8 for new hybrids compared to M. × giganteus '1993‐1780' with 71.3. The genetic diversity of new hybrids, compared with existing cultivars, was characterized using whole genome sequences. Based on pair‐wise distances, cluster analysis clearly showed increased diversity of new hybrids compared with earlier selections. These results document new triploid hybrids of M. × giganteus with enhanced biomass and theoretical ethanol yields in combination with broader genetic diversity and lowreproductive fertility.

Touchell, Darren H.↗

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↗