Acessibilidade / Reportar erro

A nuclear ribosomal DNA pseudogene in triatomines opens a new research field of fundamental and applied implications in Chagas disease

Abstract

A pseudogene, designated as "ps(5.8S+ITS-2)", paralogous to the 5.8S gene and internal transcribed spacer (ITS)-2 of the nuclear ribosomal DNA (rDNA), has been recently found in many triatomine species distributed throughout North America, Central America and northern South America. Among characteristics used as criteria for pseudogene verification, secondary structures and free energy are highlighted, showing a lower fit between minimum free energy, partition function and centroid structures, although in given cases the fit only appeared to be slightly lower. The unique characteristics of "ps(5.8S+ITS-2)" as a processed or retrotransposed pseudogenic unit of the ghost type are reviewed, with emphasis on its potential functionality compared to the functionality of genes and spacers of the normal rDNA operon. Besides the technical problem of the risk for erroneous sequence results, the usefulness of "ps(5.8S+ITS-2)" for specimen classification, phylogenetic analyses and systematic/taxonomic studies should be highlighted, based on consistence and retention index values, which in pseudogenic sequence trees were higher than in functional sequence trees. Additionally, intraindividual, interpopulational and interspecific differences in pseudogene amount and the fact that it is a pseudogene in the nuclear rDNA suggests a potential relationships with fitness, behaviour and adaptability of triatomine vectors and consequently its potential utility in Chagas disease epidemiology and control.

triatomines; Chagas disease; rDNA pseudogene; secondary structures; free energy; functionality


The genes, basic units for transferring hereditary information, provide structure, function and regulation to a biological system. A gene must go through several steps from a genetic DNA sequence to a fully-functional protein. These steps include transcription, pre-mRNA processing, translation and protein folding. The sequence may be considered nonfunctional if any of the aforementioned steps fails. The most commonly identified disablements are stop codons and frame shifts, which almost universally stop the translation of a functional protein product. Pseudogenes are homologous sequences arising from currently or evolutionarily active genes that have lost their ability to function as a result of disrupted transcription or translation. They may contain stop codons, repetitive elements, have frame shifts and/or lack of transcription. However, they might retain gene-like features. Pseudogenes are of particular interest to biologists since they can interfere with gene centric studies [such as de novo gene prediction and polymerase chain reaction (PCR) amplification] and also to evolutionary biologists because of the possibility to study their age and mutational rates and tendencies (Rouchka & Cha 2009Rouchka EC, Cha IE 2009. Current trends in pseudogene detection and characterization. Curr Bioinform 4: 112-119.).

Although traditionally noted to be nonfunctional, the view of pseudogenes as genetic elements similar to functional genes but without functional properties has recently been put in question (Balakirev & Ayala 2003Balakirev ES, Ayala FJ 2003. Pseudogenes: are they "junk" or functional DNA? Annu Rev Genet 237: 123-151., Sasidharan & Gerstein 2008Sasidharan R, Gerstein M 2008. Genomics: protein fossils live on as RNA. Nature 453: 729-731.). It was traditionally noted that since pseudogenes do not produce protein products, they are typically not under selective evolutionary pressure and thus evolve at rates consistent with neutral drift (Friedberg & Rhoads 2000Friedberg F, Rhoads AR 2000. Calculation and verification of the ages of retroprocessed pseudogenes. Mol Phylogenet Evol 16: 127-130.). Nevertheless, there is recent evidence indicating that a number of pseudogenes may be actively transcribed, as seen in the human and mouse genomes (Zhang et al. 2004Zhang Z, Carriero N, Gerstein M 2004. Comparative analysis of processed pseudogenes in the mouse and human genomes. Trends Genet 20: 62-67. , Harrison et al. 2005Harrison PM, Zheng D, Zhang Z, Carriero N, Gerstein M 2005. Transcribed processed pseudogenes in the human genome: an intermediate form of expressed retrosequence lacking protein-coding ability. Nucleic Acids Res 33: 2374-2383., Zheng et al. 2005Zheng D, Gerstein MB 2007. The ambiguous boundary between genes and pseudogenes: the dead rise up or do they? Trends Genet 23: 219-224. ). Although none of these pseudogenes are translated into proteins due to sequence disablements, it has been hypothesised that they may have roles in gene regulation by using their sequences complementary to the homologous functional gene. Transcribed processed pseudogenes have an additional effect that they themselves can become duplicated, resulting in "duplicated-processed" pseudogenes (Zhang & Gerstein 2003, Zhang et al. 2008).

Many characterised pseudogenes have been implicated in regulation of gene expression, gene regulation, and provide a potential source of genetic diversity through recombination with functional genes or exon shuffling (Pavlicek et al. 2006Pavlicek A, Gentles AJ, Paces J, Paces V, Jurka J 2006. Retroposition of processed pseudogenes: the impact of RNA stability and translational control. Trends Genet 22: 69-73.). Interactions have been described between transporter genes and pseudogenes that suggest that expression of the gene is regulated, in part, by transcription of the pseudogene (Piehler et al. 2008Piehler AP, Hellum M, Wenzel JJ, Kaminsk E, Haug KBF, Kierulf P, Kaminski WE 2008. The human ABC transporter pseudogene family: evidence for transcription and genepseudogene interference. BMC Genomics 9: 165.). A subset of mammalian pseudogenes has been demonstrated to be responsible for generating small interfering RNAs (siRNAs) by forming double-stranded RNA sequences with the corresponding protein-coding messenger RNA (mRNA) (Tam et al. 2008Tam OH, Aravin AA, Stein P, Girard A, Murchison EP, Cheloufi S, Hodges E, Anger M, Sachidanandam R, Schultz RM, Hannon GJ 2008. Pseudogene-derived small interfering RNAs regulate gene expression in mouse oocytes. Nature 453: 534-538.). These pseudogene-derived siRNAs are then in turn responsible for regulating the functional gene from which the pseudogene originates by acting to degrade the functional gene's mRNA transcripts (Gray et al. 2006Gray TA, Wilson A, Fortin PJ, Nicholls RD 2006. The putatively functional Mkrn1-p1 pseudogene is neither expressed nor imprinted, nor does it regulate its source gene in trans. Proc Natl Acad Sci USA 103: 12039-12044., Sasidharan & Gerstein 2008). A classification has even been proposed to differentiate between ghost pseudogenes that have some intermediate functionality (such as a regulatory function or transcriptional activity) and dead pseudogenes that do not have any indication of functionality and therefore are subject to neutral drift (Zheng & Gerstein 2007).

