Chain motifs: The tails and handles of complex networks

A great part of the interest in complex networks has been motivated by the presence of structured, frequently nonuniform, connectivity. Because diverse connectivity patterns tend to result in distinct network dynamics, and also because they provide the means to identify and classify several types of complex network, it becomes important to obtain meaningful measurements of the local network topology. In addition to traditional features such as the node degree, clustering coefficient, and shortest path, motifs have been introduced in the literature in order to provide complementary descriptions of the network connectivity. The current work proposes a different type of motif, namely, chains of nodes, that is, sequences of connected nodes with degree 2. These chains have been subdivided into cords, tails, rings, and handles, depending on the type of their extremities (e.g., open or connected). A theoretical analysis of the density of such motifs in random and scale-free networks is described, and an algorithm for identifying these motifs in general networks is presented. The potential of considering chains for network characterization has been illustrated with respect to five categories of real-world networks including 16 cases. Several interesting findings were obtained, including the fact that several chains were observed in real-world networks, especially the world wide web, books, and the power grid. The possibility of chains resulting from incompletely sampled networks is also investigated.

Nonuniversal behavior for aperiodic interactions within a mean-field approximation

We study the spin-1/2 Ising model on a Bethe lattice in the mean-field limit, with the interaction constants following one of two deterministic aperiodic sequences, the Fibonacci or period-doubling one. New algorithms of sequence generation were implemented, which were fundamental in obtaining long sequences and, therefore, precise results. We calculate the exact critical temperature for both sequences, as well as the critical exponents beta, gamma, and delta. For the Fibonacci sequence, the exponents are classical, while for the period-doubling one they depend on the ratio between the two exchange constants. The usual relations between critical exponents are satisfied, within error bars, for the period-doubling sequence. Therefore, we show that mean-field-like procedures may lead to nonclassical critical exponents.

Molecular determinants of improved cathepsin B inhibition by new cystatins obtained by DNA shuffling

Background: Cystatins are inhibitors of cysteine proteases. The majority are only weak inhibitors of human cathepsin B, which has been associated with cancer, Alzheimer's disease and arthritis. Results: Starting from the sequences of oryzacystatin-1 and canecystatin-1, a shuffling library was designed and a hybrid clone obtained, which presented higher inhibitory activity towards cathepsin B. This clone presented two unanticipated point mutations as well as an N-terminal deletion. Reversing each point mutation independently or both simultaneously abolishes the inhibitory activity towards cathepsin B. Homology modeling together with experimental studies of the reverse mutants revealed the likely molecular determinants of the improved inhibitory activity to be related to decreased protein stability. Conclusion: A combination of experimental approaches including gene shuffling, enzyme assays and reverse mutation allied to molecular modeling has shed light upon the unexpected inhibitory properties of certain cystatin mutants against Cathepsin B. We conclude that mutations disrupting the hydrophobic core of phytocystatins increase the flexibility of the N-terminus, leading to an increase in inhibitory activity. Such mutations need not affect the inhibitory site directly but may be observed distant from it and manifest their effects via an uncoupling of its three components as a result of increased protein flexibility.

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.

TESTING STATISTICAL HYPOTHESIS ON RANDOM TREES AND APPLICATIONS TO THE PROTEIN CLASSIFICATION PROBLEM

Efficient automatic protein classification is of central importance in genomic annotation. As an independent way to check the reliability of the classification, we propose a statistical approach to test if two sets of protein domain sequences coming from two families of the Pfam database are significantly different. We model protein sequences as realizations of Variable Length Markov Chains (VLMC) and we use the context trees as a signature of each protein family. Our approach is based on a Kolmogorov-Smirnov-type goodness-of-fit test proposed by Balding et at. [Limit theorems for sequences of random trees (2008), DOI: 10.1007/s11749-008-0092-z]. The test statistic is a supremum over the space of trees of a function of the two samples; its computation grows, in principle, exponentially fast with the maximal number of nodes of the potential trees. We show how to transform this problem into a max-flow over a related graph which can be solved using a Ford-Fulkerson algorithm in polynomial time on that number. We apply the test to 10 randomly chosen protein domain families from the seed of Pfam-A database (high quality, manually curated families). The test shows that the distributions of context trees coming from different families are significantly different. We emphasize that this is a novel mathematical approach to validate the automatic clustering of sequences in any context. We also study the performance of the test via simulations on Galton-Watson related processes.

