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 127 records · Page 7

Genome-wide survey of Aux/IAA gene family members in potato (Solanum tuberosum): Identification, expression analysis, and evaluation of their roles in tuber development

The Auxin/indole-3-acetic acid (Aux/IAA) genes encode short-lived nuclear proteins that are known to be involved in the primary cellular responses to auxin. To date, systematic analysis of the Aux/IAA genes in potato (Solanum tuberosum) has not been conducted. In this study, a total of 26 potato Aux/IAA genes were identified (designated from StIAA1 to StIAA26), and the distribution of four conserved domains shared by the StIAAs were analyzed based on multiple sequence alignment and a motif-based sequence analysis. A phylogenetic analysis of the Aux/IAA gene families of potato and Arabidopsis was also conducted. In order to assess the roles of StIAA genes in tuber development, the results of RNA-seq studies were reformatted to analyze the expression patterns of StIAA genes, and then verified by quantitative real-time PCR. A large number of StIAA genes (12 genes) were highly expressed in stolon organs and in during the tuber initiation and expansion developmental stages, and most of these genes were responsive to indoleacetic acid treatment. Our results suggested that StIAA genes were involved in the process of tuber development and provided insights into functional roles of potato Aux/IAA genes. - Highlights: • A systematic analysis of the potato AUX/IAA gene family were performed. • StIAA genes were related to auxin perception and signal transduction. • Candidate StIAA genes likely related to tuber initiation and expansion were screened.

60 APPLIED LIFE SCIENCES↗

Comparative genomics analysis of drought response between obligate CAM and C 3 photosynthesis plants

Crassulacean acid metabolism (CAM) plants exhibit elevated drought and heat tolerance compared to C 3 and C 4 plants through an inverted pattern of day/night stomatal closure and opening for CO 2 assimilation. However, the molecular responses to water-deficit conditions remain unclear in obligate CAM species. In this study, we presented genome-wide transcription sequencing analysis using leaf samples of an obligate CAM species Kalanchoë fedtschenkoi under moderate and severe drought treatments at two-time points of dawn (2-h before the start of light period) and dusk (2-h before the dark period). Differentially expressed genes were identified in response to environmental drought stress and a whole genome wide co-expression network was created as well. We found that the expression of CAM-related genes was not regulated by drought stimuli in K. fedtschenkoi. Our comparative analysis revealed that CAM species (K. fedtschenkoi) and C 3 species (Arabidopsis thaliana, Populus deltoides ‘WV94’) share some common transcriptional changes in genes involved in multiple biological processes in response to drought stress, including ABA signaling and biosynthesis of secondary metabolites.

59 BASIC BIOLOGICAL SCIENCES↗

Evolution of plant virus movement proteins from the 30K superfamily and of their homologs integrated in plant genomes

Homologs of Tobacco mosaic virus 30K cell-to-cell movement protein are encoded by diverse plant viruses. Mechanisms of action and evolutionary origins of these proteins remain obscure. We expand the picture of conservation and evolution of the 30K proteins, producing sequence alignment of the 30K superfamily with the broadest phylogenetic coverage thus far and illuminating structural features of the core all-beta fold of these proteins. Integrated copies of pararetrovirus 30K movement genes are prevalent in euphyllophytes, with at least one copy intact in nearly every examined species, and mRNAs detected for most of them. Sequence analysis suggests repeated integrations, pseudogenizations, and positive selection in those provirus genes. An unannotated 30K-superfamily gene in Arabidopsis thaliana genome is likely expressed as a fusion with the At1g37113 transcript. This molecular background of endopararetrovirus gene products in plants may change our view of virus infection and pathogenesis, and perhaps of cellular homeostasis in the hosts. - Highlights: • Sequence region shared by plant virus “30K” movement proteins has an all-beta fold. • Most euphyllophyte genomes contain integrated copies of pararetroviruses. • These integrated virus genomes often include intact movement protein genes. • Molecular evidence suggests that these “30K” genes may be selected for function.

60 APPLIED LIFE SCIENCES↗

