Bat pluripotent stem cells reveal unusual entanglement between host and viruses.

Bats are distinctive among mammals due to their ability to fly, use laryngeal echolocation, and tolerate viruses. However, there are currently no reliable cellular models for studying bat biology or their response to viral infections. Here, we created induced pluripotent stem cells (iPSCs) from two species of bats: the wild greater horseshoe bat (Rhinolophus ferrumequinum) and the greater mouse-eared bat (Myotis myotis). The iPSCs from both bat species showed similar characteristics and had a gene expression profile resembling that of cells attacked by viruses. They also had a high number of endogenous viral sequences, particularly retroviruses. These results suggest that bats have evolved mechanisms to tolerate a large load of viral sequences and may have a more intertwined relationship with viruses than previously thought. Further study of bat iPSCs and their differentiated progeny will provide insights into bat biology, virus host relationships, and the molecular basis of bats' special traits.

Identification of a DNA N6-Adenine Methyltransferase Complex and Its Impact on Chromatin Organization.

DNA N6-adenine methylation (6mA) has recently been described in diverse eukaryotes, spanning unicellular organisms to metazoa. Here, we report a DNA 6mA methyltransferase complex in ciliates, termed MTA1c. It consists of two MT-A70 proteins and two homeobox-like DNA-binding proteins and specifically methylates dsDNA. Disruption of the catalytic subunit, MTA1, in the ciliate Oxytricha leads to genome-wide loss of 6mA and abolishment of the consensus ApT dimethylated motif. Mutants fail to complete the sexual cycle, which normally coincides with peak MTA1 expression. We investigate the impact of 6mA on nucleosome occupancy in vitro by reconstructing complete, full-length Oxytricha chromosomes harboring 6mA in native or ectopic positions. We show that 6mA directly disfavors nucleosomes in vitro in a local, quantitative manner, independent of DNA sequence. Furthermore, the chromatin remodeler ACF can overcome this effect. Our study identifies a diverged DNA N6-adenine methyltransferase and defines the role of 6mA in chromatin organization.

The architecture of a scrambled genome reveals massive levels of genomic rearrangement during development.

Programmed DNA rearrangements in the single-celled eukaryote Oxytricha trifallax completely rewire its germline into a somatic nucleus during development. This elaborate, RNA-mediated pathway eliminates noncoding DNA sequences that interrupt gene loci and reorganizes the remaining fragments by inversions and permutations to produce functional genes. Here, we report the Oxytricha germline genome and compare it to the somatic genome to present a global view of its massive scale of genome rearrangements. The remarkably encrypted genome architecture contains >3,500 scrambled genes, as well as >800 predicted germline-limited genes expressed, and some posttranslationally modified, during genome rearrangements. Gene segments for different somatic loci often interweave with each other. Single gene segments can contribute to multiple, distinct somatic loci. Terminal precursor segments from neighboring somatic loci map extremely close to each other, often overlapping. This genome assembly provides a draft of a scrambled genome and a powerful model for studies of genome rearrangement.

Dietary fructose improves intestinal cell survival and nutrient absorption.

Fructose consumption is linked to the rising incidence of obesity and cancer, which are two of the leading causes of morbidity and mortality globally1,2. Dietary fructose metabolism begins at the epithelium of the small intestine, where fructose is transported by glucose transporter type 5 (GLUT5; encoded by SLC2A5) and phosphorylated by ketohexokinase to form fructose 1-phosphate, which accumulates to high levels in the cell3,4. Although this pathway has been implicated in obesity and tumour promotion, the exact mechanism that drives these pathologies in the intestine remains unclear. Here we show that dietary fructose improves the survival of intestinal cells and increases intestinal villus length in several mouse models. The increase in villus length expands the surface area of the gut and increases nutrient absorption and adiposity in mice that are fed a high-fat diet. In hypoxic intestinal cells, fructose 1-phosphate inhibits the M2 isoform of pyruvate kinase to promote cell survival5-7. Genetic ablation of ketohexokinase or stimulation of pyruvate kinase prevents villus elongation and abolishes the nutrient absorption and tumour growth that are induced by feeding mice with high-fructose corn syrup. The ability of fructose to promote cell survival through an allosteric metabolite thus provides additional insights into the excess adiposity generated by a Western diet, and a compelling explanation for the promotion of tumour growth by high-fructose corn syrup.

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges.

