Search NASA⌕ Search

Engineering topics

Williams, Kelly P.

Publications and source records attributed to Williams, Kelly P..

At least 19 records

EMMA: a new method for computing multiple sequence alignments given a constraint subset alignment

Abstract Background Adding sequences into an existing (possibly user-provided) alignment has multiple applications, including updating a large alignment with new data, adding sequences into a constraint alignment constructed using biological knowledge, or computing alignments in the presence of sequence length heterogeneity. Although this is a natural problem, only a few tools have been developed to use this information with high fidelity. Results We present EMMA (Extending Multiple alignments using MAFFT--add) for the problem of adding a set of unaligned sequences into a multiple sequence alignment (i.e., a constraint alignment). EMMA builds on MAFFT--add, which is also designed to add sequences into a given constraint alignment. EMMA improves on MAFFT--add methods by using a divide-and-conquer framework to scale its most accurate version, MAFFT-linsi--add, to constraint alignments with many sequences. We show that EMMA has an accuracy advantage over other techniques for adding sequences into alignments under many realistic conditions and can scale to large datasets with high accuracy (hundreds of thousands of sequences). EMMA is available at https://github.com/c5shen/EMMA . Conclusions EMMA is a new tool that provides high accuracy and scalability for adding sequences into an existing alignment.

Shen, Chengze↗

Biotic countermeasures that rescue Nannochloropsis gaditana from a Bacillus safensis infection

The natural assemblage of a symbiotic bacterial microbiome (bacteriome) with microalgae in marine ecosystems is now being investigated as a means to increase algal productivity for industry. When algae are grown in open pond settings, biological contamination causes an estimated 30% loss of the algal crop. Therefore, new crop protection strategies that do not disrupt the native algal bacteriome are needed to produce reliable, high-yield algal biomass. Bacteriophages offer an unexplored solution to treat bacterial pathogenicity in algal cultures because they can eliminate a single species without affecting the bacteriome. To address this, we identified a highly virulent pathogen of the microalga Nannochloropsis gaditana, the bacterium Bacillus safensis, and demonstrated rescue of the microalgae from the pathogen using phage. 16S rRNA amplicon sequencing showed that phage treatment did not alter the composition of the bacteriome. It is widely suspected that the algal bacteriome could play a protective role against bacterial pathogens. To test this, we compared the susceptibility of a bacteriome-attenuated N. gaditana culture challenged with B. safensis to a N. gaditana culture carrying a growth-promoting bacteriome. We showed that the loss of the bacteriome increased the susceptibility of N. gaditana to the pathogen. Transplanting the microalgal bacteriome to the bacteriome-attenuated culture reconstituted the protective effect of the bacteriome. Finally, the success of phage treatment was dependent on the presence of beneficial bacteriome. This study introduces two synergistic countermeasures against bacterial pathogenicity in algal cultures and a tractable model for studying interactions between microalgae, phages, pathogens, and the algae microbiome.

60 APPLIED LIFE SCIENCES↗

Whole‐Cell MALDI‐ToF MS Coupled with Untargeted Metabolomics Facilitates Investigations of Microbial Chemical Interactions

Abstract The emergence of drug‐resistant pathogens necessitates the development of new countermeasures. In this regard, the introduction of probiotics to directly attack or competitively exclude pathogens presents a useful strategy. Application of this approach requires an understanding of how a probiotic and its target pathogen interact. A key means of probiotic‐pathogen interaction involves the production of small molecules called natural products (NPs). Here, we report the use of whole‐cell matrix‐assisted laser desorption/ionization time‐of‐flight (MALDI‐ToF) mass spectrometry to characterize NP production by candidate probiotics (mouse airway microbiome isolates) when co‐cultured with the respiratory pathogen Burkholderia . We found that a Bacillus velezensis strain inhibits growth of and elicits NP production by Burkholderia thailandensis . Dereplication of known NPs detected in the metabolome of this B. velezensis strain suggests that a previously unannotated bioactive compound is involved. Thus, we present the use of whole‐cell MALDI as a broadly applicable method for screening the NP composition of microbial co‐cultures; this can be combined with other ‐omics methods to characterize probiotic‐pathogen and other microbe‐microbe interactions.

Aiosa, Nicole↗

Improved Mobilome Delineation in Fragmented Genomes

