Characterization of DNA sequence properties through network and statistical approaches

In this thesis we will see that the DNA sequence is constantly shaped by the interactions with its environment at multiple levels, showing footprints of DNA methylation, of its 3D organization and, in the case of bacteria, of the interaction with the host organisms. In the first chapter, we will see that analyzing the distribution of distances between consecutive dinucleotides of the same type along the sequence, we can detect epigenetic and structural footprints. In particular, we will see that CG distance distribution allows to distinguish among organisms of different biological complexity, depending on how much CG sites are involved in DNA methylation. Moreover, we will see that CG and TA can be described by the same fitting function, suggesting a relationship between the two. We will also provide an interpretation of the observed trend, simulating a positioning process guided by the presence and absence of memory. In the end, we will focus on TA distance distribution, characterizing deviations from the trend predicted by the best fitting function, and identifying specific patterns that might be related to peculiar mechanical properties of the DNA and also to epigenetic and structural processes. In the second chapter, we will see how we can map the 3D structure of the DNA onto its sequence. In particular, we devised a network-based algorithm that produces a genome assembly starting from its 3D configuration, using as inputs Hi-C contact maps. Specifically, we will see how we can identify the different chromosomes and reconstruct their sequences by exploiting the spectral properties of the Laplacian operator of a network. In the third chapter, we will see a novel method for source clustering and source attribution, based on a network approach, that allows to identify host-bacteria interaction starting from the detection of Single-Nucleotide Polymorphisms along the sequence of bacterial genomes.

Dna Extraction From Filter-paper Spots Of Vaginal Samples Collected After Sexual Violence.

To detect the presence of male DNA in vaginal samples collected from survivors of sexual violence and stored on filter paper. A pilot study was conducted to evaluate 10 vaginal samples spotted on sterile filter paper: 6 collected at random in April 2009 and 4 in October 2010. Time between sexual assault and sample collection was 4-48hours. After drying at room temperature, the samples were placed in a sterile envelope and stored for 2-3years until processing. DNA extraction was confirmed by polymerase chain reaction for human β-globin, and the presence of prostate-specific antigen (PSA) was quantified. The presence of the Y chromosome was detected using primers for sequences in the TSPY (Y7/Y8 and DYS14) and SRY genes. β-Globin was detected in all 10 samples, while 2 samples were positive for PSA. Half of the samples amplified the Y7/Y8 and DYS14 sequences of the TSPY gene and 30% amplified the SRY gene sequence of the Y chromosome. Four male samples and 1 female sample served as controls. Filter-paper spots stored for periods of up to 3years proved adequate for preserving genetic material from vaginal samples collected following sexual violence.

Transmission Of Intra-cellular Genetic Information: A System Proposal.

One of the great challenges of the scientific community on theories of genetic information, genetic communication and genetic coding is to determine a mathematical structure related to DNA sequences. In this paper we propose a model of an intra-cellular transmission system of genetic information similar to a model of a power and bandwidth efficient digital communication system in order to identify a mathematical structure in DNA sequences where such sequences are biologically relevant. The model of a transmission system of genetic information is concerned with the identification, reproduction and mathematical classification of the nucleotide sequence of single stranded DNA by the genetic encoder. Hence, a genetic encoder is devised where labelings and cyclic codes are established. The establishment of the algebraic structure of the corresponding codes alphabets, mappings, labelings, primitive polynomials (p(x)) and code generator polynomials (g(x)) are quite important in characterizing error-correcting codes subclasses of G-linear codes. These latter codes are useful for the identification, reproduction and mathematical classification of DNA sequences. The characterization of this model may contribute to the development of a methodology that can be applied in mutational analysis and polymorphisms, production of new drugs and genetic improvement, among other things, resulting in the reduction of time and laboratory costs.

The Staircase Structure of the Southern Brazilian Continental Shelf

We show some evidences that the Southeastern Brazilian Continental Shelf (SBCS) has a devil's staircase structure, with a sequence of scarps and terraces with widths that obey fractal formation rules. Since the formation of these features is linked with the sea-level variations, we say that the sea level changes in an organized pulsating way. Although the proposed approach was applied in a particular region of the Earth, it is suitable to be applied in an integrated way to other shelves around the world, since the analyses favor the revelation of the global sea-level variations. Copyright (C) 2009 M. S. Baptista and L. A. Conti.

A simple procedure for reliability assessment of thin composite plates against buckling

