The tritryps comparative repeatome: insights on repetitive element evolution in trypanosomatid pathogens

The major human pathogens Trypanosoma cruzi, Trypanosoma brucei, and Leishmania major are collectively known as the Tritryps. The initial comparative analysis of their genomes has uncovered that Tritryps share a great number of genes, but repetitive DNA seems to be extremely variable between them. However, the in-depth characterization of repetitive DNA in these pathogens has been in part neglected, mainly due to the well-known technical challenges of studying repetitive sequences from de novo assemblies using short reads. Here, we compared the repetitive DNA repertories between the Tritryps genomes using genome-wide, low-coverage Illumina sequencing coupled to RepeatExplorer analysis. Our work demonstrates that this extensively implemented approach for studying higher eukaryote repeatomes is also useful for protozoan parasites like trypanosomatids, as we recovered previously observed differences in the presence and amount of repetitive DNA families. Additionally, our estimations of repetitive DNA abundance were comparable to those obtained from enhanced-quality assemblies using longer reads. Importantly, our methodology allowed us to describe a previously undescribed transposable element in Leishmania major (TATE element), highlighting its potential to accurately recover distinctive features from poorly characterized repeatomes. Together, our results support the application of this low-cost, low-coverage sequencing approach for the extensive characterization of repetitive DNA evolutionary dynamics in trypanosomatid and other protozoan genomes

Satellitome analysis of Rhodnius prolixus, one of the main Chagas disease vector species

The triatomine Rhodnius prolixus is the main vector of Chagas disease in countries such as Colombia and Venezuela, and the first kissing bug whose genome has been sequenced and assembled. In the repetitive genome fraction (repeatome) of this species, the transposable elements represented 19% of R. prolixus genome, being mostly DNA transposon (Class II elements). However, scarce information has been published regarding another important repeated DNA fraction, the satellite DNA (satDNA), or satellitome. Here, we offer, for the first time, extended data about satellite DNA families in the R. prolixus genome using bioinformatics pipeline based on low-coverage sequencing data. The satellitome of R. prolixus represents 8% of the total genome and it is composed by 39 satDNA families, including four satDNA families that are shared with Triatoma infestans, as well as telomeric (TTAGG)n and (GATA)n repeats, also present in the T. infestans genome. Only three of them exceed 1% of the genome. Chromosomal hybridization with these satDNA probes showed dispersed signals over the euchromatin of all chromosomes, both in autosomes and sex chromosomes. Moreover, clustering analysis revealed that most abundant satDNA families configured several superclusters, indicating that R. prolixus satellitome is complex and that the four most abundant satDNA families are composed by different subfamilies. Additionally, transcription of satDNA families was analyzed in different tissues, showing that 33 out of 39 satDNA families are transcribed in four different patterns of expression across samples

Holocentric chromosome evolution in kissing bugs (Hemiptera: Reduviidae: Triatominae): diversification of repeated sequences

Background: The analysis of the chromosomal and genome evolution in organisms with holocentric chromosomes is restricted by the lack of primary constriction or centromere. An interesting group is the hemipteran subfamily Triatominae, vectors of Chagas disease, which affects around 6 to 7 million people worldwide. This group exhibits extensive variability in the number and chromosomal location of repeated sequences such as heterochromatin and ribosomal genes. This paper tries to reveal the significant differences of the repeated sequences among Triatoma species through the use of genomic DNA probes. Methods: We analysed the chromosomal distribution and evolution of repeated sequences in Triatoma species by genomic in situ hybridization (GISH) using genomic DNA probes from two North American Triatoma species. These genomic probes were hybridized both on their own chromosomes and on other Triatoma species from North and South America, with different amounts and chromosome location of C-heterochromatin. The results were compared with those previously described using South American Triatoma genomic probes. Results: We observed two chromosomal hybridization patterns: (i) very intense hybridization signals concentrated on specific chromosomal regions or particular chromosomes; and (ii) lower intensity hybridization signals dispersed along all chromosomes. Self-GISH on T. rubrofasciata and T. dimidiata chromosomes presented strong hybridization signals on all C-heterochromatin regions. However, when we perform genomic cross-hybridizations, only strong signals are detected on the Y chromosome, leaving the C-heterochromatic autosomal regions unmarked. Conclusions: We confirm that repeated DNA of the Y chromosome is shared among Triatoma species and probably represents an ancestral character of the Triatomini tribe. On the contrary, autosomal heterochromatic regions are constituted by species-specific DNA repeats, most probably satDNA families, suggesting that Triatoma speciation involved the amplification of diverse types of autosomal repeats. Molecular characterization of principal repetitive DNAs seems to be an appropriate approach to infer evolutionary relationships in triatomines

