Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 343 records · Page 19

Conservation of Transcription Start Sites within Genes across a Bacterial Genus

Transcription start sites (TSSs) lying inside annotated genes, on the same or opposite strand, have been observed in diverse bacteria, but the function of these unexpected transcripts is unclear. Here, we use the metal-reducing bacterium Shewanella oneidensis MR-1 and its relatives to study the evolutionary conservation of unexpected TSSs. Using high-resolution tiling microarrays and 5'-end RNA sequencing, we identified 2,531 TSSs in S. oneidensis MR-1, of which 18% were located inside coding sequences (CDSs). Comparative transcriptome analysis with seven additional Shewanella species revealed that the majority (76%) of the TSSs within the upstream regions of annotated genes (gTSSs) were conserved. Thirty percent of the TSSs that were inside genes and on the sense strand (iTSSs) were also conserved. Sequence analysis around these iTSSs showed conserved promoter motifs, suggesting that many iTSS are under purifying selection. Furthermore, conserved iTSSs are enriched for regulatory motifs, suggesting that they are regulated, and they tend to eliminate polar effects, which confirms that they are functional. In contrast, the transcription of antisense TSSs located inside CDSs (aTSSs) was significantly less likely to be conserved (22%). However, aTSSs whose transcription was conserved often have conserved promoter motifs and drive the expression of nearby genes. Overall, our findings demonstrate that some internal TSSs are conserved and drive protein expression despite their unusual locations, but the majority are not conserved and may reflect noisy initiation of transcription rather than a biological function.

bacterial genus↗

MetaHCR: a web-enabled metagenome data management system for hydrocarbon resources

The ever-increasing metagenomic data necessitate appropriate cataloguing in a way that facilitates the comparison and better contextualization of the underlying investigations. To this extent, information associated with the sequencing data as well as the original sample and the environment where it was obtained from is crucial. To date, there are not any publicly available repositories able to capture environmental metadata pertaining to hydrocarbon-rich environments. As such, contextualization and comparative analysis among sequencing datasets derived from these environments is to a certain degree hindered or cannot be fully evaluated. The metagenomics data management system for hydrocarbon resources (MetaHCRs) enables the capturing of marker gene and whole metagenome sequencing data as well as over 300 contextual attributes associated with samples, organisms, environments and geological properties, among others. Moreover, MetaHCR implements the Minimum Information about any Sequence–hydrocarbon resource specification from the Genomic Standards Consortium; it integrates a user-friendly web interface and relational database model, and it enables the generation of complex custom search. MetaHCR has been tested with 36 publicly available metagenomic studies, and its modular architecture can be easily customized for other types of environmental and metagenomics studies.

59 BASIC BIOLOGICAL SCIENCES↗

A Biosequence-based Approach to Software Characterization

For many applications, it is desirable to have some process for recognizing when software binaries are closely related without relying on them to be identical or have identical segments. Some examples include monitoring utilization of high performance computing centers or service clouds, detecting freeware in licensed code, and enforcing application whitelists. But doing so in a dynamic environment is a nontrivial task because most approaches to software similarity require extensive and time-consuming analysis of a binary, or they fail to recognize executables that are similar but nonidentical. Presented herein is a novel biosequence-based method for quantifying similarity of executable binaries. Using this method, it is shown in an example application on large-scale multi-author codes that 1) the biosequence-based method has a statistical performance in recognizing and distinguishing between a collection of real-world high performance computing applications better than 90% of ideal; and 2) an example of using family tree analysis to tune identification for a code subfamily can achieve better than 99% of ideal performance.

cyber security↗

A comparison of methods used to unveil the genetic and metabolic pool in the built environment

