Smoother: on-the-fly processing of interactome data using prefix sums.

Nucleic acid interactome data, such as chromosome conformation capture data and RNA-DNA interactome data, are currently analyzed via pipelines that must be rerun for each new parameter set. A more dynamic approach is desirable since the optimal parameter set is commonly unknown ahead of time and rerunning pipelines is a time-consuming process. We have developed an approach fast enough to process interactome data on-the-fly using a sparse prefix sum index. With this index, we created Smoother, a flexible, multifeatured visualization and analysis tool that allows interactive filtering, e.g. by mapping quality, almost instant comparisons between different normalization approaches, e.g. iterative correction, and ploidy correction. Further, Smoother can overlay other sequencing data or genomic annotations, compare different samples, and perform virtual 4C analysis. Smoother permits a novel way to interact with and explore interactome data, fostering comprehensive, high-quality data analysis. Smoother is available at https://github.com/Siegel-Lab/BioSmoother under the MIT license.

Decoding the impact of nuclear organization on antigenic variation in parasites.

Antigenic variation as a strategy to evade the host adaptive immune response has evolved in divergent pathogens. Antigenic variation involves restricted, and often mutually exclusive, expression of dominant antigens and a periodic switch in antigen expression during infection. In eukaryotes, nuclear compartmentalization, including three-dimensional folding of the genome and physical separation of proteins in compartments or condensates, regulates mutually exclusive gene expression and chromosomal translocations. In this Review, we discuss the impact of nuclear organization on antigenic variation in the protozoan pathogens Trypanosoma brucei and Plasmodium falciparum. In particular, we highlight the relevance of nuclear organization in both mutually exclusive antigen expression and genome stability, which underlie antigenic variation.

H2B.V demarcates divergent strand-switch regions, some tDNA loci, and genome compartments in Trypanosoma cruzi and affects parasite differentiation and host cell invasion.

Histone variants play a crucial role in chromatin structure organization and gene expression. Trypanosomatids have an unusual H2B variant (H2B.V) that is known to dimerize with the variant H2A.Z generating unstable nucleosomes. Previously, we found that H2B.V protein is enriched in tissue-derived trypomastigote (TCT) life forms, a nonreplicative stage of Trypanosoma cruzi, suggesting that this variant may contribute to the differences in chromatin structure and global transcription rates observed among parasite life forms. Here, we performed the first genome-wide profiling of histone localization in T. cruzi using epimastigotes and TCT life forms, and we found that H2B.V was preferentially located at the edges of divergent transcriptional strand switch regions, which encompass putative transcriptional start regions; at some tDNA loci; and between the conserved and disrupted genome compartments, mainly at trans-sialidase, mucin and MASP genes. Remarkably, the chromatin of TCT forms was depleted of H2B.V-enriched peaks in comparison to epimastigote forms. Interactome assays indicated that H2B.V associated specifically with H2A.Z, bromodomain factor 2, nucleolar proteins and a histone chaperone, among others. Parasites expressing reduced H2B.V levels were associated with higher rates of parasite differentiation and mammalian cell infectivity. Taken together, H2B.V demarcates critical genomic regions and associates with regulatory chromatin proteins, suggesting a scenario wherein local chromatin structures associated with parasite differentiation and invasion are regulated during the parasite life cycle.

Allele-specific assembly of a eukaryotic genome corrects apparent frameshifts and reveals a lack of nonsense-mediated mRNA decay.

To date, most reference genomes represent a mosaic consensus sequence in which the homologous chromosomes are collapsed into one sequence. This approach produces sequence artefacts and impedes analyses of allele-specific mechanisms. Here, we report an allele-specific genome assembly of the diploid parasite Trypanosoma brucei and reveal allelic variants affecting gene expression. Using long-read sequencing and chromosome conformation capture data, we could assign 99.5% of all heterozygote variants to a specific homologous chromosome and build a 66 Mb long allele-specific genome assembly. The phasing of haplotypes allowed us to resolve hundreds of artefacts present in the previous mosaic consensus assembly. In addition, it revealed allelic recombination events, visible as regions of low allelic heterozygosity, enabling the lineage tracing of T. brucei isolates. Interestingly, analyses of transcriptome and translatome data of genes with allele-specific premature termination codons point to the absence of a nonsense-mediated decay mechanism in trypanosomes. Taken together, this study delivers a reference quality allele-specific genome assembly of T. brucei and demonstrates the importance of such assemblies for the study of gene expression control. We expect the new genome assembly will increase the awareness of allele-specific phenomena and provide a platform to investigate them.

