Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sequence Alignment”

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 217 records · Page 12

A gatekeeper chaperone complex directs translocator secretion during Type Three Secretion

Many Gram-negative bacteria use Type Three Secretion Systems (T3SS) to deliver effector proteins into host cells. These protein delivery machines are composed of cytosolic components that recognize substrates and generate the force needed for translocation, the secretion conduit, formed by a needle complex and associated membrane spanning basal body, and translocators that form the pore in the target cell. A defined order of secretion in which needle component proteins are secreted first, followed by translocators, and finally effectors, is necessary for this system to be effective. While the secreted effectors vary significantly between organisms, the ~20 individual protein components that form the T3SS are conserved in many pathogenic bacteria. One such conserved protein, referred to as either a plug or gatekeeper, is necessary to prevent unregulated effector release and to allow efficient translocator secretion. The mechanism by which translocator secretion is promoted while effector release is inhibited by gatekeepers is unknown. We present the structure of the Chlamydial gatekeeper, CopN, bound to a translocator-specific chaperone. The structure identifies a previously unknown interface between gatekeepers and translocator chaperones and reveals that in the gatekeeper-chaperone complex the canonical translocator-binding groove is free to bind translocators. Thus, structure-based mutagenesis of the homologous complex in Shigella reveals that the gatekeeper-chaperone-translocator complex is essential for translocator secretion and for the ordered secretion of translocators prior to effectors.

59 BASIC BIOLOGICAL SCIENCES↗

Chaperone-tip adhesin complex is vital for synergistic activation of CFA/I fimbriae biogenesis

Colonization factor CFA/I defines the major adhesive fimbriae of enterotoxigenic Escherichia coli and mediates bacterial attachment to host intestinal epithelial cells. The CFA/I fimbria consists of a tip-localized minor adhesive subunit, CfaE, and thousands of copies of the major subunit CfaB polymerized into an ordered helical rod. Biosynthesis of CFA/I fimbriae requires the assistance of the periplasmic chaperone CfaA and outer membrane usher CfaC. Although the CfaE subunit is proposed to initiate the assembly of CFA/I fimbriae, how it performs this function remains elusive. Here, we report the establishment of an in vitro assay for CFA/I fimbria assembly and show that stabilized CfaA-CfaB and CfaA-CfaE binary complexes together with CfaC are sufficient to drive fimbria formation. The presence of both CfaA-CfaE and CfaC accelerates fimbria formation, while the absence of either component leads to linearized CfaB polymers in vitro. We further report the crystal structure of the stabilized CfaA-CfaE complex, revealing features unique for biogenesis of Class 5 fimbriae.

59 BASIC BIOLOGICAL SCIENCES↗

Neutralization profiles of HIV-1 viruses from the VRC01 Antibody Mediated Prevention (AMP) trials

The VRC01 Antibody Mediated Prevention (AMP) efficacy trials conducted between 2016 and 2020 showed for the first time that passively administered broadly neutralizing antibodies (bnAbs) could prevent HIV-1 acquisition against bnAb-sensitive viruses. HIV-1 viruses isolated from AMP participants who acquired infection during the study in the sub-Saharan African (HVTN 703/HPTN 081) and the Americas/European (HVTN 704/HPTN 085) trials represent a panel of currently circulating strains of HIV-1 and offer a unique opportunity to investigate the sensitivity of the virus to broadly neutralizing antibodies (bnAbs) being considered for clinical development. Pseudoviruses were constructed using envelope sequences from 218 individuals. The majority of viruses identified were clade B and C; with clades A, D, F and G and recombinants AC and BF detected at lower frequencies. We tested eight bnAbs in clinical development (VRC01, VRC07-523LS, 3BNC117, CAP256.25, PGDM1400, PGT121, 10–1074 and 10E8v4) for neutralization against all AMP placebo viruses (n = 76). Compared to older clade C viruses (1998–2010), the HVTN703/HPTN081 clade C viruses showed increased resistance to VRC07-523LS and CAP256.25. At a concentration of 1μg/ml (IC80), predictive modeling identified the triple combination of V3/V2-glycan/CD4bs-targeting bnAbs (10-1074/PGDM1400/VRC07-523LS) as the best against clade C viruses and a combination of MPER/V3/CD4bs-targeting bnAbs (10E8v4/10-1074/VRC07-523LS) as the best against clade B viruses, due to low coverage of V2-glycan directed bnAbs against clade B viruses. Overall, the AMP placebo viruses represent a valuable resource for defining the sensitivity of contemporaneous circulating viral strains to bnAbs and highlight the need to update reference panels regularly. Our data also suggests that combining bnAbs in passive immunization trials would improve coverage of global viruses.

