Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consensus sequence”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Matrix Metalloproteinases as Candidate Antigenic Determinants for Anti‐Tumor Autoantibodies in Human Ovarian Cancer: A Post Hoc Analysis

Circulating antibodies in patients with cancer can facilitate the identification of accessible epitopes on autoantigens expressed by tumors. To identify previously unrecognized protein targets in ovarian cancer, we computationally assessed a heptapeptide consensus motif (VPELGHE, flanked by two cysteine residues yielding a cyclic nonapeptide under oxidizing conditions) previously discovered via phage display-based epitope mapping of autoantibodies in patients. Eight proteins associated with ovarian cancer encompass amino acid sequences similar to the consensus motif and were, therefore, considered as candidate native autoantigens. Among these candidate targets, however, matrix metalloproteinase 14 (MMP14) demonstrates gene expression that is both high and negatively correlated with survival in ovarian cancer patient cohorts. MMP14 protein levels are also stable in tumor versus non-tumor tissues. Moreover, the corresponding heptapeptide mimic in MMP14 occurs within an α-helical secondary structural element observed in its catalytic domain. These findings demonstrate that a subset of patient-derived autoantibodies may interact with a previously unknown antigenic epitope found in MMP14 and other MMPs, thereby providing opportunities for the development of new targeted agents.

Biochemistry & Molecular Biology↗
DOE OSTI2021

CONSTAX2: improved taxonomic classification of environmental DNA markers

Abstract Summary CONSTAX—the CONSensus TAXonomy classifier—was developed for accurate and reproducible taxonomic annotation of fungal rDNA amplicon sequences and is based upon a consensus approach of RDP, SINTAX and UTAX algorithms. CONSTAX2 extends these features to classify prokaryotes as well as eukaryotes and incorporates BLAST-based classifiers to reduce classification errors. Additionally, CONSTAX2 implements a conda-installable command-line tool with improved classification metrics, faster training, multithreading support, capacity to incorporate external taxonomic databases and new isolate matching and high-level taxonomy tools, replete with documentation and example tutorials. Availability and implementation CONSTAX2 is available at https://github.com/liberjul/CONSTAXv2, and is packaged for Linux and MacOS from Bioconda with use under the MIT License. A tutorial and documentation are available at https://constax.readthedocs.io/en/latest/. Data and scripts associated with the manuscript are available at https://github.com/liberjul/CONSTAXv2_ms_code. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Recombinant And Mix-Infection Finder for SARS-CoV-2 sample

The scientific and public health communities responded to the COVID-19 pandemic with sample acquisition and genome sequencing on a scale that eclipsed all prior sequencing efforts. While this can only be characterized as a resounding success story that has cemented the use of genomics for epidemiological investigations for any future infectious disease outbreak, several retrospective studies are cataloging an array of lessons learned and issues that have yet to be addressed in order to realize the full potential of genomics as a routine biosurveillance tool. We have been both developing methods to accurately assess SARS-CoV-2 genomes from complex samples, and analyzing the large volumes of international data, both at the consensus level and the raw sequencing data. During the course of our investigations and similar to other groups, we have examined COVID-19 samples with signatures from multiple lineages of SARS-CoV-2 and will describe some of our findings during the development of a novel workflow that incorporates detection and reporting of potential co-infection within samples and also highlights any evidence of within-host recombination.

Lo, Chien-Chi↗
DOE OSTI2023

U.S. Efforts in Support of Examinations at Fukushima Daiichi - November 2022 Meeting Notes and Information Request Status