The mobilome of a microbe, i.e., its set of mobile elements, has major effects on its ecology, and is important to delineate properly in each genome. This becomes more challenging for incomplete genomes, and even more so for metagenome-assembled genomes (MAGs), where misbinning of scaffolds and other losses can occur. Genomic islands (GIs), which integrate into the host chromosome, are a major component of the mobilome. Our GI-detection software TIGER, unique in its precise mapping of GI termini, was applied to 74,561 genomes from 2,473 microbial species, each species containing at least one MAG and one isolate genome. A species-normalized deficit of ∼1.6 GIs/genome was measured for MAGs relative to isolates. To test whether this undercount was due to the higher fragmentation of MAG genomes, TIGER was updated to enable detection of split GIs whose termini are on separate scaffolds or that wrap around the origin of a circular replicon. This doubled GI yields, and the new split GIs matched the quality of single-scaffold GIs, except that highly fragmented GIs may lack central portions. Cross-scaffold search is an important upgrade to GI detection as fragmented genomes increasingly dominate public databases. TIGER2 better captures MAG microdiversity, recovering niche-defining GIs and supporting microbiome research aims such as virus-host linking and ecological assessment.

54 ENVIRONMENTAL SCIENCES↗

Data Science and Machine Learning for Genome Security

This report describes research conducted to use data science and machine learning methods to distinguish targeted genome editing versus natural mutation and sequencer machine noise. Genome editing capabilities have been around for more than 20 years, and the efficiencies of these techniques has improved dramatically in the last 5+ years, notably with the rise of CRISPR-Cas technology. Whether or not a specific genome has been the target of an edit is concern for U.S. national security. The research detailed in this report provides first steps to address this concern. A large amount of data is necessary in our research, thus we invested considerable time collecting and processing it. We use an ensemble of decision tree and deep neural network machine learning methods as well as anomaly detection to detect genome edits given either whole exome or genome DNA reads. The edit detection results we obtained with our algorithms tested against samples held out during training of our methods are significantly better than random guessing, achieving high F1 and recall scores as well as with precision overall.

59 BASIC BIOLOGICAL SCIENCES↗

Expanding the Scope of Genomic Security: Targeted Genome Editing within Microbiomes through Designer Bacteriophage Vectors

The ability to engineer the genome of a bacterial strain, not as an isolate, but while present among other microbes in a microbiome, would open new technological possibilities in the areas of medicine, energy and biomanufacturing. Our approach is to develop sets of phages (bacterial viruses) active on the target strain and themselves engineered to act not as killers but as vectors for gene delivery. This approach is rooted in our bioinformatic tools that map prophages accurately within bacterial genomes. We present new bioinformatic results in cross-contig search, design of phage genome assemblies, satellites that embed within prophages, alignment of large numbers of biological sequences, and improvement of reference databases for prophage discovery. We targeted a Pseudomonas putida strain within a lignin-degrading microbiome, but were unable to obtain active phages, and turned toward a defined microbiome of the mouse gut.

59 BASIC BIOLOGICAL SCIENCES↗

CasCollect: targeted assembly of CRISPR-associated operons from high-throughput sequencing data

Abstract CRISPR arrays and CRISPR-associated (Cas) proteins comprise a widespread adaptive immune system in bacteria and archaea. These systems function as a defense against exogenous parasitic mobile genetic elements that include bacteriophages, plasmids and foreign nucleic acids. With the continuous spread of antibiotic resistance, knowledge of pathogen susceptibility to bacteriophage therapy is becoming more critical. Additionally, gene-editing applications would benefit from the discovery of new cas genes with favorable properties. While next-generation sequencing has produced staggering quantities of data, transitioning from raw sequencing reads to the identification of CRISPR/Cas systems has remained challenging. This is especially true for metagenomic data, which has the highest potential for identifying novel cas genes. We report a comprehensive computational pipeline, CasCollect, for the targeted assembly and annotation of cas genes and CRISPR arrays—even isolated arrays—from raw sequencing reads. Benchmarking our targeted assembly pipeline demonstrates significantly improved timing by almost two orders of magnitude compared with conventional assembly and annotation, while retaining the ability to detect CRISPR arrays and cas genes. CasCollect is a highly versatile pipeline and can be used for targeted assembly of any specialty gene set, reconfigurable for user provided Hidden Markov Models and/or reference nucleotide sequences.

Podlevsky, Joshua D.↗

Computational Basis for On-Demand Production of Diversified Therapeutic Phage Cocktails