60 APPLIED LIFE SCIENCES↗

Whole-Genome and Chromosome Evolution Associated with Host Adaptation and Speciation of the Wheat Pathogen Mycosphaerella graminicola

The fungus Mycosphaerella graminicola has been a pathogen of wheat since host domestication 10,000–12,000 years ago in the Fertile Crescent. The wheat-infecting lineage emerged from closely related Mycosphaerella pathogens infecting wild grasses. We use a comparative genomics approach to assess how the process of host specialization affected the genome structure of M. graminicola since divergence from the closest known progenitor species named M. graminicola S1. The genome of S1 was obtained by Illumina sequencing resulting in a 35 Mb draft genome sequence of 32X. Assembled contigs were aligned to the previously sequenced M. graminicola genome. The alignment covered .90% of the non-repetitive portion of the M. graminicola genome with an average divergence of 7%. The sequenced M. graminicola strain is known to harbor thirteen essential chromosomes plus eight dispensable chromosomes. We found evidence that structural rearrangements significantly affected the dispensable chromosomes while the essential chromosomes were syntenic. At the nucleotide level, the essential and dispensable chromosomes have evolved differently. The average synonymous substitution rate in dispensable chromosomes is considerably lower than in essential chromosomes, whereas the average non-synonymous substitution rate is three times higher. Differences in molecular evolution can be related to different transmission and recombination patterns, as well as to differences in effective population sizes of essential and dispensable chromosomes. In order to identify genes potentially involved in host specialization or speciation, we calculated ratios of synonymous and non-synonymous substitution rates in the .9,500 aligned protein coding genes. The genes are generally under strong purifying selection. We identified 43 candidate genes showing evidence of positive selection, one encoding a potential pathogen effector protein. We conclude that divergence of these pathogens was accompanied by structural rearrangements in the small dispensable chromosomes, while footprints of positive selection were present in only a small number of protein coding genes.

59 BASIC BIOLOGICAL SCIENCES↗

OrthoPhylo

This software builds on PHAME developed at LANL to generate phylogenetic trees of bacterial whole genome sequences. Where PHAME uses whole genome alignments to generate informative sites to base tress on, PHAME-OuS annotates bacterial genes, identifies orthologous sequences, aligns related proteins, uses those alignments to inform transcript alignments, then builds trees with several methods. The first is a conventional gene concatenation and ML tree estemation method. The second attempts to reconcile gene tree with a unified species tree using quartets (ASTRAL). Both methods allow filtering of gene lists on number of species represented, length, and gappiness in order to tune noise-to-signal for tree estimation

Middlebrook, Earl↗

Alignment of RNA molecules: Binding energy and statistical properties of random sequences

A new statistical approach to the problem of pairwise alignment of RNA sequences is proposed. The problem is analyzed for a pair of interacting polymers forming an RNA-like hierarchical cloverleaf structures. An alignment is characterized by the numbers of matches, mismatches, and gaps. A weight function is assigned to each alignment; this function is interpreted as a free energy taking into account both direct monomer-monomer interactions and a combinatorial contribution due to formation of various cloverleaf secondary structures. The binding free energy is determined for a pair of RNA molecules. Statistical properties are discussed, including fluctuations of the binding energy between a pair of RNA molecules and loop length distribution in a complex. Based on an analysis of the free energy per nucleotide pair complexes of random RNAs as a function of the number of nucleotide types c, a hypothesis is put forward about the exclusivity of the alphabet c = 4 used by nature.

