Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transcriptome”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Integrated Transcriptomic and Proteomic Analysis Identifies Plasma Biomarkers of Hepatocellular Failure in Alcohol-Associated Hepatitis

Alcohol-associated hepatitis (AH) is a form of liver failure with high short-term mortality. Recent results have shown that HNF4a defective function and systemic inflammation are major disease drivers of AH. Plasma biomarkers of hepatocyte function could be useful for diagnostic and prognostic purposes. Herein an integrative analysis of hepatic RNAseq and liquid chromatography-tandem mass spectrometry (LC-MS/MS) was performed to identify plasma protein signatures for mild and severe AH patients. Alcohol-related liver disease cirrhosis (ALD)(AC), non-alcoholic fatty liver disease (NALFD), and healthy subjects (HC) were used as comparator groups. Identified proteins primarily involved in hepatocellular function were decreased in AH patients which included hepatokines, clotting factors, complement cascade components, and hepatocyte growth activators. A protein signature of AH disease severity was identified including thrombin (THRB), hepatocyte growth factor alpha (HGFA), clusterin (CLUS), human serum factor H-related protein (FHR1) and kallistatin (KAIN), which exhibited large abundance shifts between severe and non-severe AH. The combination of THRB and HGFA discriminated between severe and non-severe AH with high sensitivity and specificity. These findings were correlated with the liver expression of genes encoding secreted proteins in a similar cohort, finding a highly consistent plasma protein signature reflecting HNF4A and HNF1A functions. This unbiased proteomic-transcriptome analysis identified plasma protein signatures and pathways associated with disease severity, reflecting HNF4A/1A activity useful for diagnostic assessment in AH.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Comparative transcriptomics of CAZy enzymes in white- and brown-rot agaricomycetes: Evolutionary insights into lignocellulose degradation and the relevance of GH16 glycoside hydrolase functional divergence

White-rot and brown-rot fungi (WRF and BRF, respectively) decompose lignocellulose, the main structural component of plant biomass, through distinct mechanisms. This study examines the transcriptomic responses of three WRF species (Pleurotus ostreatus, Phanerochaete chrysosporium, and Heterobasidion irregulare) and two BRF species (Fomitopsis schrenkii and Rhodonia placenta) grown on poplar wood (W) and glucose (G) as sole carbon sources. RNA-seq analysis revealed upregulation of carbohydrate-active enzymes (CAZymes) linked to lignocellulose degradation. WRF displayed a broader enzymatic repertoire, whereas BRF employed a more selective approach. Among these responses, GH16 glycoside hydrolases were consistently upregulated across all species, including BRF. Since GH16 enzymes are involved in both plant hemicellulose modification and cell wall remodeling, their wood-induced expression may reflect multiple processes rather than a single conserved wood-decay mechanism. Structural and phylogenetic analyses revealed species-specific divergence, consistent with potential functional specialization. These findings broaden our understanding of fungal enzymatic strategies and highlight GH16 enzymes as candidates for further study. Beyond the immediate context of wood decay, this work has broader implications for fungal ecology, evolutionary biology, and biotechnological applications such as biomass conversion and sustainable bioenergy.

Fungal ecological strategies↗
DOE OSTI2023

Learning perturbation-inducible cell states from observability analysis of transcriptome dynamics

Abstract A major challenge in biotechnology and biomanufacturing is the identification of a set of biomarkers for perturbations and metabolites of interest. Here, we develop a data-driven, transcriptome-wide approach to rank perturbation-inducible genes from time-series RNA sequencing data for the discovery of analyte-responsive promoters. This provides a set of biomarkers that act as a proxy for the transcriptional state referred to as cell state. We construct low-dimensional models of gene expression dynamics and rank genes by their ability to capture the perturbation-specific cell state using a novel observability analysis. Using this ranking, we extract 15 analyte-responsive promoters for the organophosphate malathion in the underutilized host organism Pseudomonas fluorescens SBW25. We develop synthetic genetic reporters from each analyte-responsive promoter and characterize their response to malathion. Furthermore, we enhance malathion reporting through the aggregation of the response of individual reporters with a synthetic consortium approach, and we exemplify the library’s ability to be useful outside the lab by detecting malathion in the environment. The engineered host cell, a living malathion sensor, can be optimized for use in environmental diagnostics while the developed machine learning tool can be applied to discover perturbation-inducible gene expression systems in the compendium of host organisms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Deciphering triterpenoid saponin biosynthesis by leveraging transcriptome response to methyl jasmonate elicitation in Saponaria vaccaria

