Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural genomics”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Seed coat transcriptomic profiling of 5-593, a genotype important for genetic studies of seed coat color and patterning in common bean ( Phaseolus vulgaris L.)

Common bean (Phaseolus vulgaris L.) market classes have distinct seed coat colors, which are directly related to the diverse flavonoids found in the mature seed coat. To understand and elucidate the molecular mechanisms underlying the regulation of seed coat color, RNA-Seq data was collected from the black bean 5-593 and used for a differential gene expression and enrichment analysis from four different seed coat color development stages. 5-593 carries dominant alleles for 10 of the 11 major genes that control seed coat color and expression and has historically been used to develop introgression lines used for seed coat genetic analysis. Pairwise comparison among the four stages identified 6,294 differentially expressed genes (DEGs) varying from 508 to 5,780 DEGs depending on the compared stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that phenylpropanoid biosynthesis, flavonoid biosynthesis, and plant hormone signal transduction comprised the principal pathways expressed during bean seed coat pigment development. Transcriptome analysis suggested that most structural genes for flavonoid biosynthesis and some potential regulatory genes were significantly differentially expressed. Further studies detected 29 DEGs as important candidate genes governing the key enzymatic flavonoid biosynthetic pathways for common bean seed coat color development. Additionally, four gene models, Pv5-593.02G016100, 593.02G078700, Pv5-593.02G090900, and Pv5-593.06G121300, encode MYB-like transcription factor family protein were identified as strong candidate regulatory genes in anthocyanin biosynthesis which could regulate the expression levels of some important structural genes in flavonoid biosynthesis pathway. These findings provide a framework to draw new insights into the molecular networks underlying common bean seed coat pigment development.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Genetically correlated leaf tensile and morphological traits are driven by growing season length in a widespread perennial grass

Leaf tensile resistance, a leaf's ability to withstand pulling forces, is an important determinant of plant ecological strategies. One potential driver of leaf tensile resistance is growing season length. When growing seasons are long, strong leaves, which often require more time and resources to construct than weak leaves, may be more advantageous than when growing seasons are short. Growing season length and other ecological conditions may also impact the morphological traits that underlie leaf tensile resistance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Visualization of conformational changes and membrane remodeling leading to genome delivery by viral class-II fusion machinery

Chikungunya virus (CHIKV) is a human pathogen that delivers its genome to the host cell cytoplasm through endocytic low pH-activated membrane fusion mediated by class-II fusion proteins. Though structures of prefusion, icosahedral CHIKV are available, structural characterization of virion interaction with membranes has been limited. Here, we have used cryo-electron tomography to visualize CHIKV’s complete membrane fusion pathway, identifying key intermediary glycoprotein conformations coupled to membrane remodeling events. Using sub-tomogram averaging, we elucidate features of the low pH-exposed virion, nucleocapsid and full-length E1-glycoprotein’s post-fusion structure. Contrary to class-I fusion systems, CHIKV achieves membrane apposition by protrusion of extended E1-glycoprotein homotrimers into the target membrane. The fusion process also features a large hemifusion diaphragm that transitions to a wide pore for intact nucleocapsid delivery. Our analyses provide comprehensive ultrastructural insights into the class-II virus fusion system function and direct mechanistic characterization of the fundamental process of protein-mediated membrane fusion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Small and mighty: adaptation of superphylum Patescibacteria to groundwater environment drives their genome simplicity