60 APPLIED LIFE SCIENCES↗

BBMap: A Fast, Accurate, Splice-Aware Aligner

Alignment of reads is one of the primary computational tasks in bioinformatics. Of paramount importance to resequencing, alignment is also crucial to other areas - quality control, scaffolding, string-graph assembly, homology detection, assembly evaluation, error-correction, expression quantification, and even as a tool to evaluate other tools. An optimal aligner would greatly improve virtually any sequencing process, but optimal alignment is prohibitively expensive for gigabases of data. Here, we will present BBMap [1], a fast splice-aware aligner for short and long reads. We will demonstrate that BBMap has superior speed, sensitivity, and specificity to alternative high-throughput aligners bowtie2 [2], bwa [3], smalt, [4] GSNAP [5], and BLASR [6].

09 BIOMASS FUELS↗

Alignment of Short Reads: A Crucial Step for Application of Next-Generation Sequencing Data in Precision Medicine

Precision medicine or personalized medicine has been proposed as a modernized and promising medical strategy. Genetic variants of patients are the key information for implementation of precision medicine. Next-generation sequencing (NGS) is an emerging technology for deciphering genetic variants. Alignment of raw reads to a reference genome is one of the key steps in NGS data analysis. Many algorithms have been developed for alignment of short read sequences since 2008. Users have to make a decision on which alignment algorithm to use in their studies. Selection of the right alignment algorithm determines not only the alignment algorithm but also the set of suitable parameters to be used by the algorithm. Understanding these algorithms helps in selecting the appropriate alignment algorithm for different applications in precision medicine. Here, we review current available algorithms and their major strategies such as seed-and-extend and q-gram filter. We also discuss the challenges in current alignment algorithms, including alignment in multiple repeated regions, long reads alignment and alignment facilitated with known genetic variants.

59 BASIC BIOLOGICAL SCIENCES↗

Methods of resolving artifacts in Hadamard-transformed data

A method of validating data produced from a multiplexing process on an analytical instrument is disclosed. In one embodiment, the method includes using a pseudorandom sequence to encode a multiplexed segment of data; applying Hadamard transform to generate a demultiplexed segment of the data; aligning the pseudorandom sequence to the multiplexed data; and calculating a score for at least one positive value in the demultiplexed segment to find a valid demultiplexed value.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Developmental Pathway of the MPER-Directed HIV-1-Neutralizing Antibody 10E8

Antibody 10E8 targets the membrane-proximal external region (MPER) of HIV-1 gp41, neutralizes >97% of HIV-1 isolates, and lacks the auto-reactivity often associated with MPER-directed antibodies. The developmental pathway of 10E8 might therefore serve as a promising template for vaccine design, but samples from time-of-infection—often used to infer the B cell record—are unavailable. In this study, we used crystallography, next-generation sequencing (NGS), and functional assessments to infer the 10E8 developmental pathway from a single time point. Mutational analysis indicated somatic hypermutation of the 2nd-heavy chain-complementarity determining region (CDR H2) to be critical for neutralization, and structures of 10E8 variants with V-gene regions reverted to genomic origin for heavy-and-light chains or heavy chain-only showed structural differences >2 Å relative to mature 10E8 in the CDR H2 and H3. To understand these developmental changes, we used bioinformatic sieving, maximum likelihood, and parsimony analyses of immunoglobulin transcripts to identify 10E8-lineage members, to infer the 10E8-unmutated common ancestor (UCA), and to calculate 10E8-developmental intermediates. We were assisted in this analysis by the preservation of a critical D-gene segment, which was unmutated in most 10E8-lineage sequences. UCA and early intermediates weakly bound a 26-residue-MPER peptide, whereas HIV-1 neutralization and epitope recognition in liposomes were only observed with late intermediates. Antibody 10E8 thus develops from a UCA with weak MPER affinity and substantial differences in CDR H2 and H3 from the mature 10E8; only after extensive somatic hypermutation do 10E8-lineage members gain recognition in the context of membrane and HIV-1 neutralization.

