Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional 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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

DYNC1H1 regulates NSCLC cell growth and metastasis by IFN-γ-JAK-STAT signaling and is associated with an aberrant immune response

Highlights: • DYNC1H1 is upregulated and its high expression predicted poor outcomes in NSCLC. • Loss of function of DYNC1H1 inhibits cell growth and metastasis of NSCLC cells in vitro. • DYNC1H1 knockdown inhibited metastasis of NSCLC cells in vivo. • DYNC1H1 loss of function suppresses IFN-γ-JAK-STAT signaling pathway in NSCLC. • DYNC1H1 expression is correlated with aberrant immune response in NSCLC indicated by analysis of TCGA immune data. It is urgent to identify new biomarkers and therapeutic targets to ameliorate the clinical prognosis of patients with lung cancer. The functional significance and molecular mechanism of dynein cytoplasmic 1 heavy chain 1 (DYNC1H1) in nonsmall cell lung cancer (NSCLC) progression is still elusive. In our current study, publicly available data and Western blotting experiments confirmed that DYNC1H1 expression was upregulated in lung cancer samples compared with noncancerous samples. Quantitative real-time PCR (qPCR) results indicated that high DYNC1H1 expression in lung cancer tissues was significantly associated with clinical tumor stage and distal metastasis; moreover, its high expression was negatively correlated with prognosis. Functional experiments demonstrated that DYNC1H1 loss of function caused a significant decrease in cell viability and cell proliferative ability, inhibition of the cell cycle, and promotion of both migration potential and invasion potential in vitro. Animal experiments by tail vein injection of lung cancer cells showed that DYNC1H1 knockdown significantly decreased lung cancer metastasis. Mechanistically, the results from a human protein array showed changes in the IFN-γ-JAK-STAT signaling pathway, and analysis of The Cancer Genome Atlas (TCGA) immune data demonstrated that disturbance of the immune microenvironment might be involved in the impaired growth and metastatic ability mediated by DYNC1H1 loss in NSCLC. DYNC1H1 might serve as a promising biological marker of prognosis and a potential clinical therapeutic target for patients with NSCLC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Genome-wide approaches delineate the additive, epistatic, and pleiotropic nature of variants controlling fatty acid composition in peanut ( Arachis hypogaea L.)

Abstract The fatty acid composition of seed oil is a major determinant of the flavor, shelf-life, and nutritional quality of peanuts. Major QTLs controlling high oil content, high oleic content, and low linoleic content have been characterized in several seed oil crop species. Here, we employ genome-wide association approaches on a recently genotyped collection of 787 plant introduction accessions in the USDA peanut core collection, plus selected improved cultivars, to discover markers associated with the natural variation in fatty acid composition, and to explain the genetic control of fatty acid composition in seed oils. Overall, 251 single nucleotide polymorphisms (SNPs) had significant trait associations with the measured fatty acid components. Twelve SNPs were associated with two or three different traits. Of these loci with apparent pleiotropic effects, 10 were associated with both oleic (C18:1) and linoleic acid (C18:2) content at different positions in the genome. In all 10 cases, the favorable allele had an opposite effect—increasing and lowering the concentration, respectively, of oleic and linoleic acid. The other traits with pleiotropic variant control were palmitic (C16:0), behenic (C22:0), lignoceric (C24:0), gadoleic (C20:1), total saturated, and total unsaturated fatty acid content. One hundred (100) of the significantly associated SNPs were located within 1000 kbp of 55 genes with fatty acid biosynthesis functional annotations. These genes encoded, among others: ACCase carboxyl transferase subunits, and several fatty acid synthase II enzymes. With the exception of gadoleic (C20:1) and lignoceric (C24:0) acid content, which occur at relatively low abundance in cultivated peanuts, all traits had significant SNP interactions exceeding a stringent Bonferroni threshold (α = 1%). We detected 7682 pairwise SNP interactions affecting the relative abundance of fatty acid components in the seed oil. Of these, 627 SNP pairs had at least one SNP within 1000 kbp of a gene with fatty acid biosynthesis functional annotation. We evaluated 168 candidate genes underlying these SNP interactions. Functional enrichment and protein-to-protein interactions supported significant interactions (P-value < 1.0E−16) among the genes evaluated. These results show the complex nature of the biology and genes underlying the variation in seed oil fatty acid composition and contribute to an improved genotype-to-phenotype map for fatty acid variation in peanut seed oil.