DNA sequencing has revolutionized medicine over recent decades. However, analysis of large structural variation and repetitive DNA, a hallmark of human genomes, has been limited by short-read technology, with read lengths of 100-300 bp. Long-read sequencing (LRS) permits routine sequencing of human DNA fragments tens to hundreds of kilobase pairs in size, using both real-time sequencing by synthesis and nanopore-based direct electronic sequencing. LRS permits analysis of large structural variation and haplotypic phasing in human genomes and has enabled the discovery and characterization of rare pathogenic structural variants and repeat expansions. It has also recently enabled the assembly of a complete, gapless human genome that includes previously intractable regions, such as highly repetitive centromeres and homologous acrocentric short arms. With the addition of protocols for targeted enrichment, direct epigenetic DNA modification detection, and long-range chromatin profiling, LRS promises to launch a new era of understanding of genetic diversity and pathogenic mutations in human populations.

FLAIRR-Seq: A Method for Single-Molecule Resolution of Near Full-Length Antibody H Chain Repertoires.

Current Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using short-read sequencing strategies resolve expressed Ab transcripts with limited resolution of the C region. In this article, we present the near-full-length AIRR-seq (FLAIRR-seq) method that uses targeted amplification by 5' RACE, combined with single-molecule, real-time sequencing to generate highly accurate (99.99%) human Ab H chain transcripts. FLAIRR-seq was benchmarked by comparing H chain V (IGHV), D (IGHD), and J (IGHJ) gene usage, complementarity-determining region 3 length, and somatic hypermutation to matched datasets generated with standard 5' RACE AIRR-seq using short-read sequencing and full-length isoform sequencing. Together, these data demonstrate robust FLAIRR-seq performance using RNA samples derived from PBMCs, purified B cells, and whole blood, which recapitulated results generated by commonly used methods, while additionally resolving H chain gene features not documented in IMGT at the time of submission. FLAIRR-seq data provide, for the first time, to our knowledge, simultaneous single-molecule characterization of IGHV, IGHD, IGHJ, and IGHC region genes and alleles, allele-resolved subisotype definition, and high-resolution identification of class switch recombination within a clonal lineage. In conjunction with genomic sequencing and genotyping of IGHC genes, FLAIRR-seq of the IgM and IgG repertoires from 10 individuals resulted in the identification of 32 unique IGHC alleles, 28 (87%) of which were previously uncharacterized. Together, these data demonstrate the capabilities of FLAIRR-seq to characterize IGHV, IGHD, IGHJ, and IGHC gene diversity for the most comprehensive view of bulk-expressed Ab repertoires to date.

Hedgehog-interacting protein acts in the habenula to regulate nicotine intake.

Hedgehog-interacting protein (HHIP) sequesters Hedgehog ligands to repress Smoothened (SMO)-mediated recruitment of the GLI family of transcription factors. Allelic variation in HHIP confers risk of chronic obstructive pulmonary disease and other smoking-related lung diseases, but underlying mechanisms are unclear. Using single-cell and cell-type-specific translational profiling, we show that HHIP expression is highly enriched in medial habenula (MHb) neurons, particularly MHb cholinergic neurons that regulate aversive behavioral responses to nicotine. HHIP deficiency dysregulated the expression of genes involved in cholinergic signaling in the MHb and disrupted the function of nicotinic acetylcholine receptors (nAChRs) through a PTCH-1/cholesterol-dependent mechanism. Further, CRISPR/Cas9-mediated genomic cleavage of the Hhip gene in MHb neurons enhanced the motivational properties of nicotine in mice. These findings suggest that HHIP influences vulnerability to smoking-related lung diseases in part by regulating the actions of nicotine on habenular aversion circuits.

