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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Generating a high-confidence reference genome map of the Greater Duckweed by integration of cytogenomic, optical mapping, and Oxford Nanopore technologies

Duckweeds are the fastest growing angiosperms and have the potential to become a new generation of sustainable crops. Although a seed plant, Spirodela polyrhiza clones rarely flower and multiply mainly through vegetative propagation. Whole-genome sequencing using different approaches and clones yielded two reference maps. One for clone 9509, supported in its assembly by optical mapping of single DNA molecules, and one for clone 7498, supported by cytogenetic assignment of 96 fingerprinted bacterial artificial chromosomes (BACs) to its 20 chromosomes. However, these maps differ in the composition of several individual chromosome models. We validated both maps further to resolve these differences and addressed whether they could be due to chromosome rearrangements in different clones. For this purpose, we applied sequential multicolor fluorescence in situ hybridization (mcFISH) to seven S. polyrhiza clones, using 106 BACs that were mapped onto the 39 pseudomolecules for clone 7498. Furthermore we integrated high-depth Oxford Nanopore (ON) sequence data for clone 9509 to validate and revise the previously assembled chromosome models. We found no major structural rearrangements between these seven clones, identified seven chimeric pseudomolecules and Illumina assembly errors in the previous maps, respectively. A new S. polyrhiza genome map with high contiguity was produced with the ON sequence data and genome-wide synteny analysis supported the occurrence of two Whole Genome Duplication events during its evolution. This work generated a high confidence genome map for S. polyrhiza at the chromosome scale, and illustrates the complementarity of independent approaches to produce whole-genome assemblies in the absence of a genetic map.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Crewmember microbiome may influence microbial composition of ISS habitable surfaces

The International Space Station (ISS) is a complex built environment physically isolated from Earth. Assessing the interplay between the microbial community of the ISS and its crew is important for preventing biomedical and structural complications for long term human spaceflight missions. In this study, we describe one crewmember’s microbial profile from body swabs of mouth, nose, ear, skin and saliva that were collected at eight different time points pre-, during and post-flight. Additionally, environmental surface samples from eight different habitable locations in the ISS were collected from two flights. Environmental samples from one flight were collected by the crewmember and samples from the next flight were collected after the crewmember departed. The microbial composition in both environment and crewmember samples was measured using shotgun metagenomic sequencing and processed using the Livermore Metagenomics Analysis Toolkit. Ordination of sample to sample distances showed that of the eight crew body sites analyzed, skin, nostril, and ear samples are more similar in microbial composition to the ISS surfaces than mouth and saliva samples; and that the microbial composition of the crewmember’s skin samples are more closely related to the ISS surface samples collected by the crewmember on the same flight than ISS surface samples collected by other crewmembers on different flights. In these collections, species alpha diversity in saliva samples appears to decrease during flight and rebound after returning to Earth. This is the first study to compare the ISS microbiome to a crewmember’s microbiome via shotgun metagenomic sequencing. We observed that the microbiome of the surfaces inside the ISS resemble those of the crew’s skin. These data support future crew and ISS microbial surveillance efforts and the design of preventive measures to maintain crew habitat onboard spacecraft destined for long term space travel.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1991

Hierarchic models for laminated plates

Structural plates and shells are three-dimensional bodies, one dimension of which happens to be much smaller than the other two. Thus, the quality of a plate or shell model must be judged on the basis of how well its exact solution approximates the corresponding three-dimensional problem. Of course, the exact solution depends not only on the choice of the model but also on the topology, material properties, loading and constraints. The desired degree of approximation depends on the analyst's goals in performing the analysis. For these reasons models have to be chosen adaptively. Hierarchic sequences of models make adaptive selection of the model which is best suited for the purposes of a particular analysis possible. The principles governing the formulation of hierarchic models for laminated plates are presented. The essential features of the hierarchic models described models are: (1) the exact solutions corresponding to the hierarchic sequence of models converge to the exact solution of the corresponding problem of elasticity for a fixed laminate thickness; and (2) the exact solution of each model converges to the same limit as the exact solution of the corresponding problem of elasticity with respect to the laminate thickness approaching zero. The formulation is based on one parameter (beta) which characterizes the hierarchic sequence of models, and a set of constants whose influence was assessed by a numerical sensitivity study. The recommended selection of these constants results in the number of fields increasing by three for each increment in the power of beta. Numerical examples analyzed with the proposed sequence of models are included and good correlation with the reference solutions was found. Results were obtained for laminated strips (plates in cylindrical bending) and for square and rectangular plates with uniform loading and with homogeneous boundary conditions. Cross-ply and angle-ply laminates were evaluated and the results compared with those of MSC/PROBE. Hierarchic models make the computation of any engineering data possible to an arbitrary level of precision within the framework of the theory of elasticity.

