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

Mono-allelic epigenetic regulation of bi-directional polycistronic transcription initiation by RNA Polymerase II in Trypanosoma brucei

Unique for a eukaryote, protein-coding genes in trypanosomes are arranged in polycistronic units (PTUs). This genome arrangement has led to a model where Pol II transcription of PTUs is unregulated. The initial step in trypanosome lytic factor (TLF) mediated lysis of Trypanosoma brucei requires high affinity haptoglobin/hemoglobin receptor (HpHbR) binding. Here we demonstrate that by in vitro selection with TLF, resistance is obtained in a stepwise process correlating with loss of HpHbR expression at an allelic level. RNA-seq, Pol II ChIP and run-on analysis indicate HpHbR silencing is at the transcriptional level, where loss of Pol II binding at the promoter region specifically shuts down transcription of the HpHbR containing gene cluster and the adjacent opposing gene cluster. Reversible transcriptional silencing of the divergent PTUs correlates with DNA base J modification of the shared promoter region. Therefore, epigenetic mechanisms exist to regulate gene expression via Pol II transcription initiation of gene clusters in a mono-allelic fashion. These findings suggest epigenetic chromatin-based regulation of gene expression is deeply conserved among eukaryotes, including primitive eukaryotes that rely on polycistronic transcription.

Milletdb: a multi-omics database to accelerate the research of functional genomics and molecular breeding of millets

Millets are a class of nutrient-rich coarse cereals with high resistance to abiotic stress; thus, they guarantee food security for people living in areas with extreme climatic conditions and provide stress-related genetic resources for other crops. However, no platform is available to provide a comprehensive and systematic multi-omics analysis for millets, which seriously hinders the mining of stress-related genes and the molecular breeding of millets. Here, a free, web-accessible, user-friendly millets multi-omics database platform (Milletdb, http://milletdb.novogene.com) has been developed. The Milletdb contains six millets and their one related species genomes, graph-based pan-genomics of pearl millet, and stress-related multi-omics data, which enable Milletdb to be the most complete millets multi-omics database available. We stored GWAS (genome-wide association study) results of 20 yield-related trait data obtained under three environmental conditions [field (no stress), early drought and late drought] for 2 years in the database, allowing users to identify stress-related genes that support yield improvement. Milletdb can simplify the functional genomics analysis of millets by providing users with 20 different tools (e.g., 'Gene mapping', 'Co-expression', 'KEGG/GO Enrichment' analysis, etc.). On the Milletdb platform, a gene PMA1G03779.1 was identified through 'GWAS', which has the potential to modulate yield and respond to different environmental stresses. Using the tools provided by Milletdb, we found that the stress-related PLATZs TFs (transcription factors) family expands in 87.5% of millet accessions and contributes to vegetative growth and abiotic stress responses. Milletdb can effectively serve researchers in the mining of key genes, genome editing and molecular breeding of millets.

Identification of sequence-specific promoters driving polycistronic transcription initiation by RNA polymerase II in trypanosomes

Protein-coding genes in trypanosomes occur in polycistronic transcription units (PTUs). How RNA polymerase II (Pol II) initiates transcription of PTUs has not been resolved; the current model favors chromatin modifications inducing transcription rather than sequence-specific promoters. Here, we uncover core promoters by functional characterization of Pol II peaks identified by chromatin immunoprecipitation sequencing (ChIP-seq). Two distinct promoters are located between divergent PTUs, each driving unidirectional transcription. Detailed analysis identifies a 75-bp promoter that is necessary and sufficient to drive full reporter expression and contains functional motifs. Analysis of further promoters suggests transcription initiation is regulated and promoters are either focused or dispersed. In contrast to the previous model of unregulated and promoter-independent transcription initiation, we find that sequence-specific promoters determine the initiation of Pol II transcription of protein-coding genes PTUs. These findings in Trypanosoma brucei suggest that in addition of chromatin modifications, promoter motifs-based regulation of gene expression is deeply conserved among eukaryotes.

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.

Histone modification and deposition are key regulators of transcription. Here, Kraus et al. provide a quantitative histone acetylome for Trypanosoma brucei, identify histone modifications enriched at transcription start sites, and show how H4 and H2A.Z acetylation affect histone deposition and transcription in trypanosomes.

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.

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.

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.