Otyama, Paul I.↗
DOE OSTI2022

The IMG/M data management and analysis system v.7: content updates and new features

The Integrated Microbial Genomes & Microbiomes system at the Department of Energy (DOE) Joint Genome Institute (JGI) continues to provide support for users to perform comparative analysis of isolate and single cell genomes, metagenomes, and metatranscriptomes. In addition to datasets produced by the JGI, IMG v.7 also includes datasets imported from public sources such as NCBI Genbank, SRA, and the DOE National Microbiome Data Collaborative (NMDC), or submitted by external users. In the past couple years, we have continued our effort to help the user community by improving the annotation pipeline, upgrading the contents with new reference database versions, and adding new analysis functionalities such as advanced scaffold search, Average Nucleotide Identity (ANI) for high-quality metagenome bins, new cassette search, improved gene neighborhood display, and improvements to metatranscriptome data display and analysis. Here, we also extended the collaboration and integration efforts with other DOE-funded projects such as NMDC and DOE Biology Knowledgebase (KBase).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

siRNA high throughput screening identifies regulators of chloropicrin and hydrogen fluoride injury in human corneal epithelial cell models

Corneal injuries induced by various toxicants result in similar clinical presentations such as corneal opacity and neovascularization. Many studies suggest that several weeks post-exposure a convergence of the molecular mechanisms drives these progressive pathologies. However, chemical agents vary in toxicological properties, and early molecular responses are anticipated to be somewhat dissimilar for different toxicants. In this report we chose 3120 targets from the Dharmacon Human Druggable genome to screen for chloropicrin (CP) and hydrogen fluoride (HF) corneal injury as we hypothesized that targets identified in vitro may be effective as therapeutic targets in future studies. Human immortalized corneal epithelial cells (SV40-HCEC) were used for screening. Cell viability and IL-8 were analyzed to down-select hits into validation studies, where multiplex cytokine analysis and high content analysis were performed to understand toxicant effect and target function. Some endpoints were also evaluated in a second human immortalized corneal epithelial cell line, TCEpi. Over 20 targets entered validation studies for CP and HF; of these, only three targets were shared: NR3C1, RELA, and KMT5A. These findings suggest that early molecular responses to different toxicants may be somewhat distinctive and present dissimilar targets for possible early intervention.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Crabtree/Warburg-like aerobic xylose fermentation by engineered Saccharomyces cerevisiae

Bottlenecks in the efficient conversion of xylose into cost-effective biofuels have limited the widespread use of plant lignocellulose as a renewable feedstock. The yeast Saccharomyces cerevisiae ferments glucose into ethanol with such high metabolic flux that it ferments high concentrations of glucose aerobically, a trait called the Crabtree/Warburg Effect. In contrast to glucose, most engineered S. cerevisiae strains do not ferment xylose at economically viable rates and yields, and they require respiration to achieve sufficient xylose metabolic flux and energy return for growth aerobically. Here, we evolved respiration-deficient S. cerevisiae strains that can grow on and ferment xylose to ethanol aerobically, a trait analogous to the Crabtree/Warburg Effect for glucose. Through genome sequence comparisons and directed engineering, we determined that duplications of genes encoding engineered xylose metabolism enzymes, as well as TKL1, a gene encoding a transketolase in the pentose phosphate pathway, were the causative genetic changes for the evolved phenotype. Reengineered duplications of these enzymes, in combination with deletion mutations in HOG1, ISU1, GRE3, and IRA2, increased the rates of aerobic and anaerobic xylose fermentation. Importantly, we found that these genetic modifications function in another genetic background and increase the rate and yield of xylose-to-ethanol conversion in industrially relevant switchgrass hydrolysate, indicating that these specific genetic modifications may enable the sustainable production of industrial biofuels from yeast. We propose a model for how key regulatory mutations prime yeast for aerobic xylose fermentation by lowering the threshold for overflow metabolism, allowing mutations to increase xylose flux and to redirect it into fermentation products.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Abstract Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN)cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genesMAFB, PHOX2A, CHN1, andEBF3– as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Science & Technology - Other Topics↗
DOE OSTI2025

Orthogonal replication with optogenetic selection evolves yeast JEN1 into a mevalonate transporter

