Search NASA⌕ Search

SEARCH · Search NASA

Results for “single cell 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 127 records · Page 7

Synergistic epistasis enhances the co-operativity of mutualistic interspecies interactions

Abstract Early evolution of mutualism is characterized by big and predictable adaptive changes, including the specialization of interacting partners, such as through deleterious mutations in genes not required for metabolic cross-feeding. We sought to investigate whether these early mutations improve cooperativity by manifesting in synergistic epistasis between genomes of the mutually interacting species. Specifically, we have characterized evolutionary trajectories of syntrophic interactions of Desulfovibrio vulgaris (Dv) with Methanococcus maripaludis (Mm) by longitudinally monitoring mutations accumulated over 1000 generations of nine independently evolved communities with analysis of the genotypic structure of one community down to the single-cell level. We discovered extensive parallelism across communities despite considerable variance in their evolutionary trajectories and the perseverance within many evolution lines of a rare lineage of Dv that retained sulfate-respiration (SR+) capability, which is not required for metabolic cross-feeding. An in-depth investigation revealed that synergistic epistasis across pairings of Dv and Mm genotypes had enhanced cooperativity within SR− and SR+ assemblages, enabling their coexistence within the same community. Thus, our findings demonstrate that cooperativity of a mutualism can improve through synergistic epistasis between genomes of the interacting species, enabling the coexistence of mutualistic assemblages of generalists and their specialized variants.

59 BASIC BIOLOGICAL SCIENCES↗

Mapping of flumioxazin tolerance in a snap bean diversity panel leads to the discovery of a master genomic region controlling multiple stress resistance genes

Effective weed management tools are crucial for maintaining the profitable production of snap bean (Phaseolus vulgaris L.). Preemergence herbicides help the crop to gain a size advantage over the weeds, but the few preemergence herbicides registered in snap bean have poor waterhemp (Amaranthus tuberculatus) control, a major pest in snap bean production. Waterhemp and other difficult-to-control weeds can be managed by flumioxazin, an herbicide that inhibits protoporphyrinogen oxidase (PPO). However, there is limited knowledge about crop tolerance to this herbicide. We aimed to quantify the degree of snap bean tolerance to flumioxazin and explore the underlying mechanisms. We investigated the genetic basis of herbicide tolerance using genome-wide association mapping approach utilizing field-collected data from a snap bean diversity panel, combined with gene expression data of cultivars with contrasting response. The response to a preemergence application of flumioxazin was measured by assessing plant population density and shoot biomass variables. Snap bean tolerance to flumioxazin is associated with a single genomic location in chromosome 02. Tolerance is influenced by several factors, including those that are indirectly affected by seed size/weight and those that directly impact the herbicide's metabolism and protect the cell from reactive oxygen species-induced damage. Transcriptional profiling and co-expression network analysis identified biological pathways likely involved in flumioxazin tolerance, including oxidoreductase processes and programmed cell death. Transcriptional regulation of genes involved in those processes is possibly orchestrated by a transcription factor located in the region identified in the GWAS analysis. Several entries belonging to the Romano class, including Bush Romano 350, Roma II, and Romano Purpiat presented high levels of tolerance in this study. The alleles identified in the diversity panel that condition snap bean tolerance to flumioxazin shed light on a novel mechanism of herbicide tolerance and can be used in crop improvement.

60 APPLIED LIFE SCIENCES↗

Radiation induces genomic instability and mammary ductal dysplasia in Atm heterozygous mice

Ataxia-telangiectasia (AT) is a genetic syndrome resulting from the inheritance of two defective copies of the ATM gene that includes among its stigmata radiosensitivity and cancer susceptibility. Epidemiological studies have demonstrated that although women with a single defective copy of ATM (AT heterozygotes) appear clinically normal, they may never the less have an increased relative risk of developing breast cancer. Whether they are at increased risk for radiation-induced breast cancer from medical exposures to ionizing radiation is unknown. We have used a murine model of AT to investigate the effect of a single defective Atm allele, the murine homologue of ATM, on the susceptibility of mammary epithelial cells to radiation-induced transformation. Here we report that mammary epithelial cells from irradiated mice with one copy of Atm truncated in the PI-3 kinase domain were susceptible to radiation-induced genomic instability and generated a 10% incidence of dysplastic mammary ducts when transplanted into syngenic recipients, whereas cells from Atm(+/+) mice were stable and formed only normal ducts. Since radiation-induced ductal dysplasia is a precursor to mammary cancer, the results indicate that AT heterozygosity increases susceptibility to radiogenic breast cancer in this murine model system.