Genome organization and DNA accessibility control antigenic variation in trypanosomes.

Many evolutionarily distant pathogenic organisms have evolved similar survival strategies to evade the immune responses of their hosts. These include antigenic variation, through which an infecting organism prevents clearance by periodically altering the identity of proteins that are visible to the immune system of the host1. Antigenic variation requires large reservoirs of immunologically diverse antigen genes, which are often generated through homologous recombination, as well as mechanisms to ensure the expression of one or very few antigens at any given time. Both homologous recombination and gene expression are affected by three-dimensional genome architecture and local DNA accessibility2,3. Factors that link three-dimensional genome architecture, local chromatin conformation and antigenic variation have, to our knowledge, not yet been identified in any organism. One of the major obstacles to studying the role of genome architecture in antigenic variation has been the highly repetitive nature and heterozygosity of antigen-gene arrays, which has precluded complete genome assembly in many pathogens. Here we report the de novo haplotype-specific assembly and scaffolding of the long antigen-gene arrays of the model protozoan parasite Trypanosoma brucei, using long-read sequencing technology and conserved features of chromosome folding4. Genome-wide chromosome conformation capture (Hi-C) reveals a distinct partitioning of the genome, with antigen-encoding subtelomeric regions that are folded into distinct, highly compact compartments. In addition, we performed a range of analyses-Hi-C, fluorescence in situ hybridization, assays for transposase-accessible chromatin using sequencing and single-cell RNA sequencing-that showed that deletion of the histone variants H3.V and H4.V increases antigen-gene clustering, DNA accessibility across sites of antigen expression and switching of the expressed antigen isoform, via homologous recombination. Our analyses identify histone variants as a molecular link between global genome architecture, local chromatin conformation and antigenic variation.

Capturing pairwise and multi-way chromosomal conformations using chromosomal walks.

Chromosomes are folded into highly compacted structures to accommodate physical constraints within nuclei and to regulate access to genomic information. Recently, global mapping of pairwise contacts showed that loops anchoring topological domains (TADs) are highly conserved between cell types and species. Whether pairwise loops synergize to form higher-order structures is still unclear. Here we develop a conformation capture assay to study higher-order organization using chromosomal walks (C-walks) that link multiple genomic loci together into proximity chains in human and mouse cells. This approach captures chromosomal structure at varying scales. Inter-chromosomal contacts constitute only 7-10% of the pairs and are restricted by interfacing TADs. About half of the C-walks stay within one chromosome, and almost half of those are restricted to intra-TAD spaces. C-walks that couple 2-4 TADs indicate stochastic associations between transcriptionally active, early replicating loci. Targeted analysis of thousands of 3-walks anchored at highly expressed genes support pairwise, rather than hub-like, chromosomal topology at active loci. Polycomb-repressed Hox domains are shown by the same approach to enrich for synergistic hubs. Together, the data indicate that chromosomal territories, TADs, and intra-TAD loops are primarily driven by nested, possibly dynamic, pairwise contacts.

Elucidation of de novo small insertion/deletion biology with parent-of-origin phasing.

