Alfonso Méndez-Tenorio

Phylogenetic and Structural Analysis of Translationally Controlled Tumor Proteins

The translationally controlled tumor protein (TCTP) is conserved in all eukaryotes studied thus far. Recent evidence points to an important role for TCTP in the induction of cell proliferation in animals through an interaction with G proteins. TCTP may also constitute an intercellular secreted signal that modulates the immune response in the vertebrates. Because of its sequence conservation and ubiquity, the analysis of its amino acid sequence divergence between different taxa may provide insight into the structural constraints on the evolution of this protein. In the present study, we analyzed the phylogeny of TCTP sequences from a wide range of organisms and found that, with some exceptions, the groupings formed were consistent with the evolutionary history. Indeed, at the level of lower-order taxa, the groupings are in agreement with their established phylogeny, thus indicating that the substitution rates of the TCTP residues varied evenly between members of the same clade. Predicted three-dimensional structures of representative TCTPs, based on the reported 3D structure of Schizosaccharomyces pombe, indicated that these proteins are highly conserved among diverse taxonomic groups. However, analysis of the primary structure indicated subtle differences in the domain-forming pocket that potentially interacts with G proteins, particularly among Diplomonadidae, Apicomplexa, and other parasites of vertebrates. These differences support the notion that these specific TCTPs could block the normal immune response by acting as dominant negative mutants. Structural differences were also observed in a reported sequence of TCTP from Plasmodium knowlesi, in which the presence of an extra α-helix could also interfere in the interaction with G proteins.

Geographical variation in human papillomavirus prevalence in Mexican women with normal cytology

OBJECTIVE To determine the prevalence of human papillomavirus (HPV) infection and genotype distribution in Mexican women with similar lifestyles from two geographical regions who receive medical care from the Mexican Navy Health System, and to identify the associated sociodemographic and reproductive characteristics. METHODS Cervical swabs from 671 women, beneficiaries of the Mexican Navy Health System, from two distinct southern coast regions of Mexico, were analyzed. Data were obtained regarding sociodemographic variables and sexual and reproductive history. For HPV detection and typing, PCR with general primers and direct sequencing were performed on extracted DNA. Association with clinical variables was evaluated. RESULTS Most patients had a normal cytology or low-grade intraepithelial neoplasia. A high prevalence of HPV was found (43.6%), with a significant difference between the two regions studied from the southwest Pacific coast of Mexico (37.6% in Acapulco, Guerrero vs. 49.7% in Lázaro Cárdenas, Michoacán). Some differences were also found associated to HPV type distribution, particularly related to genotypes 18, 58, and 53. Factors influencing these differences could not be identified with the analysis of typical risk factors linked to the acquisition of an HPV infection. CONCLUSIONS Regional differences in HPV prevalence and distribution show an apparent geographic boundary between the studied populations that deserves further analysis, taking into account other factors such as those related to the sexual partners.

Whole Genome Sequence and Phylogenetic Analysis Show Helicobacter pylori Strains from Latin America Have Followed a Unique Evolution Pathway

Helicobacter pylori (HP) genetics may determine its clinical outcomes. Despite high prevalence of HP infection in Latin America (LA), there have been no phylogenetic studies in the region. We aimed to understand the structure of HP populations in LA mestizo individuals, where gastric cancer incidence remains high. The genome of 107 HP strains from Mexico, Nicaragua and Colombia were analyzed with 59 publicly available worldwide genomes. To study bacterial relationship on whole genome level we propose a virtual hybridization technique using thousands of high-entropy 13 bp DNA probes to generate fingerprints. Phylogenetic virtual genome fingerprint (VGF) was compared with Multi Locus Sequence Analysis (MLST) and with phylogenetic analyses of cagPAI virulence island sequences. With MLST some Nicaraguan and Mexican strains clustered close to Africa isolates, whereas European isolates were spread without clustering and intermingled with LA isolates. VGF analysis resulted in increased resolution of populations, separating European from LA strains. Furthermore, clusters with exclusively Colombian, Mexican, or Nicaraguan strains were observed, where the Colombian cluster separated from Europe, Asia, and Africa, while Nicaraguan and Mexican clades grouped close to Africa. In addition, a mixed large LA cluster including Mexican, Colombian, Nicaraguan, Peruvian, and Salvadorian strains was observed; all LA clusters separated from the Amerind clade. With cagPAI sequence analyses LA clades clearly separated from Europe, Asia and Amerind, and Colombian strains formed a single cluster. A NeighborNet analyses suggested frequent and recent recombination events particularly among LA strains. Results suggests that in the new world, H. pylori has evolved to fit mestizo LA populations, already 500 years after the Spanish colonization. This co-adaption may account for regional variability in gastric cancer risk.

Low density DNA microarray for detection of most frequent TP53 missense point mutations

