Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “comparative 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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Phenome‐to‐genome insights for evaluating root system architecture in field studies of maize

Abstract Understanding the genetic basis of root system architecture (RSA) in crops requires innovative approaches that enable both high‐throughput and precise phenotyping in field conditions. In this study, we evaluated multiple phenotyping and analytical frameworks for quantifying RSA in mature, field‐grown maize in three field experiments. We used forward and reverse genetic approaches to evaluate >1700 maize root crowns, including a diversity panel, a biparental mapping population, and maize mutant and wild‐type alleles at two known RSA genes,DEEPER ROOTING 1(DRO1) andRootless1(Rt1). We show the utility of increasing the dimensionality of traditional two‐dimensional (2D) techniques, referred to as the “2D multi‐view” method, to improve the capture of whole root system information for mapping genetic variation influencing RSA. Comparison of univariate and multivariate genome‐wide association study (GWAS) approaches revealed that multivariate traits were effective at dissecting complex RSA phenotypes and identifying pleiotropic quantitative trait loci (QTLs). Overall, three‐dimensional (3D) root models generated from X‐ray computed tomography and digital phenotyping captured a larger proportion of RSA trait variations compared to other methods of root phenotyping, as evidenced by both genome‐wide and single‐gene analyses. Among the individual root traits, root pulling force emerged as a highly heritable estimate of RSA that identified the largest number of shared QTLs with 3D phenotypes. Our study shows that integrating complementary phenotyping technologies helps to provide a more comprehensive understanding of the genetic architecture of RSA in field‐grown maize.

Genetics & Heredity↗
DOE OSTI2024

Adaptation strategies of giant viruses to low-temperature marine ecosystems

Abstract Microbes in marine ecosystems have evolved their gene content to thrive successfully in the cold. Although this process has been reasonably well studied in bacteria and selected eukaryotes, less is known about the impact of cold environments on the genomes of viruses that infect eukaryotes. Here, we analyzed cold adaptations in giant viruses (Nucleocytoviricota and Mirusviricota) from austral marine environments and compared them with their Arctic and temperate counterparts. We recovered giant virus metagenome-assembled genomes (98 Nucleocytoviricota and 12 Mirusviricota MAGs) from 61 newly sequenced metagenomes and metaviromes from sub-Antarctic Patagonian fjords and Antarctic seawater samples. When analyzing our data set alongside Antarctic and Arctic giant viruses MAGs already deposited in the Global Ocean Eukaryotic Viral database, we found that Antarctic and Arctic giant viruses predominantly inhabit sub-10°C environments, featuring a high proportion of unique phylotypes in each ecosystem. In contrast, giant viruses in Patagonian fjords were subject to broader temperature ranges and showed a lower degree of endemicity. However, despite differences in their distribution, giant viruses inhabiting low-temperature marine ecosystems evolved genomic cold-adaptation strategies that led to changes in genetic functions and amino acid frequencies that ultimately affect both gene content and protein structure. Such changes seem to be absent in their mesophilic counterparts. The uniqueness of these cold-adapted marine giant viruses may now be threatened by climate change, leading to a potential reduction in their biodiversity.

Environmental Sciences & Ecology↗
DOE OSTI2024

Genome-resolved analysis of Serratia marcescens strain SMTT infers niche specialization as a hydrocarbon-degrader

Abstract Bacteria that are chronically exposed to high levels of pollutants demonstrate genomic and corresponding metabolic diversity that complement their strategies for adaptation to hydrocarbon-rich environments. Whole genome sequencing was carried out to infer functional traits of Serratia marcescens strain SMTT recovered from soil contaminated with crude oil. The genome size (Mb) was 5,013,981 with a total gene count of 4,842. Comparative analyses with carefully selected S. marcescens strains, 2 of which are associated with contaminated soil, show conservation of central metabolic pathways in addition to intra-specific genetic diversity and metabolic flexibility. Genome comparisons also indicated an enrichment of genes associated with multidrug resistance and efflux pumps for SMTT. The SMTT genome contained genes that enable the catabolism of aromatic compounds via the protocatechuate para-degradation pathway, in addition to meta-cleavage of catechol (meta-cleavage pathway II); gene enrichment for aromatic compound degradation was markedly higher for SMTT compared to the other S. marcescens strains analysed. Our data presents a valuable genetic inventory for future studies on strains of S. marcescens and provides insights into those genomic features of SMTT with industrial potential.

