Phylogeny and evolution of the aspartyl protease family from clinically relevant Candida species

Aspartyl proteases are a class of enzymes that include the yeast aspartyl proteases and secreted aspartyl protease (Sap) superfamilies. Several Sap superfamily members have been demonstrated or suggested as virulence factors in opportunistic pathogens of the genus Candida. Candida albicans, Candida tropicalis, Candida dubliniensis and Candida parapsilosis harbour 10, four, eight and three SAP genes, respectively. In this work, genome mining and phylogenetic analyses revealed the presence of new members of the Sap superfamily in C. tropicalis (8), Candida guilliermondii (8), C. parapsilosis(11) and Candida lusitaniae (3). A total of 12 Sap families, containing proteins with at least 50% similarity, were discovered in opportunistic, pathogenic Candida spp. In several Sap families, at least two subfamilies or orthologous groups were identified, each defined by > 90% sequence similitude, functional similarity and synteny among its members. No new members of previously described Sap families were found in a Candida spp. clinical strain collection; however, the universality of SAPT gene distribution among C. tropicalis strains was demonstrated. In addition, several features of opportunistic pathogenic Candida species, such as gene duplications and inversions, similitude, synteny, putative transcription factor binding sites and genome traits of SAP gene superfamily were described in a molecular evolutionary context.

Characterisation of pvmdr1 and pvdhfr genes associated with chemoresistance in Brazilian Plasmodium vivax isolates

Plasmodium vivax control is now being hampered by drug resistance. Orthologous Plasmodium falciparum genes linked to chloroquine or sulfadoxine-pyrimethamine chemoresistance have been identified in P. vivax parasites, but few studies have been performed. The goal of the present work is to characterise pvmdr1 and pvdhfr genes in parasite isolates from a Brazilian endemic area where no molecular investigation had been previously conducted. The pvmdr1 analysis revealed the existence of single (85.7%) and double (14.3%) mutant haplotypes, while the pvdhfr examination showed the presence of double (57.2%) and triple (42.8%) mutant haplotypes. The implications of these findings are discussed.

Comparative genomics of ParaHox clusters of teleost fishes: gene cluster breakup and the retention of gene sets following whole genome duplications.

BACKGROUND: The evolutionary lineage leading to the teleost fish underwent a whole genome duplication termed FSGD or 3R in addition to two prior genome duplications that took place earlier during vertebrate evolution (termed 1R and 2R). Resulting from the FSGD, additional copies of genes are present in fish, compared to tetrapods whose lineage did not experience the 3R genome duplication. Interestingly, we find that ParaHox genes do not differ in number in extant teleost fishes despite their additional genome duplication from the genomic situation in mammals, but they are distributed over twice as many paralogous regions in fish genomes. RESULTS: We determined the DNA sequence of the entire ParaHox C1 paralogon in the East African cichlid fish Astatotilapia burtoni, and compared it to orthologous regions in other vertebrate genomes as well as to the paralogous vertebrate ParaHox D paralogons. Evolutionary relationships among genes from these four chromosomal regions were studied with several phylogenetic algorithms. We provide evidence that the genes of the ParaHox C paralogous cluster are duplicated in teleosts, just as it had been shown previously for the D paralogon genes. Overall, however, synteny and cluster integrity seems to be less conserved in ParaHox gene clusters than in Hox gene clusters. Comparative analyses of non-coding sequences uncovered conserved, possibly co-regulatory elements, which are likely to contain promoter motives of the genes belonging to the ParaHox paralogons. CONCLUSION: There seems to be strong stabilizing selection for gene order as well as gene orientation in the ParaHox C paralogon, since with a few exceptions, only the lengths of the introns and intergenic regions differ between the distantly related species examined. The high degree of evolutionary conservation of this gene cluster's architecture in particular - but possibly clusters of genes more generally - might be linked to the presence of promoter, enhancer or inhibitor motifs that serve to regulate more than just one gene. Therefore, deletions, inversions or relocations of individual genes could destroy the regulation of the clustered genes in this region. The existence of such a regulation network might explain the evolutionary conservation of gene order and orientation over the course of hundreds of millions of years of vertebrate evolution. Another possible explanation for the highly conserved gene order might be the existence of a regulator not located immediately next to its corresponding gene but further away since a relocation or inversion would possibly interrupt this interaction. Different ParaHox clusters were found to have experienced differential gene loss in teleosts. Yet the complete set of these homeobox genes was maintained, albeit distributed over almost twice the number of chromosomes. Selection due to dosage effects and/or stoichiometric disturbance might act more strongly to maintain a modal number of homeobox genes (and possibly transcription factors more generally) per genome, yet permit the accumulation of other (non regulatory) genes associated with these homeobox gene clusters.