Abstract The in vivo continuous evolution system OrthoRep (orthogonal replication) is a powerful strategy for rapid enzyme evolution inSaccharomyces cerevisiaethat diversifies genes at a rate exceeding the endogenous genome mutagenesis rate by several orders of magnitude. However, it is difficult to neofunctionalize genes using OrthoRep partly because of the way selection pressures are applied. Here we combine OrthoRep with optogenetics in a selection strategy we call OptoRep, which allows fine-tuning of selection pressure with light. With this capability, we evolved a truncated form of the endogenous monocarboxylate transporterJEN1 (JEN1t)into a de novo mevalonate importer. We demonstrate the functionality of the evolvedJEN1t(JEN1t Y180C/G ) in the production of farnesene, a renewable aviation biofuel, from mevalonate fed to fermentation media or produced by microbial consortia. This study shows that the light-induced complementation of OptoRep may improve the ability to evolve functions not currently accessible for selection, while its fine tunability of selection pressure may allow the continuous evolution of genes whose desired function has a restrictive range between providing effective selection and cellular viability.

Biochemistry & Molecular Biology↗
DOE OSTI2021

The Human Proteoform Project: Defining the human proteome

Proteins are the primary effectors of function in biology, and thus, complete knowledge of their structure and properties is fundamental to deciphering function in basic and translational research. The chemical diversity of proteins is expressed in their many proteoforms, which result from combinations of genetic polymorphisms, RNA splice variants, and posttranslational modifications. This knowledge is foundational for the biological complexes and networks that control biology yet remains largely unknown. We propose here an ambitious initiative to define the human proteome, that is, to generate a definitive reference set of the proteoforms produced from the genome. Several examples of the power and importance of proteoform-level knowledge in disease-based research are presented along with a call for improved technologies in a two-pronged strategy to the Human Proteoform Project.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRSDate not supplied

A Nonlinear Viscoelastic-Viscoplastic Constitutive Model for Epoxy Polymers

The objective of this paper is to develop a nonlinear viscoelastic–viscoplastic constitutive model for epoxy polymers. The classic nonlinear viscoelasticity model is reformulated to yield a closed-form incremental constitutive relation, which relates the stress increments to the viscoelastic strain increments. A viscoplasticity model, which consists of the Drucker–Prager yield function, nonlinear isotropic and kinematic hardening laws, and the Perzyna viscosity function, is subsequently developed. The nonlinear viscoelasticity and the vis-coplasticity models are then implemented in a radial return algorithm. A closed-form algorithmic tangent operator is derived to facilitate model implementation. The nonlinear viscoelasticity model is validated by reproducing the test data for polymethyl methacrylate(PMMA). The present constitutive model’s capabilities are demonstrated through modeling viscoelastic–viscoplastic PMMA loaded at different strain rates. The mechanics of structure genome-based micromechanics approach, embedded with the present constitutive model, is used to homogenize a unidirectional fiber-reinforced composite with a PMMA matrix, subjected to uniaxial and shear loading at different strain rates. The viscoelastic characterization of PMT-F7 epoxy is evaluated by comparing the predicted creep–recovery responses by the present constitutive model with the experimental ones.

Liang Zhang↗
DOE OSTI2024

High-resolution RNA tertiary structures in Zika virus stem-loop A for the development of inhibitory small molecules

Flaviviruses such as Zika (ZIKV) and dengue virus (DENV) are positive-sense RNA viruses belonging to Flaviviridae. The flavivirus genome contains a 5’ end stem-loop promoter sequence known as stem-loop A (SLA) that is recognized by the flavivirus polymerase NS5 during viral RNA synthesis and 5’ guanosine cap methylation. The crystal structures of ZIKV and DENV SLAs show a well-defined fold, consisting of a bottom stem, side loop, and top stem-loop, providing unique interaction sites for small-molecule inhibitors to disrupt the promoter function. To facilitate identification of small molecule binding sites in flavivirus SLA, we determined high-resolution structures of bottom and top stems of ZIKV SLA, which contain a single U- or G-bulge, respectively. Both bulge nucleotides exhibit multiple orientations, from folded back on the adjacent nucleotide to flipped out of the helix, and are stabilized by stacking or base triple interactions. These structures suggest that even a single unpaired nucleotide can provide flexibility to RNA structures and its conformation is mainly determined by the stabilizing chemical environment. To facilitate discovery of small molecule inhibitors that interfere with the functions of ZIKV SLA, we screened and identified compounds that bind to the bottom and top stems of ZIKV SLA.

