Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gene prediction”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Propagation of viral genomes by replicating ammonia-oxidising archaea during soil nitrification

Ammonia-oxidising archaea (AOA) are a ubiquitous component of microbial communities and dominate the first stage of nitrification in some soils. While we are beginning to understand soil virus dynamics, we have no knowledge of the composition or activity of those infecting nitrifiers or their potential to influence processes. This study aimed to characterise viruses having infected autotrophic AOA in two nitrifying soils of contrasting pH by following transfer of assimilated CO 2 -derived 13 C from host to virus via DNA stable-isotope probing and metagenomic analysis. Incorporation of 13 C into low GC mol% AOA and virus genomes increased DNA buoyant density in CsCl gradients but resulted in co-migration with dominant non-enriched high GC mol% genomes, reducing sequencing depth and contig assembly. We therefore developed a hybrid approach where AOA and virus genomes were assembled from low buoyant density DNA with subsequent mapping of 13 C isotopically enriched high buoyant density DNA reads to identify activity of AOA. Metagenome-assembled genomes were different between the two soils and represented a broad diversity of active populations. Sixty-four AOA-infecting viral operational taxonomic units (vOTUs) were identified with no clear relatedness to previously characterised prokaryote viruses. These vOTUs were also distinct between soils, with 42% enriched in 13 C derived from hosts. The majority were predicted as capable of lysogeny and auxiliary metabolic genes included an AOA-specific multicopper oxidase suggesting infection may augment copper uptake essential for central metabolic functioning. These findings indicate virus infection of AOA may be a frequent process during nitrification with potential to influence host physiology and activity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

PKA regulatory subunit Bcy1 couples growth, lipid metabolism, and fermentation during anaerobic xylose growth in Saccharomyces cerevisiae

Past work from our lab, focused on engineering the budding yeast Saccharomyces cerevisiae for fermentation of the non-native pentose sugar xylose, discovered that hyperactivation of the RAS/Protein Kinase A (PKA) pathway was needed for rapid anaerobic xylose fermentation. Interestingly, the mechanism of PKA hyperactivation has a dramatic impact on growth and metabolism on xylose; deletion of the RAS inhibitor IRA2 permits rapid growth and fermentation, while deletion of the PKA regulatory subunit BCY1 allows for fermentation without growth on xylose. To understand how a single deletion in the PKA pathway can decouple growth and metabolism, we performed transcriptomic analysis of these strains, predicting that altered PKA activity would impact global gene expression and identify pathways important for growth and metabolism coordination.

Gasch, Audrey↗
NASA NTRSDate not supplied

Investigating the Relationship between the Cell Wall Integrity Pathway and Unfolded Protein Response

Plants have made significant contributions to astronaut health in spaceflight missions. To further spaceflight research in optimizing plant viability, this study aims to understand the factors involved in maintaining cell wall integrity, which is vital to plant morphology and structural stability. Spaceflight can negatively impact the cell wall; thus, it is crucial to investigate how to mitigate spaceflight stressors to maintain the integrity of the cell wall. The structural integrity of plants’ cell walls depends on secondary cell wall biogenesis, which enables the repair and architectural support of plants like A. thaliana. This biogenesis is triggered by a signal transduction cascade: first initiated by cell wall stress, the CWI (cell wall integrity) pathway is activated, followed by the UPR (unfolded protein response), then the cell wall’s integrity is maintained through secondary cell wall biogenesis. Through a re-analysis of GeneLab Dataset 321 (GLDS-321), a study from NASA’s Open Science Data Repository that investigates the effects of spaceflight on the UPR, several genes were found to be associated with the cell wall. This proposal postulates a relationship between the UPR and the CWI pathway and their direct effect on secondary cell wall biogenesis by investigating IRX7, a gene associated with secondary cell wall biogenesis. The predicted outcome of overexpressing IRX7 is increased resilience of the cell wall by upregulating both the UPR and the CWI pathway, while silencing IRX7 is predicted to compromise the cell wall integrity by downregulating the UPR and the CWI pathway. This study will give insight into the needed measures to increase cell wall resilience in stressful environments: As spaceflight durations increase and uncertain climate change events progress on Earth, understanding how to optimize cell wall resilience – a fundamental pillar of plant health – can effectively enhance mass crop production and quality and ensure the physical and psychological health of astronauts in long-term space missions.