Cell-to-Cell Heterogeneity in Trypanosomes.

Recent developments in single-cell and single-molecule techniques have revealed surprising levels of heterogeneity among isogenic cells. These advances have transformed the study of cell-to-cell heterogeneity into a major area of biomedical research, revealing that it can confer essential advantages, such as priming populations of unicellular organisms for future environmental stresses. Protozoan parasites, such as trypanosomes, face multiple and often hostile environments, and to survive, they undergo multiple changes, including changes in morphology, gene expression, and metabolism. But why does only a subset of proliferative cells differentiate to the next life cycle stage? Why do only some bloodstream parasites undergo antigenic switching while others stably express one variant surface glycoprotein? And why do some parasites invade an organ while others remain in the bloodstream? Building on extensive research performed in bacteria, here we suggest that biological noise can contribute to the fitness of eukaryotic pathogens and discuss the importance of cell-to-cell heterogeneity in trypanosome infections.

Spatial integration of transcription and splicing in a dedicated compartment sustains monogenic antigen expression in African trypanosomes.

Highly selective gene expression is a key requirement for antigenic variation in several pathogens, allowing evasion of host immune responses and maintenance of persistent infections1. African trypanosomes-parasites that cause lethal diseases in humans and livestock-employ an antigenic variation mechanism that involves monogenic antigen expression from a pool of >2,600 antigen-coding genes2. In other eukaryotes, the expression of individual genes can be enhanced by mechanisms involving the juxtaposition of otherwise distal chromosomal loci in the three-dimensional nuclear space3-5. However, trypanosomes lack classical enhancer sequences or regulated transcription initiation6,7. In this context, it has remained unclear how genome architecture contributes to monogenic transcription elongation and transcript processing. Here, we show that the single expressed antigen-coding gene displays a specific inter-chromosomal interaction with a major messenger RNA splicing locus. Chromosome conformation capture (Hi-C) revealed a dynamic reconfiguration of this inter-chromosomal interaction upon activation of another antigen. Super-resolution microscopy showed the interaction to be heritable and splicing dependent. We found a specific association of the two genomic loci with the antigen exclusion complex, whereby VSG exclusion 1 (VEX1) occupied the splicing locus and VEX2 occupied the antigen-coding locus. Following VEX2 depletion, loss of monogenic antigen expression was accompanied by increased interactions between previously silent antigen genes and the splicing locus. Our results reveal a mechanism to ensure monogenic expression, where antigen transcription and messenger RNA splicing occur in a specific nuclear compartment. These findings suggest a new means of post-transcriptional gene regulation.

Distinct roles for H4 and H2A.Z acetylation in RNA transcription in African trypanosomes.

Despite histone H2A variants and acetylation of histones occurring in almost every eukaryotic organism, it has been difficult to establish direct functional links between canonical histones or H2A variant acetylation, deposition of H2A variants and transcription. To disentangle these complex interdependent processes, we devised a highly sensitive strategy for quantifying histone acetylation levels at specific genomic loci. Taking advantage of the unusual genome organization in Trypanosoma brucei, we identified 58 histone modifications enriched at transcription start sites (TSSs). Furthermore, we found TSS-associated H4 and H2A.Z acetylation to be mediated by two different histone acetyltransferases, HAT2 and HAT1, respectively. Whereas depletion of HAT2 decreases H2A.Z deposition and shifts the site of transcription initiation, depletion of HAT1 does not affect H2A.Z deposition but reduces total mRNA levels by 50%. Thus, specifically reducing H4 or H2A.Z acetylation levels enabled us to reveal distinct roles for these modifications in H2A.Z deposition and RNA transcription.

Photoemission from non-polar aromatic molecules in the gas and liquid phase.

The photoelectron spectra of both liquid and gas phase aromatic molecules are reported. The spectra were obtained using a 34.1 eV source produced by high harmonic generation and analysed with the help of high-level ab initio simulations using the reflection principle combined with path integral molecular dynamics simulations accounting for nuclear quantum effects for the gas phase. We demonstrate the suitability of three trimethylbenzenes (1,3,5-trimethylbenzene, 1,2,3-trimethylbenzene and 1,2,4-trimethylbenzene) as a solvent for liquid photoelectron spectroscopy of solute species. We also discuss the electrokinetic charging of a non-polar liquid jet.

Efficient and specific oligo-based depletion of rRNA.

