The role of the Cer1 transposon in horizontal transfer of transgenerational memory.

Animals face both external and internal dangers: pathogens threaten from the environment, and unstable genomic elements threaten from within. C. elegans protects itself from pathogens by "reading" bacterial small RNAs, using this information to both induce avoidance and transmit memories for four generations. Here, we found that memories can be transferred from either lysed animals or from conditioned media to naive animals via Cer1 retrotransposon-encoded virus-like particles. Moreover, Cer1 functions internally at the step of transmission of information from the germline to neurons and is required for learned avoidance. The presence of the Cer1 retrotransposon in wild C. elegans strains correlates with the ability to learn and inherit small-RNA-induced pathogen avoidance. Together, these results suggest that C. elegans has co-opted a potentially dangerous retrotransposon to instead protect itself and its progeny from a common pathogen through its inter-tissue signaling ability, hijacking this genomic element for its own adaptive immunity benefit.

Improving prime editing with an endogenous small RNA-binding protein.

Prime editing enables the precise modification of genomes through reverse transcription of template sequences appended to the 3' ends of CRISPR-Cas guide RNAs1. To identify cellular determinants of prime editing, we developed scalable prime editing reporters and performed genome-scale CRISPR-interference screens. From these screens, a single factor emerged as the strongest mediator of prime editing: the small RNA-binding exonuclease protection factor La. Further investigation revealed that La promotes prime editing across approaches (PE2, PE3, PE4 and PE5), edit types (substitutions, insertions and deletions), endogenous loci and cell types but has no consistent effect on genome-editing approaches that rely on standard, unextended guide RNAs. Previous work has shown that La binds polyuridine tracts at the 3' ends of RNA polymerase III transcripts2. We found that La functionally interacts with the 3' ends of polyuridylated prime editing guide RNAs (pegRNAs). Guided by these results, we developed a prime editor protein (PE7) fused to the RNA-binding, N-terminal domain of La. This editor improved prime editing with expressed pegRNAs and engineered pegRNAs (epegRNAs), as well as with synthetic pegRNAs optimized for La binding. Together, our results provide key insights into how prime editing components interact with the cellular environment and suggest general strategies for stabilizing exogenous small RNAs therein.

C. elegans interprets bacterial non-coding RNAs to learn pathogenic avoidance.

Caenorhabditis elegans must distinguish pathogens from nutritious food sources among the many bacteria to which it is exposed in its environment1. Here we show that a single exposure to purified small RNAs isolated from pathogenic Pseudomonas aeruginosa (PA14) is sufficient to induce pathogen avoidance in the treated worms and in four subsequent generations of progeny. The RNA interference (RNAi) and PIWI-interacting RNA (piRNA) pathways, the germline and the ASI neuron are all required for avoidance behaviour induced by bacterial small RNAs, and for the transgenerational inheritance of this behaviour. A single P. aeruginosa non-coding RNA, P11, is both necessary and sufficient to convey learned avoidance of PA14, and its C. elegans target, maco-1, is required for avoidance. Our results suggest that this non-coding-RNA-dependent mechanism evolved to survey the microbial environment of the worm, use this information to make appropriate behavioural decisions and pass this information on to its progeny.

Comprehensive single-cell transcriptome lineages of a proto-vertebrate.

Ascidian embryos highlight the importance of cell lineages in animal development. As simple proto-vertebrates, they also provide insights into the evolutionary origins of cell types such as cranial placodes and neural crest cells. Here we have determined single-cell transcriptomes for more than 90,000 cells that span the entirety of development-from the onset of gastrulation to swimming tadpoles-in Ciona intestinalis. Owing to the small numbers of cells in ascidian embryos, this represents an average of over 12-fold coverage for every cell at every stage of development. We used single-cell transcriptome trajectories to construct virtual cell-lineage maps and provisional gene networks for 41 neural subtypes that comprise the larval nervous system. We summarize several applications of these datasets, including annotating the synaptome of swimming tadpoles and tracing the evolutionary origin of cell types such as the vertebrate telencephalon.

AccuCor2: isotope natural abundance correction for dual-isotope tracer experiments.

