Search NASA⌕ Search

SEARCH · Search NASA

Results for “comparative transcriptomics”

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 235 records · Page 13

Multiple origins, one evolutionary trajectory: gradual evolution characterizes distinct lineages of allotetraploid Brachypodium

The “genomic shock” hypothesis posits that unusual challenges to genome integrity such as whole genome duplication may induce chaotic genome restructuring. Decades of research on polyploid genomes have revealed that this is often, but not always the case. While some polyploids show major chromosomal rearrangements and derepression of transposable elements in the immediate aftermath of whole genome duplication, others do not. Nonetheless, all polyploids show gradual diploidization over evolutionary time. To evaluate these hypotheses, we produced a chromosome-scale reference genome for the natural allotetraploid grass Brachypodium hybridum, accession “Bhyb26.” We compared 2 independently derived accessions of B. hybridum and their deeply diverged diploid progenitor species Brachypodium stacei and Brachypodium distachyon. The 2 B. hybridum lineages provide a natural timecourse in genome evolution because one formed 1.4 million years ago, and the other formed 140 thousand years ago. The genome of the older lineage reveals signs of gradual post-whole genome duplication genome evolution including minor gene loss and genome rearrangement that are missing from the younger lineage. In neither B. hybridum lineage do we find signs of homeologous recombination or pronounced transposable element activation, though we find evidence supporting steady post-whole genome duplication transposable element activity in the older lineage. Gene loss in the older lineage was slightly biased toward 1 subgenome, but genome dominance was not observed at the transcriptomic level. We propose that relaxed selection, rather than an abrupt genomic shock, drives evolutionary novelty in B. hybridum, and that the progenitor species’ similarity in transposable element load may account for the subtlety of the observed genome dominance.

59 BASIC BIOLOGICAL SCIENCES↗

Multilevel analysis between Physcomitrium patens and Mortierellaceae endophytes explores potential long‐standing interaction among land plants and fungi

SUMMARY The model moss species Physcomitrium patens has long been used for studying divergence of land plants spanning from bryophytes to angiosperms. In addition to its phylogenetic relationships, the limited number of differential tissues, and comparable morphology to the earliest embryophytes provide a system to represent basic plant architecture. Based on plant–fungal interactions today, it is hypothesized these kingdoms have a long‐standing relationship, predating plant terrestrialization. Mortierellaceae have origins diverging from other land fungi paralleling bryophyte divergence, are related to arbuscular mycorrhizal fungi but are free‐living, observed to interact with plants, and can be found in moss microbiomes globally. Due to their parallel origins, we assess here how two Mortierellaceae species, Linnemannia elongata and Benniella erionia , interact with P. patens in coculture. We also assess how Mollicute ‐related or Burkholderia ‐related endobacterial symbionts (MRE or BRE) of these fungi impact plant response. Coculture interactions are investigated through high‐throughput phenomics, microscopy, RNA‐sequencing, differential expression profiling, gene ontology enrichment, and comparisons among 99 other P. patens transcriptomic studies. Here we present new high‐throughput approaches for measuring P. patens growth, identify novel expression of over 800 genes that are not expressed on traditional agar media, identify subtle interactions between P. patens and Mortierellaceae, and observe changes to plant–fungal interactions dependent on whether MRE or BRE are present. Our study provides insights into how plants and fungal partners may have interacted based on their communications observed today as well as identifying L. elongata and B. erionia as modern fungal endophytes with P. patens .

59 BASIC BIOLOGICAL SCIENCES↗

Characterization of a Unique Pair of Ferredoxin and Ferredoxin NADP+ Reductase Isoforms That Operates in Non-Photosynthetic Glandular Trichomes

