Transcriptome profiling and Calreticulin expression in Zika virus -infected Aedes aegypti.

Zika virus (ZIKV) may cause febrile illness and neurological damage, such as microcephaly in fetuses. ZIKV is transmitted to humans by Aedes aegypti, a nearly cosmopolitan mosquito. Understanding the virus-vector molecular interactions has been promising to enhance the knowledge towards disease mitigation. Since ZIKV infection alters gene physiology of mosquitoes, we examined the expression profile of ZIKV-infected Ae. aegypti by several approaches to identify genes altered by viral infection. Transcriptomics were performed by comparing between ZIKV-infected and uninfected Ae. aegypti females, which revealed some differentially expressed genes. Most of these genes appear to be involved with immune response as evidenced by an interactome analysis, and a prominent finding was a calreticulin-like (CRT) gene, which was upregulated during the infection. Expression of CRT was also experimentally quantified by qPCR, however, it revealed no significant differences between infected and uninfected females. Instead, expression levels were highly variable among individuals and negatively correlated to viral load. We also tested the possibility of this gene to be silenced, but the double-stranded RNA did not reduce CRT expression, and actually increased the inter-individuals' expressional variability. Present results differed from our original hypothesis of upregulation by infection. They also diverged between them (comparing qPCR to Transcriptomics) and from the literature which reported augmented CRT levels in Aedes species during viral infection. Present case probably underlies a more complex virus-host interaction system than we expected. Regulation of this gene seems not to be a linear correlation between expression and viremy. As infection takes place, a complex homeostatic mechanism may act to prevent expression and other cellular tasks from drifting. It is also possible that CRT expression is simply randomly disturbed by viral infection. Taken together, results show that CRT expression profile during ZIKV infection is complex and requires different investigative approaches to be understood. Studies focused on the biochemical function of CRT protein and on its role in the native mosquito metabolic network could unravel how it is actually influenced by ZIKV. Current work contributes more by getting incidental findings and by posing new hypotheses than by answering the original questions.

Characterization of a Gene Coding for the Complement System Component FB from Loxosceles laeta Spider Venom Glands.

The human complement system is composed of more than 30 proteins and many of these have conserved domains that allow tracing the phylogenetic evolution. The complement system seems to be initiated with the appearance of C3 and factor B (FB), the only components found in some protostomes and cnidarians, suggesting that the alternative pathway is the most ancient. Here, we present the characterization of an arachnid homologue of the human complement component FB from the spider Loxosceles laeta. This homologue, named Lox-FB, was identified from a total RNA L. laeta spider venom gland library and was amplified using RACE-PCR techniques and specific primers. Analysis of the deduced amino acid sequence and the domain structure showed significant similarity to the vertebrate and invertebrate FB/C2 family proteins. Lox-FB has a classical domain organization composed of a control complement protein domain (CCP), a von Willebrand Factor domain (vWFA), and a serine protease domain (SP). The amino acids involved in Mg2+ metal ion dependent adhesion site (MIDAS) found in the vWFA domain in the vertebrate C2/FB proteins are well conserved; however, the classic catalytic triad present in the serine protease domain is not conserved in Lox-FB. Similarity and phylogenetic analyses indicated that Lox-FB shares a major identity (43%) and has a close evolutionary relationship with the third isoform of FB-like protein (FB-3) from the jumping spider Hasarius adansoni belonging to the Family Salcitidae.

Adaptive evolution in the toxicity of a spider's venom enzymes.

BACKGROUND: Sphingomyelinase D is the main toxin present in the venom of Loxosceles spiders. Several isoforms present in these venoms can be structurally classified in two groups. Class I Sphingomyelinase D contains a single disulphide bridge and variable loop. Class II Sphingomyelinase D presents an additional intrachain disulphide bridge that links a flexible loop with a catalytic loop. These classes exhibit differences in their toxic potential. In this paper we address the distribution of the structural classes of SMase D within and among species of spiders and also their evolutionary origin by means of phylogenetic analyses. We also conducted tests to assess the action of natural selection in their evolution combined to structural modelling of the affected sites. RESULTS: The majority of the Class I enzymes belong to the same clade, which indicates a recent evolution from a single common ancestor. Positively selected sites are located on the catalytic interface, which contributes to a distinct surface charge distribution between the classes. Sites that may prevent the formation of an additional bridge were found in Class I enzymes. CONCLUSIONS: The evolution of Sphingomyelinase D has been driven by natural selection toward an increase in noxiousness, and this might help explain the toxic variation between classes.