In most organisms, ribosomal RNA (rRNA) contributes to >85% of total RNA. Thus, to obtain useful information from RNA-sequencing (RNA-seq) analyses at reasonable sequencing depth, typically, mature polyadenylated transcripts are enriched or rRNA molecules are depleted. Targeted depletion of rRNA is particularly useful when studying transcripts lacking a poly(A) tail, such as some non-coding RNAs (ncRNAs), most bacterial RNAs and partially degraded or immature transcripts. While several commercially available kits allow effective rRNA depletion, their efficiency relies on a high degree of sequence homology between oligonucleotide probes and the target RNA. This restricts the use of such kits to a limited number of organisms with conserved rRNA sequences. In this study we describe the use of biotinylated oligos and streptavidin-coated paramagnetic beads for the efficient and specific depletion of trypanosomal rRNA. Our approach reduces the levels of the most abundant rRNA transcripts to less than 5% with minimal off-target effects. By adjusting the sequence of the oligonucleotide probes, our approach can be used to deplete rRNAs or other abundant transcripts independent of species. Thus, our protocol provides a useful alternative for rRNA removal where enrichment of polyadenylated transcripts is not an option and commercial kits for rRNA are not available.

Rapid Identification of Ischemic Injury in Renal Tissue by Mass-Spectrometry Imaging.

The increasing analytical speed of mass-spectrometry imaging (MSI) has led to growing interest in the medical field. Acute kidney injury is a severe disease with high morbidity and mortality. No reliable cut-offs are known to estimate the severity of acute kidney injury. Thus, there is a need for new tools to rapidly and accurately assess acute ischemia, which is of clinical importance in intensive care and in kidney transplantation. We investigated the value of MSI to assess acute ischemic kidney tissue in a porcine model. A perfusion model was developed where paired kidneys received warm (severe) or cold (minor) ischemia ( n = 8 per group). First, ischemic tissue damage was systematically assessed by two blinded pathologists. Second, MALDI-MSI of kidney tissues was performed to study the spatial distributions and compositions of lipids in the tissues. Histopathological examination revealed no significant difference between kidneys, whereas MALDI-MSI was capable of a detailed discrimination of severe and mild ischemia by differential expression of characteristic lipid-degradation products throughout the tissue within 2 h. In particular, lysolipids, including lysocardiolipins, lysophosphatidylcholines, and lysophosphatidylinositol, were dramatically elevated after severe ischemia. This study demonstrates the significant potential of MSI to differentiate and identify molecular patterns of early ischemic injury in a clinically acceptable time frame. The observed changes highlight the underlying biochemical processes of acute ischemic kidney injury and provide a molecular classification tool that can be deployed in assessment of acute ischemic kidney injury.

Ribosome Profiling Reveals HSP90 Inhibitor Effects on Stage-Specific Protein Synthesis in Leishmania donovani.

The 90-kDa heat shock protein (HSP90) of eukaryotes is a highly abundant and essential chaperone required for the maturation of regulatory and signal proteins. In the protozoan parasite Leishmania donovani, causative agent of the fatal visceral leishmaniasis, HSP90 activity is essential for cell proliferation and survival. Even more importantly, its inhibition causes life cycle progression from the insect stage to the pathogenic, mammalian stage. To unravel the molecular impact of HSP90 activity on the parasites' gene expression, we performed a ribosome profiling analysis of L. donovani, comparing genome-wide protein synthesis patterns in the presence and absence of the HSP90-specific inhibitor radicicol and an ectopically expressed radicicol-resistant HSP90 variant. We find that ribosome-protected RNA faithfully maps open reading frames and represents 97% of the annotated protein-coding genes of L. donovani. Protein synthesis was found to correlate poorly with RNA steady-state levels, indicating a regulated translation as primary mechanism for HSP90-dependent gene expression. The results confirm inhibitory effects of HSP90 on the synthesis of Leishmania proteins that are associated with the pathogenic, intracellular stage of the parasite. Those include heat shock proteins, redox enzymes, virulence-enhancing surface proteins, proteolytic pathways, and a complete set of histones. Conversely, HSP90 promotes fatty acid synthesis enzymes. Complementing radicicol treatment with the radicicol-resistant HSP90rr variant revealed important off-target radicicol effects that control a large number of the above-listed proteins. Leishmania lacks gene-specific transcription regulation and relies on regulated translation instead. Our ribosome footprinting analysis demonstrates a controlling function of HSP90 in stage-specific protein synthesis but also significant, HSP90-independent effects of the inhibitor radicicol. IMPORTANCE Leishmania parasites cause severe illness in humans and animals. They exist in two developmental stages, insect form and mammalian form, which differ in shape and gene expression. By mapping and quantifying RNA fragments protected by protein synthesis complexes, we determined the rates of protein synthesis for >90% of all Leishmania proteins in response to the inhibition of a key regulatory protein, the 90-kDa heat shock protein. We find that Leishmania depends on a regulation of protein synthesis for controlling its gene expression and that heat shock protein 90 inhibition can trigger the developmental program from insect form to mammalian form of the pathogen.

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.

