SEARCH · Engineering Papers

Results for “Gene expression”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Characterization of the Ghd8 Flowering Time Gene in a Mini-Core Collection of Miscanthus sinensis

The optimal flowering time for bioenergy crop Miscanthus is essential for environmental adaptability and biomass accumulation. However, little is known about how genes controlling flowering in other grasses contribute to flowering regulation in Miscanthus. Here, we report on the sequence characterization and gene expression of Miscanthus sinensisGhd8, a transcription factor encoding a HAP3/NF-YB DNA-binding domain, which has been identified as a major quantitative trait locus in rice, with pleiotropic effects on grain yield, heading date and plant height. In M. sinensis, we identified two homoeologous loci, MsiGhd8A located on chromosome 13 and MsiGhd8B on chromosome 7, with one on each of this paleo-allotetraploid species’ subgenomes. A total of 46 alleles and 28 predicted protein sequence types were identified in 12 wild-collected accessions. Several variants of MsiGhd8 showed a geographic and latitudinal distribution. Quantitative real-time PCR revealed that MsiGhd8 expressed under both long days and short days, and MsiGhd8B showed a significantly higher expression than MsiGhd8A. The comparison between flowering time and gene expression indicated that MsiGhd8B affected flowering time in response to day length for some accessions. This study provides insight into the conserved function of Ghd8 in the Poaceae, and is an important initial step in elucidating the flowering regulatory network of Miscanthus.

geographic distribution↗

DOE OSTI2024

Transcriptomic and metabolomic changes in postharvest sugarbeet roots reveal widespread metabolic changes in storage and identify genes potentially responsible for respiratory sucrose loss

Endogenous metabolism is primarily responsible for losses in sucrose content and processing quality in postharvest sugarbeet roots. The genes responsible for this metabolism and the transcriptional changes that regulate it, however, are largely unknown. To identify genes and metabolic pathways that participate in postharvest sugarbeet root metabolism and the transcriptional changes that contribute to their regulation, transcriptomic and metabolomic profiles were generated for sugarbeet roots at harvest and after 12, 40 and 120 d storage at 5 and 12°C and gene expression and metabolite concentration changes related to storage duration or temperature were identified. During storage, 8656 genes, or 34% of all expressed genes, and 225 metabolites, equivalent to 59% of detected metabolites, were altered in expression or concentration, indicating extensive transcriptional and metabolic changes in stored roots. These genes and metabolites contributed to a wide range of cellular and molecular functions, with carbohydrate metabolism being the function to which the greatest number of genes and metabolites classified. Because respiration has a central role in postharvest metabolism and is largely responsible for sucrose loss in sugarbeet roots, genes and metabolites involved in and correlated to respiration were identified. Seventy-five genes participating in respiration were differentially expressed during storage, including two bidirectional sugar transporter SWEET17 genes that highly correlated with respiration rate. Weighted gene co-expression network analysis identified 1896 additional genes that positively correlated with respiration rate and predicted a pyruvate kinase gene to be a central regulator or biomarker for respiration rate. Overall, these results reveal the extensive and diverse physiological and metabolic changes that occur in stored sugarbeet roots and identify genes with potential roles as regulators or biomarkers for respiratory sucrose loss.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Bioactivity Profiling of Chemical Mixtures for Hazard Characterization

Abstract The assessment and regulation of chemical toxicity to protect human health and the environment are done one chemical at a time and seldom at environmentally relevant concentrations. However, chemicals are found in the environment as mixtures, and their toxicity is largely unknown. Understanding the hazard posed by chemicals within the mixture is critical to enforce protective measures. Here, we demonstrate the application of bioactivity profiling of environmental water samples using the sentinel and ecotoxicology model species Daphnia to reveal the biomolecular response induced by exposure to real-world mixtures. We exposed a Daphnia strain to 30 sampled waters of the Chaobai River and measured the gene expression response profiles. Using a multiblock correlation analysis, we establish correlations between chemical mixtures identified in 30 water samples with gene expression patterns induced by these chemical mixtures. We identified 80 metabolic pathways putatively activated by mixtures of inorganic ions, heavy metals, polycyclic aromatic hydrocarbons, industrial chemicals, and a set of biocides, pesticides, and pharmacologically active substances. Our data-driven approach discovered both known bioactivity signatures with previously described modes of action and new pathways linked to undiscovered potential hazards. This study demonstrates the feasibility of reducing the complexity of real-world mixture toxicity to characterize the biomolecular effects of a defined number of chemical components based on gene expression monitoring of the sentinel species Daphnia.