Abstract Methyl jasmonate (MeJA) is a known elicitor of plant specialized metabolism, including triterpenoid saponins. Saponaria vaccaria is an annual herb used in traditional Chinese medicine, containing large quantities of oleanane-type triterpenoid saponins with anticancer properties and structural similarities to the vaccine adjuvant QS-21. Leveraging the MeJA-elicited saponin biosynthesis, we identify multiple enzymes catalyzing the oxidation and glycosylation of triterpenoids in S. vaccaria . This exploration is aided by Pacbio full-length transcriptome sequencing and gene expression analysis. A cellulose synthase-like enzyme can not only glucuronidate triterpenoid aglycones but also alter the product profile of a cytochrome P450 monooxygenase via preference for the aldehyde intermediate. Furthermore, the discovery of a UDP-glucose 4,6-dehydratase and a UDP-4-keto-6-deoxy-glucose reductase reveals the biosynthetic pathway for the rare nucleotide sugar UDP- d -fucose, a likely sugar donor for fucosylation of plant natural products. Our work enables the production and optimization of high-value saponins in microorganisms and plants through synthetic biology approaches.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A multi-scale pipeline linking drug transcriptomics with pharmacokinetics predicts in vivo interactions of tuberculosis drugs

Tuberculosis (TB) is the deadliest infectious disease worldwide. The design of new treatments for TB is hindered by the large number of candidate drugs, drug combinations, dosing choices, and complex pharmaco-kinetics/dynamics (PK/PD). Here we study the interplay of these factors in designing combination therapies by linking a machine-learning model, INDIGO-MTB, which predicts in vitro drug interactions using drug transcriptomics, with a multi-scale model of drug PK/PD and pathogen-immune interactions called GranSim. We calculate an in vivo drug interaction score (iDIS) from dynamics of drug diffusion, spatial distribution, and activity within lesions against various pathogen sub-populations. The iDIS of drug regimens evaluated against non-replicating bacteria significantly correlates with efficacy metrics from clinical trials. Our approach identifies mechanisms that can amplify synergistic or mitigate antagonistic drug interactions in vivo by modulating the relative distribution of drugs. Our mechanistic framework enables efficient evaluation of in vivo drug interactions and optimization of combination therapies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Comparative transcriptome analysis of cells from different areas reveals ROS responsive mechanism at sclerotial initiation stage in Morchella importuna

Morels are some of the most highly prized edible and medicinal mushrooms, with great economic and scientific value. Outdoor cultivation has been achieved and expanded on a large scale in China in recent years. Sclerotial formation is one of the most important phases during the morel life cycle, and previous reports indicated that reactive oxygen species (ROS) play an important role. However, ROS response mechanisms at sclerotial initiation (SI) stage are poorly understood. In this study, comparative transcriptome analyses were performed with sclerotial and hyphal cells at different areas in the same plate at SI stage. Gene expression was significantly different at SI stage between sclerotial formation and mycelia growth areas. GO and KEGG analyses indicated more vigorous metabolic characteristics in the hyphae area, while transcription process, DNA repair, and protein processing were enriched in sclerotial cells. Gene expression related to H 2 O 2 production was high in the hyphae area, while expression of H 2 O 2 -scavenging genes was high in sclerotial cells, leading to a higher H 2 O 2 concentration in the hyphal region than in the sclerotium. Minor differences were observed in gene expression of H 2 O 2 -induced signaling pathway in sclerotial and hyphal cells; however, expression levels of the target genes of transcription factor MSN2, important in the H 2 O 2 -induced signaling pathways, were significantly different. MSN2 enhanced stress response regulation in sclerotia by regulating these target genes. Small molecular HSPs were also found upregulated in sclerotial cells. This study indicated that sclerotial cells are more resistant to ROS stress than hyphal cells through transcriptional regulation of related genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Integration of high-density genetic mapping with transcriptome analysis uncovers numerous agronomic QTL and reveals candidate genes for the control of tillering in sorghum

