Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DNA sequencing”

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 577 records · Page 32

CRISPR/Cas9-Mediated Multi-Allelic Gene Targeting in Sugarcane Confers Herbicide Tolerance

Sugarcane is the source of 80% of the sugar and 26% of the bioethanol produced globally. However, its complex, highly polyploid genome (2 n = 100 – 120) impedes crop improvement. Here, we report efficient and reproducible gene targeting (GT) in sugarcane, enabling precise co-editing of multiple alleles via template-mediated and homology-directed repair (HDR) of DNA double strand breaks induced by the programmable nuclease CRISPR/Cas9. The evaluation of 146 independently transformed plants from five independent experiments revealed a targeted nucleotide replacement that resulted in both targeted amino acid substitutions W574L and S653I in the acetolactate synthase (ALS) in 11 lines in addition to single, targeted amino acid substitutions W574L or S653I in 25 or 18 lines, respectively. Co-editing of up to three ALS copies/alleles that confer herbicide tolerance was confirmed by Sanger sequencing of cloned long polymerase chain reaction (PCR) amplicons. This work will enable crop improvement by conversion of inferior alleles to superior alleles through targeted nucleotide substitutions.

(Saccharum hybrid complex)↗

ADEPT: a domain independent sequence alignment strategy for gpu architectures

Bioinformatic workflows frequently make use of automated genome assembly and protein clustering tools. At the core of most of these tools, a significant portion of execution time is spent in determining optimal local alignment between two sequences. This task is performed with the Smith-Waterman algorithm, which is a dynamic programming based method. With the advent of modern sequencing technologies and increasing size of both genome and protein databases, a need for faster Smith-Waterman implementations has emerged. Multiple SIMD strategies for the Smith-Waterman algorithm are available for CPUs. However, with the move of HPC facilities towards accelerator based architectures, a need for an efficient GPU accelerated strategy has emerged. Existing GPU based strategies have either been optimized for a specific type of characters (Nucleotides or Amino Acids) or for only a handful of application use-cases. In this paper, we present ADEPT, a new sequence alignment strategy for GPU architectures that is domain independent, supporting alignment of sequences from both genomes and proteins. Our proposed strategy uses GPU specific optimizations that do not rely on the nature of sequence. We demonstrate the feasibility of this strategy by implementing the Smith-Waterman algorithm and comparing it to similar CPU strategies as well as the fastest known GPU methods for each domain. ADEPT’s driver enables it to scale across multiple GPUs and allows easy integration into software pipelines which utilize large scale computational systems. We have shown that the ADEPT based Smith-Waterman algorithm demonstrates a peak performance of 360 GCUPS and 497 GCUPs for protein based and DNA based datasets respectively on a single GPU node (8 GPUs) of the Cori Supercomputer. Overall ADEPT shows 10x faster performance in a node-to-node comparison against a corresponding SIMD CPU implementation. ADEPT demonstrates a performance that is either comparable or better than existing GPU strategies. We demonstrated the efficacy of ADEPT in supporting existing bionformatics software pipelines by integrating ADEPT in MetaHipMer a high-performance denovo metagenome assembler and PASTIS a high-performance protein similarity graph construction pipeline. Our results show 10% and 30% boost of performance in MetaHipMer and PASTIS respectively.

59 BASIC BIOLOGICAL SCIENCES↗

RCSB protein data Bank: Next‐generation advanced search for exploration of experimental structures and computed structure models