Molecular Analysis of Asymptomatic Bacteriuria Escherichia coli Strain VR50 Reveals Adaptation to the Urinary Tract by Gene Acquisition

Urinary tract infections (UTIs) are among the most common infectious diseases of humans, with Escherichia coli for >80% of all cases. One extreme of UTI is asymptomatic bacteriuria (ABU), which occurs as an asymptomatic carrier state that resembles commensalism. Here, to understand the evolution and molecular mechanisms that underpin ABU, the genome of the ABU E. coli strain VR50 was sequenced. Analysis of the complete genome indicated that it most resembles E. coli K-12, with the addition of a 94-kb genomic island (GI-VR50-pheV), eight prophages, and multiple plasmids. GI-VR50- pheV has a mosaic structure and contains genes encoding a number of UTI-associated virulence factors, namely, Afa (afimbrial adhesin), two autotransporter proteins (Ag43 and Sat), and aerobactin. We demonstrated that the presence of this island in VR50 confers its ability to colonize the murine bladder, as a VR50 mutant with GI-VR50- pheV deleted was attenuated in a mouse model of UTI in vivo . We established that Afa is the island-encoded factor responsible for this phenotype using two independent deletion (Afa operon and AfaE adhesin) mutants. E. coli VR50 afa and VR50 afaE displayed significantly decreased ability to adhere to human bladder epithelial cells. In the mouse model of UTI, VR50 afa and VR50 afaE displayed reduced bladder colonization compared to wild-type VR50, similar to the colonization level of the GI-VR50- pheV mutant. In conlusion, our study suggests that E. coli VR50 is a commensal-like strain that has acquired fitness factors that facilitate colonization of the human bladder.

59 BASIC BIOLOGICAL SCIENCES↗

Identification of a Latin American-specific BabA adhesin variant through whole genome sequencing of Helicobacter pylori patient isolates from Nicaragua

In this study, Helicobacter pylori (H. pylori) is one of the most common bacterial infections in humans and this infection can lead to gastric ulcers and gastric cancer. H. pylori is one of the most genetically variable human pathogens and the ability of the bacterium to bind to the host epithelium as well as the presence of different virulence factors and genetic variants within these genes have been associated with disease severity. Nicaragua has particularly high gastric cancer incidence and we therefore studied Nicaraguan clinical H. pylori isolates for factors that could contribute to cancer risk. The complete genomes of fifty-two Nicaraguan H. pylori i isolates were sequenced and assembled de novo, and phylogenetic and virulence factor analyses were performed. The Nicaraguan isolates showed phylogenetic relationship with West African isolates in whole-genome sequence comparisons and with Western and urban South-and Central American isolates using MLSA (Multi-locus sequence analysis). A majority, 77 % of the isolates carried the cancer-associated virulence gene cagA and also the s1/i1/m1 vacuolating cytotoxin, vacA allele combination, which is linked to increased severity of disease. Specifically, we also found that Nicaraguan isolates have a blood group-binding adhesin (BabA) variant highly similar to previously reported BabA sequences from Latin America, including from isolates belonging to other phylogenetic groups. These BabA sequences were found to be under positive selection at several amino acid positions that differed from the global collection of isolates. In conclusion, the discovery of a Latin American BabA variant, independent of overall phylogenetic background, suggests hitherto unknown host or environmental factors within the Latin American population giving H. pylori isolates carrying this adhesin variant a selective advantage, which could affect pathogenesis and risk for sequelae through specific adherence properties.

59 BASIC BIOLOGICAL SCIENCES↗

Deep functional analysis of synII, a 770-kilobase synthetic yeast chromosome