Characterization of anti-silencing factor 1 in Leishmania major

Anti-silencing factor 1 (ASF1) is a histone chaperone that contributes to the histone deposition during nucleosome assembly in newly replicated DNA. It is involved in chromatin disassembly, transcription activation and in the cellular response to DNA damage. In Leishmania major the ASF1 gene (LmASF1) is located in chromosome 20 and codes for a protein showing 67% of identity with the Trypanosoma brucei TbASF1a. Compared to orthologous proteins, LmASF1 conserves the main residues relevant for its various biological functions. To study ASF1 in Leishmania we generated a mutant overexpressing LmASF1 in L. major. We observed that the excess of LmASF1 impaired promastigotes growth rates and had no impact on cell cycle progress. Differently from yeast, ASF1 overproduction in Leishmania did not affect expression levels of genes located on telomeres, but led to an upregulation of proteins involved in chromatin remodelling and physiological stress, such as heat shock proteins, oxidoreductase activity and proteolysis. In addition, we observed that LmASF1 mutant is more susceptible to the DNA damaging agent, methyl methane sulphonate, than the control line. Therefore, our study suggests that ASF1 from Leishmania pertains to the chromatin remodelling machinery of the parasite and acts on its response to DNA damage.

The level of ascorbate peroxidase is enhanced in benznidazole-resistant populations of Trypanosoma cruzi and its expression is modulated by stress generated by hydrogen peroxide

Ascorbate peroxidases (APX) are class I heme-containing enzymes that convert hydrogen peroxide into water molecules. The gene encoding APX has been characterized in 11 strains of Trypanosoma cruzi that are sensitive or resistant to benznidazole (BZ). Bioinformatic analysis revealed the presence of two complete copies of the T. cruzi APX (TcAPX) gene in the genome of the parasite, while karyotype analysis showed that the gene was present in the 2.000-kb chromosome of all of the strains analyzed. The sequence of TcAPX exhibited greater levels of similarity to those of orthologous enzymes from Leishmania spp than to APXs from the higher plant Arabidopsis thaliana. Northern blot and real-time reverse transcriptase polymerase chain reaction (RT-PCR) analyses revealed no significant differences in TcAPX mRNA levels between the T. cruzi strains analyzed. On the other hand, Western blots showed that the expression levels of TcAPX protein were, respectively, two and three-fold higher in T. cruzi populations with in vitro induced (17 LER) and in vivo selected (BZR) resistance to BZ, in comparison with their corresponding susceptible counterparts. Moreover, the two BZ-resistant populations exhibited higher tolerances to exogenous hydrogen peroxide than their susceptible counterparts and showed TcAPX levels that increased in a dose-dependent manner following exposure to 100 and 200 µM hydrogen peroxide.

Trypanosoma rangeli protein tyrosine phosphatase is associated with the parasite's flagellum

Protein tyrosine phosphatases (PTPs) play an essential role in the regulation of cell differentiation in pathogenic trypanosomatids. In this study, we describe a PTP expressed by the non-pathogenic protozoan Trypanosoma rangeli (TrPTP2). The gene for this PTP is orthologous to the T. brucei TbPTP1 and Trypanosoma cruzi (TcPTP2) genes. Cloning and expression of the TrPTP2 and TcPTP2 proteins allowed anti-PTP2 monoclonal antibodies to be generated in BALB/c mice. When expressed by T. rangeli epimastigotes and trypomastigotes, native TrPTP2 is detected as a ~65 kDa protein associated with the parasite's flagellum. Given that the flagellum is an important structure for cell differentiation in trypanosomatids, the presence of a protein responsible for tyrosine dephosphorylation in the T. rangeli flagellum could represent an interesting mechanism of regulation in this structure.