Our recent investigations indicated that isoforms of ferredoxin (Fd) and ferredoxin NADP+ reductase (FNR) play essential roles for the reductive steps of the 2C-methyl-D-erythritol 4-phosphate (MEP) pathway of terpenoid biosynthesis in peppermint glandular trichomes (GTs). Based on an analysis of several transcriptome data sets, we demonstrated the presence of transcripts for a leaf-type FNR (L-FNR), a leaf-type Fd (Fd I), a root-type FNR (R-FNR), and two root-type Fds (Fd II and Fd III) in several members of the mint family (Lamiaceae). The present study reports on the biochemical characterization of all Fd and FNR isoforms of peppermint (Mentha × piperita L.). The redox potentials of Fd and FNR isoforms were determined using photoreduction methods. Based on a diaphorase assay, peppermint R-FNR had a substantially higher specificity constant (kcat/Km) for NADPH than L-FNR. Similar results were obtained with ferricyanide as an electron acceptor. When assayed for NADPH–cytochrome c reductase activity, the specificity constant with the Fd II and Fd III isoforms (when compared to Fd I) was slightly higher for L-FNR and substantially higher for R-FNR. Based on real-time quantitative PCR assays with samples representing various peppermint organs and cell types, the Fd II gene was expressed very highly in metabolically active GTs (but also present at lower levels in roots), whereas Fd III was expressed at low levels in both roots and GTs. Our data provide evidence that high transcript levels of Fd II, and not differences in the biochemical properties of the encoded enzyme when compared to those of Fd III, are likely to support the formation of copious amounts of monoterpene via the MEP pathway in peppermint GTs. This work has laid the foundation for follow-up studies to further investigate the roles of a unique R-FNR–Fd II pair in non-photosynthetic GTs of the Lamiaceae.

59 BASIC BIOLOGICAL SCIENCES↗

Differential transcriptomic alterations in nasal versus lung tissue of acrolein-exposed rats

Introduction: Acrolein is a significant component of anthropogenic and wildfire emissions, as well as cigarette smoke. Although acrolein primarily deposits in the upper respiratory tract upon inhalation, patterns of site-specific injury in nasal versus pulmonary tissues are not well characterized. This assessment is critical in the design of in vitro and in vivo studies performed for assessing health risk of irritant air pollutants. Methods: In this study, male and female Wistar-Kyoto rats were exposed nose-only to air or acrolein. Rats in the acrolein exposure group were exposed to incremental concentrations of acrolein (0, 0.1, 0.316, 1 ppm) for the first 30 min, followed by a 3.5 h exposure at 3.16 ppm. In the first cohort of male and female rats, nasal and bronchoalveolar lavage fluids were analyzed for markers of inflammation, and in a second cohort of males, nasal airway and left lung tissues were used for mRNA sequencing. Results: Protein leakage in nasal airways of acrolein-exposed rats was similar in both sexes; however, inflammatory cells and cytokine increases were more pronounced in males when compared to females. No consistent changes were noted in bronchoalveolar lavage fluid of males or females except for increases in total cells and IL-6. Acrolein-exposed male rats had 452 differentially expressed genes (DEGs) in nasal tissue versus only 95 in the lung. Pathway analysis of DEGs in the nose indicated acute phase response signaling, Nrf2-mediated oxidative stress, unfolded protein response, and other inflammatory pathways, whereas in the lung, xenobiotic metabolism pathways were changed. Genes associated with glucocorticoid and GPCR signaling were also changed in the nose but not in the lung. Discussion: These data provide insights into inhaled acrolein-mediated sex-specific injury/inflammation in the nasal and pulmonary airways. The transcriptional response in the nose reflects acrolein-induced acute oxidative and cytokine signaling changes, which might have implications for upper airway inflammatory disease susceptibility.

Alewel, Devin I.↗

Single-Cell RNA-Seq Reveals Transcriptomic Heterogeneity and Post-Traumatic Osteoarthritis-Associated Early Molecular Changes in Mouse Articular Chondrocytes