GL4HS↗
DOE OSTI2021

Discovery of widespread transcription initiation at microsatellites predictable by sequence-based deep neural network

Using the Cap Analysis of Gene Expression (CAGE) technology, the FANTOM5 consortium provided one of the most comprehensive maps of transcription start sites (TSSs) in several species. Strikingly, ~72% of them could not be assigned to a specific gene and initiate at unconventional regions, outside promoters or enhancers. Here, we probe these unassigned TSSs and show that, in all species studied, a significant fraction of CAGE peaks initiate at microsatellites, also called short tandem repeats (STRs). To confirm this transcription, we develop Cap Trap RNA-seq, a technology which combines cap trapping and long read MinION sequencing. We train sequence-based deep learning models able to predict CAGE signal at STRs with high accuracy. These models unveil the importance of STR surrounding sequences not only to distinguish STR classes, but also to predict the level of transcription initiation. Importantly, genetic variants linked to human diseases are preferentially found at STRs with high transcription initiation level, supporting the biological and clinical relevance of transcription initiation at STRs. Together, our results extend the repertoire of non-coding transcription associated with DNA tandem repeats and complexify STR polymorphism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Investigating genomic prediction strategies for grain carotenoid traits in a tropical/subtropical maize panel

Abstract Vitamin A deficiency remains prevalent on a global scale, including in regions where maize constitutes a high percentage of human diets. One solution for alleviating this deficiency has been to increase grain concentrations of provitamin A carotenoids in maize (Zea mays ssp. mays L.)—an example of biofortification. The International Maize and Wheat Improvement Center (CIMMYT) developed a Carotenoid Association Mapping panel of 380 inbred lines adapted to tropical and subtropical environments that have varying grain concentrations of provitamin A and other health-beneficial carotenoids. Several major genes have been identified for these traits, 2 of which have particularly been leveraged in marker-assisted selection. This project assesses the predictive ability of several genomic prediction strategies for maize grain carotenoid traits within and between 4 environments in Mexico. Ridge Regression-Best Linear Unbiased Prediction, Elastic Net, and Reproducing Kernel Hilbert Spaces had high predictive abilities for all tested traits (β-carotene, β-cryptoxanthin, provitamin A, lutein, and zeaxanthin) and outperformed Least Absolute Shrinkage and Selection Operator. Furthermore, predictive abilities were higher when using genome-wide markers rather than only the markers proximal to 2 or 13 genes. These findings suggest that genomic prediction models using genome-wide markers (and assuming equal variance of marker effects) are worthwhile for these traits even though key genes have already been identified, especially if breeding for additional grain carotenoid traits alongside β-carotene. Predictive ability was maintained for all traits except lutein in between-environment prediction. The TASSEL (Trait Analysis by aSSociation, Evolution, and Linkage) Genomic Selection plugin performed as well as other more computationally intensive methods for within-environment prediction. The findings observed herein indicate the utility of genomic prediction methods for these traits and could inform their resource-efficient implementation in biofortification breeding programs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A call for caution in the biological interpretation of viral auxiliary metabolic genes

Virus-encoded auxiliary metabolic genes (AMGs) are non-essential genes that increase viral fitness by maintaining or manipulating host metabolism during infection. AMGs are intriguing from an evolutionary perspective, as most viral genomes are highly compact and have limited coding capacity for accessory genes. Advances in viral (meta)genomics have expanded the detection of putative AMGs from viruses in diverse environments. However, this has also led to many instances of misannotation due to the limitations of annotation tools, resulting in misinterpretations about the roles of some viral genes. Here, we highlight studies that support claims about AMGs with more than just function predictions for guidance on best practices. We then propose the adoption of an expanded, inclusive view of all genes auxiliary to core viral functions with the term ‘auxiliary viral genes’ (AVGs), alongside an associated eco-evolutionary framework for considering the types of analyses that can better support claims made about AVGs.

Environmental microbiology↗
DOE OSTI2025

Testing for the Genomic Footprint of Conflict Between Life Stages in an Angiosperm and Moss Species