Comparative repeatome analysis on Triatoma infestans Andean and Non-Andean lineages, main vector of Chagas disease

Triatoma infestans is the most important Chagas disease vector in South America. Two main evolutionary lineages, named Andean and non-Andean, have been recognized by geographical distribution, phenetic and genetic characteristics. One of the main differences is the genomic size, varying over 30% in their haploid DNA content. Here we realize a genome wide analysis to compare the repetitive genome fraction (repeatome) between both lineages in order to identify the main repetitive DNA changes occurred during T. infestans differentiation process. RepeatExplorer analysis using Illumina reads showed that both lineages exhibit the same amount of non-repeat sequences, and that satellite DNA is by far the major component of repetitive DNA and the main responsible for the genome size differentiation between both lineages. We characterize 42 satellite DNA families, which are virtually all present in both lineages but with different amount in each lineage. Furthermore, chromosomal location of satellite DNA by fluorescence in situ hybridization showed that genomic variations in T. infestans are mainly due to satellite DNA families located on the heterochromatic regions. The results also show that many satDNA families are located on the euchromatic regions of the chromosomes

Biological attributes of the kissing bug Triatoma rubrofasciata from Vietnam

Background Triatoma rubrofasciata is the only kissing bug species distributed globally. In the Americas, this species transmits the parasite Trypanosoma cruzi, responsible for Chagas disease. The presence of T. rubrofasciata in several Asian countries has greatly increased recently. In Vietnam, it is found in large numbers, closely associated with human environments. Although T. rubrofasciata from Asia is not infected with Tryp. cruzi, it carries other parasites such as Trypanosoma lewisi and Trypanosoma conorhini. Reports of bites by T. rubrofasciata have increased significantly in several places of Vietnam, becoming a public health problem as it produces severe anaphylactic reactions. Methods Specimens of T. rubrofasciata were collected from seven provinces in central Vietnam. We analyzed different biological attributes (life-cycle, starvation resistance, feeding and reproductive capacities) and genetic characteristics (chromosomes and DNA sequences) of T. rubrofasciata from Vietnam and compared them with Brazilian specimens. Natural infection with Tryp. conorhini and Tryp. lewisi were analyzed in a sample of 100 collected insects. Results Species identification of T. rubrofasciata from central Vietnam was corroborated by genetic markers. Cytogenetic analyses showed that T. rubrofasciata from central Vietnam share the same chromosomal characteristics with individuals from Brazil and Hanoi. DNA sequence analyses of a mitochondrial cytochrome b gene fragment showed little variation between Old and New World specimens. Our study sample, compared with Brazilian individuals, showed a higher survival capacity revealed by a higher hatching rate (98% compared with 80.5%), a larger amount of blood taken in single meal and long-term starvation resistance. Furthermore, this species had a high natural rate of infection with Tryp. conorhini (46%) and Tryp. lewisi (27%). Conclusions For T. rubrofasciata of Vietnam, a high rate of fecundity throughout the year, a high capacity for starvation, and its occurrence in synanthropic environments of urban areas with a high availability of food sources are risk factors to be taken into account by vector control campaigns. The several allergic reactions caused by their bites and their high infection with Tryp. lewisi highlight the need to implement specific control programmes for T. rubrofasciata in Vietnam

Expanding an expanded genome: long-read sequencing of Trypanosoma cruzi

Although the genome of Trypanosoma cruzi, the causative agent of Chagas disease, was first made available in 2005, with additional strains reported later, the intrinsic genome complexity of this parasite (the abundance of repetitive sequences and genes organized in tandem) has traditionally hindered high-quality genome assembly and annotation. This also limits diverse types of analyses that require high degrees of precision. Long reads generated by third-generation sequencing technologies are particularly suitable to address the challenges associated with T. cruzi's genome since they permit direct determination of the full sequence of large clusters of repetitive sequences without collapsing them. This, in turn, not only allows accurate estimation of gene copy numbers but also circumvents assembly fragmentation. Here, we present the analysis of the genome sequences of two T. cruzi clones: the hybrid TCC (TcVI) and the non-hybrid Dm28c (TcI), determined by PacBio Single Molecular Real-Time (SMRT) technology. The improved assemblies herein obtained permitted us to accurately estimate gene copy numbers, abundance and distribution of repetitive sequences (including satellites and retroelements). We found that the genome of T. cruzi is composed of a 'core compartment' and a 'disruptive compartment' which exhibit opposite GC content and gene composition. Novel tandem and dispersed repetitive sequences were identified, including some located inside coding sequences. Additionally, homologous chromosomes were separately assembled, allowing us to retrieve haplotypes as separate contigs instead of a unique mosaic sequence. Finally, manual annotation of surface multigene families, mucins and trans-sialidases allows now a better overview of these complex groups of genes

