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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Selective expansion of myeloid and NK cells in humanized mice yields human-like vaccine responses.

Mice engrafted with components of a human immune system have become widely-used models for studying aspects of human immunity and disease. However, a defined methodology to objectively measure and compare the quality of the human immune response in different models is lacking. Here, by taking advantage of the highly immunogenic live-attenuated yellow fever virus vaccine YFV-17D, we provide an in-depth comparison of immune responses in human vaccinees, conventional humanized mice, and second generation humanized mice. We demonstrate that selective expansion of human myeloid and natural killer cells promotes transcriptomic responses akin to those of human vaccinees. These enhanced transcriptomic profiles correlate with the development of an antigen-specific cellular and humoral response to YFV-17D. Altogether, our approach provides a robust scoring of the quality of the human immune response in humanized mice and highlights a rational path towards developing better pre-clinical models for studying the human immune response and disease.

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.

Whole-Genome Sequence and Variant Analysis of W303, a Widely-Used Strain of Saccharomyces cerevisiae.

The yeast Saccharomyces cerevisiae has emerged as a superior model organism. Selection of distinct laboratory strains of S. cerevisiae with unique phenotypic properties, such as superior mating or sporulation efficiencies, has facilitated advancements in research. W303 is one such laboratory strain that is closely related to the first completely sequenced yeast strain, S288C. In this work, we provide a high-quality, annotated genome sequence for W303 for utilization in comparative analyses and genome-wide studies. Approximately 9500 variations exist between S288C and W303, affecting the protein sequences of ∼700 genes. A listing of the polymorphisms and divergent genes is provided for researchers interested in identifying the genetic basis for phenotypic differences between W303 and S288C. Several divergent functional gene families were identified, including flocculation and sporulation genes, likely representing selection for desirable laboratory phenotypes. Interestingly, remnants of ancestor wine strains were found on several chromosomes. Finally, as a test of the utility of the high-quality reference genome, variant mapping revealed more accurate identification of accumulated mutations in passaged mismatch repair-defective strains.

Genome-wide Purification of Extrachromosomal Circular DNA from Eukaryotic Cells.

Extrachromosomal circular DNAs (eccDNAs) are common genetic elements in Saccharomyces cerevisiae and are reported in other eukaryotes as well. EccDNAs contribute to genetic variation among somatic cells in multicellular organisms and to evolution of unicellular eukaryotes. Sensitive methods for detecting eccDNA are needed to clarify how these elements affect genome stability and how environmental and biological factors induce their formation in eukaryotic cells. This video presents a sensitive eccDNA-purification method called Circle-Seq. The method encompasses column purification of circular DNA, removal of remaining linear chromosomal DNA, rolling-circle amplification of eccDNA, deep sequencing, and mapping. Extensive exonuclease treatment was required for sufficient linear chromosomal DNA degradation. The rolling-circle amplification step by φ29 polymerase enriched for circular DNA over linear DNA. Validation of the Circle-Seq method on three S. cerevisiae CEN.PK populations of 10(10) cells detected hundreds of eccDNA profiles in sizes larger than 1 kilobase. Repeated findings of ASP3-1, COS111, CUP1, RSC30, HXT6, HXT7 genes on circular DNA in both S288c and CEN.PK suggests that DNA circularization is conserved between strains at these loci. In sum, the Circle-Seq method has broad applicability for genome-scale screening for eccDNA in eukaryotes as well as for detecting specific eccDNA types.

Formation of Extrachromosomal Circular DNA from Long Terminal Repeats of Retrotransposons in Saccharomyces cerevisiae.

Extrachromosomal circular DNA (eccDNA) derived from chromosomal Ty retrotransposons in yeast can be generated in multiple ways. Ty eccDNA can arise from the circularization of extrachromosomal linear DNA during the transpositional life cycle of retrotransposons, or from circularization of genomic Ty DNA. Circularization may happen through nonhomologous end-joining (NHEJ) of long terminal repeats (LTRs) flanking Ty elements, by Ty autointegration, or by LTR-LTR recombination. By performing an in-depth investigation of sequence reads stemming from Ty eccDNAs obtained from populations of Saccharomyces cerevisiae S288c, we find that eccDNAs predominantly correspond to full-length Ty1 elements. Analyses of sequence junctions reveal no signs of NHEJ or autointegration events. We detect recombination junctions that are consistent with yeast Ty eccDNAs being generated through recombination events within the genome. This opens the possibility that retrotransposable elements could move around in the genome without an RNA intermediate directly through DNA circularization.

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.

Rapid genome assembly and comparison decode intrastrain variation in human alphaherpesviruses.

UNLABELLED: Herpes simplex virus (HSV) is a widespread pathogen that causes epithelial lesions with recurrent disease that manifests over a lifetime. The lifelong aspect of infection results from latent viral infection of neurons, a reservoir from which the virus reactivates periodically. Recent work has demonstrated the breadth of genetic variation in globally distributed HSV strains. However, the amount of variation or capacity for mutation within one strain has not been well studied. Here we developed and applied a streamlined new approach for assembly and comparison of large DNA viral genomes such as HSV-1. This viral genome assembly (VirGA) workflow incorporates a combination of de novo assembly, alignment, and annotation strategies to automate the generation of draft genomes for large viruses. We applied this approach to quantify the amount of variation between clonal derivatives of a common parental virus stock. In addition, we examined the genetic basis for syncytial plaque phenotypes displayed by a subset of these strains. In each of the syncytial strains, we found an identical DNA change, affecting one residue in the gB (UL27) fusion protein. Since these identical mutations could have appeared after extensive in vitro passaging, we applied the VirGA sequencing and comparison approach to two clinical HSV-1 strains isolated from the same patient. One of these strains was syncytial upon first culturing; its sequence revealed the same gB mutation. These data provide insight into the extent and origin of genome-wide intrastrain HSV-1 variation and present useful methods for expansion to in vivo patient infection studies. IMPORTANCE: Herpes simplex virus (HSV) infects more than 70% of adults worldwide, causing epithelial lesions and recurrent disease that manifests over a lifetime. Prior work has demonstrated that HSV strains vary from country to country and between individuals. However, the amount of variation within one strain has not been well studied. To address this, we developed a new approach for viral genome assembly (VirGA) and analysis. We used this approach to quantify the amount of variation between sister clones of a common parental virus stock and to determine the basis of a unique fusion phenotype displayed by several variants. These data revealed that while sister clones of one HSV stock are more than 98% identical, these variants harbor enough genetic differences to change their observed characteristics. Comparative genomics approaches will allow us to explore the impacts of viral inter- and intrastrain diversity on drug and vaccine efficacy.

Genome Sequence of the Anterograde-Spread-Defective Herpes Simplex Virus 1 Strain MacIntyre.

We used paired-end Illumina deep sequencing and de novo assembly to determine the genome sequence of herpes simplex virus 1 (HSV-1) strain MacIntyre (aka McIntyre). The MacIntyre strain originated from the brain of a patient with lethal HSV encephalitis and has a unique limitation in its neuronal spread, moving solely in the retrograde direction.