New therapies are necessary to combat increasingly antibiotic-resistant bacterial pathogens. We have developed a technology platform of computational, molecular biology, and microbiology tools which together enable on-demand production of phages that target virtually any given bacterial isolate. Two complementary computational tools that identify and precisely map prophages and other integrative genetic elements in bacterial genomes are used to identify prophage-laden bacteria that are close relatives of the target strain. Phage genomes are engineered to disable lysogeny, through use of long amplicon PCR and Gibson assembly. Finally, the engineered phage genomes are introduced into host bacteria for phage production. As an initial demonstration, we used this approach to produce a phage cocktail against the opportunistic pathogen Pseudomonas aeruginosa PAO1. Two prophageladen P. aeruginosa strains closely related to PAO1 were identified, ATCC 39324 and ATCC 27853. Deep sequencing revealed that mitomycin C treatment of these strains induced seven phages that grow on P. aeruginosa PAO1. The most diverse five phages were engineered for nonlysogeny by deleting the integrase gene (int), which is readily identifiable and typically conveniently located at one end of the prophage. The Δint phages, individually and in cocktails, killed P. aeruginosa PAO1 in liquid culture as well as in a waxworm (Galleria mellonella) model of infection.

59 BASIC BIOLOGICAL SCIENCES↗

Genome Sequences of Burkholderia thailandensis Strains E421, E426, and DW503

ABSTRACT We present the draft genome sequences of three Burkholderia thailandensis strains, E421, E426, and DW503. E421 consists of 90 contigs of 6,639,935 bp and 67.73% GC content. E426 consists of 106 contigs of 6,587,853 bp and 67.73% GC content. DW503 consists of 102 contigs of 6,458,767 bp and 67.64% GC content.

59 BASIC BIOLOGICAL SCIENCES↗

New candidates for regulated gene integrity revealed through precise mapping of integrative genetic elements

Abstract Integrative genetic elements (IGEs) are mobile multigene DNA units that integrate into and excise from host bacterial genomes. Each IGE usually targets a specific site within a conserved host gene, integrating in a manner that preserves target gene function. However, a small number of bacterial genes are known to be inactivated upon IGE integration and reactivated upon excision, regulating phenotypes of virulence, mutation rate, and terminal differentiation in multicellular bacteria. The list of regulated gene integrity (RGI) cases has been slow-growing because IGEs have been challenging to precisely and comprehensively locate in genomes. We present software (TIGER) that maps IGEs with unprecedented precision and without attB site bias. TIGER uses a comparative genomic, ping-pong BLAST approach, based on the principle that the IGE integration module (i.e. its int-attP region) is cohesive. The resultant IGEs from 2168 genomes, along with integrase phylogenetic analysis and gene inactivation tests, revealed 19 new cases of genes whose integrity is regulated by IGEs (including dut, eccCa1, gntT, hrpB, merA, ompN, prkA, tqsA, traG, yifB, yfaT and ynfE), as well as recovering previously known cases (in sigK, spsM, comK, mlrA and hlb genes). It also recovered known clades of site-promiscuous integrases and identified possible new ones.

59 BASIC BIOLOGICAL SCIENCES↗

Precise mapping of new group I introns in tRNA genes

Bacterial tRNA have been found interrupted at various positions in the anticodon loop by group I introns, in four types. The primary bioinformatic tool for group I intron discovery is a covariance model that can identify conserved features in the catalytic core and can sometimes identify the typical uridine residue at the -1 position, preceding the 5-prime splice site, but cannot identify the typical guanidine residue at the omega position, preceding the 3-prime splice site, to achieve precise mapping. One approach to complete the automation of group I intron mapping is to focus instead on the exons, which is enabled by the regularity of tRNAs. We develop a software module, within a larger package (tFind) aimed at mapping bacterial tRNA and tmRNA genes precisely, that expands this list of four known classes of intron-interrupted tRNAs to 21 cases. A new covariance model for these introns is presented. The wobble base pair formed by the -1 uridine is considered a determinant of the 5-prime splice site, yet one reasonably large new type bears a cytidine nucleotide at that position.

Williams, Kelly P.↗

RNAcentral: a hub of information for non-coding RNA sequences

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences, collating information on ncRNA sequences of all types from a broad range of organisms. We have recently added a new genome mapping pipeline that identifies genomic locations for ncRNA sequences in 296 species. We have also added several new types of functional annotations, such as tRNA secondary structures, Gene Ontology annotations, and miRNA-target interactions. A new quality control mechanism based on Rfam family assignments identifies potential contamination, incomplete sequences, and more. The RNAcentral database has become a vital component of many workflows in the RNA community, serving as both the primary source of sequence data for academic and commercial groups, as well as a source of stable accessions for the annotation of genomic and functional features. These examples are facilitated by an improved RNAcentral web interface, which features an updated genome browser, a new sequence feature viewer, and improved text search functionality.

Sweeney, Blake A.↗

Extraction and biomolecular analysis of dermal interstitial fluid collected with hollow microneedles