Engineering↗

DOE OSTI2026

Genetic and epigenetic signatures of Populus trichocarpa in response to abiotic stresses

Plants respond to abiotic stresses such as drought, heat, and salinity through both shared and stress-specific regulatory pathways. However, the role of epigenetic mechanisms, particularly DNA methylation, in modulating these responses is still underexplored. Here, we integrated transcriptome profiling with whole-genome bisulfite sequencing (WGBS) to investigate the dynamic relationship between gene expression and DNA methylation in Populus trichocarpa during brief early-stage exposure to abiotic stress. Each stress elicited distinct transcriptional and methylation signatures; however, a subset of conserved stress-responsive genes was commonly regulated across treatments. We found that ∼5-6% of differentially expressed genes also displayed differential methylation patterns, suggesting a coordinated role for DNA methylation in regulating gene expression. Motif enrichment analysis of differentially methylated regions revealed binding sites of key transcription factor families, including ERF, bHLH, and ABF, highlighting potential role for methylation in modulating transcription factor targeting. Furthermore, we identified stress-inducible molecular markers with potential applications in early stress detection and functional dissection of gene regulatory pathways. Together, these findings provide new insights into the coordinated genetic and epigenetic responses to abiotic stress and provide a foundation for developing biosensors and breeding strategies to enhance stress resilience woody plants.

DNA methylation↗

DOE OSTI2022

The MuvB complex binds and stabilizes nucleosomes downstream of the transcription start site of cell-cycle dependent genes

The chromatin architecture in promoters is thought to regulate gene expression, but it remains uncertain how most transcription factors (TFs) impact nucleosome position. The MuvB TF complex regulates cell-cycle dependent gene-expression and is critical for differentiation and proliferation during development and cancer. MuvB can both positively and negatively regulate expression, but the structure of MuvB and its biochemical function are poorly understood. Here we determine the overall architecture of MuvB assembly and the crystal structure of a subcomplex critical for MuvB function in gene repression. We find that the MuvB subunits LIN9 and LIN37 function as scaffolding proteins that arrange the other subunits LIN52, LIN54 and RBAP48 for TF, DNA, and histone binding, respectively. Biochemical and structural data demonstrate that MuvB binds nucleosomes through an interface that is distinct from LIN54-DNA consensus site recognition and that MuvB increases nucleosome occupancy in a reconstituted promoter. We find in arrested cells that MuvB primarily associates with a tightly positioned +1 nucleosome near the transcription start site (TSS) of MuvB-regulated genes. These results support a model that MuvB binds and stabilizes nucleosomes just downstream of the TSS on its target promoters to repress gene expression.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Quantitative trait locus (QTL) mapping and transcriptome profiling identify QTLs and candidate genes associated with heat stress response during reproductive development in Camelina sativa

Camelina sativa (L.) Crantz is a low-input oilseed crop that has great potential in providing sustainable feedstock for biofuels and bioproducts. Climate change is threatening production of camelina with rising global temperatures. Elucidating the genetic response to high temperatures is essential for successful breeding of heat-tolerant camelina varieties. Here, we report a combinatorial approach to identifying candidate genes associated with heat stress by quantitative trait locus (QTL) mapping and comparative transcriptome profiling. A population of recombinant inbred lines (RILs) was grown in a controlled growth chamber under the high-temperature regimes for 14 days beginning at the onset of the reproductive stage. Several traits related to seed production were evaluated at maturity. The QTL analysis identified several regions with co-located traits on chromosomes 8, 10, and 12. Two RILs with contrasting phenotypic responses to heat stress were chosen for gene expression profiling via RNA sequencing. Multiple pathways and genes were found to be strongly affected by heat stress, and many genes expressed differently between the two RILs. Several genes identified within the QTL regions were considered strong candidates that may control heat tolerance during reproduction in camelina. These studies provide resources for future studies that may assist in improving the heat tolerance of camelina.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

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↗

