A two-phase strategy for detecting recombination in nucleotide sequences

Genetic recombination can produce heterogeneous phylogenetic histories within a set of homologous genes. Delineating recombination events is important in the study of molecular evolution, as inference of such events provides a clearer picture of the phylogenetic relationships among different gene sequences or genomes. Nevertheless, detecting recombination events can be a daunting task, as the performance of different recombination-detecting approaches can vary, depending on evolutionary events that take place after recombination. We previously evaluated the effects of post-recombination events on the prediction accuracy of recombination-detecting approaches using simulated nucleotide sequence data. The main conclusion, supported by other studies, is that one should not depend on a single method when searching for recombination events. In this paper, we introduce a two-phase strategy, applying three statistical measures to detect the occurrence of recombination events, and a Bayesian phylogenetic approach to delineate breakpoints of such events in nucleotide sequences. We evaluate the performance of these approaches using simulated data, and demonstrate the applicability of this strategy to empirical data. The two-phase strategy proves to be time-efficient when applied to large datasets, and yields high-confidence results.

Knowledge discovery and data mining in biological databases

The new technologies for Knowledge Discovery from Databases (KDD) and data mining promise to bring new insights into a voluminous growing amount of biological data. KDD technology is complementary to laboratory experimentation and helps speed up biological research. This article contains an introduction to KDD, a review of data mining tools, and their biological applications. We discuss the domain concepts related to biological data and databases, as well as current KDD and data mining developments in biology.

Analysis of Mouse Keratin 6a Regulatory Sequences in Transgenic Mice Reveals Constitutive, Tissue-Specific Expression by a Keratin 6a Minigene

The analysis of keratin 6 expression is complicated by the presence of multiple isoforms that are expressed constitutively in a number of internal stratified epithelia, in palmoplantar epidermis, and in the companion cell layer of the hair follicle. In addition, keratin 6 expression is inducible in interfollicular epidermis and the outer root sheath of the follicle, in response to wounding stimuli, phorbol esters, or retinoic acid. In order to establish the critical regions involved in the regulation of keratin 6a (the dominant isoform in mice), we generated transgenic mice with two different-sized mouse keratin 6a constructs containing either 1.3 kb or 0.12 kb of 5' flanking sequence linked to the lacZ reporter gene. Both constructs also contained the first intron and the 3' flanking sequence of mouse keratin 6a. Ectopic expression of either transgene was not observed. Double-label immunofluorescence analyses demonstrated expression of the reporter gene in keratin 6 expressing tissues, including the hair follicle, tongue, footpad, and nail bed, showing that both transgenes retained keratinocyte-specific expression. Quantitative analysis of beta -galactosidase activity verified that both the 1.3 and 0.12 kb keratin 6a promoter constructs produced similar levels of the reporter. Notably, both constructs were constitutively expressed in the outer root sheath and interfollicular epidermis in the absence of any activating stimulus, suggesting that they lack the regulatory elements that normally silence transcription in these cells. This study has revealed that a keratin 6a minigene contains critical cis elements that mediate tissue-specific expression and that the elements regulating keratin 6 induction lie distal to the 1.3 kb promoter region.

Identification of "pathologs" (disease-related genes) from the RIKEN mouse cDNA dataset using human curation plus FACTS, a new biological information extraction system

Background: A major goal in the post-genomic era is to identify and characterise disease susceptibility genes and to apply this knowledge to disease prevention and treatment. Rodents and humans have remarkably similar genomes and share closely related biochemical, physiological and pathological pathways. In this work we utilised the latest information on the mouse transcriptome as revealed by the RIKEN FANTOM2 project to identify novel human disease-related candidate genes. We define a new term patholog to mean a homolog of a human disease-related gene encoding a product ( transcript, anti-sense or protein) potentially relevant to disease. Rather than just focus on Mendelian inheritance, we applied the analysis to all potential pathologs regardless of their inheritance pattern. Results: Bioinformatic analysis and human curation of 60,770 RIKEN full-length mouse cDNA clones produced 2,578 sequences that showed similarity ( 70 - 85% identity) to known human-disease genes. Using a newly developed biological information extraction and annotation tool ( FACTS) in parallel with human expert analysis of 17,051 MEDLINE scientific abstracts we identified 182 novel potential pathologs. Of these, 36 were identified by computational tools only, 49 by human expert analysis only and 97 by both methods. These pathologs were related to neoplastic ( 53%), hereditary ( 24%), immunological ( 5%), cardio-vascular (4%), or other (14%), disorders. Conclusions: Large scale genome projects continue to produce a vast amount of data with potential application to the study of human disease. For this potential to be realised we need intelligent strategies for data categorisation and the ability to link sequence data with relevant literature. This paper demonstrates the power of combining human expert annotation with FACTS, a newly developed bioinformatics tool, to identify novel pathologs from within large-scale mouse transcript datasets.