The zinc finger protein TcZFP2 binds target mRNAs enriched during Trypanosoma cruzi metacyclogenesis

Trypanosomes are parasitic protozoa in which gene expression is primarily controlled through the regulation of mRNA stability and translation. This post-transcriptional control is mediated by various families of RNA-binding proteins, including those with zinc finger CCCH motifs. CCCH zinc finger proteins have been shown to be essential to differentiation events in trypanosomatid parasites. Here, we functionally characterise TcZFP2 as a predicted post-transcriptional regulator of differentiation in Trypanosoma cruzi. This protein was detected in cell culture-derived amastigotes and trypomastigotes, but it was present in smaller amounts in metacyclic trypomastigote forms of T. cruzi. We use an optimised recombinant RNA immunopreciptation followed by microarray analysis assay to identify TcZFP2 target mRNAs. We further demonstrate that TcZFP2 binds an A-rich sequence in which the adenosine residue repeats are essential for high-affinity recognition. An analysis of the expression profiles of the genes encoding the TcZFP2-associated mRNAs throughout the parasite life cycle by microarray hybridisation showed that most of the associated mRNAs were upregulated in the metacyclic trypomastigote forms, also suggesting a role for TcZFP2 in metacyclic trypomastigote differentiation. Knockdown of the orthologous Trypanosoma brucei protein levels showed ZFP2 to be a positive regulator of specific target mRNA abundance.

Hominoid fission of chromosome 14/15 and the role of segmental duplications.

Ape chromosomes homologous to human chromosomes 14 and 15 were generated by a fission event of an ancestral submetacentric chromosome, where the two chromosomes were joined head-to-tail. The hominoid ancestral chromosome most closely resembles the macaque chromosome 7. In this work, we provide insights into the evolution of human chromosomes 14 and 15, performing a comparative study between macaque boundary region 14/15 and the orthologous human regions. We construct a 1.6-Mb contig of macaque BAC clones in the region orthologous to the ancestral hominoid fission site and use it to define the structural changes that occurred on human 14q pericentromeric and 15q subtelomeric regions. We characterize the novel euchromatin-heterochromatin transition region (∼20 Mb) acquired during the neocentromere establishment on chromosome 14, and find it was mainly derived through pericentromeric duplications from ancestral hominoid chromosomes homologous to human 2q14-qter and 10. Further, we show a relationship between evolutionary hotspots and low-copy repeat loci for chromosome 15, revealing a possible role of segmental duplications not only in mediating but also in "stitching" together rearrangement breakpoints.

Plasmodium vivax antigen discovery based on alpha-helical coiled coil protein motif.

Protein α-helical coiled coil structures that elicit antibody responses, which block critical functions of medically important microorganisms, represent a means for vaccine development. By using bioinformatics algorithms, a total of 50 antigens with α-helical coiled coil motifs orthologous to Plasmodium falciparum were identified in the P. vivax genome. The peptides identified in silico were chemically synthesized; circular dichroism studies indicated partial or high α-helical content. Antigenicity was evaluated using human sera samples from malaria-endemic areas of Colombia and Papua New Guinea. Eight of these fragments were selected and used to assess immunogenicity in BALB/c mice. ELISA assays indicated strong reactivity of serum samples from individuals residing in malaria-endemic regions and sera of immunized mice, with the α-helical coiled coil structures. In addition, ex vivo production of IFN-γ by murine mononuclear cells confirmed the immunogenicity of these structures and the presence of T-cell epitopes in the peptide sequences. Moreover, sera of mice immunized with four of the eight antigens recognized native proteins on blood-stage P. vivax parasites, and antigenic cross-reactivity with three of the peptides was observed when reacted with both the P. falciparum orthologous fragments and whole parasites. Results here point to the α-helical coiled coil peptides as possible P. vivax malaria vaccine candidates as were observed for P. falciparum. Fragments selected here warrant further study in humans and non-human primate models to assess their protective efficacy as single components or assembled as hybrid linear epitopes.