Actis, Ricardo Luis↗
DOE OSTI2016

Full core LOCA analysis for BWR/6 - Methodology and first results

As part of collaboration between the Paul Scherrer Institute (PSI) and the Swiss regulator ENSI, a study aimed at core-wide estimates of the fuel rod behaviour during large-break loss-of-coolant-accidents (LB-LOCA) in BWR/6 plants was recently conducted on the basis of a Swiss reactor. To perform such full core LOCA analysis, a methodology consisting of a multi-physics sequence of interfaced best-estimate models and codes was established: CASMO-4E/SIMULATE-3 (C4/S3) for the core modelling, TRACE for the plant system analysis and FALCON for the fuel thermomechanical calculations during base irradiation (BI) as well as for the transient analyses (TA). As starting point, reference C4/S3 models validated for all operated cycles of the given plant were used to select as situation target, a recent cycle representative of modern BWR core loadings and fuel assembly designs. Based on the selected C4/S3 cycle model, PSI in-house computational modules were developed or updated in order to set up the TRACE full core model with regards to assembly geometries, fuel rods thermophysical properties and initial 3-D power distributions. Also, an interface to FALCON BI models was established for the transfer of the irradiation history, including primarily linear heat generation rates (LHGR), for an equivalent fuel pin. For each assembly, the power of the equivalent pin is defined with the maximum power of all pins at each axial position. However, in order to reduce the computational time of the fuel performance analyses, the whole core was subdivided into groups of assemblies, each represented by the most powerful assembly within the group. On this basis, the FALCON BI analyses were conducted with coupling to the PSI advanced model for gas release and swelling (GRWS-A) and an interface was developed to transfer the calculated fuel rod characteristics at end of irradiation to initialize the core fuel rods heat structures in the TRACE model. Subsequently, double-ended recirculation line rupture scenarios at hot-full-power (HFP) conditions were investigated with TRACE in order to determine a conservative transient case as function of cycle burnup, single failure criterion and availability of safety systems. For the most conservative case of this study, an interface was set-up to transfer the TRACE results as boundary conditions to the FALCON TA models, including thus primarily time and axial distributions of the clad outer temperatures. As final step, the FALCON TA were conducted for the representative assemblies in order to infer core-wide estimates of the fuel rod behaviour with specific focus on equivalent cladding reacted levels, hoop strains as well as margins to failure criteria. The objective of this paper is to present this methodology along with the obtained first results. Regarding the latter, it is found that principally, even for the most conservative LOCA scenario of this study, the clad temperature excursion remains mild enough at all core locations such that no significant rod ballooning and no cladding failure is predicted, something in line with corresponding analyses conducted earlier for BWR/4 plants. (authors)

42 ENGINEERING↗
DOE OSTI2017

Stable Isotope Labeling for Improved Comparative Analysis of RNA Digests by Mass Spectrometry