Phylogenetic analysis of the Monocotylidae (Monogenea) inferred from 28S rDNA sequences

The current classification of the Monocotylidae (Monogenea) is based on a phylogeny generated from morphological characters. The present study tests the morphological phylogenetic hypothesis using molecular methods. Sequences from domains C2 and D1 and the partial domains C1 and D2 from the 28S rDNA gene for 26 species of monocotylids from six of the seven subfamilies were used. Trees were generated using maximum parsimony, neighbour joining and maximum likelihood algorithms. The maximum parsimony tree, with branches showing less than 70% bootstrap support collapsed, had a topology identical to that obtained using the maximum likelihood analysis. The neighbour joining tree, with branches showing less than 70% support collapsed. differed only in its placement of Heterocotyle capricornensis as the sister group to the Decacotylinae clade. The molecular tree largely supports the subfamilies established using morphological characters. Differences are primarily how the subfamilies are related to each other. The monophyly of the Calicotylinae and Merizocotylinae and their sister group relationship is supported by high bootstrap values in all three methods, but relationships within the Merizocotylinae are unclear. Merizocotyle is paraphyletic and our data suggest that Mycteronastes and Thaumatocotyle, which were synonymized with Merizocotyle after the morphological cladistic analysis, should perhaps be resurrected as valid genera. The monophyly of the Monocotylinae and Decacotylinae is also supported by high bootstrap values. The Decacotylinae, which was considered previously to be the sister group to the Calicotylinae plus Merizocotylinae, is grouped in an unresolved polychotomy with the Monocotylinae and members of the Heterocotylinae. According to our molecular data, the Heterocotylinae is paraphyletic. Molecular data support a sister group relationship between Troglocephalus rhinobatidis and Neoheterocotyle rhinobatidis to the exclusion of the other species of Neoheterocotyle and recognition of Troglocephalus renders Neoheterocotyle,le paraphyletic. We propose Troglocephalus incertae sedis. An updated classification and full species list of the Monocotylidae is provided. (C) 2001 Australian Society for Parasitology Inc. Published by Elsevier Science Ltd. All rights reserved.

Phylogenetic analysis of evolutionary relationships of the planctomycete division of the domain bacteria based on amino acid sequences of elongation factor Tu

Sequences from the tuf gene coding for the elongation factor EF-Tu were amplified and sequenced from the genomic DNA of Pirellula marina and Isosphaera pallida, two species of bacteria within the order Planctomycetales. A near-complete (1140-bp) sequence was obtained from Pi. marina and a partial (759-bp) sequence was obtained for I. pallida. Alignment of the deduced Pi. marina EF-Tu amino acid sequence against reference sequences demonstrated the presence of a unique Il-amino acid sequence motif not present in any other division of the domain Bacteria. Pi. marina shared the highest percentage amino acid sequence identity with I. pallida but showed only a low percentage identity with other members of the domain Bacteria. This is consistent with the concept of the planctomycetes as a unique division of the Bacteria. Neither primary sequence comparison of EF-Tu nor phylogenetic analysis supports any close relationship between planctomycetes and the chlamydiae, which has previously been postulated on the basis of 16S rRNA. Phylogenetic analysis of aligned EF-Tu amino acid sequences performed using distance, maximum-parsimony, and maximum likelihood approaches yielded contradictory results with respect to the position of planctomycetes relative to other bacteria, It is hypothesized that long-branch attraction effects due to unequal evolutionary rates and mutational saturation effects may account for some of the contradictions.

Molecular phylogenetics of Lentibulariaceae inferred from plastid rps16 Intron and trnL-F DNA sequences: implications for character evolution and biogeography

