Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequence 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Recommendations for Uniform Variant Calling of SARS-CoV-2 Genome Sequence across Bioinformatic Workflows

Genomic sequencing of clinical samples to identify emerging variants of SARS-CoV-2 has been a key public health tool for curbing the spread of the virus. As a result, an unprecedented number of SARS-CoV-2 genomes were sequenced during the COVID-19 pandemic, which allowed for rapid identification of genetic variants, enabling the timely design and testing of therapies and deployment of new vaccine formulations to combat the new variants. However, despite the technological advances of deep sequencing, the analysis of the raw sequence data generated globally is neither standardized nor consistent, leading to vastly disparate sequences that may impact identification of variants. Here, we show that for both Illumina and Oxford Nanopore sequencing platforms, downstream bioinformatic protocols used by industry, government, and academic groups resulted in different virus sequences from same sample. These bioinformatic workflows produced consensus genomes with differences in single nucleotide polymorphisms, inclusion and exclusion of insertions, and/or deletions, despite using the same raw sequence as input datasets. Here, we compared and characterized such discrepancies and propose a specific suite of parameters and protocols that should be adopted across the field. Consistent results from bioinformatic workflows are fundamental to SARS-CoV-2 and future pathogen surveillance efforts, including pandemic preparation, to allow for a data-driven and timely public health response.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Attosecond Time Delay Trends across the Isoelectronic Noble Gas Sequence

The analysis and measurement of Wigner time delays can provide detailed information about the electronic environment within and around atomic and molecular systems, with one the key differences being the lack of a long-range potential after a halogen ion undergoes photoionization. In this work, we use relativistic random-phase approximation to calculate the average Wigner delay from the highest occupied subshells of the atomic pairings (2p, 2s in Fluorine, Neon), (3p, 3s in Chlorine, Argon), (4p, 4s, 3d, in Bromine, Krypton), and (5p, 5s, 4d in Iodine, Xenon). The qualitative behaviors of the Wigner delays between the isoelectronic pairings were found to be similar in nature, with the only large differences occurring at photoelectron energies less than 20 eV and around Cooper minima. Interestingly, the relative shift in Wigner time delays between negatively charged halogens and noble gases decreases as atomic mass increases. All atomic pairings show large differences at low energies, with noble gas atoms showing large positive Wigner delays, while negatively charged halogen ions show negative delays. The implications for photoionization studies in halide-containing molecules is also discussed.

Grafstrom, Brock↗
DOE OSTI2017

Anomaly Detection in Moving-Camera Video Sequences Using Principal Subspace Analysis

This study presents a family of algorithms based on sparse decompositions that detect anomalies in video sequences obtained from slow moving cameras. These algorithms start by computing the union of subspaces that best represents all the frames from a reference (anomaly free) video as a low-rank projection plus a sparse residue. Then, they perform a low-rank representation of a target (possibly anomalous) video by taking advantage of both the union of subspaces and the sparse residue computed from the reference video. Such algorithms provide good detection results while at the same time obviating the need for previous video synchronization. However, this is obtained at the cost of a large computational complexity, which hinders their applicability. Another contribution of this paper approaches this problem by using intrinsic properties of the obtained data representation in order to restrict the search space to the most relevant subspaces, providing computational complexity gains of up to two orders of magnitude. The developed algorithms are shown to cope well with videos acquired in challenging scenarios, as verified by the analysis of 59 videos from the VDAO database that comprises videos with abandoned objects in a cluttered industrial scenario.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Next-generation sequencing-based miRNA expression analysis in Parp1-deficient embryonic stem cell-derived exosomes