INTRODUCTION Although much effort has been devoted to studying yeast in the past few decades, our understanding of this model organism is still limited. Rapidly developing DNA synthesis techniques have made a “build-to-understand” approach feasible to reengineer on the genome scale. Here, we report on the completion of a 770-kilobase synthetic yeast chromosome II (synII). SynII was characterized using extensive Trans-Omics tests. Despite considerable sequence alterations, synII is virtually indistinguishable from wild type. However, an up-regulation of translational machinery was observed and can be reversed by restoring the transfer RNA (tRNA) gene copy number. RATIONALE Following the “design-build-test-debug” working loop, synII was successfully designed and constructed in vivo. Extensive Trans-Omics tests were conducted, including phenomics, transcriptomics, proteomics, metabolomics, chromosome segregation, and replication analyses. By both complementation assays and SCRaMbLE (synthetic chromosome rearrangement and modification by loxP -mediated evolution), we targeted and debugged the origin of a growth defect at 37°C in glycerol medium. RESULTS To efficiently construct megabase-long chromosomes, we developed an I- Sce I–mediated strategy, which enables parallel integration of synthetic chromosome arms and reduced the overall integration time by 50% for synII. An I- Sce I site is introduced for generating a double-strand break to promote targeted homologous recombination during mitotic growth. Despite hundreds of modifications introduced, there are still regions sharing substantial sequence similarity that might lead to undesirable meiotic recombinations when intercrossing the two semisynthetic chromosome arm strains. Induction of the I- Sce I–mediated double-strand break is otherwise lethal and thus introduced a strong selective pressure for targeted homologous recombination. Since our strategy is designed to generate a markerless synII and leave the URA3 marker on the wild-type chromosome, we observed a tenfold increase in URA3 -deficient colonies upon I- Sce I induction, meaning that our strategy can greatly bias the crossover events toward the designated regions. By incorporating comprehensive phenotyping approaches at multiple levels, we demonstrated that synII was capable of powering the growth of yeast indistinguishably from wild-type cells (see the figure), showing highly consistent biological processes comparable to the native strain. Meanwhile, we also noticed modest but potentially significant up-regulation of the translational machinery. The main alteration underlying this change in expression is the deletion of 13 tRNA genes. A growth defect was observed in one very specific condition—high temperature (37°C) in medium with glycerol as a carbon source—where colony size was reduced significantly. We targeted and debugged this defect by two distinct approaches. The first approach involved phenotype screening of all intermediate strains followed by a complementation assay with wild-type sequences in the synthetic strain. By doing so, we identified a modification resulting from PCRTag recoding in TSC10 , which is involved in regulation of the yeast high-osmolarity glycerol (HOG) response pathway. After replacement with wild-type TSC10 , the defect was greatly mitigated. The other approach, debugging by SCRaMbLE, showed rearrangements in regions containing HOG regulation genes. Both approaches indicated that the defect is related to HOG response dysregulation. Thus, the phenotypic defect can be pinpointed and debugged through multiple alternative routes in the complex cellular interactome network. CONCLUSION We have demonstrated that synII segregates, replicates, and functions in a highly similar fashion compared with its wild-type counterpart. Furthermore, we believe that the iterative “design-build-test-debug” cycle methodology, established here, will facilitate progression of the Sc2.0 project in the face of the increasing synthetic genome complexity. SynII characterization. ( A ) Cell cycle comparison between synII and BY4741 revealed by the percentage of cells with separated CEN2-GFP dots, metaphase spindles, and anaphase spindles. ( B ) Replication profiling of synII (red) and BY4741 (black) expressed as relative copy number by deep sequencing. ( C ) RNA sequencing analysis revealed that the significant up-regulation of translational machinery in synII is induced by the deletion of tRNA genes in synII.

Shen, Yue↗

Metagenomics-Guided Survey, Isolation, and Characterization of Uranium Resistant Microbiota from the Savannah River Site, USA