DOE OSTI2022

Harmonizing model organism data in the Alliance of Genome Resources

The Alliance of Genome Resources (the Alliance) is a combined effort of 7 knowledgebase projects: Saccharomyces Genome Database, WormBase, FlyBase, Mouse Genome Database, the Zebrafish Information Network, Rat Genome Database, and the Gene Ontology Resource. The Alliance seeks to provide several benefits: better service to the various communities served by these projects; a harmonized view of data for all biomedical researchers, bioinformaticians, clinicians, and students; and a more sustainable infrastructure. The Alliance has harmonized cross-organism data to provide useful comparative views of gene function, gene expression, and human disease relevance. The basis of the comparative views is shared calls of orthology relationships and the use of common ontologies. The key types of data are alleles and variants, gene function based on gene ontology annotations, phenotypes, association to human disease, gene expression, protein–protein and genetic interactions, and participation in pathways. The information is presented on uniform gene pages that allow facile summarization of information about each gene in each of the 7 organisms covered (budding yeast, roundworm Caenorhabditis elegans, fruit fly, house mouse, zebrafish, brown rat, and human). The harmonized knowledge is freely available on the alliancegenome.org portal, as downloadable files, and by APIs. We expect other existing and emerging knowledge bases to join in the effort to provide the union of useful data and features that each knowledge base currently provides.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Mef2c regulates bone mass through Sost -dependent and -independent mechanisms

We report Mef2c is a transcription factor that mediates key cellular behaviors that promote endochondral ossification and bone formation. Previously, Mef2c has been shown to regulate Sost transcription via its osteocyte-specific enhancer, ECR5, and conditional deletions of Mef2 cfl/fl with either Col1-Cre or Dmp1-Cre produced generalized high bone mass (HBM) consistent with Van Buchem Disease phenotypes. However, Sost -/- ; Mef2 cfl/fl ; Dmp1-Cre mice produced a significantly higher bone mass phenotype that Sost-/- alone suggesting that Mef2c modulates bone mass through additional mechanisms, independent of Sost. To identify new Mef2c transcriptional targets important in bone metabolism, we profiled gene expression by single-cell RNA sequencing in subpopulations of cells isolated from Mef2 cfl/fl ; Dmp1-Cre and Mef2 cfl/fl ; Bglap-Cre femurs, both strains exhibiting similar high bone mass phenotypes. However, we found Mef2cfl/fl; Bglap-Cre to also display a growth plate defect characterized by an expansion of several osteoprogenitor subpopulations. Differential gene expression analysis identified a total of 96 up- and 2434 down- regulated genes in Mef2 cfl/fl ; Bglap-Cre and 176 up- and 1041 down- regulated genes in Mef2 cfl/fl ; Dmp1-Cre bone cell subpopulations compared to wildtype mice. Mef2c deletion affected the transcriptomes across several cell types including mesenchymal progenitors (MP), osteoprogenitors (OSP), osteoblast (OB), and osteocyte (OCY) subpopulations. Several energy metabolism genes such as Uqcrb, Ndufv2, Ndufs3, Ndufa13, Ndufb9, Ndufb5, Cox6a1, Cox5a, Atp5o, Atp5g2, Atp5b, Atp5 were significantly down regulated in Mef2c-deficient OBs and OCYs, in both strains. Binding motif analysis of promoter regions of differentially expressed genes identified Mef2c binding in Bone Sialoprotein (BSP/Ibsp), a gene known to cause increased trabecular BV/TV in the femurs of Ibsp -/- mice. Immunohistochemical analysis confirmed the absence of Ibsp protein in OBs and OCYs. These findings suggests that the HBM in Sost -/- ; Mef2 cfl/fl ; Dmp1-Cre is caused by a multitude of transcriptional changes in genes that regulate bone formation, two of which are Sost and Ibsp.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Systemic immunological responses are dependent on sex and ovarian hormone presence following acute inhaled woodsmoke exposure

