Search NASA⌕ Search

SEARCH · Search NASA

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 253 records · Page 14

Capturing an Early Gene Induction Event during Wood Decay by the Brown Rot Fungus Rhodonia placenta

Brown rot fungi dominate wood decomposition in coniferous forests, and their carbohydrate-selective mechanisms are of commercial interest. Brown rot was recently described as a two-step, sequential mechanism orchestrated by fungi using differentially expressed genes (DEGs) and consisting of oxidation via reactive oxygen species (ROS) followed by enzymatic saccharification. There have been indications, however, that the initial oxidation step itself might require induction. To capture this early gene regulation event, here, we integrated fine-scale cryosectioning with whole-transcriptome sequencing to dissect gene expression at the single-hyphal-cell scale (tens of micrometers). This improved the spatial resolution 50-fold, relative to previous work, and we were able to capture the activity of the first 100 μm of hyphal front growth by Rhodonia placenta in aspen wood. This early decay period was dominated by delayed gene expression patterns as the fungus ramped up its mechanism. These delayed DEGs included many genes implicated in ROS pathways (lignocellulose oxidation [LOX]) that were previously and incorrectly assumed to be constitutively expressed. These delayed DEGs, which include those with and without predicted functions, also create a focused subset of target genes for functional genomics. However, this delayed pattern was not universal, with a few genes being upregulated immediately at the hyphal front. Most notably, this included a gene commonly implicated in hydroquinone and iron redox cycling: benzoquinone reductase.

59 BASIC BIOLOGICAL SCIENCES↗

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

ABSTRACT The Shox2 homeodomain transcriptional regulator is known for its critical functions during mouse embryogenesis, enabling accurate development of limbs, craniofacial structures, neural populations and the cardiac conduction system. At the genomic level, the Shox2 gene is flanked by an extensive gene desert, a continuous non-coding genomic region spanning over 500 kilobases that contains a multitude of evolutionarily conserved elements with predicted cis -regulatory activities. However, the transcriptional enhancer potential of the vast majority of these elements in combination with the biological necessity of the gene desert have not yet been explored. Using transgenic reporter assays in mouse embryos to validate an extensive set of stringent epigenomic enhancer predictions, we identify several novel gene desert enhancers with distinct tissue-specific activities in Shox2 expressing tissues. 4C-seq chromatin conformation capture further uncovers a repertoire of gene desert enhancers with overlapping activities in the proximal limb, in a compartment essential for Shox2 -mediated stylopod formation. Leveraging CRISPR/Cas9 to delete the gene desert region contained in the Shox2 topologically associated domain (TAD), we demonstrate that this complex cis -regulatory platform is essential for embryonic survival and required for control of region-specific Shox2 expression in multiple developing tissues. While transcription of Shox2 in the embryonic limb is only moderately affected by gene desert loss, Shox2 expression in craniofacial and cardiac domains is nearly abolished. In particular, Shox2 transcripts in the sinus venosus (SV) encompassing the sinoatrial node (SAN) were depleted in embryos lacking the gene desert, likely accounting for the embryonic lethality due to Shox2 -dependency of the SAN pacemaker. Finally, we discover a 1.5kb SV enhancer within the deleted gene desert region, which may act as a genomic module controlling the development of the cardiac conduction system. In summary, our results identify a gene desert indispensable for pleiotropic patterning and highlight the importance of these extensive regulatory landscapes for embryonic development and viability.

Abassah-Oppong, Samuel↗

Predicting metabolic modules in incomplete bacterial genomes with MetaPathPredict

The reconstruction of complete microbial metabolic pathways using ‘omics data from environmental samples remains challenging. Computational pipelines for pathway reconstruction that utilize machine learning methods to predict the presence or absence of KEGG modules in incomplete genomes are lacking. Here, we present MetaPathPredict, a software tool that incorporates machine learning models to predict the presence of complete KEGG modules within bacterial genomic datasets. Using gene annotation data and information from the KEGG module database, MetaPathPredict employs deep learning models to predict the presence of KEGG modules in a genome. MetaPathPredict can be used as a command line tool or as a Python module, and both options are designed to be run locally or on a compute cluster. Benchmarks show that MetaPathPredict makes robust predictions of KEGG module presence within highly incomplete genomes.