In triatomines vectors of Chagas disease, all pseudogenes so far described have been found in the mitochondrial DNA (mtDNA) genome and proved to be nonfunctional due to the stop codons they present (Mas-Coma & Bargues 2009Mas-Coma S, Bargues MD 2009. Populations, hybrids and the systematic concepts of species and subspecies in Chagas disease triatomine vectors inferred from nuclear ribosomal and mitochondrial DNA. Acta Trop 110: 112-136. ). Recently, however, a new paralogous sequence of the 5.8S gene and the internal transcribed spacer (ITS)-2 of the nuclear ribosomal DNA (rDNA) have been detected to be present in many triatomine species distributed throughout North, Central and northern South America and which has revealed to be a pseudogene with surprising characteristics suggesting potential functionality (Bargues et al. 2014). This appears to be the fourth report of a pseudogene in the nuclear ribosomal operon of an animal, after two in insects such as the grasshopper Podisma pedestris (Keller et al. 2006Keller I, Chintauan-Marquier IC, Veltsos, P, Nichols RA 2006. Ribosomal DNA in the grasshopper Podisma pedestris: escape from concerted evolution. Genetics 174: 863-874.) and the pseudogenic ITS-2 sequences in the malaria vector Anopheles albitarsis (Li & Wilkerson 2007Li C, Wilkerson RC 2007. Intragenomic rDNA ITS2 variation in the Neotropical Anopheles (Nyssorhynchus) albitarsis complex (Diptera: Culicidae). J Hered 98: 51-59. ) and in a vertebrate such as the stone flounder fish Kareius bicoloratus(Xu et al. 2009Xu J, Zhang Q, Xu X, Wanga Z, Qi J 2009. Intragenomic variability and pseudogenes of ribosomal DNA in stone flounder Kareius bicoloratus. Mol Phylogenet Evol 52: 157-166. ).

Pseudogene verification

The coexistence of two different 5.8S+ITS-2 sequences within the same triatomine specimen, one corresponding to the normal functional rDNA operon and another paralogous one corresponding to the pseudogenic, was confirmed by (i) double signal in the sequencing chromatograms (Fig. 1), (ii) cloning and subsequent clone sequencing, (iii) different specific primer sequencing for each functional and paralogous sequences and (iv) relative quantification by real-time PCR (Bargues et al. 2014).

Fig. 1:
detection of double signal in the sequencing chromatogram fragment obtained when working with a specimen of Triatoma phyllosoma longipennis showing superposition of internal transcribed spacer-2 [functional (H1) and pseudogenic (psH1)]. See corresponding ClustalW2 alignment of the psH1 and H1 paralogous sequences of the haplotype 1 of this species at the bottom.

A combination of different criteria was used to identify the pseudogene. In the 5.8S gene, criteria for pseudogene identification included length variability, lower guanine-cytosine content, mutations regarding the functional uniform sequence and relatively high base substitutions in evolutionary conserved sites (Table I). At ITS-2 level, criteria were the shorter sequence and large proportion of insertions and deletions (indels) (Table I) (Bargues et al. 2014).

TABLE I.
Comparative data (extreme values and means in parentheses) furnishing criteria for the identification of the "ps(5.8S+ITS-2)" pseudogene in triatomine vectors

DNA secondary structures - The tendency of complementary strands of DNA to form double helices is well known since long time ago. Single stranded nucleic acid sequences contain many complementary regions that have the potential to form double helices when the molecule folds back on itself. The resulting pattern of double helical stretches interspersed with loops is what is called the secondary structure of an RNA or DNA. Secondary structure elements may in turn be arranged in space to form three-dimensional tertiary structure, leading to additional noncovalent interactions. Tertiary interactions are weaker than secondary structure and, consequently, RNA folding is regarded as a hierarchical process in which secondary structure forms before tertiary structure (Thirumalai et al. 2001Thirumalai D, Lee N, Woodson SA, Klimov D 2001. Early events in RNA folding. Annu Rev Phys Chem 52: 751-762.). Since formation of tertiary structure usually does not induce changes in secondary structure, the two processes can be described independently.

Functional RNA molecules [transfer RNA, ribosomal RNAs (rRNAs) etc., as opposed to pure coding sequences], usually have characteristic spatial structures and therefore also characteristic secondary structures, that are prerequisites for their function. Their secondary structures are therefore highly conserved in evolution (Hofacker & Stadler 2008Hofacker IL, Fontana W, Stadler PF, Bonhoeffer SL, Tacker M, Schuster P 1994. Fast folding and comparison of RNA secondary structures. Monatsh Chem 125: 167-188. ). Kinetics of RNA secondary structure formation plays an important role in many biological functions. An RNA secondary structure is considered to be locally optimal if there is no lower energy structure that can be obtained by the addition or removal of a single base pair (bp), where energy is defined according to the widely accepted Turner nearest neighbour model (Lorenz & Clote 2011Lorenz WA, Clote P 2011. Computing the partition function for kinetically trapped RNA secondary structures. PLoS ONE 6: e16178.).