The newly defined superphylum Patescibacteria such as Parcubacteria (OD1) and Microgenomates (OP11) has been found to be prevalent in groundwater, sediment, lake, and other aquifer environments. Recently increasing attention has been paid to this diverse superphylum including > 20 candidate phyla (a large part of the candidate phylum radiation, CPR) because it refreshed our view of the tree of life. However, adaptive traits contributing to its prevalence are still not well known. Results Here, we investigated the genomic features and metabolic pathways of Patescibacteria in groundwater through genome-resolved metagenomics analysis of >600 Gbp sequence data. We observed that, while the members of Patescibacteria have reduced genomes (~ 1 Mbp) exclusively, functions essential to growth and reproduction such as genetic information processing were retained. Surprisingly, they have sharply reduced redundant and nonessential functions, including specific metabolic activities and stress response systems. The Patescibacteria have ultra-small cells and simplified membrane structures, including flagellar assembly, transporters, and two-component systems. Despite the lack of CRISPR viral defense, the bacteria may evade predation through deletion of common membrane phage receptors and other alternative strategies, which may explain the low representation of prophage proteins in their genomes and lack of CRISPR. By establishing the linkages between bacterial features and the groundwater environmental conditions, our results provide important insights into the functions and evolution of this CPR group.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Postglacial migration shaped the genomic diversity and global distribution of the wild ancestor of lager-brewing hybrids

The wild, cold-adapted parent of hybrid lager-brewing yeasts, Saccharomyces eubayanus, has a complex and understudied natural history. The exploration of this diversity can be used both to develop new brewing applications and to enlighten our understanding of the dynamics of yeast evolution in the wild. Here, we integrate whole genome sequence and phenotypic data of 200 S. eubayanus strains, the largest collection known to date. S. eubayanus has a multilayered population structure, consisting of two major populations that are further structured into six subpopulations. Four of these subpopulations are found exclusively in the Patagonian region of South America; one is found predominantly in Patagonia and sparsely in Oceania and North America; and one is specific to the Holarctic ecozone. Plant host associations differed between subpopulations and between S. eubayanus and its sister species, Saccharomyces uvarum. S. eubayanus is most abundant and genetically diverse in northern Patagonia, where some locations harbor more genetic diversity than is found outside of South America, suggesting that northern Patagonia east of the Andes was a glacial refugium for this species. All but one subpopulation shows isolation-by-distance, and gene flow between subpopulations is low. However, there are strong signals of ancient and recent outcrossing, including two admixed lineages, one that is sympatric with and one that is mostly isolated from its parental populations. Using our extensive biogeographical data, we build a robust model that predicts all known and a handful of additional regions of the globe that are climatically suitable for S. eubayanus, including Europe where host accessibility and competitive exclusion by other Saccharomyces species may explain its continued elusiveness. We conclude that this industrially relevant species has rich natural diversity with many factors contributing to its complex distribution and natural history.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Structure of RADX and mechanism for regulation of RAD51 nucleofilaments

Replication fork reversal is a fundamental process required for resolution of encounters with DNA damage. A key step in the stabilization and eventual resolution of reversed forks is formation of RAD51 nucleoprotein filaments on exposed single strand DNA (ssDNA). To avoid genome instability, RAD51 filaments are tightly controlled by a variety of positive and negative regulators. RADX (RPA-related RAD51-antagonist on the X chromosome) is a recently discovered negative regulator that binds tightly to ssDNA, directly interacts with RAD51, and regulates replication fork reversal and stabilization in a context-dependent manner. Here, we present a structure-based investigation of RADX’s mechanism of action. Mass photometry experiments showed that RADX forms multiple oligomeric states in a concentration-dependent manner, with a predominance of trimers in the presence of ssDNA. The structure of RADX, which has no structurally characterized orthologs, was determined ab initio by cryo-electron microscopy (cryo-EM) from maps in the 2 to 4 Å range. The structure reveals the molecular basis for RADX oligomerization and the coupled multi-valent binding of ssDNA binding. The interaction of RADX with RAD51 filaments was imaged by negative stain EM, which showed a RADX oligomer at the end of filaments. Based on these results, we propose a model in which RADX functions by capping and restricting the end of RAD51 filaments.

Science & Technology - Other Topics↗
DOE OSTI2020

Human XPG nuclease structure, assembly, and activities with insights for neurodegeneration and cancer from pathogenic mutations