Stable isotope labeling techniques have been widely applied in the field of metabolomics and proteomics. Before the measured mass spectral data can be used for quantitative analysis, it must be accurately corrected for isotope natural abundance and tracer isotopic impurity. Despite the increasing popularity of dual-isotope tracing strategy such as 13C-15N or 13C-2H, there are no accurate tools for correcting isotope natural abundance for such experiments in a resolution-dependent manner. Here, we present AccuCor2 as an R-based tool to perform the correction for 13C-15N or 13C-2H labeling experiments. Our method uses a newly designed algorithm to construct the correction matrices that link labeling pattern and measured mass fractions, then use non-negative least-squares to solve the labeling patterns. Our results show that the dual-isotope experiments often require a mass resolution that is high enough to resolve 13C and 15N or 13C and 2H. Otherwise, the labeling pattern is not solvable. However, this mass resolution may not be sufficiently high to resolve other non-tracer elements such as oxygen or sulfur from the tracer elements. Therefore, we design AccuCor2 to perform the correction based on the actual mass resolution of the measurements. Using both simulated and experimental data, we show that AccuCor2 performs accurate and resolution-dependent correction for dual-isotope tracer data.

Analysis of Host Responses to Hepatitis B and Delta Viral Infections in a Micro-scalable Hepatic Co-culture System.

Hepatitis B virus (HBV) remains a major global health problem with 257 million chronically infected individuals worldwide, of whom approximately 20 million are co-infected with hepatitis delta virus (HDV). Progress toward a better understanding of the complex interplay between these two viruses and the development of novel therapies have been hampered by the scarcity of suitable cell culture models that mimic the natural environment of the liver. Here, we established HBV and HBV/HDV co-infections and super-infections in self-assembling co-cultured primary human hepatocytes (SACC-PHHs) for up to 28 days in a 384-well format and highlight the suitability of this platform for high-throughput drug testing. We performed RNA sequencing at days 8 and 28 on SACC-PHHs, either HBV mono-infected or HBV/HDV co-infected. Our transcriptomic analysis demonstrates that hepatocytes in SACC-PHHs maintain a mature hepatic phenotype over time, regardless of infection condition. We confirm that HBV is a stealth virus, as it does not induce a strong innate immune response; rather, oxidative phosphorylation and extracellular matrix-receptor interactions are dysregulated to create an environment that promotes persistence. Notably, HDV co-infection also did not lead to statistically significant transcriptional changes across multiple donors and replicates. The lack of innate immune activation is not due to SACC-PHHs being impaired in their ability to induce interferon stimulated genes (ISGs). Rather, polyinosinic:polycytidylic acid exposure activates ISGs, and this stimulation significantly inhibits HBV infection, yet only minimally affects the ability of HDV to infect and persist. Conclusion: These data demonstrate that the SACC-PHH system is a versatile platform for studying HBV/HDV co-infections and holds promise for performing chemical library screens and improving our understanding of the host response to such infections.

Extrachromosomal circular DNA is common in yeast.

Examples of extrachromosomal circular DNAs (eccDNAs) are found in many organisms, but their impact on genetic variation at the genome scale has not been investigated. We mapped 1,756 eccDNAs in the Saccharomyces cerevisiae genome using Circle-Seq, a highly sensitive eccDNA purification method. Yeast eccDNAs ranged from an arbitrary lower limit of 1 kb up to 38 kb and covered 23% of the genome, representing thousands of genes. EccDNA arose both from genomic regions with repetitive sequences ≥ 15 bases long and from regions with short or no repetitive sequences. Some eccDNAs were identified in several yeast populations. These eccDNAs contained ribosomal genes, transposon remnants, and tandemly repeated genes (HXT6/7, ENA1/2/5, and CUP1-1/-2) that were generally enriched on eccDNAs. EccDNAs seemed to be replicated and 80% contained consensus sequences for autonomous replication origins that could explain their maintenance. Our data suggest that eccDNAs are common in S. cerevisiae, where they might contribute substantially to genetic variation and evolution.