Exploiting CRISPR-Cas9 technology to investigate individual histone modifications.

Despite their importance for most DNA-templated processes, the function of individual histone modifications has remained largely unknown because in vivo mutational analyses are lacking. The reason for this is that histone genes are encoded by multigene families and that tools to simultaneously edit multiple genomic loci with high efficiency are only now becoming available. To overcome these challenges, we have taken advantage of the power of CRISPR-Cas9 for precise genome editing and of the fact that most DNA repair in the protozoan parasite Trypanosoma brucei occurs via homologous recombination. By establishing an episome-based CRISPR-Cas9 system for T. brucei, we have edited wild type cells without inserting selectable markers, inserted a GFP tag between an ORF and its 3'UTR, deleted both alleles of a gene in a single transfection, and performed precise editing of genes that exist in multicopy arrays, replacing histone H4K4 with H4R4 in the absence of detectable off-target effects. The newly established genome editing toolbox allows for the generation of precise mutants without needing to change other regions of the genome, opening up opportunities to study the role of individual histone modifications, catalytic sites of enzymes or the regulatory potential of UTRs in their endogenous environments.

Comparative ribosome profiling uncovers a dominant role for translational control in Toxoplasma gondii.

BACKGROUND: The lytic cycle of the protozoan parasite Toxoplasma gondii, which involves a brief sojourn in the extracellular space, is characterized by defined transcriptional profiles. For an obligate intracellular parasite that is shielded from the cytosolic host immune factors by a parasitophorous vacuole, the brief entry into the extracellular space is likely to exert enormous stress. Due to its role in cellular stress response, we hypothesize that translational control plays an important role in regulating gene expression in Toxoplasma during the lytic cycle. Unlike transcriptional profiles, insights into genome-wide translational profiles of Toxoplasma gondii are lacking. METHODS: We have performed genome-wide ribosome profiling, coupled with high throughput RNA sequencing, in intracellular and extracellular Toxoplasma gondii parasites to investigate translational control during the lytic cycle. RESULTS: Although differences in transcript abundance were mostly mirrored at the translational level, we observed significant differences in the abundance of ribosome footprints between the two parasite stages. Furthermore, our data suggest that mRNA translation in the parasite is potentially regulated by mRNA secondary structure and upstream open reading frames. CONCLUSION: We show that most of the Toxoplasma genes that are dysregulated during the lytic cycle are translationally regulated.

GT-rich promoters can drive RNA pol II transcription and deposition of H2A.Z in African trypanosomes.

Genome-wide transcription studies are revealing an increasing number of "dispersed promoters" that, unlike "focused promoters", lack well-conserved sequence motifs and tight regulation. Dispersed promoters are nevertheless marked by well-defined chromatin structures, suggesting that specific sequence elements must exist in these unregulated promoters. Here, we have analyzed regions of transcription initiation in the eukaryotic parasite Trypanosoma brucei, in which RNA polymerase II transcription initiation occurs over broad regions without distinct promoter motifs and lacks regulation. Using a combination of site-specific and genome-wide assays, we identified GT-rich promoters that can drive transcription and promote the targeted deposition of the histone variant H2A.Z in a genomic context-dependent manner. In addition, upon mapping nucleosome occupancy at high resolution, we find nucleosome positioning to correlate with RNA pol II enrichment and gene expression, pointing to a role in RNA maturation. Nucleosome positioning may thus represent a previously unrecognized layer of gene regulation in trypanosomes. Our findings show that even highly dispersed, unregulated promoters contain specific DNA elements that are able to induce transcription and changes in chromatin structure.

Genome-wide analysis of chromatin structures in Trypanosoma brucei using high-resolution MNase-ChIP-seq.

Specific DNA-protein interactions are the basis for many important cellular mechanisms like the regulation of gene expression or replication. Knowledge about the precise genomic locations of DNA-protein interactions is important because it provides insight into the regulation of these processes. Recently, we have adapted an approach that combines micrococcal nuclease (MNase) digestion of chromatin with chromatin immunoprecipitation in Trypanosoma brucei. Here, we describe in detail how this method can be used to map the genome-wide distribution of nucleosomes or other DNA-binding proteins at high resolution in T. brucei.