Non-NASA Center↗

Construction and characterization of a genome-scale ordered mutant collection of Bacteroides thetaiotaomicron

Ordered transposon-insertion collections, in which specific transposon-insertion mutants are stored as monocultures in a genome-scale collection, represent a promising tool for genetic dissection of human gut microbiota members. However, publicly available collections are scarce and the construction methodology remains in early stages of development. Here, we describe the assembly of a genome-scale ordered collection of transposon-insertion mutants in the model gut anaerobe Bacteroides thetaiotaomicron VPI-5482 that we created as a resource for the research community. We used flow cytometry to sort single cells from a pooled library, located mutants within this initial progenitor collection by applying a pooling strategy with barcode sequencing, and re-arrayed specific mutants to create a condensed collection with single-insertion strains covering >2500 genes. To demonstrate the potential of the condensed collection for phenotypic screening, we analyzed growth dynamics and cell morphology. We identified both growth defects and altered cell shape in mutants disrupting sphingolipid synthesis and thiamine scavenging. Finally, we analyzed the process of assembling the B. theta condensed collection to identify inefficiencies that limited coverage. We demonstrate as part of this analysis that the process of assembling an ordered collection can be accurately modeled using barcode sequencing data. We expect that utilization of this ordered collection will accelerate research into B. theta physiology and that lessons learned while assembling the collection will inform future efforts to assemble ordered mutant collections for an increasing number of gut microbiota members.

59 BASIC BIOLOGICAL SCIENCES↗

Complete genome of Polaromonas vacuolata KCTC 22033 T isolated from beneath Antarctic Sea ice

Polaromonas vacuolata KCTC 22033 T is an obligate aerobic, Gram-negative, psychrophilic and rod-shaped bacterium isolated from beneath the sea ice off the coast of the Palmer Peninsula, Anvers Islands, Antarctica. P. vacuolata is the type species of Polaromonas genus and the first example of gas vacuolate Betaproteobacteria isolated from marine habitats. In this study, we report a complete genome of P. vacuolata KCTC 22033 T , which consists of 3,837,686 bp (G + C content of 52.07%) with a single chromosome, 3461 protein-coding genes, 56 tRNAs and 6 rRNA operons. Genomic analysis revealed the presence of genes involved in bacterial adaptation under saline conditions, cold adaptation via the production of gas vesicles and cell adhesion proteins, and a photoheterotrophic lifestyle when challenged by starvation. Intriguingly, several of these genes were likely acquired from species outside the Polaromonas genus. The genomic information therefore describes the unique evolution and adaptation of P. vacuolata to its extraordinary habitat, i.e., beneath the Antarctic sea ice.

59 BASIC BIOLOGICAL SCIENCES↗

Discovery of FoTO1 and Taxol genes enables biosynthesis of baccatin III

Abstract Plants make complex and potent therapeutic molecules 1,2 , but sourcing these molecules from natural producers or through chemical synthesis is difficult, which limits their use in the clinic. A prominent example is the anti-cancer therapeutic paclitaxel (sold under the brand name Taxol), which is derived from yew trees (Taxusspecies) 3 . Identifying the full paclitaxel biosynthetic pathway would enable heterologous production of the drug, but this has yet to be achieved despite half a century of research 4 . WithinTaxus’ large, enzyme-rich genome 5 , we suspected that the paclitaxel pathway would be difficult to resolve using conventional RNA-sequencing and co-expression analyses. Here, to improve the resolution of transcriptional analysis for pathway identification, we developed a strategy we term multiplexed perturbation × single nuclei (mpXsn) to transcriptionally profile cell states spanning tissues, cell types, developmental stages and elicitation conditions. Our data show that paclitaxel biosynthetic genes segregate into distinct expression modules that suggest consecutive subpathways. These modules resolved seven new genes, allowing a de novo 17-gene biosynthesis and isolation of baccatin III, the industrial precursor to Taxol, inNicotiana benthamianaleaves, at levels comparable with the natural abundance inTaxusneedles. Notably, we found that a nuclear transport factor 2 (NTF2)-like protein, FoTO1, is crucial for promoting the formation of the desired product during the first oxidation, resolving a long-standing bottleneck in paclitaxel pathway reconstitution. Together with a new β-phenylalanine-CoA ligase, the eight genes discovered here enable the de novo biosynthesis of 3’-N-debenzoyl-2’-deoxypaclitaxel. More broadly, we establish a generalizable approach to efficiently scale the power of co-expression analysis to match the complexity of large, uncharacterized genomes, facilitating the discovery of high-value gene sets.