Genetics & Heredity↗
DOE OSTI2024

Extraordinary preservation of gene collinearity over three hundred million years revealed in homosporous lycophytes

Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes ( Selaginella and Isoetes ) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum . Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Comparative phylogenetics of repetitive elements in a diverse order of flowering plants (Brassicales)

Genome sizes of plants have long piqued the interest of researchers due to the vast differences among organisms. However, the mechanisms that drive size differences have yet to be fully understood. Two important contributing factors to genome size are expansions of repetitive elements, such as transposable elements (TEs), and whole-genome duplications (WGD). Although studies have found correlations between genome size and both TE abundance and polyploidy, these studies typically test for these patterns within a genus or species. The plant order Brassicales provides an excellent system to further test if genome size evolution patterns are consistent across larger time scales, as there are numerous WGDs. This order is also home to one of the smallest plant genomes, Arabidopsis thaliana—chosen as the model plant system for this reason—as well as to species with very large genomes. With new methods that allow for TE characterization from low-coverage genome shotgun data and 71 taxa across the Brassicales, we confirm the correlation between genome size and TE content, however, we are unable to reconstruct phylogenetic relationships and do not detect any shift in TE abundance associated with WGD.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A genome-informed higher rank classification of the biotechnologically important fungal subphylum Saccharomycotina

The subphylum Saccharomycotina is a lineage in the fungal phylum Ascomycota that exhibits levels of genomic diversity similar to those of plants and animals. The Saccharomycotina consist of more than 1 200 known species currently divided into 16 families, one order, and one class. Species in this subphylum are ecologically and metabolically diverse and include important opportunistic human pathogens, as well as species important in biotechnological applications. Many traits of biotechnological interest are found in closely related species and often restricted to single phylogenetic clades. However, the biotechnological potential of most yeast species remains unexplored. Although the subphylum Saccharomycotina has much higher rates of genome sequence evolution than its sister subphylum, Pezizomycotina, it contains only one class compared to the 16 classes in Pezizomycotina. The third subphylum of Ascomycota, the Taphrinomycotina, consists of six classes and has approximately 10 times fewer species than the Saccharomycotina. These data indicate that the current classification of all these yeasts into a single class and a single order is an underappreciation of their diversity. Our previous genome-scale phylogenetic analyses showed that the Saccharomycotina contains 12 major and robustly supported phylogenetic clades; seven of these are current families (Lipomycetaceae, Trigonopsidaceae, Alloascoideaceae, Pichiaceae, Phaffomycetaceae, Saccharomycodaceae, and Saccharomycetaceae), one comprises two current families (Dipodascaceae and Trichomonascaceae), one represents the genus Sporopachydermia, and three represent lineages that differ in their translation of the CUG codon (CUG-Ala, CUG-Ser1, and CUG-Ser2). Using these analyses in combination with relative evolutionary divergence and genome content analyses, we propose an updated classification for the Saccharomycotina, including seven classes and 12 orders that can be diagnosed by genome content. This updated classification is consistent with the high levels of genomic diversity within this subphylum and is necessary to make the higher rank classification of the Saccharomycotina more comparable to that of other fungi, as well as to communicate efficiently on lineages that are not yet formally named.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

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

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

Kulyashov, M. A.↗
DOE OSTI2022

Chapter 20: Advances and Application of CRISPR-Cas Systems