Phylogenetic relationships among 75 species of Lentibulariaceae, representing the three recognized genera, were assessed by cladistic analysis of DNA sequences from the plastid rps16 intron and the trnL-F region. Sequence data from the two loci were analyzed both separately and in combination. Consensus trees from all analyses are congruent, and parsimony jackknife results demonstrate strong support for relationships both between and within each of the three demonstrably monophyletic genera. The genus Pinguicula is sister to a Genlisea-Utricularia clade, the phylogenetic structure within this clade closely follows Taylor's recent sectional delimitations based on morphology. Three principal clades are shown within Utricularia, with the basal sections Polypoinpholyx and Pleiochasia together forming the sister lineage of the remaining Utricularia species. Of the fundamental morphological specializations, the stoloniferous growth form apparently arose independently within Genlisea and Utricularia three times, and within Utricularia itself, perhaps more than once. The epiphytic habit has evolved independently at least three times, in Pinguicula, in Utricularia section Phyllaria, and within the two sections Orchidioides and Iperua (in the latter as bromeliad tank-epiphytes). The suspended aquatic habit may have evolved independently within sections Utricularia and Vesiculina. Biogeographic optimization on the phylogeny demonstrates patterns commonly associated with the boreotropics hypothesis and limits the spatial origin of Lentibulariaceae to temperate Eurasia or tropical America.

C-terminal sequences in R-Ras are involved in integrin regulation and in plasma membrane microdomain distribution

The small GTPases R-Ras and H-Ras are highly homologous proteins with contrasting biological properties, for example, they differentially modulate integrin affinity: H-Ras suppresses integrin activation in fibroblasts whereas R-Ras can reverse this effect of H-Ras. To gain insight into the sequences directing this divergent phenotype, we investigated a panel of H-Ras/R-Ras chimeras and found that sequences in the R-Ras hypervariable C-terminal region including amino acids 175-203 are required for the R-Ras ability to increase integrin activation in CHO cells; however, the proline-rich site in this region, previously reported to bind the adaptor protein Nck, was not essential for this effect. In addition, we found that the GTPase TC21 behaved similarly to R-Ras. Because the C-termini of Ras proteins can control their subcellular localization, we compared the localization of H-Ras and R-Ras. In contrast to H-Ras, which migrates out of lipid rafts upon activation, we found that activated R-Ras remained localized to lipid rafts. However, functionally distinct H-Ras/R-Ras chimeras containing different C-terminal R-Ras segments localized to lipid rafts irrespective of their integrin phenotype. (C) 2003 Elsevier Inc. All rights reserved.

The value of idiosyncratic markers and changes to conserved tRNA sequences from the mitochondrial genome of hard ticks (Acari: Ixodida: Ixodidae) for phylogenetic inference

Idiosyncratic markers are features of genes and genomes that are so unusual that it is unlikely that they evolved more than once in a lineage of organisms. Here we explore further the potential of idiosyncratic markers and changes to typically conserved tRNA sequences for phylogenetic inference. Hard ticks were chosen as the model group because their phylogeny has been studied extensively. Fifty-eight candidate markers from hard ticks ( family Ixodidae) and 22 markers from the subfamily Rhipicephalinae sensu lato were mapped onto phylogenies of these groups. Two of the most interesting markers, features of the secondary structure of two different tRNAs, gave strong support to the hypothesis that species of the Prostriata ( Ixodes spp.) are monophyletic. Previous analyses of genes and morphology did not strongly support this relationship, instead suggesting that the Prostriata is paraphyletic with respect to the Metastriata ( the rest of the hard ticks). Parallel or convergent evolution was not found in the arrangements of mitochondrial genes in ticks nor were there any reversals to the ancestral arthropod character state. Many of the markers identified were phylogenetically informative, whereas others should be informative with study of additional taxa. Idiosyncratic markers and changes to typically conserved nucleotides in tRNAs that are phylogenetically informative were common in this data set, and thus these types of markers might be found in other organisms.

Phylogeny of the Sporobolus indicus complex, based on internal transcribed spacer (ITS) sequences