Articular cartilage is a connective tissue lining the surfaces of synovial joints. When the cartilage severely wears down, it leads to osteoarthritis (OA), a debilitating disease that affects millions of people globally. The articular cartilage is composed of a dense extracellular matrix (ECM) with a sparse distribution of chondrocytes with varying morphology and potentially different functions. Elucidating the molecular and functional profiles of various chondrocyte subtypes and understanding the interplay between these chondrocyte subtypes and other cell types in the joint will greatly expand our understanding of joint biology and OA pathology. Although recent advances in high-throughput OMICS technologies have enabled molecular-level characterization of tissues and organs at an unprecedented resolution, thorough molecular profiling of articular chondrocytes has not yet been undertaken, which may be in part due to the technical difficulties in isolating chondrocytes from dense cartilage ECM. In this study, we profiled articular cartilage from healthy and injured mouse knee joints at a single-cell resolution and identified nine chondrocyte subtypes with distinct molecular profiles and injury-induced early molecular changes in these chondrocytes. We also compared mouse chondrocyte subpopulations to human chondrocytes and evaluated the extent of molecular similarity between mice and humans. This work expands our view of chondrocyte heterogeneity and rapid molecular changes in chondrocyte populations in response to joint trauma and highlights potential mechanisms that trigger cartilage degeneration.

59 BASIC BIOLOGICAL SCIENCES↗

Machine learning analysis of RB-TnSeq fitness data predicts functional gene modules in Pseudomonas putida KT2440

ABSTRACT There is growing interest in engineering Pseudomonas putida KT2440 as a microbial chassis for the conversion of renewable and waste-based feedstocks, and metabolic engineering of P. putida relies on the understanding of the functional relationships between genes. In this work, independent component analysis (ICA) was applied to a compendium of existing fitness data from randomly barcoded transposon insertion sequencing (RB-TnSeq) of P. putida KT2440 grown in 179 unique experimental conditions. ICA identified 84 independent groups of genes, which we call fModules (“functional modules”), where gene members displayed shared functional influence in a specific cellular process. This machine learning-based approach both successfully recapitulated previously characterized functional relationships and established hitherto unknown associations between genes. Selected gene members from fModules for hydroxycinnamate metabolism and stress resistance, acetyl coenzyme A assimilation, and nitrogen metabolism were validated with engineered mutants of P. putida . Additionally, functional gene clusters from ICA of RB-TnSeq data sets were compared with regulatory gene clusters from prior ICA of RNAseq data sets to draw connections between gene regulation and function. Because ICA profiles the functional role of several distinct gene networks simultaneously, it can reduce the time required to annotate gene function relative to manual curation of RB-TnSeq data sets. IMPORTANCE This study demonstrates a rapid, automated approach for elucidating functional modules within complex genetic networks. While Pseudomonas putida randomly barcoded transposon insertion sequencing data were used as a proof of concept, this approach is applicable to any organism with existing functional genomics data sets and may serve as a useful tool for many valuable applications, such as guiding metabolic engineering efforts in other microbes or understanding functional relationships between virulence-associated genes in pathogenic microbes. Furthermore, this work demonstrates that comparison of data obtained from independent component analysis of transcriptomics and gene fitness datasets can elucidate regulatory-functional relationships between genes, which may have utility in a variety of applications, such as metabolic modeling, strain engineering, or identification of antimicrobial drug targets.

09 BIOMASS FUELS↗

Impact of Porcine Arterivirus, Influenza B, and Their Coinfection on Antiviral Response in the Porcine Lung