Abstract The maintenance of genetic variation by balancing selection is of considerable interest to evolutionary biologists. An important but understudied potential driver of balancing selection is antagonistic pleiotropy between diploid and haploid stages of the plant life cycle. Despite sharing a common genome, sporophytes (2n) and gametophytes (n) may undergo differential or even opposing selection. Theoretical work suggests antagonistic pleiotropy between life stages can generate balancing selection and maintain genetic variation. Despite the potential for far-reaching consequences of gametophytic selection, empirical tests of its pleiotropic effects (neutral, synergistic, or antagonistic) on sporophytes are generally lacking. Here, we examined the population genomic signals of selection across life stages in the angiosperm Rumex hastatulus and the moss Ceratodon purpureus. We compared gene expression between life stages and sexes, combined with neutral diversity statistics and the analysis of the distribution of fitness effects. In contrast to what would be predicted under balancing selection due to antagonistic pleiotropy, we found that unbiased genes between life stages were under stronger purifying selection, likely explained by a predominance of synergistic pleiotropy between life stages and strong purifying selection on broadly expressed genes. In addition, we found that 30% of candidate genes under balancing selection in R. hastatulus were located within inversion polymorphisms. Our findings provide novel insights into the genome-wide characteristics and consequences of plant gametophytic selection.

Evolutionary Biology↗
DOE OSTI2021

Multi-Omics Driven Metabolic Network Reconstruction and Analysis of Lignocellulosic Carbon Utilization in Rhodosporidium toruloides

An oleaginous yeast Rhodosporidium toruloides is a promising host for converting lignocellulosic biomass to bioproducts and biofuels. In this work, we performed multi-omics analysis of lignocellulosic carbon utilization in R. toruloides and reconstructed the genome-scale metabolic network of R. toruloides . High-quality metabolic network models for model organisms and orthologous protein mapping were used to build a draft metabolic network reconstruction. The reconstruction was manually curated to build a metabolic model using functional annotation and multi-omics data including transcriptomics, proteomics, metabolomics, and RB-TDNA sequencing. The multi-omics data and metabolic model were used to investigate R. toruloides metabolism including lipid accumulation and lignocellulosic carbon utilization. The developed metabolic model was validated against high-throughput growth phenotyping and gene fitness data, and further refined to resolve the inconsistencies between prediction and data. We believe that this is the most complete and accurate metabolic network model available for R. toruloides to date.

09 BIOMASS FUELS↗
NASA NTRSDate not supplied

The maize rough endosperm6 (rgh6) mutant encodes a predicted DEAD-box RNA helicase

Maize rough endosperm (rgh) mutants have defective seeds with a rough, etched, or pitted endosperm surface. Molecular genetic analysis of this mutant class has identified multiple RNA processing proteins critical to endosperm development. We isolated the rgh6 locus from the UniformMu transposon tagging population. Mutant kernels have reduced grain-fill with defective embryos that fail to germinate. Self-pollination of rgh6 heterozygotes produces mutant seeds at frequency consistent with a single recessive mutation. Molecular mapping localized rgh6 to a 60 kbp interval on chromosome 5. PCR analysis of gene models within the fine-map interval identified a Mutator (Mu) transposon insertion within a predicted DEAD-box RNA helicase gene. A second allele of rgh6 was identified from UniformMu population. Reciprocal crosses between rgh6-umu1 and rgh6-umu2 heterozygotes produces mutant seeds at a frequency consistent with a 3:1 ratio indicating these mutants are allelic. Both rgh6-umu1 and rgh6-umu2 have Mu element insertions at the same insertion site within the DEAD-box RNA helicase gene. However, the terminal inverted repeat sequences are polymorphic between the alleles showing the alleles are independent mutations. We conclude that mutation of the DEAD-box RNA helicase causes the rgh6 phenotype. Further characterization of the rgh6 mutant phenotype and biochemical investigation of the RGH6 protein is expected to give further insights on the roles of RNA processing in endosperm development.

maize↗
DOE OSTI2022

Temporal regulation of cold transcriptional response in switchgrass