60 APPLIED LIFE SCIENCES↗

Evolutionary genomics suggests that CheV is an additional adaptor for accommodating specific chemoreceptors within the chemotaxis signaling complex

Escherichia coli and Salmonella enterica are models for many experiments in molecular biology including chemotaxis, and most of the results obtained with one organism have been generalized to another. While most components of the chemotaxis pathway are strongly conserved between the two species, Salmonella genomes contain some chemoreceptors and an additional protein, CheV, that are not found in E. coli. The role of CheV was examined in distantly related species Bacillus subtilis and Helicobacter pylori , but its role in bacterial chemotaxis is still not well understood. We tested a hypothesis that in enterobacteria CheV functions as an additional adaptor linking the CheA kinase to certain types of chemoreceptors that cannot be effectively accommodated by the universal adaptor CheW. Phylogenetic profiling, genomic context and comparative protein sequence analyses suggested that CheV interacts with specific domains of CheA and chemoreceptors from an orthologous group exemplified by the Salmonella McpC protein. Structural consideration of the conservation patterns suggests that CheV and CheW share the same binding spot on the chemoreceptor structure, but have some affinity bias towards chemoreceptors from different orthologous groups. Finally, published experimental results and data newly obtained via comparative genomics support the idea that CheV functions as a "phosphate sink" possibly to off-set the over-stimulation of the kinase by certain types of chemoreceptors. Altogether, our results strongly suggest that CheV is an additional adaptor for accommodating specific chemoreceptors within the chemotaxis signaling complex.

59 BASIC BIOLOGICAL SCIENCES↗

MicroPheno: predicting environments and host phenotypes from 16S rRNA gene sequencing using a k-mer based representation of shallow sub-samples

Microbial communities play important roles in the function and maintenance of various biosystems, ranging from the human body to the environment. A major challenge in microbiome research is the classification of microbial communities of different environments or host phenotypes. The most common and cost-effective approach for such studies to date is 16S rRNA gene sequencing. Recent falls in sequencing costs have increased the demand for simple, efficient and accurate methods for rapid detection or diagnosis with proved applications in medicine, agriculture and forensic science. We describe a reference- and alignment-free approach for predicting environments and host phenotypes from 16S rRNA gene sequencing based on k-mer representations that benefits from a bootstrapping framework for investigating the sufficiency of shallow sub-samples. Deep learning methods as well as classical approaches were explored for predicting environments and host phenotypes. A k-mer distribution of shallow sub-samples outperformed Operational Taxonomic Unit (OTU) features in the tasks of body-site identification and Crohn’s disease prediction. Aside from being more accurate, using k-mer features in shallow sub-samples allows (i) skipping computationally costly sequence alignments required in OTU-picking and (ii) provided a proof of concept for the sufficiency of shallow and short-length 16S rRNA sequencing for phenotype prediction. In addition, k-mer features predicted representative 16S rRNA gene sequences of 18 ecological environments, and 5 organismal environments with high macro-F1 scores of 0.88 and 0.87. For large datasets, deep learning outperformed classical methods such as Random Forest and Support Vector Machine.

59 BASIC BIOLOGICAL SCIENCES↗

Impedance Characterization of Utility-Scale Renewable Energy and Storage Systems