Background: Rural regions of the western United States have experienced a noticeable surge in both the frequency and severity of acute wildfire events, which brings significant challenges to both public safety and environmental conservation efforts, with impacts felt globally. Identifying factors contributing to immune dysfunction, including endocrinological phenotypes, is essential to understanding how hormones may influence toxicological susceptibility. Methods: This exploratory study utilized male and female C57BL/6 mice as in vivo models to investigate distinct responses to acute woodsmoke (WS) exposure with a focus on sex-based differences. In a second set of investigations, two groups were established within the female mouse cohort. In one group, mice experienced ovariectomy (OVX) to simulate an ovarian hormone-deficient state similar to surgical menopause, while the other group received Sham surgery as controls, to investigate the mechanistic role of ovarian hormone presence in driving immune dysregulation following acute WS exposure. Each experimental cohort followed a consecutive 2-day protocol with daily 4-h exposure intervals under two conditions: control HEPA-filtered air (FA) and acute WS to simulate an acute wildfire episode. Results: Metals analysis of WS particulate matter (PM) revealed significantly increased levels of 63 Cu, 182 W, 208 Pb, and 238 U, compared to filtered air (FA) controls, providing insights into the specific metal components most impacted by the changing dynamics of wildfire occurrences in the region. Male and female mice exhibited diverse patterns in lung mRNA cytokine expression following WS exposure, with males showing downregulation and females displaying upregulation, notably for IL-1β, TNF-α, CXCL-1, CCL-5, TGF-β, and IL-6. After acute WS exposure, there were notable differences in the responses of macrophages, neutrophils, and bronchoalveolar lavage (BAL) cytokines IL-10, IL-6, IL-1β, and TNF-α. Significant diverse alterations were observed in BAL cytokines, specifically IL-1β, IL-10, IL-6, and TNF-α, as well as in the populations of immune cells, such as macrophages and polymorphonuclear leukocytes, in both Sham and OVX mice, following acute WS exposure. These findings elucidated the profound influence of hormonal changes on inflammatory outcomes, delineating substantial sex-related differences in immune activation and revealing altered immune responses in OVX mice due to ovarian hormone deficiency. In addition, the flow cytometry analysis highlighted the complex interaction between OVX surgery, acute WS exposure, and their collective impact on immune cell populations within the hematopoietic bone marrow niche. Conclusions: In summary, both male and female mice, alongside females subjected to OVX and those who had sham surgery, exhibit significant variations in the expression of proinflammatory cytokines, chemokines, lung mRNA gene expression, and related functional networks linked to signaling pathways. These differences potentially act as mediators of sex-specific and hormonal influences in the systemic inflammatory response to acute WS exposure during a wildfire event. Understanding the regulatory roles of genes expressed differentially under environmental stressors holds considerable implications, aiding in identifying sex-specific therapeutic targets for addressing acute lung inflammation and injury.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Symmetric adenine methylation is an essential DNA modification in the early-diverging fungus Rhizopus microsporus

The discovery of N6-methyladenine (6mA) in eukaryotic genomes, typically found in prokaryotic DNA, has revolutionized epigenetics. Here, we show that symmetric 6mA is essential in the early diverging fungus Rhizopus microsporus, as the absence of the MT-A70 complex (MTA1c) responsible for this modification results in a lethal phenotype. 6mA is present in 70% of the genes, correlating with the presence of H3K4me3 and H2A.Z in open euchromatic regions. This modification is found predominantly in nucleosome linker regions, influencing the nucleosome positioning around the transcription start sites of highly expressed genes. Controlled downregulation of MTA1c reduces symmetric 6mA sites affecting nucleosome positioning and histone modifications, leading to altered gene expression, which is likely the cause of the severe phenotypic changes observed. Our study highlights the indispensable role of the DNA 6mA in a multicellular organism and delineates the mechanisms through which this epigenetic mark regulates gene expression in a eukaryotic genome.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Building customizable auto-luminescent luciferase-based reporters in plants