In some circumstances ice floes may be modeled as beams. In general this modeling supposes constant thickness, which contradicts field observations. Action of currents, wind and the sequence of contacts, causes thickness to vary. Here this effect is taken into consideration on the modeling of the behavior of ice hitting inclined walls of offshore platforms. For this purpose, the boundary value problem is first equated. The set of equations so obtained is then transformed into a system of equations, that is then solved numerically. For this sake an implicit solution is developed, using a shooting method, with the accompanying Jacobian. In-plane coupling and the dependency of the boundary terms on deformation, make the problem non-linear and the development particular. Deformation and internal resultants are then computed for harmonic forms of beam profile. Forms of giving some additional generality to the problem are discussed.

Identification of a Highly Antigenic Linear B Cell Epitope within Plasmodium vivax Apical Membrane Antigen 1 (AMA-1)

Apical membrane antigen 1 (AMA-1) is considered to be a major candidate antigen for a malaria vaccine. Previous immunoepidemiological studies of naturally acquired immunity to Plasmodium vivax AMA-1 (PvAMA-1) have shown a higher prevalence of specific antibodies to domain II (DII) of AMA-1. In the present study, we confirmed that specific antibody responses from naturally infected individuals were highly reactive to both full-length AMA-1 and DII. Also, we demonstrated a strong association between AMA-1 and DII IgG and IgG subclass responses. We analyzed the primary sequence of PvAMA-1 for B cell linear epitopes co-occurring with intrinsically unstructured/ disordered regions (IURs). The B cell epitope comprising the amino acid sequence 290-307 of PvAMA-1 (SASDQPTQYEEEMTDYQK), with the highest prediction scores, was identified in domain II and further selected for chemical synthesis and immunological testing. The antigenicity of the synthetic peptide was identified by serological analysis using sera from P. vivax-infected individuals who were knowingly reactive to the PvAMA-1 ectodomain only, domain II only, or reactive to both antigens. Although the synthetic peptide was recognized by all serum samples specific to domain II, serum with reactivity only to the full-length protein presented 58.3% positivity. Moreover, IgG reactivity against PvAMA-1 and domain II after depletion of specific synthetic peptide antibodies was reduced by 18% and 33% (P = 0.0001 for both), respectively. These results suggest that the linear epitope SASDQPTQYEEEMTDYQK is highly antigenic during natural human infections and is an important antigenic region of the domain II of PvAMA-1, suggesting its possible future use in pre-clinical studies.

Unexpected Diversity of Cellular Immune Responses against Nef and Vif in HIV-1-Infected Patients Who Spontaneously Control Viral Replication

Background: HIV-1-infected individuals who spontaneously control viral replication represent an example of successful containment of the AIDS virus. Understanding the anti-viral immune responses in these individuals may help in vaccine design. However, immune responses against HIV-1 are normally analyzed using HIV-1 consensus B 15-mers that overlap by 11 amino acids. Unfortunately, this method may underestimate the real breadth of the cellular immune responses against the autologous sequence of the infecting virus. Methodology and Principal Findings: Here we compared cellular immune responses against nef and vif-encoded consensus B 15-mer peptides to responses against HLA class I-predicted minimal optimal epitopes from consensus B and autologous sequences in six patients who have controlled HIV-1 replication. Interestingly, our analysis revealed that three of our patients had broader cellular immune responses against HLA class I-predicted minimal optimal epitopes from either autologous viruses or from the HIV-1 consensus B sequence, when compared to responses against the 15-mer HIV-1 type B consensus peptides. Conclusion and Significance: This suggests that the cellular immune responses against HIV-1 in controller patients may be broader than we had previously anticipated.

Report of a chimeric origin of transposable elements in a bovine-coding gene

Despite the wide distribution of transposable elements (TEs) in mammalian genomes, part of their evolutionary significance remains to be discovered. Today there is a substantial amount of evidence showing that TEs are involved in the generation of new exons in different species. In the present study, we searched 22,805 genes and reported the occurrence of TE-cassettes in coding sequences of 542 cow genes using the RepeatMasker program. Despite the significant number (542) of genes with TE insertions in exons only 14 (2.6%) of them were translated into protein, which we characterized as chimeric genes. From these chimeric genes, only the FAST kinase domains 3 (FASTKD3) gene, present on chromosome BTA 20, is a functional gene and showed evidence of the exaptation event. The genome sequence analysis showed that the last exon coding sequence of bovine FASTKD3 is similar to 85% similar to the ART2A retrotransposon sequence. In addition, comparison among FASTKD3 proteins shows that the last exon is very divergent from those of Homo sapiens, Pan troglodytes and Canis familiares. We suggest that the gene structure of bovine FASTKD3 gene could have originated by several ectopic recombinations between TE copies. Additionally, the absence of TE sequences in all other species analyzed suggests that the TE insertion is clade-specific, mainly in the ruminant lineage.

