Molecular Modeling Approaches for Determining Gene Function: application to a Putative Poly-A Binding Protein from Leishmania amazonensis (LaPABP)

The great expansion in the number of genome sequencing projects has revealed the importance of computational methods to speed up the characterization of unknown genes. These studies have been improved by the use of three dimensional information from the predicted proteins generated by molecular modeling techniques. In this work, we disclose the structure-function relationship of a gene product from Leishmania amazonensis by applying molecular modeling and bioinformatics techniques. The analyzed sequence encodes a 159 aminoacids polypeptide (estimated 18 kDa) and was denoted LaPABP for its high homology with poly-A binding proteins from trypanosomatids. The domain structure, clustering analysis and a three dimensional model of LaPABP, basically obtained by homology modeling on the structure of the human poly-A binding protein, are described. Based on the analysis of the electrostatic potential mapped on the model's surface and conservation of intramolecular contacts responsible for folding stabilization we hypothesize that this protein may have less avidity to RNA than it's L. major counterpart but still account for a significant functional activity in the parasite. The model obtained will help in the design of mutagenesis experiments aimed to elucidate the mechanism of gene expression in trypanosomatids and serve as a starting point for its exploration as a potential source of targets for a rational chemotherapy.

Classification of human astrocytic gliomas on the basis of gene expression: a correlated group of genes with angiogenic activity emerges as a strong predictor of subtypes.

The development of targeted treatment strategies adapted to individual patients requires identification of the different tumor classes according to their biology and prognosis. We focus here on the molecular aspects underlying these differences, in terms of sets of genes that control pathogenesis of the different subtypes of astrocytic glioma. By performing cDNA-array analysis of 53 patient biopsies, comprising low-grade astrocytoma, secondary glioblastoma (respective recurrent high-grade tumors), and newly diagnosed primary glioblastoma, we demonstrate that human gliomas can be differentiated according to their gene expression. We found that low-grade astrocytoma have the most specific and similar expression profiles, whereas primary glioblastoma exhibit much larger variation between tumors. Secondary glioblastoma display features of both other groups. We identified several sets of genes with relatively highly correlated expression within groups that: (a). can be associated with specific biological functions; and (b). effectively differentiate tumor class. One prominent gene cluster discriminating primary versus nonprimary glioblastoma comprises mostly genes involved in angiogenesis, including VEGF fms-related tyrosine kinase 1 but also IGFBP2, that has not yet been directly linked to angiogenesis. In situ hybridization demonstrating coexpression of IGFBP2 and VEGF in pseudopalisading cells surrounding tumor necrosis provided further evidence for a possible involvement of IGFBP2 in angiogenesis. The separating groups of genes were found by the unsupervised coupled two-way clustering method, and their classification power was validated by a supervised construction of a nearly perfect glioma classifier.

Drug resistance and genotypes of strains of Mycobacterium tuberculosis isolated from human immunodeficiency virus-infected and non-infected tuberculosis patients in Bauru, São Paulo, Brazil

Little is known about transmission and drug resistance of tuberculosis (TB) in Bauru, State of São Paulo. The objective of this study was to evaluate risk factors for transmission of Mycobacterium tuberculosis strains in this area. Strains were collected from patients attended at ambulatory services in the region and susceptibility towards the main first line antibiotics was determined and fingerprinting performed. A total of 57 strains were submitted to susceptibility testing: 23 (42.6%) were resistant to at least one drug while 3 (13%) were resistant against both rifampicin and isoniazide. Resistant strains had been isolated from patients that had not (n = 13) or had (n = 9) previously been submitted to anti-TB treatment, demonstrating a preoccupying high level of primary resistance in the context of the study. All strains were submitted to IS6110 restriction fragment length polymorphism (IS6110-RFLP) and double repetitive element PCR (DRE-PCR). Using IS6110-RFLP, 26.3% of the strains were clustered and one cluster of 3 patients included 2 HIV-infected individuals that had been hospitalized together during 16 days; clustering of strains of patients from the hospital was however not higher than that of patients attended at health posts. According to DRE-PCR, 55.3% belonged to a cluster, confirming the larger discriminatory power of IS6110-RFLP when compared to DRE-PCR, that should therefore be used as a screening procedure only. No clinical, epidemiological or microbiological characteristics were associated with clustering so risk factors for transmission of TB could not be defined in the present study.

Genetic analysis of phenoxyalkanoic acid degradation in Sphingomonas herbicidovorans MH.