Biochemistry & Molecular Biology↗
DOE OSTI2020

Clades of huge phages from across Earth’s ecosystems

Bacteriophages typically have small genomes and depend on their bacterial hosts for replication. Here we sequenced DNA from diverse ecosystems and found hundreds of phage genomes with lengths of more than 200 kilobases (kb), including a genome of 735 kb, which is—to our knowledge—the largest phage genome to be described to date. Thirty-five genomes were manually curated to completion (circular and no gaps). Expanded genetic repertoires include diverse and previously undescribed CRISPR–Cas systems, transfer RNAs (tRNAs), tRNA synthetases, tRNA-modification enzymes, translation-initiation and elongation factors, and ribosomal proteins. The CRISPR–Cas systems of phages have the capacity to silence host transcription factors and translational genes, potentially as part of a larger interaction network that intercepts translation to redirect biosynthesis to phage-encoded functions. In addition, some phages may repurpose bacterial CRISPR–Cas systems to eliminate competing phages. We phylogenetically define the major clades of huge phages from human and other animal microbiomes, as well as from oceans, lakes, sediments, soils and the built environment. We conclude that the large gene inventories of huge phages reflect a conserved biological strategy, and that the phages are distributed across a broad bacterial host range and across Earth’s ecosystems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A minimal complex of KHNYN and zinc-finger antiviral protein binds and degrades single-stranded RNA

Detecting viral infection is a key role of the innate immune system. The genomes of some RNA viruses have a high CpG dinucleotide content relative to most vertebrate cell RNAs, making CpGs a molecular marker of infection. The human zinc-finger antiviral protein (ZAP) recognizes CpG, mediates clearance of the foreign CpG-rich RNA, and causes attenuation of CpG-rich RNA viruses. While ZAP binds RNA, it lacks enzymatic activity that might be responsible for RNA degradation and thus requires interacting cofactors for its function. One of these cofactors, KHNYN, has a predicted nuclease domain. Using biochemical approaches, we found that the KHNYN NYN domain is a single-stranded RNA ribonuclease that does not have sequence specificity and digests RNA with or without CpG dinucleotides equivalently in vitro. We show that unlike most KH domains, the KHNYN KH domain does not bind RNA. Indeed, a crystal structure of the KH region revealed a double-KH domain with a negatively charged surface that accounts for the lack of RNA binding. Rather, the KHNYN C-terminal domain (CTD) interacts with the ZAP RNA-binding domain (RBD) to provide target RNA specificity. We define a minimal complex composed of the ZAP RBD and the KHNYN NYN-CTD and use a fluorescence polarization assay to propose a model for how this complex interacts with a CpG dinucleotide-containing RNA. In the context of the cell, this module would represent the minimum ZAP and KHNYN domains required for CpG-recognition and ribonuclease activity essential for attenuation of viruses with clusters of CpG dinucleotides.

Yeoh, Zoe C. (ORCID:0000000226949068)↗
DOE OSTI2022

Chronic Cardiotoxicity Assays Using Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes (hiPSC-CMs)

Cardiomyocytes (CMs) differentiated from human induced pluripotent stem cells (hiPSCs) are increasingly used in cardiac safety assessment, disease modeling and regenerative medicine. A vast majority of cardiotoxicity studies in the past have tested acute effects of compounds and drugs; however, these studies lack information on the morphological or physiological responses that may occur after prolonged exposure to a cardiotoxic compound. In this review, we focus on recent advances in chronic cardiotoxicity assays using hiPSC-CMs. We summarize recently published literature on hiPSC-CMs assays applied to chronic cardiotoxicity induced by anticancer agents, as well as non-cancer classes of drugs, including antibiotics, anti-hepatitis C virus (HCV) and antidiabetic drugs. We then review publications on the implementation of hiPSC-CMs-based assays to investigate the effects of non-pharmaceutical cardiotoxicants, such as environmental chemicals or chronic alcohol consumption. We also highlight studies demonstrating the chronic effects of smoking and implementation of hiPSC-CMs to perform genomic screens and metabolomics-based biomarker assay development. The acceptance and wide implementation of hiPSC-CMs-based assays for chronic cardiotoxicity assessment will require multi-site standardization of assay protocols, chronic cardiac maturity marker reproducibility, time points optimization, minimal cellular variation (commercial vs. lab reprogrammed), stringent and matched controls and close clinical setting resemblance. A comprehensive investigation of long-term repeated exposure-induced effects on both the structure and function of cardiomyocytes can provide mechanistic insights and recapitulate drug and environmental cardiotoxicity.

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

