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.

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.

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.

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.

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.