Cross-species and cross-platform gene expression studies with the Bioconductor-compliant R package 'annotationTools'.

Background: The variety of DNA microarray formats and datasets presently available offers an unprecedented opportunity to perform insightful comparisons of heterogeneous data. Cross-species studies, in particular, have the power of identifying conserved, functionally important molecular processes. Validation of discoveries can now often be performed in readily available public data which frequently requires cross-platform studies.Cross-platform and cross-species analyses require matching probes on different microarray formats. This can be achieved using the information in microarray annotations and additional molecular biology databases, such as orthology databases. Although annotations and other biological information are stored using modern database models ( e. g. relational), they are very often distributed and shared as tables in text files, i.e. flat file databases. This common flat database format thus provides a simple and robust solution to flexibly integrate various sources of information and a basis for the combined analysis of heterogeneous gene expression profiles.Results: We provide annotationTools, a Bioconductor-compliant R package to annotate microarray experiments and integrate heterogeneous gene expression profiles using annotation and other molecular biology information available as flat file databases. First, annotationTools contains a specialized set of functions for mining this widely used database format in a systematic manner. It thus offers a straightforward solution for annotating microarray experiments. Second, building on these basic functions and relying on the combination of information from several databases, it provides tools to easily perform cross-species analyses of gene expression data.Here, we present two example applications of annotationTools that are of direct relevance for the analysis of heterogeneous gene expression profiles, namely a cross-platform mapping of probes and a cross-species mapping of orthologous probes using different orthology databases. We also show how to perform an explorative comparison of disease-related transcriptional changes in human patients and in a genetic mouse model.Conclusion: The R package annotationTools provides a simple solution to handle microarray annotation and orthology tables, as well as other flat molecular biology databases. Thereby, it allows easy integration and analysis of heterogeneous microarray experiments across different technological platforms or species.

Comparison of splice sites in mammals and chicken

We have carried out an initial analysis of the dynamics of the recent evolution of the splice-sites sequences on a large collection of human, rodent (mouse and rat), and chicken introns. Our results indicate that the sequences of splice sites are largely homogeneous within tetrapoda. We have also found that orthologous splice signals between human and rodents and within rodents are more conserved than unrelated splice sites, but the additional conservation can be explained mostly by background intron conservation. In contrast, additional conservation over background is detectable in orthologous mammalian and chicken splice sites. Our results also indicate that the U2 and U12 intron classes seem to have evolved independently since the split of mammals and birds; we have not been able to find a convincing case of interconversion between these two classes in our collections of orthologous introns. Similarly, we have not found a single case of switching between AT-AC and GT-AG subtypes within U12 introns, suggesting that this event has been a rare occurrence in recent evolutionary times. Switching between GT-AG and the noncanonical GC-AG U2 subtypes, on the contrary, does not appear to be unusual; in particular, T to C mutations appear to be relatively well tolerated in GT-AG introns with very strong donor sites.

Association of NTRK3 and its interaction with NGF suggest an altered cross-regulation of the neurotrophin signaling pathway in eating disorders

Eating disorders (EDs) are complex psychiatric diseases that include anorexia nervosa and bulimia nervosa, and have higher than 50% heritability. Previous studies have found association of BDNF and NTRK2 to ED, while animal models suggest that other neurotrophin genes might also be involved in eating behavior. We have performed a family-based association study with 151 TagSNPs covering 10 neurotrophin signaling genes: NGFB, BDNF, NTRK1, NGFR/p75, NTF4/5, NTRK2, NTF3, NTRK3, CNTF and CNTFR in 371 ED trios of Spanish, French and German origin. Besides several nominal associations, we found a strong significant association after correcting for multiple testing (P = 1.04 × 10−4) between ED and rs7180942, located in the NTRK3 gene, which followed an overdominant model of inheritance. Interestingly, HapMap unrelated individuals carrying the rs7180942 risk genotypes for ED showed higher levels of expression of NTRK3 in lymphoblastoid cell lines. Furthermore, higher expression of the orthologous murine Ntrk3 gene was also detected in the hypothalamus of the anx/anx mouse model of anorexia. Finally, variants in NGFB gene appear to modify the risk conferred by the NTRK3 rs7180942 risk genotypes (P = 4.0 × 10−5) showing a synergistic epistatic interaction. The reported data, in addition to the previous reported findings for BDNF and NTRK2, point neurotrophin signaling genes as key regulators of eating behavior and their altered cross-regulation as susceptibility factors for EDs.