THE LOWER CENTRAL AND DERIVED SERIES OF THE BRAID GROUPS OF THE SPHERE

In this paper, we determine the lower central and derived series for the braid groups of the sphere. We are motivated in part by the study of Fadell-Neuwirth short exact sequences, but the problem is important in its own right. The braid groups of the 2-sphere S(2) were studied by Fadell, Van Buskirk and Gillette during the 1960s, and are of particular interest due to the fact that they have torsion elements (which were characterised by Murasugi). We first prove that for all n epsilon N, the lower central series of the n-string braid group B(n)(S(2)) is constant from the commutator subgroup onwards. We obtain a presentation of Gamma(2)(Bn(S(2))), from which we observe that Gamma(2)(B(4)(S(2))) is a semi-direct product of the quaternion group Q(8) of order 8 by a free group F(2) of rank 2. As for the derived series of Bn(S(2)), we show that for all n >= 5, it is constant from the derived subgroup onwards. The group Bn(S(2)) being finite and soluble for n <= 3, the critical case is n = 4 for which the derived subgroup is the above semi-direct product Q(8) (sic) F(2). By proving a general result concerning the structure of the derived subgroup of a semi-direct product, we are able to determine completely the derived series of B(4)(S(2)) which from (B(4)(S(2)))(4) onwards coincides with that of the free group of rank 2, as well as its successive derived series quotients.

Environmental stresses inhibit splicing in the aquatic fungus Blastocladiella emersonii

Background: Exposure of cells to environmental stress conditions can lead to the interruption of several intracellular processes, in particular those performed by macromolecular complexes such as the spliceosome. Results: During nucleotide sequencing of cDNA libraries constructed using RNA isolated from B. emersonii cells submitted to heat shock and cadmium stress, a large number of ESTs with retained introns was observed. Among the 6,350 ESTs obtained through sequencing of stress cDNA libraries, 181 ESTs presented putative introns (2.9%), while sequencing of cDNA libraries from unstressed B. emersonii cells revealed only 0.2% of ESTs containing introns. These data indicate an enrichment of ESTs with introns in B. emersonii stress cDNA libraries. Among the 85 genes corresponding to the ESTs that retained introns, 19 showed more than one intron and three showed three introns, with intron length ranging from 55 to 333 nucleotides. Canonical splicing junctions were observed in most of these introns, junction sequences being very similar to those found in introns from genes previously characterized in B. emersonii, suggesting that inhibition of splicing during stress is apparently a random process. Confirming our observations, analyses of gpx3 and hsp70 mRNAs by Northern blot and S1 protection assays revealed a strong inhibition of intron splicing in cells submitted to cadmium stress. Conclusion: In conclusion, data indicate that environmental stresses, particularly cadmium treatment, inhibit intron processing in B. emersonii, revealing a new adaptive response to cellular exposure to this heavy metal.

Cross-Talk Between Mitochondria and NADPH Oxidase: Effects of Mild Mitochondrial Dysfunction on Angiotensin II-Mediated Increase in Nox Isoform Expression and Activity in Vascular Smooth Muscle Cells

