Lossless filter for multiple repeats with bounded edit distance

Background: Identifying local similarity between two or more sequences, or identifying repeats occurring at least twice in a sequence, is an essential part in the analysis of biological sequences and of their phylogenetic relationship. Finding such fragments while allowing for a certain number of insertions, deletions, and substitutions, is however known to be a computationally expensive task, and consequently exact methods can usually not be applied in practice. Results: The filter TUIUIU that we introduce in this paper provides a possible solution to this problem. It can be used as a preprocessing step to any multiple alignment or repeats inference method, eliminating a possibly large fraction of the input that is guaranteed not to contain any approximate repeat. It consists in the verification of several strong necessary conditions that can be checked in a fast way. We implemented three versions of the filter. The first is simply a straightforward extension to the case of multiple sequences of an application of conditions already existing in the literature. The second uses a stronger condition which, as our results show, enable to filter sensibly more with negligible (if any) additional time. The third version uses an additional condition and pushes the sensibility of the filter even further with a non negligible additional time in many circumstances; our experiments show that it is particularly useful with large error rates. The latter version was applied as a preprocessing of a multiple alignment tool, obtaining an overall time (filter plus alignment) on average 63 and at best 530 times smaller than before (direct alignment), with in most cases a better quality alignment. Conclusion: To the best of our knowledge, TUIUIU is the first filter designed for multiple repeats and for dealing with error rates greater than 10% of the repeats length.

Survival time of random walk in random environment among soft obstacles

We consider a Random Walk in Random Environment (RWRE) moving in an i.i.d. random field of obstacles. When the particle hits an obstacle, it disappears with a positive probability. We obtain quenched and annealed bounds on the tails of the survival time in the general d-dimensional case. We then consider a simplified one-dimensional model (where transition probabilities and obstacles are independent and the RWRE only moves to neighbour sites), and obtain finer results for the tail of the survival time. In addition, we study also the ""mixed"" probability measures (quenched with respect to the obstacles and annealed with respect to the transition probabilities and vice-versa) and give results for tails of the survival time with respect to these probability measures. Further, we apply the same methods to obtain bounds for the tails of hitting times of Branching Random Walks in Random Environment (BRWRE).

KeaA, a Dictyostelium kelch-domain protein that regulates the response to stress and development

Background: The protein kinase YakA is responsible for the growth arrest and induction of developmental processes that occur upon starvation of Dictyostelium cells. yakA-cells are aggregation deficient, have a faster cell cycle and are hypersensitive to oxidative and nitrosoative stress. With the aim of isolating members of the YakA pathway, suppressors of the death induced by nitrosoative stress in the yakA-cells were identified. One of the suppressor mutations occurred in keaA, a gene identical to DG1106 and similar to Keap1 from mice and the Kelch protein from Drosophila, among others that contain Kelch domains. Results: A mutation in keaA suppresses the hypersensitivity to oxidative and nitrosoative stresses but not the faster growth phenotype of yakA-cells. The growth profile of keaA deficient cells indicates that this gene is necessary for growth. keaA deficient cells are more resistant to nitrosoative and oxidative stress and keaA is necessary for the production and detection of cAMP. A morphological analysis of keaA deficient cells during multicellular development indicated that, although the mutant is not absolutely deficient in aggregation, cells do not efficiently participate in the process. Gene expression analysis using cDNA microarrays of wild-type and keaA deficient cells indicated a role for KeaA in the regulation of the cell cycle and pre-starvation responses. Conclusions: KeaA is required for cAMP signaling following stress. Our studies indicate a role for kelch proteins in the signaling that regulates the cell cycle and development in response to changes in the environmental conditions.

Disparity between Multilocus Enzyme Electrophoresis, Microsatellite Markers and Pulsed-Field Gel Electrophoresis in epidemiological tracking of Candida albicans