Network genealogy of 195-bp satellite DNA supports the superimposed hybridization hypothesis of Trypanosoma cruzi evolutionary pattern.

Trypanosoma cruzi is highly diverse genetically and has been partitioned into six discrete typing units (DTUs), recently re-named T. cruzi I-VI. Although T. cruzi reproduces predominantly by binary division, accumulating evidence indicates that particular DTUs are the result of hybridization events. Two major scenarios for the origin of the hybrid lineages have been proposed. It is accepted widely that the most heterozygous TcV and TcVI DTUs are the result of genetic exchange between TcII and TcIII strains. On the other hand, the participation of a TcI parental in the current genome structure of these hybrid strains is a matter of debate. Here, sequences of the T. cruzi-specific 195-bp satellite DNA of TcI, TcII, TcIII, TcV, and TcVI strains have been used for inferring network genealogies. The resulting genealogy showed a high degree of reticulation, which is consistent with more than one event of hybridization between the Tc DTUs. The data also strongly suggest that TcIII is a hybrid with two distinct sets of satellite sequences, and that genetic exchange between TcI and TcII parentals occurred within the pedigree of the TcV and TcVI DTUs. Although satellite DNAs belong to the fast-evolving portion of eukaryotic genomes, in >100 satellite units of nine T. cruzi strains we found regions that display 100% identity. No DTU-specific consensus motifs were identified, inferring species-wide conservation.

Trypanosoma cruzi: exploring the nuclear genome of zymodeme 3 stocks by chromosome size polymorphism.

Chagas disease is emerging in the Brazilian Amazon. We evaluated the position of eight zymodeme 3 isolates from Amazonian sylvatic vectors and one human case in relation to Trypanosoma cruzi I and II major groups and hybrid strains by chromosome size polymorphism. Nineteen isolates were analyzed by mapping nine coding sequences on chromosomal bands (0.6-3.3Mbp). Numerical analysis was based on the absolute chromosomal size difference index (aCSDI). A dendrogram was obtained applying the minimum evolution criterion and considering the aCSDI values to estimate the branch lengths. The isolates were distributed in four groups. Group A clustered hybrid isolates; Groups B and C, T. cruzi II and T. cruzi I isolates, respectively. Seven Z3 stocks were clustered in Group D, which showed low intra-group diversity and was the most divergent. The proportion of two different-sized homologous chromosomes was determined. Wild vectors harboring Z3 stocks constitute a potential reservoir of human infection in the Amazon.

Comparative analysis of genomic sequences suggests that Trypanosoma cruzi CL Brener contains two sets of non-intercalated repeats of satellite DNA that correspond to T. cruzi I and T. cruzi II types.

Approximately 10% of the Trypanosoma cruzi genome is formed by a satellite DNA, composed by 195-bp repeats organized in 30+/-10 kb clusters in some, but not all chromosomes. Here, the satellite DNA of six representative T. cruzi strains was sequenced and used for phylogenetic inference. The results show that CL Brener contains satellite repeats from T. cruzi I and T. cruzi II strains, although type II sequences are more abundant. The presence of types I and II sequences extends previous propositions that genetic exchange between the two major T. cruzi lineages have occurred in CL Brener, although our data accommodate alternative scenarios of hybridization within T. cruzi II, as proposed by others. Altogether, present data suggest a complex origin for CL Brener. Sequence analysis of satellites isolated from chromosomal bands indicates that satellite DNA sequences are not chromosome specific. Neighbor analysis of in tandem satellite DNAs containing up to five repeats shows that each cluster contains only one type of sequence. Consequently, clusters with intercalated types I and II repeats were not found. We propose that the CL Brener genome contains large pieces of satellite DNA originated mainly from chromosomes of T. cruzi II with introgression of T. cruzi I lineage.

Chromosomal polymorphism, gene synteny and genome size in T. cruzi I and T. cruzi II groups.