BackgroundWe have developed an oligonucleotide microarray (genosensor) utilizing a double tandem hybridization technique to search for 9 point mutations located in the most frequently altered codons of the TP53 gene. Isolated and multiplexed PCR products, 108 and 92 bp long, from exons 7 and 8, respectively, were obtained from 24 different samples. Single-stranded target DNA was then prepared from isolated or multiplexed PCR products, through cyclic DNA synthesis. Independent ssDNAs were annealed with the corresponding pairs of labeled stacking oligonucleotides to create partially duplex DNA having a 7-nt gap, which contains the sequence that will be interrogated by the capture probes forming double tandem hybridization. In the case of multiplexed ssPCR products, only two stacking oligonucleotides were added per target, therefore the gap for the PCR products having two consecutive codons to be interrogated in exon 7 was 12 nt long, so only single tandem hybridization was produced with these respective probes.Results18 codon substitutions were found by DNA sequencing. In 13 of them a perfect correlation with the pattern of hybridization was seen (In 5 no signal was seen with the wt probe while a new signal was seen with the appropriate mutant probe, and in 8 more, as expected, no signal was seen with any probe due to the absence of the corresponding probe in the array). In 3 other cases a mutation was falsely suggested by the combination of the absence of the wild type signal along with a false signal in the other probe. In the other 2 cases the presence of the mutation was not detected due to the production of a false hybridization signal with the wild type probe. In both cases (false mutation or no mutation detected) relatively stable mismatched target/probe duplexes should be formed. These problems could be avoided by the addition of probes to improve the performance of the array.ConclusionOur results demonstrate that a simple TP53 microarray employing short (7-mer) probes, used in combination with single or double tandem hybridization approach and a simple or multiplex target preparation method, can identify common TP53 missense mutations from a variety of DNA sources with good specificity.

Detection of mutations in RET proto-oncogene codon 634 through double tandem hybridization

We developed a procedure to detect the 7 point mutations at Cys634 of the proto-oncogene RET, which is responsible for medullary thyroid carcinoma (MTC). Genomic DNA was prepared from blood samples obtained from normal and MTC-affected individuals belonging to a family with a history of the disease. The RET genotype for each individual was first established by performing restriction and sequencing analyses. Single-stranded target DNA was prepared by asymmetric polymerase chain reaction (PCR) amplification of a 93-bp fragment containing Cys634. The target was annealed with pairs of prelabeled stacking oligonucleotides designed to create appropriate 7-nucleotide gaps, which served as the sites of subsequent hybridization with glass-immobilized 7-mer probes. The target-stacking oligonucleotide duplexes were hybridized with DNA chips containing a set of eight 7-mer probes designed to detect the wild-type sequence and the seven point mutations described. We tested two sets of immobilized probes containing internal or 5′-terminal codon-634 single-base variations. Both groups of probes were able to discriminatively identify the mutations. The hybridization patterns indicated that the disease in this family was due to the C634Y mutation, in accord with the original sequence analysis. The hybridization-based mutation assignment was additionally supported by determination of the control homozygous and heterozygous hybridization patterns produced with synthetic targets having the normal or codon 634 mutant sequences. The effects of mismatch type and nearest-neighbor sequences on the occurrence of false-positive (mismatched) hybridizations are discussed.

Bioinformatics prediction of siRNAs as potential antiviral agents against dengue viruses.

Dengue virus (DENV 1-4) represents the major emerging arthropod-borne viral infection in the world. Currently, there is neither an available vaccine nor a specific treatment. Hence, there is a need of antiviral drugs for these viral infections; we describe the prediction of short interfering RNA (siRNA) as potential therapeutic agents against the four DENV serotypes. Our strategy was to carry out a series of multiple alignments using ClustalX program to find conserved sequences among the four DENV serotype genomes to obtain a consensus sequence for siRNAs design. A highly conserved sequence among the four DENV serotypes, located in the encoding sequence for NS4B and NS5 proteins was found. A total of 2,893 complete DENV genomes were downloaded from the NCBI, and after a depuration procedure to identify identical sequences, 220 complete DENV genomes were left. They were edited to select the NS4B and NS5 sequences, which were aligned to obtain a consensus sequence. Three different servers were used for siRNA design, and the resulting siRNAs were aligned to identify the most prevalent sequences. Three siRNAs were chosen, one targeted the genome region that codifies for NS4B protein and the other two; the region for NS5 protein. Predicted secondary structure for DENV genomes was used to demonstrate that the siRNAs were able to target the viral genome forming double stranded structures, necessary to activate the RNA silencing machinery.

LifePrint: a novel k-tuple distance method for construction of phylogenetic trees