Maxicircle architecture and evolutionary insights into Trypanosoma cruzi complex

We sequenced maxicircles from T. cruzi strains representative of the species evolutionary diversity by using long-read sequencing, which allowed us to uncollapse their repetitive regions, finding that their real lengths range from 35 to 50 kb. T. cruzi maxicircles have a common architecture composed of four regions: coding region (CR), AT-rich region, short (SR) and long repeats (LR). Distribution of genes, both in order and in strand orientation are conserved, being the main differences the presence of deletions affecting genes coding for NADH dehydrogenase subunits, reinforcing biochemical findings that indicate that complex I is not functional in T. cruzi. Moreover, the presence of complete minicircles into maxicircles of some strains lead us to think about the origin of minicircles. Finally, a careful phylogenetic analysis was conducted using coding regions of maxicircles from up to 29 strains, and 1108 single copy nuclear genes from all of the DTUs, clearly establishing that taxonomically T. cruzi is a complex of species composed by group 1 that contains clades A (TcI), B (TcIII) and D (TcIV), and group 2 (1 and 2 do not coincide with groups I and II described decades ago) containing clade C (TcII), being all hybrid strains of the BC type. Three variants of maxicircles exist in T. cruzi: a, b and c, in correspondence with clades A, B, and C from mitochondrial phylogenies. While A and C carry maxicircles a and c respectively, both clades B and D carry b maxicircle variant; hybrid strains also carry the b- variant. We then propose a new nomenclature that is self-descriptive and makes use of both the phylogenetic relationships and the maxicircle variants present in T. cruzi

Multidisciplinary approach detects speciation within the kissing bug Panstrongylus rufotuberculatus populations (Hemiptera, Heteroptera, Reduviidae)

BACKGROUND Panstrongylus rufotuberculatus (Hemiptera-Reduviidae) is a triatomine species with a wide geographic distribution and a broad phenotypic variability. In some countries, this species is found infesting and colonising domiciliary ecotopes representing an epidemiological risk factor as a vector of Trypanosoma cruzi, etiological agent of Chagas disease. In spite of this, little is known about P. rufotuberculatus genetic diversity. METHODS Cytogenetic studies and DNA sequence analyses of one nuclear (ITS-2) and two mitochondrial DNA sequences (cyt b and coI) were carried out in P. rufotuberculatus individuals collected in Bolivia, Colombia, Ecuador and Mexico. Moreover, a geometric morphometrics study was applied to Bolivian, Colombian, Ecuadorian and French Guiana samples. OBJECTIVES To explore the genetic and phenetic diversity of P. rufotuberculatus from different countries, combining chromosomal studies, DNA sequence analyses and geometric morphometric comparisons. FINDINGS We found two chromosomal groups differentiated by the number of X chromosomes and the chromosomal position of the ribosomal DNA clusters. In concordance, two main morphometric profiles were detected, clearly separating the Bolivian sample from the other ones. Phylogenetic DNA analyses showed that both chromosomal groups were closely related to each other and clearly separated from the remaining Panstrongylus species. High nucleotide divergence of cyt b and coI fragments were observed among P. rufotuberculatus samples from Bolivia, Colombia, Ecuador and Mexico (Kimura 2-parameter distances higher than 9%). MAIN CONCLUSIONS Chromosomal and molecular analyses supported that the two chromosomal groups could represent different closely related species. We propose that Bolivian individuals constitute a new Panstrongylus species, being necessary a detailed morphological study for its formal description. The clear morphometric discrimination based on the wing venation pattern suggests such morphological description might be conclusive.MEC-DICYT: II/FVF/2019/054CSIC: No. 16

Evolución cariotípica en insectos vectores de la enfermedad de chagas (Hemiptera, Reduviidae, Triatominae) :análisis del Cluster Ribosomal 45S

Orientadora: Dra. Yanina Panzera Crespo.Co-orientador: Francisco Panzera Arball