The mechanisms underlying de novo insertion/deletion (indel) genesis, such as polymerase slippage, have been hypothesized but not well characterized in the human genome. We implemented two methodological improvements, which were leveraged to dissect indel mutagenesis. We assigned de novo variants to parent-of-origin (i.e., phasing) with low-coverage long-read whole-genome sequencing, achieving better phasing compared to short-read sequencing (medians of 84% and 23%, respectively). We then wrote an application programming interface to classify indels into three subtypes according to sequence context. Across three cohorts with different phasing methods (Ntrios = 540, all cohorts), we observed that one de novo indel subtype, change in copy count (CCC), was significantly correlated with father's (p = 7.1 × 10-4 ) but not mother's (p = .45) age at conception. We replicated this effect in three cohorts without de novo phasing (ppaternal = 1.9 × 10-9 , pmaternal = .61; Ntrios = 3,391, all cohorts). Although this is consistent with polymerase slippage during spermatogenesis, the percentage of variance explained by paternal age was low, and we did not observe an association with replication timing. These results suggest that spermatogenesis-specific events have a minor role in CCC indel mutagenesis, one not observed for other indel subtypes nor for maternal age in general. These results have implications for indel modeling in evolution and disease.

Microevolution in response to transient heme-iron restriction enhances intracellular bacterial community development and persistence.

Bacterial pathogens must sense, respond and adapt to a myriad of dynamic microenvironmental stressors to survive. Adaptation is key for colonization and long-term ability to endure fluctuations in nutrient availability and inflammatory processes. We hypothesize that strains adapted to survive nutrient deprivation are more adept for colonization and establishment of chronic infection. In this study, we detected microevolution in response to transient nutrient limitation through mutation of icc. The mutation results in decreased 3',5'-cyclic adenosine monophosphate phosphodiesterase activity in nontypeable Haemophilus influenzae (NTHI). In a preclinical model of NTHI-induced otitis media (OM), we observed a significant decrease in the recovery of effusion from ears infected with the icc mutant strain. Clinically, resolution of OM coincides with the clearance of middle ear fluid. In contrast to this clinical paradigm, we observed that the icc mutant strain formed significantly more intracellular bacterial communities (IBCs) than the parental strain early during experimental OM. Although the number of IBCs formed by the parental strain was low at early stages of OM, we observed a significant increase at later stages that coincided with absence of recoverable effusion, suggesting the presence of a mucosal reservoir following resolution of clinical disease. These data provide the first insight into NTHI microevolution during nutritional limitation and provide the first demonstration of IBCs in a preclinical model of chronic OM.

Capture of complete ciliate chromosomes in single sequencing reads reveals widespread chromosome isoforms.

BACKGROUND: Whole-genome shotgun sequencing, which stitches together millions of short sequencing reads into a single genome, ushered in the era of modern genomics and led to a rapid expansion of the number of genome sequences available. Nevertheless, assembly of short reads remains difficult, resulting in fragmented genome sequences. Ultimately, only a sequencing technology capable of capturing complete chromosomes in a single run could resolve all ambiguities. Even "third generation" sequencing technologies produce reads far shorter than most eukaryotic chromosomes. However, the ciliate Oxytricha trifallax has a somatic genome with thousands of chromosomes averaging only 3.2 kbp, making it an ideal candidate for exploring the benefits of sequencing whole chromosomes without assembly. RESULTS: We used single-molecule real-time sequencing to capture thousands of complete chromosomes in single reads and to update the published Oxytricha trifallax JRB310 genome assembly. In this version, over 50% of the completed chromosomes with two telomeres derive from single reads. The improved assembly includes over 12,000 new chromosome isoforms, and demonstrates that somatic chromosomes derive from variable rearrangements between somatic segments encoded up to 191,000 base pairs away. However, while long reads reduce the need for assembly, a hybrid approach that supplements long-read sequencing with short reads for error correction produced the most complete and accurate assembly, overall. CONCLUSIONS: This assembly provides the first example of complete eukaryotic chromosomes captured by single sequencing reads and demonstrates that traditional approaches to genome assembly can mask considerable structural variation.

Identification of Differentially Expressed Splice Variants by the Proteogenomic Pipeline Splicify.