Interferon (IFN) cytokines induce an autonomous antiviral state in cells of the infected site to restrict virus spreading and critically regulate overall antiviral response. The antiviral state leads to host protection through expression of hundreds of IFN-stimulated genes that restrict viral infection through multiple mechanisms, for example, directly in viral genome degradation and indirectly through cellular metabolic inhibition. Young pigs were split into four treatment groups: control, porcine reproductive and respiratory syndrome virus (PRRSV, also known as porcine arterivirus) infected, influenza B virus (IBV) infected, and IBV/PRRSV coinfection. Lung tissue was collected at 3, 5, and 7 days post infection (dpi) for control, PRRSV and IBV/PRRSV coinfection, and at 3 and 5 dpi for IBV. Transcriptomic analysis, using usegalaxy.org tools, was performed against the S.scrofa 11.1 reference genome. Differentially expressed gene (DEG) analysis was carried out using DeSeq2 based on the model treatment + dpi + treatment:dpi + E. Downstream analysis examined the interaction of DEG at each dpi for over-enriched gene ontology (G.O.) terms and pathways. Comparisons of the infected groups vs. the controls yielded a total of (n = 1412) DEGs for the PRRSV group and (n = 1578) for the IBV/PRRSV group across all timepoints. The IBV group had (n = 64) total DEGs across 3 and 5 dpi. Expression data were considered statistically significant based on false discovery rate (FDR) ⫹ 0.1. Venn diagram comparisons of the DEGs across dpi showed that groups shared only 16 DEGs at 3 dpi, no DEGs were shared at 5 dpi, and for 7 dpi, only the PRRSV and IBV/PRRSV groups were compared and shared a total of 43 DEGs. Across the comparisons, differential expression was observed in antiviral genes such as IRF1, MX1, and OAS2. The IBV and IBV/PRRSV groups showed higher expression of antiviral genes at earlier dpi than the PRRSV group. Additionally, downregulated genes from the comparisons clustered around Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways effecting lung development and cellular integrity. Early expression of host IFN and antiviral genes may lead to viral RNA degradation, and assembly and transcription inhibition in the IBV infections. In comparison, expression of antiviral genes in the PRRSV group decreased across time. The decrease may explain why PRRSV infections persist, while IBV clears. Moreover, all infected groups showed prolonged upregulation in neutrophil degranulation pathway activity, possibly exacerbating symptomatic lung lesion pathology seen in these respiratory infections.

59 BASIC BIOLOGICAL SCIENCES↗

Brachypodium Phenylalanine Ammonia Lyase (PAL) Promotes Antiviral Defenses against Panicum mosaic virus and Its Satellites

Brachypodium distachyon has recently emerged as a premier model plant for monocot biology, akin to Arabidopsis thaliana. We previously reported genome-wide transcriptomic and alternative splicing changes occurring in Brachypodium during compatible infections with Panicum mosaic virus (PMV) and its satellite virus (SPMV). Here, we dissected the role of Brachypodium phenylalanine ammonia lyase 1 (PAL1), a key enzyme for phenylpropanoid and salicylic acid (SA) biosynthesis and the induction of plant defenses. Targeted metabolomics profiling of PMV-infected and PMV- plus SPMV-infected (PMV/SPMV) Brachypodium plants revealed enhanced levels of multiple defense-related hormones and metabolites such as cinnamic acid, SA, and fatty acids and lignin precursors during disease progression. The virus-induced accumulation of SA and lignin was significantly suppressed upon knockdown of B. distachyon PAL1 (BdPAL1) using RNA interference (RNAi). The compromised SA accumulation in PMV/SPMV-infected BdPAL1 RNAi plants correlated with weaker induction of multiple SA-related defense gene markers (pathogenesis related 1 [PR-1], PR-3, PR-5, and WRKY75) and enhanced susceptibility to PMV/SPMV compared to that of wild-type (WT) plants. Furthermore, exogenous application of SA alleviated the PMV/SPMV necrotic disease phenotypes and delayed plant death caused by single and mixed infections. Together, our results support an antiviral role for BdPAL1 during compatible host-virus interaction, perhaps as a last resort attempt to rescue the infected plant.

59 BASIC BIOLOGICAL SCIENCES↗