Xeroderma pigmentosum group G (XPG) protein is both a functional partner in multiple DNA damage responses (DDR) and a pathway coordinator and structure-specific endonuclease in nucleotide excision repair (NER). Different mutations in the XPG gene ERCC5 lead to either of two distinct human diseases: Cancer-prone xeroderma pigmentosum (XP-G) or the fatal neurodevelopmental disorder Cockayne syndrome (XP-G/CS). To address the enigmatic structural mechanism for these differing disease phenotypes and for XPG’s role in multiple DDRs, here we determined the crystal structure of human XPG catalytic domain (XPGcat), revealing XPG-specific features for its activities and regulation. Furthermore, XPG DNA binding elements conserved with FEN1 superfamily members enable insights on DNA interactions. Notably, all but one of the known pathogenic point mutations map to XPGcat, and both XP-G and XP-G/CS mutations destabilize XPG and reduce its cellular protein levels. Mapping the distinct mutation classes provides structure-based predictions for disease phenotypes: Residues mutated in XP-G are positioned to reduce local stability and NER activity, whereas residues mutated in XP-G/CS have implied long-range structural defects that would likely disrupt stability of the whole protein, and thus interfere with its functional interactions. Combined data from crystallography, biochemistry, small angle X-ray scattering, and electron microscopy unveil an XPG homodimer that binds, unstacks, and sculpts duplex DNA at internal unpaired regions (bubbles) into strongly bent structures, and suggest how XPG complexes may bind both NER bubble junctions and replication forks. Collective results support XPG scaffolding and DNA sculpting functions in multiple DDR processes to maintain genome stability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Phylogenomic analysis of 589 metagenome-assembled genomes encompassing all major prokaryotic lineages from the gut of higher termites

“Higher” termites have been able to colonize all tropical and subtropical regions because of their ability to digest lignocellulose with the aid of their prokaryotic gut microbiota. Over the last decade, numerous studies based on 16S rRNA gene amplicon libraries have largely described both the taxonomy and structure of the prokaryotic communities associated with termite guts. Host diet and microenvironmental conditions have emerged as the main factors structuring the microbial assemblages in the different gut compartments. Additionally, these molecular inventories have revealed the existence of termite-specific clusters that indicate coevolutionary processes in numerous prokaryotic lineages. However, for lack of representative isolates, the functional role of most lineages remains unclear. We reconstructed 589 metagenome-assembled genomes (MAGs) from the different gut compartments of eight higher termite species that encompass 17 prokaryotic phyla. By iteratively building genome trees for each clade, we significantly improved the initial automated assignment, frequently up to the genus level. We recovered MAGs from most of the termite-specific clusters in the radiation of, for example, Planctomycetes, Fibrobacteres, Bacteroidetes, Euryarchaeota, Bathyarchaeota, Spirochaetes, Saccharibacteria, and Firmicutes, which to date contained only few or no representative genomes. Moreover, the MAGs included abundant members of the termite gut microbiota. This dataset represents the largest genomic resource for arthropod-associated microorganisms available to date and contributes substantially to populating the tree of life. More importantly, it provides a backbone for studying the metabolic potential of the termite gut microbiota, including the key members involved in carbon and nitrogen biogeochemical cycles, and important clues that may help cultivating representatives of these understudied clades.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Multi-layered heterochromatin interaction as a switch for DIM2-mediated DNA methylation

Functional crosstalk between DNA methylation, histone H3 lysine-9 trimethylation (H3K9me3) and heterochromatin protein 1 (HP1) is essential for proper heterochromatin assembly and genome stability. However, how repressive chromatin cues guide DNA methyltransferases for region-specific DNA methylation remains largely unknown. Here, we report structure-function characterizations of DNA methyltransferase Defective-In-Methylation-2 (DIM2) in Neurospora . The DNA methylation activity of DIM2 requires the presence of both H3K9me3 and HP1. Our structural study reveals a bipartite DIM2-HP1 interaction, leading to a disorder-to-order transition of the DIM2 target-recognition domain that is essential for substrate binding. Furthermore, the structure of DIM2-HP1-H3K9me3-DNA complex reveals a substrate-binding mechanism distinct from that for its mammalian orthologue DNMT1. In addition, the dual recognition of H3K9me3 peptide by the DIM2 RFTS and BAH1 domains allosterically impacts the DIM2-substrate binding, thereby controlling DIM2-mediated DNA methylation. Together, this study uncovers how multiple heterochromatin factors coordinately orchestrate an activity-switching mechanism for region-specific DNA methylation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