Even with the advent of high throughput methods to detect modified ribonucleic acids (RNAs), mass spectrometry remains a reliable method to detect, characterize, and place post-transcriptional modifications within an RNA sequence. Here we have developed a stable isotope labeling comparative analysis of RNA digests (SIL-CARD) approach, which improves upon the original {sup 18}O/{sup 16}O labeling CARD method. Like the original, SIL-CARD allows sequence or modification information from a previously uncharacterized in vivo RNA sample to be obtained by direct comparison with a reference RNA, the sequence of which is known. This reference is in vitro transcribed using a {sup 13}C/{sup 15}N isotopically enriched nucleoside triphosphate (NTP). The two RNAs are digested with an endonuclease, the specificity of which matches the labeled NTP used for transcription. As proof of concept, several transfer RNAs (tRNAs) were characterized by SIL-CARD, where labeled guanosine triphosphate was used for the reference in vitro transcription. RNase T1 digestion products from the in vitro transcript will be 15 Da higher in mass than the same digestion products from the in vivo tRNA that are unmodified, leading to a doublet in the mass spectrum. Singlets, rather than doublets, arise if a sequence variation or a post-transcriptional modification is present that results in a relative mass shift different from 15 Da. Moreover, the use of the in vitro synthesized tRNA transcript allows for quantitative measurement of RNA abundance. Overall, SIL-CARD simplifies data analysis and enhances quantitative RNA modification mapping by mass spectrometry. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Optimizing metaproteomics database construction: lessons from a study of the vaginal microbiome

Metaproteomics, a method for untargeted, high-throughput identification of proteins in complex samples, provides functional information about microbial communities and can tie functions to specific taxa. Metaproteomics often generates less data than other omics techniques, but analytical workflows can be improved to increase usable data in metaproteomic outputs. Identification of peptides in the metaproteomic analysis is performed by comparing mass spectra of sample peptides to a reference database of protein sequences. Although these protein databases are an integral part of the metaproteomic analysis, few studies have explored how database composition impacts peptide identification. Here, we used cervicovaginal lavage (CVL) samples from a study of bacterial vaginosis (BV) to compare the performance of databases built using six different strategies. We evaluated broad versus sample-matched databases, as well as databases populated with proteins translated from metagenomic sequencing of the same samples versus sequences from public repositories. Smaller sample-matched databases performed significantly better, driven by the statistical constraints on large databases. Additionally, large databases attributed up to 34% of significant bacterial hits to taxa absent from the sample, as determined orthogonally by 16S rRNA gene sequencing. We also tested a set of hybrid databases which included bacterial proteins from NCBI RefSeq and translated bacterial genes from the samples. These hybrid databases had the best overall performance, identifying 1,068 unique human and 1,418 unique bacterial proteins, ~30% more than a database populated with proteins from typical vaginal bacteria and fungi. Our findings can help guide the optimal identification of proteins while maintaining statistical power for reaching biological conclusions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

A hybrid computational strategy to address WGS variant analysis in >5000 samples

Abstract Background The decreasing costs of sequencing are driving the need for cost effective and real time variant calling of whole genome sequencing data. The scale of these projects are far beyond the capacity of typical computing resources available with most research labs. Other infrastructures like the cloud AWS environment and supercomputers also have limitations due to which large scale joint variant calling becomes infeasible, and infrastructure specific variant calling strategies either fail to scale up to large datasets or abandon joint calling strategies. Results We present a high throughput framework including multiple variant callers for single nucleotide variant (SNV) calling, which leverages hybrid computing infrastructure consisting of cloud AWS, supercomputers and local high performance computing infrastructures. We present a novel binning approach for large scale joint variant calling and imputation which can scale up to over 10,000 samples while producing SNV callsets with high sensitivity and specificity. As a proof of principle, we present results of analysis on Cohorts for Heart And Aging Research in Genomic Epidemiology (CHARGE) WGS freeze 3 dataset in which joint calling, imputation and phasing of over 5300 whole genome samples was produced in under 6 weeks using four state-of-the-art callers. The callers used were SNPTools, GATK-HaplotypeCaller, GATK-UnifiedGenotyper and GotCloud. We used Amazon AWS, a 4000-core in-house cluster at Baylor College of Medicine, IBM power PC Blue BioU at Rice and Rhea at Oak Ridge National Laboratory (ORNL) for the computation. AWS was used for joint calling of 180 TB of BAM files, and ORNL and Rice supercomputers were used for the imputation and phasing step. All other steps were carried out on the local compute cluster. The entire operation used 5.2 million core hours and only transferred a total of 6 TB of data across the platforms. Conclusions Even with increasing sizes of whole genome datasets, ensemble joint calling of SNVs for low coverage data can be accomplished in a scalable, cost effective and fast manner by using heterogeneous computing platforms without compromising on the quality of variants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