Dermal interstitial fluid (ISF) is an underutilized information-rich biofluid potentially useful in health status monitoring applications whose contents remain challenging to characterize. Here, we present a facile microneedle approach for dermal ISF extraction with minimal pain and no blistering for human subjects and rats. Extracted ISF volumes were sufficient for determining transcriptome, and proteome signatures. We noted similar profiles in ISF, serum, and plasma samples, suggesting that ISF can be a proxy for direct blood sampling. Dynamic changes in RNA-seq were recorded in ISF from induced hypoxia conditions. Finally, we report the first isolation and characterization, to our knowledge, of exosomes from dermal ISF. The ISF exosome concentration is 12–13 times more enriched when compared to plasma and serum and represents a previously unexplored biofluid for exosome isolation. This minimally invasive extraction approach can enable mechanistic studies of ISF and demonstrates the potential of ISF for real-time health monitoring applications.

59 BASIC BIOLOGICAL SCIENCES↗

Predictive Pathogen Biology: Genome-Based Prediction of Pathogenic Potential and Countermeasures Targets

Bacterial pathogens have numerous processes by which their genomic DNA is acquired or rearranged as part of their normal physiology (e.g., exchange of plasmids through conjugation) or by bacteriophage that parasitize bacteria and often insert into the bacterial genome as prophages. These processes occur with relatively high probability/frequency, and may lead to sudden changes in virulence, as new genetic material is added to the chromosome, or structural changes in the chromosome affect gene expression. We set out to devise methods to measure the rates of these processes in bacteria using next generation DNA sequencing. Using very deep sequencing on genomes we had assembled, using library preparation methods and bioinformatics tools designed to help find mobile elements and signs of their insertion, we were able to find numerous examples of attempted novel genome arrangements, revealing data that can be used to calculate rates of different mechanisms of genome change.

59 BASIC BIOLOGICAL SCIENCES↗

IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets

IslandViewer is a widely-used webserver for the prediction and interactive visualization of genomic islands (GIs, regions of probable horizontal origin) in bacterial and archaeal genomes. GIs disproportionately encode factors that enhance the adaptability and competitiveness of the microbe within a niche, including virulence factors and other medically or environmentally important adaptations. We report here the release of IslandViewer 4, with novel features to accommodate the needs of larger-scale microbial genomics analysis, while expanding GI predictions and improving its flexible visualization interface. A user management web interface as well as an HTTP API for batch analyses are now provided with a secured authentication to facilitate the submission of larger numbers of genomes and the retrieval of results. In addition, IslandViewer’s integrated GI predictions from multiple methods have been improved and expanded by integrating the precise Islander method for pre-computed genomes, as well as an updated IslandPath-DIMOB for both pre-computed and user-supplied custom genome analysis. Finally, pre-computed predictions including virulence factors and antimicrobial resistance are now available for 6193 complete bacterial and archaeal strains publicly available in RefSeq. IslandViewer 4 provides key enhancements to facilitate the analysis of GIs and better understand their role in the evolution of successful environmental microbes and pathogens.

59 BASIC BIOLOGICAL SCIENCES↗

Experimental single-strain mobilomics reveals events that shape pathogen emergence

Virulence and resistance genes carried on mobile DNAs such as genomic islands (GIs) and plasmids promote bacterial pathogen emergence. An early step in the mobilization of GIs is their excision, which produces both a circular form of the GI and a deletion site in the chromosome; circular forms have also been described for some bacterial insertion sequences (ISs). We demonstrate that the recombinant sequence produced at the junction of such circles, and their corresponding deletion sites, can be detected sensitively in high throughput sequencing data, using new computational methods that enable empirical discovery of new mobile DNAs. Applied to the rich mobilome of a single strain (Kpn2146) of the emerging multidrug-resistant pathogen Klebsiella pneumoniae, our approach detected circular junctions for six GIs and seven IS types (several of the latter not previously known to circularize). Our methods further revealed differential biology of multiple mobile DNAs, imprecision of integrases and transposases, and differential activity among identical IS copies for IS26, ISKpn18 and ISKpn21. Exonuclease was used to enrich for circular dsDNA molecules, and internal calibration with the native Kpn2146 plasmids showed that not all molecules bearing GI and IS circular junctions were circular dsDNAs. Transposition events were also detected, revealing replicon preference (ISKpn18 preferring a conjugative IncA/C2 plasmid), local action (IS26), regional preferences, selection (against capsule synthesis), and left-right IS end swapping. Efficient discovery and global characterization of numerous mobile elements per experiment will allow detailed accounting of bacterial evolution, explaining the new gene combinations that arise in emerging pathogens.

60 APPLIED LIFE SCIENCES↗