Information obtained from Fukushima Daiichi Nuclear Power Station (Daiichi) is required to inform future Decontamination and Decommissioning (D&D) activities, improving the ability of the Tokyo Electric Power Company Holdings, Incorporated (TEPCO Holdings) to characterize potential hazards and to ensure the safety of workers involved with cleanup activities. This information also has important implications for the safety and operation of U.S. commercial nuclear power plants. This document summarizes results from the Fiscal Year 2023 (FY2023) U.S. effort to review Daiichi information and extract insights to enhance the safety of existing and future nuclear power plant designs. This U.S. effort, which was initiated in 2014 by the Department of Energy Office of Nuclear Energy, is completed by a group of experts in reactor safety and plant operations that identify examination needs and evaluate recent Daiichi examination data to address these needs. Fukushima-related information and associated discussions during these meetings benefit operating, new, and advanced reactors. Significant safety insights have been and are continuing to be obtained in several areas: system and component performance, radionuclide surveys and sampling, debris end-state location, combustible gas effects, and plant operations and maintenance. In addition to reducing uncertainties related to severe accident modeling progression, these insights have and continue to be used to update guidance for severe accident prevention, mitigation, and emergency planning. Furthermore, Daiichi-related activities, such as code modeling improvements and analysis, testing, and new technology deployment efforts, have the potential to offer additional benefits to the operating fleet and new LWR and non-LWR designs. U.S. evaluations of obtained examination information and input regarding future Daiichi examinations are of interest to several organizations within Japan. Since its inception, the U.S. has provided consensus input for high priority time-sequenced examination tasks and supporting research activities. In their Mid-to-Long-term Examination Plan for 1F investigations, TEPCO included all remaining U.S. consensus information requests and additional information requests they identified. TEPCO periodically provides reports on the status of these requests (reflecting D&D priorities, new insights from investigations, and new technologies that become available). Hence, U.S. experts agreed that it was appropriate for TEPCO to track and prioritize these information requests as D&D progresses. U.S. experts will continue to review and comment on the information obtained from examinations and, as needed, provide additional details and relevant background material to support future examinations. As documented in this report, several other items, such as additional details on information requests pertaining to ex-vessel examinations, relevant references from prior research, additional documents to provide insights regarding recent investigation findings, and reviews of recently released documents, were agreed to during the FY2023 meeting.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

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↗
DOE OSTI2024

Recommendations for Uniform Variant Calling of SARS-CoV-2 Genome Sequence across Bioinformatic Workflows

Genomic sequencing of clinical samples to identify emerging variants of SARS-CoV-2 has been a key public health tool for curbing the spread of the virus. As a result, an unprecedented number of SARS-CoV-2 genomes were sequenced during the COVID-19 pandemic, which allowed for rapid identification of genetic variants, enabling the timely design and testing of therapies and deployment of new vaccine formulations to combat the new variants. However, despite the technological advances of deep sequencing, the analysis of the raw sequence data generated globally is neither standardized nor consistent, leading to vastly disparate sequences that may impact identification of variants. Here, we show that for both Illumina and Oxford Nanopore sequencing platforms, downstream bioinformatic protocols used by industry, government, and academic groups resulted in different virus sequences from same sample. These bioinformatic workflows produced consensus genomes with differences in single nucleotide polymorphisms, inclusion and exclusion of insertions, and/or deletions, despite using the same raw sequence as input datasets. Here, we compared and characterized such discrepancies and propose a specific suite of parameters and protocols that should be adopted across the field. Consistent results from bioinformatic workflows are fundamental to SARS-CoV-2 and future pathogen surveillance efforts, including pandemic preparation, to allow for a data-driven and timely public health response.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Parallel String Graph Construction and Transitive Reduction for De Novo Genome Assembly

One of the most computationally intensive tasks in computational biology is de novo genome assembly, the decoding of the sequence of an unknown genome from redundant and erroneous short sequences. A common assembly paradigm identifies overlapping sequences, simplifies their layout, and creates consensus. Despite many algorithms developed in the literature, the efficient assembly of large genomes is still an open problem. In this work, we introduce new distributed-memory parallel algorithms for overlap detection and layout simplification steps of de novo genome assembly, and implement them in the diBELLA 2D pipeline. Our distributed memory algorithms for both overlap detection and layout simplification are based on linear-algebra operations over semirings using 2D distributed sparse matrices. Our layout step consists of performing a transitive reduction from the overlap graph to a string graph. We provide a detailed communication analysis of the main stages of our new algorithms. diBELLA 2D achieves near linear scaling with over 80% parallel efficiency for the human genome, reducing the runtime for overlap detection by 1.2-1.3× for the human genome and 1.5-1.9× for C.elegans compared to the state-of-the-art. Our transitive reduction algorithm outperforms an existing distributed-memory implementation by 10.5-13.3× for the human genome and 18-29× for the C. elegans. Our work paves the way for efficient de novo assembly of large genomes using long reads in distributed memory.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A minimal motif for sequence recognition by mitochondrial transcription factor A (TFAM)