Homology Modeling and Analysis of Structure Predictions of the Bovine Rhinitis B Virus RNA Dependent RNA Polymerase (RdRp)

Bovine Rhinitis B Virus (BRBV) is a picornavirus responsible for mild respiratory infection of cattle. It is probably the least characterized among the aphthoviruses. BRBV is the closest relative known to Foot and Mouth Disease virus (FMDV) with a ~43% identical polyprotein sequence and as much as 67% identical sequence for the RNA dependent RNA polymerase (RdRp), which is also known as 3D polymerase (3D pol ). In the present study we carried out phylogenetic analysis, structure based sequence alignment and prediction of three-dimensional structure of BRBV 3D pol using a combination of different computational tools. Model structures of BRBV 3D pol were verified for their stereochemical quality and accuracy. The BRBV 3D pol structure predicted by SWISS-MODEL exhibited highest scores in terms of stereochemical quality and accuracy, which were in the range of 2Å resolution crystal structures. The active site, nucleic acid binding site and overall structure were observed to be in agreement with the crystal structure of unliganded as well as template/primer (T/P), nucleotide tri-phosphate (NTP) and pyrophosphate (PPi) bound FMDV 3D pol (PDB, 1U09 and 2E9Z). The closest proximity of BRBV and FMDV 3D pol as compared to human rhinovirus type 16 (HRV-16) and rabbit hemorrhagic disease virus (RHDV) 3D pols is also substantiated by phylogeny analysis and root-mean square deviation (RMSD) between C-α traces of the polymerase structures. The absence of positively charged α-helix at C terminal, significant differences in non-covalent interactions especially salt bridges and CH-pi interactions around T/P channel of BRBV 3D pol compared to FMDV 3D pol , indicate that despite a very high homology to FMDV 3D pol , BRBV 3D pol may adopt a different mechanism for handling its substrates and adapting to physiological requirements. Our findings will be valuable in the design of structure-function interventions and identification of molecular targets for drug design applicable to Aphthovirus RdRps.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Tailoring Hierarchical Structure and Rare Earth Affinity of Compositionally Identical Polymers via Sequence Control

Macromolecule sequence, structure, and function are inherently intertwined. While well-established relationships exist in proteins, they are more challenging to define for synthetic polymer nanoparticles due to their molecular weight, sequence, and conformational dispersities. Furthermore, to explore the impact of sequence on nanoparticle structure, we synthesized a set of 16 compositionally identical, sequence-controlled polymers with distinct monomer patterning of dimethyl acrylamide and a bioinspired, structure-driving di(phenylalanine) acrylamide (FF). Sequence control was achieved through multiblock polymerizations, yielding unique ensembles of polymer sequences which were simulated by kinetic Monte Carlo simulations. Systematic analysis of the global (tertiary- and quaternary-like) structure in this amphiphilic copolymer series revealed the effect of multiple sequence descriptors: the number of domains, the hydropathy of terminal domains, and the patchiness (density) of FF within a domain, each of which impacted both chain collapse and the distribution of single- and multichain assemblies. Furthermore, both the conformational freedom of chain segments and local-scale, β-sheet-like interactions were sensitive to the patchiness of FF. To connect sequence, structure, and target function, we evaluated an additional series of nine sequence-controlled copolymers as sequestrants for rare earth elements (REEs) by incorporating a functional acrylic acid monomer into select polymer scaffolds. We identified key sequence variables that influence the binding affinity, capacity, and selectivity of the polymers for REEs. Collectively, these results highlight the potential of and boundaries of sequence control via multiblock polymerizations to drive primary sequence ensembles hierarchical structures, and ultimately the functionality of compositionally identical polymeric materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Molecular phylogenetic analysis of the Paguristes tortugae Schmitt, 1933 complex and selected other Paguroidea (Crustacea: Decapoda: Anomura)