BACKGROUND:A majority of indoor residential microbes originate from humans, pets, and outdoor air and are not adapted to the built environment (BE). Consequently, a large portion of the microbes identified by DNA-based methods are either dead or metabolically inactive. Although many exceptions have been noted, the ribosomal RNA fraction of the sample is more likely to represent either viable or metabolically active cells. We examined methodological variations in sample processing using a defined, mock BE microbial community to better understand the scope of technique-based vs. biological-based differences in both ribosomal transcript (rRNA) and gene (DNA) sequence community analysis. Based on in vitro tests, a protocol was adopted for the analysis of the genetic and metabolic pool (DNA vs. rRNA) of air and surface microbiomes within a residential setting. RESULTS:We observed differences in DNA/RNA co-extraction efficiency for individual microbes, but overall, a greater recovery of rRNA using FastPrep (> 50%). Samples stored with various preservation methods at - 80°C experienced a rapid decline in nucleic acid recovery starting within the first week, although post-extraction rRNA had no significant degradation when treated with RNAStable. We recommend that co-extraction samples be processed as quickly as possible after collection. The in vivo analysis revealed significant differences in the two components (genetic and metabolic pool) in terms of taxonomy, community structure, and microbial association networks. Rare taxa present in the genetic pool showed higher metabolic potential (RNA:DNA ratio), whereas commonly detected taxa of outdoor origins based on DNA sequencing, especially taxa of the Sphingomonadales order, were present in lower relative abundances in the viable community. CONCLUSIONS:Although methodological variations in sample preparations are high, large differences between the DNA and RNA fractions of the total microbial community demonstrate that direct examination of rRNA isolated from a residential BE microbiome has the potential to identify the more likely viable or active portion of the microbial community. In an environment that has primarily dead and metabolically inactive cells, we suggest that the rRNA fraction of BE samples is capable of providing a more ecologically relevant insight into the factors that drive indoor microbial community dynamics.

59 BASIC BIOLOGICAL SCIENCES↗

Mechanism of H 2 S Oxidation by the Dissimilatory Perchlorate-Reducing Microorganism Azospira suillum PS