Abstract Mitochondrial transcription factor A (TFAM) plays a critical role in mitochondrial transcription initiation and mitochondrial DNA (mtDNA) packaging. Both functions require DNA binding, but in one case TFAM must recognize a specific promoter sequence, while packaging requires coating of mtDNA by association with non sequence-specific regions. The mechanisms by which TFAM achieves both sequence-specific and non sequence-specific recognition have not yet been determined. Existing crystal structures of TFAM bound to DNA allowed us to identify two guanine-specific interactions that are established between TFAM and the bound DNA. These interactions are observed when TFAM is bound to both specific promoter sequences and non-sequence specific DNA. These interactions are established with two guanine bases separated by 10 random nucleotides (GN10G). Our biochemical results demonstrate that the GN10G consensus is essential for transcriptional initiation and contributes to facilitating TFAM binding to DNA substrates. Furthermore, we report a crystal structure of TFAM in complex with a non sequence-specific sequence containing a GN10G consensus. The structure reveals a unique arrangement in which TFAM bridges two DNA substrates while maintaining the GN10G interactions. We propose that the GN10G consensus is key to facilitate the interaction of TFAM with DNA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

PCP consensus protein/peptide alphavirus antigens stimulate broad spectrum neutralizing antibodies

Vaccines based on proteins and peptides may be safer and if calculated based on many sequences, more broad-spectrum than those designed based on single strains. Physicochemical Property Consensus (PCP con ) alphavirus (AV) antigens from the B-domain of the E2 envelope protein were designed, synthesized recombinantly and shown to be immunogenic (i.e. sera after inoculation detected the antigen in dotspots and ELISA). Antibodies in sera after inoculation with B-region antigens based on individual AV species (eastern or Venezuelan equine encephalitis (EEEV con , VEEV con ), or chikungunya (CHIKV con ) bound only their cognate protein, while those designed against multiple species (Mosaik con and EVC con ) recognized all three serotype specific antigens. The VEEV con and EEEV con sera only showed antiviral activity against their related strains (in plaque reduction neutralization assays (PRNT 50/80 ). Peptides designed to surface exposed areas of the E2-A-domain of CHIKV con were added to CHIKV con inocula to provide anti-CHIKV antibodies. EVC con , based on three different alphavirus species, combined with E2-A-domain peptides from AllAV con , a PCPcon of 24 diverse AV, generated broad spectrum, antiviral antibodies against VEEV, EEEV and CHIKV, AV with less than 35% amino acid identity to each other (>65% diversity). This is a promising start to a molecularly defined vaccine against all AV. A further study with these antigens can illuminate what areas are most important for a robust immune response, resistant to mutations in rapidly evolving viruses. The validated computational methods can also be used to design broad spectrum antigens against many other pathogen families.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

NASA GeneLab RNA-seq consensus pipeline: Standardized processing of short-read RNA-seq data

With the development of transcriptomic technologies, we are able to quantify precise changes in gene expression profiles from astronauts and other organisms exposed to spaceflight. Members of NASA GeneLab and GeneLab-associated analysis working groups (AWGs) have developed a consensus pipeline for analyzing short-read RNA-sequencing data from spaceflight-associated experiments. The pipeline includes quality control, read trimming, mapping, and gene quantification steps, culminating in the detection of differentially expressed genes. This data analysis pipeline and the results of its execution using data submitted to GeneLab are now all publicly available through the GeneLab database. We present here the full details and rationale for the construction of this pipeline in order to promote transparency, reproducibility, and reusability of pipeline data; to provide a template for data processing of future spaceflight-relevant datasets; and to encourage cross-analysis of data from other databases with the data available in GeneLab.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

An expanded role for the transcription factor WRINKLED1 in the biosynthesis of triacylglycerols during seed development

The transcription factor WRINKLED1 ( WRI1 ) is known as a master regulator of fatty acid synthesis in developing oilseeds of Arabidopsis thaliana and other species. WRI1 is known to directly stimulate the expression of many fatty acid biosynthetic enzymes and a few targets in the lower part of the glycolytic pathway. However, it remains unclear to what extent and how the conversion of sugars into fatty acid biosynthetic precursors is controlled by WRI 1. To shortlist possible gene targets for future in-planta experimental validation, here we present a strategy that combines phylogenetic foot printing of cis-regulatory elements with additional layers of evidence. Upstream regions of protein-encoding genes in A. thaliana were searched for the previously described DNA-binding consensus for WRI1, the ASML1/WRI1 (AW)-box. For about 900 genes, AW-box sites were found to be conserved across orthologous upstream regions in 11 related species of the crucifer family. For 145 select potential target genes identified this way, affinity of upstream AW-box sequences to WRI1 was assayed by Microscale Thermophoresis. This allowed definition of a refined WRI1 DNA-binding consensus. We find that known WRI1 gene targets are predictable with good confidence when upstream AW-sites are phylogenetically conserved, specifically binding WRI1 in the in vitro assay, positioned in proximity to the transcriptional start site, and if the gene is co-expressed with WRI1 during seed development. When targets predicted in this way are mapped to central metabolism, a conserved regulatory blueprint emerges that infers concerted control of contiguous pathway sections in glycolysis and fatty acid biosynthesis by WRI1. Several of the newly predicted targets are in the upper glycolysis pathway and the pentose phosphate pathway. Of these, plastidic isoforms of fructokinase ( FRK 3) and of phosphoglucose isomerase ( PGI 1) are particularly corroborated by previously reported seed phenotypes of respective null mutations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Intestinal single-cell atlas reveals novel lymphocytes in pigs with similarities to human cells

Lymphocytes can heavily influence intestinal health, but resolving intestinal lymphocyte function is challenging as the intestine contains a vastly heterogeneous mixture of cells. Pigs are an advantageous biomedical model, but deeper understanding of intestinal lymphocytes is warranted to improve model utility. Twenty-six cell types were identified in the porcine ileum by single-cell RNA sequencing and further compared with cells in human and murine ileum. Though general consensus of cell subsets across species was revealed, some porcine-specific lymphocyte subsets were identified. Differential tissue dissection and in situ analyses conferred spatial context, revealing similar locations of lymphocyte subsets in Peyer’s patches and epithelium in pig-to-human comparisons. Like humans, activated and effector lymphocytes were abundant in the ileum but not periphery of pigs, suggesting tissue-specific and/or activation-associated gene expression. Gene signatures for peripheral and ileal innate lymphoid cells newly discovered in pigs were defined and highlighted similarities to human innate lymphoid cells. Overall, we reveal novel lymphocyte subsets in pigs and highlight utility of pigs for intestinal research applications.

Life Sciences & Biomedicine - Other Topics↗
DOE OSTI2024

The atomic gas sequence and mass–metallicity relation from dwarfs to massive galaxies

Galaxy scaling relations provide insights into the processes that drive galaxy evolution. The extension of these scaling relations into the dwarf galaxy regime is of particular interest. This is because dwarf galaxies represent a crucial stage in galaxy evolution, and understanding them could also shed light on their role in reionizing the early Universe. There is currently no consensus on the processes that dominate the evolution of dwarfs. In this work, we constrain the atomic gas sequence (stellar mass versus atomic gas fraction) and mass–metallicity relation (stellar mass versus gas-phase metallicity) from dwarf (⁠10 6.5 M ⊙ ⁠) to massive (⁠10 11.5 M ⊙ ) galaxies in the local Universe. The combined optical and 21-cm spectroscopic observations of the Dark Energy Spectroscopic Instrument and Arecibo Legacy Fast ALFA surveys allow us to constrain both scaling relations simultaneously. We find a slope change of the atomic gas sequence at a stellar mass of ~10 9 M ⊙ ⁠. We also find that the shape and scatter of the atomic gas sequence and mass–metallicity relation are strongly linked for both dwarfs and more massive galaxies. Consequently, the low-mass slope change of the atomic gas sequence is imprinted onto the mass–metallicity relation of dwarf galaxies. The mass scale of the measured slope change is consistent with a predicted escape velocity threshold below which low-mass galaxies experience significant supernova-driven gas loss, as well as with a reduction in cold gas accretion onto more massive galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Tousled-like kinase 2 targets ASF1 histone chaperones through client mimicry

Tousled-like kinases (TLKs) are nuclear serine-threonine kinases essential for genome maintenance and proper cell division in animals and plants. A major function of TLKs is to phosphorylate the histone chaperone proteins ASF1a and ASF1b to facilitate DNA replication-coupled nucleosome assembly, but how TLKs selectively target these critical substrates is unknown. Here, we show that TLK2 selectivity towards ASF1 substrates is achieved in two ways. First, the TLK2 catalytic domain recognizes consensus phosphorylation site motifs in the ASF1 C-terminal tail. Second, a short sequence at the TLK2 N-terminus docks onto the ASF1a globular N-terminal domain in a manner that mimics its histone H3 client. Disrupting either catalytic or non-catalytic interactions through mutagenesis hampers ASF1 phosphorylation by TLK2 and cell growth. Our results suggest that the stringent selectivity of TLKs for ASF1 is enforced by an unusual interaction mode involving mutual recognition of a short sequence motifs by both kinase and substrate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Detection of Candida DNA in peritoneal fluids by PCR assay optimizing the diagnosis and treatment for intra-abdominal candidiasis in high-risk ICU patients: A prospective cohort study

Background Intra-abdominal candidiasis (IAC) is the predominant type of invasive candidiasis with high mortality in critically ill patients. This study aimed to investigate whether the polymerase chain reaction (PCR) assay for detecting Candida DNA in peritoneal fluids (PF) is useful in diagnosing and management of IAC in high-risk patients in intensive care unit (ICU). Methods A prospective single-center cohort study of surgical patients at high risk for IAC was conducted in the ICU. PF was collected from the abdominal drainage tubes (within 24 h) or by percutaneous puncture. Direct PF smear microscopy, PF culture, blood culture, and serum (1–3)-β-D-glucan were performed in all patients. For Candida PCR assay, the ITS1/ITS4 primers that targeted the ITS1-5.8 s-ITS2 regions were used for PCR, and sequencing analysis was used to identify the pathogen at the species level. IAC was defined according to the 2013 European consensus criteria. Results Among 83 patients at high risk for IAC, the IAC criteria were present in 17 (20.5%). The sensitivity and specificity of the Candida PCR assay were 64.7 and 89.4%, respectively, and the area under the receiver operating characteristic curve was 0.77 (95% CI: 0.63–0.91). In this cohort, the positive predictive value and negative predictive value were 90.8% (95% CI: 80.3–96.2%) and 61.1% (95% CI: 36.1–81.7%), respectively. Diagnostic consistency was moderate (kappa 0.529, p < 0.001) according to the 2013 European consensus criteria. Conclusion Detection of Candida DNA in PF using PCR can be considered an adjunct to existing routine diagnostic tools which may optimize the diagnosis and antifungal treatment of IAC in high-risk patients in the ICU.

Xie, Min↗
DOE OSTI2024

Watching a signaling protein function: What has been learned over four decades of time-resolved studies of photoactive yellow protein

Photoactive yellow protein (PYP) is a signaling protein whose internal p-coumaric acid chromophore undergoes reversible, light-induced trans-to-cis isomerization, which triggers a sequence of structural changes that ultimately lead to a signaling state. Since its discovery nearly 40 years ago, PYP has attracted much interest and has become one of the most extensively studied proteins found in nature. The method of time-resolved crystallography, pioneered by Keith Moffat, has successfully characterized intermediates in the PYP photocycle at near atomic resolution over 12 decades of time down to the sub-picosecond time scale, allowing one to stitch together a movie and literally watch a protein as it functions. But how close to reality is this movie? To address this question, results from numerous complementary time-resolved techniques including x-ray crystallography, x-ray scattering, and spectroscopy are discussed. Emerging from spectroscopic studies is a general consensus that three time constants are required to model the excited state relaxation, with a highly strained ground-state cis intermediate formed in less than 2.4 ps. Persistent strain drives the sequence of structural transitions that ultimately produce the signaling state. Crystal packing forces produce a restoring force that slows somewhat the rates of interconversion between the intermediates. Moreover, the solvent composition surrounding PYP can influence the number and structures of intermediates as well as the rates at which they interconvert. When chloride is present, the PYP photocycle in a crystal closely tracks that in solution, which suggests the epic movie of the PYP photocycle is indeed based in reality.

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

Structural variant analysis of a cancer reference cell line sample using multiple sequencing technologies

The cancer genome is commonly altered with thousands of structural rearrangements including insertions, deletions, translocation, inversions, duplications, and copy number variations. Thus, structural variant (SV) characterization plays a paramount role in cancer target identification, oncology diagnostics, and personalized medicine. As part of the SEQC2 Consortium effort, the present study established and evaluated a consensus SV call set using a breast cancer reference cell line and matched normal control derived from the same donor, which were used in our companion benchmarking studies as reference samples. We systematically investigated somatic SVs in the reference cancer cell line by comparing to a matched normal cell line using multiple NGS platforms including Illumina short-read, 10X Genomics linked reads, PacBio long reads, Oxford Nanopore long reads, and high-throughput chromosome conformation capture (Hi-C). We established a consensus SV call set of a total of 1788 SVs including 717 deletions, 230 duplications, 551 insertions, 133 inversions, 146 translocations, and 11 breakends for the reference cancer cell line. To independently evaluate and cross-validate the accuracy of our consensus SV call set, we used orthogonal methods including PCR-based validation, Affymetrix arrays, Bionano optical mapping, and identification of fusion genes detected from RNA-seq. We evaluated the strengths and weaknesses of each NGS technology for SV determination, and our findings provide an actionable guide to improve cancer genome SV detection sensitivity and accuracy. A high-confidence consensus SV call set was established for the reference cancer cell line. A large subset of the variants identified was validated by multiple orthogonal methods.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Colour gradients of low-redshift galaxies in the DESI Legacy Imaging Survey

ABSTRACT Radial colour gradients within galaxies arise from gradients of stellar age, metallicity, and dust reddening. Large samples of colour gradients from wide-area imaging surveys can complement smaller integral-field spectroscopy data sets and can be used to constrain galaxy formation models. Here, we measure colour gradients for low-redshift galaxies (z < 0.1) using photometry from the DESI Legacy Imaging Survey DR9. Our sample comprises ∼93 000 galaxies with spectroscopic redshifts and ∼574 000 galaxies with photometric redshifts. We focus on gradients across a radial range 0.5Reff to Reff, which corresponds to the inner disc of typical late-type systems at low redshift. This region has been the focus of previous statistical studies of colour gradients and has recently been explored by spectroscopic surveys such as MaNGA. We find that the colour gradients of most galaxies in our sample are negative (redder towards the centre), consistent with the literature. We investigate empirical relationships between colour gradient, average g − r and r − z colour, Mr, M⋆, and sSFR. Trends of gradient strength with Mr (M⋆) show an inflection around Mr ∼ −21 ($\log _{10} \, M_\star /\mathrm{M_\odot }\sim 10.5$). Below this mass, colour gradients become steeper with increasing M⋆, whereas colour gradients in more massive galaxies become shallower. We find that positive gradients (bluer stars at smaller radii) are typical for galaxies of $M_{\star }\sim 10^{8}\, \mathrm{M_\odot }$. We compare our results to age and metallicity gradients in two data sets derived from fits of different stellar population libraries to MaNGA spectra, but find no clear consensus explanation for the trends we observe. Both MaNGA data sets seem to imply a significant contribution from dust reddening, in particular, to explain the flatness of colour gradients along the red sequence.

Astronomy & Astrophysics↗