Abstract The Protein Data Bank (PDB), established in 1971, is the primary global, open‐access archive for experimentally determined 3D macromolecular structures (proteins, RNA, DNA). The research‐focused RCSB.org web‐portal provides access to these data alongside more than one million machine‐learning‐predicted structure models, greatly expanding the available structural landscape. Rapid growth of both experimental and computational structures has increased the need for powerful yet accessible search tools that serve a broad and diverse scientific community. Herein, we describe a redesigned RCSB Protein Data Bank RCSB.org Advanced Search capability that supports intuitive discovery of 3D structures through a unified interface. This interface integrates annotation‐, sequence‐, and 3D structure‐based searches, embeds an interactive 3D viewer, and incorporates curated biological knowledge, such as catalytic site definitions from Mechanism and Catalytic Site Atlas and ligand‐guided structural motifs, for constructing geometry‐driven queries. A new Chemical Search tool allows definition of chemical queries via an integrated drawing tool or standard identifiers, seamlessly combining them with annotation filters. By allowing query definition directly within spatial and chemical contexts, these search interfaces reduce the need for detailed knowledge of residue numbering, chain identifiers, or external cheminformatics software. This capability enables efficient exploration of structures, chemical diversity, and structure–function relationships across all life domains. The redesigned interfaces can be accessed directly at rcsb.org/search/advanced for Advanced Search and rcsb.org/search/chemical for Chemical Search.