Science & Technology - Other Topics↗

Conservation of endo -glucanase 16 (EG16) activity across highly divergent plant lineages

Plant cell walls are highly dynamic structures that are composed predominately of polysaccharides. As such, endogenous carbohydrate active enzymes (CAZymes) are central to the synthesis and subsequent modification of plant cells during morphogenesis. The endo-glucanase 16 (EG16) members constitute a distinct group of plant CAZymes, angiosperm orthologs of which were recently shown to have dual β-glucan/xyloglucan hydrolase activity. Molecular phylogeny indicates that EG16 members comprise a sister clade with a deep evolutionary relationship to the widely studied apoplastic xyloglucan endo-transglycosylases/hydrolases (XTH). A cross-genome survey indicated that EG16 members occur as a single ortholog across species and are widespread in early diverging plants, including the non-vascular bryophytes, for which functional data were previously lacking. Remarkably, enzymological characterization of an EG16 ortholog from the model moss Physcomitrella patens (PpEG16) revealed that EG16 activity and sequence/structure are highly conserved across 500 million years of plant evolution, vis-à-vis orthologs from grapevine and poplar. Ex vivo biomechanical assays demonstrated that the application of EG16 gene products caused abrupt breakage of etiolated hypocotyls rather than slow extension, thereby indicating a mode-of-action distinct from endogenous expansins and microbial endo-glucanases. As a result, the biochemical data presented here will inform future genomic, genetic, and physiological studies of EG16 enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

Spatially structured bacterial interactions alter algal carbon flow to bacteria

Phytoplankton account for nearly half of global photosynthetic carbon fixation, and the fate of that carbon is regulated in large part by microbial food web processing. We currently lack a mechanistic understanding of how interactions among heterotrophic bacteria impact the fate of photosynthetically fixed carbon. Here, we used a set of bacterial isolates capable of growing on exudates from the diatom Phaeodactylum tricornutum to investigate how bacteria-bacteria interactions affect the balance between exudate remineralization and incorporation into biomass. With exometabolomics and genome-scale metabolic modeling, we estimated the degree of resource competition between bacterial pairs. In a sequential spent media experiment, we found that pairwise interactions were more beneficial than predicted based on resource competition alone, and 30% exhibited facilitative interactions. To link this to carbon fate, we used single-cell isotope tracing in a custom cultivation system to compare the impact of different "primary" bacterial strains in close proximity to live P. tricornutum on a distal "secondary" strain. We found that a primary strain with a high degree of competition decreased secondary strain carbon drawdown by 51% at the single-cell level, providing a quantitative metric for the "cost" of competition on algal carbon fate. Additionally, a primary strain classified as facilitative based on sequential interactions increased total algal-derived carbon assimilation by 7.6 times, integrated over all members, compared to the competitive primary strain. Our findings suggest that the degree of interaction between bacteria along a spectrum from competitive to facilitative is directly linked to algal carbon drawdown.

genome-scale metabolic model↗

Gene conversion is strongly induced in human cells by double-strand breaks and is modulated by the expression of BCL-x(L)

Homology-directed repair (HDR) of DNA double-strand breaks (DSBs) contributes to the maintenance of genomic stability in rodent cells, and it has been assumed that HDR is of similar importance in DSB repair in human cells. However, some outcomes of homologous recombination can be deleterious, suggesting that factors exist to regulate HDR. We demonstrated previously that overexpression of BCL-2 or BCL-x(L) enhanced the frequency of X-ray-induced TK1 mutations, including loss of heterozygosity events presumed to arise by mitotic recombination. The present study was designed to test whether HDR is a prominent DSB repair pathway in human cells and to determine whether ectopic expression of BCL-x(L) affects HDR. Using TK6-neo cells, we find that a single DSB in an integrated HDR reporter stimulates gene conversion 40-50-fold, demonstrating efficient DSB repair by gene conversion in human cells. Significantly, DSB-induced gene conversion events are 3-4-fold more frequent in TK6 cells that stably overexpress the antiapoptotic protein BCL-X(L). Thus, HDR plays an important role in maintaining genomic integrity in human cells, and ectopic expression of BCL-x(L) enhances HDR of DSBs. This is the first study to highlight a function for BCL-x(L) in modulating DSB repair in human cells.