59 BASIC BIOLOGICAL SCIENCES↗

Deletion analysis of the itaconic acid biosynthesis gene cluster components in Aspergillus pseudoterreus ATCC32359

The filamentous fungus Aspergillus terreus has been successfully used for industrial production of itaconic acid (IA) for many years. The IA biosynthesis pathway has recently been characterized at a molecular genetic level as an IA gene cluster (IA cluster) by a clone-based transcriptomic approach. The cluster consists of four genes, including genes for cis-aconitic acid decarboxylase (cadA), a predicted transcription factor (tf), mitochondrial organic acid transporter (mttA), and MFS (Major Facilitator Superfamily) type transporter (mfsA). In this research, we performed EST (Expressed Sequence Tag) analysis and systematic gene deletions to investigate the role of those genes in IA biosynthesis in A. pseudoterreus ATCC32359. EST analysis showed a similar expression pattern among those four genes that was distinct from neighbouring genes and further confirmed that they belong to the same biosynthesis cluster. Systematic gene deletion analysis demonstrated that tf, cadA, mttA, and mfsA genes in the cluster are essential for IA production. Deletion of any of them will either completely abolish the IA production or dramatically decrease amount of IA produced. More importantly, a significant amount of aconitic acid was detected in the cadA deletion strain but not in the other deletion strains. Therefore, the cadA deletion strain is a novel microorganism host for production of aconitic acid and other value-added chemicals from lignocellulosic biomass.

Deng, Shuang↗

Synthetic genetic circuits as a means of reprogramming plant roots

We report the shape of a plant’s root system influences its ability to reach essential nutrients in the soil and to acquire water during drought. Progress in engineering plant roots to optimize water and nutrient acquisition has been limited by our capacity to design and build genetic programs that alter root growth in a predictable manner. We developed a collection of synthetic transcriptional regulators for plants that can be compiled to create genetic circuits. These circuits control gene expression by performing Boolean logic operations and can be used to predictably alter root structure. This work demonstrates the potential of synthetic genetic circuits to control gene expression across tissues and reprogram plant growth.

59 BASIC BIOLOGICAL SCIENCES↗

Modification and analysis of context-specific genome-scale metabolic models: methane-utilizing microbial chassis as a case study

ABSTRACT Context-specific genome-scale model (CS-GSM) reconstruction is becoming an efficient strategy for integrating and cross-comparing experimental multi-scale data to explore the relationship between cellular genotypes, facilitating fundamental or applied research discoveries. However, the application of CS modeling for non-conventional microbes is still challenging. Here, we present a graphical user interface that integrates COBRApy, EscherPy, and RIPTiDe, Python-based tools within the BioUML platform, and streamlines the reconstruction and interrogation of the CS genome-scale metabolic frameworks via Jupyter Notebook. The approach was tested using -omics data collected for Methylotuvimicrobium alcaliphilum 20Z R , a prominent microbial chassis for methane capturing and valorization. We optimized the previously reconstructed whole genome-scale metabolic network by adjusting the flux distribution using gene expression data. The outputs of the automatically reconstructed CS metabolic network were comparable to manually optimized i IA409 models for Ca-growth conditions. However, the CS model questions the reversibility of the phosphoketolase pathway and suggests higher flux via primary oxidation pathways. The model also highlighted unresolved carbon partitioning between assimilatory and catabolic pathways at the formaldehyde-formate node. Only a very few genes and only one enzyme with a predicted function in C1 metabolism, a homolog of the formaldehyde oxidation enzyme ( fae1-2 ), showed a significant change in expression in La-growth conditions. The CS-GSM predictions agreed with the experimental measurements under the assumption that the Fae1-2 is a part of the tetrahydrofolate-linked pathway. The cellular roles of the tungsten (W)-dependent formate dehydrogenase ( fdhAB ) and fae homologs ( fae1-2 and fae3 ) were investigated via mutagenesis. The phenotype of the f dhAB mutant followed the model prediction. Furthermore, a more significant reduction of the biomass yield was observed during growth in La-supplemented media, confirming a higher flux through formate. M. alcaliphilum 20Z R mutants lacking fae1-2 did not display any significant defects in methane or methanol-dependent growth. However, contrary to fae1, the fae1-2 homolog failed to restore the formaldehyde-activating enzyme function in complementation tests. Overall, the presented data suggest that the developed computational workflow supports the reconstruction and validation of CS-GSM networks of non-model microbes. IMPORTANCE The interrogation of various types of data is a routine strategy to explore the relationship between genotype and phenotype. An efficient approach for integrating and cross-comparing experimental multi-scale data in the context of whole-genome-based metabolic network reconstruction becomes a powerful tool that facilitates fundamental and applied research discoveries. The present study describes the reconstruction of a context-specific (CS) model for the methane-utilizing bacterium, Methylotuvimicrobium alcaliphilum 20Z R . M. alcaliphilum 20Z R is becoming an attractive microbial platform for the production of biofuels, chemicals, pharmaceuticals, and bio-sorbents for capturing atmospheric methane. We demonstrate that this pipeline can help reconstruct metabolic models that are similar to manually curated networks. Furthermore, the model is able to highlight previously overlooked pathways, thus advancing fundamental knowledge of non-model microbial systems or promoting their development toward biotechnological or environmental implementations.