Phenoxyalkanoic acid degradation is well studied in Beta- and Gammaproteobacteria, but the genetic background has not been elucidated so far in Alphaproteobacteria. We report the isolation of several genes involved in dichlor- and mecoprop degradation from the alphaproteobacterium Sphingomonas herbicidovorans MH and propose that the degradation proceeds analogously to that previously reported for 2,4-dichlorophenoxyacetic acid (2,4-D). Two genes for alpha-ketoglutarate-dependent dioxygenases, sdpA(MH) and rdpA(MH), were found, both of which were adjacent to sequences with potential insertion elements. Furthermore, a gene for a dichlorophenol hydroxylase (tfdB), a putative regulatory gene (cadR), two genes for dichlorocatechol 1,2-dioxygenases (dccA(I/II)), two for dienelactone hydrolases (dccD(I/II)), part of a gene for maleylacetate reductase (dccE), and one gene for a potential phenoxyalkanoic acid permease were isolated. In contrast to other 2,4-D degraders, the sdp, rdp, and dcc genes were scattered over the genome and their expression was not tightly regulated. No coherent pattern was derived on the possible origin of the sdp, rdp, and dcc pathway genes. rdpA(MH) was 99% identical to rdpA(MC1), an (R)-dichlorprop/alpha-ketoglutarate dioxygenase from Delftia acidovorans MC1, which is evidence for a recent gene exchange between Alpha- and Betaproteobacteria. Conversely, DccA(I) and DccA(II) did not group within the known chlorocatechol 1,2-dioxygenases, but formed a separate branch in clustering analysis. This suggests a different reservoir and reduced transfer for the genes of the modified ortho-cleavage pathway in Alphaproteobacteria compared with the ones in Beta- and Gammaproteobacteria.

IS6110 restriction fragment length polymorphism of Mycobacterium tuberculosis isolated from patients with pulmonary tuberculosis in Campinas, Brazil: evidence of intercontinental distribution of strains

Tuberculosis (TB) is a major concern in developing countries. In Brazil, few genotyping studies have been conducted to verify the number of IS6110 copies present in local prevalent strains of Mycobacterium tuberculosis, the distribution and clustering of strains. IS6110 DNA fingerprinting was performed on a sample of M. tuberculosis isolates from patients with AFB smear-positive pulmonary TB, at a hospital in Brazil. The IS6110 profiles were analyzed and compared to a M. tuberculosis database of the Houston Tuberculosis Initiative, Houston, US. Seventy-six fingerprints were obtained from 98 patients. All M. tuberculosis strains had an IS6110 copy number between 5-21 allowing for differentiation of the isolates. Human immunodeficiency virus infection was confirmed in nearly half the patients of whom data was available. Fifty-eight strains had unique patterns, while 17 strains were grouped in 7 clusters (2 to 6 strains). When compared to the HTI database, 6 strains matched isolates from El Paso, Ciudad de Juarez, Houston, and New York. Recently acquired infections were documented in 19% of cases. The community transmission of infection is intense, since some clustered strains were recovered during the four-year study period. The intercontinental dissemination of M. tuberculosis strains is suspected by demonstration of identical fingerprints in a distant country.

A comprehensive characterization of genome-wide copy number aberrations in colorectal cancer reveals novel oncogenes and patterns of alterations.

To develop a comprehensive overview of copy number aberrations (CNAs) in stage-II/III colorectal cancer (CRC), we characterized 302 tumors from the PETACC-3 clinical trial. Microsatellite-stable (MSS) samples (n = 269) had 66 minimal common CNA regions, with frequent gains on 20 q (72.5%), 7 (41.8%), 8 q (33.1%) and 13 q (51.0%) and losses on 18 (58.6%), 4 q (26%) and 21 q (21.6%). MSS tumors have significantly more CNAs than microsatellite-instable (MSI) tumors: within the MSI tumors a novel deletion of the tumor suppressor WWOX at 16 q23.1 was identified (p<0.01). Focal aberrations identified by the GISTIC method confirmed amplifications of oncogenes including EGFR, ERBB2, CCND1, MET, and MYC, and deletions of tumor suppressors including TP53, APC, and SMAD4, and gene expression was highly concordant with copy number aberration for these genes. Novel amplicons included putative oncogenes such as WNK1 and HNF4A, which also showed high concordance between copy number and expression. Survival analysis associated a specific patient segment featured by chromosome 20 q gains to an improved overall survival, which might be due to higher expression of genes such as EEF1B2 and PTK6. The CNA clustering also grouped tumors characterized by a poor prognosis BRAF-mutant-like signature derived from mRNA data from this cohort. We further revealed non-random correlation between CNAs among unlinked loci, including positive correlation between 20 q gain and 8 q gain, and 20 q gain and chromosome 18 loss, consistent with co-selection of these CNAs. These results reinforce the non-random nature of somatic CNAs in stage-II/III CRC and highlight loci and genes that may play an important role in driving the development and outcome of this disease.