Phenotypes such as branching, photoperiod sensitivity, and height were modified during plant domestication and crop improvement. Here, we perform quantitative trait locus (QTL) mapping of these and other agronomic traits in a recombinant inbred line (RIL) population derived from an interspecific cross between Sorghum propinquum and Sorghum bicolor inbred Tx7000. Using low-coverage Illumina sequencing and a bin-mapping approach, we generated ~1920 bin markers spanning ~875 cM. Phenotyping data were collected and analyzed from two field locations and one greenhouse experiment for six agronomic traits, thereby identifying a total of 30 QTL. Many of these QTL were penetrant across environments and co-mapped with major QTL identified in other studies. Other QTL uncovered new genomic regions associated with these traits, and some of these were environment-specific in their action. To further dissect the genetic underpinnings of tillering, we complemented QTL analysis with transcriptomics, identifying 6189 genes that were differentially expressed during tiller bud elongation. We identified genes such as Dormancy Associated Protein 1 (DRM1) in addition to various transcription factors that are differentially expressed in comparisons of dormant to elongating tiller buds and lie within tillering QTL, suggesting that these genes are key regulators of tiller elongation in sorghum. Our study demonstrates the usefulness of this RIL population in detecting domestication and improvement-associated genes in sorghum, thus providing a valuable resource for genetic investigation and improvement to the sorghum community.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Cuticle development and the underlying transcriptome–metabolome associations during early seedling establishment

Abstract The plant cuticle is a complex extracellular lipid barrier that has multiple protective functions. This study investigated cuticle deposition by integrating metabolomics and transcriptomics data gathered from six different maize seedling organs of four genotypes, the inbred lines B73 and Mo17, and their reciprocal hybrids. These datasets captured the developmental transition of the seedling from heterotrophic skotomorphogenic growth to autotrophic photomorphogenic growth, a transition that is highly vulnerable to environmental stresses. Statistical interrogation of these data revealed that the predominant determinant of cuticle composition is seedling organ type, whereas the seedling genotype has a smaller effect on this phenotype. Gene-to-metabolite associations assessed by integrated statistical analyses identified three gene networks associated with the deposition of different elements of the cuticle: cuticular waxes; monomers of lipidized cell wall biopolymers, including cutin and suberin; and both of these elements. These gene networks reveal three metabolic programs that appear to support cuticle deposition, including processes of chloroplast biogenesis, lipid metabolism, and molecular regulation (e.g. transcription factors, post-translational regulators, and phytohormones). This study demonstrates the wider physiological metabolic context that can determine cuticle deposition and lays the groundwork for new targets for modulating the properties of this protective barrier.

Plant Sciences↗
DOE OSTI2023

Canalization of Phenotypes—When the Transcriptome is Constantly but Weakly Perturbed

Abstract Recent studies have increasingly pointed to microRNAs (miRNAs) as the agent of gene regulatory network (GRN) stabilization as well as developmental canalization against constant but small environmental perturbations. To analyze mild perturbations, we construct a Dicer-1 knockdown line (dcr-1 KD) in Drosophila that modestly reduces all miRNAs by, on average, ∼20%. The defining characteristic of stabilizers is that, when their capacity is compromised, GRNs do not change their short-term behaviors. Indeed, even with such broad reductions across all miRNAs, the changes in the transcriptome are very modest during development in stable environment. By comparison, broad knockdowns of other regulatory genes (esp. transcription factors) by the same method should lead to drastic changes in the GRNs. The consequence of destabilization may thus be in long-term development as postulated by the theory of canalization. Flies with modest miRNA reductions may gradually deviate from the developmental norm, resulting in late-stage failures such as shortened longevity. In the optimal culture condition, the survival to adulthood is indeed normal in the dcr-1 KD line but, importantly, adult longevity is reduced by ∼90%. When flies are stressed by high temperature, dcr-1 KD induces lethality earlier in late pupation and, as the perturbations are shifted earlier, the affected stages are shifted correspondingly. Hence, in late stages of development with deviations piling up, GRN would be increasingly in need of stabilization. In conclusion, miRNAs appear to be a solution to weak but constant environmental perturbations.