EVOLUTION OF FLORAL MORPHOLOGY AND POLLINATION SYSTEM IN BIGNONIEAE (BIGNONIACEAE)

The radiation of angiosperms is associated with shifts among pollination modes that are thought to have driven the diversification of floral forms. However, the exact sequence of evolutionary events that led to such great diversity in floral traits is unknown for most plant groups. Here, we characterize the patterns of evolution of individual floral traits and overall floral morphologies in the tribe Bignonieae (Bignoniaceae). We identified 12 discrete traits that are associated with seven floral types previously described for the group and used a penalized likelihood tree of the tribe to reconstruct the ancestral states of those traits at all nodes of the phylogeny of Bignonieae. In addition, evolutionary correlations among traits were conducted using a maximum likelihood approach to test whether the evolution of individual floral traits followed the correlated patterns of evolution expected under the ""pollination syndrome"" concept. The ancestral Bignonieae flower presented an Anemopaegma-type morphology, which was followed by several parallel shifts in floral morphologies. Those shifts occurred through intermediate stages resulting in mixed floral morphologies as well as directly from the Anemopaegma-type morphology to other floral types. Positive and negative evolutionary correlations among traits fit patterns expected under the pollination syndrome perspective, suggesting that interactions between Bignonieae flowers and pollinators likely played important roles in the diversification of the group as a whole.

LIFE CYCLE AND POPULATION STRUCTURE OF AEGLA PAULENSIS (DECAPODA: ANOMURA: AEGLIDAE)

We describe growth, longevity, sex ratio, reproductive period, and recruitment of Aegla paulensis from Jaragua Stale Park, Sao Paulo, Brazil (23 degrees 27'27.9 '' S; 46 degrees 45'32.3 '' W). The population was sampled monthly (September 2007 through August 2009) with the aid of traps. Over five thousand individuals were captured, sexed, measured (carapace length = CL) and inspected for reproductive traits (females only), and then released back to the sampling site. The pattern of the reproductive cycle was strongly seasonal (austral mid autumn through late winter), with a single recruitment pulse per year. The obtained von Bertalanffy growth equations were CL = 21.25[1-e(-0.041(t + 1.250))] and CL = 16.52[1-e(-0.049(t + 1.823))] for males and females, respectively. Males (mean CL +/- SD = 11.86 +/- 2.79 mm) attain larger sizes than females (mean CL +/- SD = 10.84 +/- 2.36 mm). Aegla paulensis reproduces twice during an estimated life span of 40.2 months for females and 33.9 months for males. Temporal variation of sex ratio showed a distinctive pattern characterized by a sequence of three distinct periods that repeated from one year to another, and which suggested that a behavioral component influence the proportion of sex in adult specimens sampled with traps during reproductive and non-reproductive periods.

Musical genres: beating to the rhythms of different drums

Online music databases have increased significantly as a consequence of the rapid growth of the Internet and digital audio, requiring the development of faster and more efficient tools for music content analysis. Musical genres are widely used to organize music collections. In this paper, the problem of automatic single and multi-label music genre classification is addressed by exploring rhythm-based features obtained from a respective complex network representation. A Markov model is built in order to analyse the temporal sequence of rhythmic notation events. Feature analysis is performed by using two multi-variate statistical approaches: principal components analysis (unsupervised) and linear discriminant analysis (supervised). Similarly, two classifiers are applied in order to identify the category of rhythms: parametric Bayesian classifier under the Gaussian hypothesis (supervised) and agglomerative hierarchical clustering (unsupervised). Qualitative results obtained by using the kappa coefficient and the obtained clusters corroborated the effectiveness of the proposed method.

The mitochondrial genome of the phytopathogenic basidiomycete Moniliophthora perniciosa is 109 kb in size and contains a stable integrated plasmid