Habitat specialization and edge effects of soil microbial communities in a fragmented landscape

Abstract Soil microorganisms play outsized roles in nutrient cycling, plant health, and climate regulation. Despite their importance, we have a limited understanding of how soil microbes are affected by habitat fragmentation, including their responses to conditions at fragment edges, or “edge effects.” To understand the responses of soil communities to edge effects, we analyzed the distributions of soil bacteria, archaea, and fungi in an experimentally fragmented system of open patches embedded within a forest matrix. In addition, we identified taxa that consistently differed among patch, edge, or matrix habitats (“specialists”) and taxa that showed no habitat preference (“nonspecialists”). We hypothesized that microbial community turnover would be most pronounced at the edge between habitats. We also hypothesized that specialist fungi would be more likely to be mycorrhizal than nonspecialist fungi because mycorrhizae should be affected more by different plant hosts among habitats, whereas specialist prokaryotes would have smaller genomes (indicating reduced metabolic versatility) and be less likely to be able to sporulate than nonspecialist prokaryotes. Across all replicate sites, the matrix and patch soils harbored distinct microbial communities. However, sites where the contrasts in vegetation and pH between the patch and matrix were most pronounced exhibited larger differences between patch and matrix communities and tended to have edge communities that differed from those in the patch and forest. There were similar numbers of patch and matrix specialists, but very few edge specialist taxa. Acidobacteria and ectomycorrhizae were more likely to be forest specialists, while Chloroflexi, Ascomycota, and Glomeromycota (i.e., arbuscular mycorrhizae) were more likely to be patch specialists. Contrary to our hypotheses, nonspecialist bacteria were not more likely than specialist bacteria to have larger genomes or to be spore‐formers. We found partial support for our mycorrhizal hypothesis: arbuscular mycorrhizae, but not ectomycorrhizae, were more likely to be specialists. Overall, our results indicate that soil microbial communities are sensitive to edges, but not all taxa are equally affected, with arbuscular mycorrhizae in particular showing a strong response to habitat edges. In the context of increasing habitat fragmentation worldwide, our results can help inform efforts to maintain the structure and functioning of the soil microbiome.