Transcription factor map alignment of promoter regions

We address the problem of comparing and characterizing the promoter regions of genes with similar expression patterns. This remains a challenging problem in sequence analysis, because often the promoter regions of co-expressed genes do not show discernible sequence conservation. In our approach, thus, we have not directly compared the nucleotide sequence of promoters. Instead, we have obtained predictions of transcription factor binding sites, annotated the predicted sites with the labels of the corresponding binding factors, and aligned the resulting sequences of labels—to which we refer here as transcription factor maps (TF-maps). To obtain the global pairwise alignment of two TF-maps, we have adapted an algorithm initially developed to align restriction enzyme maps. We have optimized the parameters of the algorithm in a small, but well-curated, collection of human–mouse orthologous gene pairs. Results in this dataset, as well as in an independent much larger dataset from the CISRED database, indicate that TF-map alignments are able to uncover conserved regulatory elements, which cannot be detected by the typical sequence alignments.

African signatures of recent positive selection in human FOXI1

Background: The human FOXI1 gene codes for a transcription factor involved in the physiology of the inner ear, testis, and kidney. Using three interspecies comparisons, it has been suggested that this may be a gene underhuman-specific selection. We sought to confirm this finding by using an extended set of orthologous sequences.Additionally, we explored for signals of natural selection within humans by sequencing the gene in 20 Europeans,20 East Asians and 20 Yorubas and by analysing SNP variation in a 2 Mb region centered on FOXI1 in 39worldwide human populations from the HGDP-CEPH diversity panel.Results: The genome sequences recently available from other primate and non-primate species showed that FOXI1divergence patterns are compatible with neutral evolution. Sequence-based neutrality tests were not significant inEuropeans, East Asians or Yorubas. However, the Long Range Haplotype (LRH) test, as well as the iHS and XP-Rsbstatistics revealed significantly extended tracks of homozygosity around FOXI1 in Africa, suggesting a recentepisode of positive selection acting on this gene. A functionally relevant SNP, as well as several SNPs either on theputatively selected core haplotypes or with significant iHS or XP-Rsb values, displayed allele frequencies stronglycorrelated with the absolute geographical latitude of the populations sampled.Conclusions: We present evidence for recent positive selection in the FOXI1 gene region in Africa. Climate mightbe related to this recent adaptive event in humans. Of the multiple functions of FOXI1, its role in kidney-mediatedwater-electrolyte homeostasis is the most obvious candidate for explaining a climate-related adaptation.

Comparative analysis of cancer genes in the human and chimpanzee genomes

Background: Cancer is a major medical problem in modern societies. However, the incidence of this disease in non-human primates is very low. To study whether genetic differences between human and chimpanzee could contribute to their distinct cancer susceptibility, we have examined in the chimpanzee genome the orthologous genes of a set of 333 human cancer genes. Results: This analysis has revealed that all examined human cancer genes are present in chimpanzee, contain intact open reading frames and show a high degree of conservation between both species. However, detailed analysis of this set of genes has shown some differences in genes of special relevance for human cancer. Thus, the chimpanzee gene encoding p53 contains a Pro residue at codon 72, while this codon is polymorphic in humans and can code for Arg or Pro, generating isoforms with different ability to induce apoptosis or interact with p73. Moreover, sequencing of the BRCA1 gene has shown an 8 Kb deletion in the chimpanzee sequence that prematurely truncates the co-regulated NBR2 gene. Conclusion: These data suggest that small differences in cancer genes, as those found in tumor suppressor genes, might influence the differences in cancer susceptibility between human and chimpanzee. Nevertheless, further analysis will be required to determine the exact contribution of the genetic changes identified in this study to the different cancer incidence in non-human primates.