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 775 records · Page 43

Detailed analysis of putative genes encoding small proteins in legume genomes

Diverse plant genome sequencing projects coupled with powerful bioinformatics tools have facilitated massive data analysis to construct specialized databases classified according to cellular function. However, there are still a considerable number of genes encoding proteins whose function has not yet been characterized. Included in this category are small proteins (SPs, 30–150 amino acids) encoded by short open reading frames (sORFs). SPs play important roles in plant physiology, growth, and development. Unfortunately, protocols focused on the genome-wide identification and characterization of sORFs are scarce or remain poorly implemented. As a result, these genes are underrepresented in many genome annotations. In this work, we exploited publicly available genome sequences of Phaseolus vulgaris, Medicago truncatula, Glycine max, and Lotus japonicus to analyze the abundance of annotated SPs in plant legumes. Our strategy to uncover bona fide sORFs at the genome level was centered in bioinformatics analysis of characteristics such as evidence of expression (transcription), presence of known protein regions or domains, and identification of orthologous genes in the genomes explored. We collected 6170, 10,461, 30,521, and 23,599 putative sORFs from P. vulgaris, G. max, M. truncatula, and L. japonicus genomes, respectively. Expressed sequence tags (ESTs) available in the DFCI Gene Index database provided evidence that ~one-third of the predicted legume sORFs are expressed. Most potential SPs have a counterpart in a different plant species and counterpart regions or domains in larger proteins. Potential functional sORFs were also classified according to a reduced set of GO categories, and the expression of 13 of them during P. vulgaris nodule ontogeny was confirmed by qPCR. This analysis provides a collection of sORFs that potentially encode for meaningful SPs, and offers the possibility of their further functional evaluation.

59 BASIC BIOLOGICAL SCIENCES↗

Identification of new marker genes from plant single‐cell RNA‐seq data using interpretable machine learning methods

Summary An essential step in the analysis of single‐cell RNA sequencing data is to classify cells into specific cell types using marker genes. In this study, we have developed a machine learning pipeline called single‐cell predictive marker (SPmarker) to identify novel cell‐type marker genes in the Arabidopsis root. Unlike traditional approaches, our method uses interpretable machine learning models to select marker genes. We have demonstrated that our method can: assign cell types based on cells that were labelled using published methods; project cell types identified by trajectory analysis from one data set to other data sets; and assign cell types based on internal GFP markers. Using SPmarker, we have identified hundreds of new marker genes that were not identified before. As compared to known marker genes, the new marker genes have more orthologous genes identifiable in the corresponding rice single‐cell clusters. The new root hair marker genes also include 172 genes with orthologs expressed in root hair cells in five non‐ Arabidopsis species, which expands the number of marker genes for this cell type by 35–154%. Our results represent a new approach to identifying cell‐type marker genes from scRNA‐seq data and pave the way for cross‐species mapping of scRNA‐seq data in plants.

54 ENVIRONMENTAL SCIENCES↗

Description of Pegethrix niliensis sp. nov., a Novel Cyanobacterium from the Nile River Basin, Egypt: A Polyphasic Analysis and Comparative Study of Related Genera in the Oculatellales Order

In this paper, we examine the filamentous cyanobacterial strain NILCB16 and describe it as a new species within the genus Pegethrix. The original population was sampled from a mat growing in an irrigation canal in the Nile River, Egypt. Initially classified under Plectonema or Planktolyngbya, the strain is a potential producer of the toxins microcystin and β-N-Methylamino-L-Alanine (BMAA). Additionally, we reviewed the taxonomic relationships between the Oculatellales genera. To describe the new species, we conducted a polyphasic study, encompassing 16S rRNA gene phylogenetic analyses performed using both Maximum Likelihood and Bayesian methods, sequence identity (p-distance) analysis, 16S-23S ITS secondary structures, and morphological and habitat comparisons. The phylogenetic analysis revealed that strain NILCB16 clustered within the Pegethrix clade with strong phylogenetic support, but in a distinct position from other species in the genus. The strain shared a maximum 16S rRNA gene identity of 97.3% with P. qiandaoensis and 96.1% with the type species, P. bostrychoides. Morphologically, NILCB16 can be differentiated from other species in the genus by its lack of false branching. Our phylogenetic analyses also show that Pegethrix, Cartusia, Elainella, and Maricoleus are clustered with strong phylogenetic support. They exhibit high 16S rRNA gene identity and are morphologically indistinguishable, suggesting they could potentially be merged into a single genus in the future.