Bioluminescence is a powerful biological signal that scientists have repurposed as a reporter for gene expression in plants and animals. However, there are downsides associated with the need to provide a substrate to these reporters, including its high cost and non-uniform tissue penetration. In this work we reconstitute a fungal bioluminescence pathway (FBP) in planta using a composable toolbox of parts. We demonstrate that the FBP can create luminescence across various tissues in a broad range of plants without external substrate addition. We also show how our toolbox can be used to deploy the FBP in planta to build auto-luminescent reporters for the study of gene-expression and hormone fluxes. A low-cost imaging platform for gene expression profiling is also described. These experiments lay the groundwork for future construction of programmable auto-luminescent plant traits, such as light driven plant-pollinator interactions or light emitting plant-based sensors.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Transcriptomics reveal a mechanism of niche defense: two beneficial root endophytes deploy an antimicrobial GH18-CBM5 chitinase to protect their hosts

Effector secretion is crucial for root endophytes to establish and protect their ecological niche. We used time-resolved transcriptomics to monitor effector gene expression dynamics in two closely related Sebacinales, Serendipita indica and Serendipita vermifera, during symbiosis with three plant species, competition with the phytopathogenic fungus Bipolaris sorokiniana, and cooperation with root-associated bacteria. We observed increased effector gene expression in response to biotic interactions, particularly with plants, indicating their importance in host colonization. Some effectors responded to both plants and microbes, suggesting dual roles in intermicrobial competition and plant–microbe interactions. A subset of putative antimicrobial effectors, including a GH18-CBM5 chitinase, was induced exclusively by microbes. Functional analyses of this chitinase revealed its antimicrobial and plant-protective properties. We conclude that dynamic effector gene expression underpins the ability of Sebacinales to thrive in diverse ecological niches with a single fungal chitinase contributing substantially to niche defense.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Perimatrix of middle ear cholesteatoma: A granulation tissue with a specific transcriptomic signature

Objectives/Hypothesis: To establish comprehensive transcriptomic profiles of cholesteatoma perimatrix tissue and granulation tissue from chronic otitis media (COM) that did not develop cholesteatoma, which can indicate molecular pathways involved in the cholesteatoma perimatrix pathology and invasiveness. Study Design: Retrospective Case Series. Methods: Transcriptome data were obtained from cholesteatoma perimatrix tissue and COM granulation tissue by an Illumina iScan microarray. Differentially expressed genes (DEGs) were subsequently analyzed using both bioinformatical functional annotation and network analysis. Expression of candidate genes (MMP9 and LCN2) was validated by quantitative reverse transcription-polymerase chain reaction (qRT-PCR) on a larger group of samples. Results: Analysis of the transcriptome led to the identification of 169 differentially expressed genes between investigated tissues. Bioinformatic analysis suggested that most significant biological processes involving DEGs were previously described in cholesteatoma pathology. Network analysis identified ERBB2, TFAP2A, and TP63 as major hubs of the DEGs molecular network. Furthermore, it was observed that the cellular component most significantly enriched in DEGs was extracellular space containing 47 DEGs. Using qRT-PCR, it was confirmed that mRNA levels of the major extracellular hub (MMP9) are increased, whereas its interacting molecule (LCN2) mRNA levels were decreased in cholesteatoma perimatrix tissue compared to COM granulation tissue. Conclusions: The current study approach offers an overall look at molecular mechanisms that describe the cholesteatoma entity by focusing exclusively on the perimatrix processes in comparison to COM granulation tissue. The observed differences in gene expression between cholesteatoma perimatrix and COM granulation tissue could suggest novel markers potentially influenced by the perimatrix–matrix molecular interplay, which is not present in COM without cholesteatoma. Level of Evidence: NA. Laryngoscope, 130:E220–E227, 2020. © 2019 The American Laryngological, Rhinological and Otological Society, Inc.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

