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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1993

Sequence and structural implications of a bovine corneal keratan sulfate proteoglycan core protein. Protein 37B represents bovine lumican and proteins 37A and 25 are unique

Amino acid sequence from tryptic peptides of three different bovine corneal keratan sulfate proteoglycan (KSPG) core proteins (designated 37A, 37B, and 25) showed similarities to the sequence of a chicken KSPG core protein lumican. Bovine lumican cDNA was isolated from a bovine corneal expression library by screening with chicken lumican cDNA. The bovine cDNA codes for a 342-amino acid protein, M(r) 38,712, containing amino acid sequences identified in the 37B KSPG core protein. The bovine lumican is 68% identical to chicken lumican, with an 83% identity excluding the N-terminal 40 amino acids. Location of 6 cysteine and 4 consensus N-glycosylation sites in the bovine sequence were identical to those in chicken lumican. Bovine lumican had about 50% identity to bovine fibromodulin and 20% identity to bovine decorin and biglycan. About two-thirds of the lumican protein consists of a series of 10 amino acid leucine-rich repeats that occur in regions of calculated high beta-hydrophobic moment, suggesting that the leucine-rich repeats contribute to beta-sheet formation in these proteins. Sequences obtained from 37A and 25 core proteins were absent in bovine lumican, thus predicting a unique primary structure and separate mRNA for each of the three bovine KSPG core proteins.

NASA Discipline Cell Biology↗
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↗
NASA NTRS2004

cWINNOWER algorithm for finding fuzzy dna motifs

The cWINNOWER algorithm detects fuzzy motifs in DNA sequences rich in protein-binding signals. A signal is defined as any short nucleotide pattern having up to d mutations differing from a motif of length l. The algorithm finds such motifs if a clique consisting of a sufficiently large number of mutated copies of the motif (i.e., the signals) is present in the DNA sequence. The cWINNOWER algorithm substantially improves the sensitivity of the winnower method of Pevzner and Sze by imposing a consensus constraint, enabling it to detect much weaker signals. We studied the minimum detectable clique size qc as a function of sequence length N for random sequences. We found that qc increases linearly with N for a fast version of the algorithm based on counting three-member sub-cliques. Imposing consensus constraints reduces qc by a factor of three in this case, which makes the algorithm dramatically more sensitive. Our most sensitive algorithm, which counts four-member sub-cliques, needs a minimum of only 13 signals to detect motifs in a sequence of length N = 12,000 for (l, d) = (15, 4). Copyright Imperial College Press.

Evaluation Studies↗
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↗
NASA NTRS2003

cWINNOWER Algorithm for Finding Fuzzy DNA Motifs

The cWINNOWER algorithm detects fuzzy motifs in DNA sequences rich in protein-binding signals. A signal is defined as any short nucleotide pattern having up to d mutations differing from a motif of length l. The algorithm finds such motifs if multiple mutated copies of the motif (i.e., the signals) are present in the DNA sequence in sufficient abundance. The cWINNOWER algorithm substantially improves the sensitivity of the winnower method of Pevzner and Sze by imposing a consensus constraint, enabling it to detect much weaker signals. We studied the minimum number of detectable motifs qc as a function of sequence length N for random sequences. We found that qc increases linearly with N for a fast version of the algorithm based on counting three-member sub-cliques. Imposing consensus constraints reduces qc, by a factor of three in this case, which makes the algorithm dramatically more sensitive. Our most sensitive algorithm, which counts four-member sub-cliques, needs a minimum of only 13 signals to detect motifs in a sequence of length N = 12000 for (l,d) = (15,4).

Liang, Shoudan↗
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↗
NASA NTRS2021

PURE mRNA display and cDNA display provide rapid detection of core epitope motif via high‐throughput sequencing

The reconstructed in vitro translation system known as the PURE system has been used in a variety of cell‐free experiments such as the expression of native and de novo proteins as well as various display methods to select for functional polypeptides. We developed a refined PURE‐based display method for the preparation of stable messenger RNA (mRNA) and complementary DNA (cDNA)‐peptide conjugates and validated its utility for in vitro selection. Our conjugate formation efficiency exceeded 40%, followed by gel purification to allow minimum carry‐over of components from the translation system to the downstream assay enabling clean and efficient random peptide sequence screening. We chose the commercially available anti‐FLAG M2 antibody as a target molecule for validation. Starting from approximately 1.7 × 10(exp 12) random sequences, a round‐by‐round high‐throughput sequencing showed clear enrichment of the FLAG epitope DYKDDD as well as revealing consensus FLAG epitope motif DYK(D/L/N)(L/Y/D/N/F)D. Enrichment of core FLAG motifs lacking one of the four key residues (DYKxxD) indicates that Tyr(Y) and Lys (K) appear as the two key residues essential for binding. Furthermore, the comparison between mRNA display and cDNA display method resulted in overall similar performance with slightly higher enrichment for mRNA display. We also show that gel purification steps in the refined PURE‐based display method improve conjugate formation efficiency and enhance the enrichment rate of FLAG epitope motifs in later rounds of selection especially for mRNA display. Overall, the generalized procedure and consistent performance of two different display methods achieved by the commercially available PURE system will be useful for future studies to explore the sequence and functional space of diverse polypeptides.

cDNA display, FLAG epitope, mRNA display, peptide ↗
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↗
NASA NTRS2011

Local Estimators for Spacecraft Formation Flying

A formation estimation architecture for formation flying builds upon the local information exchange among multiple local estimators. Spacecraft formation flying involves the coordination of states among multiple spacecraft through relative sensing, inter-spacecraft communication, and control. Most existing formation flying estimation algorithms can only be supported via highly centralized, all-to-all, static relative sensing. New algorithms are needed that are scalable, modular, and robust to variations in the topology and link characteristics of the formation exchange network. These distributed algorithms should rely on a local information-exchange network, relaxing the assumptions on existing algorithms. In this research, it was shown that only local observability is required to design a formation estimator and control law. The approach relies on breaking up the overall information-exchange network into sequence of local subnetworks, and invoking an agreement-type filter to reach consensus among local estimators within each local network. State estimates were obtained by a set of local measurements that were passed through a set of communicating Kalman filters to reach an overall state estimation for the formation. An optimization approach was also presented by means of which diffused estimates over the network can be incorporated in the local estimates obtained by each estimator via local measurements. This approach compares favorably with that obtained by a centralized Kalman filter, which requires complete knowledge of the raw measurement available to each estimator.

Fathpour, Nanaz↗
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↗