Morphological characters, as presently applied to describe members of the Paguristes tortugae Schmitt, 1933 species complex, appear to be of limited value in inferring phylogenetic relationships within the genus, and may have similarly misinformed understanding of relationships between members of this complex and those presently assigned to the related genera Areopaguristes Rahayu & McLaughlin, 2010 and Pseudopaguristes McLaughlin, 2002. Previously undocumented observations of similarities and differences in color patterns among populations additionally suggest genetic divergences within some species, or alternatively seem to support phylogenetic groupings of some species. In the present study, a Maximum Likelihood (ML) phylogenetic analysis was undertaken based on the H3, 12S mtDNA, and 16S mtDNA sequences of 148 individuals, primarily representatives of paguroid species from the western Atlantic. This molecular analysis supported a polyphyletic Diogenidae Ortmann, 1892, although incomplete taxonomic sampling among the genera of Diogenidae limits the utility of this finding for resolving family level relationships. Several hypotheses regarding the evolutionary relationships among hermit crab genera were refuted by the Kishino-Hasegawa (KH). ShimodairaHasegawa (SH) and Approximately Unbiased (AU) tree topology tests, among them the hypothesis that Areopaguristes is monophyletic. A lack of support for the monophyly of Areopaguristes calls into question the phylogenetic validity of gill number for the differentiation of Paguristes, Areopaguristes, and Pseudopaguristes. The study was inconclusive with regard to the relationships among these three genera, but previously unknown diversity within both Paguristes and Areopaguristes was demonstrated. Existence of an undescribed species confounded under the name Paguristes tortugae Schmitt, 1933 was supported by genetics, morphology, and coloration. A second undescribed species with remarkable similarity to Areopaguristes hummi Wass, 1955 was discovered based on genetics and coloration.

12S↗
NASA NTRS2002

Hemichordates and the Origin of Chordates

At the start of the period of the NASA grant three years ago, we had no information on the organization and development of the body axis of the hemichordate, Saccoglossus kowalevskii. Now we have substantial findings about the anteroposterior axis and dorsoventral axis, and based on this information, we have new insights about the origin of chordates from ancestral deuterostomes. We found ways to obtain and preserve large numbers of embryos and hatched juveniles. We can now collect about 40,000 embryos in the month of September, the time of S. kowalevskii spawning at Woods Hole. Excellent cDNA libraries were prepared from three developmental stages. From these libraries, we directly isolated about 30 gene ortholog sequences by screening and pcr techniques, all of these sequences of interest in the inquiry about the animal's organization and development. We also performed a mid-sized EST project (60,000 randomly picked clones, many of these arrayed). About half of these have been analyzed so far by blastx and are suitable for direct use of clones. We have obtained about 50 interesting sequences from this set. The rest still await analysis. Thus, at this time we have isolated orthologs of 80 genes that are known to be expressed in chordates in conserved domains and known to have interesting roles in chordate organization and development. The orthology of the S. kowalevskii sequences has been verified by neighbor joining and parsimony methods, with bootstrap estimates of validity. The S. kowalevskii sequences cluster with other deuterostome sequences, namely, other hemichordates, echinoderms, ascidians, amphioxus, or vertebrates, depending on what sequences are available in the database for comparison. We have used these sequences to do high quality in situ hybridization on S. kowalevskii embryos, and the results can be divided into three sections-those concerning the anteroposterior axis of S. kowalevskii in comparison to the same axis of chordates, those concerning the dorsoventral axis of S. kowalevskii in comparison to the same axis of chordates, and those concerning the signals and transcription factors found in the endoderm, of S. kowalevskii compared to the signals and transcription factors in the endo-mesodermal cells of Spemann's organizer of chordates.

Gerhart, John↗
DOE OSTI2019

Sequenceserver: A Modern Graphical User Interface for Custom BLAST Databases