Mitochondria and NADPH oxidase activation are concomitantly involved in pathogenesis of many vascular diseases. However, possible cross-talk between those ROS-generating systems is unclear. We induced mild mitochondrial dysfunction due to mitochondrial DNA damage after 24 h incubation of rabbit aortic smooth muscle (VSMC) with 250 ng/mL ethidium bromide (EtBr). VSMC remained viable and had 29% less oxygen consumption, 16% greater baseline hydrogen peroxide, and unchanged glutathione levels. Serum-stimulated proliferation was unaltered at 24 h. Although PCR amplification of several mtDNA sequences was preserved, D-Loop mtDNA region showed distinct amplification of shorter products after EtBr. Such evidence for DNA damage was further enhanced after angiotensin-II (AngII) incubation. Remarkably, the normally observed increase in VSMC membrane fraction NADPH oxidase activity after AngII was completely abrogated after EtBr, together with failure to upregulate Nox1 mRNA expression. Conversely, basal Nox4 mRNA expression increased 1.6-fold, while being unresponsive to AngII. Similar loss in AngII redox response occurred after 24 h antimycin-A incubation. Enhanced Nox4 expression was unassociated with endoplasmic reticulum stress markers. Protein disulfide isomerase, an NADPH oxidase regulator, exhibited increased expression and inverted pattern of migration to membrane fraction after EtBr. These results unravel functionally relevant cross-talk between mitochondria and NADPH oxidase, which markedly affects redox responses to AngII. Antioxid Redox Signal 11, 1265-1278.

Haplotype distribution of five nuclear genes based on network genealogies and Bayesian inference indicates that Trypanosoma cruzi hybrid strains are polyphyletic

Chagas disease is still a major public health problem in Latin America. Its causative agent, Trypanosoma cruzi, can be typed into three major groups, T. cruzi I, T. cruzi II and hybrids. These groups each have specific genetic characteristics and epidemiological distributions. Several highly virulent strains are found in the hybrid group; their origin is still a matter of debate. The null hypothesis is that the hybrids are of polyphyletic origin, evolving independently from various hybridization events. The alternative hypothesis is that all extant hybrid strains originated from a single hybridization event. We sequenced both alleles of genes encoding EF-1 alpha, actin and SSU rDNA of 26 T. cruzi strains and DHFR-TS and TR of 12 strains. This information was used for network genealogy analysis and Bayesian phylogenies. We found T. cruzi I and T. cruzi II to be monophyletic and that all hybrids had different combinations of T. cruzi I and T. cruzi II haplotypes plus hybrid-specific haplotypes. Bootstrap values (networks) and posterior probabilities (Bayesian phylogenies) of clades supporting the monophyly of hybrids were far below the 95% confidence interval, indicating that the hybrid group is polyphyletic. We hypothesize that T. cruzi I and T. cruzi II are two different species and that the hybrids are extant representatives of independent events of genome hybridization, which sporadically have sufficient fitness to impact on the epidemiology of Chagas disease.

GPIomics: global analysis of glycosylphosphatidylinositol-anchored molecules of Trypanosoma cruzi

Glycosylphosphatidylinositol (GPI) anchoring is a common, relevant posttranslational modification of eukaryotic surface proteins. Here, we developed a fast, simple, and highly sensitive (high attomole-low femtomole range) method that uses liquid chromatography-tandem mass spectrometry (LC-MS(n)) for the first large-scale analysis of GPI-anchored molecules (i.e., the GPIome) of a eukaryote, Trypanosoma cruzi, the etiologic agent of Chagas disease. Our genome-wise prediction analysis revealed that approximately 12% of T. cruzi genes possibly encode GPI-anchored proteins. By analyzing the GPIome of T. cruzi insect-dwelling epimastigote stage using LC-MS(n), we identified 90 GPI species, of which 79 were novel. Moreover, we determined that mucins coded by the T. cruzi small mucin-like gene (TcSMUG S) family are the major GPI-anchored proteins expressed on the epimastigote cell surface. TcSMUG S mucin mature sequences are short (56-85 amino acids) and highly O-glycosylated, and contain few proteolytic sites, therefore, less likely susceptible to proteases of the midgut of the insect vector. We propose that our approach could be used for the high throughput GPIomic analysis of other lower and higher eukaryotes. Molecular Systems Biology 7 April 2009; doi:10.1038/msb.2009.13