Kulyashov, M. A.↗

Tunturi virus isolates and metagenome-assembled viral genomes provide insights into the virome of Acidobacteriota in Arctic tundra soils

Arctic soils are climate-critical areas, where microorganisms play crucial roles in nutrient cycling processes. Acidobacteriota are phylogenetically and physiologically diverse bacteria that are abundant and active in Arctic tundra soils. Still, surprisingly little is known about acidobacterial viruses in general and those residing in the Arctic in particular. Here, we applied both culture-dependent and -independent methods to study the virome of Acidobacteriota in Arctic soils. Five virus isolates, Tunturi 1–5, were obtained from Arctic tundra soils, Kilpisjärvi, Finland (69°N), using Tunturiibacter spp. strains originating from the same area as hosts. The new virus isolates have tailed particles with podo- (Tunturi 1, 2, 3), sipho- (Tunturi 4), or myovirus-like (Tunturi 5) morphologies. The dsDNA genomes of the viral isolates are 63–98 kbp long, except Tunturi 5, which is a jumbo phage with a 309-kbp genome. Tunturi 1 and Tunturi 2 share 88% overall nucleotide identity, while the other three are not related to one another. For over half of the open reading frames in Tunturi genomes, no functions could be predicted. To further assess the Acidobacteriota-associated viral diversity in Kilpisjärvi soils, bulk metagenomes from the same soils were explored and a total of 1881 viral operational taxonomic units (vOTUs) were bioinformatically predicted. Almost all vOTUs (98%) were assigned to the class Caudoviricetes. For 125 vOTUs, including five (near-)complete ones, Acidobacteriota hosts were predicted. Acidobacteriota-linked vOTUs were abundant across sites, especially in fens. Terriglobia-associated proviruses were observed in Kilpisjärvi soils, being related to proviruses from distant soils and other biomes. Approximately genus- or higher-level similarities were found between the Tunturi viruses, Kilpisjärvi vOTUs, and other soil vOTUs, suggesting some shared groups of Acidobacteriota viruses across soils. This study provides acidobacterial virus isolates as laboratory models for future research and adds insights into the diversity of viral communities associated with Acidobacteriota in tundra soils. Predicted virus-host links and viral gene functions suggest various interactions between viruses and their host microorganisms. Largely unknown sequences in the isolates and metagenome-assembled viral genomes highlight a need for more extensive sampling of Arctic soils to better understand viral functions and contributions to ecosystem-wide cycling processes in the Arctic.

54 ENVIRONMENTAL SCIENCES↗

Targeted curation of the gut microbial gene content modulating human cardiovascular disease

Despite the promise of the gut microbiome to predict human health, few studies expose the molecular-scale processes underpinning such forecasts. We mined over 200,000 gut-derived genomes from cultivated and uncultivated microbial lineages to inventory the gut microorganisms and their gene content that control trimethylamine-induced cardiovascular disease. We assigned an atherosclerotic profile to the 6,341 microbial genomes that encoded metabolisms associated with heart disease, creating the Methylated Amine Gene Inventory of Catabolism database (MAGICdb). From microbiome gene expression data sets, we demonstrate that MAGICdb enhanced the recovery of disease-relevant genes and identified the most active microorganisms, unveiling future therapeutic targets. From the feces of healthy and diseased subjects, we show that MAGICdb predicted cardiovascular disease status as effectively as traditional lipid blood tests. This functional microbiome catalog is a public, exploitable resource, designed to enable a new era of microbiota-based therapeutics and diagnostics