NASA Discipline Radiation Health↗

Proposed minimal standards for description of methanogenic archaea

Methanogenic archaea are a diverse, polyphyletic group of strictly anaerobic prokaryotes capable of producing methane as their primary metabolic product. It has been over three decades since minimal standards for their taxonomic description have been proposed. In light of advancements in technology and amendments in systematic microbiology, revision of the older criteria for taxonomic description is essential. Most of the previously recommended minimum standards regarding phenotypic characterization of pure cultures are maintained. Electron microscopy and chemotaxonomic methods like whole-cell protein and lipid analysis are desirable but not required. Because of advancements in DNA sequencing technologies, obtaining a complete or draft whole genome sequence for type strains and its deposition in a public database are now mandatory. Genomic data should be used for rigorous comparison to close relatives using overall genome related indices such as average nucleotide identity and digital DNA–DNA hybridization. Phylogenetic analysis of the 16S rRNA gene is also required and can be supplemented by phylogenies of the mcrA gene and phylogenomic analysis using multiple conserved, single-copy marker genes. Additionally, it is now established that culture purity is not essential for studying prokaryotes, and description of Candidatus methanogenic taxa using single-cell or metagenomics along with other appropriate criteria is a viable alternative. The revisions to the minimal criteria proposed here by the members of the Subcommittee on the Taxonomy of Methanogenic Archaea of the International Committee on Systematics of Prokaryotes should allow for rigorous yet practical taxonomic description of these important and diverse microbes.

Microbiology↗

An in silico assessment of gene function and organization of the phenylpropanoid pathway metabolic networks in Arabidopsis thaliana and limitations thereof

The Arabidopsis genome sequencing in 2000 gave to science the first blueprint of a vascular plant. Its successful completion also prompted the US National Science Foundation to launch the Arabidopsis 2010 initiative, the goal of which is to identify the function of each gene by 2010. In this study, an exhaustive analysis of The Institute for Genomic Research (TIGR) and The Arabidopsis Information Resource (TAIR) databases, together with all currently compiled EST sequence data, was carried out in order to determine to what extent the various metabolic networks from phenylalanine ammonia lyase (PAL) to the monolignols were organized and/or could be predicted. In these databases, there are some 65 genes which have been annotated as encoding putative enzymatic steps in monolignol biosynthesis, although many of them have only very low homology to monolignol pathway genes of known function in other plant systems. Our detailed analysis revealed that presently only 13 genes (two PALs, a cinnamate-4-hydroxylase, a p-coumarate-3-hydroxylase, a ferulate-5-hydroxylase, three 4-coumarate-CoA ligases, a cinnamic acid O-methyl transferase, two cinnamoyl-CoA reductases) and two cinnamyl alcohol dehydrogenases can be classified as having a bona fide (definitive) function; the remaining 52 genes currently have undetermined physiological roles. The EST database entries for this particular set of genes also provided little new insight into how the monolignol pathway was organized in the different tissues and organs, this being perhaps a consequence of both limitations in how tissue samples were collected and in the incomplete nature of the EST collections. This analysis thus underscores the fact that even with genomic sequencing, presumed to provide the entire suite of putative genes in the monolignol-forming pathway, a very large effort needs to be conducted to establish actual catalytic roles (including enzyme versatility), as well as the physiological function(s) for each member of the (multi)gene families present and the metabolic networks that are operative. Additionally, one key to identifying physiological functions for many of these (and other) unknown genes, and their corresponding metabolic networks, awaits the development of technologies to comprehensively study molecular processes at the single cell level in particular tissues and organs, in order to establish the actual metabolic context.

NASA Program Fundamental Space Biology↗

A maize semi-dwarf mutant reveals a GRAS transcription factor involved in brassinosteroid signaling