Switchgrass low-land ecotypes have significantly higher biomass but lower cold tolerance compared to up-land ecotypes. Understanding the molecular mechanisms underlying cold response, including the ones at transcriptional level, can contribute to improving tolerance of high-yield switchgrass under chilling and freezing environmental conditions. Here, by analyzing an existing switchgrass transcriptome dataset, the temporal cis- regulatory basis of switchgrass transcriptional response to cold is dissected computationally. We found that the number of cold-responsive genes and enriched Gene Ontology terms increased as duration of cold treatment increased from 30 min to 24 hours, suggesting an amplified response/cascading effect in cold-responsive gene expression. To identify genomic sequences likely important for regulating cold response, machine learning models predictive of cold response were established using k -mer sequences enriched in the genic and flanking regions of cold-responsive genes but not non-responsive genes. These k -mers, referred to as putative cis -regulatory elements (pCREs) are likely regulatory sequences of cold response in switchgrass. There are in total 655 pCREs where 54 are important in all cold treatment time points. Consistent with this, eight of 35 known cold-responsive CREs were similar to top-ranked pCREs in the models and only these eight were important for predicting temporal cold response. More importantly, most of the top-ranked pCREs were novel sequences in cold regulation. Our findings suggest additional sequence elements important for cold-responsive regulation previously not known that warrant further studies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Telomere-to-telomere assemblies of chromosome 10 reveal complex adaptive variation of 3-ketoacyl-CoA-synthases in Populus trichocarpa likely driven by Helitrons

The model woody plant Populus trichocarpa displays an atypical alkene-diverse wax cuticle likely driven by copy number variation (CNV) of 3-ketoacyl-CoA synthases ( KCS ), which has been difficult to confirm with short-read assemblies. Long-read sequencing enables the development of telomere-to-telomere resources to detect cryptic variation, including CNVs, which are currently missed. Integrating this information can improve genomic prediction for breeding and provide insights into the evolutionary basis of important traits. Our analysis of 78 long-read haplotypes from chromosome 10 identified more than twice as many KCS genes as previously reported, and numerous intragenic non-synonymous substitutions. Random Forest predictive models highlighted the importance of Potri.010G079500 in producing very long chain alkenes; however, its absence did not predict previously reported alkene-deficient phenotypes. Instead, alkene levels are best predicted by the combinations of KCS copies. Additionally, amino acid substitutions clustered around ligand and donor binding pockets, suggesting they contribute to differing wax cuticle composition. Finally, each KCS gene and copy was linked to a Helitron transposon. A phylogenetic analysis suggests Helitrons are the evolutionary mechanism for generating KCS tandem arrays. Long-read generated telomere-to-telomere assemblies of P. trichocarpa chromosome 10 revealed large-effect loci critical to genetic studies that are unattainable from short-reads. This new resource produced novel insights into genome structure and function, and a novel mechanism for generating tandem gene duplication. Our results highlight that, given current challenges in annotation and assembly, detailed and focused long-read sequences are key to interpreting complex genomic regions that contain tandem copy number variants.

09 BIOMASS FUELS↗
NASA NTRS2005

Human Adaptation Genetic Response Suites: Toward New Interventions and Countermeasures for Spaceflight

Genetic response suites in human lymphocytes in response to microgravity are important to identify and further study in order to augment human physiological adaptation to novel environments. Emerging technologies, such as DNA micro array profiling, have the potential to identify novel genes that are involved in mediating adaptation to these environments. These genes may prove to be therapeutically valuable as new targets for countermeasures, or as predictive biomarkers of response to these new environments. Human lymphocytes cultured in lg and microgravity analog culture were analyzed for their differential gene expression response. Different groups of genes related to the immune response, cardiovascular system and stress response were then analyzed. Analysis of cells from multiple donors reveals a small shared set that are likely to be essential to adaptation. These three groups focus on human adaptation to new environments. The shared set contains genes related to T cell activation, immune response and stress response to analog microgravity.

Sundaresan, A.↗
DOE OSTI2021

Soil enzymes as indicators of soil function: A step toward greater realism in microbial ecological modeling