The genetic and biochemical basis of perchlorate-dependent H 2 S oxidation (PSOX) was investigated in the dissimilatory perchlorate-reducing microorganism (DPRM) Azospira suillum PS (PS). Previously, it was shown that all known DPRMs innately oxidize H 2 S, producing elemental sulfur (S o ). Although the process involving PSOX is thermodynamically favorable (ΔG ° ' = -206 kJ ∙ mol -1 H 2 S), the underlying biochemical and genetic mechanisms are currently unknown. Interestingly, H 2 S is preferentially utilized over physiological electron donors such as lactate or acetate although no growth benefit is obtained from the metabolism. Here, we determined that PSOX is due to a combination of enzymatic and abiotic interactions involving reactive intermediates of perchlorate respiration. Using various approaches, including barcode analysis by sequencing (Bar-seq), transcriptome sequencing (RNA-seq), and proteomics, along with targeted mutagenesis and biochemical characterization, we identified all facets of PSOX in PS. In support of our proposed model, deletion of identified upregulated PS genes traditionally known to be involved in sulfur redox cycling (e.g., Sox, sulfide:quinone reductase [SQR]) showed no defect in PSOX activity. Proteomic analysis revealed differential abundances of a variety of stress response metal efflux pumps and divalent heavy-metal transporter proteins, suggesting a general toxicity response. Furthermore, in vitro biochemical studies demonstrated direct PSOX mediated by purified perchlorate reductase (PcrAB) in the absence of other electron transfer proteins. The results of these studies support a model in which H 2 S oxidation is mediated by electron transport chain short-circuiting in the periplasmic space where the PcrAB directly oxidizes H 2 S to S o . The biogenically formed reactive intermediates (ClO 2 - and O 2 ) subsequently react with additional H 2 S, producing polysulfide and S o as end products.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

Multiple Lines of Evidences Reveal Mechanisms Underpinning Mercury Resistance and Volatilization by Stenotrophomonas sp. MA5 Isolated from the Savannah River Site (SRS), USA

A largely understudied microbially mediated mercury (Hg) bioremediative pathway includes the volatilization of Hg 2+ to Hg 0 . Therefore, studies on Hg resistant bacteria (HgR), isolated from historically long-term contaminated environments, can serve as models to understand mechanisms underpinning Hg cycling. Towards this end, a mercury resistant bacterial strain, identified as Stenotrophomonas sp., strain MA5, was isolated from Mill Branch on the Savannah River Site (SRS); an Hg-impacted ecosystem. Minimum inhibitory concentration (MIC) analysis showed Hg resistance of up to 20 µg/mL by MA5 with 95% of cells retaining viability. Microcosm studies showed that the strain depleted more than 90% of spiked Hg 2+ within the first 24 h of growth and the detection of volatilized mercury indicated that the strain was able to reduce Hg 2+ to Hg 0 . To understand molecular mechanisms of Hg volatilization, a draft whole genome sequence was obtained, annotated and analyzed, which revealed the presence of a transposon-derived mer operon (merRTPADE) in MA5, known to transport and reduce Hg 2+ into Hg 0 . Based on the whole genome sequence of strain MA5, qRT-PCR assays were designed on merRTPADE, we found a ~40-fold higher transcription of mer T, P, A, D and E when cells were exposed to 5 µg/mL Hg 2+ . Interestingly, strain MA5 increased cellular size as a function of increasing Hg concentrations, which is likely an evolutionary response mechanism to cope with Hg stress. Moreover, metal contaminated environments are shown to co-select for antibiotic resistance. When MA5 was screened for antibiotic resistance, broad resistance against penicillin, streptomycin, tetracycline, ampicillin, rifampicin, and erythromycin was found; this correlated with the presence of multiple gene determinants for antibiotic resistance within the whole genome sequence of MA5. Overall, this study provides an in-depth understanding of the underpinnings of Stenotrophomonas-mercury interactions that facilitate cellular survival in a contaminated soil habitat.

59 BASIC BIOLOGICAL SCIENCES↗

Deletion analysis of the itaconic acid biosynthesis gene cluster components in Aspergillus pseudoterreus ATCC32359

The filamentous fungus Aspergillus terreus has been successfully used for industrial production of itaconic acid (IA) for many years. The IA biosynthesis pathway has recently been characterized at a molecular genetic level as an IA gene cluster (IA cluster) by a clone-based transcriptomic approach. The cluster consists of four genes, including genes for cis-aconitic acid decarboxylase (cadA), a predicted transcription factor (tf), mitochondrial organic acid transporter (mttA), and MFS (Major Facilitator Superfamily) type transporter (mfsA). In this research, we performed EST (Expressed Sequence Tag) analysis and systematic gene deletions to investigate the role of those genes in IA biosynthesis in A. pseudoterreus ATCC32359. EST analysis showed a similar expression pattern among those four genes that was distinct from neighbouring genes and further confirmed that they belong to the same biosynthesis cluster. Systematic gene deletion analysis demonstrated that tf, cadA, mttA, and mfsA genes in the cluster are essential for IA production. Deletion of any of them will either completely abolish the IA production or dramatically decrease amount of IA produced. More importantly, a significant amount of aconitic acid was detected in the cadA deletion strain but not in the other deletion strains. Therefore, the cadA deletion strain is a novel microorganism host for production of aconitic acid and other value-added chemicals from lignocellulosic biomass.

Deng, Shuang↗

ClassyFlu: Classification of Influenza A Viruses with Discriminatively Trained Profile-HMMs

Accurate and rapid characterization of influenza A virus (IAV) hemagglutinin (HA) and neuraminidase (NA) sequences with respect to subtype and clade is at the basis of extended diagnostic services and implicit to molecular epidemiologic studies. ClassyFlu is a new tool and web service for the classification of IAV sequences of the HA and NA gene into subtypes and phylogenetic clades using discriminatively trained profile hidden Markov models (HMMs), one for each subtype or clade. ClassyFlu merely requires as input unaligned, full-length or partial HA or NA DNA sequences. It enables rapid and highly accurate assignment of HA sequences to subtypes H1–H17 but particularly focusses on the finer grained assignment of sequences of highly pathogenic avian influenza viruses of subtype H5N1 according to the cladistics proposed by the H5N1 Evolution Working Group. NA sequences are classified into subtypes N1–N10. ClassyFlu was compared to semiautomatic classification approaches using BLAST and phylogenetics and additionally for H5 sequences to the new ‘‘Highly Pathogenic H5N1 Clade Classification Tool’’ (IRD-CT) proposed by the Influenza Research Database. Our results show that both web tools (ClassyFlu and IRD-CT), although based on different methods, are nearly equivalent in performance and both are more accurate and faster than semiautomatic classification. A retraining of ClassyFlu to altered cladistics as well as an extension of ClassyFlu to other IAV genome segments or fragments thereof is undemanding. This is exemplified by unambiguous assignment to a distinct cluster within subtype H7 of sequences of H7N9 viruses which emerged in China early in 2013 and caused more than 130 human infections.

59 BASIC BIOLOGICAL SCIENCES↗

IMG-ABC: An Atlas of Biosynthetic Gene Clusters to Fuel the Discovery of Novel Secondary Metabolites

In the discovery of secondary metabolites (SMs), large-scale analysis of sequence data is a promising exploration path that remains largely underutilized due to the lack of relevant computational resources. We present IMG-ABC (https://img.jgi.doe.gov/abc/) -- An Atlas of Biosynthetic gene Clusters within the Integrated Microbial Genomes (IMG) system1. IMG-ABC is a rich repository of both validated and predicted biosynthetic clusters (BCs) in cultured isolates, single-cells and metagenomes linked with the SM chemicals they produce and enhanced with focused analysis tools within IMG. The underlying scalable framework enables traversal of phylogenetic dark matter and chemical structure space -- serving as a doorway to a new era in the discovery of novel molecules.

99 GENERAL AND MISCELLANEOUS↗

Metatranscriptomic insights on gene expression and regulatory controls in Candidatus Accumulibacter phosphatis

Previous studies on enhanced biological phosphorus removal (EBPR) have focused on reconstructing genomic blueprints for the model polyphosphate-accumulating organism Candidatus Accumulibacter phosphatis. Here, a time series metatranscriptome generated from enrichment cultures of Accumulibacter was used to gain insight into anerobic/aerobic metabolism and regulatory mechanisms within an EBPR cycle. Co-expressed gene clusters were identified displaying ecologically relevant trends consistent with batch cycle phases. Transcripts displaying increased abundance during anerobic acetate contact were functionally enriched in energy production and conversion, including upregulation of both cytoplasmic and membrane-bound hydrogenases demonstrating the importance of transcriptional regulation to manage energy and electron flux during anerobic acetate contact. We hypothesized and demonstrated hydrogen production after anerobic acetate contact, a previously unknown strategy for Accumulibacter to maintain redox balance. Genes involved in anerobic glycine utilization were identified and phosphorus release after anerobic glycine contact demonstrated, suggesting that Accumulibacter routes diverse carbon sources to acetyl-CoA formation via previously unrecognized pathways. A comparative genomics analysis of sequences upstream of co-expressed genes identified two statistically significant putative regulatory motifs. One palindromic motif was identified upstream of genes involved in PHA synthesis and acetate activation and is hypothesized to be a phaR binding site, hence representing a hypothetical PHA modulon. A second motif was identified ~35 base pairs (bp) upstream of a large and diverse array of genes and hence may represent a sigma factor binding site. As a result, this analysis provides a basis and framework for further investigations into Accumulibacter metabolism and the reconstruction of regulatory networks in uncultured organisms.

54 ENVIRONMENTAL SCIENCES↗

Combined theoretical and experimental study of the valence, Rydberg and ionic states of fluorobenzene

New photoelectron spectra (PES) and ultra violet (UV) and vacuum UV (VUV) absorption spectra of fluorobenzene recorded at higher resolution than previously, have been combined with mass-resolved (2 + 1) and (3 + 1) resonance enhanced multiphoton ionization (REMPI) spectra; this has led to the identification of numerous Rydberg states. The PES have been compared with earlier mass-analyzed threshold ionization and photoinduced Rydberg ionization (PIRI) spectra to give an overall picture of the ionic state sequence. The analysis of these spectra using both equations of motion with coupled cluster singles and doubles (EOM-CCSD) configuration interaction and time dependent density functional theory (TDDFT) calculations have been combined with vibrational analysis of both the hot and cold bands of the spectra, in considerable detail. The results extend several earlier studies on the vibronic coupling leading to conical intersections between the X{sup 2}B{sub 1} and A{sup 2}A{sub 2} states, and a further trio (B, C, and D) of states. The conical intersection of the X and A states has been explicitly identified, and its structure and energetics evaluated. The energy sequence of the last group is only acceptable to the present study if given as B{sup 2}B{sub 2}<C{sup 2}B{sub 1}<D{sup 2}A{sub 1}, a conclusion which is in agreement with most previous EOM-CCSD and other calculations. However, this symmetry ordering of the B and C states forces reconsideration of the nature of the PIRI spectrum. The coupling between these two states is induced by the a{sub 2} modes, ν{sub 12} and ν{sub 14} and we propose that the 14{sup 1} band is observed in the B{sup 2}B{sub 2} band in the PES for the first time, because of the improved resolution. This same assignment is given to the lowest energy band in the PIRI spectrum which was previously assigned as the origin band and further conclude that the entire PIRI spectrum is induced by ν{sub 12} and ν{sub 14}. The relative intensities of the various Rydberg state peaks in the VUV absorption and REMPI spectra of fluorobenzene are very similar to those observed in the equivalent spectra of benzene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A deep learning approach to real-time HIV outbreak detection using genetic data

Pathogen genomic sequence data are increasingly made available for epidemiological monitoring. A main interest is to identify and assess the potential of infectious disease outbreaks. While popular methods to analyze sequence data often involve phylogenetic tree inference, they are vulnerable to errors from recombination and impose a high computational cost, making it difficult to obtain real-time results when the number of sequences is in or above the thousands. Here, we propose an alternative strategy to outbreak detection using genomic data based on deep learning methods developed for image classification. The key idea is to use a pairwise genetic distance matrix calculated from viral sequences as an image, and develop convolutional neutral network (CNN) models to classify areas of the images that show signatures of active outbreak, leading to identification of subsets of sequences taken from an active outbreak. We showed that our method is efficient in finding HIV-1 outbreaks with R0 ≥ 2.5, and overall a specificity exceeding 98% and sensitivity better than 92%. We validated our approach using data from HIV-1 CRF01 in Europe, containing both endemic sequences and a well-known dual outbreak in intravenous drug users. Our model accurately identified known outbreak sequences in the background of slower spreading HIV. Importantly, we detected both outbreaks early on, before they were over, implying that had this method been applied in real-time as data became available, one would have been able to intervene and possibly prevent the extent of these outbreaks. This approach is scalable to processing hundreds of thousands of sequences, making it useful for current and future real-time epidemiological investigations, including public health monitoring using large databases and especially for rapid outbreak identification.

59 BASIC BIOLOGICAL SCIENCES↗

A cost-benefit analysis to select the most effective method for positional cloning: genotyping by sequencing versus allele-specific PCR

Next-generation sequencing technologies have enabled the rapid generation of high-resolution genetic maps to enable the identification of genomic regions associated with traits of interest. In certain cases these regions of interest need to be narrowed down via fine mapping to increase resolution to the level of the individual gene. In those cases, there is a choice of genotyping options, whereby cost, time, and information content need to be considered. We developed detailed, customizable models to compare the cost of genotyping by sequencing and allele-specific PCR for fine-mapping genomic regions of interest. The models were validated experimentally with data from a fine-mapping experiment designed to identify candidate genes within a 3.2-Mb disease resistance locus in the cereal crop sorghum (Sorghum bicolor (L.). Moench), which has a sequenced genome. The size of the mapping population, genetic diversity of the parental lines, choice of reagents, and labor cost were shown to influence the overall cost of the two methods, and, consequently, which method would be most cost-effective for a given experiment. Here, the model can be easily customized to reflect experiments with different species, genetic populations, experimental and personnel costs to determine the most cost-effective procedure for a new positional cloning project.

Cost analysis↗

Impacts of Glutaraldehyde on Microbial Community Structure and Degradation Potential in Streams Impacted by Hydraulic Fracturing

The environmental impacts of hydraulic fracturing, particularly those of surface spills in aquatic ecosystems, are not fully understood. The goals of this study were to (1) understand the effect of previous exposure to hydraulic fracturing fluids on aquatic microbial community structure and (2) examine the impacts exposure has on biodegradation potential of the biocide glutaraldehyde. Microcosms were constructed from hydraulic fracturing-impacted and nonhydraulic fracturing-impacted streamwater within the Marcellus shale region in Pennsylvania. Microcosms were amended with glutaraldehyde and incubated aerobically for 56 days. Microbial community adaptation to glutaraldehyde was monitored using 16S rRNA gene amplicon sequencing and quantification by qPCR. Abiotic and biotic glutaraldehyde degradation was measured using ultra-performance liquid chromatography--high resolution mass spectrometry and total organic carbon. It was found that nonhydraulic fracturing-impacted microcosms biodegraded glutaraldehyde faster than the hydraulic fracturing-impacted microcosms, showing a decrease in degradation potential after exposure to hydraulic fracturing activity. Hydraulic fracturing-impacted microcosms showed higher richness after glutaraldehyde exposure compared to unimpacted streams, indicating an increased tolerance to glutaraldehyde in hydraulic fracturing impacted streams. Beta diversity and differential abundance analysis of sequence count data showed different bacterial enrichment for hydraulic fracturing-impacted and nonhydraulic fracturing-impacted microcosms after glutaraldehyde addition. These findings demonstrated a lasting effect on microbial community structure and glutaraldehyde degradation potential in streams impacted by hydraulic fracturing operations.

54 ENVIRONMENTAL SCIENCES↗

Genomic, Transcriptional, and Phenotypic Analysis of the Glucose Derepressed Clostridium beijerinckii Mutant Exhibiting Acid Crash Phenotype

Clostridium beijerinckii is a predominant solventogenic bacterium that is used for the ABE fermentation. Various C. beijerinckii mutants are constructed for desirable phenotypes. The C. beijerinckii mutant BA105 harboring a glucose derepression phenotype was previously isolated and demonstrated the enhanced amylolytic activity in the presence of glucose. Despite its potential use, BA105 is not further characterized and utilized. Therefore, the authors investigate fermentation phenotypes of BA105 in this study. Under the typical batch fermentation conditions, BA105 consistently exhibits acid crash phenotype resulting in limited glucose uptake and cell growth. However, when the culture pH is maintained above 5.5, BA105 exhibits the increased glucose uptake and butanol production than did the wild‐type. To further analyze BA105, the authors perform genome sequencing and RNA sequencing. Genome analysis identifies two SNPs unique to BA105, in the upstream region of AbrB regulator (Cbei_4885) and the ROK family glucokinase (Cbei_4895) which are involved in catabolite repression and regulation of sugar metabolism. Transcriptional analysis of BA105 reveals significant differential expression of the genes associated with the PTS sugar transport system and acid production. This study improves understanding of the acid crash phenomenon and provides the genetic basis underlying the catabolite derepression phenotype of C. beijericnkii .

Seo, Seung‐Oh↗

Developing an Automated Uncertainty Quantification Tool to Improve Watershed-Scale Predictions of Water and Nutrient Cycling

Managing the flow of water, nutrients, and contaminants in watersheds is vital to addressing pressing issues related to water scarcity, access to clean drinking water, energy production, resilience to natural and anthropogenic perturbations, and ecological restoration. Decisions about the management of watersheds critically depend on the accuracy with which the flow of water and chemicals through the watershed can be predicted by computer models. Prediction uncertainty can be reduced by matching the model to data, which are collected in the field at great expense. The contribution of watershed characterization data to reducing uncertainty of relevant model predictions can be evaluated in a so-called data-worth analysis, which provides transparent, quantitative metrics about a data set’s value for the support of relevant watershed management objectives. To achieve this goal, we developed a software package that implements the data-worth analysis approach for use with state-of-the-art watershed models. The purpose of the proposed data-worth analysis is to help decision-makers allocate resources for watershed characterization such that the uncertainty in model predictions can be significantly reduced, which leads to better, more effective management decisions. At the same time, watershed characterization costs can be reduced. The specific technical objectives of this SBIR/STTR Phase II project were to develop a framework and associated software toolsets that implement the uncertainty quantification and data-worth analysis approach for use with state-of-the-art watershed models. This goal was achieved by (A) developing a user-friendly, robust software package that is accessible to a wide audience, including watershed managers, policy-makers, and public stakeholders; (B) by demonstrating application of the prototype on several use cases that are representative of complex watershed management challenges spanning a range of scales and that consider different open-source, DOE-based codes and other modeling platforms; and (C) by gathering information about the needs and requirements from potential users to help guide future developments, ensuring that the final product will be commercially viable. The developed software consists of a graphical user interface that guides the user through a sequence of analysis steps, supported by toolsets that leverage state-of-the-art computational simulation-optimization capabilities. A prototype of the software runs on multiple platforms (PC, Mac, multi-processor Linux environment), is linked to diverse watershed simulators (e.g., ECOSYS, TOUGH2, TOUGHREACT, Amanzi-ATS), performs multiple analysis tasks (predictive simulations, sensitivity analysis, uncertainty analysis, automatic parameter estimation, and data-worth analysis, multicomponent geothermometry), and is readily extensible to include external simulators and analysis tools. The software is being commercialized and will be continually updated to address user needs.

58 GEOSCIENCES↗

Structure and mutagenic analysis of the lipid II flippase MurJ from Escherichia coli

The peptidoglycan cell wall provides an essential protective barrier in almost all bacteria, defining cellular morphology and conferring resistance to osmotic stress and other environmental hazards. The precursor to peptidoglycan, lipid II, is assembled on the inner leaflet of the plasma membrane. However, peptidoglycan polymerization occurs on the outer face of the plasma membrane, and lipid II must be flipped across the membrane by the MurJ protein before its use in peptidoglycan synthesis. Due to its central role in cell wall assembly, MurJ is of fundamental importance in microbial cell biology and is a prime target for novel antibiotic development. However, relatively little is known regarding the mechanisms of MurJ function, and structural data for MurJ are available only from the extremophile Thermosipho africanus. Here, we report the crystal structure of substrate-free MurJ from the gram-negative model organism Escherichia coli, revealing an inward-open conformation. Taking advantage of the genetic tractability of E. coli, we performed high-throughput mutagenesis and next-generation sequencing to assess mutational tolerance at every amino acid in the protein, providing a detailed functional and structural map for the enzyme and identifying sites for inhibitor development. Finally, through the use of sequence coevolution analysis, we identify functionally important interactions in the outward-open state of the protein, supporting a rocker-switch model for lipid II transport.

59 BASIC BIOLOGICAL SCIENCES↗