This paper presents the impedance characterization of utility-scale wind turbines and PV/storage inverters using a multimegawatt grid simulator-based grid integration testbed. The paper describes the testbed components including a 7 MVA grid simulator, 5 MW dynamometer, medium-voltage measurement system, and a real-time control system used for automated impedance characterization of wind turbines and inverters. Different aspects of impedance measurement are discussed, including the selection of the magnitude and type of perturbation, measurement of frequency cross-coupling between sequence impedance responses, and the alignment of the reference frame for sequence impedance measurement. The paper also demonstrates high-fidelity validation of EMT simulation models from turbine manufacturers and characterization of grid-forming inverters using impedance measurements.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Impedance Characterization of Utility-Scale Renewable Energy and Storage Systems: Preprint

This paper presents the impedance characterization of utility-scale wind turbines and PV/storage inverters using a multimegawatt grid simulator-based grid integration testbed. The paper describes testbed components including a 7 MVA grid simulator, 5 MW dynamometer, medium-voltage measurement system, and real-time control system used for automated impedance characterization of wind turbines and inverters. Different aspects of impedance measurement are discussed, including the selection of the magnitude and type of perturbation, measurement of frequency cross-coupling between sequence impedance responses, and the alignment of the reference frame for sequence impedance measurement. The paper also demonstrates high-fidelity validation of EMT simulation models and characterization of grid-forming inverters using impedance measurements.

17 WIND ENERGY↗

Impedance Measurement of Wind Turbines using a Multimegawatt Grid Simulator: Preprint

This paper presents the impedance measurements of a 4 MW wind turbine conducted using a multimegawatt grid simulator at the National Renewable Energy Laboratory in Colorado, United States. This paper describes different systems, including a 7 MVA grid simulator, a 5 MW dynamometer, a medium-voltage measurement system, and a real-time control system used for the automated impedance measurement of a utility-scale wind turbine. Different aspects of impedance measurement are discussed, including the selection of the magnitude of injected perturbation, measurement of frequency cross-coupling between the sequence impedance responses, and the alignment of the reference frame for sequence impedance measurement. The paper also demonstrates high-fidelity validation of an EMT simulation model of the wind turbine using impedance measurements.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The Release 6 reference sequence of the Drosophila melanogaster genome

Drosophila melanogaster plays an important role in molecular, genetic, and genomic studies of heredity, development, metabolism, behavior, and human disease. The initial reference genome sequence reported more than a decade ago had a profound impact on progress in Drosophila research, and improving the accuracy and completeness of this sequence continues to be important to further progress. We previously described improvement of the 117-Mb sequence in the euchromatic portion of the genome and 21 Mb in the heterochromatic portion, using a whole-genome shotgun assembly, BAC physical mapping, and clone-based finishing. Here, we report an improved reference sequence of the single-copy and middle-repetitive regions of the genome, produced using cytogenetic mapping to mitotic and polytene chromosomes, clone-based finishing and BAC fingerprint verification, ordering of scaffolds by alignment to cDNA sequences, incorporation of other map and sequence data, and validation by whole-genome optical restriction mapping. These data substantially improve the accuracy and completeness of the reference sequence and the order and orientation of sequence scaffolds into chromosome arm assemblies. Representation of the Y chromosome and other heterochromatic regions is particularly improved. The new 143.9-Mb reference sequence, designated Release 6, effectively exhausts clone-based technologies for mapping and sequencing. Highly repeat-rich regions, including large satellite blocks and functional elements such as the ribosomal RNA genes and the centromeres, are largely inaccessible to current sequencing and assembly methods and remain poorly represented. In conclusion, further significant improvements will require sequencing technologies that do not depend on molecular cloning and that produce very long reads.

59 BASIC BIOLOGICAL SCIENCES↗

Automated registration of multispectral MR vessel wall images of the carotid artery