A new gene-editing technique called clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) proteins have revolutionized genome engineering because of its high efficiency, relatively low cost, and ease of use compared with other techniques such as zinc-finger nucleases and transcription activator-like effector nucleases. CRISPR-Cas systems have transformed biological research, quickly becoming the preferred method for engineering specific genome sequences in industrially relevant microbes, important food crops, and human cells. This powerful new tool has the potential for limitless applications ranging from the treatment of human diseases to improving food security to the generation of renewable and sustainable bioproducts and biofuels. In this chapter, we will discuss advances and applications for these CRISPR-Cas systems in biological engineering.

bacterial↗
DOE OSTI2023

Comparative feedstock analysis in Setaria viridis L. as a model for C 4 bioenergy grasses and panicoid crop species

Second generation feedstocks for bioethanol will likely include a sizable proportion of perennial C4 grasses, principally in the Panicoideae clade. The Panicoideae contain agronomically important annual grasses including Zea mays L. (maize), Sorghum bicolor (L.) Moench (sorghum), and Saccharum officinarum L. (sugar cane) as well as promising second generation perennial feedstocks including Miscanthus × giganteus and Panicum virgatum L. (switchgrass). The underlying complexity of these polyploid grass genomes is a major limitation for their direct manipulation and thus driving a need for rapidly cycling comparative model. Setaria viridis (green millet) is a rapid cycling C4 panicoid grass with a relatively small and sequenced diploid genome and abundant seed production. Stable, transient, and protoplast transformation technologies have also been developed for Setaria viridis making it a potentially excellent model for other C4 bioenergy grasses. Here, the lignocellulosic feedstock composition, cellulose biosynthesis inhibitor response and saccharification dynamics of Setaria viridis are compared with the annual sorghum and maize and the perennial switchgrass bioenergy crops as a baseline study into the applicability for translational research. A genome-wide systematic investigation of the cellulose synthase-A genes was performed identifying eight candidate sequences. Two developmental stages; (a) metabolically active young tissue and (b) metabolically plateaued (mature) material are examined to compare biomass performance metrics.

09 BIOMASS FUELS↗
DOE OSTI2022

Accelerating Scientific Workflows on HPC Platforms with In Situ Processing

Scientific workflows drive most modern large-scale science breakthroughs by allowing scientists to define their computations as a set of jobs executed in a given order based on their data dependencies. Workflow management systems (WMSs) have become key to automating scientific workflows-executing computational jobs and orchestrating data transfers between those jobs running on complex high-performance computing (HPC) platforms. Traditionally, WMSs use files to communicate between jobs: a job writes out files that are read by other jobs. However, HPC machines face a growing gap between their storage and compute capabilities. To address that concern, the scientific community has adopted a new approach called in situ, which bypasses costly parallel filesystem I/O operations with faster in-memory or in-network communications. When using in situ approaches, communication and computations can be interleaved. In this work, we leverage the Decaf in situ dataflow framework to accelerate task-based scientific workflows managed by the Pegasus WMS, by replacing file communications with faster MPI messaging. We propose a new execution engine that uses Decaf to manage communications within a sub-workflow (i.e., set of jobs) to optimize inter-job communications. We consider two workflows in this study: (i) a synthetic workflow that benchmarks and compares file- and MPI-based communication; and (ii) a realistic bioinformatics workflow that computes mu-tational overlaps in the human genome. Experiments show that in situ communication can improve the bioinformatics workflow execution time by 22% to 30% compared with file communication. Our results motivate further opportunities and challenges for bridging traditional WMSs with in situ frameworks.

Decaf↗
DOE OSTI2023

Comparison of wild-type KT2440 and genome-reduced EM42 Pseudomonas putida strains for muconate production from aromatic compounds and glucose