Pulsed-field gel electrophoresis and DNA hybridization were used to establish and compare some parameters of the molecular karyotype of nine stocks classified into Trypanosoma cruzi I and T. cruzi II groups. The isolates showed a variable number of chromosomal bands (17-22) comprised between 0.4 and 3.3 Mbp. The total number of chromosomes and the genome size were estimated based on the fluorescence intensity of SYBR Green I-stained chromosomal bands. Differences in the length of the telomeric regions among the stocks and between chromosomes of the same stock were observed. No correlation was found between the length of the telomeric region and the group to which the isolate belongs. Hybridization of 54 genetic markers revealed extensive chromosome size polymorphism. Nevertheless, the most represented pattern was the hybridization of the probes in larger chromosomes in stocks of T. cruzi II as compared to T. cruzi I. Eight putative syntenic groups, encompassing 29 non-redundant genetic markers and distributed in 11 CL Brener chromosomal bands were disclosed. The syntenic groups were conserved in all the stocks. The relative abundance of repetitive DNA sequences was determined. C6, B11/L1Tc and E12 elements presented maximum 1.7-fold variation in copy number, whereas 195-bp satellite DNA (120,000 copies in Y strain) was four- to nine-fold more abundant in T. cruzi II stocks. The novel aspects of T. cruzi karyotype here presented contribute to the comprehension of the genome organization of this parasite and will assist the assignment of scaffold to the CL Brener chromosomal bands.

Evaluation of Trypanosoma cruzi hybrid stocks based on chromosomal size variation.

Although all classical lines of evidence point to the fact that Trypanosoma cruzi has a predominantly clonal evolution, accumulating data show that some T. cruzi stocks are the result of hybridisation events. We evaluated whether chromosomal polymorphism would give evolutionary information on hybrid isolates. Twenty-three coding sequences were mapped on the chromosomes of nine parasite stocks, four of which are putative hybrids (CL Brener and rDNA group 1/2). Phenetic analyses of karyotype data were based on the absolute chromosomal size difference index (aCSDI), a method that assumes that the genomic distance between two organisms is the sum of the size differences between their homologous chromosomes. aCSDI-based dendrograms obtained from a variable number of probes (3-18 probes) defined in all the cases three clusters: two corresponding, respectively, to T. cruzi I and T. cruzi II groups; and a third one, to rDNA group 1/2. CL Brener was alternatively positioned in T. cruzi II or rDNA group 1/2 clusters. Three clusters were also observed in the dendrogram constructed with restriction fragment length polymorphism (RFLP) data from 18 probes. The topology of the chromosome and RFLP dendrograms is similar, with a significant correlation coefficient (r=0.86062; P<0.0001), supporting a strong structuring of the clusters. This study also revealed that hybrid stocks have a larger proportion of two different-sized homologous chromosomes, as compared with non-hybrid strains. Overall, our results show that chromosomes are valuable characters for identification of evolutionary groups, in particular, T. cruzi hybrid organisms.

Two types of ribosomal RNA genes in hybrid Trypanosoma cruzi strains.

Trypanosoma cruzi isolates can be divided into two major phylogenetic lineages-T. cruzi I and T. cruzi II. The population structure is predominantly clonal, with sexuality having no or limited influence on the evolution of the parasite. Isoenzymes and nuclear gene sequences have provided evidence that some T. cruzi strains are hybrids. Previous work of our group has shown that the putative hybrid strains designated as group 1/2 contain two types of rDNA units, corresponding to those found in T. cruzi I and T. cruzi II. In this study, the presence and transcription of the two types of ribosomal RNA (rRNA) cistrons were investigated in epimastigotes, metacyclic and tissue culture trypomastigotes of group 1/2 isolates. PCR and RT-PCR assays indicate that both types of cistrons are present in group 1/2 strains, but only type-2 genes are transcribed in all developmental stages. The structure of the promoter regions of group 1/2 was compared to reference T. cruzi I and T. cruzi II strains. In all cases, the transcription start point was mapped to a conserved A residue located approximately 1800 bp upstream the 18S rRNA gene. The distribution of rDNA clusters in chromosomal bands of group 1/2 was evaluated by pulsed-field gel electrophoresis (PFGE). The majority of type-2 rDNA genes are localized in a 1.5 Mbp band, whereas type-1 cistrons are mostly concentrated in a 1.1 Mbp band. The structural and functional studies of group 1/2 ribosomal cistrons described here may shed light on the evolutionary processes that took place during the generation of such hybrid organisms.