Comparing newly obtained and previously known nucleotide and amino-acid sequences underpins modern biological research. BLAST is a well-established tool for such comparisons but is challenging to use on new data sets. We combined a user-centric design philosophy with sustainable software development approaches to create Sequenceserver, a tool for running BLAST and visually inspecting BLAST results for biological interpretation. Sequenceserver uses simple algorithms to prevent potential analysis errors and provides flexible text-based and visual outputs to support researcher productivity. Our software can be rapidly installed for use by individuals or on shared servers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Identifying RNA contacts from SHAPE-MaP by partial correlation analysis

In a recent paper Siegfried et al. published a new sequence-based structural RNA assay that utilizes mutational profiling to detect base pairing (MaP). Output from MaP provides information about both pairing (via reactivities) and contact (via correlations). Reactivities can be coupled to partition function folding models for structural inference, while correlations can reveal pairs of sites that may be in structural proximity. The possibility for inference of 3D contacts via MaP suggests a novel approach to structural prediction for RNA analogous to covariance structural prediction for proteins. We explore this approach and show that partial correlation analysis outperforms na\"ive correlation analysis. Our results should be applicable to a wide range of high-throughput sequencing based RNA structural assays that are under development.

Tambe, Akshay↗
DOE OSTI2023

Simulation-Based Recovery Action Analysis Using the EMRALD Dynamic Risk Assessment Tool

A recovery action is defined as the action that prevents deviant conditions from producing unwanted effects. It generally indicates a kind of countermeasure performed in response to a failure of human action. The recovery actions especially play an important role in complex systems like nuclear power plants (NPPs), which consist of highly sophisticated controllers to ensure that desired performance and safety must be achieved and maintained. This is because a combination of human error and its recovery failure may be able to cause a catastrophic effect on a system. Analyzing recovery actions has been a critical part of HRA, which is a technique to evaluate human errors and provide human error probabilities (HEPs) for application in probabilistic safety assessment (PSA). If recovery actions are not adequately analyzed and applied to PSA models, the PSA results may be under-estimated or be not able to reasonably account for the failure of human actions in the context of PSA. For this reason, some regulatory documents such as ASME/ANS RA-Sb-2013 by the American Society for Mechanical Engineers and the American Nuclear Society and NUREG-1792 by U.S. Nuclear Regulatory Commission have emphasized the importance of recovery analysis within the HRA. A couple of existing HRA methods, such as the Technique for Human Error-Rate Prediction (THERP), the Cause-Based Decision Tree (CBDT), and the Korean Standard HRA (K-HRA), have respectively suggested their own approaches to the HRA recovery analysis. However, there are a couple of limitations to treating recovery actions using only the current HRA methods available. The biggest limitation is that the existing recovery analysis does not explicitly consider a variety of recovery action types and recovery sequences as they occur in actual NPPs. To handle the limitations of existing recovery analysis, this study proposes a simulation-based recovery analysis method using the Event Modeling Risk Assessment Using Linked Diagram (EMRALD) software. The EMRALD software is a dynamic simulation tool for PSA. It supports realistic and dynamic modeling of human actions as they would be performed at NPPs. It is also favorable to simultaneously model the specific moment at which an action is performed, the time it takes to perform the action, and the failure probability of that action. In this paper, a detailed methodology for modeling recovery actions in the simulation platform is proposed with a couple of examples. Then, outputs from the simulation are discussed as reviewing if this novel approach can complement the challenges of existing recovery analyses.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

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

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

Touchell, Darren H.↗
DOE OSTI2020

Accurate and complete genomes from metagenomes