Brassinosteroids (BR) and gibberellins (GA) regulate plant height and leaf angle in maize (Zea mays). Mutants with defects in BR or GA biosynthesis or signaling identify components of these pathways and enhance our knowledge about plant growth and development. In this study, we characterized three recessive mutant alleles of GRAS transcription factor 42 (gras42) in maize, a GRAS transcription factor gene orthologous to the DWARF AND LOW TILLERING (DLT) gene of rice (Oryza sativa). These maize mutants exhibited semi-dwarf stature, shorter and wider leaves, and more upright leaf angle. Transcriptome analysis revealed a role for GRAS42 as a determinant of BR signaling. Analysis of the expression consequences from loss of GRAS42 in the gras42-mu1021149 mutant indicated a weak loss of BR signaling in the mutant, consistent with its previously demonstrated role in BR signaling in rice. Loss of BR signaling was also evident by the enhancement of weak BR biosynthetic mutant alleles in double mutants of nana plant1-1 and gras42-mu1021149. The gras42-mu1021149 mutant had little effect on GA-regulated gene expression, suggesting that GRAS42 is not a regulator of core GA signaling genes in maize. Single-cell expression data identified gras42 expressed among cells in the G2/M phase of the cell cycle consistent with its previously demonstrated role in cell cycle gene expression in Arabidopsis (Arabidopsis thaliana). Cis-acting natural variation controlling GRAS42 transcript accumulation was identified by expression genome-wide association study (eGWAS) in maize. Our results demonstrate a conserved role for GRAS42/SCARECROW-LIKE 28 (SCL28)/DLT in BR signaling, clarify the role of this gene in GA signaling, and suggest mechanisms of tillering and leaf angle control by BR.

54 ENVIRONMENTAL SCIENCES↗

Simulation of Radiation-Induced DNA Damage With the Code RITRACKS

INTRODUCTION DNA damage is one of the most physiologically important effects of ionizing radiation. Clustered DNA damage events, like double-strand breaks (DSBs), have the most notable biological consequences. DNA damage types depend on both the track structure of the radiation and the spatial organization of the DNA. High linear energy transfer (LET) charged nuclei, found in galactic cosmic rays (GCR), are known to produce large numbers of complex DNA damage events. The human genome is packaged into chromatin, which can take on locus-dependent and cell type-dependent spatial conformations that correspond to epigenetic states, such as more open, extended structures in transcriptionally active chromatin. These epigenetic differences can affect DNA break patterns in response to ionizing radiation, potentially creating distinct DNA repair and signaling outcomes across the genome in different cells. MATERIAL AND METHODS The code RITRACKS (Relativistic Ion Tracks), which simulates stochastic radiation track structures and radiation chemistry, was used to model damage on isolated and histone-bound DNA by various types of ions and photons. The changes made to the code to perform radiation-induced DNA damage, and simulation results on single nucleosomes are given in our recent paper. In this work, the DNA building capabilities of RITRACKS have been extended to simulate more complex DNA structures build on the coarse-grain simulation framework meso-WLCsim. This code can sample generic chromatin fiber conformation ensembles based on the geometry of nucleosomes and mechanical properties of DNA. Using RITRACKS, we simulated the fragment length distributions (FLD) of irradiated DNA structures built using the chromatin conformations of WLCsim and obtained results representative of those obtained with Radiation-Induced Correlated Cleavage with sequencing (RICC-Seq) experiments [6]. We have also performed Fe ion and photon irradiations of K562, IMR90, BJ and RPE-1 cells at NSRL to experimentally validate results. Sample processing and data analysis are in progress and any available preliminary results will be discussed. DISCUSSION The recent updates in the code RITRACKS allow the calculation of several quantities such as the DNA damage yield and the FLD. This approach can be used to model epigenetic state-specific chromatin structure parameters to leverage the epigenetic state data available for many human cell types to infer relative DNA damage sensitivity among genomic loci.

I Plante↗

Non-photosynthetic lineages sibling to Cyanobacteria associate with eukaryotes in the open ocean

Margulisbacteria are elusive uncultivated bacteria that have illuminated evolutionary transitions in the progenitor of Cyanobacteria, the latter being a critically important phylum that underpins oxygenic photosynthesis. The non-photosynthetic Margulisbacteria were discovered in a sulfidic spring and later in other habitats. Currently, this candidate phylum partitions into the Riflemargulisbacteria, primarily from sediments and groundwater, the Termititenax from insect gut microbiomes, and the Marinamargulisbacteria, from marine samples. We found that Marinamargulisbacteria amplicons were unusually distributed in size-fractionated samples from the sunlit photic and dark twilight zones of the ocean. Further, sequencing of wild marine protists rendered genomic information for distinct marinamargulisbacterial clades co-associated with uncultivated, non-photosynthetic Stramenopila and Opisthokonta protists. Phylogenomic analyses combining these data and available metagenome-assembled genomes (MAGs) and single-amplified genomes (SAGs) from sorted bacteria revealed new Marinamargulisbacteria lineages. The lineages delineate by their environment, forming clades comprising freshwater, marine pelagic, or sediment/hypoxic taxa. In conclusion, the remarkable diversity of Margulisbacteria indicates success in colonizing various habitats, potentially in a conserved strategy involving eukaryotic cells.