‘ Candidatus Liberibacter asiaticus’-Encoded BCP Peroxiredoxin Suppresses Lipopolysaccharide-Mediated Defense Signaling and Nitrosative Stress In Planta

The lipopolysaccharides (LPS) of gram-negative bacteria trigger a nitrosative and oxidative burst in both animals and plants during pathogen invasion. Liberibacter crescens strain BT-1 is a surrogate for functional genomic studies of the uncultured pathogenic ‘Candidatus Liberibacter’ spp. that are associated with severe diseases such as citrus greening and potato zebra chip. Structural determination of L. crescens LPS revealed the presence of a very long chain fatty acid modification. L. crescens LPS pretreatment suppressed growth of Xanthomonas perforans on nonhost tobacco (Nicotiana benthamiana) and X. citri subsp. citri on host orange (Citrus sinensis), confirming bioactivity of L. crescens LPS in activation of systemic acquired resistance (SAR). L. crescens LPS elicited a rapid burst of nitric oxide (NO) in suspension cultured tobacco cells. Pharmacological inhibitor assays confirmed that arginine-utilizing NO synthase (NOS) activity was the primary source of NO generation elicited by L. crescens LPS. LPS treatment also resulted in biological markers of NO-mediated SAR activation, including an increase in the glutathione pool, callose deposition, and activation of the salicylic acid and azelaic acid (AzA) signaling networks. Transient expression of ‘Ca. L. asiaticus’ bacterioferritin comigratory protein (BCP) peroxiredoxin in tobacco compromised AzA signaling, a prerequisite for LPS-triggered SAR. Western blot analyses revealed that ‘Ca. L. asiaticus’ BCP peroxiredoxin prevented peroxynitrite-mediated tyrosine nitration in tobacco. ‘Ca. L. asiaticus’ BCP peroxiredoxin (i) attenuates NO-mediated SAR signaling and (ii) scavenges peroxynitrite radicals, which would facilitate repetitive cycles of ‘Ca. L. asiaticus’ acquisition and transmission by fecund psyllids throughout the limited flush period in citrus.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Novel anaerobic selenium oxyanion reducers native to FGD wastewater for enhanced selenium removal

Biological treatment is a recognized approach for removing selenate and selenite oxyanions present in flue gas desulfurization (FGD) wastewater. However, the knowledge of the specific microbial species or communities responsible for reducing water-soluble selenium oxyanions to insoluble elemental selenium remains limited. In addition, the selenium oxyanion reduction genes and pathways have yet to be understood in these wastewaters. This study characterizes selenium oxyanion-reducing bacteria (SeRB) native to FGD wastewater, and the resulting elemental selenium particles formed. By selecting native SeRB microbes in a defined media, a novel resolution of these organisms has been achieved. This research identifies previously unrecognized selenium oxyanion-reducing capabilities in Anaerosolibacter, alongside predominant SeRB from Mesobacillus and Tepidibacillus genera. This work encompasses both 16S and metagenomic techniques to recover novel metagenome-assembled genomes, distinct to this environment. The biogenic selenium produced by these organisms was predominantly of elemental selenium, either amorphous or with a hexagonal structure. This study identifies the SeRB present in FGD wastewater and characterizes their selenium products, offering crucial insights to enhance the efficiency of biological treatment strategies and the potential of selenium recovery from this industrial waste.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Genomic analysis and identification of a novel superantigen, SargEY, in Staphylococcus argenteus isolated from atopic dermatitis lesions