Lu, Guang-An↗
DOE OSTI2023

Transcriptomic Analysis Reveals Novel Regulators of the Scots Pine Stilbene Pathway

Stilbenes accumulate in Scots pine heartwood where they have important roles in protecting wood from decaying fungi. They are also part of active defense responses, and their production is induced by different (a)biotic stressors. The specific transcriptional regulators as well as the enzyme responsible for activating the stilbene precursor cinnamate in the pathway are still unknown. UV-C radiation was the first discovered artificial stress activator of the pathway. Here, we describe a large-scale transcriptomic analysis of pine needles in response to UV-C and treatment with translational inhibitors, both activating the transcription of stilbene pathway genes. We used the data to identify putative candidates for the missing CoA ligase and for pathway regulators. We further showed that the pathway is transcriptionally activated by phosphatase inhibitor, ethylene and jasmonate treatments, as in grapevine, and that the stilbene synthase promoter retains its inducibility in some of the tested conditions in Arabidopsis, a species that normally does not synthesize stilbenes. Shared features between gymnosperm and angiosperm regulation and partially retained inducibility in Arabidopsis suggest that pathway regulation occurs not only via ancient stress-response pathway(s) but also via species-specific regulators. Understanding which genes control the biosynthesis of stilbenes in Scots pine aids breeding of more resistant trees.

4-Coumarate-CoA ligase↗
DOE OSTI2024

Signature analysis of high-throughput transcriptomics screening data for mechanistic inference and chemical grouping

Abstract High-throughput transcriptomics (HTTr) uses gene expression profiling to characterize the biological activity of chemicals in in vitro cell-based test systems. As an extension of a previous study testing 44 chemicals, HTTr was used to screen an additional 1,751 unique chemicals from the EPA’s ToxCast collection in MCF7 cells using 8 concentrations and an exposure duration of 6 h. We hypothesized that concentration-response modeling of signature scores could be used to identify putative molecular targets and cluster chemicals with similar bioactivity. Clustering and enrichment analyses were conducted based on signature catalog annotations and ToxPrint chemotypes to facilitate molecular target prediction and grouping of chemicals with similar bioactivity profiles. Enrichment analysis based on signature catalog annotation identified known mechanisms of action (MeOAs) associated with well-studied chemicals and generated putative MeOAs for other active chemicals. Chemicals with predicted MeOAs included those targeting estrogen receptor (ER), glucocorticoid receptor (GR), retinoic acid receptor (RAR), the NRF2/KEAP/ARE pathway, AP-1 activation, and others. Using reference chemicals for ER modulation, the study demonstrated that HTTr in MCF7 cells was able to stratify chemicals in terms of agonist potency, distinguish ER agonists from antagonists, and cluster chemicals with similar activities as predicted by the ToxCast ER Pathway model. Uniform manifold approximation and projection (UMAP) embedding of signature-level results identified novel ER modulators with no ToxCast ER Pathway model predictions. Finally, UMAP combined with ToxPrint chemotype enrichment was used to explore the biological activity of structurally related chemicals. The study demonstrates that HTTr can be used to inform chemical risk assessment by determining in vitro points of departure, predicting chemicals’ MeOA and grouping chemicals with similar bioactivity profiles.

Toxicology↗
DOE OSTI2021

Serendipita Fungi Modulate the Switchgrass Root Transcriptome to Circumvent Host Defenses and Establish a Symbiotic Relationship