Purpose: Atherosclerosis is the primary cause of heart disease and stroke. The detailed assessment of atherosclerosis of the carotid artery requires high resolution imaging of the vessel wall using multiple MR sequences with different contrast weightings. These images allow manual or automated classification of plaque components inside the vessel wall. Automated classification requires all sequences to be in alignment, which is hampered by patient motion. In clinical practice, correction of this motion is performed manually. Previous studies applied automated image registration to correct for motion using only nondeformable transformation models and did not perform a detailed quantitative validation. The purpose of this study is to develop an automated accurate 3D registration method, and to extensively validate this method on a large set of patient data. In addition, the authors quantified patient motion during scanning to investigate the need for correction. Methods: MR imaging studies (1.5T, dedicated carotid surface coil, Philips) from 55 TIA/stroke patients with ipsilateral <70% carotid artery stenosis were randomly selected from a larger cohort. Five MR pulse sequences were acquired around the carotid bifurcation, each containing nine transverse slices: T1-weighted turbo field echo, time of flight, T2-weighted turbo spin-echo, and pre- and postcontrast T1-weighted turbo spin-echo images (T1W TSE). The images were manually segmented by delineating the lumen contour in each vessel wall sequence and were manually aligned by applying throughplane and inplane translations to the images. To find the optimal automatic image registration method, different masks, choice of the fixed image, different types of the mutual information image similarity metric, and transformation models including 3D deformable transformation models, were evaluated. Evaluation of the automatic registration results was performed by comparing the lumen segmentations of the fixed image and moving image after registration. Results: The average required manual translation per image slice was 1.33 mm. Translations were larger as the patient was longer inside the scanner. Manual alignment took 187.5 s per patient resulting in a mean surface distance of 0.271 ± 0.127 mm. After minimal user interaction to generate the mask in the fixed image, the remaining sequences are automatically registered with a computation time of 52.0 s per patient. The optimal registration strategy used a circular mask with a diameter of 10 mm, a 3D B-spline transformation model with a control point spacing of 15 mm, mutual information as image similarity metric, and the precontrast T1W TSE as fixed image. A mean surface distance of 0.288 ± 0.128 mm was obtained with these settings, which is very close to the accuracy of the manual alignment procedure. The exact registration parameters and software were made publicly available. Conclusions: An automated registration method was developed and optimized, only needing two mouse clicks to mark the start and end point of the artery. Validation on a large group of patients showed that automated image registration has similar accuracy as the manual alignment procedure, substantially reducing the amount of user interactions needed, and is multiple times faster. In conclusion, the authors believe that the proposed automated method can replace the current manual procedure, thereby reducing the time to analyze the images.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Genome-Wide Transcription Factor DNA Binding Sites and Gene Regulatory Networks in Clostridium thermocellum

Clostridium thermocellum is a thermophilic bacterium recognized for its natural ability to effectively deconstruct cellulosic biomass. While there is a large body of studies on the genetic engineering of this bacterium and its physiology to-date, there is limited knowledge in the transcriptional regulation in this organism and thermophilic bacteria in general. The study herein is the first report of a large-scale application of DNA-affinity purification sequencing (DAP-seq) to transcription factors (TFs) from a bacterium. We applied DAP-seq to > 90 TFs in C. thermocellum and detected genome-wide binding sites for 11 of them. We then compiled and aligned DNA binding sequences from these TFs to deduce the primary DNA-binding sequence motifs for each TF. These binding motifs are further validated with electrophoretic mobility shift assay (EMSA) and are used to identify individual TFs’ regulatory targets in C. thermocellum . Our results led to the discovery of novel, uncharacterized TFs as well as homologues of previously studied TFs including RexA-, LexA-, and LacI-type TFs. We then used these data to reconstruct gene regulatory networks for the 11 TFs individually, which resulted in a global network encompassing the TFs with some interconnections. As gene regulation governs and constrains how bacteria behave, our findings shed light on the roles of TFs delineated by their regulons, and potentially provides a means to enable rational, advanced genetic engineering of C. thermocellum and other organisms alike toward a desired phenotype.

59 BASIC BIOLOGICAL SCIENCES↗