Pseudogenic 5.8S and ITS-2 secondary structuresand free energy - They were different from the functional foldings, different one another, showing less negative values for minimum free energy (mfe) and centroid predictions (Table II). A prediction of the secondary structures was made for both functional and paralogous sequences independently for 5.8S and ITS-2. Although the complete sequence was obtained for the functional 5.8S gene, a long fragment (only lacking the 21 bp at the 5' end) was used for secondary structure prediction and free energy evaluation, to be strict in comparisons between functional and pseudogenic sequences. In the case of ITS-2, sequences used for these comparisons were the complete spacer sequences for both functional and pseudogenic paralogues. Associated free energy values were evaluated using the mfe algorithm (Zuker & Stiegler 1981Zuker M 1989. On finding all suboptimal foldings of an RNA molecule. Science 244: 48-52. , Zuker 1989) and the latest free energy rules (Mathews et al. 2004Mathews DH, Disney MD, Childs JL, Schroeder S, Zuker M, Turner DH 2004. Incorporating chemical modification constraints into a dynamic programming algorithm for prediction of RNA secondary structure. Proc Natl Acad Sci USA 101: 7287-7292.). Fold predictions and bp probabilities were made using the RNAfold program (Hofacker et al. 1994). For comparative purposes, the stability of the secondary structure of the putative pseudogenic sequences vs. that of functional sequences was assessed, for each single sequence, by (i) calculating their mfe structure, (ii) computing their partition function (pf) structure and (iii) analysing the centroid structure, in mountain plots of their positional height values.

TABLE II.
Minimum free energy (mfe) prediction for the optimal secondary structure and for the thermodynamic ensemble of the functional (FuH) and pseudogenic (PsH) paralogous sequences of the 5.8S ribosomal RNA gene and the ribosomal DNA internal transcribed spacer (ITS)-2 of selected haplotypes representing 13 Triatominae taxa as comparison examples

The mfe structure of an RNA sequence is the secondary structure that contributes a minimum of free energy. The pf sums all Boltzmann weighted free energies of each secondary structure that is possible given an RNA sequence, thus providing the possibility to calculate base pairing probabilities for each possible pair of bases and the obtaining of an ensemble structure depicting the bp probabilities when pf folding is selected. The centroid structure of an RNA sequence is the secondary structure with minimal bp distance to all other secondary structures in the Boltzmann ensemble [with regard to Boltzmann thermodynamic entropy concepts, see recent reappraisal by Kalinin and Kononogov (2005)Kalinin MI, Kononogov SA 2005. Boltzmann's constant, the energy meaning of temperature and thermodynamic irreversibility. Meas Tech 48: 632-636.]. A mountain plot represents a secondary structure in a plot of height vs. position, where the height mk is given by the number of bp enclosing the base at position k, i.e. loops correspond to plateaus (hairpin loops are peaks) and helices to slopes.

In the pseudogenic 5.8S gene, the corresponding mfe (delta G) values were less negative, ranging from -32.90 to -38.20 kcal/mol and from -18.91 to -26.91 kcal/mol for optimal and centroid secondary structure prediction, respectively (Table II). The lower fit between the three secondary structure prediction approaches assayed (mfe, pf and centroid) became evident (Fig. 2B-F) when compared with the functional uniform 5.8S secondary structure (Fig. 2A).

Fig. 2:
secondary structure mountain plots of height vs. position (height mk = number of base pairs enclosing the base at position k) for the functional and pseudogenic sequences of the 5.8S gene of representative Triatominae taxa: A: functional conserved gene in Triatomini; B: pseudogenic sequence in Triatoma dimidiata dimidiata from Guatemala; C: pseudogenic sequence in T. dimidiata capitata from Colombia; D: pseudogenic sequence in T. sp. aff. dimidiata sensu (Bargues et al. 2008) from Yucatan, Mexico; E: pseudogenic sequence in Triatoma phyllosoma phyllosoma from Mexico; F: pseudogenic sequence in Triatoma sanguisuga from the United States of America. Centroid: centroid structure; mfe: minimum free energy structure; pf: partition function structure.

At ITS-2, the secondary structures predicted were showing considerably less negative values for mfe and centroid secondary structure predictions (-80.10 to -95.70 kcal/mol and from -37.85 to -79.80 kcal/mol, respectively). Large discrepancy was also observed in the values of ensemble diversity with regard to the same parameters for the functional ITS-2 (Table II). The lower fit between the three secondary structure prediction approaches assayed (mfe, pf and centroid) became very evident (Fig. 3F-J) when compared with their perfect fit in the respective functional ITS-2 secondary structure according to species (Fig. 3A-E), although in given cases the result did not appear bad at all, as in the case of Triatoma di- midiata dimidiata H1 from Guatemala (Fig. 3G).

Fig. 3:
secondary structure mountain plots of height vs. position (height mk = number of base pairs enclosing the base at position k) for the functional and pseudogenic sequences of the internal transcribed spacer-2 of representative Triatominae taxa. Functional sequences (A-E) and pseudogenic sequences (F-J) in Triatoma phyllosoma phyllosoma from Mexico (A, F), Triatoma dimidiata dimidiata from Guatemala (B, G), Triatoma mexicana from Mexico (C, H), Triatoma nitida from Guatemala (D, I), Triatoma sanguisuga from the United States of America (E, J). Centroid: centroid structure; mfe: minimum free energy structure; pf: partition function structure.

Unique characteristics of the "ps(5.8S+ITS-2)" pseudogene

The combination of several characteristics of the "ps(5.8S+ITS-2)" paralogous sequence detected in triatomine vector species proved it to be a unique pseudogene (Bargues et al. 2014): found in the nuclear DNA instead of in the mtDNA genome; correspond to variants of parts of the nuclear rDNA operon instead of mitochondrial sequences inserted in the nuclear genome (Numts); sequence copies identical inside the same host individual; sequence copies varying inside the same host species; found in an animal, in which such nuclear rDNA pseudogenes appear to be very rare when compared to their detection in plants (Márquez et al. 2003, Manen 2004Manen JF 2004. Are both sympatric species Ilex perado and Ilex canariensis secretly hybridizing? Indication from nuclear markers collected in Tenerife. BMC Evol Biol 4: 46., Razafimandimbison et al. 2004Razafimandimbison SG, Kellogg E, Bremer B 2004. Recent origin and phylogenetic utility of divergent ITS putative pseudogenes: a case study from Naucleeae (Rubiaceae). Syst Biol 53: 177-192., Roselló et al. 2007, Thornill et al. 2007Thornill DJ, Lajeunesse T, Santos SR 2007. Measuring rDNA diversity in eukaryotic microbial systems: how intragenomic variation, pseudogenes and PCR artifacts confound biodiversity estimates. Mol Ecol 16: 5326-5340., Zheng et al. 2008, Xu et al. 2009, Xiao et al. 2010Xiao LQ, Möller M, Zhu H 2010. High nrDNA ITS polymorphism in the ancient extant seed plant Cycas: incomplete concerted evolution and the origin of pseudogenes. Mol Phylogenet Evol 55: 168-177., Glass et al. 2013Glass DJ, Takebayashi N, Olson LE, Taylor DL 2013. Evaluation of the authenticity of a highly novel environmental sequence from boreal forest soil using ribosomal RNA secondary structure modeling. Mol Phylogenet Evol 67: 234-245.); fourth report of a pseudogene in the nuclear ribosomal operon of an animal, after two in insects and one in a vertebrate; first time a pseudogene is found in many different species (13 triatomine taxa) belonging to the same taxonomic group (and not in only 1 species, as in all pseudogenes found in other animal organisms and the majority of those found in plants - only a very few found in a few closely related plant species); not present in all species of Triatominae, but only in given species of the Triatomini tribe; interspecific variability following coherent evolutionary rates; present in phylogenetically closely related triatomine species sharing a geographical distribution, allowing for the assessment of a set of evolutionary characteristics in a pseudogene for the first time.

The discovery of a pseudogene in many phylogenetically related species is unique in animals and allowed for an estimation of its palaeobiogeographical origin based on molecular clock data, inheritance pathways, evolutionary rate and pattern and geographical spread (Bargues et al. 2014).

Pseudogene distribution - The presence of both functional and pseudogenic variants of the 5.8S+ITS-2 has been confirmed at least in the following triatomine species (Bargues et al. 2014): Meccus dimidiata dimidiata, M. d. capitata, M. d. maculipennis, M. d. hegneri, M. sp. aff. dimidiata, Meccus phyllosoma phyllosoma, M. p. longipennis, M. p. pallidipennis, M. p. picturata, M. p. mazzottii, Triatoma mexicana, Triatoma nitida and Triatoma sanguisuga.

This pseudogene, paralogous to the functional rDNA 5.8S and ITS-2, has been shown to be absent in the following triatomine species (Bargues et al. 2014): Triatoma gerstaeckeri, Triatoma barberi, Triatoma rubida, Triatoma ryckmani, Triatoma maculata, Triatoma infestans, Dipetalogaster maxima, Panstrongylus megistus and Rhodnius prolixus. The species Triatoma bassolsae, in a specimen of which the pseudogene was not found, should be reassessed due to problems posed in the systematic classification of the specimen (Bargues et al. 2014).

This means that the "ps(5.8S+ITS-2)" pseudogene has followed an evolution enabling it to cover a geographical distribution throughout North America, Central America and northern South America [see Carcavallo et al. (1999)Carcavallo RU, de Casas SIC, Sherlock IA, Giron IG, Jurberg J, Galvão C, Segura CAM, Noireau F 1999. Geographical distribution and alti-latitudinal dispersion. In RU Carcavallo, IG Girón, J Jurberg, H Lent, Atlas of Chagas disease vectors in the Americas, Editora Fiocruz, Rio de Janeiro, p. 747-792.regarding the distribution of the triatomine species presenting the pseudogene].

The molecular clock calculations, based on the evolutionary rates in triatomines (Bargues et al. 2000) indicated a relict pseudogene of a very ancient origin with an estimation for the appearance of the pseudogenic sequence of around 11-25.7 million years. Consequently, this additionally suggests that species with older origins or species not related to those in which the pseudogene has been found, should not present the pseudogene. The present geographical distribution of the triatomine taxa presenting the pseudogene suggests, moreover, that the origin of this pseudogene should have taken place in a triatomine ancestor living in the aforementioned old period in Mexico.

Functionality and applications - In eukaryotic cells, the vital role of rRNA molecules in protein synthesis leads to strong selection pressure to maintain functional rRNA molecules. The 18S, 5.8S and 28S of the rRNA cistron are produced by a series of nucleolytic reactions that liberate the mature rRNAs from the large primary precursor transcript synthesised by RNA polymerase 1. The three genes are transcribed as a unit by the RNA polymerase 1 into a single precursor pre-rRNA molecule, which undergoes a series of processing steps resulting in mature and fully functional rRNA. These processes include excision of the transcribed spacer regions, nucleotide modifications such as methylation and pseudoridylation, terminal additions of nucleotides and further cleavages and trimming of the precursor molecule (Venema & Tollervey 1999Venema J, Tollervey D 1999. Ribosome synthesis in Saccharomyces cerevisiae. Annu Rev Genet 33: 261-311.). The excision of the ITS regions is necessary for the production of mature rRNA molecules. Post-transcriptional processes split the cistron, removing two ITS. Thus, ITS-1 and ITS-2 are not incorporated into the structure of the ribosome and are relatively more variable and have less homology. Consequently, ITS-1 and ITS-2 are not subject to the same functional constraints as the rRNA genes and are therefore subject to higher evolutionary rates leading to greater variability in both nucleotide sequence and length.

Spacers play a role in the processing reactions: ITS-1 is required for the processing of the 3' end of the 18S and the 5' end of the 5.8S molecules, while ITS-2 is required for the processing of the 3' end of the 5.8S and the 5' end of the 28S molecules. For this processing of rRNA transcript to be correct and efficient, the proper folding of the ITSs appears to be essential (Côté & Peculis 2001). The importance of ITS secondary structure in rRNA processing is supported by the fact that while the primary sequence of these regions can be quite variable, the secondary structures they form are relatively well conserved across eukaryotes. Even minor modifications to these spacers can inhibit or prevent the formation of mature rRNA products (van Nues et al. 1994van Nues RW, Rientjes JM, van der Sande CAFM, Zerp SF, Sluiter C, Venema J, Planta RJ, Raue HA 1994. Separate structural elements within internal transcribed spacer 1 of Saccharomyces cerevisiae precursor ribosomal RNA direct the formation of 17S and 26S rRNA. Nucleic Acids Res 22: 912-919.).

ITS-2 plays an essential role in the maturation of the pre-RNA as its secondary structures, acquired shortly after transcription, contain the cleavage sites and secondary structure motifs recognised by the enzymatic complexes that act in the processing of the pre-RNA. Factors as biochemical events that disrupt the correct assemblage of its secondary structure have been shown to cause a decrease in the amounts or complete absence of mature 28S rRNA (Hunter et al. 2007Hunter RL, Lajeunese TC, Santos SR 2007. Structure and evolution of the rDNA internal transcribed spacer (ITS) region 2 in the symbiotic dinoflagellates (Symbiodinium, Dinophyta). J Phycol 43: 120-128.). Specific positions within ITS-2 have been identified as sites where processing reactions occur that generate the 7S pre-rRNA and the mature 28S rRNA (Venema & Tollervey 1999). Nucleotides surrounding processing sites are important as substitutions either decrease or eliminate 28S rRNA maturation or lead to decreased cell growth rates. Furthermore, mutations disrupting the secondary structures in the ITS regions may also reduce or eliminate the production of precursor molecules and mature rRNA products (van Nues et al. 1994, Côté & Peculis 2001).

In triatomines, the complete characterisation of the paralogous sequence in question indicated a processed or retrotransposed pseudogenic unit, evolving independently from the concerted evolution acting on the functional normal rDNA operon and following its own concerted evolution. Aspects indicating that this pseudogene follows its own concerted evolution are: (i) quantification demonstrated that there are many copies of the pseudogene inside a triatomine specimen, (ii) cloning allowed to verify that all pseudogene copies have an identical sequence inside the same specimen (no intraindividual pseudogene variability was ever detected), (iii) comparison of different triatomines showed that the pseudogene sequence evolves coherently according to the species phylogeny and (iv) sequence analyses demonstrated that the pseudogene follows an evolutionary rate similar to that of the functional ITS-2 (Bargues et al. 2014). Processed or retrotransposed pseudogenes present distinguishable characteristics indicating RNA processing, including: (i) lack of noncoding intervening sequences (intronic regions and promoters), (ii) presence of poly-A tracts at the carboxyl (3') end and (iii) homological extensions, i.e., flanking repeat regions which are associated with insertion sites of transposable elements (Gibson 1994Gibson LJ 1994. Pseudogenes and origins. Origins 21: 91-108., Cooper 1999Cooper DN 1999. Pseudogenes and their formation. In DN Cooper, Human gene evolution, BIOS Scientific Publishers, Oxford, p. 265-285. ).

The analyses of the triatomine paralogous rDNA sequences also indicated it to be a ghost pseudogene, with some intermediate functionality of the type of a regulatory function or transcriptional activity, at least with regard to the processing of the mature 28S rRNA (Bargues et al. 2014). This would be the case of course if the pseudogene does not expand also including the totality of the large subunit or a part sufficiently long as for truncation and consequent process abolishing. If the pseudogenic ITS-2 would somehow allow for the processing of the mature 28S, this would perhaps be at the origin of changes in behaviour and adaptation capacity, as seen with T. infestansITS-1 minisatellites (Bargues et al. 2006). The fact that divergence differences between functional ITS-2 and pseudogenic ITS-2 did not become higher, despite the long evolutionary period elapsed since the origin of the pseudogene, may suggest that the pseudogene is subject to some constraints instead of evolving free by neutral drift (Bargues et al. 2014).

Fundamental implications in triatomine research

The existence of two paralogous complete sequences of the ITS-2 inside the same triatomine specimen may be a technical problem easily overlooked in direct sequencing, because the lower stability of the secondary structures of pseudogenic sequences causes them to be preferentially amplified under standard PCR conditions (Harpke & Peterson 2006Harpke D, Peterson A 2006. Non-concerted ITS evolution in Mammillaria (Cactaceae). Mol Phylogenet Evol 41: 579-593., Zheng et al. 2008). This may happen when applying somewhat lower annealing temperatures, as widely applied when the same primers are used for a large set of different species. The risk for erroneous sequences represented by the presence of a pseudogenic sequence becomes an important potential problem, given that the ITS-2 is a crucial molecular tool for specimen classification purposes (Mas-Coma & Bargues 2009).

Additional to the aforementioned technical risk to be considered henceforth, this relict pseudogene proves to be a valuable marker for specimen classification, phylogenetic analyses and systematic/taxonomic studies. In the few organisms in which a pseudogene has proved to be relict and to show a substitution rate not outpacing speciation, they have shown great usefulness in reconstructing well supported phylogenetic trees, as verified in pseudogenes with a sequence divergence between them and their corresponding functional copies of about 15% (Zhang et al. 2008). In triatomines, the "ps(5.8S+ITS-2)" pseudogene fulfils the needed requirements: (i) an old origin, (ii) divergences observed between functional sequences and between pseudogenic sequences similar in average and (iii) sufficient resolution (Bargues et al. 2014).

The "ps(5.8S+ITS-2)" pseudogene sequence shows the highest resolution, even at the lowest taxonomic level of subspecies (Bargues et al. 2014): when the same pseudogene was found in specimens presenting different functional ITS-2 haplotype, they proved to belong to the same subspecies; when the same pseudogene was found in different subspecies, in all cases these were subspecies taxa whose validity was already previously put in question by the results obtained in previous functional ITS-2 analyses.

Regarding utility for phylogenetic analyses, the monophyletism of the "ps(5.8S+ITS-2)" pseudogene sequence in a together tree with functional sequences suggests that it no longer interacts genetically with functional copies and thus could be used as another completely new resource to infer phylogeny. Moreover, exhaustive comparative studies made showed that topologies obtained by pseudogenic sequences are more congruent than topologies obtained by functional sequences (Bargues et al. 2014): values of both the consistence index (CI) and the retention index (RI) proved to be very high in the trees reconstructed with pseudogenic sequences, indicating very high congruence of the topologies obtained; for the same paralogous sequences (whether 5.8S+ITS-2 or ITS-2 alone), values of both CI and RI in the trees of the pseudogenic sequences were always higher than those calculated for the trees of the functional sequences; support values for the evaluation of the reliability of the nodes appear to be similar when using the pseudogenic 5.8S+ITS-2 sequence and the functional ITS-2 sequence, sometimes even somewhat higher supports when using the pseudogene.

The comparison of the phylogenetic trees obtained when only using the functional ITS-2 sequence with the tree obtained with the pseudogenic 5.8S+ITS-2 sequence proved that both furnish pronouncedly similar topologies. However, the compared analysis showed some small differences at the level of the nodes which furnished interesting additional information which may be valuable to assess the relationships between species, allowing for the clarification of taxonomic arrangements. The value of the "ps(5.8S+ITS-2)" pseudogene for such purposes should be consequently highlighted. Despite a high indel number, low mutation number and an evolutionary rate similar to the functional ITS-2, that pseudogene distinguishes different taxa and furnishes coherent phylogenetic topologies with resolution similar to the functional ITS-2 (Bargues et al. 2014).

Examples of the taxonomic usefulness of the "ps(5.8S+ITS-2)" pseudogene sequence include (Bargues et al. 2014): (i) the monophyletism it shows for the Phyllosoma complex species and, hence, the support it gives to the validity of the genus Meccus, (ii) possibility for a new assessment of the present arrangement of triatomine species in different complexes mainly based on phenotypic characteristics (Schofield & Galvão 2009Schofield CJ, Diotaiuti L, Dujardin JP 1999. The process of domestication in Triatominae. Mem Inst Oswaldo Cruz 94 (Suppl. I): 375-378. ), (iii) confirmation of taxonomic entity at species level, as for instance in the case of M. sp. aff. dimidiata (Bargues et al. 2008) and (iv) assessment capacity at the level of subspecies, as for instance within the original Phyllosoma subcomplex taxa M. p. phyllosoma, M. p. longipennis, M. p. pallidipennis, M. p. picturata and M. p. mazzottii (Mas-Coma & Bargues 2009).

Potential applications for Chagas disease

All data obtained suggest that the "ps(5.8S+ITS-2)" pseudogene should have functionality, as to understand why it has been kept by inheritance throughout triatomine lineages following coherent evolutionary rates and patterns during such a long evolutionary period since its old origin, instead of being lost. Additionally, intraindividual and interpopulational as well as interspecific differences in the amount of the pseudogene, both relative levels and functional ITS-2/pseudogenic ITS ratios assessed by real time PCR, pose a question mark (Bargues et al. 2014). The fact that it is a pseudogene in the nuclear rDNA suggests possible relationships with fitness, behaviour and adaptability of the triatomine vectors, given the crucial role of rRNA molecules in protein synthesis.

It should be considered that nuclear rDNA spacer sequences are noncoding regions typically containing mini and microsatellite repeats (Mas-Coma & Bargues 2009). In that sense, recent studies refer to the involvement of microsatellite polymorphism in the social behaviour of animals (Hammock & Young 2005Hammock EAD, Young LJ 2005. Microsatellite instability generates diversity in brain and sociobehavioral traits. Science 308: 1630-1634.). Moreover, mini and microsatellites were early suspected to form secondary structures that may play an important role in the mutational process due to their repetitive nature and highly biased nucleotide composition. Expanded repeats in noncoding regions interfere with the metabolism of several cellular pathways, such as methylation, transcription, splicing, RNA processing, nuclear export and translation and the resulting expanded mRNAs often acquire an altered function (Richard et al. 2008Richard GF, Kerrest A, Dujon B 2008. Comparative genomics and molecular dynamics of DNA repeats in eukaryotes. Microbiol Mol Biol Rev 72: 686-727.).

The transmission of Chagas disease is mainly related to triatomine species adapted to live within human dwellings. An increasing number of species seems to be following a similar adaptive route from sylvatic to domestic habitats (Schofield et al. 1999). The understanding of the capacity of triatomines to colonise the domicile and inherent biological and genetic processes related to this domiciliation process from the wild is of considerable importance in relation to epidemiological surveillance and control of Chagas disease vectors (WHO 1991WHO - World Health Organization 1991. Control of Chagas disease: report of a WHO expert committee. Available from: whqlibdoc.who.int/trs/WHO_TRS_811.pdf.
whqlibdoc.who.int/trs/WHO_TRS_811.pdf...
, Dujardin 1998Dujardin JP 1998. Population genetics and the natural history of domestication in Triatominae. Mem Inst Oswaldo Cruz 93 (Suppl. II): 34-36. ).

DNA techniques are genetic tools which have been used to help in the understanding of the domiciliation capacity of triatomine species (Bargues et al. 2002), not only from de novo phenomena in triatomine species nowadays beginning to show a domiciliation trend in different areas, but also in the endeavour of forecasting human habitat re-colonisation by the same previously existing or other sylvatic triatomine species once the effect of insecticide spraying disappears, after control campaigns (Dujardin et al. 2000).

Studies on the main South American vector species T. infestans showed a large variation of the haploid DNA content, including a strikingly high DNA content difference between Andean sylvatic populations and non-Andean intradomiciliated specimens (mean reduction of 30%, with a maximum of up to 40%) and a correlation between presence/absence of minisatellites and larger/smaller genome size (Bargues et al. 2006). A relationship between total DNA content per cell, C-heterochromatin revealed by C-banding of chromosomes and highly repetitive DNA sequences has been suggested (Panzera et al. 2004Panzera F, Dujardin JP, Nicolini P, Coraccio MN, Rose V, Tellez T, Bermudez H, Bargues MD, Mas-Coma S, O'Connor JE, Perez R 2004. Genomic changes of Chagas disease vector, South America. Emerg Infect Dis 10: 438-446.).

The correlation which was observed between the presence or absence of minisatellites in the ITS-1 spacer and genome size is worth mentioning. In Bolivia, both the higher number of minisatellite repeats and the intraspecific variability in the ITS-1 length related to minisatellites agreed with the higher total DNA content and the intraspecific variability of genome size detected in Bolivian populations, which appeared to be larger than the variability found in the numerous populations studied from all other countries (Bargues et al. 2006). Thus, minisatellite repeats in ITS-1 proved to have a parallelism with sylvatic or intradomiciliary populations.

Concluding remarks

The results of the aforementioned analyses indicate that the "ps(5.8S+ITS-2)" pseudogene constitutes a new valuable marker useful in fundamental and applied studies, for specimen classification, phylogenetic analyses and systematic/taxonomic studies in triatomine vectors, as well as for Chagas disease epidemiology and control, thus opening a new broad research field on the disease throughout North, Central and northern South America.

ACKNOWLEDGEMENTS

This work benefited from international collaboration through the ECLAT network. All numerous specialists having collaborated in triatomine collections or having furnished triatomine specimens from different countries are greatly acknowledged.

REFERENCES

  • Balakirev ES, Ayala FJ 2003. Pseudogenes: are they "junk" or functional DNA? Annu Rev Genet 237: 123-151.
  • Bargues MD, Klisiowicz DR, Gonzalez-Candelas F, Ramsey JM, Monroy C, Ponce C, Salazar-Schettino PM, Panzera F, Abad-Franch F, Sousa OE, Schofield CJ, Dujardin JP, Guhl F, Mas-Coma S 2008. Phylogeography and genetic variation of Triatoma dimidiata, the main Chagas disease vector in Central America and its position within the genus Triatoma. PLoS Negl Trop Dis 2: e233.
  • Bargues MD, Klisiowicz DR, Panzera F, Noireau F, Marcilla A, Perez R, Rojas MG, O'Connor JE, Gonzalez-Candelas F, Galvão C, Jurberg J, Carcavallo RU, Dujardin JP, Mas-Coma S 2006. Origin and phylogeography of the Chagas disease main vector Triatoma infestans based on nuclear rDNA sequences and genome size. Infect Genet Evol 6: 46-62.
  • Bargues MD, Marcilla A, Dujardin JP, Mas-Coma S 2002. Triatomine vectors of Trypanosoma cruzi: a molecular perspective based on nuclear ribosomal DNA markers. Trans R Soc Trop Med Hyg 96: 159-164.
  • Bargues MD, Marcilla A, Ramsey JM, Dujardin JP, Schofield CJ, Mas-Coma S 2000. Nuclear rDNA-based molecular clock of the evolution of Triatominae (Hemiptera: Reduviidae), vectors of Chagas disease. Mem Inst Oswaldo Cruz 95: 567-573.
  • Bargues MD, Zuiaga MA, Mas-Coma S 2014. Nuclear rDNA pseudogenes in Chagas disease vectors: evolutionary implications of a new 5.8S+ITS-2 paralogous sequence marker in triatomines of North, Central and Northern South America. Infect Genet Evol 21: 134-156.
  • Carcavallo RU, de Casas SIC, Sherlock IA, Giron IG, Jurberg J, Galvão C, Segura CAM, Noireau F 1999. Geographical distribution and alti-latitudinal dispersion. In RU Carcavallo, IG Girón, J Jurberg, H Lent, Atlas of Chagas disease vectors in the Americas, Editora Fiocruz, Rio de Janeiro, p. 747-792.
  • Cooper DN 1999. Pseudogenes and their formation. In DN Cooper, Human gene evolution, BIOS Scientific Publishers, Oxford, p. 265-285.
  • Côté CA, Peculis BA 2001. Role of the ITS2-proximal stem and evidence for indirect recognition of processing sites in pre-rRNA processing in yeast. Nucleic Acids Res 29: 2106-2116.
  • Dujardin JP 1998. Population genetics and the natural history of domestication in Triatominae. Mem Inst Oswaldo Cruz 93 (Suppl. II): 34-36.
  • Dujardin JP, Schofield CJ, Panzera F 2000. Les vecteurs de la maladie de Chagas. Recherches taxonomiques, biologiques et génétiques, Académie Royale des Sciences d'Outre Mer, Buxelles, 189 pp.
  • Friedberg F, Rhoads AR 2000. Calculation and verification of the ages of retroprocessed pseudogenes. Mol Phylogenet Evol 16: 127-130.
  • Gibson LJ 1994. Pseudogenes and origins. Origins 21: 91-108.
  • Glass DJ, Takebayashi N, Olson LE, Taylor DL 2013. Evaluation of the authenticity of a highly novel environmental sequence from boreal forest soil using ribosomal RNA secondary structure modeling. Mol Phylogenet Evol 67: 234-245.
  • Gray TA, Wilson A, Fortin PJ, Nicholls RD 2006. The putatively functional Mkrn1-p1 pseudogene is neither expressed nor imprinted, nor does it regulate its source gene in trans. Proc Natl Acad Sci USA 103: 12039-12044.
  • Hammock EAD, Young LJ 2005. Microsatellite instability generates diversity in brain and sociobehavioral traits. Science 308: 1630-1634.
  • Harpke D, Peterson A 2006. Non-concerted ITS evolution in Mammillaria (Cactaceae). Mol Phylogenet Evol 41: 579-593.
  • Harrison PM, Zheng D, Zhang Z, Carriero N, Gerstein M 2005. Transcribed processed pseudogenes in the human genome: an intermediate form of expressed retrosequence lacking protein-coding ability. Nucleic Acids Res 33: 2374-2383.
  • Hofacker IL, Fontana W, Stadler PF, Bonhoeffer SL, Tacker M, Schuster P 1994. Fast folding and comparison of RNA secondary structures. Monatsh Chem 125: 167-188.
  • Hofacker IL, Stadler PF 2008. RNA secondary structures. In T Lengauer, Bioinformatics - from genomes to therapies, Wiley-VCH Verlag GmbH, Weinheim, p. 439-489.
  • Hunter RL, Lajeunese TC, Santos SR 2007. Structure and evolution of the rDNA internal transcribed spacer (ITS) region 2 in the symbiotic dinoflagellates (Symbiodinium, Dinophyta). J Phycol 43: 120-128.
  • Kalinin MI, Kononogov SA 2005. Boltzmann's constant, the energy meaning of temperature and thermodynamic irreversibility. Meas Tech 48: 632-636.
  • Keller I, Chintauan-Marquier IC, Veltsos, P, Nichols RA 2006. Ribosomal DNA in the grasshopper Podisma pedestris: escape from concerted evolution. Genetics 174: 863-874.
  • Li C, Wilkerson RC 2007. Intragenomic rDNA ITS2 variation in the Neotropical Anopheles (Nyssorhynchus) albitarsis complex (Diptera: Culicidae). J Hered 98: 51-59.
  • Lorenz WA, Clote P 2011. Computing the partition function for kinetically trapped RNA secondary structures. PLoS ONE 6: e16178.
  • Manen JF 2004. Are both sympatric species Ilex perado and Ilex canariensis secretly hybridizing? Indication from nuclear markers collected in Tenerife. BMC Evol Biol 4: 46.
  • Márquez LM, Miller DJ, Mackenzie JB, Van Oppen MJM 2003. Pseudogenes contribute to extreme diversity of nuclear ribosomal DNA in the hard coral Acropora. Mol Biol Evol 20: 1077-1086.
  • Mas-Coma S, Bargues MD 2009. Populations, hybrids and the systematic concepts of species and subspecies in Chagas disease triatomine vectors inferred from nuclear ribosomal and mitochondrial DNA. Acta Trop 110: 112-136.
  • Mathews DH, Disney MD, Childs JL, Schroeder S, Zuker M, Turner DH 2004. Incorporating chemical modification constraints into a dynamic programming algorithm for prediction of RNA secondary structure. Proc Natl Acad Sci USA 101: 7287-7292.
  • Panzera F, Dujardin JP, Nicolini P, Coraccio MN, Rose V, Tellez T, Bermudez H, Bargues MD, Mas-Coma S, O'Connor JE, Perez R 2004. Genomic changes of Chagas disease vector, South America. Emerg Infect Dis 10: 438-446.
  • Pavlicek A, Gentles AJ, Paces J, Paces V, Jurka J 2006. Retroposition of processed pseudogenes: the impact of RNA stability and translational control. Trends Genet 22: 69-73.
  • Piehler AP, Hellum M, Wenzel JJ, Kaminsk E, Haug KBF, Kierulf P, Kaminski WE 2008. The human ABC transporter pseudogene family: evidence for transcription and genepseudogene interference. BMC Genomics 9: 165.
  • Razafimandimbison SG, Kellogg E, Bremer B 2004. Recent origin and phylogenetic utility of divergent ITS putative pseudogenes: a case study from Naucleeae (Rubiaceae). Syst Biol 53: 177-192.
  • Richard GF, Kerrest A, Dujon B 2008. Comparative genomics and molecular dynamics of DNA repeats in eukaryotes. Microbiol Mol Biol Rev 72: 686-727.
  • Roselló JA, Lázaro A, Cosín R, Molins A 2007. Phylogeographic split in Buxus balearica (Buxaceae) as evidenced by nuclear ribosomal markers: when ITS paralogues are welcome. J Mol Evol 64: 143-157.
  • Rouchka EC, Cha IE 2009. Current trends in pseudogene detection and characterization. Curr Bioinform 4: 112-119.
  • Sasidharan R, Gerstein M 2008. Genomics: protein fossils live on as RNA. Nature 453: 729-731.
  • Schofield CJ, Diotaiuti L, Dujardin JP 1999. The process of domestication in Triatominae. Mem Inst Oswaldo Cruz 94 (Suppl. I): 375-378.
  • Schofield CJ, Galvão C 2009. Classification, evolution and species groups within the Triatominae. Acta Trop 110: 88-100.
  • Tam OH, Aravin AA, Stein P, Girard A, Murchison EP, Cheloufi S, Hodges E, Anger M, Sachidanandam R, Schultz RM, Hannon GJ 2008. Pseudogene-derived small interfering RNAs regulate gene expression in mouse oocytes. Nature 453: 534-538.
  • Thirumalai D, Lee N, Woodson SA, Klimov D 2001. Early events in RNA folding. Annu Rev Phys Chem 52: 751-762.
  • Thornill DJ, Lajeunesse T, Santos SR 2007. Measuring rDNA diversity in eukaryotic microbial systems: how intragenomic variation, pseudogenes and PCR artifacts confound biodiversity estimates. Mol Ecol 16: 5326-5340.
  • van Nues RW, Rientjes JM, van der Sande CAFM, Zerp SF, Sluiter C, Venema J, Planta RJ, Raue HA 1994. Separate structural elements within internal transcribed spacer 1 of Saccharomyces cerevisiae precursor ribosomal RNA direct the formation of 17S and 26S rRNA. Nucleic Acids Res 22: 912-919.
  • Venema J, Tollervey D 1999. Ribosome synthesis in Saccharomyces cerevisiae. Annu Rev Genet 33: 261-311.
  • WHO - World Health Organization 1991. Control of Chagas disease: report of a WHO expert committee. Available from: whqlibdoc.who.int/trs/WHO_TRS_811.pdf.
    » whqlibdoc.who.int/trs/WHO_TRS_811.pdf
  • Xiao LQ, Möller M, Zhu H 2010. High nrDNA ITS polymorphism in the ancient extant seed plant Cycas: incomplete concerted evolution and the origin of pseudogenes. Mol Phylogenet Evol 55: 168-177.
  • Xu J, Zhang Q, Xu X, Wanga Z, Qi J 2009. Intragenomic variability and pseudogenes of ribosomal DNA in stone flounder Kareius bicoloratus. Mol Phylogenet Evol 52: 157-166.
  • Zhang Z, Carriero N, Gerstein M 2004. Comparative analysis of processed pseudogenes in the mouse and human genomes. Trends Genet 20: 62-67.
  • Zhang Z, Gerstein M 2003. The human genome has 49 cytochrome c pseudogenes, including a relic of a primordial gene that still functions in mouse. Gene 312: 61-72.
  • Zhang ZD, Cayting P, Weinstock G, Gerstein M 2008. Analysis of nuclear receptor pseudogenes in vertebrates: how the silent tell their stories. Mol Biol Evol 25: 131-143.
  • Zheng D, Gerstein MB 2007. The ambiguous boundary between genes and pseudogenes: the dead rise up or do they? Trends Genet 23: 219-224.
  • Zheng D, Zhang Z, Harrison PM, Karro J, Carriero N, Gerstein M 2005. Integrated pseudogene annotation for human chromosome 22: evidence for transcription. J Mol Biol 349: 27-45.
  • Zheng X, Cai D, Yao L, Teng Y 2008. Non-concerted ITS evolution, early origin and phylogenetic utility of ITS pseudogenes in Pyrus. Mol Phylogenet Evol 48: 892-903.
  • Zuker M 1989. On finding all suboptimal foldings of an RNA molecule. Science 244: 48-52.
  • Zuker M, Stiegler P 1981. Optimal computer folding of large RNA sequences using thermodynamics and auxiliary information. Nucleic Acids Res 9: 133-148.
  • Financial support: RICET (RD12/0018/0013 of RETICS), VI National Plan of I+D+I 2008-2011, ISCIII - Subdirección General de Redes y Centros de Investigación Coopetariva and FEDER, Ministry of Health, Spain, PROMETEO (2012/042), Programa de Ayudas para Grupos de Investigación de Excelencia, Generalitat Valenciana, Spain

Publication Dates

  • Publication in this collection
    06 Mar 2015
  • Date of issue
    May 2015

History

  • Received
    28 Oct 2014
  • Accepted
    10 Feb 2015
Instituto Oswaldo Cruz, Ministério da Saúde Av. Brasil, 4365 - Pavilhão Mourisco, Manguinhos, 21040-900 Rio de Janeiro RJ Brazil, Tel.: (55 21) 2562-1222, Fax: (55 21) 2562 1220 - Rio de Janeiro - RJ - Brazil
E-mail: memorias@fiocruz.br