Despite the recent advancements in culturomics, isolation of the majority of environmental microbiota performing critical ecosystem services, such as bioremediation of contaminants, remains elusive. Towards this end, we conducted a metagenomics-guided comparative assessment of soil microbial diversity and functions present in uraniferous soils relative to those that grew in diffusion chambers (DC) or microbial traps (MT), followed by isolation of uranium (U) resistant microbiota. Shotgun metagenomic analysis performed on the soils used to establish the DC/MT chambers revealed Proteobacterial phyla and Burkholderia genus to be the most abundant among bacteria. The chamber-associated growth conditions further increased their abundances relative to the soils. Ascomycota was the most abundant fungal phylum in the chambers relative to the soils, with Penicillium as the most dominant genus. Metagenomics-based taxonomic findings completely mirrored the taxonomic composition of the retrieved isolates such that the U-resistant bacteria and fungi mainly belonged to Burkholderia and Penicillium species, thus confirming that the chambers facilitated proliferation and subsequent isolation of specific microbiota with environmentally relevant functions. Furthermore, shotgun metagenomic analysis also revealed that the gene classes for carbohydrate metabolism, virulence, and respiration predominated with functions related to stress response, membrane transport, and metabolism of aromatic compounds were also identified, albeit at lower levels. Of major note was the successful isolation of a potentially novel Penicillium species using the MT approach, as evidenced by whole genome sequence analysis and comparative genomic analysis, thus enhancing our overall understanding on the uranium cycling microbiota within the tested uraniferous soils.

Jaswal, Rajneesh↗

Optical follow-up study of 32 high-redshift galaxy cluster candidates from Planck with the William Herschel Telescope

Abstract The Planck satellite has detected cluster candidates via the Sunyaev Zel’dovich (SZ) effect, but the optical follow-up required to confirm these candidates is still incomplete, especially at high redshifts and for SZ detections at low significance. In this work, we present our analysis of optical observations obtained for 32 Planck cluster candidates using ACAM on the 4.2-m William Herschel Telescope. These cluster candidates were pre-selected using SDSS, WISE, and Pan-STARRS images to likely represent distant clusters at redshifts z ≳ 0.7. We obtain photometric redshift and richness estimates for all of the cluster candidates from a red-sequence analysis of r-, i-, and z-band imaging data. In addition, long-slit observations allow us to measure the redshifts of a subset of the clusters spectroscopically. The optical richness is often lower than expected from the inferred SZ mass when compared to scaling relations previously calibrated at low redshifts. This likely indicates the impact of Eddington bias and projection effects or noise-induced detections, especially at low-SZ significance. Thus, optical follow-up not only provides redshift measurements, but also an important independent verification method. We find that 18 (7) of the candidates at redshifts z > 0.5 (z > 0.8) are at least half as rich as expected from scaling relations, thereby clearly confirming these candidates as massive clusters. While the complex selection function of our sample due to our pre-selection hampers its use for cosmological studies, we do provide a validation of massive high-redshift clusters particularly suitable for further astrophysical investigations.

Zohren, Hannah↗

Detection of 224 candidate structured RNAs by comparative analysis of specific subsets of intergenic regions

The discovery of structured non-coding RNAs (ncRNAs) in bacteria can reveal new facets of biology and biochemistry. Comparative genomics analyses executed by powerful computer algorithms have successfully been used to uncover many novel bacterial ncRNA classes in recent years. However, this general search strategy favors the discovery of more common ncRNA classes, whereas progressively rarer classes are correspondingly more difficult to identify. In the current study, we confront this problem by devising several methods to select subsets of intergenic regions that can concentrate these rare RNA classes, thereby increasing the probability that comparative sequence analysis approaches will reveal their existence. By implementing these methods, we discovered 224 novel ncRNA classes, which include ROOL RNA, an RNA class averaging 581 nt and present in multiple phyla, several highly conserved and widespread ncRNA classes with properties that suggest sophisticated biochemical functions and a multitude of putative cis-regulatory RNA classes involved in a variety of biological processes. We expect that further research on these newly found RNA classes will reveal additional aspects of novel biology, and allow for greater insights into the biochemistry performed by ncRNAs.

59 BASIC BIOLOGICAL SCIENCES↗

An optimized FM-index library for nucleotide and amino acid search