The Oxytricha trifallax macronuclear genome: a complex eukaryotic genome with 16,000 tiny chromosomes.

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor "silent" germline micronuclear genome by a process of "unscrambling" and fragmentation. The tiny macronuclear "nanochromosomes" typically encode single, protein-coding genes (a small portion, 10%, encode 2-8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.

Mapping the pericentric heterochromatin by comparative genomic hybridization analysis and chromosome deletions in Drosophila melanogaster.

Heterochromatin represents a significant portion of eukaryotic genomes and has essential structural and regulatory functions. Its molecular organization is largely unknown due to difficulties in sequencing through and assembling repetitive sequences enriched in the heterochromatin. Here we developed a novel strategy using chromosomal rearrangements and embryonic phenotypes to position unmapped Drosophila melanogaster heterochromatic sequence to specific chromosomal regions. By excluding sequences that can be mapped to the assembled euchromatic arms, we identified sequences that are specific to heterochromatin and used them to design heterochromatin specific probes ("H-probes") for microarray. By comparative genomic hybridization (CGH) analyses of embryos deficient for each chromosome or chromosome arm, we were able to map most of our H-probes to specific chromosome arms. We also positioned sequences mapped to the second and X chromosomes to finer intervals by analyzing smaller deletions with breakpoints in heterochromatin. Using this approach, we were able to map >40% (13.9 Mb) of the previously unmapped heterochromatin sequences assembled by the whole-genome sequencing effort on arm U and arm Uextra to specific locations. We also identified and mapped 110 kb of novel heterochromatic sequences. Subsequent analyses revealed that sequences located within different heterochromatic regions have distinct properties, such as sequence composition, degree of repetitiveness, and level of underreplication in polytenized tissues. Surprisingly, although heterochromatin is generally considered to be transcriptionally silent, we detected region-specific temporal patterns of transcription in heterochromatin during oogenesis and early embryonic development. Our study provides a useful approach to elucidate the molecular organization and function of heterochromatin and reveals region-specific variation of heterochromatin.

Pseudomonas fluorescens 15 small RNA Pfs1 mediates transgenerational epigenetic inheritance of pathogen avoidance in C. elegans through the Ephrin receptor VAB-1.

C. elegans are exposed to a variety of pathogenic and non-pathogenic bacteria species in their natural environment. Correspondingly, C. elegans has evolved an ability to discern between nutritive and infectious bacterial food sources. Here we show that C. elegans can learn to avoid the pathogenic bacteria Pseudomonas fluorescens 15 (PF15), and that this learned avoidance behavior is passed on to progeny for four generations, as we previously demonstrated for Pseudomonas aeruginosa (PA14) and Pseudomonas vranovensis, using similar mechanisms, including the involvement of both the TGF-ß ligand DAF-7 and Cer1 retrotransposon-encoded virus-like particles. PF15 small RNAs are both necessary and sufficient to induce this transgenerational avoidance behavior. Unlike PA14 or P. vranovensis, PF15 does not use P11, Pv1, or a small RNA with maco-1 homology for this avoidance; instead, an unrelated PF15 small RNA, Pfs1, that targets the C. elegans vab-1 Ephrin receptor gene is necessary and sufficient for learned avoidance, suggesting the evolution of yet another bacterial sRNA/C. elegans gene target pair involved in transgenerational inheritance of pathogen avoidance. As VAB-2 Ephrin receptor ligand and MACO-1 knockdown also induce PF15 avoidance, we have begun to understand the genetic pathway involved in small RNA targeted pathogenic avoidance. Moreover, these data show that axon guidance pathway genes (VAB-1 and VAB-2) have previously unknown adult roles in regulating neuronal function. C. elegans may have evolved multiple bacterial specificity-encoded small RNA-dependent mechanisms to avoid different pathogenic bacteria species, thereby providing progeny with a survival advantage in a dynamic environment.

A wide extent of inter-strain diversity in virulent and vaccine strains of alphaherpesviruses.

Alphaherpesviruses are widespread in the human population, and include herpes simplex virus 1 (HSV-1) and 2, and varicella zoster virus (VZV). These viral pathogens cause epithelial lesions, and then infect the nervous system to cause lifelong latency, reactivation, and spread. A related veterinary herpesvirus, pseudorabies (PRV), causes similar disease in livestock that result in significant economic losses. Vaccines developed for VZV and PRV serve as useful models for the development of an HSV-1 vaccine. We present full genome sequence comparisons of the PRV vaccine strain Bartha, and two virulent PRV isolates, Kaplan and Becker. These genome sequences were determined by high-throughput sequencing and assembly, and present new insights into the attenuation of a mammalian alphaherpesvirus vaccine strain. We find many previously unknown coding differences between PRV Bartha and the virulent strains, including changes to the fusion proteins gH and gB, and over forty other viral proteins. Inter-strain variation in PRV protein sequences is much closer to levels previously observed for HSV-1 than for the highly stable VZV proteome. Almost 20% of the PRV genome contains tandem short sequence repeats (SSRs), a class of nucleic acids motifs whose length-variation has been associated with changes in DNA binding site efficiency, transcriptional regulation, and protein interactions. We find SSRs throughout the herpesvirus family, and provide the first global characterization of SSRs in viruses, both within and between strains. We find SSR length variation between different isolates of PRV and HSV-1, which may provide a new mechanism for phenotypic variation between strains. Finally, we detected a small number of polymorphic bases within each plaque-purified PRV strain, and we characterize the effect of passage and plaque-purification on these polymorphisms. These data add to growing evidence that even plaque-purified stocks of stable DNA viruses exhibit limited sequence heterogeneity, which likely seeds future strain evolution.

Metabolic cycling in single yeast cells from unsynchronized steady-state populations limited on glucose or phosphate.

Oscillations in patterns of expression of a large fraction of yeast genes are associated with the "metabolic cycle," usually seen only in prestarved, continuous cultures of yeast. We used FISH of mRNA in individual cells to test the hypothesis that these oscillations happen in single cells drawn from unsynchronized cultures growing exponentially in chemostats. Gene-expression data from synchronized cultures were used to predict coincident appearance of mRNAs from pairs of genes in the unsynchronized cells. Quantitative analysis of the FISH results shows that individual unsynchronized cells growing slowly because of glucose limitation or phosphate limitation show the predicted oscillations. We conclude that the yeast metabolic cycle is an intrinsic property of yeast metabolism and does not depend on either synchronization or external limitation of growth by the carbon source.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees.

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

Herpes simplex virus 1 pUL34 plays a critical role in cell-to-cell spread of virus in addition to its role in virus replication.

Herpes simplex virus (HSV) pUL34 plays a critical role in virus replication by mediating egress of nucleocapsids from the infected cell nucleus. We have identified a mutation in pUL34 (Y68A) that produces a major defect in virus replication and impaired nuclear egress but also profoundly inhibits cell-to-cell spread and trafficking of gE. Virion release to the extracellular medium is not affected by the Y68A mutation, indicating that the mutation specifically inhibits cell-to-cell spread. We isolated extragenic suppressors of the Y68A plaque formation defect and mapped them by a combination of high-throughput Illumina sequencing and PCR-based screening. We found that suppression is highly correlated with a nonsense mutation in the US9 gene, which plays a critical role in cell-to-cell spread of HSV-1 in neurons. The US9 mutation alone is not sufficient to suppress the Y68A spread phenotype, indicating a likely role for multiple viral factors.

GCN2 adapts protein synthesis to scavenging-dependent growth.

Pancreatic cancer cells with limited access to free amino acids can grow by scavenging extracellular protein. In a murine model of pancreatic cancer, we performed a genome-wide CRISPR screen for genes required for scavenging-dependent growth. The screen identified key mediators of macropinocytosis, peripheral lysosome positioning, endosome-lysosome fusion, lysosomal protein catabolism, and translational control. The top hit was GCN2, a kinase that suppresses translation initiation upon amino acid depletion. Using isotope tracers, we show that GCN2 is not required for protein scavenging. Instead, GCN2 prevents ribosome stalling but without slowing protein synthesis; cells still use all of the limiting amino acids as they emerge from lysosomes. GCN2 also adapts gene expression to the nutrient-poor environment, reorienting protein synthesis away from ribosomes and toward lysosomal hydrolases, such as cathepsin L. GCN2, cathepsin L, and the other genes identified in the screen are potential therapeutic targets in pancreatic cancer.

Sequence variability in clinical and laboratory isolates of herpes simplex virus 1 reveals new mutations.

Herpes simplex virus 1 (HSV-1) is a well-adapted human pathogen that can invade the peripheral nervous system and persist there as a lifelong latent infection. Despite their ubiquity, only one natural isolate of HSV-1 (strain 17) has been sequenced. Using Illumina high-throughput sequencing of viral DNA, we obtained the genome sequences of both a laboratory strain (F) and a low-passage clinical isolate (H129). These data demonstrated the extent of interstrain variation across the entire genome of HSV-1 in both coding and noncoding regions. We found many amino acid differences distributed across the proteome of the new strain F sequence and the previously known strain 17, demonstrating the spectrum of variability among wild-type HSV-1 proteins. The clinical isolate, strain H129, displays a unique anterograde spread phenotype for which the causal mutations were completely unknown. We have defined the sequence differences in H129 and propose a number of potentially causal genes, including the neurovirulence protein ICP34.5 (RL1). Further studies will be required to demonstrate which change(s) is sufficient to recapitulate the spread defect of strain H129. Unexpectedly, these data also revealed a frameshift mutation in the UL13 kinase in our strain F isolate, demonstrating how deep genome sequencing can reveal the full complement of background mutations in any given strain, particularly those passaged or plaque purified in a laboratory setting. These data increase our knowledge of sequence variation in large DNA viruses and demonstrate the potential of deep sequencing to yield insight into DNA genome evolution and the variation among different pathogen isolates.

Monitoring mammalian mitochondrial translation with MitoRiboSeq.

Several essential components of the electron transport chain, the major producer of ATP in mammalian cells, are encoded in the mitochondrial genome. These 13 proteins are translated within mitochondria by 'mitoribosomes'. Defective mitochondrial translation underlies multiple inborn errors of metabolism and has been implicated in pathologies such as aging, metabolic syndrome and cancer. Here, we provide a detailed ribosome profiling protocol optimized to interrogate mitochondrial translation in mammalian cells (MitoRiboSeq), wherein mitoribosome footprints are generated with micrococcal nuclease and mitoribosomes are separated from cytosolic ribosomes and other RNAs by ultracentrifugation in a single straightforward step. We highlight critical steps during library preparation and provide a step-by-step guide to data analysis accompanied by open-source bioinformatic code. Our method outputs mitoribosome footprints at single-codon resolution. Codons with high footprint densities are sites of mitoribosome stalling. We recently applied this approach to demonstrate that defects in mitochondrial serine catabolism or in mitochondrial tRNA methylation cause stalling of mitoribosomes at specific codons. Our method can be applied to study basic mitochondrial biology or to characterize abnormalities in mitochondrial translation in patients with mitochondrial disorders.

Genome Sequence of the Virulent Model Herpes Simplex Virus 1 Strain McKrae Demonstrates the Presence of at Least Two Widely Used Variant Strains.

Herpes simplex virus 1 (HSV-1) strain McKrae was isolated in 1965 and has been utilized by many laboratories. Three HSV-1 strain McKrae stocks have been sequenced previously, revealing discrepancies in key genes. We sequenced the genome of HSV-1 strain McKrae from the laboratory of James M. Hill to better understand the genetic differences between isolates.

Conservation of cell-intrinsic immune responses in diverse nonhuman primate species.

Differences in immune responses across species can contribute to the varying permissivity of species to the same viral pathogen. Understanding how our closest evolutionary relatives, nonhuman primates (NHPs), confront pathogens and how these responses have evolved over time could shed light on host range barriers, especially for zoonotic infections. Here, we analyzed cell-intrinsic immunity of primary cells from the broadest panel of NHP species interrogated to date, including humans, great apes, and Old and New World monkeys. Our analysis of their transcriptomes after poly(I:C) transfection revealed conservation in the functional consequences of their response. In mapping reads to either the human or the species-specific genomes, we observed that with the current state of NHP annotations, the percent of reads assigned to a genetic feature was largely similar regardless of the method. Together, these data provide a baseline for the cell-intrinsic responses elicited by a potent immune stimulus across multiple NHP donors, including endangered species, and serve as a resource for refining and furthering the existing annotations of NHP genomes.

Four Key Steps Control Glycolytic Flux in Mammalian Cells.

Altered glycolysis is a hallmark of diseases including diabetes and cancer. Despite intensive study of the contributions of individual glycolytic enzymes, systems-level analyses of flux control through glycolysis remain limited. Here, we overexpress in two mammalian cell lines the individual enzymes catalyzing each of the 12 steps linking extracellular glucose to excreted lactate, and find substantial flux control at four steps: glucose import, hexokinase, phosphofructokinase, and lactate export (and not at any steps of lower glycolysis). The four flux-controlling steps are specifically upregulated by the Ras oncogene: optogenetic Ras activation rapidly induces the transcription of isozymes catalyzing these four steps and enhances glycolysis. At least one isozyme catalyzing each of these four steps is consistently elevated in human tumors. Thus, in the studied contexts, flux control in glycolysis is concentrated in four key enzymatic steps. Upregulation of these steps in tumors likely underlies the Warburg effect.