Highlights: • This is the first report for the extracellular role of Parp1 via exosomal miRNAs. • Parp1{sup −/−} embryonic stem cell-derived exosomes were examined for miRNA regulation. • Next-generation sequencing revealed upregulated and downregulated miRNAs. • Upregulated miRNAs were associated with five signaling pathways. • Gene ontology annotation for growth/stress-related signaling and cell communication. Poly (ADP-ribose) polymerase family, member 1 (Parp1) has pleiotropic and disparate functions in multiple cellular signaling pathways through post-translational protein modification. It contributes to the regulation of various cellular processes, including DNA damage repair, cell death, and cell differentiation, genetically or epigenetically. Meanwhile, the functions of Parp1 in intercellular signaling remain to be established. To examine the functions of Parp1 in intercellular signaling, we examined microRNA (miRNA) regulation in exosomes derived from Parp1-deficient (Parp1{sup −/−}) embryonic stem (ES) cells. The percentages of miRNAs among total RNAs, including small RNAs such as miRNAs, snRNAs, snoRNAs, tRNAs, exonic RNAs, and intronic RNAs, in Parp1{sup +/+} and Parp1{sup −/−} ES cell-derived exosomes were 8.2% and 3.5%, respectively. Overall, 329 distinct miRNAs exhibited ≥2-fold changes (118 upregulated; 211 downregulated). The upregulated miRNAs targeted 810 candidate genes, and the downregulated miRNAs targeted 716 candidate genes. Pathway analyses revealed that the upregulated miRNAs were significantly associated with five pathways including MAPK signaling cascades (p Parp1{sup −/−} ES cells. In quantitative analyses of miRNA expression, miR365-3p, let-7a-5p, miR196b-5p, miR203-3p, miR98-5p, and miR146a-5p were increased by ≥ 2-fold in Parp1{sup −/−} ES cell-derived exosomes. Gene ontology enrichment analyses revealed that the upregulated miRNAs were significantly annotated for growth and stress-related cell signaling and cell communication (p Parp1 deficiency in ES cells led to inhibition of cell–cell communication, possibly by intercellular signal transduction, suggesting that Parp1 functions extracellularly by regulating exosomal miRNAs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Benchmarking second and third-generation sequencing platforms for microbial metagenomics

Shotgun metagenomic sequencing is a common approach for studying the taxonomic diversity and metabolic potential of complex microbial communities. Current methods primarily use second generation short read sequencing, yet advances in third generation long read technologies provide opportunities to overcome some of the limitations of short read sequencing. Here, we compared seven platforms, encompassing second generation sequencers (Illumina HiSeq 300, MGI DNBSEQ-G400 and DNBSEQ-T7, ThermoFisher Ion GeneStudio S5 and Ion Proton P1) and third generation sequencers (Oxford Nanopore Technologies MinION R9 and Pacific Biosciences Sequel II). We constructed three uneven synthetic microbial communities composed of up to 87 genomic microbial strains DNAs per mock, spanning 29 bacterial and archaeal phyla, and representing the most complex and diverse synthetic communities used for sequencing technology comparisons. Our results demonstrate that third generation sequencing have advantages over second generation platforms in analyzing complex microbial communities, but require careful sequencing library preparation for optimal quantitative metagenomic analysis. Our sequencing data also provides a valuable resource for testing and benchmarking bioinformatics software for metagenomics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Identification of a new circulating recombinant form of human immunodeficiency virus type 1, CRF124_cpx involving subtypes A, G, H, and CRF27_cpx in Angola

Angola, located in Central Africa, has around 320,000 (270,000–380,000) people living with human immunodeficiency virus (HIV)/AIDS, equivalent to 1% of the country’s population at the end of 2021. A previous study conducted in 2012, using Angolan samples collected between 2008 and 2010 revealed a high prevalence of HIV-1 recombinants, around 42% of sequences, with 21% showing the same UH profile in partial pol region which were grouped into a monophyletic cluster with high bootstrap support. Thus, the objective of the present work was to obtain complete genomes of those sequences and characterize them, aiming at a description of a new circulating recombinant form (CRF). Whole blood from nine HIV-1 UH pol -infected individuals had their genomic DNA extracted, and nested PCR was used to amplify seven overlapping fragments targeting the full-length HIV-1 genome. The final classification was based on maximum likelihood trees, and recombination analyses were performed using a bootscan from the Simplot program. BLAST and Los Alamos Database inspections were used to search other similar H-like pol sequences. Complete genome amplification was possible for three samples, partial genomes were obtained for the other three, and only pol was available for the remaining three sequences. Bootscan analysis of the two whole-genome and three partial genome sequences retrieved from people living with HIV/AIDS (PLHIVA) without epidemiological linkage showed the same complex recombination profile involving HIV-1 subtypes A/G/H/CRF27_cpx, with a total of six recombinant breakpoints, aiming to classify a new HIV-1 CRF124_cpx. We found no other full-length HIV-1 genomes with the same mosaic profile; however, we identified 33 partial pol sequences, mainly sampled from Angola between 2001 to 2019, with the same H-like profile. Bayesian analysis of H and H-like pol sequences indicates that CRF124_cpx probably originated in Angola at mid-1970s, indicating that this CRF has been circulating in the country for a long time. In summary, our study describes a new CRF circulating principally in Angola and highlights the importance of continuing molecular surveillance studies, especially in countries with high molecular diversity of HIV.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Comprehensive High Throughput One‐Pot Analysis of Glycoproteins through Tool for Rapid Analysis of Glycopeptide by Permethylation (TRAP) Method.

Glycoproteins plays key roles in various cellular functions such as protein folding, protein targeting, receptor function, cell–cell recognition, protein stability, signaling, and other cellular functions. Currently a comprehensive and universal tool for the automated analysis of glycopeptides through their tandem mass spectrometry and data processing is a tedious process. We recently developed a method, Tool for Rapid Analysis of glycopeptide by Permethylation (TRAP), that simplifies glycoprotein analysis by integrating glycan sequencing and glycopeptide analysis in a single experiment. We employed this approach for the detailed analysis of glycosidic linkages, composition, glycan sequencing, glycosylation site mapping and sequencing of peptide moiety of glycopeptides of glycoproteins from both prokaryotes and eukaryotes. Here we are presenting an easy to use, robust, suitable for large as well as small sample sets and also cost‐effective glycoprotein analysis tool through TRAP approach. A permethylation protocol which can be automated and suitable for reliable permethylation of glycopeptides from the protease digest of glycoproteins was developed. Initial studies on glycopeptides from a model glycoprotein bovine RNAse B indicates that improved sensitivity and throughput of permethylation reactions can be achieved in conducting the reaction in microscale. This approach would enable convenient, rapid, and reliable high‐throughput screening of protein glycosylation and can be employed as a cost‐effective kit for glycosylated biomarker discoveries. In order to improvise the comprehensive glycosylation analysis and throughput, we have also developed an LC‐MS/MS platform for the analysis of glycopeptides processed through TRAP. We have conducted the analysis of standard glycoproteins such as bovine RNase B, human transferrin, and bovine fetuin and our preliminary results indicate that TRAP hyphened with LC‐MS/MS enables characterization of isomers of each glycoform at each glycosylation sites of complex glycoprotein mixtures. Support or Funding Information This research was supported by the National Institutes of Health (NIH) funded research, Grant R21GM122633; resource for biomedical glycomics, Grant P41GM103490; instrument grant, Grant S10OD018530; and U.S. Department of Energy grant, Grant DE‐SC0015662, to the Complex Carbohydrate Research Center. This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal .

Shajahan, Asif↗
DOE OSTI2019

Diverse tumour susceptibility in Collaborative Cross mice: identification of a new mouse model for human gastric tumourigenesis

Objective: The Collaborative Cross (CC) is a mouse population model with diverse and reproducible genetic backgrounds used to identify novel disease models and genes that contribute to human disease. Since spontaneous tumour susceptibility in CC mice remains unexplored, we assessed tumour incidence and spectrum. Design: We monitored 293 mice from 18 CC strains for tumour development. Genetic association analysis and RNA sequencing were used to identify susceptibility loci and candidate genes. We analysed genomes of patients with gastric cancer to evaluate the relevance of genes identified in the CC mouse model and measured the expression levels of ISG15 by immunohistochemical staining using a gastric adenocarcinoma tissue microarray. Association of gene expression with overall survival (OS) was assessed by Kaplan-Meier analysis. Results: CC mice displayed a wide range in the incidence and types of spontaneous tumours. More than 40% of CC036 mice developed gastric tumours within 1 year. Genetic association analysis identified Nfκb1 as a candidate susceptibility gene, while RNA sequencing analysis of non-tumour gastric tissues from CC036 mice showed significantly higher expression of inflammatory response genes. In human gastric cancers, the majority of human orthologues of the 166 mouse genes were preferentially altered by amplification or deletion and were significantly associated with OS. Higher expression of the CC036 inflammatory response gene signature is associated with poor OS. Finally, ISG15 protein is elevated in gastric adenocarcinomas and correlated with shortened patient OS. Conclusions: CC strains exhibit tremendous variation in tumour susceptibility, and we present CC036 as a spontaneous laboratory mouse model for studying human gastric tumourigenesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

A1 Signature pattern and phylogenetic analysis of full-length env genes in 20 hemophiliacs infected with Korean subclade of HIV-1 subtype B

Twenty patients with hemophilia (HP) were diagnosed with HIV-1, one to two years after exposure to domestic clotting factor IX (DCF), in 1989–90. Previous analysis of pol and vif genes confirmed that HIV-1 transmission to 20 HPs occurred through infusion of DCF. In this study, we determined full-length env gene sequences (about 2,550 bp) in 21 HPs, three plasma donors whose plasma were used in DCF production, and 45 local controls. The env gene from frozen stored sera obtained 1–3 years after diagnosis as well as from samples collected years after infection, was amplified by RT-PCR, and subjected to direct sequencing. Phylogenetic analysis revealed that all sequences were subtype B, with 109 sequences from 64 patients (20 HPs, 3 serum donors, 41 local controls, LCs) belonging to the Korean subclade (KSB) and 13 sequences from 5 patients (1HP infected outside Korea, 4 LCs) not in the KSB. Sequences of the env gene from donors O and P plus the 20 HPs comprised two subclusters within KSB. In addition, signature pattern analysis revealed signature nucleotides at 45 positions between the HPs and LCs (P < 0.05). Surprisingly, specific signature nucleotides positions were conserved 100% in clusters O and P only (at 4 and 1 positions, respectively) with none in LCs. Within both clusters, sequence identity was high. These results are consistent with our previous conclusion that these 20 HPs were infected with viruses in the DCF used for treatment. In addition, we found that there are 25 signature pattern residues in amino acids in env gene of KSB distinct from subtype B.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Unusual modifications of protein biomarkers expressed by plasmid, prophage, and bacterial host of pathogenic Escherichia coli identified using top‐down proteomic analysis

Rationale Pathogenic bacteria often carry prophage (bacterial viruses) and plasmids (small circular pieces of DNA) that may harbor toxin, antibacterial, and antibiotic resistance genes. Proteomic characterization of pathogenic bacteria should include the identification of host proteins and proteins produced by prophage and plasmid genomes. Methods Protein biomarkers of two strains of Shiga toxin–producingEscherichia coli(STEC) were identified using antibiotic induction, matrix‐assisted laser desorption/ionization tandem time‐of‐flight (MALDI‐TOF‐TOF) tandem mass spectrometry (MS/MS) with post‐source decay (PSD), top‐down proteomic (TDP) analysis, and plasmid sequencing. Alphafold2 was also used to compare predicted in silico structures of the identified proteins to prominent fragment ions generated using MS/MS‐PSD. Strain samples were also analyzed with and without chemical reduction treatment to detect the attachment of pendant groups bound by thioester or disulfide bonds. Results Shiga toxin was detected and/or identified in both STEC strains. For the first time, we also identified the osmotically inducible protein (OsmY) whose sequence unexpectedly had two forms: a full and a truncated sequence. The truncated OsmY terminates in the middle of an α‐helix as determined by Alphafold2. A plasmid‐encoded colicin immunity protein was also identified with and without attachment of an unidentified cysteine‐bound pendant group (~307 Da). Plasmid sequencing confirmed top‐down analysis and the identification of a promoter upstream of the immunity gene that is activated by antibiotic induction, that is, SOS box. Conclusions TDP analysis, coupled with other techniques (e.g., antibiotic induction, chemical reduction, plasmid sequencing, and in silico protein modeling), is a powerful tool to identify proteins (and their modifications), including prophage‐ and plasmid‐encoded proteins, produced by pathogenic microorganisms.

Biochemistry & Molecular Biology↗
DOE OSTI2019

Genetic Analysis of the Full-Length gag Gene from the Earliest Korean Subclade B of HIV-1: An Outbreak among Korean Hemophiliacs

We determined the earliest full-length HIV-1 gag gene sequences in 110 patients with HIV-1, including 20 hemophiliacs (HPs) and 90 local controls (LCs). The gag gene from stored sera was amplified using RT-PCR, and was subjected to direct sequencing. Phylogenetic analysis indicated that 94 and 16 sequences belonged to the Korean subclade of HIV-1 subtype B (KSB) and subtype B, respectively. A total of 12 signature pattern amino acids were found within the KSB, distinct from the worldwide consensus of subtype B. Within the KSB, the gag gene sequences from donors O and P and those from the 20 HPs comprised two subclusters. In particular, sequences from donor O strongly clustered with those of eight HPs. Moreover, signature pattern analysis indicated that 14 signature nucleotides were shared between the HPs and LCs within KSB (p < 0.01). Among the 14 nucleotides, positions 9 and 5 belonged to clusters O and P, respectively. In conclusion, signature pattern analysis for the gag gene revealed 12 signature pattern residues within the KSB and also confirmed the previous conclusion that the 20 HPs were infected with viruses due to incompletely inactivated clotting factor IX. This study is the first genetic analysis of the HIV-1 gag gene in Korea.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

MerCat: a versatile k-mer counter and diversity estimator for database-independent property analysis obtained from metagenomic and/or metatranscriptomic sequencing data

MerCat is a parallel, highly scalable and modular property software package for robust analysis of features in next-generation sequencing data. MerCat inputs include assembled contigs and raw sequence reads from any platform resulting in feature abundance counts tables. MerCat allows for direct analysis of data properties without reference sequence database dependency commonly used by search tools such as BLAST and/or DIAMOND for compositional analysis of whole community shotgun sequencing (e.g. metagenomes and metatranscriptomes).

metabolomics↗
DOE OSTI2023

National Prevalence of Caprine Prion Protein Genetic Variability at Codons 146, 211, and 222 in Goat Herds in the United States

Scrapie is a neurodegenerative disease that impacts sheep and goats, characterized by gradual and progressive changes in neurological function. Recent research shows that the scrapie incubation period is significantly influenced by specific variations in amino acids within the prion protein gene (PRNP). The objective of this study was to estimate the national prevalence of caprine PRNP genetic variability at codons 146, 211, and 222 in goat populations across the United States. A total of 3052 blood, ear tissue, and brain tissue samples were collected from goats from 50 states. The participating states were categorized into four Veterinary Service (VS) district regions. The samples underwent DNA extraction, and the PRNP variants corresponding to codons 146, 211, and 222 were amplified and sequenced. The analysis of PRNP variants, when compared to the PRNP reference sequence, revealed seven alleles in twelve genotypes. The homozygous 146NN, 211RR, and 222QQ alleles, which have been linked to an increased risk of scrapie, were found to be the most prevalent among all the goats. The heterozygous 222QK, 211RQ, 146SD, 146ND, and 146NS alleles and the homozygous 222KK, 146SS, and 146DD alleles, known to be associated with reduced scrapie susceptibility and a prolonged incubation period after experimental challenge, were found in 1.098% (222QK), 2.33% (211RQ), 0.58% (146SD), 3.13% (146ND), 20.68% (146NS), 0.005% (222KK), 3.31% (146SS), and 0.67% (146DD) of goats, respectively. The 222QK allele was found most frequently in goats tested from the east (VS District 1, 1.59%) and southwest (VS District 4, 1.08%) regions, whereas the 211RQ allele was found most often in goats tested from the Midwest (VS District 2, 8.03%) and east (VS District 1, 6.53%) regions. The 146NS allele was found most frequently in goats tested from the northwest (VS District 3, 29.02%) and southwest (VS District 4, 20.69%) regions. Our results showed that the prevalence of less susceptible genotypes at PRNP codon 146 may be sufficient to use genetic susceptibility testing in some herds. This may reduce the number of goats removed as part of a herd clean-up plan and may promote the selective breeding goats for less susceptible alleles in high-risk herds at the national level.

59 BASIC BIOLOGICAL SCIENCES↗