We present here the sequence of the mitochondrial genome of the basidiomycete phytopathogenic hemibiotrophic fungus Moniliophthora perniciosa, causal agent of the Witches` Broom Disease in Theobroma cacao. The DNA is a circular molecule of 109103 base pairs, with 31.9 % GC, and is the largest sequenced so far. This size is due essentially to the presence of numerous non-conserved hypothetical ORFs. It contains the 14 genes coding for proteins involved in the oxidative phosphorylation, the two rRNA genes, one ORF coding for a ribosomal protein (rps3), and a set of 26 tRNA genes that recognize codons for all amino acids. Seven homing endonucleases are located inside introns. Except atp8, all conserved known genes are in the same orientation. Phylogenetic analysis based on the cox genes agrees with the commonly accepted fungal taxonomy. An uncommon feature of this mitochondrial genome is the presence of a region that contains a set of four, relatively small, nested, inverted repeats enclosing two genes coding for polymerases with an invertron-type structure and three conserved hypothetical genes interpreted as the stable integration of a mitochondrial linear plasmid. The integration of this plasmid seems to be a recent evolutionary event that could have implications in fungal biology. This sequence is available under GenBank accession number AY376688. (c) 2008 The British Mycological Society. Published by Elsevier Ltd. All rights reserved.

Nickel sorption capacity of ground xylem of Quercus ilex trees and effects of selected ligands present in the xylem sap

In this work the influence of four different ligands present in the xylem sap of Quercus ilex (histidine, citric, oxalic and aspartic acids) on Ni(II) adsorption by xylem was investigated. Grinded xylem was trapped in acrylic columns and solutions of Ni(II), in the absence and presence of the four ligands prepared in KNO(3) 0-1 mol L(-1) at pH 5.5, were percolated through the column. Aliquots of solutions were recovered in the column end for Ni determination by graphite furnace atomic absorption spectrometry (GFAAS). The experimental. data to describe Ni sorption by xylem in both the presence and absence of ligands was better explained by the Freundlich isotherm model. The decreasing affinity order of ligands for Ni was: oxalic acid > citric acid > histidine > aspartic acid. On the other hand, the Ni(II) adsorption by xylem increased following the inverse sequence of ligands. Potentiometric titrations of acidic groups were carried out to elucidate the sorption site groups available in Q. ilex xylem. The potentiometric titration has shown three sorption sites: pK(a) 2.6 (57.7% of the sorption sites), related to monobasic aliphatic carboxylic acids or nitrogen aromatic bases, pK(a) 8.1 (9.6%) and pK(a) 9.9 (32.7%), related to phenolic groups. (C) 2008 Elsevier GmbH. All rights reserved.

Effect of hemicellulose and lignin on enzymatic hydrolysis of cellulose from brewer`s spent grain

Enzymatic hydrolysis of brewer`s spent grain in three different forms: original (untreated), pretreated by dilute acid (cellulignin), and pretreated by a sequence of dilute acid and dilute alkali (cellulose pulp), was studied to verify the effect of hemicellulose and lignin on cellulose conversion into glucose. The hydrolysis was carried out using a commercial cellulase concentrate (Celluclast 1.5 L) in an enzyme/substrate ratio of 45 FPU/g, 2% (w/v) substrate concentration, 45 degrees C for 96 h. According to the results, the cellulose hydrolysis was affected by the presence of hemicellulose and/or lignin in the sample. The cellulose conversion ratio (defined as glucose yield + cellobiose yield) from cellulignin was 3.5-times higher than that from untreated sample, whereas from cellulose pulp such value was 4-times higher, correspondent to 91.8% (glucose yield of 85.6%). This best result was probably due to the strong modification in the material structure caused by the hemicellulose and lignin removal from the sample. As a consequence, the cellulose fibers were separated being more susceptible to the enzymatic attack. It was concluded that the lower the hemicellulose and lignin contents in the sample, the higher the efficiency of cellulose hydrolysis. (C) 2007 Elsevier Inc. All rights reserved.

A Constructive Genetic Algorithm for permutation flowshop scheduling

The general flowshop scheduling problem is a production problem where a set of n jobs have to be processed with identical flow pattern on in machines. In permutation flowshops the sequence of jobs is the same on all machines. A significant research effort has been devoted for sequencing jobs in a flowshop minimizing the makespan. This paper describes the application of a Constructive Genetic Algorithm (CGA) to makespan minimization on flowshop scheduling. The CGA was proposed recently as an alternative to traditional GA approaches, particularly, for evaluating schemata directly. The population initially formed only by schemata, evolves controlled by recombination to a population of well-adapted structures (schemata instantiation). The CGA implemented is based on the NEH classic heuristic and a local search heuristic used to define the fitness functions. The parameters of the CGA are calibrated using a Design of Experiments (DOE) approach. The computational results are compared against some other successful algorithms from the literature on Taillard`s well-known standard benchmark. The computational experience shows that this innovative CGA approach provides competitive results for flowshop scheduling; problems. (C) 2007 Elsevier Ltd. All rights reserved.