Soil carbon (C) and nitrogen (N) cycles and their complex responses to environmental changes have received increasing attention. However, large uncertainties in model predictions remain, partially due to the lack of explicit representation and parameterization of microbial processes. One great challenge is to effectively integrate rich microbial functional traits into ecosystem modeling for better predictions. Here, using soil enzymes as indicators of soil function, we developed a competitive dynamic enzyme allocation scheme and detailed enzyme-mediated soil inorganic N processes in the Microbial-ENzyme Decomposition (MEND) model. We conducted a rigorous calibration and validation of MEND with diverse soil C-N fluxes, microbial C:N ratios, and functional gene abundances from a 12-year CO 2 × N grassland experiment (BioCON) in Minnesota, USA. In addition to accurately simulating soil CO 2 fluxes and multiple N variables, the model correctly predicted microbial C:N ratios and their negative response to enriched N supply. Model validation further showed that, compared to the changes in simulated enzyme concentrations and decomposition rates, the changes in simulated activities of eight C-N-associated enzymes were better explained by the measured gene abundances in responses to elevated atmospheric CO 2 concentration. In conclusion, our results demonstrated that using enzymes as indicators of soil function and validating model predictions with functional gene abundances in ecosystem modeling can provide a basis for testing hypotheses about microbially mediated biogeochemical processes in response to environmental changes. Further development and applications of the modeling framework presented here will enable microbial ecologists to address ecosystem-level questions beyond empirical observations, toward more predictive understanding, an ultimate goal of microbial ecology.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2017

Systemic Response to Microgravity: Utilizing GeneLab Datasets to Identify Molecular Targets for Future Hypotheses-Driven Spaceflight Studies

Biological risks associated with microgravity are a major concern for long-term space travel. Although determination of risk has been a focus for NASA research, data examining systemic (i.e., multi- or pan-tissue) responses to space flight are sparse. To perform our analysis, we utilized the NASA GeneLab database which is a publicly available repository containing a wide array of omics results from experiments conducted with: i) with different flight conditions (space shuttle (STS) missions vs. International Space Station (ISS); ii) a variety of tissues; and 3) assays that measure epigenetic, transcriptional, and protein expression changes. Meta-analysis of the transcriptomic data from 7 different murine and rat data sets, examining tissues such as liver, kidney, adrenal gland, thymus, mammary gland, skin, and skeletal muscle (soleus, extensor digitorum longus, tibialis anterior, quadriceps, and gastrocnemius) revealed for the first time, the existence of potential master regulators coordinating systemic responses to microgravity in rodents. We identified p53, TGF(beta)1 and immune related pathways as the highly prevalent pan-tissue signaling pathways that are affected by microgravity. Some variability in the degree of change in their expression across species, strain and time of flight was also observed. Interestingly, while certain skeletal muscle (gastrocnemius and soleus) exhibited an overall down-regulation of these genes, some other muscle types such as the extensor digitorum longus, tibialis anterior and quadriceps, showed an up-regulated expression, indicative of potential compensatory mechanisms to prevent microgravity-induced atrophy. Key genes isolated by unbiased systems analyses displayed a major overlap between tissue types and flight conditions and established TGF(beta)1 to be the most connected gene across all data sets. Finally, a set of microgravity responsive miRNA signature was identified and based on their predicted functional state and subsequent impact on health, a theoretical health risk score was calculated. The genes and miRNAs identified from our analyses can be targeted for future research involving efficient countermeasure design. Our study thus exemplifies the utility of GeneLab data repository to aid in the process of performing novel hypothesis based spaceflight research aimed at elucidating the global impact of environmental stressors at multiple biological scales.

GeneLab↗
NASA NTRS2017

Systemic Microgravity Response: Utilizing GeneLab to Develop Hypotheses for Spaceflight Risks

Biological risks associated with microgravity are a major concern for long-term space travel. Although determination of risk has been a focus for NASA research, data examining systemic (i.e., multi- or pan-tissue) responses to space flight are sparse. To perform our analysis, we utilized the NASA GeneLab database which is a publicly available repository containing a wide array of omics results from experiments conducted with: i) with different flight conditions (space shuttle (STS) missions vs. International Space Station (ISS); ii) a variety of tissues; and 3) assays that measure epigenetic, transcriptional, and protein expression changes. Meta-analysis of the transcriptomic data from 7 different murine and rat data sets, examining tissues such as liver, kidney, adrenal gland, thymus, mammary gland, skin, and skeletal muscle (soleus, extensor digitorum longus, tibialis anterior, quadriceps, and gastrocnemius) revealed for the first time, the existence of potential master regulators coordinating systemic responses to microgravity in rodents. We identified p53, TGF1 and immune related pathways as the highly prevalent pan-tissue signaling pathways that are affected by microgravity. Some variability in the degree of change in their expression across species, strain and time of flight was also observed. Interestingly, while certain skeletal muscle (gastrocnemius and soleus) exhibited an overall down-regulation of these genes, some other muscle types such as the extensor digitorum longus, tibialis anterior and quadriceps, showed an up-regulated expression, indicative of potential compensatory mechanisms to prevent microgravity-induced atrophy. Key genes isolated by unbiased systems analyses displayed a major overlap between tissue types and flight conditions and established TGF1 to be the most connected gene across all data sets. Finally, a set of microgravity responsive miRNA signature was identified and based on their predicted functional state and subsequent impact on health, a theoretical health risk score was calculated. The genes and miRNAs identified from our analyses can be targeted for future research involving efficient countermeasure design. Our study thus exemplifies the utility of GeneLab data repository to aid in the process of performing novel hypothesis based spaceflight research aimed at elucidating the global impact of environmental stressors at multiple biological scales.