During surveillance of Staphylococcus aureus in lesions from patients with atopic dermatitis (AD), we isolated Staphylococcus argenteus, a species registered in 2011 as a new member of the genus Staphylococcus and previously considered a lineage of S. aureus. Genome sequence comparisons between S. argenteus isolates and representative S. aureus clinical isolates from various origins revealed that the S. argenteus genome from AD patients closely resembles that of S. aureus causing skin infections. We previously reported that 17%–22% of S. aureus isolated from skin infections produce staphylococcal enterotoxin Y (SEY), which predominantly induces T-cell proliferation via the T-cell receptor (TCR) Vα pathway. Complete genome sequencing of S. argenteus isolates revealed a gene encoding a protein similar to superantigen SEY, designated as SargEY, on its chromosome. Population structure analysis of S. argenteus revealed that these isolates are ST2250 lineage, which was the only lineage positive for the SEY-like gene among S. argenteus. Recombinant SargEY demonstrated immunological cross-reactivity with anti-SEY serum. SargEY could induce proliferation of human CD4 + and CD8 + T cells, as well as production of TNF-α and IFN-γ. SargEY showed emetic activity in a marmoset monkey model. S arg EY and SET (a phylogenetically close but uncharacterized SE) revealed their dependency on TCR Vα in inducing human T-cell proliferation. Additionally, TCR sequencing revealed other previously undescribed Vα repertoires induced by SEH. S arg EY and SEY may play roles in exacerbating the respective toxin-producing strains in AD.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Combinatorial Biosynthesis of Sulfated Benzenediol Lactones with a Phenolic Sulfotransferase from Fusarium graminearum PH-1

Total biosynthesis or whole-cell biocatalytic production of sulfated small molecules relies on the discovery and implementation of appropriate sulfotransferase enzymes. Although fungi are prominent biocatalysts and have been used to sulfate drug-like phenolics, no gene encoding a sulfotransferase enzyme has been functionally characterized from these organisms. Here, we identify a phenolic sulfotransferase, FgSULT1, by genome mining from the plant-pathogenic fungus Fusarium graminearum PH-1. We expressed FgSULT1 in a Saccharomyces cerevisiae chassis to modify a broad range of benzenediol lactones and their nonmacrocyclic congeners, together with an anthraquinone, with the resulting unnatural natural product (uNP) sulfates displaying increased solubility. FgSULT1 shares low similarity with known animal and plant sulfotransferases. Instead, it forms a sulfotransferase family with putative bacterial and fungal enzymes for phase II detoxification of xenobiotics and allelochemicals. Among fungi, putative FgSULT1 homologues are encoded in the genomes of Fusarium spp. and a few other genera in nonsyntenic regions, some of which may be related to catabolic sulfur recycling. Computational structure modeling combined with site-directed mutagenesis revealed that FgSULT1 retains the key catalytic residues and the typical fold of characterized animal and plant sulfotransferases. Our work opens the way for the discovery of hitherto unknown fungal sulfotransferases and provides a synthetic biological and enzymatic platform that can be adapted to produce bioactive sulfates, together with sulfate ester standards and probes for masked mycotoxins, precarcinogenic toxins, and xenobiotics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Topological structure of complex predictions

Abstract Current complex prediction models are the result of fitting deep neural networks, graph convolutional networks or transducers to a set of training data. A key challenge with these models is that they are highly parameterized, which makes describing and interpreting the prediction strategies difficult. We use topological data analysis to transform these complex prediction models into a simplified topological view of the prediction landscape. The result is a map of the predictions that enables inspection of the model results with more specificity than dimensionality-reduction methods such as tSNE and UMAP. The methods scale up to large datasets across different domains. We present a case study of a transformer-based model previously designed to predict expression levels of a piece of DNA in thousands of genomic tracks. When the model is used to study mutations in the BRCA1 gene, our topological analysis shows that it is sensitive to the location of a mutation and the exon structure of BRCA1 in ways that cannot be found with tools based on dimensionality reduction. Moreover, the topological framework offers multiple ways to inspect results, including an error estimate that is more accurate than model uncertainty. Further studies show how these ideas produce useful results in graph-based learning and image classification.

Computer Science↗
DOE OSTI2024

African Swine Fever Virus Protein–Protein Interaction Prediction