A k-mer based approach for classifying viruses without taxonomy identifies viral associations in human autism and plant microbiomes

Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an essential role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus’ genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples, even if they do not have a taxonomy. To identify and quantify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses to create a database of metagenomic viruses. The resulting Kraken2 database of the metagenomic viruses can be found here: https://www.osti.gov/biblio/1615774 and is compatible with Kraken2. We then integrated the viral classification database with databases created with genomes from NCBI for use with ParaKraken (a parallelized version of Kraken provided in Supplemental Zip 1), a metagenomic/transcriptomic classifier. To illustrate the breadth of our utility for classifying metagenome viruses, we analyzed data from a plant metagenome study identifying genotypic and compartment specific differences between two Populus genotypes in three different compartments. We also identified a significant increase in abundance of eight viral sequences in post mortem brains in a human metatranscriptome study comparing Autism Spectrum Disorder patients and controls. We also show the potential accuracy for classifying viruses by utilizing both the JGI and NCBI viral databases to identify the uniqueness of viral sequences. Finally, we validate the accuracy of viral classification with NCBI databases containing viruses with taxonomy to identify pathogenic viruses in known COVID-19 and cassava brown streak virus infection samples. Our method represents the compulsory first step in better understanding the role of viruses in the microbiome by allowing for a more complete identification of sequences without taxonomy. Better classification of viruses will improve identifying associations between viruses and their hosts as well as viruses and other microbiome members. Despite the lack of taxonomy, this database of metagenomic viruses can be used with any tool that utilizes a taxonomy, such as Kraken, for accurate classification of viruses.

59 BASIC BIOLOGICAL SCIENCES↗

Transcriptional, metabolic, physiological and developmental responses of switchgrass to phosphorus limitation

Abstract Knowing how switchgrass ( Panicum virgatum L . ) responds and adapts to phosphorus (P)‐limitation will aid efforts to optimize P acquisition and use in this species for sustainable biomass production. This integrative study investigated the impacts of mild, moderate, and severe P‐stress on genome transcription and whole‐plant metabolism, physiology and development in switchgrass. P‐limitation reduced overall plant growth, increased root/shoot ratio, increased root branching at moderate P‐stress, and decreased root diameter with increased density and length of root hairs at severe P‐stress. RNA‐seq analysis revealed thousands of genes that were differentially expressed under moderate and severe P‐stress in roots and/or shoots compared to P‐replete plants, with many stress‐induced genes involved in transcriptional and other forms of regulation, primary and secondary metabolism, transport, and other processes involved in P‐acquisition and homeostasis. Amongst the latter were multiple miRNA399 genes and putative targets of these. Metabolite profiling showed that levels of most sugars and sugar alcohols decreased with increasing P stress, while organic and amino acids increased under mild and moderate P‐stress in shoots and roots, although this trend reversed under severe P‐stress, especially in shoots.

59 BASIC BIOLOGICAL SCIENCES↗

Coupling Environmental Whole Mixture Toxicity Screening with Unbiased RNA-Seq Reveals Site-Specific Biological Responses in Zebrafish

Passive sampling device (PSD) extracts paired with developmental toxicity assays in Danio Rerio (zebrafish) are excellent sensors for whole mixture toxicity associated with the bioavailable non-polar organics at environmental sites. We expand this concept by incorporating RNA-Seq in 48-h post fertilization zebrafish statically exposed to PSD extracts from two Portland Harbor Superfund Site locations: river mile 6.5W (RM 6.5W) and river mile 7W (RM 7W). RM 6.5W contained higher concentrations of polycyclic aromatic hydrocarbons (PAHs), but the diagnostic ratios of both extracts indicated similar PAH sourcing and composition. Developmental screens determined RM 6.5W to be more toxic with the most sensitive endpoint being a “wavy” notochord malformation. Differential gene expression from exposure to both extracts was largely parallel, although more pronounced for RM 6.5W. When compared to the gene expression associated with individual chemical exposures, PSD extracts produced some gene signatures parallel to PAHs but were more closely matched by oxygenated-PAHs. Additionally, differential expression, reminiscent of the wavy notochord phenotype, was not accounted for by either class of chemical, indicating the potential of other contaminants driving mixture toxicity. These techniques offer a compelling method for non-targeted hazard characterization of whole mixtures in an in vivo vertebrate system without requiring complete chemical characterization.

59 BASIC BIOLOGICAL SCIENCES↗

Hemiptera phylogenomic resources: Tree‐based orthology prediction and conserved exon identification