Abstract Background Pattern matching is a key step in a variety of biological sequence analysis pipelines. The FM-index is a compressed data structure for pattern matching, with search run time that is independent of the length of the database text. Implementation of the FM-index is reasonably complicated, so that increased adoption will be aided by the availability of a fast and flexible FM-index library. Results We present AvxWindowedFMindex (AWFM-index), a lightweight, open-source, thread-parallel FM-index library written in C that is optimized for indexing nucleotide and amino acid sequences. AWFM-index introduces a new approach to storing FM-index data in a strided bit-vector format that enables extremely efficient computation of the FM-index occurrence function via AVX2 bitwise instructions, and combines this with optional on-disk storage of the index’s suffix array and a cache-efficient lookup table for partial k-mer searches. The AWFM-index performs exact match count and locate queries faster than SeqAn3’s FM-index implementation across a range of comparable memory footprints. When optimized for speed, AWFM-index is $$\sim $$ ∼ 2–4x faster than SeqAn3 for nucleotide search, and $$\sim $$ ∼ 2–6x faster for amino acid search; it is also $$\sim $$ ∼ 4x faster with similar memory footprint when storing the suffix array in on-disk SSD storage. Conclusions AWFM-index is easy to incorporate into bioinformatics software, offers run-time performance parameterization, and provides clients with FM-index functionality at both a high-level (count or locate all instances of a query string) and low-level (step-wise control of the FM-index backward-search process). The open-source library is available for download at https://github.com/TravisWheelerLab/AvxWindowFmIndex.

59 BASIC BIOLOGICAL SCIENCES↗

Organochlorine pesticides (OCPs) and polychlorinated biphenyls (PCBs) in Cyprinidae fish: Towards hints of their arrangements using advanced classification methods

Highlights: • SOM and DT classifiers searched for relation among 25 POPs in three Cyprinidae fish. • The concentrations of indicator congeners dominated over dioxin-like PCBs. • Based on SOM, both dioxin-like and indicator PCBs differentiated the fish species. • DT assorted the samples on the levels of PCB-126, -156, -170, -28, -105, and -114. To tackle the ever-present global concern regarding human exposure to persistent organic pollutants (POPs) via food products, this study strived to indicate associations between organochlorine pesticides (OCPs) and polychlorinated biphenyls (PCBs) in lake-fish tissue depending on the species and sampling season. Apart from the monitoring initiatives recommended in the Global Monitoring Plan for POPs, the study discussed 7 OCPs and 18 PCB congeners determined in three Cyprinidae species (rudd, carp, and Prussian carp) from Vransko Lake (Croatia), which are widely domesticated and reared as food fish across Europe and Asia. We exploit advanced classification algorithms, the Kohonen self-organizing maps (SOM) and Decision Trees (DT), to search for POP patterns typical for the investigated species. As indicated by SOM, some of the dioxin-like and non-dioxin-like PCBs (PCB-28, PCB-74, PCB-52, PCB-101, PCB-105, PCB-114, PCB-118, PCB-156 and PCB-157), α-HCH and β-HCH caused dissimilarities among fish species, but regardless of their weight and length. To support these suggestions, DT analysis sequenced the fish species and seasons based on the concentration of heavier congeners. The presented assumptions indicated that the supplemental application of SOM and DT offers advantageous features over the usually rough interpretation of POPs pattern and over the single use of the methods.

54 ENVIRONMENTAL SCIENCES↗

Transcriptomic Analyses Elucidate Adaptive Differences of Closely Related Strains of Pseudomonas aeruginosa in Fuel

