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.

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.

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.

ddPCRclust: an R package and Shiny app for automated analysis of multiplexed ddPCR data.

Motivation: Droplet digital PCR (ddPCR) is an emerging technology for quantifying DNA. By partitioning the target DNA into ∼20 000 droplets, each serving as its own PCR reaction compartment, a very high sensitivity of DNA quantification can be achieved. However, manual analysis of the data is time consuming and algorithms for automated analysis of non-orthogonal, multiplexed ddPCR data are unavailable, presenting a major bottleneck for the advancement of ddPCR transitioning from low-throughput to high-throughput. Results: ddPCRclust is an R package for automated analysis of data from Bio-Rad's droplet digital PCR systems (QX100 and QX200). It can automatically analyze and visualize multiplexed ddPCR experiments with up to four targets per reaction. Results are on par with manual analysis, but only take minutes to compute instead of hours. The accompanying Shiny app ddPCRvis provides easy access to the functionalities of ddPCRclust through a web-browser based GUI. Availability and implementation: R package: https://github.com/bgbrink/ddPCRclust; Interface: https://github.com/bgbrink/ddPCRvis/; Web: https://bibiserv.cebitec.uni-bielefeld.de/ddPCRvis/. Supplementary information: Supplementary data are available at Bioinformatics online.

Omics Fusion - A Platform for Integrative Analysis of Omics Data.

We present Omics Fusion, a new web-based platform for integrative analysis of omics data. Omics Fusion provides a collection of new and established tools and visualization methods to support researchers in exploring omics data, validating results or understanding how to adjust experiments in order to make new discoveries. It is easily extendible and new visualization methods are added continuously. It is available for free under: https://fusion.cebitec.uni-bielefeld.de/.

Dynamics of protein and polar lipid recruitment during lipid droplet assembly in Chlamydomonas reinhardtii.

In plants, neutral lipids are frequently synthesized and stored in seed tissues, where the assembly of lipid droplets (LDs) coincides with the accumulation of triacylglycerols (TAGs). In addition, photosynthetic, vegetative cells can form cytosolic LDs and much less information is known about the makeup and biogenesis of these LDs. Here we focus on Chlamydomonas reinhardtii as a reference model for LDs in a photosynthetic cell, because in this unicellular green alga LD dynamics can be readily manipulated by nitrogen availability. Nitrogen deprivation leads to cellular quiescence during which cell divisions cease and TAGs accumulate. The major lipid droplet protein (MLDP) forms a proteinaceous coat surrounding mature LDs. Reducing the amount of MLDP affects LD size and number, TAG breakdown and timely progression out of cellular quiescence following nitrogen resupply. Depending on nitrogen availability, MLDP recruits different proteins to LDs, tubulins in particular. Conversely, depolymerization of microtubules drastically alters the association of MLDP with LDs. LDs also contain select chloroplast envelope membrane proteins hinting at an origin of LDs, at least in part, from chloroplast membranes. Moreover, LD surface lipids are rich in de novo synthesized fatty acids, and are mainly composed of galactolipids which are typical components of chloroplast membranes. The composition of the LD membrane is altered in the absence of MLDP. Collectively, our results suggest a mechanism for LD formation in C. reinhardtii involving chloroplast envelope membranes by which specific proteins are recruited to LDs and a specialized polar lipid monolayer surrounding the LD is formed.