Winfrey, Claire C. [Department of Ecology and Evol↗
DOE OSTI2020

Quality-Control Mechanism for Telomerase RNA Folding in the Cell

Long non-coding RNAs can often fold into different conformations. Telomerase RNA, an essential component of the telomerase ribonucleoprotein (RNP) enzyme, must fold into a defined structure to fulfill its function with the protein catalytic subunit (TERT) and other accessory factors. However, the mechanism by which the correct folding of telomerase RNA is warranted in a cell is still unknown. Here we show that La-related protein Pof8 specifically recognizes the conserved pseudoknot region of telomerase RNA and instructs the binding of the Lsm2-8 complex to its mature 3' end, thus selectively protecting the correctly folded RNA from exonucleolytic degradation. In the absence of Pof8, TERT assembles with misfolded RNA and produces little telomerase activity. Therefore, Pof8 plays a key role in telomerase RNA folding quality control, ensuring that TERT only assembles with functional telomerase RNA to form active telomerase. Our finding reveals a mechanism for non-coding RNA folding quality control.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Population ecology and biogeochemical implications of ssDNA and dsDNA viruses along a permafrost thaw gradient

Anthropogenic-driven climate change is accelerating permafrost thaw, threatening to release vast carbon stores through increased microbial activity. While microbial roles are increasingly studied, the contributions of viruses remain largely unexplored, in part due to soil-associated technical challenges that have hindered their detection and characterization. Here, we applied an optimized virion enrichment workflow along a permafrost thaw gradient, identifying 9,963 viral populations (vOTUs), including single- and double-stranded DNA viruses, with 99.9% novelty compared to other soils. Hosts were predicted for 38% of vOTUs, spanning nine archaeal, and 36 bacterial phyla, 22% of which were linked to metagenome-assembled genomes, including key carbon-cycling taxa. Genomic analyses revealed 811 putative auxiliary metabolic genes (AMGs) from 658 vOTUs, nearly half involved in carbon processing. These included 59 glycoside hydrolases (GH) across nine GH families, 45 for monosaccharide degradation, and seven involved in short-chain fatty acid and C1 metabolism, linking viruses to both early and late stages of carbon turnover. Additionally, six vOTUs carried racD, which may stabilize microbial necromass and promote long-term carbon storage. Viral and AMG functional diversity increased with thaw stage, indicating that viruses might participate in a broadening range of microbial metabolic processes as permafrost thaws. These findings expand our understanding of virus contributions in microbial carbon processing and suggest their important role in deciphering soil carbon fate under changing climate conditions.

Biological and medical sciences↗
DOE OSTI2021

Atypical Divergence of SARS-CoV-2 Orf8 from Orf7a within the Coronavirus Lineage Suggests Potential Stealthy Viral Strategies in Immune Evasion

Orf8, one of the most puzzling genes in the SARS lineage of coronaviruses, marks a unique and striking difference in genome organization between SARS-CoV-2 and SARS-CoV-1. Here, using sequence comparisons, we unequivocally reveal the distant sequence similarities between SARS-CoV-2 Orf8 with its SARS-CoV-1 counterparts and the X4-like genes of coronaviruses, including its highly divergent “paralog” gene Orf7a, whose product is a potential immune antagonist of known structure. Supervised sequence space walks unravel identity levels that drop below 10% and yet exhibit subtle conservation patterns in this novel superfamily, characterized by an immunoglobulin-like beta sandwich topology. We document the high accuracy of the sequence space walk process in detail and characterize the subgroups of the superfamily in sequence space by systematic annotation of gene and taxon groups. While SARS-CoV-1 Orf7a and Orf8 genes are most similar to bat virus sequences, their SARS-CoV-2 counterparts are closer to pangolin virus homologs, reflecting the fine structure of conservation patterns within the SARS-CoV-2 genomes. The divergence between Orf7a and Orf8 is exceptionally idiosyncratic, since Orf7a is more constrained, whereas Orf8 is subject to rampant change, a peculiar feature that may be related to hitherto-unknown viral infection strategies. Despite their common origin, the Orf7a and Orf8 protein families exhibit different modes of evolutionary trajectories within the coronavirus lineage, which might be partly attributable to their complex interactions with the mammalian host cell, reflected by a multitude of functional associations of Orf8 in SARS-CoV-2 compared to a very small number of interactions discovered for Orf7a.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Rapid Cis–Trans Coevolution Driven by a Novel Gene Retroposed from a Eukaryotic Conserved CCR4–NOT Component in Drosophila

Young, or newly evolved, genes arise ubiquitously across the tree of life, and they can rapidly acquire novel functions that influence a diverse array of biological processes. Previous work identified a young regulatory duplicate gene in Drosophila, Zeus that unexpectedly diverged rapidly from its parent, Caf40, an extremely conserved component in the CCR4–NOT machinery in post-transcriptional and post-translational regulation of eukaryotic cells, and took on roles in the male reproductive system. This neofunctionalization was accompanied by differential binding of the Zeus protein to loci throughout the Drosophila melanogaster genome. However, the way in which new DNA-binding proteins acquire and coevolve with their targets in the genome is not understood. Here, by comparing Zeus ChIP-Seq data from D. melanogaster and D. simulans to the ancestral Caf40 binding events from D. yakuba, a species that diverged before the duplication event, we found a dynamic pattern in which Zeus binding rapidly coevolved with a previously unknown DNA motif, which we term Caf40 and Zeus-Associated Motif (CAZAM), under the influence of positive selection. Interestingly, while both copies of Zeus acquired targets at male-biased and testis-specific genes, D. melanogaster and D. simulans proteins have specialized binding on different chromosomes, a pattern echoed in the evolution of the associated motif. Using CRISPR-Cas9-mediated gene knockout of Zeus and RNA-Seq, we found that Zeus regulated the expression of 661 differentially expressed genes (DEGs). Our results suggest that the evolution of young regulatory genes can be coupled to substantial rewiring of the transcriptional networks into which they integrate, even over short evolutionary timescales. Our results thus uncover dynamic genome-wide evolutionary processes associated with new genes.

59 BASIC BIOLOGICAL SCIENCES↗