metatranscriptomics↗

Peatland Microbial Community Composition Is Driven by a Natural Climate Gradient

Peatlands are important players in climate change–biosphere feedbacks via long-term net carbon (C) accumulation in soil organic matter and as potential net C sources including the potent greenhouse gas methane (CH 4 ). Interactions of climate, site-hydrology, plant community, and groundwater chemical factors influence peatland development and functioning, including C dioxide (CO 2 ) and CH 4 fluxes, but the role of microbial community composition is not well understood. To assess microbial functional and taxonomic dissimilarities, we used high throughput sequencing of the small subunit ribosomal DNA (SSU rDNA) to determine bacterial and archaeal community composition in soils from twenty North American peatlands. Targeted DNA metabarcoding showed that although Proteobacteria, Acidobacteria, and Actinobacteria were the dominant phyla on average, intermediate and rich fens hosted greater diversity and taxonomic richness, as well as an array of candidate phyla when compared with acidic and nutrient-poor poor fens and bogs. Moreover, pH was revealed to be the strongest predictor of microbial community structure across sites. Predictive metagenome content (PICRUSt) showed increases in specific genes, such as purine/pyrimidine and amino-acid metabolism in mid-latitude peatlands from 38 to 45° N, suggesting a shift toward utilization of microbial biomass over utilization of initial plant biomass in these microbial communities. Overall, there appears to be noticeable differences in community structure between peatland classes, as well as differences in microbial metabolic activity between latitudes. Finally, these findings are in line with a predicted increase in the decomposition and accelerated C turnover, and suggest that peatlands north of 37° latitude may be particularly vulnerable to climate change.

59 BASIC BIOLOGICAL SCIENCES↗

In silico λ-dynamics predicts protein binding specificities to modified RNAs

Abstract RNA modifications shape gene expression through a variety of chemical changes to canonical RNA bases. Although numbering in the hundreds, only a few RNA modifications are well characterized, in part due to the absence of methods to identify modification sites. Antibodies remain a common tool to identify modified RNA and infer modification sites through straightforward applications. However, specificity issues can result in off-target binding and confound conclusions. This work utilizes in silico λ-dynamics to efficiently estimate binding free energy differences of modification-targeting antibodies between a variety of naturally occurring RNA modifications. Crystal structures of inosine and N6-methyladenosine (m6A) targeting antibodies bound to their modified ribonucleosides were determined and served as structural starting points. λ-Dynamics was utilized to predict RNA modifications that permit or inhibit binding to these antibodies. In vitro RNA-antibody binding assays supported the accuracy of these in silico results. High agreement between experimental and computed binding propensities demonstrated that λ-dynamics can serve as a predictive screen for antibody specificity against libraries of RNA modifications. More importantly, this strategy is an innovative way to elucidate how hundreds of known RNA modifications interact with biological molecules without the limitations imposed by in vitro or in vivo methodologies.

Biochemistry & Molecular Biology↗

Functional characterization of glycosyltransferases in duckweed to enable predictive biology

Glycosyltransferases (GTs) catalyze the formation of glycosidic linkages to produce almost all complex carbohydrates. This project used a multi-disciplinary, high-throughput (HTP) biochemical and computational biology approach focused on duckweed as a model energy crop, to study carbohydrate metabolic processes. To achieve this, developed and carried out out high-throughput (HTP) functional characterization of plant glycosyltransferases (GTs) role of enzymatic microenvironments be assessed through a combined proteomic and computational biology approach, and the combined data was used to populate deep-learning frameworks to predict plant GT function. Functional validation achieved through this research is being used to assign gene function and study plant processes at the systems level to efficiently link the genome sequence with gene function. Together, the combined approaches used within this study provide a foundation for how computational prediction, in combination with high-throughput functional validation, can be used to study plant processes at the systems level and translate knowledge gained to efficiently link genome sequence with gene function in a species agnostic manner.