Abstract High‐throughput sequencing of transcriptomes and targeted genomic regions are advancing our knowledge of The Tree of Life. Building phylogenies with regions of the genome requires 1‐to‐1 orthologue resources of genes and noncoding loci. One organismal group that has received little attention in this area is the Hemiptera, the fifth largest insect order represented by ~103,590 named species. Here, we present a set of 3,872 Hemiptera 1‐to‐1 orthogroups based on tree‐based orthology inference of eight Hemiptera species with publicly available genome sequences. We also estimate a set of 406 orthologous exons with similar mRNA splice sites that can be used for Sanger sequencing and develop enrichment probes for targeted genome sequencing for phylogenomic inference. We show this novel set of orthologues is informative at the protein, coding sequence and exon molecular levels and provides robust branch support in both gene tree–species tree methods and concatenated sequence phylogenies. In addition, we demonstrate the utility of these loci to resolve relationships in whiteflies, Bemisia tabaci , a large species complex with few phylogenomic resources. Last, we compare our Hemiptera phylogeny with previously published phylogenies and other orthologue databases, while providing suggestions on further improvement to this phylogenomic resource.

Owen, Christopher L.↗

Protein N -Glycans in Healthy and Sclerotic Glomeruli in Diabetic Kidney Disease

Diabetes is expected to directly affect renal glycosylation; yet to date, there has not been a comprehensive evaluation of alterations in N-glycan composition in the glomeruli of patients with diabetic kidney disease (DKD). Here, we used untargeted mass spectrometry imaging to identify N-glycan structures in healthy and sclerotic glomeruli in formalin-fixed paraffin-embedded sections from needle biopsies of five patients with DKD and three healthy kidney samples. Regional proteomics was performed on glomeruli from additional biopsies from the same patients to compare the abundances of enzymes involved in glycosylation. Secondary analysis of single-nucleus RNA sequencing (snRNAseq) data were used to inform on transcript levels of glycosylation machinery in different cell types and states. We detected 120 N-glycans, and among them, we identified 12 of these protein post-translated modifications that were significantly increased in glomeruli. All glomeruli-specific N-glycans contained an N-acetyllactosamine epitope. Five N-glycan structures were highly discriminant between sclerotic and healthy glomeruli. Sclerotic glomeruli had an additional set of glycans lacking fucose linked to their core, and they did not show tetra-antennary structures that were common in healthy glomeruli. Orthogonal omics analyses revealed lower protein abundance and lower gene expression involved in synthesizing fucosylated and branched N-glycans in sclerotic podocytes. In snRNAseq and regional proteomics analyses, we observed that genes and/or proteins involved in sialylation and N-acetyllactosamine synthesis were also downregulated in DKD glomeruli, but this alteration remained undetectable by our spatial N-glycomics assay. Integrative spatial glycomics, proteomics, and transcriptomics revealed protein N-glycosylation characteristic of sclerotic glomeruli in DKD.

60 APPLIED LIFE SCIENCES↗

Metapangenomics reveals depth-dependent shifts in metabolic potential for the ubiquitous marine bacterial SAR324 lineage

