Intrinsic DNA curvature in trypanosomes
© The Author(s) 2017
Received: 18 January 2017
Accepted: 1 November 2017
Published: 9 November 2017
Trypanosoma cruzi and Trypanosoma brucei are protozoan parasites causing Chagas disease and African sleeping sickness, displaying unique features of cellular and molecular biology. Remarkably, no canonical signals for RNA polymerase II promoters, which drive protein coding genes transcription, have been identified so far. The secondary structure of DNA has long been recognized as a signal in biological processes and more recently, its involvement in transcription initiation in Leishmania was proposed. In order to study whether this feature is conserved in trypanosomatids, we undertook a genome wide search for intrinsic DNA curvature in T. cruzi and T. brucei.
Using a region integrated intrinsic curvature (RIIC) scoring that we previously developed, a non-random distribution of sequence-dependent curvature was observed. High RIIC scores were found to be significantly correlated with transcription start sites in T. cruzi, which have been mapped in divergent switch regions, whereas in T. brucei, the high RIIC scores correlated with sites that have been involved not only in RNA polymerase II initiation but also in termination. In addition, we observed regions with high RIIC score presenting in-phase tracts of Adenines, in the subtelomeric regions of the T. brucei chromosomes that harbor the variable surface glycoproteins genes.
In both T. cruzi and T. brucei genomes, a link between DNA conformational signals and gene expression was found. High sequence dependent curvature is associated with transcriptional regulation regions. High intrinsic curvature also occurs at the T. brucei chromosome subtelomeric regions where the recombination processes involved in the evasion of the immune host system take place. These findings underscore the relevance of indirect DNA readout in these ancient eukaryotes.
Trypanosoma cruzi and Trypanosoma brucei (family Trypanosomatidae, order Kinetoplastida) are flagellated protozoan parasites that cause Chagas disease and African sleeping sickness, respectively. They infect the poorest rural populations in developing countries in tropical and subtropical regions leading to tens of thousands of human deaths every year .
Although trypanosomes share many characteristics at the molecular and biochemical levels, they exhibit different life cycles. T. cruzi has a complex life cycle alternating between two extracellular forms in the triatomine insect: the epimastigote and the infective metacyclic trypomastigote, and two forms in the mammalian host: the intracellular amastigote and the infective bloodstream trypomastigote. In contrast, T. brucei is exclusively extracellular, alternating between the procyclic form in the tsetse fly and the bloodstream form in the mammalian host . In order to evade the immune system, T. brucei forms a dense coat of variant surface glycoproteins (VSG) that are expressed one at a time from telomeric expression sites and are derived from a repertoire of up to 2000 genes [3, 4]. The high expression of a single VSG gene is carried out by an RNA polymerase I (RNAPI) capable of high transcription initiation rates . This mRNA synthesis property by RNAPI is one of the molecular characteristics that distinguish this parasite from other eukaryotic organisms. Recently, we determined that the eukaryotic conserved intrinsic curvature, which is a characteristic of the RNAPI core promoters, is present not only at the rDNA loci of T. brucei, T. cruzi and Leishmania but also within RNAPI promoters involved in transcribing protein coding genes in T. brucei (VSG and procyclins) .
Trypanosomes have additional unique gene organization and expression features, including the organization of genes in directional clusters, the constitutive transcription of large polycistronic gene units, the amplification of genes in response to environmental stimuli, the trans-splicing of mRNAs, the extensive editing of mitochondrial transcripts [7, 8] and the dependence on post transcriptional regulation mechanisms to coordinate gene expression [2, 9]. Subtelomeric regions are composed of variable repetitive elements and contain genes involved in antigenic variation in T. brucei or genes encoding surface antigens in T. cruzi [10, 11]. The hyper-modified base J (β-d-Glucopyranosyloxymethyluracil), predominantly present in repetitive DNA sequences in telomeres and subtelomeres of trypanosomatids, has been more recently localized in RNA polymerase II transcription initiation and termination sites . Canonical signals for RNA polymerase II promoters have only been described for the genes encoding the spliced leader (SL), a small RNA added by trans-splicing to all the protein coding genes . Nonetheless, transcription starting sites (TSS) [14, 15] and histone variants implicated in the initiation process [16, 17] have been described mainly at the strand switch regions (SSRs) that separate the heads of the polycistronic gene units, named divergent SSRs (DSSRs). On the contrary, the SSRs that separate the tails of the polycistronic gene units, named convergent SSRs (CSSRs), have been shown to preferentially contain sites of transcription termination as well as polymerase III transcribed tRNA genes [15, 18]. Transcription also terminates and initiates at internal TSSs, where transcription of an upstream polycistronic unit terminates and transcription of a downstream polycistronic unit on the same strand initiates [17, 19, 20]. A bias in poly-dinucleotides abundance has also been reported for those regions . A link between transcription and DNA replication has been recently described in T. brucei . In spite of the important insights achieved (reviewed in ), the molecular signals underlying these processes remain mostly under investigation [24, 25].
Intrinsically curved DNA structures are often found around origins of DNA replication, DNA recombination loci and in regions that regulate transcription. In eukaryotes, this feature is common not only to the promoters of genes transcribed by RNA polymerase I, but also to many promoters of RNA polymerase II (reviewed in ). We have previously reported an association of transcription start sites with regions of high regional integrated intrinsic curvature (RIIC) score in Leishmania . Though differences between trypanosomes and Leishmania, due to different base composition content and derived intrinsic curvature (IC) may exist, here we studied whether high RIIC regions could also be associated to a biological phenomenon in T. cruzi and T. brucei. We found that high RIIC scores were indeed associated with regions involved in transcriptional regulation, such as DSSRs in T. cruzi and regions enriched in base J in T. brucei. A concentration of high curvature regions in the T. brucei subtelomeres, overlapping regions of silent VSGs, was also observed, and a canonical signal for DNA bending consisting on A tracts repeated in phase was discovered therein. Our findings suggest a link between DNA conformational signals and gene expression in trypanosomes.
Here we present data demonstrating that regions of high curvature are not randomly distributed along T. brucei and T. cruzi chromosomes. In both cases the high curvature signals seem more widespread than the profile observed for Leish, particularly for T. brucei. It is worth noting that the Leishmania intrinsic curvature distribution, compared to other organisms, ranging from prokaryotes to humans, including trypanosomes, showed fewer regions with high curvature and a higher density of regions with lower curvature. So the T. cruzi and T. brucei profile could be attributed to the presence of more regions of high curvature leading to higher median IC (3.32° and 3.52° per helical turn, respectively) compared to Leishmania (2.63 per helical turn).
Considering the additive contribution to the DNA curvature of each sequential nucleotides in a given region, we compiled the data by using the region integrated intrinsic curvature score, we have previously developed. The enrichment of high RIIC in DSSRs in T. cruzi is suggestive of a putative association of DNA intrinsic curvature with transcription initiation as previously observed in Leishmania. Such curvatures may either facilitate the binding of the RNA polymerase and/or help enable the formation of the open DNA complex. Although an enrichment of histone acetylation in the DSSR has been described in T. cruzi , no studies have been carried out to investigate the presence of TSS markers inside the polycistronic units which prevents a genome wide association study. Considering that high curvature regions are also associated with internal TSSs in Leishmania, one can only speculate that a high percentage of the 172 internal high RIIC scoring regions found in T. cruzi may also serve as internal TSSs. If this consideration holds true, and some of the 172 internal sites are actual TSS regions, then the number of true positives would increase while the number of false positives will decrease resulting in a more significant association (with an overall increase in the Mathews correlation coefficient).
A difference between the pattern observed here in T. brucei and the one previously seen in Leishmania is evident. In the latter, the association of high RIIC with transcription initiation is clear . Besides, for this parasite, the location of base J in the core genome corresponds to transcription termination sites [33–35]. In contrast, in T. brucei, a high of RIIC is seen with base J, which in this organism is considered a marker of transcription boundaries , while the correlation with H4K10ac (marker of transcription initiation) is not significant. Indeed, the coincidence in some regions of high RIIC and H4K10ac enriched regions, might be just the consequence of high association of base J with these regions  (68.8%, 86 out of 125). However, since base J is not present in the procyclic stage of the parasite , these coincidences may be only restricted to the bloodstream stage. Recently, an attenuation role in gene expression at specific sites within the polycistronic gene units has been proposed for base J in T. brucei . It is tempting to speculate that the bulky glucosyl moiety of the base J affects the intrinsic curvature and bendability of regulatory regions involved in facilitating the transcription initiation or processivity. This modification being stage-specific, the location of stage-specific genes downstream to the regions with base J would be expected. Interestingly, in addition to the dispersed internal location of base J, its presence is especially abundant at telomeric and subtelomeric regions in T. brucei  where stage-specific genes are located.
From our analysis, in phase A runs emerge as a common motif in the subtelomeric regions in T. brucei. It has been well established that runs of 4 or more adenines with a 10 bp phasing cause DNA bending [37–40]. The molecular structures of these DNA tracts are unusual and vary depending on the genome context . Multiple roles for the A-tract curvature have been proposed . Among those, A tracts have been implicated in organizing DNA architecture, enhancing the recombination process and assisting chromatin structure. It is also worthwhile noting that curved DNA may vary its shape depending on temperature and physiological changes encompassing the VSG expression silencing in the insect form  and tight control in the bloodstream stage.
Although drugs that bind DNA through the recognition and distortion of DNA structure are known to act as trypanocides, their use has been mostly discarded due to the collateral effects on the host . Nonetheless, the work presented here, draws attention to the potential relevance of refining the design of chemotherapeutic agents focusing on the high curvature occurring at regions important for regulation of gene expression and in the case of T. brucei, for host immune evasion.
Highly curved DNA has been recognized as a signal in transcriptional processes both in prokaryotes and eukaryotes (reviewed in ). Using the regional integrated intrinsic curvature scoring, we have previously shown the association of regions of high intrinsic curvature with those related to the transcription initiation in Leishmania . Here we show that in T. cruzi, the divergent strand switches, which are the only regions up to now known to be associated to transcription initiation, are also distinctively characterized by high RIIC scores, supporting the extrapolation of the same conclusion. Meanwhile in T. brucei, the association of regions of high intrinsic curvature with the locations of base J was found to be significant. Since base J has been involved in transcription initiation and termination, this result may indicate a role for DNA intrinsic curvature in these processes.
These findings point out to the existence of putative conserved process- and species specific-DNA architectural signals in kinetoplastids. Further studies would be necessary to deeply analyze the recognition signal commonalities and differences that may be involved in different steps of the transcriptional mechanisms.
In addition, in T. brucei, the remarkable concentration of regions with high intrinsic curvature at the chromosome subtelomeric regions, where the species-specific genes for the gene family of the highly variable surface glycoproteins are located, may suggest a putative involvement of this structural signal either facilitating the recombinational process or mediating the chromatin silencing and/or granting the vast antigenic variability needed for the efficient evasion of the immune host system that T. brucei has developed.
Globally, the data here presented, while establishing particularities within kinetoplastids, underscore the relevance of indirect DNA readout in these ancient eukaryotes.
The genome data for T. brucei strain TREU 927 and the T. cruzi CL Brener Esmeraldo-like contigs were downloaded from TritrypDB (version 2.1). T. brucei chromosome regions were classified as ‘core’ and ‘subtelomeric’ regions as described in .
Genome-wide intrinsic curvature calculation
The bend.it algorithm , kindly provided by Dr. S. Pongor for local runs, was used to obtain the IC values (°/helical turn) for each base on the individual chromosomes which were binned into 200 Kbase segments. The default window size (31 bp), bendability values from nucleosome binding and DNaseI parameters were used. In-house scripts either in Python or in R programming languages were developed to filter and/or further analyze results. For visualization, wiggle (WIG) files were generated for each chromosome.
Genome-wide regional integrated intrinsic curvature calculation
The RIIC score was calculated as the area under the curvature plot using the Riemann sum. To assess if the T. brucei subtelomeric and VSG array regions have a significantly higher RIIC score than the rest of the genome, their RIIC score was compared to a density function representing the population of RIIC scores for equal-length regions in the genome as in . A region was considered highly curved if the RIIC score was bigger than the 95% confidence interval for the population. The genome wide search for high RIIC regions was performed as in  with minor modifications. Briefly, every chromosome was scanned for 600 bp regions with a RIIC score greater than the 80th percentile value for that chromosome. For the counting of T. cruzi SSRs associated with high RIIC, two criteria were established. Namely, a SSR was not considered if it presented internal sequencing gaps and/or was defined by polycistronic units of less than 6 genes.
To test the association between high RIIC containing and H4K10ac containing regions a contingency matrix was built for each chromosome. The matrix was constructed by classifying genomic regions as containing both signals, only one signal or none and counting each instance. The independence of both variables was tested using the Fisher exact test in R using a p value of 0.01 to define significance. The statistical correlation of these regions was tested using the Mathews correlation coefficient working on the previously generated contingency matrix.
For the analysis of motif-based sequence associated to the regions of high intrinsic curvature, the MEME suite was used. The search was performed on a fragment of T. brucei chromosome 9 spanning bases from 1 to 319,439 (the VSG array region). A total of 103 sequences surrounding 30 bp each high IC peak were selected for motif discovery. The randomly shuffled sequences were also submitted to the MEME suite as control.
PS participated in the design of the study, carried out the data analyses, discussed the results and drafted the manuscript. NME participated in the design of the study and discussed the results. BG conceived and participated in the design of the study, discussed the results and drafted the manuscript. All authors read and approved the final manuscript.
We acknowledge Dr. R. Sabatini and Dr. N. Siegel for providing T. brucei Chip-seq data for H4K10ac and base J.
The authors declare that they have no competing interests.
Availability of data and materials
All data are available upon request.
Consent for publication
Ethics approval and consent to participate
This work was supported by Fondo Clemente Estable (Agencia Nacional de Investigación e Innovación); Comisión Sectorial de Investigación Científica (Universidad de la República) and Programa de Desarrollo de Ciencias Básicas. PS received a Ph.D. fellowship from Agencia Nacional de Investigación e Innovación in Uruguay.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Field V, Gautret P, Schlagenhauf P, Burchard GD, Caumes E, Jensenius M, Castelli F, Gkrania-Klotsas E, Weld L, Lopez-Velez R, et al. Travel and migration associated infectious diseases morbidity in Europe, 2008. BMC Infect Dis. 2010;10:330.View ArticlePubMedPubMed CentralGoogle Scholar
- Kramer S. Developmental regulation of gene expression in the absence of transcriptional control: the case of kinetoplastids. Mol Biochem Parasitol. 2012;181:61–72.View ArticlePubMedGoogle Scholar
- Boothroyd JC. Antigenic variation in African trypanosomes. Annu Rev Microbiol. 1985;39:475–502.View ArticlePubMedGoogle Scholar
- Berriman M, Hall N, Sheader K, Bringaud F, Tiwari B, Isobe T, Bowman S, Corton C, Clark L, Cross GA, et al. The architecture of variant surface glycoprotein gene expression sites in Trypanosoma brucei. Mol Biochem Parasitol. 2002;122:131–40.View ArticlePubMedGoogle Scholar
- Rudenko G, Lee MG, Van der Ploeg LH. The PARP and VSG genes of Trypanosoma brucei do not resemble RNA polymerase II transcription units in sensitivity to Sarkosyl in nuclear run-on assays. Nucleic Acids Res. 1992;20:303–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Smircich P, Duhagon MA, Garat B. Conserved curvature of RNA polymerase I core promoter beyond rRNA genes: the case of the Tritryps. Genomics Proteomics Bioinformatics. 2015;13:355–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Campbell DA, Thomas S, Sturm NR. Transcription in kinetoplastid protozoa: why be normal? Microbes Infect. 2003;5:1231–40.View ArticlePubMedGoogle Scholar
- Stuart KD, Schnaufer A, Ernst NL, Panigrahi AK. Complex management: RNA editing in trypanosomes. Trends Biochem Sci. 2005;30:97–105.View ArticlePubMedGoogle Scholar
- Martinez-Calvillo S, Vizuet-de-Rueda JC, Florencio-Martinez LE, Manning-Cela RG, Figueroa-Angulo EE. Gene expression in trypanosomatid parasites. J Biomed Biotechnol. 2010;2010:525241.View ArticlePubMedPubMed CentralGoogle Scholar
- Dreesen O, Cross GA. Consequences of telomere shortening at an active VSG expression site in telomerase-deficient Trypanosoma brucei. Eukaryot Cell. 2006;5:2114–9.View ArticlePubMedPubMed CentralGoogle Scholar
- El-Sayed NM, Myler PJ, Bartholomeu DC, Nilsson D, Aggarwal G, Tran AN, Ghedin E, Worthey EA, Delcher AL, Blandin G, et al. The genome sequence of Trypanosoma cruzi, etiologic agent of Chagas disease. Science. 2005;309:409–15.View ArticlePubMedGoogle Scholar
- Cliffe LJ, Siegel TN, Marshall M, Cross GA, Sabatini R. Two thymidine hydroxylases differentially regulate the formation of glucosylated DNA at regions flanking polymerase II polycistronic transcription units throughout the genome of Trypanosoma brucei. Nucleic Acids Res. 2010;38:3923–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Gilinger G, Bellofatto V. Trypanosome spliced leader RNA genes contain the first identified RNA polymerase II gene promoter in these organisms. Nucleic Acids Res. 2001;29:1556–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Martinez-Calvillo S, Yan S, Nguyen D, Fox M, Stuart K, Myler PJ. Transcription of Leishmania major Friedlin chromosome 1 initiates in both directions within a single region. Mol Cell. 2003;11:1291–9.View ArticlePubMedGoogle Scholar
- Martinez-Calvillo S, Nguyen D, Stuart K, Myler PJ. Transcription initiation and termination on Leishmania major chromosome 3. Eukaryot Cell. 2004;3:506–17.View ArticlePubMedPubMed CentralGoogle Scholar
- Respuela P, Ferella M, Rada-Iglesias A, Aslund L. Histone acetylation and methylation at sites initiating divergent polycistronic transcription in Trypanosoma cruzi. J Biol Chem. 2008;283:15884–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Siegel TN, Hekstra DR, Kemp LE, Figueiredo LM, Lowell JE, Fenyo D, Wang X, Dewell S, Cross GA. Four histone variants mark the boundaries of polycistronic transcription units in Trypanosoma brucei. Genes Dev. 2009;23:1063–76.View ArticlePubMedPubMed CentralGoogle Scholar
- Padilla-Mejia NE, Florencio-Martinez LE, Figueroa-Angulo EE, Manning-Cela RG, Hernandez-Rivas R, Myler PJ, Martinez-Calvillo S. Gene organization and sequence analyses of transfer RNA genes in Trypanosomatid parasites. BMC Genomics. 2009;10:232.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomas S, Green A, Sturm NR, Campbell DA, Myler PJ. Histone acetylations mark origins of polycistronic transcription in Leishmania major. BMC Genomics. 2009;10:152.View ArticlePubMedPubMed CentralGoogle Scholar
- Kolev NG, Franklin JB, Carmi S, Shi H, Michaeli S, Tschudi C. The transcriptome of the human pathogen Trypanosoma brucei at single-nucleotide resolution. PLoS Pathog. 2010;6:e1001090.View ArticlePubMedPubMed CentralGoogle Scholar
- Duhagon MA, Smircich P, Forteza D, Naya H, Williams N, Garat B. Comparative genomic analysis of dinucleotide repeats in Tritryps. Gene. 2011;487:29–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Tiengwe C, Marcello L, Farr H, Dickens N, Kelly S, Swiderski M, Vaughan D, Gull K, Barry JD, Bell SD, McCulloch R. Genome-wide analysis reveals extensive functional interaction between DNA replication initiation and transcription in the genome of Trypanosoma brucei. Cell Rep. 2012;2:185–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Choi J, El-Sayed NM. Functional genomics of trypanosomatids. Parasite Immunol. 2012;34:72–9.View ArticlePubMedGoogle Scholar
- Mao Y, Najafabadi HS, Salavati R. Genome-wide computational identification of functional RNA elements in Trypanosoma brucei. BMC Genomics. 2009;10:355.View ArticlePubMedPubMed CentralGoogle Scholar
- Kelly S, Wickstead B, Maini PK, Gull K. Ab initio identification of novel regulatory elements in the genome of Trypanosoma brucei by Bayesian inference on sequence segmentation. PLoS ONE. 2011;6:e25666.View ArticlePubMedPubMed CentralGoogle Scholar
- Ohyama T. DNA conformation and transcription. Georgetown: Landes Bioscience: Springer Science Business Media; c2005.Google Scholar
- Smircich P, Forteza D, El-Sayed NM, Garat B. Genomic analysis of sequence-dependent DNA curvature in Leishmania. PLoS ONE. 2013;8:e63068.View ArticlePubMedPubMed CentralGoogle Scholar
- Weatherly DB, Boehlke C, Tarleton RL. Chromosome level assembly of the hybrid Trypanosoma cruzi genome. BMC Genomics. 2009;10:255.View ArticlePubMedPubMed CentralGoogle Scholar
- Vlahovicek K, Kajan L, Pongor S. DNA analysis servers: plot.it, bend.it, model.it and IS. Nucleic Acids Res. 2003;31:3686–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Tosato V, Ciarloni L, Ivens AC, Rajandream MA, Barrell BG, Bruschi CV. Secondary DNA structure analysis of the coding strand switch regions of five Leishmania major Friedlin chromosomes. Curr Genet. 2001;40:186–94.View ArticlePubMedGoogle Scholar
- Borst P, Sabatini R. Base J: discovery, biosynthesis, and possible functions. Annu Rev Microbiol. 2008;62:235–51.View ArticlePubMedGoogle Scholar
- Hertz-Fowler C, Figueiredo LM, Quail MA, Becker M, Jackson A, Bason N, Brooks K, Churcher C, Fahkro S, Goodhead I, et al. Telomeric expression sites are highly conserved in Trypanosoma brucei. PLoS ONE. 2008;3:e3527.View ArticlePubMedPubMed CentralGoogle Scholar
- van Luenen HG, Farris C, Jan S, Genest PA, Tripathi P, Velds A, Kerkhoven RM, Nieuwland M, Haydock A, Ramasamy G, et al. Glucosylated hydroxymethyluracil, DNA base J, prevents transcriptional readthrough in Leishmania. Cell. 2012;150:909–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Reynolds D, Cliffe L, Forstner KU, Hon CC, Siegel TN, Sabatini R. Regulation of transcription termination by glucosylated hydroxymethyluracil, base J, in Leishmania major and Trypanosoma brucei. Nucleic Acids Res. 2014;42:9717–29.View ArticlePubMedPubMed CentralGoogle Scholar
- Genest PA, Baugh L, Taipale A, Zhao W, Jan S, van Luenen HG, Korlach J, Clark T, Luong K, Boitano M, et al. Defining the sequence requirements for the positioning of base J in DNA using SMRT sequencing. Nucleic Acids Res. 2015;43:2102–15.View ArticlePubMedPubMed CentralGoogle Scholar
- van Leeuwen F, Wijsman ER, Kieft R, van der Marel GA, van Boom JH, Borst P. Localization of the modified base J in telomeric VSG gene expression sites of Trypanosoma brucei. Genes Dev. 1997;11:3232–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Trifonov EN. Sequence-dependent deformational anisotropy of chromatin DNA. Nucleic Acids Res. 1980;8:4041–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Marini JC, Levene SD, Crothers DM, Englund PT. Bent helical structure in kinetoplast DNA. Proc Natl Acad Sci USA. 1982;79:7664–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Nadeau JG, Crothers DM. Structural basis for DNA bending. Proc Natl Acad Sci USA. 1989;86:2622–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Haran TE, Mohanty U. The unique structure of A-tracts and intrinsic DNA bending. Q Rev Biophys. 2009;42:41–81.View ArticlePubMedGoogle Scholar
- De Santis P, Scipioni A. Sequence-dependent collective properties of DNAs and their role in biological systems. Phys Life Rev. 2013;10:41–67.View ArticlePubMedGoogle Scholar
- Pandya UM, Sandhu R, Li B. Silencing subtelomeric VSGs by Trypanosoma brucei RAP1 at the insect stage involves chromatin structure changes. Nucleic Acids Res. 2013;41:7673–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Wainwright M. Dyes, trypanosomiasis and DNA: a historical and critical review. Biotech Histochem. 2010;85:341–54.View ArticlePubMedGoogle Scholar