Ku-deficient yeast strains exhibit alternative states of silencing competence.

In Saccharomyces cerevisiae, efficient silencer function requires telomere proximity, i.e. compartments of the nucleoplasm enriched in silencing factors. Accordingly, silencers located far from telomeres function inefficiently. We show here that cells lacking yKu balance between two mitotically stable states of silencing competence. In one, a partial delocalization of telomeres and silencing factors throughout the nucleoplasm correlates with enhanced silencing at a non-telomeric locus, while in the other, telomeres retain their focal pattern of distribution and there is no repression at the non-telomeric locus, as observed in wild-type cells. The two states also differ in their level of residual telomeric silencing. These findings indicate the existence of a yKu-independent pathway of telomere clustering and Sir localization. Interestingly, this pathway appears to be under epigenetic control.

Spatial analysis of the distribution of leprosy in the State of Ceará, Northeast Brazil

The aim of this study was to describe spatial patterns of the distribution of leprosy and to investigate spatial clustering of incidence rates in the state of Ceará, Northeast Brazil. The average incidence rate of leprosy for the period of 1991 to 1999 was calculated for each municipality of Ceará. Maps were used to describe the spatial distribution of the disease, and spatial statistics were applied to explore large- and small-scale variations of incidence rates. Three regions were identified in which the incidence of leprosy was particularly high. A spatial gradient in the incidence rates was identified, with a tendency of high rates to be concentrated on the North-South axis in the middle region of the state. Moran's I statistic indicated that a significant spatial autocorrelation also existed. The spatial distribution of leprosy in Ceará is heterogeneous. The reasons for spatial clustering of disease rates are not known, but might be related to an heterogeneous distribution of other factors such as crowding, social inequality, and environmental characteristics which by themselves determine the transmission of Mycobacterium leprae.

Soil fungal communities of grasslands are environmentally structured at a regional scale in the Alps.

Studying patterns of species distributions along elevation gradients is frequently used to identify the primary factors that determine the distribution, diversity and assembly of species. However, despite their crucial role in ecosystem functioning, our understanding of the distribution of below-ground fungi is still limited, calling for more comprehensive studies of fungal biogeography along environmental gradients at various scales (from regional to global). Here, we investigated the richness of taxa of soil fungi and their phylogenetic diversity across a wide range of grassland types along a 2800 m elevation gradient at a large number of sites (213), stratified across a region of the Western Swiss Alps (700 km(2)). We used 454 pyrosequencing to obtain fungal sequences that were clustered into operational taxonomic units (OTUs). The OTU diversity-area relationship revealed uneven distribution of fungal taxa across the study area (i.e. not all taxa are everywhere) and fine-scale spatial clustering. Fungal richness and phylogenetic diversity were found to be higher in lower temperatures and higher moisture conditions. Climatic and soil characteristics as well as plant community composition were related to OTU alpha, beta and phylogenetic diversity, with distinct fungal lineages suggesting distinct ecological tolerances. Soil fungi, thus, show lineage-specific biogeographic patterns, even at a regional scale, and follow environmental determinism, mediated by interactions with plants.

IS6110 restriction fragment length polymorphism typing of Mycobacterium tuberculosis isolates from East Azerbaijan Province of Iran

To investigate the genetic variation among Mycobacterium tuberculosis isolates in the East Azerbaijan Province of Iran and to evaluate the level of and risk factors for recent transmission of tuberculosis (TB), we performed IS6110-based restriction fragment length polymorphism analysis of strains, isolated from 105 patients during the period of September 2002 to March 2003 in TB centers and university hospitals of the province. Among 105 isolates, 81 different IS6110 patterns were found, of which 70 were observed only once and 11 were shared by two to eight isolates. Ninety-six isolates (91.4%) were found to have more than five copies of IS6110 and together with high patterns polymorphism, shows that IS6110-RFLP typing could be useful for studying the epidemiology of TB in Azerbaijan. The minimum estimated rate of recent transmission was 23%, suggesting that the degree of recent transmission in East Azerbaijan Province is relatively low. Clustering was not associated with age, sex or site of infection of TB but drug-resistant isolates were less likely to be clustered than sensitive isolates (p < 0.05).

The spatial distribution of Schistosoma mansoni infection before and after chemotherapy in the Jequitinhonha Valley in Brazil