Rose, Yana [Research Collaboratory for Structural ↗

Worldwide Population Genomics Reveal Long-Term Stability of the Mitochondrial Genome Architecture in a Keystone Marine Plant

Mitochondrial genomes (mitogenomes) of flowering plants are composed of multiple chromosomes. Recombination within and between the mitochondrial chromosomes may generate diverse DNA molecules termed isoforms. The isoform copy number and composition can be dynamic within and among individual plants due to uneven replication and homologous recombination. Nonetheless, despite their functional importance, the level of mitogenome conservation within species remains understudied. Whether the ontogenetic variation translates to evolution of mitogenome composition over generations is currently unknown. Here we show that the mitogenome composition of the seagrass Zostera marina is conserved among worldwide populations that diverged ca. 350,000 years ago. Using long-read sequencing, we characterized the Z. marina mitochondrial genome and inferred the repertoire of recombination-induced configurations. To characterize the mitochondrial genome architecture worldwide and study its evolution, we examined the mitogenome in Z. marina meristematic region sampled in 16 populations from the Pacific and Atlantic oceans. Our results reveal a striking similarity in the isoform relative copy number, indicating a high conservation of the mitogenome composition among distantly related populations and within the plant germline, despite a notable variability during individual ontogenesis. Our study supplies a link between observations of dynamic mitogenomes at the level of plant individuals and long-term mitochondrial evolution.

(sub)stoichiometric shifting↗

Reekeekee- and roodoodooviruses, two different Microviridae clades constituted by the smallest DNA phages

Small circular single-stranded DNA viruses of the Microviridae family are both prevalent and diverse in all ecosystems. They usually harbor a genome between 4.3 and 6.3 kb, with a microvirus recently isolated from a marine Alphaproteobacteria being the smallest known genome of a DNA phage (4.248 kb). A subfamily, Amoyvirinae, has been proposed to classify this virus and other related small Alphaproteobacteria-infecting phages. Here, we report the discovery, in meta-omics data sets from various aquatic ecosystems, of sixteen complete microvirus genomes significantly smaller (2.991–3.692 kb) than known ones. Phylogenetic analysis reveals that these sixteen genomes represent two related, yet distinct and diverse, novel groups of microviruses—amoyviruses being their closest known relatives. We propose that these small microviruses are members of two tentatively named subfamilies Reekeekeevirinae and Roodoodoovirinae. As known microvirus genomes encode many overlapping and overprinted genes that are not identified by gene prediction software, we developed a new methodology to identify all genes based on protein conservation, amino acid composition, and selection pressure estimations. Surprisingly, only four to five genes could be identified per genome, with the number of overprinted genes lower than that in phiX174. These small genomes thus tend to have both a lower number of genes and a shorter length for each gene, leaving no place for variable gene regions that could harbor overprinted genes. Even more surprisingly, these two Microviridae groups had specific and different gene content, and major differences in their conserved protein sequences, highlighting that these two related groups of small genome microviruses use very different strategies to fulfill their lifecycle with such a small number of genes. The discovery of these genomes and the detailed prediction and annotation of their genome content expand our understanding of ssDNA phages in nature and are further evidence that these viruses have explored a wide range of possibilities during their long evolution.

59 BASIC BIOLOGICAL SCIENCES↗

Sequencing and Analysis of the Sex Determination Region of Populus trichocarpa

The ages and sizes of a sex-determination region (SDR) are difficult to determine in non-model species. Due to the lack of recombination and enrichment of repetitive elements in SDRs, the quality of assembly with short sequencing reads is universally low. Unique features present in the SDRs help provide clues about how SDRs are established and how they evolve in the absence of recombination. Several Populus species have been reported with a male heterogametic configuration of sex (XX/XY system) mapped on chromosome 19, but the exact location of the SDR has been inconsistent among species, and thus far, none of these SDRs has been fully assembled in a genomic context. Here we identify the Y-SDR from a Y-linked contig directly from a long-read PacBio assembly of a Populus trichocarpa male individual. We also identified homologous gene sequences in the SDR of P. trichocarpa and the SDR of the W chromosome in Salix purpurea. We show that inverted repeats (IRs) found in the Y-SDR and the W-SDR are lineage-specific. We hypothesize that, although the two IRs are derived from the same orthologous gene within each species, they likely have independent evolutionary histories. Furthermore, the truncated inverted repeats in P. trichocarpa may code for small RNAs that target the homologous gene for RNA-directed DNA methylation. These findings support the hypothesis that diverse sex-determining systems may be achieved through similar evolutionary pathways, thereby providing a possible mechanism to explain the lability of sex-determination systems in plants in general.

59 BASIC BIOLOGICAL SCIENCES↗

Integration of Complete Plasmids Containing Bont Genes into Chromosomes of Clostridium parabotulinum , Clostridium sporogenes , and Clostridium argentinense

At least 40 toxin subtypes of botulinum neurotoxins (BoNTs), a heterogenous group of bacterial proteins, are produced by seven different clostridial species. A key factor that drives the diversity of neurotoxigenic clostridia is the association of bont gene clusters with various genomic locations including plasmids, phages and the chromosome. Analysis of Clostridium sporogenes BoNT/B1 strain CDC 1632, C. argentinense BoNT/G strain CDC 2741, and Clostridium parabotulinum BoNT/B1 strain DFPST0006 genomes revealed bont gene clusters within plasmid-like sequences within the chromosome or nested in large contigs, with no evidence of extrachromosomal elements. A nucleotide sequence (255,474 bp) identified in CDC 1632 shared 99.5% identity (88% coverage) with bont/B1-containing plasmid pNPD7 of C. sporogenes CDC 67071; CDC 2741 contig AYSO01000020 (1.1 MB) contained a ~140 kb region which shared 99.99% identity (100% coverage) with plasmid pRSJ17_1 of C. argentinense BoNT/G strain 89G; and DFPST0006 contig JACBDK0100002 (573 kb) contained a region that shared 100% identity (99%) coverage with the bont/B1-containing plasmid pCLD of C. parabotulinum Okra. This is the first report of full-length plasmid DNA-carrying complete neurotoxin gene clusters integrated in three distinct neurotoxigenic species: C. parabotulinum, C. sporogenes and C. argentinense.

59 BASIC BIOLOGICAL SCIENCES↗

Nickase fidelity drives EvolvR-mediated diversification in mammalian cells

Abstract In vivo genetic diversifiers have previously enabled efficient searches of genetic variant fitness landscapes for continuous directed evolution. However, existing genomic diversification modalities for mammalian genomic loci exclusively rely on deaminases to generate transition mutations within target loci, forfeiting access to most missense mutations. Here, we engineer CRISPR-guided error-prone DNA polymerases (EvolvR) to diversify all four nucleotides within genomic loci in mammalian cells. We demonstrate that EvolvR generates both transition and transversion mutations throughout a mutation window of at least 40 bp and implement EvolvR to evolve previously unreported drug-resistantMAP2K1variants via substitutions not achievable with deaminases. Moreover, we discover that the nickase’s mismatch tolerance limits EvolvR’s mutation window and substitution biases in a gRNA-specific fashion. To compensate for gRNA-to-gRNA variability in mutagenesis, we maximize the number of gRNA target sequences by incorporating a PAM-flexible nickase into EvolvR. Finally, we find a strong correlation between predicted free energy changes underlying R-loop formation and EvolvR’s performance using a given gRNA. The EvolvR system diversifies all four nucleotides to enable the evolution of mammalian cells, while nuclease and gRNA-specific properties underlying nickase fidelity can be engineered to further enhance EvolvR’s mutation rates.

Science & Technology - Other Topics↗

Pseudodesulfovibrio mercurii sp. nov., a mercury-methylating bacterium isolated from sediment

The sulfate-reducing, mercury-methylating strain ND132 T was isolated from the brackish anaerobic bottom sediments of Chesapeake Bay, USA. Capable of high levels of mercury (Hg) methylation, ND132 T has been widely used as a model strain to study the process and to determine the genetic basis of Hg methylation. Originally called Desulfovibrio desulfuricans ND132 T on the basis of an early partial 16S rRNA sequence, the strain has never been formally described. Phylogenetic and physiological traits place this strain within the genus Pseudodesulfovibrio, in the recently reclassified phylum Desulfobacterota (formerly Deltaproteobacteria ). ND132 T is most closely related to Pseudodesulfovibrio hydrargyri BerOc1 T and Pseudodesulfovibrio indicus J2T. Analysis of average nucleotide identity (ANI) of whole-genome sequences showed roughly 88 % ANI between P. hydrargyri BerOc1T and ND132T, and 84 % similarity between ND132T and P. indicus J2 T . These cut-off scores <95 %, along with a multi-gene phylogenetic analysis of members of the family Desulfovibrionacea, and differences in physiology indicate that all three strains represent separate species. The Gram-stain-negative cells are vibrio-shaped, motile and not sporulated. ND132 T is a salt-tolerant mesophile with optimal growth in the laboratory at 32 °C, 2 % salinity, and pH 7.8. The DNA G+C content of the genomic DNA is 65.2 %. It is an incomplete oxidizer of short chain fatty acids, using lactate, pyruvate and fumarate with sulfate or sulfite as the terminal electron acceptors. ND132 T can respire fumarate using pyruvate as an electron donor. The major fatty acids are iso-C 15 :0 , anteiso-C 15 : 0 , iso-C 17 : 0 , iso-C 17 : 1 ω9c and anteiso-C 17 : 0 . We propose the classification of strain ND132 T (DSM 110689, ATCC TSD-224) as the type strain Pseudodesulfovibrio mercurii sp. nov.

59 BASIC BIOLOGICAL SCIENCES↗

Simple biochemical features underlie transcriptional activation domain diversity and dynamic, fuzzy binding to Mediator

Gene activator proteins comprise distinct DNA-binding and transcriptional activation domains (ADs). Because few ADs have been described, we tested domains tiling all yeast transcription factors for activation in vivo and identified 150 ADs. By mRNA display, we showed that 73% of ADs bound the Med15 subunit of Mediator, and that binding strength was correlated with activation. AD-Mediator interaction in vitro was unaffected by a large excess of free activator protein, pointing to a dynamic mechanism of interaction. Structural modeling showed that ADs interact with Med15 without shape complementarity (‘fuzzy’ binding). ADs shared no sequence motifs, but mutagenesis revealed biochemical and structural constraints. Finally, a neural network trained on AD sequences accurately predicted ADs in human proteins and in other yeast proteins, including chromosomal proteins and chromatin remodeling complexes. These findings solve the longstanding enigma of AD structure and function and provide a rationale for their role in biology.

60 APPLIED LIFE SCIENCES↗

Viruses interact with hosts that span distantly related microbial domains in dense hydrothermal mats

Abstract Many microbes in nature reside in dense, metabolically interdependent communities. We investigated the nature and extent of microbe-virus interactions in relation to microbial density and syntrophy by examining microbe-virus interactions in a biomass dense, deep-sea hydrothermal mat. Using metagenomic sequencing, we find numerous instances where phylogenetically distant (up to domain level) microbes encode CRISPR-based immunity against the same viruses in the mat. Evidence of viral interactions with hosts cross-cutting microbial domains is particularly striking between known syntrophic partners, for example those engaged in anaerobic methanotrophy. These patterns are corroborated by proximity-ligation-based (Hi-C) inference. Surveys of public datasets reveal additional viruses interacting with hosts across domains in diverse ecosystems known to harbour syntrophic biofilms. We propose that the entry of viral particles and/or DNA to non-primary host cells may be a common phenomenon in densely populated ecosystems, with eco-evolutionary implications for syntrophic microbes and CRISPR-mediated inter-population augmentation of resilience against viruses.

59 BASIC BIOLOGICAL SCIENCES↗

Ancient co-option of LTR retrotransposons as yeast centromeres

Centromeres ensure accurate chromosome segregation, yet their DNA evolves rapidly across eukaryotes leaving the origins of new centromere architectures unclear. The brewer’s yeast Saccharomyces cerevisiae exemplifies this long-standing puzzle. Its centromeres shifted ancestrally from large, repeat-rich, epigenetically specified forms to the compact, genetically defined ‘point’ centromeres. How this transition occurred has remained unresolved6. Here we identify evolutionarily related ‘proto-point’ centromeres that provide a resolution to the evolutionary origins of point centromeres. Proto-point centromeres contain a single centromeric nucleosome positioned over an AT-rich core, accompanied by relaxed organization and sequence variability of flanking cis-elements. In two species, these proto-point centromeres lie within retrotransposon-derived repeat clusters, linking ancestral repeat-rich centromeres to genetically encoded ones. Comparative and phylogenetic analyses indicate that proto-point and point centromeres evolved in an ancestor with retrotransposon-rich centromeres. These results identify long-terminal-repeat retrotransposons, specifically Ty5 sequences, as the genetic substrate for point-centromere evolution and provide a mechanistic route by which an epigenetic centromere can become genetically specified. More broadly, they show how selfish elements can be co-opted to perform essential chromosomal functions.

Haase, Max A. B. [Max Planck Institute of Molecula↗

The fate of 35S rRNA genes in the allotetraploid grass Brachypodium hybridum

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

35S rDNA evolution↗

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↗

KBase Narrative - Atacama Trench & Sn Clemente 2018 cruise Metagenomes -Sequence quality filtration & Assembly-

Deep-sea microbiomes sampled with traditional or pressure-retaining instrumentation A fundamental property of life in the sea is that with increasing depth must come increased adaptation to elevated hydrostatic pressures. This project aims to assess technology for the accurate characterization of microbial biodiversity and biochemistry at greater depths, and validate its use in some of the deepest and least-explored waters in the Ocean: Hadal Trench. Compared to the 0.1 megapascal (MPa) of surface waters, organisms inhabiting at the average depth of the Ocean (3800 m) have to deal with pressures of 38 MPa (5,500 psi). The effects of these high pressure have been most thoroughly investigated in microorganisms, where it has been found to affect membrane integrity and transport, cytoskeletal assembly and both DNA and protein synthesis. We will explore the diversity and activity of microorganisms from 2000 m to >8,000 m depths within the highly productive water column of hadal trenches, comparing seawater samples obtained with and without substantial pressure changes during recovery. A pressure retaining device was used to sample the microbial communities of the watercolumn with reduced decompression during retreival (Peoples et. al., 2019). References: Peoples LM, Norenberg M, Price D, McGoldrick M, Novotny M, Bochdansky A, Bartlett DH. A full-ocean-depth rated modular lander and pressure-retaining sampler capable of collecting hadal-endemic microbes under in situ conditions. Deep Sea Research Part I: Oceanographic Research Papers. 2019 Jan 1;143:50-7.

Plominsky, Alvaro M.↗

A new genus and species of xanthoid crab (Decapoda: Brachyura) from offshore hard bank habitats in the Gulf of Mexico

Speciose populations of small xanthoid crabs on offshore banks and reefs of the northern Gulf of Mexico include a new species that is not assignable to presently named genera. Morphological diagnoses of the new genus and species are underpinned by previously published gene sequence analyses, originally misattributed to another species but now known to apply to this taxon. Herein named Guinope tiara n. gen., n. sp., the species shows molecular phylogenetic affinities with the family Linnaeoxanthidae Stevcic, 2005, an ally of panopeid and pseudorhombilid crabs. Specimens from Occulina banks off the Florida Atlantic coast, previously regarded to represent Garthiope barbadensis (Rathbun, 1921), are not that species but instead morphologically assignable to Guinope n. gen. Whether they represent variants of Guinope tiara n. gen, n. sp. or a second species of the genus awaits the collection of fresh materials for DNA analyses.

Zoology↗

Bacterial genome editing by coupling Cre-lox and CRISPR-Cas9 systems

The past decade has been a golden age for microbiology, marked by the discovery of an unprecedented increase in the number of novel bacterial species. Yet gaining biological knowledge of those organisms has not kept pace with sequencing efforts. To unlock this genetic potential there is an urgent need for generic (i.e. non-species specific) genetic toolboxes. Recently, we developed a method, termed chassis-independent recombinase-assisted genome engineering (CRAGE), enabling the integration and expression of large complex gene clusters directly into the chromosomes of diverse bacteria. Here we expand upon this technology by incorporating CRISPR-Cas9 allowing precise genome editing across multiple bacterial species. To do that we have developed a landing pad that carries one wild-type and two mutant lox sites to allow integration of foreign DNA at two locations through Cre-lox recombinase-mediated cassette exchange (RMCE). The first RMCE event is to integrate the Cas9 and the DNA repair protein genes RecET, and the second RMCE event enables the integration of customized sgRNA and a repair template. Following this workflow, we achieved precise genome editing in four different gammaproteobacterial species. We also show that the inserted landing pad and the entire editing machinery can be removed scarlessly after editing. We report here the construction of a single landing pad transposon and demonstrate its functionality across multiple species. The modular design of the landing pad and accessory vectors allows design and assembly of genome editing platforms for other organisms in a similar way. We believe this approach will greatly expand the list of bacteria amenable to genetic manipulation and provides the means to advance our understanding of the microbial world.

59 BASIC BIOLOGICAL SCIENCES↗

The protist Aurantiochytrium has universal subtelomeric rDNAs and is a host for mirusviruses

Viruses are the most abundant biological entities in the world's oceans, where they play important ecological and biogeochemical roles. Metagenomics is revealing new groups of eukaryotic viruses, although disconnected from known hosts. Among these are the recently described mirusviruses, which share some similarities with herpesviruses.50 years ago, "herpes-type" viral particles were found in a thraustochytrid member of the labyrinthulomycetes, a diverse group of abundant and ecologically important marine eukaryotes, but could not be further characterized by methods then available. Long-read sequencing has allowed us to connect the biology of mirusviruses and thraustochytrids. We sequenced the genome of the genetically tractable model thraustochytrid Aurantiochytrium limacinum ATCC MYA-1381 and found that its 26 linear chromosomes have an extraordinary configuration. Subtelomeric ribosomal DNAs (rDNAs) found at all chromosome ends are interspersed with long repeated sequence elements denoted as long repeated-telomere and rDNA spacers (LORE-TEARS). We identified two genomic elements that are related to mirusvirus genomes. The first is a ~300-kbp episome (circular element 1 [CE1]) present at a high copy number. Strikingly, the second, distinct, mirusvirus-like element is integrated between two sets of rDNAs and LORE-TEARS at the left end of chromosome 15 (LE-Chr15). Similar to metagenomically derived mirusviruses, these putative A. limacinum mirusviruses have a virion module related to that of herpesviruses along with an informational module related to nucleocytoplasmic large DNA viruses (NCLDVs). CE1 and LE-Chr15 bear striking similarities to episomal and endogenous latent forms of herpesviruses, respectively, and open new avenues of research into marine virus-host interactions.

59 BASIC BIOLOGICAL SCIENCES↗