Molecular and morphological characterization of Amathia distans Busk and Amathia brasiliensis Busk (Bryozoa: Ctenostomata) from the tropical and subtropical Western Atlantic

Morphological and molecular analyses have proven to be complementary tools of taxonomic information for the redescription of the ctenostome bryozoans Amathia brasiliensis Busk, 1886 and Amathia distans Busk, 1886. The two species, originally described from material collected by the `Challenger` expedition but synonymized by later authors, now have their status fixed by means of the selection of lectotypes, morphological observations and analyses of DNA sequences described here. The morphological characters allowing the identification of living and/or preserved specimens are (1) A. brasiliensis: whitish-pale pigment spots in the frontal surface of stolons and zooids, and a wide stolon with biserial zooid clusters growing in clockwise and anti-clockwise spirals along it, the spirality direction being maintained from maternal to daughter stolons; and (2) A. distans: bright yellow pigment spots in stolonal and zooidal surfaces including lophophores, and a slender stolon, thickly cuticularized, with biserial zooid clusters growing in clockwise and anti-clockwise spirals along it and the spirality direction not maintained from maternal to daughter stolons. Pairwise comparisons of DNA sequences of the mitochondrial genes cytochrome c oxidase subunit I and large ribosomal RNA subunit revealed deep genetic divergence between A. brasiliensis and A. distans. Finally, analyses of those sequences within a Bayesian phylogenetic context recovered their genealogical species status.

Diversity of free-living marine nematodes (Enoplida) from Baja California assessed by integrative taxonomy

We used morphological and molecular approaches to evaluate the diversity of free-living marine nematodes (order Enoplida) at four coastal sites in the Gulf of California and three on the Pacific coast of Baja California, Mexico. We identified 22 morphological species belonging to six families, of which Thoracostomopsidae and Oncholaimidae were the most diverse. The genus Mesacanthion (Thoracostomopsidae) was the most widespread and diverse. Five allopatric species, genetically and morphologically differentiated, were found in two localities in the Gulf of California (M. sp1 and M. sp2) and three in the Pacific coast (M. sp3, M. sp4 and M. sp5). Overall, we produced 19 and 20 sequences for the 18S and 28S genes, respectively. Neither gene displayed intraspecific polymorphisms, which allowed us to establish that some morphological variation was likely either ontogenetic or due to phenotypic plasticity. Although 18S and 28S phylogenies were topologically congruent (incongruence length difference test, P > 0.05), divergences between species were much higher in the 28S gene. Moreover, this gene possessed a stronger phylogenetic signal to resolve relationships involving Rhabdodemania and Bathylaimus. On the other hand, the close relationship of Pareurystomina (Enchilidiidae) with oncholaimids warrants further study. The 28S sequences (D2D3 domain) may be better suited for DNA barcoding of marine nematodes than those from the 18S rDNA, particularly for differentiating closely related or cryptic species. Finally, our results underline the relevance of adopting an integrative approach encompassing morphological and molecular analyses to improve the assessment of marine nematode diversity and advance their taxonomy.

Influence of the bacterioplankton community of a tropical eutrophic lagoon on the bacterial community of its neighbouring ocean