Abstract Background Oceanic microbiomes play a pivotal role in the global carbon cycle and are central to the transformation and recycling of carbon and energy in the ocean’s interior. SAR324 is a ubiquitous but poorly understood uncultivated clade of Deltaproteobacteria that inhabits the entire water column, from ocean surface waters to its deep interior. Although some progress has been made in elucidating potential metabolic traits of SAR324 in the dark ocean, very little is known about the ecology and the metabolic capabilities of this group in the euphotic and twilight zones. To investigate the comparative genomics, ecology, and physiological potential of the SAR324 clade, we examined the distribution and variability of key genomic features and metabolic pathways in this group from surface waters to the abyss in the North Pacific Subtropical Gyre, one of the largest biomes on Earth. Results We leveraged a pangenomic ecological approach, combining spatio-temporally resolved single-amplified genome, metagenomic, and metatranscriptomic datasets. The data revealed substantial genomic diversity throughout the SAR324 clade, with distinct depth and temporal distributions that clearly differentiated ecotypes. Phylogenomic subclade delineation, environmental distributions, genomic feature similarities, and metabolic capacities revealed strong congruence. The four SAR324 ecotypes delineated in this study revealed striking divergence from one another with respect to their habitat-specific metabolic potentials. The ecotypes living in the dark or twilight oceans shared genomic features and metabolic capabilities consistent with a sulfur-based chemolithoautotrophic lifestyle. In contrast, those inhabiting the sunlit ocean displayed higher plasticity energy-related metabolic pathways, supporting a presumptive photoheterotrophic lifestyle. In epipelagic SAR324 ecotypes, we observed the presence of two types of proton-pumping rhodopsins, as well as genomic, transcriptomic, and ecological evidence for active photoheterotrophy, based on xanthorhodopsin-like light-harvesting proteins. Conclusions Combining pangenomic and both metagenomic and metatranscriptomic profiling revealed a striking divergence in the vertical distribution, genomic composition, metabolic potential, and predicted lifestyle strategies of geographically co-located members of the SAR324 bacterial clade. The results highlight the utility of metapangenomic approaches employed across environmental gradients, to decipher the properties and variation in function and ecological traits of specific phylogenetic clades within complex microbiomes.

59 BASIC BIOLOGICAL SCIENCES↗

Oncogenic Ras and ΔNp63α cooperate to recruit immunosuppressive polymorphonuclear myeloid-derived suppressor cells in a mouse model of squamous cancer pathogenesis

Introduction Amplification of human chromosome 3q26-29, which encodes oncoprotein ΔNp63 among other isoforms of the p63 family, is a feature common to squamous cell carcinomas (SCCs) of multiple tissue origins. Along with overexpression of ΔNp63, activation of the protooncogene, RAS , whether by overexpression or oncogenic mutation, is frequently observed in many cancers. In this study, analysis of transcriptome data from The Cancer Genome Atlas (TCGA) demonstrated that expression of TP63 mRNA , particularly ΔNp63 isoforms, and HRAS are significantly elevated in advanced squamous cell carcinomas of the head and neck (HNSCCs), suggesting pathological significance. However, how co-overexpressed ΔNp63 and HRAS affect the immunosuppressive tumor microenvironment (TME) is incompletely understood. Methods Here, we established and characterized an immune competent mouse model using primary keratinocytes with retroviral-mediated overexpression of ΔNp63α and constitutively activated HRAS (v-ras Ha G12R) to evaluate the role of these oncogenes in the immune TME. Results In this model, orthotopic grafting of wildtype syngeneic keratinocytes expressing both v-ras Ha and elevated levels of ΔNp63α consistently yield carcinomas in syngeneic hosts, while cells expressing v-ras Ha alone yield predominantly papillomas. We found that polymorphonuclear (PMN) myeloid cells, experimentally validated to be immunosuppressive and thus representing myeloid-derived suppressor cells (PMN-MDSCs), were significantly recruited into the TME of carcinomas arising early following orthotopic grafting of ΔNp63α/v-ras Ha -expressing keratinocytes. ΔNp63α/v-ras Ha -driven carcinomas expressed higher levels of chemokines implicated in recruitment of MDSCs compared to v-ras Ha -initiated tumors, providing a heretofore undescribed link between ΔNp63α/HRAS-driven carcinomas and the development of an immunosuppressive TME. Conclusion These results support the utilization of a genetic carcinogenesis model harboring specific genomic drivers of malignancy to study mechanisms underlying the development of local immunosuppression.

60 APPLIED LIFE SCIENCES↗

OperonSEQer: A set of machine-learning algorithms with threshold voting for detection of operon pairs using short-read RNA-sequencing data