The fungal family Serendipitaceae encompasses root-associated lineages with endophytic, ericoid, orchid, and ectomycorrhizal lifestyles. Switchgrass is an important bioenergy crop for cellulosic ethanol production owing to high biomass production on marginal soils otherwise unfit for food crop cultivation. The aim of this study was to investigate the host plant responses to Serendipita spp. colonization by characterizing the switchgrass root transcriptome during different stages of symbiosis in vitro. For this, we included a native switchgrass strain, Serendipita bescii, and a related strain, S. vermifera, isolated from Australian orchids. Serendipita colonization progresses from thin hyphae that grow between root cells to, finally, the production of large, bulbous hyphae that fill root cells during the later stages of colonization. We report that switchgrass seems to perceive both fungi prior to physical contact, leading to the activation of chemical and structural defense responses and putative host disease resistance genes. Subsequently, the host defense system appears to be quenched and carbohydrate metabolism adjusted, potentially to accommodate the fungal symbiont. In addition, prior to contact, switchgrass exhibited significant increases in root hair density and root surface area. Furthermore, genes involved in phytohormone metabolism such as gibberellin, jasmonic acid, and salicylic acid were activated during different stages of colonization. Both fungal strains induced plant gene expression in a similar manner, indicating a conserved plant response to members of this fungal order. Understanding plant responsiveness to Serendipita spp. will inform our efforts to integrate them into forages and row crops for optimal plant–microbe functioning, thus facilitating low-input, sustainable agricultural practices.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Comprehensive analyses of 723 transcriptomes enhance genetic and biological interpretations for complex traits in cattle

By uniformly analyzing 723 RNA-seq data from 91 tissues and cell types, we built a comprehensive gene atlas and studied tissue specificity of genes in cattle. We demonstrated that tissue-specific genes significantly reflected the tissue-relevant biology, showing distinct promoter methylation and evolution patterns (e.g., brain-specific genes evolve slowest, whereas testisspecific genes evolve fastest). Through integrative analyses of those tissue-specific genes with large-scale genome-wide association studies, we detected relevant tissues/cell types and candidate genes for 45 economically important traits in cattle, including blood/immune system (e.g., CCDC88C) for male fertility, brain (e.g., TRIM46 and RAB6A) for milk production, and multiple growth-related tissues (e.g., FGF6 and CCND2) for body conformation. We validated these findings by using epigenomic data across major somatic tissues and sperm. Collectively, our findings provided novel insights into the genetic and biological mechanisms underlying complex traits in cattle, and our transcriptome atlas can serve as a primary source for biological interpretation, functional validation, studies of adaptive evolution, and genomic improvement in livestock.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Serendipita bescii promotes winter wheat growth and modulates the host root transcriptome under phosphorus and nitrogen starvation

Summary Serendipita vermifera ssp. bescii , hereafter referred to as S. bescii , is a root‐associated fungus that promotes plant growth in both its native switchgrass host and a variety of monocots and dicots. Winter wheat ( Triticum aestivum L.), a dual‐purpose crop, used for both forage and grain production, significantly contributes to the agricultural economies of the Southern Great Plains, USA. In this study, we investigated the influence of S. bescii on growth and transcriptome regulation of nitrogen (N) and phosphorus (P) metabolism in winter wheat. Serendipita bescii significantly improved lateral root growth and forage biomass under a limited N or P regime. Further, S. bescii activated sets of host genes regulating N and P starvation responses. These genes include, root‐specific auxin transport, strigolactone and gibberellin biosynthesis, degradation of phospholipids and biosynthesis of glycerolipid, downregulation of ammonium transport and nitrate assimilation, restriction of protein degradation by autophagy and subsequent N remobilization. All these genes are hypothesized to regulate acquisition, assimilation and remobilization of N and P. Based on transcriptional level gene regulation and physiological responses to N or P limitation, we suggest S. bescii plays a critical role in modulating stress imposed by limitation of these two critical nutrients in winter wheat.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Gene family expansions and transcriptome signatures uncover fungal adaptations to wood decay