Base J represses genes at the end of polycistronic gene clusters in Leishmania major by promoting RNAP II termination.

The genomes of kinetoplastids are organized into polycistronic gene clusters that are flanked by the modified DNA base J. Previous work has established a role of base J in promoting RNA polymerase II termination in Leishmania spp. where the loss of J leads to termination defects and transcription into adjacent gene clusters. It remains unclear whether these termination defects affect gene expression and whether read through transcription is detrimental to cell growth, thus explaining the essential nature of J. We now demonstrate that reduction of base J at specific sites within polycistronic gene clusters in L. major leads to read through transcription and increased expression of downstream genes in the cluster. Interestingly, subsequent transcription into the opposing polycistronic gene cluster does not lead to downregulation of sense mRNAs. These findings indicate a conserved role for J regulating transcription termination and expression of genes within polycistronic gene clusters in trypanosomatids. In contrast to the expectations often attributed to opposing transcription, the essential nature of J in Leishmania spp. is related to its role in gene repression rather than preventing transcriptional interference resulting from read through and dual strand transcription.

Histone H3 Variant Regulates RNA Polymerase II Transcription Termination and Dual Strand Transcription of siRNA Loci in Trypanosoma brucei.

Base J, ß-D-glucosyl-hydroxymethyluracil, is a chromatin modification of thymine in the nuclear DNA of flagellated protozoa of the order Kinetoplastida. In Trypanosoma brucei, J is enriched, along with histone H3 variant (H3.V), at sites involved in RNA Polymerase (RNAP) II termination and telomeric sites involved in regulating variant surface glycoprotein gene (VSG) transcription by RNAP I. Reduction of J in T. brucei indicated a role of J in the regulation of RNAP II termination, where the loss of J at specific sites within polycistronic gene clusters led to read-through transcription and increased expression of downstream genes. We now demonstrate that the loss of H3.V leads to similar defects in RNAP II termination within gene clusters and increased expression of downstream genes. Gene derepression is intensified upon the subsequent loss of J in the H3.V knockout. mRNA-seq indicates gene derepression includes VSG genes within the silent RNAP I transcribed telomeric gene clusters, suggesting an important role for H3.V in telomeric gene repression and antigenic variation. Furthermore, the loss of H3.V at regions of overlapping transcription at the end of convergent gene clusters leads to increased nascent RNA and siRNA production. Our results suggest base J and H3.V can act independently as well as synergistically to regulate transcription termination and expression of coding and non-coding RNAs in T. brucei, depending on chromatin context (and transcribing polymerase). As such these studies provide the first direct evidence for histone H3.V negatively influencing transcription elongation to promote termination.

Fragment ion patchwork quantification for measuring site-specific acetylation degrees.

We introduce fragment ion patchwork quantification as a new mass spectrometry-based approach for the highly accurate quantification of site-specific acetylation degrees. This method combines (13)C1-acetyl derivatization on the protein level, proteolysis by low-specificity proteases and quantification on the fragment ion level. Acetylation degrees are determined from the isotope patterns of acetylated b and y ions. We show that this approach allows to determine site-specific acetylation degrees of all lysine residues for all core histones of Trypanosoma brucei. In addition, we demonstrate how this approach can be used to identify substrate sites of histone acetyltransferases.

Novel insights into RNP granules by employing the trypanosome's microtubule skeleton as a molecular sieve.

RNP granules are ribonucleoprotein assemblies that regulate the post-transcriptional fate of mRNAs in all eukaryotes. Their exact function remains poorly understood, one reason for this is that RNP granule purification has not yet been achieved. We have exploited a unique feature of trypanosomes to prepare a cellular fraction highly enriched in starvation stress granules. First, granules remain trapped within the cage-like, subpellicular microtubule array of the trypanosome cytoskeleton while soluble proteins are washed away. Second, the microtubules are depolymerized and the granules are released. RNA sequencing combined with single molecule mRNA FISH identified the short and highly abundant mRNAs encoding ribosomal mRNAs as being excluded from granules. By mass spectrometry we have identified 463 stress granule candidate proteins. For 17/49 proteins tested by eYFP tagging we have confirmed the localization to granules, including one phosphatase, one methyltransferase and two proteins with a function in trypanosome life-cycle regulation. The novel method presented here enables the unbiased identification of novel RNP granule components, paving the way towards an understanding of RNP granule function.