Hentschke, Guilherme Scotta (ORCID:000000034396024↗

STELLAR AGES AND CONVECTIVE CORES IN FIELD MAIN-SEQUENCE STARS: FIRST ASTEROSEISMIC APPLICATION TO TWO KEPLER TARGETS

Using asteroseismic data and stellar evolution models we obtain the first detection of a convective core in a Kepler field main-sequence star, putting a stringent constraint on the total size of the mixed zone and showing that extra mixing beyond the formal convective boundary exists. In a slightly less massive target the presence of a convective core cannot be conclusively discarded, and thus its remaining main-sequence lifetime is uncertain. Our results reveal that best-fit models found solely by matching individual frequencies of oscillations corrected for surface effects do not always properly reproduce frequency combinations. Moreover, slightly different criteria to define what the best-fit model is can lead to solutions with similar global properties but very different interior structures. We argue that the use of frequency ratios is a more reliable way to obtain accurate stellar parameters, and show that our analysis in field main-sequence stars can yield an overall precision of 1.5%, 4%, and 10% in radius, mass, and age, respectively. We compare our results with those obtained from global oscillation properties, and discuss the possible sources of uncertainties in asteroseismic stellar modeling where further studies are still needed.

79 ASTRONOMY AND ASTROPHYSICS↗

Syntgen: a system to generate temporal networks with user-specified topology

Abstract In the last few years, the study of temporal networks has progressed markedly. The evolution of clusters of nodes (or communities) is one of the major focus of these studies. However, the time dimension increases complexity, introducing new constructs and requiring novel and enhanced algorithms. In spite of recent improvements, the relative scarcity of timestamped representations of empiric networks, with known ground truth, hinders algorithm validation. A few approaches have been proposed to generate synthetic temporal networks that conform to static topological specifications while in general adopting an ad hoc approach to temporal evolution. We believe there is still a need for a principled synthetic network generator that conforms to problem domain topological specifications from a static as well as temporal perspective. Here, we present such a system. The unique attributes of our system include accepting arbitrary node degree and cluster size distributions and temporal evolution under user control, while supporting tunable joint distribution and temporal correlation of node degrees. Theoretical contributions include the analysis of conditions for graphic sequences of inter- and intracluster node degrees and cluster sizes and the development of a heuristic to search for the cluster membership of nodes that minimizes the shared information distance between clusterings. Our work shows that this system is capable of generating networks under user controlled topology with up to thousands of nodes and hundreds of clusters with strong topology adherence. Much larger networks are possible with relaxed requirements. The generated networks support algorithm validation as well as problem domain analysis.

Pereira, Luis Ramada↗

Complete genome sequence of Anabaena variabilis ATCC 29413

Anabaena variabilis ATCC 29413 is a filamentous, heterocyst-forming cyanobacterium that has served as a model organism, with an extensive literature extending over 40 years. The strain has three distinct nitrogenases that function under different environmental conditions and is capable of photoautotrophic growth in the light and true heterotrophic growth in the dark using fructose as both carbon and energy source. While this strain was first isolated in 1964 in Mississippi and named Anabaena flos-aquae MSU A-37, it clusters phylogenetically with cyanobacteria of the genus Nostoc. The strain is a moderate thermophile, growing well at approximately 40° C. Here we provide some additional characteristics of the strain, and an analysis of the complete genome sequence

59 BASIC BIOLOGICAL SCIENCES↗

Complete genome sequence of Anabaena variabilis ATCC 29413

Anabaena variabilis ATCC 29413 is a filamentous, heterocyst-forming cyanobacterium that has served as a model organism, with an extensive literature extending over 40 years. The strain has three distinct nitrogenases that function under different environmental conditions and is capable of photoautotrophic growth in the light and true heterotrophic growth in the dark using fructose as both carbon and energy source. While this strain was first isolated in 1964 in Mississippi and named Ana-baena flos-aquae MSU A-37, it clusters phylogenetically with cyanobacteria of the genus Nostoc. The strain is a moderate thermophile, growing well at approximately 40 C. Here we provide some additional characteristics of the strain, and an analysis of the complete genome sequence.

Thiel, Teresa↗

Protocol for applying a network-enabled gene discovery pipeline to non-model plant species

Identifying upstream regulators of key genes is essential for understanding gene regulatory mechanisms and translating these insights into functional targets. Here, we present a protocol for applying the network-enabled gene discovery pipeline (NEEDLE) to non-model plant species. We describe steps for environment setup, data preparation, computational analysis, expected outputs, and parameter considerations. NEEDLE integrates RNA sequencing (RNA-seq) processing, weighted gene co-expression analysis (WGCNA), Gene Network Inference with Ensemble of trees (GENIE3), and promoter conservation analysis to prioritize candidate transcriptional regulators.

Plant Sciences↗

Corn Storage Protein - A Molecular Genetic Model

Corn is the highest yielding crop on earth and probably the most valuable agricultural product of the United States. Because it converts sun energy through photosynthesis into starch and proteins, we addressed energy savings by focusing on protein quality. People and animals require essential amino acids derived from the digestion of proteins. If proteins are relatively low in certain essential amino acids, the crop becomes nutritionally defective and has to be supplemented. Such deficiency affects meat and fish production and countries where corn is a staple. Because corn seed proteins have relatively low levels of lysine and methionine, a diet has to be supplemented with soybeans for the missing lysine and with chemically synthesized methionine. We therefore have studied genes expressed during maize seed development and their chromosomal organization. A critical technical requirement for the understanding of the molecular structure of genes and their positional information was DNA sequencing. Because of the length of sequences, DNA sequencing methods themselves were insufficient for this type of analysis. We therefore developed the so-called “DNA shotgun sequencing” strategy, where overlapping DNA fragments were sequenced in parallel and used to reconstruct large DNA molecules via overlaps. Our publications became the most frequently cited ones during the decade of 1981-1990 and former Associate Director of Science for the Office of Basic Energy Sciences Patricia M. Dehmer presented our work as one of the great successes of this program. A major component of the sequencing strategy was the development of bacterial strains and vectors, which were also used to develop the first biotechnology crops. These crops possessed new traits thanks to the expression of foreign genes in plants. To enable such expression, chimeric genes had to be constructed using our materials and methods by the industry. Because we made our materials and methods freely available to academia and industry, progress in plant research and new crop development could accelerate and benefit the public.

60 APPLIED LIFE SCIENCES↗

Archaeal and bacterial communities across a chronosequence of drained lake basins in arctic alaska

We examined patterns in soil microbial community composition across a successional gradient of drained lake basins in the Arctic Coastal Plain. Analysis of 16S rRNA gene sequences revealed that methanogens closely related to Candidatus ‘Methanoflorens stordalenmirensis’ were the dominant archaea, comprising >50% of the total archaea at most sites, with particularly high levels in the oldest basins and in the top 57 cm of soil (active and transition layers). Bacterial community composition was more diverse, with lineages from OP11, Actinobacteria, Bacteroidetes, and Proteobacteria found in high relative abundance across all sites. Notably, microbial composition appeared to converge in the active layer, but transition and permafrost layer communities across the sites were significantly different to one another. Microbial biomass using fatty acid-based analysis indicated that the youngest basins had increased abundances of gram-positive bacteria and saprotrophic fungi at higher soil organic carbon levels, while the oldest basins displayed an increase in only the gram-positive bacteria. While this study showed differences in microbial populations across the sites relevant to basin age, the dominance of Candidatus ‘M. stordalenmirensis’ across the chronosequence indicates the potential for changes in local carbon cycling, depending on how these methanogens and associated microbial communities respond to warming temperatures.

59 BASIC BIOLOGICAL SCIENCES↗

Diversity and Structure of Diazotrophic Communities in Mangrove Rhizosphere, Revealed by High-Throughput Sequencing

Diazotrophic communities make an essential contribution to the productivity through providing new nitrogen. However, knowledge of the roles that both mangrove tree species and geochemical parameters play in shaping mangove rhizosphere diazotrophic communities is still elusive. Here, a comprehensive examination of the diversity and structure of microbial communities in the rhizospheres of three mangrove species, Rhizophora apiculata , Avicennia marina , and Ceriops tagal , was undertaken using high-throughput sequencing of the 16S rRNA and nifH genes. Our results revealed a great diversity of both the total microbial composition and the diazotrophic composition specifically in the mangrove rhizosphere. Deltaproteobacteria and Gammaproteobacteria were both ubiquitous and dominant, comprising an average of 45.87 and 86.66% of total microbial and diazotrophic communities, respectively. Sulfate-reducing bacteria belonging to the Desulfobacteraceae and Desulfovibrionaceae were the dominant diazotrophs. Community statistical analyses suggested that both mangrove tree species and additional environmental variables played important roles in shaping total microbial and potential diazotroph communities in mangrove rhizospheres. In contrast to the total microbial community investigated by analysis of 16S rRNA gene sequences, most of the dominant diazotrophic groups identified by nifH gene sequences were significantly different among mangrove species. The dominant diazotrophs of the family Desulfobacteraceae were positively correlated with total phosphorus, but negatively correlated with the nitrogen to phosphorus ratio. Finally, the Pseudomonadaceae were positively correlated with the concentration of available potassium, suggesting that diazotrophs potentially play an important role in biogeochemical cycles, such as those of nitrogen, phosphorus, sulfur, and potassium, in the mangrove ecosystem.

59 BASIC BIOLOGICAL SCIENCES↗

Standardized phylogenetic and molecular evolutionary analysis applied to species across the microbial tree of life

There is growing interest in reconstructing phylogenies from the copious amounts of genome sequencing projects that target related viral, bacterial or eukaryotic organisms. To facilitate the construction of standardized and robust phylogenies for disparate types of projects, we have developed a complete bioinformatic workflow, with a web-based component to perform phylogenetic and molecular evolutionary (PhaME) analysis from sequencing reads, draft assemblies or completed genomes of closely related organisms. Furthermore, the ability to incorporate raw data, including some metagenomic samples containing a target organism (e.g. from clinical samples with suspected infectious agents), shows promise for the rapid phylogenetic characterization of organisms within complex samples without the need for prior assembly.

59 BASIC BIOLOGICAL SCIENCES↗

A draft genome of field pennycress ( Thlaspi arvense ) provides tools for the domestication of a new winter biofuel crop

Field pennycress (Thlaspi arvense L.) is being domesticated as a new winter cover crop and biofuel species for the Midwestern United States that can be double-cropped between corn and soybeans. A genome sequence will enable the use of new technologies to make improvements in pennycress. To generate a draft genome, a hybrid sequencing approach was used to generate 47 Gb of DNA sequencing reads from both the Illumina and PacBio platforms. These reads were used to assemble 6,768 genomic scaffolds. The draft genome was annotated using the MAKER pipeline, which identified 27,390 predicted protein-coding genes, with almost all of these predicted peptides having significant sequence similarity to Arabidopsis proteins. A comprehensive analysis of pennycress gene homologues involved in glucosinolate biosynthesis, metabolism, and transport pathways revealed high sequence conservation compared with other Brassicaceae species, and helps validate the assembly of the pennycress gene space in this draft genome. Additional comparative genomic analyses indicate that the knowledge gained from years of basic Brassicaceae research will serve as a powerful tool for identifying gene targets whose manipulation can be predicted to result in improvements for pennycress.

09 BIOMASS FUELS↗

Revisiting N 2 fixation in Guerrero Negro intertidal microbial mats with a functional single-cell approach

Photosynthetic microbial mats are complex, stratified ecosystems in which high rates of primary production create a demand for nitrogen, met partially by N 2 fixation. Dinitrogenase reductase ( nifH ) genes and transcripts from Cyanobacteria and heterotrophic bacteria (for example, Deltaproteobacteria ) were detected in these mats, yet their contribution to N 2 fixation is poorly understood. We used a combined approach of manipulation experiments with inhibitors, nifH sequencing and single-cell isotope analysis to investigate the active diazotrophic community in intertidal microbial mats at Laguna Ojo de Liebre near Guerrero Negro, Mexico. Acetylene reduction assays with specific metabolic inhibitors suggested that both sulfate reducers and members of the Cyanobacteria contributed to N 2 fixation, whereas 15 N 2 tracer experiments at the bulk level only supported a contribution of Cyanobacteria . Cyanobacterial and nifH Cluster III (including deltaproteobacterial sulfate reducers) sequences dominated the nifH gene pool, whereas the nifH transcript pool was dominated by sequences related to Lyngbya spp. Single-cell isotope analysis of 15 N 2 -incubated mat samples via high-resolution secondary ion mass spectrometry (NanoSIMS) revealed that Cyanobacteria were enriched in 15 N, with the highest enrichment being detected in Lyngbya spp. filaments (on average 4.4 at% 15 N), whereas the Deltaproteobacteria (identified by CARD-FISH) were not significantly enriched. We investigated the potential dilution effect from CARD-FISH on the isotopic composition and concluded that the dilution bias was not substantial enough to influence our conclusions. As a result, our combined data provide evidence that members of the Cyanobacteria , especially Lyngbya spp., actively contributed to N 2 fixation in the intertidal mats, whereas support for significant N 2 fixation activity of the targeted deltaproteobacterial sulfate reducers could not be found.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

REFLECTION OF PROPAGATING SLOW MAGNETO-ACOUSTIC WAVES IN HOT CORONAL LOOPS: MULTI-INSTRUMENT OBSERVATIONS AND NUMERICAL MODELING

Slow MHD waves are important tools for understanding coronal structures and dynamics. In this paper, we report a number of observations from the X-Ray Telescope (XRT) on board HINODE and Solar Dynamic Observatory /Atmospheric Imaging Assembly (AIA) of reflecting longitudinal waves in hot coronal loops. To our knowledge, this is the first report of this kind as seen from the XRT and simultaneously with the AIA. The wave appears after a micro-flare occurs at one of the footpoints. We estimate the density and temperature of the loop plasma by performing differential emission measure (DEM) analysis on the AIA image sequence. The estimated speed of propagation is comparable to or lower than the local sound speed, suggesting it to be a propagating slow wave. The intensity perturbation amplitude, in every case, falls very rapidly as the perturbation moves along the loop and eventually vanishes after one or more reflections. To check the consistency of such reflection signatures with the obtained loop parameters, we perform a 2.5D MHD simulation, which uses the parameters obtained from our observation as inputs, and perform forward modeling to synthesize AIA 94 Å images. Analyzing the synthesized images, we obtain the same properties of the observables as for the real observation. From the analysis we conclude that a footpoint heating can generate a slow wave which then reflects back and forth in the coronal loop before fading. Our analysis of the simulated data shows that the main agent for this damping is anisotropic thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗

Prevalence of gp160 polymorphisms known to be related to decreased susceptibility to temsavir in different subtypes of HIV-1 in the Los Alamos National Laboratory HIV Sequence Database

Fostemsavir, a prodrug of the gp120-directed attachment inhibitor temsavir, is indicated for use in heavily treatment-experienced individuals with MDR HIV-1. Reduced susceptibility to temsavir in the clinic maps to discrete changes at amino acid positions in gp160: S375, M426, M434 and M475.To query the Los Alamos National Laboratory (LANL) HIV Sequence Database for the prevalence of polymorphisms at gp160 positions of interest. Full-length gp160 sequences (N = 7560) were queried for amino acid polymorphisms relative to the subtype B consensus at positions of interest; frequencies were reported for all sequences and among subtypes/circulating recombinant forms (CRFs) with ≥10 isolates in the database. Among 239 subtypes in the database, the 5 most prevalent were B (n = 2651, 35.1%), C (n = 1626, 21.5%), CRF01_AE (n = 674, 8.9%), A1 (n = 273, 3.6%) and CRF02_AG (n = 199, 2.6%). Among all 7560 sequences, the most prevalent amino acids at positions of interest (S 375 , 73.5%; M 426 , 82.1%; M 434 , 88.2%; M 475 , 89.9%) were the same as the subtype B consensus. Specific polymorphisms with the potential to decrease temsavir susceptibility (S 375 H/I/M/N/T/Y, M 426 L/P, M 434 I/K and M475I) were found in <10% of isolates of subtypes D, G, A6, BC, F1, CRF07_BC, CRF08_BC, 02A, CRF06_cpx, F2, 02G and 02B. S 375 H and M 475 I were predominant among CRF01_AE (S375H, 99.3%; M 475 I, 76.3%; consistent with previously reported low temsavir susceptibility of this CRF) and 01B (S 375 H, 71.7%; M 475 I, 49.5%). Analysis of the LANL HIV Sequence Database found a low prevalence of gp160 amino acid polymorphisms with the potential to reduce temsavir susceptibility overall and among most of the common subtypes.

59 BASIC BIOLOGICAL SCIENCES↗

Risk-Informed Safety Margin Characterization Case Study: Selection of Electrical Equipment to Be Subjected to Environmental Qualification

The Risk-Informed Safety Margin Characterization (RISMC) pathway of the DOE’s Light Water Reactor Sustainability (LWRS) program focuses on advancing the state of the art in safety analysis and risk assessment to support decision-making on nuclear power plant operation well beyond the originally designed lifetime of the plants (i.e., beyond 60 years). Among the issues being addressed in RISMC is the significance of SSC aging and how confident we are about our understanding of its impact on the margin between the loads SSCs are expected to see during normal operation and accident conditions, and the SSC capacities (their ability to resist those loads) as the SSCs age. In this paper, a summary is provided of a case study that examines SSC aging from an environmental qualification (EQ) perspective. The case study illustrates how the state of knowledge regarding SSC margin can be characterized given the overall integrated plant design, and was developed to demonstrate a method for deciding on which cables to focus, which cables are not so important from an environmental qualification margin standpoint, and what plant design features or operating characteristics determine the role that environmental qualification plays in establishing a safety case on which decisions regarding margin can be made. The selection of cables for which demonstration of margin with respect to aging and environmental challenges uses a technique known as Prevention Analysis. Prevention Analysis is a Boolean method for optimal selection of SSCs (that is, those combinations of SSCs both necessary and sufficient to meet a predetermined selection criterion) in a manner that allows demonstration that plant-level safety can be demonstrated by the collection of selected SSCs alone. Choosing the set of SSCs that is necessary and sufficient to satisfy the safety objectives, and demonstrating that the safety objectives can be met effectively, determines where resources are best allocated to assure SSC performance margin. The paper describes the resulting component types that were selected by Prevention Analysis and identifies the accident sequence characteristics that cause these component types to be important from an EQ and aging perspective (and, hence, worthwhile evaluating the extent of safety margin). In addition, component types not selected as needing significant margin from an EQ and aging perspective are discussed and an engineering rationale is developed justifying the lack of need to apply resources to demonstrating margin for these component types. This rationale is in terms of design features of the plant and operating characteristics that make these component types less important from an EQ and aging perspective. While the case study focuses on EQ and aging of equipment and cables located inside the containment of this PWR, the prevention analysis method is demonstrated to be an effective technique for identification of minimal collections of components that would be effective in managing safety for a variety of issues associated with aging and long-term operation of the fleet of plants.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Machine-Learning-Based Analysis of the Guy-Greenbrier, Arkansas Earthquakes: A Tale of Two Sequences

Here, we revisited the June 2010 to October 2011 Guy-Greenbrier earthquake sequence in central Arkansas using PhaseNet, a deep neural network trained to pick P and S arrival times. We applied PhaseNet to continuous waveform data and used phase association and hypocenter relocation to locate nearly 90,000 events. Our catalog suggests that the sequence consists of two adjacent earthquake sequences on the same fault and that the second sequence may be associated with the wastewater disposal well to the west of the Guy-Greenbrier Fault, rather than the wells to the north and the east that were previously implicated. We find that each sequence is composed of many small clusters that exhibit diffusion along the fault at shorter timescales. Our study demonstrates that machine-learning-based earthquake catalog development is now feasible and will yield new insights into earthquake behavior.

58 GEOSCIENCES↗