Pseudomonas aeruginosa can utilize hydrocarbons, but different strains have various degrees of adaptation despite their highly conserved genome. P. aeruginosa ATCC 33988 is highly adapted to hydrocarbons, while P. aeruginosa strain PAO1, a human pathogen, is less adapted and degrades jet fuel at a lower rate than does ATCC 33988. We investigated fuel-specific transcriptomic differences between these strains in order to ascertain the underlying mechanisms utilized by the adapted strain to proliferate in fuel. During growth in fuel, the genes related to alkane degradation, heat shock response, membrane proteins, efflux pumps, and several novel genes were upregulated in ATCC 33988. Overexpression of alk genes in PAO1 provided some improvement in growth, but it was not as robust as that of ATCC 33988, suggesting the role of other genes in adaptation. Expression of the function unknown gene PA5359 from ATCC 33988 in PAO1 increased the growth in fuel. Bioinformatic analysis revealed that PA5359 is a predicted lipoprotein with a conserved Yx(FWY)xxD motif, which is shared among bacterial adhesins. Overexpression of the putative resistance-nodulation-division (RND) efflux pump PA3521 to PA3523 increased the growth of the ATCC 33988 strain, suggesting a possible role in fuel tolerance. Interestingly, the PAO1 strain cannot utilize n -C 8 and n -C 10 . The expression of green fluorescent protein (GFP) under the control of alkB promoters confirmed that alk gene promoter polymorphism affects the expression of alk genes. Promoter fusion assays further confirmed that the regulation of alk genes was different in the two strains. Protein sequence analysis showed low amino acid differences for many of the upregulated genes, further supporting transcriptional control as the main mechanism for enhanced adaptation. Lastly, these results support that specific signal transduction, gene regulation, and coordination of multiple biological responses are required to improve the survival, growth, and metabolism of fuel in adapted strains. This study provides new insight into the mechanistic differences between strains and helpful information that may be applied in the improvement of bacterial strains for resistance to biotic and abiotic factors encountered during bioremediation and industrial biotechnological processes.

54 ENVIRONMENTAL SCIENCES↗

Myotube-derived factor promotes oligodendrocyte precursor cell proliferation

Highlights: • Myotube-derived factors promote proliferation of oligodendrocyte precursor cells (OPCs). • Myotube-derived factors do not promote OPC migration and differentiation. • Candidate genes which promote OPC proliferation were identified by RNA sequence analysis. • Molecules released by myotubes may be involved in OPC proliferation. Muscle cells secrete numerous molecules that function as endocrine hormones and regulate the functions of distant organs. Myelination in the central nervous system (CNS) is regulated by peripheral hormones. However, the effects of muscle-derived molecules on myelination have not been sufficiently analyzed. In this study, we show that muscle-releasing factors promote proliferation of oligodendrocyte precursor cells (OPCs), which is an element of myelination process. Supernatants of mouse myotube cultures stimulated bromodeoxyuridine (BrdU) incorporation into mouse OPCs. Mouse myotube supernatants did not enhance mouse OPC transmigration and myelin basic protein (MBP) expression. RNA sequencing identified candidate genes with hormonal functions that were expressed in mouse myotubes. These data support the possibility that hormonal molecules secreted by myotubes contribute to OPC proliferation and myelination.

60 APPLIED LIFE SCIENCES↗

Advances in molecular serotyping and subtyping of Escherichia coli

Escherichia coli plays an important role as a member of the gut microbiota; however, pathogenic strains also exist, including various diarrheagenic E. coli pathotypes and extraintestinal pathogenic E. coli that cause illness outside of the GI-tract. E. coli have traditionally been serotyped using antisera against the ca. 186 O-antigens and 53 H-flagellar antigens. Phenotypic methods, including bacteriophage typing and O- and H- serotyping for differentiating and characterizing E. coli have been used for many years; however, these methods are generally time consuming and not always accurate. Advances in next generation sequencing technologies have made it possible to develop genetic-based subtyping and molecular serotyping methods for E. coli , which are more discriminatory compared to phenotypic typing methods. Furthermore, whole genome sequencing (WGS) of E. coli is replacing established subtyping methods such as pulsedfield gel electrophoresis, providing a major advancement in the ability to investigate food-borne disease outbreaks and for trace-back to sources. Furthermore, a variety of sequence analysis tools and bioinformatic pipelines are being developed to analyze the vast amount of data generated by WGS and to obtain specific information such as O- and H-group determination and the presence of virulence genes and other genetic markers.

59 BASIC BIOLOGICAL SCIENCES↗

Natural Allelic Variations in Highly Polyploidy Saccharum Complex

