Experimental and Bioinformatic upgrade for genome-wide mapping of nucleosomes in Trypanosoma cruzi

1In Trypanosoma cruzi, as in every eukaryotic cell, DNA is packaged into chromatin by octamers of histone proteins that constitute nucleosomes. Besides compacting DNA, nucleosomes control DNA dependent processes by modulating the access of DNA binding proteins to regulatory elements on the DNA; or by providing the platform for additional layers of regulation given by histone variants and histone post-translational modifications. In trypanosomes, protein coding genes are constitutively transcribed as polycistronic units. Therefore, gene expression is controlled mainly post transcriptionally. However, chromatin organization and the histone code influence transcription, cell cycle progression, replication and DNA repair. Hence, determining nucleosome position is of uppermost importance to understand the peculiarities of these processes in trypanosomes. Digestion of chromatin with micrococcal nuclease followed by deep sequencing has been widely applied for genome-wide mapping of nucleosomes in several organisms. Nonetheless, this parasite presents numerous singularities. On one hand, special growth conditions and cell manipulation are required. On the other hand, chromatin organization shows some uniqueness that demands a specially designed analytical approach. An additional entanglement is given by the nature of its genome harboring a large content of repetitive sequences and the poor quality of the genome assembly and annotation of many strains. Here, we adapted this broadly used method to the hybrid reference strain, CL Brener. Particularly, we developed an exhaustive and thorough computational workflow for data analysis, highlighting the relevance of using its whole genome as a reference instead of the commonly used Esmeraldo-like haplotype. Moreover, the performance of two aligners, Bowtie2 and HISAT2 was tested to find the most appropriate tool to map any genomic read to reference genomes bearing this complexity.