09 BIOMASS FUELS↗

YEASTRACT+: a portal for the exploitation of global transcription regulation and metabolic model data in yeast biotechnology and pathogenesis

Abstract YEASTRACT+ (http://yeastract-plus.org/) is a tool for the analysis, prediction and modelling of transcription regulatory data at the gene and genomic levels in yeasts. It incorporates three integrated databases: YEASTRACT (http://yeastract-plus.org/yeastract/), PathoYeastract (http://yeastract-plus.org/pathoyeastract/) and NCYeastract (http://yeastract-plus.org/ncyeastract/), focused on Saccharomyces cerevisiae, pathogenic yeasts of the Candida genus, and non-conventional yeasts of biotechnological relevance. In this release, YEASTRACT+ offers upgraded information on transcription regulation for the ten previously incorporated yeast species, while extending the database to another pathogenic yeast, Candida auris. Since the last release of YEASTRACT+ (January 2020), a fourth database has been integrated. CommunityYeastract (http://yeastract-plus.org/community/) offers a platform for the creation, use, and future update of YEASTRACT-like databases for any yeast of the users’ choice. CommunityYeastract currently provides information for two Saccharomyces boulardii strains, Rhodotorula toruloides NP11 oleaginous yeast, and Schizosaccharomyces pombe 972h-. In addition, YEASTRACT+ portal currently gathers 304 547 documented regulatory associations between transcription factors (TF) and target genes and 480 DNA binding sites, considering 2771 TFs from 11 yeast species. A new set of tools, currently implemented for S. cerevisiae and C. albicans, is further offered, combining regulatory information with genome-scale metabolic models to provide predictions on the most promising transcription factors to be exploited in cell factory optimisation or to be used as novel drug targets. The expansion of these new tools to the remaining YEASTRACT+ species is ongoing.

Teixeira, Miguel Cacho (ORCID:0000000256766174)↗

Peripheral Blood Transcript Signatures after Internal 131 I-mIBG Therapy in Relapsed and Refractory Neuroblastoma Patients Identifies Early and Late Biomarkers of Internal 131 I Exposures

131 I-metaiodobenzylguanidine ( 131 I-mIBG) is a targeted radiation therapy developed for the treatment of advanced neuroblastoma. We have previously shown that this patient cohort can be used to predict absorbed dose associated with early 131 I exposure, 72 h after treatment. We now expand these studies to identify gene expression differences associated with 131 I-mIBG exposure 15 days after treatment. Total RNA from peripheral blood lymphocytes was isolated from 288 whole blood samples representing 59 relapsed or refractory neuroblastoma patients before and after 131 I-mIBG treatment. We found that several transcripts predictive of early exposure returned to baseline levels by day 15, however, selected transcripts did not return to baseline. At 72 h, all 17 selected pathway-specific transcripts were differentially expressed. Transcripts CDKN1A (P < 0.000001), FDXR (P < 0.000001), DDB2 (P < 0.000001), and BBC3 (P < 0.000001) showed the highest up-regulation at 72 h after 131 I-mIBG exposure, with mean log2 fold changes of 2.55, 2.93, 1.86 and 1.85, respectively. At day 15 after 131I-mIBG, 11 of the 17 selected transcripts were differentially expressed, with XPC, STAT5B, PRKDC, MDM2, POLH, IGF1R, and SGK1 displaying significant up-regulation at 72 h and significant down-regulation at day 15. Interestingly, transcripts FDXR (P = 0.01), DDB2 (P = 0.03), BCL2 (P = 0.003), and SESN1 (P < 0.0003) maintained differential expression 15 days after 131 I-mIBG treatment. These results suggest that transcript levels for DNA repair, apoptosis, and ionizing radiation-induced cellular stress are still changing by 15 days after 131 I-mIBG treatment. Our studies showcase the use of biodosimetry gene expression panels as predictive biomarkers following early (72 h) and late (15 days) internal 131 I exposure. Our findings also demonstrate the utility of our transcript panel to differentiate exposed from non-exposed individuals up to 15 days after exposure from internal 131 I.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

The Architecture of Metabolism Maximizes Biosynthetic Diversity in the Largest Class of Fungi

Abstract Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or “specialized” metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered the most phylogenetically diverse class of fungi and known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation among sets of BGCs in individual genomes is due to nonoverlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

59 BASIC BIOLOGICAL SCIENCES↗

Novosphingobium aromaticivorans LigR coordinates transcription of genes involved in metabolism of multiple types of aromatics

Aromatic compounds are a ubiquitous and diverse family of chemicals with functions as biomolecules, natural products, industrial chemicals, and pollutants. Novosphingobium aromaticivorans DSM 12444 uses multiple inducible pathways to catabolize H-, G-, and S-type aromatics that contain zero, one, or two methoxy groups, respectively. Here, we obtain a systems-level view of the transcriptional control of its aromatic metabolic pathways. Several in vitro analyses found that a N. aromaticivorans homolog of the Sphingobium lignivorans SYK-6 transcription factor LigR bound genomic DNA upstream of genes involved in metabolism of multiple aromatic types. We found that a ΔLigR mutant had growth defects on all three types of aromatics as sole carbon sources. Transcriptomic analysis revealed that LigR was required to increase expression of gene products that function in metabolism of all three aromatic types. We also found that, in media containing both glucose and an aromatic carbon source, the ΔLigR mutant directed intermediates through alternative aromatic metabolic pathways. Protein-DNA binding assays showed that N. aromaticivorans LigR binds immediately upstream of promoters of genes involved in aromatic metabolism. We found that N. aromaticivorans LigR coordinates the expression of enzymes that function in the catabolism of H-, G-, and S-type aromatics, and that there are differences in the role of LigR in N. aromaticivorans and S. lignivorans. A comparative genomic analysis predicted that LigR homologs and the aromatic-metabolizing genes that it directly regulates are often co-localized in the genomes of Sphingomonadales, but often not found in this arrangement in many other known aromatic metabolizing bacteria.

Aromatic Compound Degradation↗

Evaluating E. coli genome‐scale metabolic model accuracy with high‐throughput mutant fitness data

Abstract The Escherichia coli genome‐scale metabolic model (GEM) is an exemplar systems biology model for the simulation of cellular metabolism. Experimental validation of model predictions is essential to pinpoint uncertainty and ensure continued development of accurate models. Here, we quantified the accuracy of four subsequent E. coli GEMs using published mutant fitness data across thousands of genes and 25 different carbon sources. This evaluation demonstrated the utility of the area under a precision–recall curve relative to alternative accuracy metrics. An analysis of errors in the latest (iML1515) model identified several vitamins/cofactors that are likely available to mutants despite being absent from the experimental growth medium and highlighted isoenzyme gene‐protein‐reaction mapping as a key source of inaccurate predictions. A machine learning approach further identified metabolic fluxes through hydrogen ion exchange and specific central metabolism branch points as important determinants of model accuracy. This work outlines improved practices for the assessment of GEM accuracy with high‐throughput mutant fitness data and highlights promising areas for future model refinement in E. coli and beyond.

59 BASIC BIOLOGICAL SCIENCES↗

Subnetwork representation learning for discovering network biomarkers in predicting lymph node metastasis in early oral cancer

Cervical lymph node metastasis is the leading cause of poor prognosis in oral tongue squamous cell carcinoma and also occurs in the early stages. The current clinical diagnosis depends on a physical examination that is not enough to determine whether micrometastasis remains. The transcriptome profiling technique has shown great potential for predicting micrometastasis by capturing the dynamic activation state of genes. However, there are several technical challenges in using transcriptome data to model patient conditions: (1) An Insufficient number of samples compared to the number of genes, (2) Complex dependence between genes that govern the cancer phenotype, and (3) Heterogeneity between patients between cohorts that differ geographically and racially. We developed a computational framework to learn the subnetwork representation of the transcriptome to discover network biomarkers and determine the potential of metastasis in early oral tongue squamous cell carcinoma. Our method achieved high accuracy in predicting the potential of metastasis in two geographically and racially different groups of patients. The robustness of the model and the reproducibility of the discovered network biomarkers show great potential as a tool to diagnose lymph node metastasis in early oral cancer.

59 BASIC BIOLOGICAL SCIENCES↗