The Rodrigo de Freitas lagoon (RFL) is a tropical eutrophic coastal ecosystem located in the urban area of Rio de Janeiro, Brazil. This environment consists of freshwater but has communication with the ocean through a channel (Jardim de Alah`s Channel). The aim of this study was to evaluate the influence of lagoon water on the nearby ocean using molecular and traditional microbiological methods. We hypothesised that due to the eutrophic low-salinity environment, the bacterioplankton community from the RFL would have a native ""brackish"" composition influenced by both freshwater and marine phylotypes, and that bacterial phylotypes of this community would be detected in oceanic samples closer to the channel between the lagoon and the ocean. The cultivation and microscopy experiments clearly showed this influence. Bacterial cell counts revealed that the greater amounts of bacterial cells present in the lagoon increased the observed values seen at oceanic stations near the channel. The Denaturing gradient gel eletrophoresis community profiles also showed a clear influence of Rodrigo de Freitas lagoon waters on the adjacent beaches. The band patterns found for the stations near the channel showed that these communities were mixtures of the communities of the lagoon and sea, and as the distance from the channel increased, the samples became more similar to ocean bacterial communities. A 16S rRNA gene clone library was constructed using a sample acquired from the connection point between the lagoon and the ocean. Around 52% of the sequences in the library showed similarity to the genus Proteobacteria (1% Alpha, 21% Beta, 19% Gamma and 29% unclassified Proteobacteria), and the second most abundant genus was Bacteroidetes, with 15% of the total clones. The results showed that the structure of the bacterial community had both freshwater and marine characteristics.

A NOVEL EPIPHYTIC CYANOBACTERIAL SPECIES FROM THE GENUS BRASILONEMA CAUSING DAMAGE TO EUCALYPTUS LEAVES

A cyanobacterial mat colonizing the leaves of Eucalyptus grandis was determined to be responsible for serious damage affecting the growth and development of whole plants under the clonal hybrid nursery conditions. The dominant cyanobacterial species was isolated in BG-11 medium lacking a source of combined nitrogen and identified by cell morphology characters and molecular phylogenetic analysis (16S rRNA gene and cpcBA-IGS sequences). The isolated strain represents a novel species of the genus Brasilonema and is designated Brasilonema octagenarum strain UFV-E1. Thin sections of E. grandis leaves analyzed by light and electron microscopy showed that the B. octagenarum UFV-E1 filaments penetrate into the leaf mesophyll. The depth of infection and the mechanism by which the cyanobacterium invades leaf tissue were not determined. A major consequence of colonization by this cyanobacterium is a reduction in photosynthesis in the host since the cyanobacterial mats decrease the amount of light incident on leaf surfaces. Moreover, the cyanobacteria also interfere with stomatal gas exchange, decreasing CO2 assimilation. To our knowledge, this is the first report of an epiphytic cyanobacterial species causing damage to E. grandis leaves.

Highly Toxic Microcystis aeruginosa Strain, Isolated from So Paulo-Brazil, Produce Hepatotoxins and Paralytic Shellfish Poison Neurotoxins

While evaluating several laboratory-cultured cyanobacteria strains for the presence of paralytic shellfish poison neurotoxins, the hydrophilic extract of Microcystis aeruginosa strain SPC777-isolated from Billings`s reservoir, So Paulo, Brazil-was found to exhibit lethal neurotoxic effect in mouse bioassay. The in vivo test showed symptoms that unambiguously were those produced by PSP. In order to identify the presence of neurotoxins, cells were lyophilized, and the extracts were analyzed by HPLC-FLD and HPLC-MS. HPLC-FLD analysis revealed four main Gonyautoxins: GTX4(47.6%), GTX2(29.5%), GTX1(21.9%), and GTX3(1.0%). HPLC-MS analysis, on other hand, confirmed both epimers, with positive Zwitterions M(+) 395.9 m/z for GTX3/GTX2 and M(+) 411 m/z for GTX4/GTX1 epimers. The hepatotoxins (Microcystins) were also evaluated by ELISA and HPLC-MS analyses. Positive immunoreaction was observed by ELISA assay. Alongside, the HPLC-MS analyses revealed the presence of [l-ser(7)] MCYST-RR. The N-methyltransferase (NMT) domain of the microcystin synthetase gene mcyA was chosen as the target sequence to detect the presence of the mcy gene cluster. PCR amplification of the NMT domain, using the genomic DNA of the SPC777 strain and the MSF/MSR primer set, resulted in the expected 1,369 bp product. The phylogenetic analyses grouped the NMT sequence with the NMT sequences of other known Microcystis with high bootstrap support. The taxonomical position of M. aeruginosa SPC777 was confirmed by a detailed morphological description and a phylogenetic analysis of 16S rRNA gene sequence. Therefore, co-production of PSP neurotoxins and microcystins by an isolated M. aeruginosa strain is hereby reported for the first time.