Genomes are an integral component of the biological information about an organism; thus, the more complete the genome, the more informative it is. Historically, bacterial and archaeal genomes were reconstructed from pure (monoclonal) cultures, and the first reported sequences were manually curated to completion. However, the bottleneck imposed by the requirement for isolates precluded genomic insights for the vast majority of microbial life. Shotgun sequencing of microbial communities, referred to initially as community genomics and subsequently as genome-resolved metagenomics, can circumvent this limitation by obtaining metagenome-assembled genomes (MAGs); but gaps, local assembly errors, chimeras, and contamination by fragments from other genomes limit the value of these genomes. Here, we discuss genome curation to improve and, in some cases, achieve complete (circularized, no gaps) MAGs (CMAGs). To date, few CMAGs have been generated, although notably some are from very complex systems such as soil and sediment. Through analysis of about 7000 published complete bacterial isolate genomes, we verify the value of cumulative GC skew in combination with other metrics to establish bacterial genome sequence accuracy. The analysis of cumulative GC skew identified potential misassemblies in some reference genomes of isolated bacteria and the repeat sequences that likely gave rise to them. We discuss methods that could be implemented in bioinformatic approaches for curation to ensure that metabolic and evolutionary analyses can be based on very high-quality genomes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

KBase Narrative - Supplemental Data for "Mixed waste contamination selects for a mobile genetic element population enriched in multiple heavy metal resistance genes"

Here we provide the complete set of circular elements produced by SCAPP from metagenomic data that may represent mobile genetic elements (MGEs). For the publication, these sequences were de-replicated within the high [U] and low [U] sample sets. Additionally, 98 of the circular elements were removed for downstream analysis due to high suspicion of being (1) chloroplast or mitochondrial sequences or (2) erroneously circularized long repeat regions. The complete list of circular elements used for the analysis along with metadata can be found in Tables S5 and S6.

Goff, Jennifer↗
DOE OSTI2020

Single-Cell Genomics of Novel Actinobacteria With the Wood–Ljungdahl Pathway Discovered in a Serpentinizing System

Serpentinite-hosted systems represent modern-day analogs of early Earth environments. In these systems, water-rock interactions generate highly alkaline and reducing fluids that can contain hydrogen, methane, and low-molecular-weight hydrocarbons-potent reductants capable of fueling microbial metabolism. In this study, we investigated the microbiota of Hakuba Happo hot springs (~50°C; pH~10.5–11), located in Nagano (Japan), which are impacted by the serpentinization process. Analysis of the 16S rRNA gene amplicon sequences revealed that the bacterial community comprises Nitrospirae (47%), “Parcubacteria” (19%), Deinococcus-Thermus (16%), and Actinobacteria (9%), among others. Notably, only 57 amplicon sequence variants (ASV) were detected, and fifteen of these accounted for 90% of the amplicons. Among the abundant ASVs, an early-branching, uncultivated actinobacterial clade identified as RBG-16-55-12 in the SILVA database was detected. Ten single-cell genomes (average pairwise nucleotide identity: 0.98–1.00; estimated completeness: 33–93%; estimated genome size: ~2.3 Mb) that affiliated with this clade were obtained. Taxonomic classification using single copy genes indicates that the genomes belong to the actinobacterial class-level clade UBA1414 in the Genome Taxonomy Database. Based on metabolic pathway predictions, these actinobacteria are anaerobes, capable of glycolysis, dissimilatory nitrate reduction and CO 2 fixation via the Wood–Ljungdahl (WL) pathway. Several other genomes within UBA1414 and two related class-level clades also encode the WL pathway, which has not yet been reported for the Actinobacteria phylum. For the Hakuba actinobacterium, the energy metabolism related to the WL pathway is likely supported by a combination of the Rnf complex, group 3b and 3d [NiFe]-hydrogenases, [FeFe]-hydrogenases, and V-type (H + /Na + pump) ATPase. The genomes also harbor a form IV ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) complex, also known as a RubisCO-like protein, and contain signatures of interactions with viruses, including clustered regularly interspaced short palindromic repeat (CRISPR) regions and several phage integrases. This is the first report and detailed genome analysis of a bacterium within the Actinobacteria phylum capable of utilizing the WL pathway. The Hakuba actinobacterium is a member of the clade UBA1414/RBG-16-55-12, formerly within the group “OPB41.” We propose to name this bacterium ‘Candidatus Hakubanella thermoalkaliphilus.'

59 BASIC BIOLOGICAL SCIENCES↗