The entire internal transcribed spacer ( ITS) region, including the 5.8S subunit of the nuclear ribosomal DNA ( rDNA), was sequenced by direct double-stranded sequencing of polymerase chain reaction (PCR) amplified fragments. The study included 40 Sporobolus ( Family Poaceae, subfamily Chloridoideae) seed collections from 14 putative species ( all 11 species from the S. indicus complex and three Australian native species). These sequences, along with those from two out-group species [ Pennisetum alopecuroides ( L.) Spreng. and Heteropogon contortus ( L.) P. Beauv. ex Roemer & Schultes, Poaceae, subfamily Panicoideae], were analysed by the parsimony method (PAUP; version 4.0b4a) to infer phylogenetic relationships among these species. The length of the ITS1, 5.8S subunit and ITS2 region were 222, 164 and 218 base pairs ( bp), respectively, in all species of the S. indicus complex, except for the ITS2 region of S. diandrus P. Beauv. individuals, which was 217 bp long. Of the 624 characters included in the analysis, 245 ( 39.3%) of the 330 variable sites contained potential phylogenetic information. Differences in sequences among the members of the S. pyramidalis P. Beauv., S. natalensis (Steud.) Dur & Schinz and S. jacquemontii Kunth. collections were 0%, while differences ranged from 0 to 2% between these and other species of the complex. Similarly, differences in sequences among collections of S. laxus B. K. Simon, S. sessilis B. K. Simon, S. elongatus R. Br. and S. creber De Nardi were 0%, compared with differences of 1-2% between these four species and the rest of the complex. When comparing S. fertilis ( Steud.) Clayton and S. africanus (Poir.) Robyns & Tourney, differences between collections ranged from 0 to 1%. Parsimony analysis grouped all 11 species of the S. indicus complex together, indicating a monophyletic origin. For the entire data set, pair-wise distances among members of the S. indicus complex varied from 0.00 to 1.58%, compared with a range of 20.08-21.44% among species in the complex and the Australian native species studied. A strict consensus phylogenetic tree separated 11 species of the S. indicus complex into five major clades. The phylogeny, based on ITS sequences, was found to be congruent with an earlier study on the taxonomic relationship of the weedy Sporobolus grasses revealed from random amplified polymorphic DNA ( RAPD). However, this cladistic analysis of the complex was not in agreement with that created on past morphological analyses and therefore gives a new insight into the phylogeny of the S. indicus complex.

Spatio-temporal alignment of pedobarographic image sequences

O documento em anexo encontra-se na versão post-print (versão corrigida pelo editor).

Gold electrode modified by self-assembled monolayers of thiols to determine DNA sequences hybridization

The process of immobilization of biological molecules is one of the most important steps in the construction of a biosensor. In the case of DNA, the way it exposes its bases can result in electrochemical signals to acceptable levels. The use of self-assembled monolayer that allows a connection to the gold thiol group and DNA binding to an aldehydic ligand resulted in the possibility of determining DNA hybridization. Immobilized single strand of DNA (ssDNA) from calf thymus pre-formed from alkanethiol film was formed by incubating a solution of 2-aminoethanothiol (Cys) followed by glutaraldehyde (Glu). Cyclic voltammetry (CV) was used to characterize the self-assembled monolayer on the gold electrode and, also, to study the immobilization of ssDNA probe and hybridization with the complementary sequence (target ssDNA). The ssDNA probe presents a well-defined oxidation peak at +0.158 V. When the hybridization occurs, this peak disappears which confirms the efficacy of the annealing and the DNA double helix performing without the presence of electroactive indicators. The use of SAM resulted in a stable immobilization of the ssDNA probe, enabling the hybridization detection without labels. This study represents a promising approach for molecular biosensor with sensible and reproducible results.

Characterization and Biological Activities of Ocellatin Peptides from the Skin Secretion of the Frog Leptodactylus pustulatus

Eight new peptides were isolated from the skin secretion of the frog Leptodactylus pustulatus and their amino acid sequences determined by de novo sequencing and by cDNA cloning. Structural similarities between them and other antimicrobial peptides from the skin secretion of Leptodactylus genus frogs were found. Ocellatins-PT1 to -PT5 (25 amino acid residues) are amidated at the C-terminus, while ocellatins-PT6 to -PT8 (32 amino acid residues) have free carboxylates. Antimicrobial activity, hemolytic tests, and cytotoxicity against a murine fibroblast cell line were investigated. All peptides, except for ocellatin-PT2, have antimicrobial activity against at least one Gram negative strain. Ocellatin-PT8 inhibited the growth of Escherichia coli, Staphylococcus aureus, Klebsiella pneumoniae, and Salmonella choleraesuis strains with MICs in the 60−240 μM range. No significant effect was observed in human erythrocytes and in a murine fibroblast cell line after exposure to the peptides at MICs. A comparison between sequences obtained by both direct HPLC-MS de novo sequencing and cDNA cloning demonstrates the secretion of mature peptides derived from a pre-pro-peptide structure.