GeneLab↗
DOE OSTI2021

Transcription factor FOXC1 positively regulates SFRP1 expression in androgenetic alopecia

Androgenetic alopecia (AGA) is the most common type of hair loss dysfunction. Secreted frizzled related protein 1 (SFRP1) is found to be associated with hair loss, but its role in AGA and the regulation mechanism of its transcription level is unclear. The aim of our study is to explore the expression of SFRP1 in AGA samples and its transcriptional mechanism. Male frontal and occipital scalp hair follicles from AGA patients were collected, and human dermal papilla cells (DPCs) were isolated and cultured. SFRP1 gene was cloned and constructed into recombinant plasmids to perform dual-luciferase reporter assay. Transcription factor binding sites were predicted through the Jaspar website and further confirmed by the chromatin immunoprecipitation (ChIP) assay. Expression of genes in DPCs was determined by immunofluorescence (IF) staining, quantitative real-time PCR (qRT-PCR) and western blotting. Our findings showed that SFRP1 was highly expressed in DPCs of AGA patients. The core promoter region of SFRP1 was from −100 to +50 bp and was found to be positively regulated by forkhead box C1 (FOXC1), a transcription factor related to hair growth, both at mRNA and protein level in DPCs. Our study suggests that FOXC1 plays an important role in regulating SFRP1 transcription, which may provide new insights into the development of therapeutic strategies for the treatment of AGA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Integration of Computational Docking into Anti-Cancer Drug Response Prediction Models

Cancer is a heterogeneous disease in that tumors of the same histology type can respond differently to a treatment. Anti-cancer drug response prediction is of paramount importance for both drug development and patient treatment design. Although various computational methods and data have been used to develop drug response prediction models, it remains a challenging problem due to the complexities of cancer mechanisms and cancer-drug interactions. To better characterize the interaction between cancer and drugs, we investigate the feasibility of integrating computationally derived features of molecular mechanisms of action into prediction models. Specifically, we add docking scores of drug molecules and target proteins in combination with cancer gene expressions and molecular drug descriptors for building response models. The results demonstrate a marginal improvement in drug response prediction performance when adding docking scores as additional features, through tests on large drug screening data. We discuss the limitations of the current approach and provide the research community with a baseline dataset of the large-scale computational docking for anti-cancer drugs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Integrating genomics and metabolomics for scalable non-ribosomal peptide discovery

Non-Ribosomal Peptides (NRPs) represent a biomedically important class of natural products that include a multitude of antibiotics and other clinically used drugs. NRPs are not directly encoded in the genome but are instead produced by metabolic pathways encoded by biosynthetic gene clusters (BGCs). Since the existing genome mining tools predict many putative NRPs synthesized by a given BGC, it remains unclear which of these putative NRPs are correct and how to identify post-assembly modifications of amino acids in these NRPs in a blind mode, without knowing which modifications exist in the sample. To address this challenge, here we report NRPminer, a modification-tolerant tool for NRP discovery from large (meta)genomic and mass spectrometry datasets. We show that NRPminer is able to identify many NRPs from different environments, including four previously unreported NRP families from soil-associated microbes and NRPs from human microbiota. Furthermore, in this work we demonstrate the anti-parasitic activities and the structure of two of these NRP families using direct bioactivity screening and nuclear magnetic resonance spectrometry, illustrating the power of NRPminer for discovering bioactive NRPs.

59 BASIC BIOLOGICAL SCIENCES↗