Pseudomonas putida KT2440 is a robust, aromatic catabolic bacterium that has been widely engineered to convert bio-based and waste-based feedstocks to target products. Towards industrial domestication of P. putida KT2440, rational genome reduction has been previously conducted, resulting in P. putida strain EM42, which exhibited characteristics that could be advantageous for production strains. Here, we compared P. putida KT2440-and EM42-derived strains for cis,cis-muconic acid production from an aromatic compound, p-coumarate, and in separate strains, from glucose. To our surprise, the EM42-derived strains did not outperform the KT2440-derived strains in muconate production from either substrate. In bioreactor cultivations, KT2440-and EM42-derived strains produced muconate from p-coumarate at titers of 45 g/L and 37 g/L, respectively, and from glucose at 20 g/L and 13 g/L, respectively. To provide additional insights about the differences in the parent strains, we analyzed growth profiles of KT2440 and EM42 on aromatic compounds as the sole carbon and energy sources. In general, the EM42 strain exhibited reduced growth rates but shorter growth lags than KT2440. We also observed that EM42-derived strains resulted in higher growth rates on glucose compared to KT2440-derived strains, but only at the lowest glucose concentrations tested. Transcriptomics revealed that genome reduction in EM42 had global effects on transcript levels and showed that the EM42-derived strains that produce muconate from glucose exhibit reduced modulation of gene expression in response to changes in glucose concentrations. Overall, our results highlight that additional studies are warranted to understand the effects of genome reduction on microbial metabolism and physiology, especially when intended for use in production strains.

09 BIOMASS FUELS↗
DOE OSTI2022

METABOLIC: high-throughput profiling of microbial genomes for functional traits, metabolism, biogeochemistry, and community-scale functional networks

Background Advances in microbiome science are being driven in large part due to our ability to study and infer microbial ecology from genomes reconstructed from mixed microbial communities using metagenomics and single-cell genomics. Such omics-based techniques allow us to read genomic blueprints of microorganisms, decipher their functional capacities and activities, and reconstruct their roles in biogeochemical processes. Currently available tools for analyses of genomic data can annotate and depict metabolic functions to some extent; however, no standardized approaches are currently available for the comprehensive characterization of metabolic predictions, metabolite exchanges, microbial interactions, and microbial contributions to biogeochemical cycling. Results We present METABOLIC (METabolic And BiogeOchemistry anaLyses In miCrobes), a scalable software to advance microbial ecology and biogeochemistry studies using genomes at the resolution of individual organisms and/or microbial communities. The genome-scale workflow includes annotation of microbial genomes, motif validation of biochemically validated conserved protein residues, metabolic pathway analyses, and calculation of contributions to individual biogeochemical transformations and cycles. The community-scale workflow supplements genome-scale analyses with determination of genome abundance in the microbiome, potential microbial metabolic handoffs and metabolite exchange, reconstruction of functional networks, and determination of microbial contributions to biogeochemical cycles. METABOLIC can take input genomes from isolates, metagenome-assembled genomes, or single-cell genomes. Results are presented in the form of tables for metabolism and a variety of visualizations including biogeochemical cycling potential, representation of sequential metabolic transformations, community-scale microbial functional networks using a newly defined metric “MW-score” (metabolic weight score), and metabolic Sankey diagrams. METABOLIC takes ~ 3 h with 40 CPU threads to process ~ 100 genomes and corresponding metagenomic reads within which the most compute-demanding part of hmmsearch takes ~ 45 min, while it takes ~ 5 h to complete hmmsearch for ~ 3600 genomes. Tests of accuracy, robustness, and consistency suggest METABOLIC provides better performance compared to other software and online servers. To highlight the utility and versatility of METABOLIC, we demonstrate its capabilities on diverse metagenomic datasets from the marine subsurface, terrestrial subsurface, meadow soil, deep sea, freshwater lakes, wastewater, and the human gut. Conclusion METABOLIC enables the consistent and reproducible study of microbial community ecology and biogeochemistry using a foundation of genome-informed microbial metabolism, and will advance the integration of uncultivated organisms into metabolic and biogeochemical models. METABOLIC is written in Perl and R and is freely available under GPLv3 at https://github.com/AnantharamanLab/METABOLIC.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Improving I/O-aware Workflow Scheduling via Data Flow Characterization and trade-off Analysis