Sugarcane ( Saccharum spp.) is an important sugar and biofuel crop with high polyploid and complex genomes. The Saccharum complex, comprised of Saccharum genus and a few related genera, are important genetic resources for sugarcane breeding. A large amount of natural variation exists within the Saccharum complex. Though understanding their allelic variation has been challenging, it is critical to dissect allelic structure and to identify the alleles controlling important traits in sugarcane. To characterize natural variations in Saccharum complex, a target enrichment sequencing approach was used to assay 12 representative germplasm accessions. In total, 55,946 highly efficient probes were designed based on the sorghum genome and sugarcane unigene set targeting a total of 6 Mb of the sugarcane genome. A pipeline specifically tailored for polyploid sequence variants and genotype calling was established. BWAmem and sorghum genome approved to be an acceptable aligner and reference for sugarcane target enrichment sequence analysis, respectively. Genetic variations including 1,166,066 non-redundant SNPs, 150,421 InDels, 919 gene copy number variations, and 1,257 gene presence/absence variations were detected. SNPs from three different callers (Samtools, Freebayes, and GATK) were compared and the validation rates were nearly 90%. Based on the SNP loci of each accession and their ploidy levels, 999,258 single dosage SNPs were identified and most loci were estimated as largely homozygotes. An average of 34,397 haplotype blocks for each accession was inferred. The highest divergence time among the Saccharum spp. was estimated as 1.2 million years ago (MYA). Saccharum spp. diverged from Erianthus and Sorghum approximately 5 and 6 MYA, respectively. Furthermore, the target enrichment sequencing approach provided an effective way to discover and catalog natural allelic variation in highly polyploid or heterozygous genomes.

59 BASIC BIOLOGICAL SCIENCES↗

ORIGEN-based Nuclear Fuel Inventory Module for Fuel Cycle Assessment: Final Project Report

The goal of this project, “ORIGEN-based Nuclear Fuel Depletion Module for Fuel Cycle Assessment" is to create a physics-based reactor depletion and decay module for the Cyclus nuclear fuel cycle simulator in order to assess nuclear fuel inventories over a broad space of reactor operating conditions. The overall goal of this approach is to facilitate evaluations of nuclear fuel inventories for a broad space of scenarios, including extended used nuclear fuel storage and cascading impacts on fuel cycle options such as actinide recovery in used nuclear fuel, particularly for multiple recycle scenarios. The advantages of a physics-based approach (compared to a recipe-based approach which has been typically employed for fuel cycle simulators) is in its inherent flexibility; such an approach can more readily accommodate the broad space of potential isotopic vectors that may be encountered under advanced fuel cycle options. In order to develop this flexible reactor analysis capability, we are leveraging the Origen nuclear fuel depletion and decay module from SCALE to produce a standalone “depletion engine” which will serve as the kernel of a Cyclus-based reactor analysis module. The ORIGEN depletion module is a rigorously benchmarked and extensively validated tool for nuclear fuel analysis and thus its incorporation into the Cyclus framework can bring these capabilities to bear on the problem of evaluating long-term impacts of fuel cycle option choices on relevant metrics of interest, including materials inventories and availability (for multiple recycle scenarios), long-term waste management and repository impacts, etc. Developing this Origen-based analysis capability for Cyclus requires the refinement of the Origen analysis sequence to the point where it can reasonably be compiled as a standalone sequence outside of SCALE; i.e., wherein all of the computational aspects of Origen (including reactor cross-section library processing and interpolation, input and output processing, and depletion/decay solvers) can be self-contained into a single executable sequence. Further, to embed this capability into other software environments (such as the Cyclus fuel cycle simulator) requires that Origen’s capabilities be encapsulated into a portable, self-contained library which other codes can then call directly through function calls, thereby directly accessing the solver and data processing capabilities of Origen. Additional components relevant to this work include modernization of the reactor data libraries used by Origen for conducting nuclear fuel depletion calculations. This work has included the development of new fuel assembly lattices not previously available (such as for CANDU heavy-water reactor assemblies) as well as validation of updated lattices for light-water reactors updated to employ modern nuclear data evaluations. The CyBORG reactor analysis module as-developed under this workscope is fully capable of dynamic calculation of depleted fuel compositions from all commercial U.S. reactor assembly types as well as a number of international fuel types, including MOX, VVER, MAGNOX, and PHWR CANDU fuel assemblies. In addition, the Origen-based depletion engine allows for CyBORG to evaluate novel fuel assembly and reactor design types via creation of Origen reactor data libraries via SCALE. The establishment of this new modeling capability affords fuel cycle modelers a substantially improved ability to model dynamically-changing fuel cycle and reactor conditions, including recycled fuel compositions from fuel cycle scenarios involving material recycle into thermal-spectrum systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Characterization of resistance to a potent D -peptide HIV entry inhibitor