Proteogenomics, i.e. comprehensive integration of genomics and proteomics data, is a powerful approach identifying novel protein biomarkers. This is especially the case for proteins that differ structurally between disease and control conditions. As tumor development is associated with aberrant splicing, we focus on this rich source of cancer specific biomarkers. To this end, we developed a proteogenomic pipeline, Splicify, which can detect differentially expressed protein isoforms. Splicify is based on integrating RNA massive parallel sequencing data and tandem mass spectrometry proteomics data to identify protein isoforms resulting from differential splicing between two conditions. Proof of concept was obtained by applying Splicify to RNA sequencing and mass spectrometry data obtained from colorectal cancer cell line SW480, before and after siRNA-mediated downmodulation of the splicing factors SF3B1 and SRSF1. These analyses revealed 2172 and 149 differentially expressed isoforms, respectively, with peptide confirmation upon knock-down of SF3B1 and SRSF1 compared with their controls. Splice variants identified included RAC1, OSBPL3, MKI67, and SYK. One additional sample was analyzed by PacBio Iso-Seq full-length transcript sequencing after SF3B1 downmodulation. This analysis verified the alternative splicing identified by Splicify and in addition identified novel splicing events that were not represented in the human reference genome annotation. Therefore, Splicify offers a validated proteogenomic data analysis pipeline for identification of disease specific protein biomarkers resulting from mRNA alternative splicing. Splicify is publicly available on GitHub (https://github.com/NKI-TGO/SPLICIFY) and suitable to address basic research questions using pre-clinical model systems as well as translational research questions using patient-derived samples, e.g. allowing to identify clinically relevant biomarkers.

Parallel Evolution of Two Clades of an Atlantic-Endemic Pathogenic Lineage of Vibrio parahaemolyticus by Independent Acquisition of Related Pathogenicity Islands.

Shellfish-transmitted Vibrio parahaemolyticus infections have recently increased from locations with historically low disease incidence, such as the Northeast United States. This change coincided with a bacterial population shift toward human-pathogenic variants occurring in part through the introduction of several Pacific native lineages (ST36, ST43, and ST636) to nearshore areas off the Atlantic coast of the Northeast United States. Concomitantly, ST631 emerged as a major endemic pathogen. Phylogenetic trees of clinical and environmental isolates indicated that two clades diverged from a common ST631 ancestor, and in each of these clades, a human-pathogenic variant evolved independently through acquisition of distinct Vibrio pathogenicity islands (VPaI). These VPaI differ from each other and bear little resemblance to hemolysin-containing VPaI from isolates of the pandemic clonal complex. Clade I ST631 isolates either harbored no hemolysins or contained a chromosome I-inserted island we call VPaIß that encodes a type 3 secretion system (T3SS2ß) typical of Trh hemolysin producers. The more clinically prevalent and clonal ST631 clade II had an island we call VPaIγ that encodes both tdh and trh and that was inserted in chromosome II. VPaIγ was derived from VPaIß but with some additional acquired elements in common with VPaI carried by pandemic isolates, exemplifying the mosaic nature of pathogenicity islands. Genomics comparisons and amplicon assays identified VPaIγ-type islands containing tdh inserted adjacent to the ure cluster in the three introduced Pacific and most other emergent lineages that collectively cause 67% of infections in the Northeast United States as of 2016.IMPORTANCE The availability of three different hemolysin genotypes in the ST631 lineage provided a unique opportunity to employ genome comparisons to further our understanding of the processes underlying pathogen evolution. The fact that two different pathogenic clades arose in parallel from the same potentially benign lineage by independent VPaI acquisition is surprising considering the historically low prevalence of community members harboring VPaI in waters along the Northeast U.S. coast that could serve as the source of this material. This illustrates a possible predisposition of some lineages to not only acquire foreign DNA but also become human pathogens. Whereas the underlying cause for the expansion of V. parahaemolyticus lineages harboring VPaIγ along the U.S. Atlantic coast and spread of this element to multiple lineages that underlies disease emergence is not known, this work underscores the need to define the environment factors that favor bacteria harboring VPaI in locations of emergent disease.

First Report of blaIMP-14 on a Plasmid Harboring Multiple Drug Resistance Genes in Escherichia coli Sequence Type 131.

The blaIMP-14 carbapenem resistance gene has largely previously been observed in Pseudomonas aeruginosa and Acinetobacter spp. As part of global surveillance and sequencing of carbapenem-resistant Escherichia coli, we identified a sequence type 131 strain harboring blaIMP-14 within a class 1 integron, itself nested within an ∼54-kb multidrug resistance region on an epidemic IncA/C2 plasmid. The emergence of blaIMP-14 in this context in the ST131 lineage is of potential clinical concern.

Complete genome sequence of Vibrio parahaemolyticus ST36 strain MAVP-26, a clinical isolate from an oyster-borne human gastric infection.

A Pacific native lineage of Vibrio parahaemolyticus ST36 serotype O4:K12 was introduced into the Atlantic, which increased local source illnesses. To identify genetic determinants of virulence and ecological resiliency and track their transfer into endemic populations, we constructed a complete genome of a 2013 Atlantic-traced clinical isolate by hybrid assembly.

A PRMT5-ZNF326 axis mediates innate immune activation upon replication stress.

DNA replication stress (RS) is a widespread phenomenon in carcinogenesis, causing genomic instability and extensive chromatin alterations. DNA damage leads to activation of innate immune signaling, but little is known about transcriptional regulators mediating such signaling upon RS. Using a chemical screen, we identified protein arginine methyltransferase 5 (PRMT5) as a key mediator of RS-dependent induction of interferon-stimulated genes (ISGs). This response is also associated with reactivation of endogenous retroviruses (ERVs). Using quantitative mass spectrometry, we identify proteins with PRMT5-dependent symmetric dimethylarginine (SDMA) modification induced upon RS. Among these, we show that PRMT5 targets and modulates the activity of ZNF326, a zinc finger protein essential for ISG response. Our data demonstrate a role for PRMT5-mediated SDMA in the context of RS-induced transcriptional induction, affecting physiological homeostasis and cancer therapy.

A master regulator of opioid reward in the ventral prefrontal cortex.

In addition to their intrinsic rewarding properties, opioids can also evoke aversive reactions that protect against misuse. Cellular mechanisms that govern the interplay between opioid reward and aversion are poorly understood. We used whole-brain activity mapping in mice to show that neurons in the dorsal peduncular nucleus (DPn) are highly responsive to the opioid oxycodone. Connectomic profiling revealed that DPn neurons innervate the parabrachial nucleus (PBn). Spatial and single-nuclei transcriptomics resolved a population of PBn-projecting pyramidal neurons in the DPn that express µ-opioid receptors (µORs). Disrupting µOR signaling in the DPn switched oxycodone from rewarding to aversive and exacerbated the severity of opioid withdrawal. These findings identify the DPn as a key substrate for the abuse liability of opioids.

Single cell view of tumor microenvironment gradients in pleural mesothelioma.

Immunotherapies have shown great promise in pleural mesothelioma (PM), yet most patients still do not achieve significant clinical response, highlighting the importance of improving understanding of the tumor microenvironment (TME). Here, we utilized high-throughput, single-cell RNA-sequencing to de novo identify 54 expression programs and construct a comprehensive cellular catalogue of the PM TME. We found four cancer-intrinsic programs associated with poor disease outcome and a novel fetal-like, endothelial cell population that likely responds to VEGF signaling and promotes angiogenesis. Throughout cellular compartments, we observe substantial difference in the TME associated with a cancer-intrinsic sarcomatoid signature, including enrichment in fetal-like endothelial cells, CXCL9+ macrophages, cytotoxic, exhausted, and regulatory T cells, which we validated using imaging and bulk deconvolution analyses on independent cohorts. Finally, we show, both computationally and experimentally, that NKG2A-HLA-E interaction between NK and tumor cells represents an important new therapeutic axis in PM, especially for epithelioid cases.