Purpose Here we describe LifePrint, a sequence alignment-independent k-tuple distance method to estimate relatedness between complete genomes. Methods We designed a representative sample of all possible DNA tuples of length 9 (9-tuples). The final sample comprises 1878 tuples (called the LifePrint set of 9-tuples; LPS9) that are distinct from each other by at least two internal and noncontiguous nucleotide differences. For validation of our k-tuple distance method, we analyzed several real and simulated viroid genomes. Using different distance metrics, we scrutinized diverse viroid genomes to estimate the k-tuple distances between these genomic sequences. Then we used the estimated genomic k-tuple distances to construct phylogenetic trees using the neighbor-joining algorithm. A comparison of the accuracy of LPS9 and the previously reported 5-tuple method was made using symmetric differences between the trees estimated from each method and a simulated “true” phylogenetic tree. Results The identified optimal search scheme for LPS9 allows only up to two nucleotide differences between each 9-tuple and the scrutinized genome. Similarity search results of simulated viroid genomes indicate that, in most cases, LPS9 is able to detect single-base substitutions between genomes efficiently. Analysis of simulated genomic variants with a high proportion of base substitutions indicates that LPS9 is able to discern relationships between genomic variants with up to 40% of nucleotide substitution. Conclusion Our LPS9 method generates more accurate phylogenetic reconstructions than the previously proposed 5-tuples strategy. LPS9-reconstructed trees show higher bootstrap proportion values than distance trees derived from the 5-tuple method.

Bacterial classification using genomic fingerprints obtained by virtual hybridization.

In silico genomic fingerprints were produced by virtual hybridization of 191 fully sequenced bacterial genomes using a set of 15,264 13-mer probes specially designed to produce universal whole genome fingerprints. A novel approach for constructing phylogenetic trees, based on comparative analysis of genomic fingerprints, was developed. The resultant bacterial phylogenetic tree had strong similarities to those produced from the alignment of conserved sequences. Notably, the trees derived from the alignment of other conserved COG genes divided the Bacillus and Corynebacterium genera into the same subgroups produced by the novel bacterial tree. A number of discrepancies between both techniques were observed for the grouping of some Lactobacillus species. However, a detailed analysis of the alignment of these genomes using other bioinformatics tools revealed that the grouping of these organisms in the novel tree was more satisfactory than the groupings from previous classifications, which used only a few conserved genes. All these data suggest that the bacterial taxonomy produced by genomic fingerprints is satisfactory, but sometimes different from classical taxonomies. Discrepancies probably arise because the fingerprinting technique analyzes genomic sequences and reveals more information than previously used approaches.

Overview of recurrent chromosomal losses in retinoblastoma detected by low coverage next generation sequencing.

Genes are frequently lost or gained in malignant tumors and the analysis of these changes can be informative about the underlying tumor biology. Retinoblastoma is a pediatric intraocular malignancy, and since deletions in chromosome 13 have been described in this tumor, we performed genome wide sequencing with the Illumina platform to test whether recurrent losses could be detected in low coverage data from DNA pools of Rb cases. An in silico reference profile for each pool was created from the human genome sequence GRCh37p5; a chromosome integrity score and a graphics 40 Kb window analysis approach, allowed us to identify with high resolution previously reported non random recurrent losses in all chromosomes of these tumors. We also found a pattern of gains and losses associated to clear and dark cytogenetic bands respectively. We further analyze a pool of medulloblastoma and found a more stable genomic profile and previously reported losses in this tumor. This approach facilitates identification of recurrent deletions from many patients that may be biological relevant for tumor development.

Conformational changes associated with L16P and T118M mutations in the membrane-embedded PMP22 protein, consequential in CMT-1A

Peripheral myelin protein 22 (PMP22) resides in the plasma membrane and is required for myelin formation in the peripheral nervous system. Excess PMP22 mutants accumulate in the endoplasmic reticulum (ER) resulting in the inherited neuropathies of Charcot–Marie–Tooth disease. However, there was no evidence of the structure of PMP22 or how mutations affect its folding. Therefore, in this study, we combined bioinformatics and homology modeling approaches to obtain three-dimensional native and mutated PMP22 models and its anchoring to a POPC membrane, submitted to .5-μs MD simulations, to determine how the L16P and T118M mutations affect the conformational behavior of PMP22. In addition, we investigated the ability of the native and mutated species to accumulate in the ER, via interaction with RER1, by combining protein–protein docking and MD simulations, taking the conformations that were most representative of the native and mutated PMP22 systems and RER1 conformations. Principal component analysis over MD simulations revealed that L16P and T118M mutations resulted in increased structural instability compared to the native form, which is consistent with previous experimental findings of increased structural fluctuations along a loop connecting transmembrane α-helix1 and α-helix2. Docking and MD simulations coupled with the MMGBSA approach allowed the identification that the binding interface for the PMP22-RER1 complex takes place through transmembrane α-helix1 and α-helix2, with higher effective binding free energy values between the mutated PMP22 systems and RER1 than for the native PMP22, mainly through van der Waals interactions.