Operon prediction in prokaryotes is critical not only for understanding the regulation of endogenous gene expression, but also for exogenous targeting of genes using newly developed tools such as CRISPR-based gene modulation. A number of methods have used transcriptomics data to predict operons, based on the premise that contiguous genes in an operon will be expressed at similar levels. While promising results have been observed using these methods, most of them do not address uncertainty caused by technical variability between experiments, which is especially relevant when the amount of data available is small. In addition, many existing methods do not provide the flexibility to determine the stringency with which genes should be evaluated for being in an operon pair. We present OperonSEQer, a set of machine learning algorithms that uses the statistic and p-value from a non-parametric analysis of variance test (Kruskal-Wallis) to determine the likelihood that two adjacent genes are expressed from the same RNA molecule. We implement a voting system to allow users to choose the stringency of operon calls depending on whether your priority is high recall or high specificity. In addition, we provide the code so that users can retrain the algorithm and re-establish hyperparameters based on any data they choose, allowing for this method to be expanded as additional data is generated. We show that our approach detects operon pairs that are missed by current methods by comparing our predictions to publicly available long-read sequencing data. OperonSEQer therefore improves on existing methods in terms of accuracy, flexibility, and adaptability.

59 BASIC BIOLOGICAL SCIENCES↗

Identifying strengths and weaknesses of methods for computational network inference from single-cell RNA-seq data

Single-cell RNA-sequencing (scRNA-seq) offers unparalleled insight into the transcriptional programs of different cellular states by measuring the transcriptome of thousands of individual cells. An emerging problem in the analysis of scRNA-seq is the inference of transcriptional gene regulatory networks and a number of methods with different learning frameworks have been developed to address this problem. Here, we present an expanded benchmarking study of eleven recent network inference methods on seven published scRNA-seq datasets in human, mouse, and yeast considering different types of gold standard networks and evaluation metrics. We evaluate methods based on their computing requirements as well as on their ability to recover the network structure. We find that, while most methods have a modest recovery of experimentally derived interactions based on global metrics such as Area Under the Precision Recall curve, methods are able to capture targets of regulators that are relevant to the system under study. Among the top performing methods that use only expression were SCENIC, PIDC, MERLIN or Correlation. Addition of prior biological knowledge and the estimation of transcription factor activities resulted in the best overall performance with the Inferelator and MERLIN methods that use prior knowledge outperforming methods that use expression alone. We found that imputation for network inference did not improve network inference accuracy and could be detrimental. Comparisons of inferred networks for comparable bulk conditions showed that the networks inferred from scRNA-seq datasets are often better or at par with the networks inferred from bulk datasets. Our analysis should be beneficial in selecting methods for network inference. At the same time, this highlights the need for improved methods and better gold standards for regulatory network inference from scRNAseq datasets.

59 BASIC BIOLOGICAL SCIENCES↗

Modulation of Plant and Fungal Gene Expression Upon Cd Exposure and Symbiosis in Ericoid Mycorrhizal Vaccinium myrtillus

The success of Ericaceae in stressful habitats enriched in heavy metals has been ascribed to the distinctive abilities of their mycorrhizal fungal partners to withstand heavy metal stress and to enhance metal tolerance in the host plant. Whereas heavy metal tolerance has been extensively investigated in some ericoid mycorrhizal (ERM) fungi, the molecular and cellular mechanisms that extend tolerance to the host plant are currently unknown. Here, we show a reduced Cd content in Cd-exposed mycorrhizal roots of Vaccinium myrtillus colonized by a metal tolerant isolate of the fungus Oidiodendron maius as compared to non-mycorrhizal roots. To better understand this phenotype, we applied Next Generation Sequencing technologies to analyze gene expression in V. myrtillus and O. maius Zn grown under normal and Cd-stressed conditions, in the free living and in the mycorrhizal status. The results clearly showed that Cd had a stronger impact on plant gene expression than symbiosis, whereas fungal gene expression was mainly regulated by symbiosis. The higher abundance of transcripts coding for stress related proteins in non-mycorrhizal roots may be related to the higher Cd content. Regulated plant metal transporters have been identified that may play a role in reducing Cd content in mycorrhizal roots exposed to this metal.

Cd tolerance, transcriptomics↗