Various molecular systems are available for epidemiological, genetic, evolutionary, taxonomic and systematic studies of innumerable fungal infections, especially those caused by the opportunistic pathogen C. albicans. A total of 75 independent oral isolates were selected in order to compare Multilocus Enzyme Electrophoresis (MLEE), Electrophoretic Karyotyping (EK) and Microsatellite Markers (Simple Sequence Repeats - SSRs), in their abilities to differentiate and group C. albicans isolates (discriminatory power), and also, to evaluate the concordance and similarity of the groups of strains determined by cluster analysis for each fingerprinting method. Isoenzyme typing was performed using eleven enzyme systems: Adh, Sdh, M1p, Mdh, Idh, Gdh, G6pdh, Asd, Cat, Po, and Lap (data previously published). The EK method consisted of chromosomal DNA separation by pulsed-field gel electrophoresis using a CHEF system. The microsatellite markers were investigated by PCR using three polymorphic loci: EF3, CDC3, and HIS3. Dendrograms were generated by the SAHN method and UPGMA algorithm based on similarity matrices (S(SM)). The discriminatory power of the three methods was over 95%, however a paired analysis among them showed a parity of 19.7-22.4% in the identification of strains. Weak correlation was also observed among the genetic similarity matrices (S(SM)(MLEE) x S(SM)(EK) x S(SM)(SSRs)). Clustering analyses showed a mean of 9 +/- 12.4 isolates per cluster (3.8 +/- 8 isolates/taxon) for MLEE, 6.2 +/- 4.9 isolates per cluster (4 +/- 4.5 isolates/taxon) for SSRs, and 4.1 +/- 2.3 isolates per cluster (2.6 +/- 2.3 isolates/taxon) for EK. A total of 45 (13%), 39(11.2%), 5 (1.4%) and 3 (0.9%) clusters pairs from 347 showed similarity (Si) of 0.1-10%, 10.1-20%, 20.1-30% and 30.1-40%, respectively. Clinical and molecular epidemiological correlation involving the opportunistic pathogen C. albicans may be attributed dependently of each method of genotyping (i.e., MLEE, EK, and SSRs) supplemented with similarity and grouping analysis. Therefore, the use of genotyping systems that give results which offer minimum disparity, or the combination of the results of these systems, can provide greater security and consistency in the determination of strains and their genetic relationships. (C) 2010 Elsevier B.V. All rights reserved.

Runoff sources and land cover change in the Amazon: an end-member mixing analysis from small watersheds

The flowpaths by which water moves from watersheds to streams has important consequences for the runoff dynamics and biogeochemistry of surface waters in the Amazon Basin. The clearing of Amazon forest to cattle pasture has the potential to change runoff sources to streams by shifting runoff to more surficial flow pathways. We applied end-member mixing analysis (EMMA) to 10 small watersheds throughout the Amazon in which solute composition of streamwater and groundwater, overland flow, soil solution, throughfall and rainwater were measured, largely as part of the Large-Scale Biosphere-Atmosphere Experiment in Amazonia. We found a range in the extent to which streamwater samples fell within the mixing space determined by potential flowpath end-members, suggesting that some water sources to streams were not sampled. The contribution of overland flow as a source of stream flow was greater in pasture watersheds than in forest watersheds of comparable size. Increases in overland flow contribution to pasture streams ranged in some cases from 0% in forest to 27-28% in pasture and were broadly consistent with results from hydrometric sampling of Amazon forest and pasture watersheds that indicate 17- to 18-fold increase in the overland flow contribution to stream flow in pastures. In forest, overland flow was an important contribution to stream flow (45-57%) in ephemeral streams where flows were dominated by stormflow. Overland flow contribution to stream flow decreased in importance with increasing watershed area, from 21 to 57% in forest and 60-89% in pasture watersheds of less than 10 ha to 0% in forest and 27-28% in pastures in watersheds greater than 100 ha. Soil solution contributions to stream flow were similar across watershed area and groundwater inputs generally increased in proportion to decreases in overland flow. Application of EMMA across multiple watersheds indicated patterns across gradients of stream size and land cover that were consistent with patterns determined by detailed hydrometric sampling.

A simulation-based evolutionary multiobjective approach to manufacturing cell formation

The purpose of this paper is to propose a multiobjective optimization approach for solving the manufacturing cell formation problem, explicitly considering the performance of this said manufacturing system. Cells are formed so as to simultaneously minimize three conflicting objectives, namely, the level of the work-in-process, the intercell moves and the total machinery investment. A genetic algorithm performs a search in the design space, in order to approximate to the Pareto optimal set. The values of the objectives for each candidate solution in a population are assigned by running a discrete-event simulation, in which the model is automatically generated according to the number of machines and their distribution among cells implied by a particular solution. The potential of this approach is evaluated via its application to an illustrative example, and a case from the relevant literature. The obtained results are analyzed and reviewed. Therefore, it is concluded that this approach is capable of generating a set of alternative manufacturing cell configurations considering the optimization of multiple performance measures, greatly improving the decision making process involved in planning and designing cellular systems. (C) 2010 Elsevier Ltd. All rights reserved.

Spatial genetic structure of Hymenaea stigonocarpa Mart. ex Hayne assessed with chloroplast microsatellite markers