The scientific computing paradigm has transitioned from compute-intensive to I/O-intensive and memory-intensive in the past decade, especially when data-driven science has become common practice. Numerous empirical I/O-aware scheduling optimizations have been developed by incorporating I/O capacity and bandwidth as constraints into scheduling. Unfortunately, there is a lack of data flow (I/O) characterization tool and an understanding of trade-offs between concurrency, locality, and I/O bandwidth. To bridge the gap, this work 1) presents a set of descriptors to characterize, organize, and visualize I/O profiles, including flow size, I/O bandwidth, and operation count, which group data flows by I/O types, tasks, and files; 2) proposes an I/O Roofline model-based trade-off analysis to find the optimal trade-off between flow operational intensity, concurrency, and flow performance. The I/O descriptors generate useful insights into complicated I/O behaviors, suggesting distinct concurrency, storage, and scheduling to be used by types, tasks, and files. The proposed trade-off analysis guides scheduling decisions that generate resource assignment with the best flow parallelism. We evaluate our I/O-aware scheduling methodology on a highly I/O-intensive workflow–1000 Genomes. The experimental results demonstrate speedups of up to 2.4× compared to the state-of-the- art methods.

Guo, Luanzheng [BATTELLE (PACIFIC NW LAB)]↗
DOE OSTI2025

Packaged delivery of CRISPR–Cas9 ribonucleoproteins accelerates genome editing

Effective genome editing requires a sufficient dose of CRISPR–Cas9 ribonucleoproteins (RNPs) to enter the target cell while minimizing immune responses, off-target editing, and cytotoxicity. Clinical use of Cas9 RNPs currently entails electroporation into cells ex vivo, but no systematic comparison of this method to packaged RNP delivery has been made. Here we compared two delivery strategies, electroporation and enveloped delivery vehicles (EDVs), to investigate the Cas9 dosage requirements for genome editing. Using fluorescence correlation spectroscopy, we determined that >1300 Cas9 RNPs per nucleus are typically required for productive genome editing. EDV-mediated editing was >30-fold more efficient than electroporation, and editing occurs at least 2-fold faster for EDV delivery at comparable total Cas9 RNP doses. We hypothesize that differences in efficacy between these methods result in part from the increased duration of RNP nuclear residence resulting from EDV delivery. Our results directly compare RNP delivery strategies, showing that packaged delivery could dramatically reduce the amount of CRISPR–Cas9 RNPs required for experimental or clinical genome editing.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Tracking the ancestry of known and ‘ghost’ homeologous subgenomes in model grass Brachypodium polyploids

Unraveling the evolution of plant polyploids is a challenge when their diploid progenitor species are extinct or unknown or when genome sequences of known progenitors are unavailable. Existing subgenome identification methods cannot adequately infer the homeologous genomes that are present in the allopolyploids if they do not take into account the potential existence of unknown progenitors. We addressed this challenge in the widely distributed dysploid grass genus Brachypodium, which is a model genus for temperate cereals and biofuel grasses. We used a transcriptome-based phylogeny and newly designed subgenome detection algorithms coupled with a comparative chromosome barcoding analysis. Our phylogenomic subgenome detection pipeline was validated in Triticum allopolyploids, which have known progenitor genomes, and then used to infer the identities of three subgenomes derived from extant diploid species and four subgenomes derived from unknown diploid progenitors (ghost subgenomes) in six Brachypodium polyploids (B. mexicanum, B. boissieri, B. retusum, B. phoenicoides, B. rupestre and B. hybridum), of which five contain undescribed homeologous subgenomes. The existence of the seven Brachypodium progenitor genomes in the polyploids was confirmed by their karyotypic barcode profiles. Comparative phylogenomics of nuclear versus plastid trees allowed us to formulate hypothetical homoploid hybridizations and allo- and autopolyploidization scenarios that could have generated the six Brachypodium polyploids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Engineering self-deliverable ribonucleoproteins for genome editing in the brain