Background: PIE12-trimer is a highly potent D-peptide HIV-1 entry inhibitor that broadly targets group M isolates. It specifically binds the three identical conserved hydrophobic pockets at the base of the gp41 N-trimer with sub-femtomolar affinity. This extremely high affinity for the transiently exposed gp41 trimer provides a reserve of binding energy (resistance capacitor) to prevent the viral resistance pathway of stepwise accumulation of modest affinity-disrupting mutations. Such modest mutations would not affect PIE12-trimer potency and therefore not confer a selective advantage. Viral passaging in the presence of escalating PIE12-trimer concentrations ultimately selected for PIE12-trimer resistant populations, but required an extremely extended timeframe (> 1 year) in comparison to other entry inhibitors. Eventually, HIV developed resistance to PIE12-trimer by mutating Q577 in the gp41 pocket. Results: Using deep sequence analysis, we identified three mutations at Q577 (R, N and K) in our two PIE12-trimer resistant pools. Each point mutant is capable of conferring the majority of PIE12-trimer resistance seen in the polyclonal pools. Surface plasmon resonance studies demonstrated substantial affinity loss between PIE12-trimer and the Q577R-mutated gp41 pocket. A high-resolution X-ray crystal structure of PIE12 bound to the Q577R pocket revealed the loss of two hydrogen bonds, the repositioning of neighboring residues, and a small decrease in buried surface area. The Q577 mutations in an NL4-3 backbone decreased viral growth rates. Fitness was ultimately rescued in resistant viral pools by a suite of compensatory mutations in gp120 and gp41, of which we identified seven candidates from our sequencing data. Conclusions: These data show that PIE12-trimer exhibits a high barrier to resistance, as extended passaging was required to develop resistant virus with normal growth rates. The primary resistance mutation, Q577R/N/K, found in the conserved gp41 pocket, substantially decreases inhibitor affinity but also damages viral fitness, and candidate compensatory mutations in gp160 have been identified.

59 BASIC BIOLOGICAL SCIENCES↗

Integration of Russian Tier-1 Grid Center with High Performance Computers at NRC-KI for LHC experiments and beyond HENP

The LHC experiments are preparing for the precision measurements and further discoveries that will be made possible by higher LHC energies from April 2015 (LHC Run2). The need for simulation, data processing and analysis would overwhelm the expected capacity of grid infrastructure computing facilities deployed by the Worldwide LHC Computing Grid (WLCG). To meet this challenge the integration of the opportunistic resources into LHC computing model is highly important. The Tier-1 facility at Kurchatov Institute (NRC-KI) in Moscow is a part of WLCG and it will process, simulate and store up to 10% of total data obtained from ALICE, ATLAS and LHCb experiments. In addition Kurchatov Institute has supercomputers with peak performance 0.12 PFLOPS. The delegation of even a fraction of supercomputing resources to the LHC Computing will notably increase total capacity. In 2014 the development a portal combining a Tier-1 and a supercomputer in Kurchatov Institute was started to provide common interfaces and storage. The portal will be used not only for HENP experiments, but also by other data- and compute-intensive sciences like biology with genome sequencing analysis; astrophysics with cosmic rays analysis, antimatter and dark matter search, etc.

97 MATHEMATICS AND COMPUTING↗