The use of chloroplast DNA markers (cpDNA) helps to elucidate questions related to ecology, evolution and genetic structure. The knowledge of inter-and intra-population genetic structure allows to design effective conservation and management strategies for tropical tree species. With the aim to help the conservation of Hymenaea stigonocarpa of the Cerrado (Brazilian savanna) in Sao Paulo State, an analysis of the spatial genetic structure (SGS) was conducted in two populations using five universal chloroplast microsatellite loci (cpSSR). The population of 68 trees of H. stigonocarpa in the Ecological Station of Itirapina (ESI) had a single haplotype, indicating a strong founder effect. In turn, the population of 47 trees of H. stigonocarpa in a contiguous area that includes the Ecological Station of Assis and the Assis State Forest (ESA), showed six haplotypes ((n) over cap (h) = 6) with a moderate haplotype diversity ((h) over cap = 0667 + 0094), revealing that it was founded by a small number of maternal lineages. The SGS analysis for the population ESA/ASF, using Moran`s I index, indicated limited seed dispersal. Considering SGS, for ex situ conservation strategies in the population ESA/ASF, seed harvesting should require a minimum distance of 750 m among seed-trees.

Distribution of Aceria guerreronis and Neoseiulus baraki among and within coconut bunches in northeast Brazil

Aceria guerreronis Keifer (Acari: Eriophyidae) is considered a major pest of coconut in many countries in the Americas, Africa and parts of Asia. Neoseiulus baraki Athias-Henriot (Acari: Phytoseiidae) is one of the predatory mites most commonly found in association with A. guerreronis in parts of northeast Brazil. The objective of this work was to study the distribution of A. guerreronis and N. baraki among and within coconut bunches. The hypothesis was tested that A. guerreronis and N. baraki are homogenously distributed over the fruits in a bunch, independent of the fruits` age and position. Five collections of bunches, each corresponding to leaves 12-16 from apex (about 2-6 month-old), were conducted in each of three fields in northeastern Brazil, from February to October, 2007. A total of 1,986 fruits were examined. The number of mites, the percentage of fruits hosting them and the level of damage caused by A. guerreronis were evaluated. The highest density of A. guerreronis was observed on fruits of bunch 4 whereas the highest density of N. baraki was observed on bunch 5. Considering all fruits together, no significant differences were observed between densities of either A. guerreronis or N. baraki among the basal, median and apical thirds of the bunches. In younger bunches, fruits of the apical region tend to have lower densities of both mites than fruits of the basal region. This pattern, in association with a similar pattern for the percentage of fruits hosting N. baraki, suggests that the predator initially reaches the basal bunch region, from where it moves to the apical region. The results of the present study suggest that the pest population reduction in bunches older than bunch 4 could be due to (1) an effect of the predator, (2) reduction of the proportion of undamaged tissues amenable to attack, and/or (3) less favorable characteristics of the fruits to attack by A. guerreronis, as indicated by their increasing lignin content as they get older.

Age-related association of rDNA and telomeres with the nuclear matrix in mouse hepatocytes

Transcribed sequences have been suggested to be associated with the nuclear matrix, differing from non-transcribing sequences, which have been reported to be contained in DNA loops. However, although a dozen of genes have their expression level affected by aging, data on chromatin-nuclear matrix interactions under this physiological condition are still scarce. In the present study, liver imprints from young, adult and old mice were subjected to FISH (fluorescence in situ hybridization) for 45S rDNA and telomeric sequences, with or without a lysis treatment to produce extended chromatin fibres. There was an increased amount of 45S rDNA sequences located in DNA loops as the animals grow older, while telomeric sequences were always observed in DNA loops irrespective of the animal age. We assume that active rRNA genes associate with the nuclear matrix, while DNA loops contain silent sequences. Transcription of each 45S rDNA repeat unit is suggested to be dependent on its interaction with the nuclear matrix.

Molecular and morphological diversity in Japanese rice germplasm