The delivery of CRISPR ribonucleoproteins (RNPs) for genome editing in vitro and in vivo has important advantages over other delivery methods, including reduced off-target and immunogenic effects. However, effective delivery of RNPs remains challenging in certain cell types due to low efficiency and cell toxicity. To address these issues, we engineer self-deliverable RNPs that can promote efficient cellular uptake and carry out robust genome editing without the need for helper materials or biomolecules. Screening of cell-penetrating peptides (CPPs) fused to CRISPR-Cas9 protein identifies potent constructs capable of efficient genome editing of neural progenitor cells. Further engineering of these fusion proteins establishes a C-terminal Cas9 fusion with three copies of A22p, a peptide derived from human semaphorin-3a, that exhibits substantially improved editing efficacy compared to other constructs. We find that self-deliverable Cas9 RNPs generate robust genome edits in clinically relevant genes when injected directly into the mouse striatum. Overall, self-deliverable Cas9 proteins provide a facile and effective platform for genome editing in vitro and in vivo.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Genomic variation across Chinook salmon populations reveals effects of a duplication on migration alleles and supports fine scale structure

Abstract The distribution of ecotypic variation in natural populations is influenced by neutral and adaptive evolutionary forces that are challenging to disentangle. This study provides a high‐resolution portrait of genomic variation in Chinook salmon ( Oncorhynchus tshawytscha ) with emphasis on a region of major effect for ecotypic variation in migration timing. With a filtered data set of ~13 million single nucleotide polymorphisms (SNPs) from low‐coverage whole genome resequencing of 53 populations (3566 barcoded individuals), we contrasted patterns of genomic structure within and among major lineages and examined the extent of a selective sweep at a major effect region underlying migration timing (GREB1L/ROCK1). Neutral variation provided support for fine‐scale structure of populations, while allele frequency variation in GREB1L/ROCK1 was highly correlated with mean return timing for early and late migrating populations within each of the lineages ( r 2 = .58–.95; p < .001). However, the extent of selection within the genomic region controlling migration timing was much narrower in one lineage (interior stream‐type) compared to the other two major lineages, which corresponded to the breadth of phenotypic variation in migration timing observed among lineages. Evidence of a duplicated block within GREB1L/ROCK1 may be responsible for reduced recombination in this portion of the genome and contributes to phenotypic variation within and across lineages. Lastly, SNP positions across GREB1L/ROCK1 were assessed for their utility in discriminating migration timing among lineages, and we recommend multiple markers nearest the duplication to provide highest accuracy in conservation applications such as those that aim to protect early migrating Chinook salmon. These results highlight the need to investigate variation throughout the genome and the effects of structural variants on ecologically relevant phenotypic variation in natural species.

Horn, Rebekah L.↗
DOE OSTI2021

Impact of short-read sequencing on the misassembly of a plant genome

Abstract Background Availability of plant genome sequences has led to significant advances. However, with few exceptions, the great majority of existing genome assemblies are derived from short read sequencing technologies with highly uneven read coverages indicative of sequencing and assembly issues that could significantly impact any downstream analysis of plant genomes. In tomato for example, 0.6% (5.1 Mb) and 9.7% (79.6 Mb) of short-read based assembly had significantly higher and lower coverage compared to background, respectively. Results To understand what the causes may be for such uneven coverage, we first established machine learning models capable of predicting genomic regions with variable coverages and found that high coverage regions tend to have higher simple sequence repeat and tandem gene densities compared to background regions. To determine if the high coverage regions were misassembled, we examined a recently available tomato long-read based assembly and found that 27.8% (1.41 Mb) of high coverage regions were potentially misassembled of duplicate sequences, compared to 1.4% in background regions. In addition, using a predictive model that can distinguish correctly and incorrectly assembled high coverage regions, we found that misassembled, high coverage regions tend to be flanked by simple sequence repeats, pseudogenes, and transposon elements. Conclusions Our study provides insights on the causes of variable coverage regions and a quantitative assessment of factors contributing to plant genome misassembly when using short reads and the generality of these causes and factors should be tested further in other species.

59 BASIC BIOLOGICAL SCIENCES↗