The African swine fever virus (ASFV) is an often deadly disease in swine and poses a threat to swine livestock and swine producers. With its complex genome containing more than 150 coding regions, developing effective vaccines for this virus remains a challenge due to a lack of basic knowledge about viral protein function and protein–protein interactions between viral proteins and between viral and host proteins. In this work, we identified ASFV-ASFV protein–protein interactions (PPIs) using artificial intelligence-powered protein structure prediction tools. We benchmarked our PPI identification workflow on the Vaccinia virus, a widely studied nucleocytoplasmic large DNA virus, and found that it could identify gold-standard PPIs that have been validated in vitro in a genome-wide computational screening. We applied this workflow to more than 18,000 pairwise combinations of ASFV proteins and were able to identify seventeen novel PPIs, many of which have corroborating experimental or bioinformatic evidence for their protein–protein interactions, further validating their relevance. Two protein–protein interactions, I267L and I8L, I267L__I8L, and B175L and DP79L, B175L__DP79L, are novel PPIs involving viral proteins known to modulate host immune response.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Zika virus RNA structure controls its unique neurotropism by bipartite binding to Musashi-1

Human RNA binding protein Musashi-1 (MSI1) plays a critical role in neural progenitor cells (NPCs) by binding to various host RNA transcripts. The canonical MSI1 binding site (MBS), A/GU (1-3) AG single-strand motif, is present in many RNA virus genomes, but only Zika virus (ZIKV) genome has been demonstrated to bind MSI1. Herein, we identified the AUAG motif and the AGAA tetraloop in the Xrn1-resistant RNA 2 (xrRNA2) as the canonical and non-canonical MBS, respectively, and both are crucial for ZIKV neurotropism. More importantly, the unique AGNN-type tetraloop is evolutionally conserved, and distinguishes ZIKV from other known viruses with putative MBSs. Integrated structural analysis showed that MSI1 binds to the AUAG motif and AGAA tetraloop of ZIKV in a bipartite fashion. Thus, our results not only identified an unusual viral RNA structure responsible for MSI recognition, but also revealed a role for the highly structured xrRNA in controlling viral neurotropism.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Disruption of histone acetylation homeostasis reveals multilayered chromatin regulation for transcriptional resiliency

Background Epigenetic modifications, nucleosome occupancy, and three-dimensional chromatin architecture collectively create a multi-layered, highly interactive regulatory system for controlling genomic functionality. Dysregulation of epigenetic processes leads to a plethora of abnormalities including disease states. Therapies focused on epigenetic modulation can alter gene expression to correct dysfunction, though the perpetuation of these states and the relationships among chromatin regulatory layers is not well understood. Results Here, we investigated global and local chromatin structural and functional responses after acute histone deacetylase inhibitor treatment (suberoylanilide hydroxamic acid) in lung cancer cells across time. Treatment substantially increased global histone acetylation resulting in a pervasive but not distinctive signature. The spread of acetylation did not significantly impact global chromatin accessibility, and nucleosome remodeling largely occurred at finer scales in functionally relevant genomic regions. Indeed, both H 3 K 4 trimethylation, a mark of active transcription, and gene expression changes were altered in a controlled locus-specific manner, suggesting aberrant acetylation indirectly leads to balanced and bidirectional gene expression profiles from tighter regulation of other chromatin features. HDACi treatment induced (13%) genomic rearrangement in chromatin compartmentalization and moderate weakening of topologically associating domains. Conclusions Continuous wavelet analysis of these features demonstrates that scale-dependent, locus-specific factors influence the relationship between chromatin architecture and functional output, suggesting that regulation of transcription and nucleosome remodeling is not entirely (nor linearly) dependent upon large scale compartment exchange. Structural and functional responses are most pronounced early after treatment with partial persistence of differential local chromatin features and expression later in time; this highlights the plasticity of chromatin regulation, which may have implications for the efficacy of epigenetic treatments. These results demonstrate the effectiveness of multi-layered regulation of transcription: in resilient systems, disruption of one chromatin feature does not distort the regulation of other features in supporting a transcriptional program that allows for survival.

59 BASIC BIOLOGICAL SCIENCES↗
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↗