TransPlatformer

We propose TransPlatformer for translating toxicogenomics from one platform to another. Transcriptomic profiling has evolved through multiple generations of technology, from microarrays (e.g., Affymetrix, CodeLink) to more recent high-throughput sequencing and targeted panels such as S1500+. Microarrays, which dominated gene expression studies in the early 2000s, provided affordable and high-throughput transcript quantification but suffered from cross-hybridization issues and limited dynamic range . RNA-Seq, introduced in the late 2000s, revolutionized transcriptomics by enabling unbiased and comprehensive gene expression analysis, albeit at higher costs and computational demands . Despite advances, many studies rely on historical microarray data, necessitating the translation of legacy data into modern platforms to ensure continuity and comparability. This translation is complicated by factors such as platform-specific probe design, differences in transcript coverage, and batch effects . Existing methods for cross-platform mapping include statistical normalization, machine learning models, and biological anchoring approaches. The ability to translate transcriptomic data between platforms has broad implications, including enhanced meta-analyses, improved toxicological modeling, and better integration of historical datasets with contemporary research. TransPlatformer seeks to contribute to this effort by evaluating translation methodologies and proposing novel strategies to improve cross-platform gene expression harmonization. In this repository there are code examples for TransPlatformer implementation

Cong, Guojing↗

DOE OSTI2026

A Fluorescence‐Based Transient Expression Assay for the Analysis of Upstream Open Reading Frames in Plants

Upstream open reading frames (uORFs) are regulatory elements present in the 5′ leaders of mRNA that can significantly impact downstream gene expression in eukaryotes. In crop engineering, editing of uORFs can provide an avenue to upregulate expression of native genes without the need to add persistent transgenic copies. Even with genome-wide methods to identify translated uORFs such as ribosome profiling, their functional characterization depends on validation through reporter gene assays and mutagenesis studies. Current screening methods for plants use luciferases or protoplasts to measure differential gene expression between wild-type and mutated transcript leaders, which requires tissue processing and/or substrate addition. Here, we present a time- and cost-efficient alternative to investigate transcript leaders by co-expression of two fluorescent proteins in Nicotiana benthamiana leaf tissue and test our assay on genes involved in photoprotection, editing of which could provide a pathway to increase CO 2 assimilation during sun–shade transitions.

Nicotiana benthamiana↗

DOE OSTI2022

Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach

Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression. Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to previous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding. Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Integrated transcriptomic and metabolomic investigation of the genes and metabolites involved in swine follicular cyst formation

Follicular cysts are a common reproductive disorder in mammals that is usually caused by stress. However, the pathogenesis of follicular cysts in sows remains unclear. To provide new insights into the mechanisms of follicular cyst formation in pigs, we conducted a combined transcriptomic and metabolomic analysis on theca interna and mural granulosa cells of follicular cysts and mature follicles. We identified 2,533 up-regulated and 1,355 down-regulated genes in follicular cysts, compared with mature follicles. These differentially expressed genes were mainly found in signaling pathways related to tumor formation and cortisol synthesis and secretion as shown by Ingenuity Pathway Analysis, which predicted 4,362 upstream regulatory factors. The combined gene expression and pathway analysis identified the following genes as potential biomarkers for porcine follicular cysts: cytochrome P450 family 2 subfamily C polypeptide 18, L-lactate dehydrogenase, carbamoyl-phosphate synthase, fibroblast growth factor 7, integrin binding sialoprotein, interleukin 23 receptor, prolactin receptor, epiregulin, interleukin 1 receptor type II, arginine vasopressin receptor 1A, fibroblast growth factor 10, claudin 7, G Protein Subunit Gamma 3, cholecystokinin B receptor and cytosolic phospholipase A2 . Metabolomics analysis found significant differences in 87 metabolites, which were enriched in unsaturated fatty acid biosynthesis, and sphingolipid signaling pathways. These results provide valuable information on the molecular mechanisms of follicular cyst formation, which may facilitate the development of new therapeutics to prevent and treat follicular cysts.

Dai, Jiage↗