Germplasm molecular and phenotypic characterization is instrumental to its utilization and to genetic variability incorporation into rice breeding programmes. The diversity within 192 Japanese rice accessions was analysed for 22 agro-morphological traits and 24 single sequence repeat markers. A total of 181 alleles were detected, 38 of which were exclusive. The number of alleles/marker ranged from 2 to 16, with an average of 7.54 alleles/locus and the H(e) value ranged from 0.01 to 0.82, with an average of 0.46. The accessions showed diversity at molecular and phenotypic level and few showed also good agronomic performance. Tocher`s method applied on a total-dissimilarity matrix was used to determine cluster formation of 13 diversity groups. Most of the accessions (81%) were clustered within a group, whereas eight accessions (Kyuushuu, Eika Ine, Ishiwari Mochi, Col/Fukui/1965, Ookuma Nishiki, Suzume Shirazu, Iwate Ryoon and Toga) did not cluster with other accessions.

Characterization of new polymorphic functional markers for sugarcane

Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP)[02/01167-1]

Podzolization as a deferralitization process: dynamics and chemistry of ground and surface waters in an Acrisol-Podzol sequence of the upper Amazon Basin

Hydrochemical processes involved in the development of hydromorphic Podzols are a major concern for the upper Amazon Basin because of the extent of the areas affected by such processes and the large amounts of organic carbon and associated metals exported to the rivers. The dynamics and chemical composition of ground and surface waters were studied along an Acrisol-Podzol sequence lying in an open depression of a plateau. Water levels were monitored along the sequence over a period of 2 years by means of piezometers. Water was sampled in zero-tension lysimeters for groundwater and for surface water in the drainage network of the depression. The pH and concentrations of organic carbon and major elements (Si, Fe and Al) were determined. The contrasted changes reported for concentrations of Si, organic carbon and metals (Fe, Al) mainly reflect the dynamics of the groundwater and the weathering conditions that prevail in the soils. Iron is released by the reductive dissolution of Fe oxides, mostly in the Bg horizons of the upslope Acrisols. It moves laterally under the control of hydraulic gradients and migrates through the iron-depleted Podzols where it is exported to the river network. Aluminium is released from the dissolution of Al-bearing minerals (gibbsite and kaolinite) at the margin of the podzolic area but is immobilized as organo-Al complexes in spodic horizons. In downslope positions, the quick recharge of the groundwater and large release of organic compounds lead to acidification and a loss of metals (mainly Al), previously stored in the Podzols.

Heteroplasmy in Hair: Study of Mitochondrial DNA Third Hypervariable Region in Hair and Blood Samples

Mitochondrial DNA (mtDNA) analysis has proved useful for forensic identification especially in cases where nuclear DNA is not available, such as with hair evidence. Heteroplasmy, the presence of more than one type of mtDNA in one individual, is a common situation often reported in the first and second mtDNA hypervariable regions (HV1/HV2), particularly in hair samples. However, there is no data about heteroplasmy frequency in the third mtDNA hypervariable region (HV3). To investigate possible heteroplasmy hotspots, HV3 from hair and blood samples of 100 individuals were sequenced and compared. No point heteroplasmy was observed, but length heteroplasmy was, both in C-stretch and CA repeat. To observe which CA ""alleles"" were present in each tissue, PCR products were cloned and re-sequenced. However, no variation among CA alleles was observed. Regarding forensic practice, we conclude that point heteroplasmy in HV3 is not as frequent as in the HV1/HV2.

Clonal transmission of ESBL-producing Klebsiella spp. at a university hospital in Brazil

The present study was designed to determine the prevalence and extended-spectrum beta-lactamase (ESBL) types in clinical isolates of Klebsiella spp. at a university hospital located in the Brazilian southern region (Ribeirao Preto, Sao Paulo) as well as their antibiotic susceptibility and genetic profiles. This study included 147 non-repeat Klebsiella spp. isolates collected from January to June 2000, of which 23 K. pneumoniae and 8 K. oxytoca were selected as ESBL producers by using the Oxoid combination disk method and Etest ESBL strip. beta-lactamases were characterized by IEF, PCR and sequencing assays using primers for ESBL genes. Antibiotic susceptibility was evaluated by MicroScan system. Dissemination of two major clones of ESBL-producing Klebsiella spp. occurred in the hospital. According to the results obtained in this study there was a clonal spread of CTX-M-producing K. oxytoca in five clinics and dissemination of ESBL-producing K. pneumoniae in the nursery and pediatrics wards.

A spectral clustering algorithm for manufacturing cell formation

A graph clustering algorithm constructs groups of closely related parts and machines separately. After they are matched for the least intercell moves, a refining process runs on the initial cell formation to decrease the number of intercell moves. A simple modification of this main approach can deal with some practical constraints, such as the popular constraint of bounding the maximum number of machines in a cell. Our approach makes a big improvement in the computational time. More importantly, improvement is seen in the number of intercell moves when the computational results were compared with best known solutions from the literature. (C) 2009 Elsevier Ltd. All rights reserved.