Schistosomiasis prevalence and egg counts remained low one year after chemotherapy in most households in a hyperendemic rural area in northern Minas Gerais but several distinct spatial patterns could be observed in relation to IgE levels and to a lesser extent to exposure risk (TBM) and type of water supply. An inverse relationship between pre-treatment household prevalence and egg counts on the one hand and post-treatment IgE levels on the other were noted in two of the five communities. Low exposure risk was associated with the low pre-treatment infection rates in the central village but did not contribute to the decline of infection rates after chemotherapy in the study area, as indicated by the significant increase in water contact during the posttreatment period (p < 0.0001). Distance between households and the streams and socioeconomic factors were also unimportant in predicting the spatial distribution of infection. These results are consistent with the production and antiparasitic effect of high levels of IgE in Schistosoma mansoni infection.

Man biting rate seasonal variation of malaria vectors in Roraima, Brazil

Malaria control has been directed towards regional actions where more detailed knowledge of local determinants of transmission is of primary importance. This is a short report on range distribution and biting indices for Anopheles darlingi and An. albitarsis during the dry and rainy season that follows river level variation in a savanna/alluvial forest malaria system area in the Northern Amazon Basin. Distribution range and adult biting indices were at their highest during the rainy season for both An. darlingi and An. albitarsis. During the rainy season the neighboring alluvial forest was extensively flooded. This coincided with highest rates in malaria transmission with case clustering near the river. As the river receded, anopheline distribution range and density decreased. This decrease in distribution and density corresponded to a malaria decrease in the near area. An exponential regression function was derived to permit estimations of An. darlingi distribution over specified distances. Anopheline spatio-temporal variations lead to uneven malaria case distribution and are of important implications for control strategies.

Typing of Mycobacterium tuberculosis strains isolated in Community Health Centers of Rio de Janeiro City, Brazil

Fingerprinting of Mycobacterium tuberculosis strains from tuberculosis (TB) patients attended in Community Health Centers (CHCs) of Rio de Janeiro was performed to verify possible risk factors for TB transmission. A prospective community-based study was performed during the period of July 1996 to December 1996 by collecting sputum samples of 489 patients in 11 different CHCs in four different planning areas (APs) of the city. Bacteriological, clinical, and epidemiological information was collected and M. tuberculosis genotypes defined after restriction fragment length polymorphism (IS6110-RFLP) and double repetitive element (DRE) fingerprinting of RFLP-clustered cases. Risk factors for TB transmission were looked for using three levels of cluster stringency. Among 349 (71%) positive cultures obtained, IS6110-RFLP typing could be performed on strains from 153 different patients. When using identity of RFLP patterns as cluster definition, 49 (32%) of the strains belonged to a cluster and none of the clinical or epidemiologic characteristics was associated with higher clustering levels. However, higher clustering level was observed in the AP including the central region of the city when compared to others. This strongly suggests that more recent transmission occurs in that area and this may be related with higher incidence of TB and HIV in this region.

A preliminary analysis of the population genetics and molecular phylogenetics of Onchocerca volvulus (Nematoda: Filarioidea) using nuclear ribosomal second internal transcribed spacer sequences

Nuclear internal transcribed spacer 2 (ITS2) rDNA sequences were used for a molecular phylogenetics analysis of five Onchocerca species. The sister species of the human parasite O. volvulus was found to be the cattle parasite O. ochengi and not O. gibsoni, contrary to chromosomal evidence. The genetic differentiation of two African populations (representing the two African strains) and a Brazilian population of O. volvulus was also studied. Phylogenetic and network reconstruction did not show any clustering of ITS2 alleles on geographic or strain grounds. Furthermore, population genetics tests showed no indication of population differentiation but suggested gene flow among the three populations.

Artificial neural network study of whole-cell bacterial bioreporter response determined using fluorescence flow cytometry.

Genetically engineered bioreporters are an excellent complement to traditional methods of chemical analysis. The application of fluorescence flow cytometry to detection of bioreporter response enables rapid and efficient characterization of bacterial bioreporter population response on a single-cell basis. In the present study, intrapopulation response variability was used to obtain higher analytical sensitivity and precision. We have analyzed flow cytometric data for an arsenic-sensitive bacterial bioreporter using an artificial neural network-based adaptive clustering approach (a single-layer perceptron model). Results for this approach are far superior to other methods that we have applied to this fluorescent bioreporter (e.g., the arsenic detection limit is 0.01 microM, substantially lower than for other detection methods/algorithms). The approach is highly efficient computationally and can be implemented on a real-time basis, thus having potential for future development of high-throughput screening applications.