59 BASIC BIOLOGICAL SCIENCES↗

Identification of functional non-coding variants associated with orofacial cleft

Oral facial cleft (OFC) comprises cleft lip with or without cleft palate (CL/P) or cleft palate only. Genome wide association studies (GWAS) of isolated OFC have identified common single nucleotide polymorphisms (SNPs) in many genomic loci where the presumed effector gene (for example, IRF6 in the 1q32 locus) is expressed in embryonic oral epithelium. To identify candidates for functional SNPs at eight such loci we conduct a massively parallel reporter assay in a fetal oral epithelial cell line, revealing SNPs with allele-specific effects on enhancer activity. We filter these SNPs against chromatin-mark evidence of enhancers and test a subset in traditional reporter assays, which support the candidacy of SNPs at loci containing FOXE1, IRF6, MAFB, TFAP2A, and TP63. For two SNPs near IRF6 and one near FOXE1, we engineer the genome of induced pluripotent stem cells, differentiate the cells into embryonic oral epithelium, and discover allele-specific effects on the levels of effector gene expression, and, in two cases, the binding affinity of transcription factors FOXE1 or ETS2. Conditional analyses of GWAS data suggest the two functional SNPs near IRF6 account for the majority of risk for CL/P at this locus. This study connects genetic variation associated with OFC to mechanisms of pathogenesis.

Kumari, Priyanka↗

Binding profiles for 961 Drosophila and C. elegans transcription factors reveal tissue-specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here, we present the culmination of the efforts of the modENCODE (model organism Encyclopedia of DNA Elements) and modERN (model organism Encyclopedia of Regulatory Networks) consortia to systematically assay TF binding events in vivo in two major model organisms,Drosophila melanogaster(fly) andCaenorhabditis elegans(worm). These data sets comprise 605 TFs identifying 3.6 M sites in the fly and 356 TFs identifying 0.9 M sites in the worm, and represent the majority of the regulatory space in each genome. We demonstrate that TFs associate with chromatin in clusters termed “metapeaks,” that larger metapeaks have characteristics of high-occupancy target (HOT) regions, and that the importance of consensus sequence motifs bound by TFs depends on metapeak size and complexity. Combining ChIP-seq data with single-cell RNA-seq data in a machine-learning model identifies TFs with a prominent role in promoting target gene expression in specific cell types, even differentiating between parent–daughter cells during embryogenesis. These data are a rich resource for the community that should fuel and guide future investigations into TF function. To facilitate data accessibility and utility, all strains expressing green fluorescent protein (GFP)-tagged TFs are available at the stock centers for each organism. The chromatin immunoprecipitation sequencing data are available through the ENCODE Data Coordinating Center, GEO, and through a direct interface that provides rapid access to processed data sets and summary analyses, as well as widgets to probe the cell-type-specific TF–target relationships.

Biochemistry & Molecular Biology↗

Genome-wide functional screens enable the prediction of high activity CRISPR-Cas9 and -Cas12a guides in Yarrowia lipolytica

Abstract Genome-wide functional genetic screens have been successful in discovering genotype-phenotype relationships and in engineering new phenotypes. While broadly applied in mammalian cell lines and in E. coli , use in non-conventional microorganisms has been limited, in part, due to the inability to accurately design high activity CRISPR guides in such species. Here, we develop an experimental-computational approach to sgRNA design that is specific to an organism of choice, in this case the oleaginous yeast Yarrowia lipolytica . A negative selection screen in the absence of non-homologous end-joining, the dominant DNA repair mechanism, was used to generate single guide RNA (sgRNA) activity profiles for both SpCas9 and LbCas12a. This genome-wide data served as input to a deep learning algorithm, DeepGuide, that is able to accurately predict guide activity. DeepGuide uses unsupervised learning to obtain a compressed representation of the genome, followed by supervised learning to map sgRNA sequence, genomic context, and epigenetic features with guide activity. Experimental validation, both genome-wide and with a subset of selected genes, confirms DeepGuide’s ability to accurately predict high activity sgRNAs. DeepGuide provides an organism specific predictor of CRISPR guide activity that with retraining could be applied to other fungal species, prokaryotes, and other non-conventional organisms.

59 BASIC BIOLOGICAL SCIENCES↗