Because they comprise some of the most efficient wood-decayers, Polyporales fungi impact carbon cycling in forest environment. Despite continuous discoveries on the enzymatic machinery involved in wood decomposition, the vision on their evolutionary adaptation to wood decay and genome diversity remains incomplete. We combined the genome sequence information from 50 Polyporales species, including 26 newly sequenced genomes and sought for genomic and functional adaptations to wood decay through the analysis of genome composition and transcriptome responses to different carbon sources. The genomes of Polyporales from different phylogenetic clades showed poor conservation in macrosynteny, indicative of genome rearrangements. We observed different gene family expansion/contraction histories for plant cell wall degrading enzymes in core polyporoids and phlebioids and captured expansions for genes involved in signalling and regulation in the lineages of white rotters. Furthermore, we identified conserved cupredoxins, thaumatin-like proteins and lytic polysaccharide monooxygenases with a yet uncharacterized appended module as new candidate players in wood decomposition. Given the current need for enzymatic toolkits dedicated to the transformation of renewable carbon sources, the observed genomic diversity among Polyporales strengthens the relevance of mining Polyporales biodiversity to understand the molecular mechanisms of wood decay.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Comparative transcriptomics uncovers poplar and fungal genetic determinants of ectomycorrhizal compatibility

Ectomycorrhizal symbiosis supports tree growth and is crucial for nutrient cycling and temperate and boreal ecosystems functioning. The establishment of functional ectomycorrhiza (ECM) first requires the association of compatible partners. However, host and fungal genetic determinants governing mycorrhizal compatibility are unknown. To identify such factors in poplar and its fungal associates, we mined existing and de novo tree and fungal transcriptional datasets. We identified co-expressed genes enabling ECM symbiosis at early and mature stages of the interaction. These sets of genes can be divided into general fungal-sensing and ECM-specific components. We highlight the importance of fungal modulation of plant JA-related defenses and the regulation of secretory pathways for ECM compatibility, including upregulation of key fungal small secreted proteins, the downregulation of plant secreted peroxidases, and the downregulation of plant cell wall remodeling proteins concomitantly with the upregulation of fungal glycosyl hydrolases acting on pectin. Not only gene regulation, but also its temporal scale and dynamics seem to play a crucial role for mycorrhizal compatibility. The expression profile of the host Common Symbiosis Pathway and nutrient transporters was also studied, revealing constitutive levels of expression and moderate upregulation in compatible ECM interactions. Overall, these results underscore the importance of novel biological functions during the establishment of ECM symbiosis, help us gain insights into the molecular events determining mycorrhiza compatibility, and serve as a data-rich transcriptomic resource to open new research questions in the field.

Marqués‐Gálvez, José Eduardo↗
DOE OSTI2022

Unified epigenomic, transcriptomic, proteomic, and metabolomic taxonomy of Alzheimer’s disease progression and heterogeneity

Alzheimer’s disease (AD) is a heterogeneous disorder with abnormalities in multiple biological domains. In an advanced machine learning analysis of postmortem brain and in vivo blood multi-omics molecular data ( N = 1863), we integrated epigenomic, transcriptomic, proteomic, and metabolomic profiles into a multilevel biological AD taxonomy. We obtained a personalized multilevel molecular index of AD dementia progression that predicts severity of neuropathologies, and identified three robust molecular-based subtypes that explain much of the pathologic and clinical heterogeneity of AD. These subtypes present distinct patterns of alteration in DNA methylation, RNA, proteins, and metabolites, identifiable in the brain and subsequently in blood. In addition, the genetic variations that predispose to the various AD subtypes in brain predict distinct spatial patterns of alteration in cell types, suggesting a unique influence of each putative AD variant on neuropathological mechanisms. These observations support that an individually tailored multi-omics molecular taxonomy of AD may represent distinct targets for preventive or treatment interventions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Transcriptomic data sets for Novosphingobium aromaticivorans grown with the β-5-linked aromatic dimer dehydrodiconiferyl alcohol and the related G-aromatic monomers vanillin and ferulic acid

ABSTRACT The transcriptomes of a 2-pyrone-4,6-dicarboxylic acid-producing strain of Novosphingobium aromaticivorans DSM12444 were determined when grown in minimal medium containing glucose alone or glucose plus vanillin, ferulic acid, or the β-5-linked aromatic dimer dehydrodiconiferyl alcohol as carbon sources. Here, we present the RNA-sequencing data we obtained.

Metz, Fletcher↗