Shannon, rényie and tsallis entropy analysis of DNA using phase plane

This paper analyzes DNA information using entropy and phase plane concepts. First, the DNA code is converted into a numerical format by means of histograms that capture DNA sequence length ranging from one up to ten bases. This strategy measures dynamical evolutions from 4 up to 410 signal states. The resulting histograms are analyzed using three distinct entropy formulations namely the Shannon, Rényie and Tsallis definitions. Charts of entropy versus sequence length are applied to a set of twenty four species, characterizing 486 chromosomes. The information is synthesized and visualized by adapting phase plane concepts leading to a categorical representation of chromosomes and species.

Analysis of the genetic polymorphism of Paracoccidioides brasiliensis and Paracoccidioides cerebriformis "Moore" by random amplified polymorphic DNA (RAPD) and 28S ribosomal DNA sequencing: Paracoccidioides cerebriformis revisited

Our purpose was to compare the genetic polymorphism of six samples of P. brasiliensis (113, 339, BAT, T1F1, T3B6, T5LN1), with four samples of P. cerebriformis (735, 741, 750, 361) from the Mycological Laboratory of the Instituto de Medicina Tropical de São Paulo, using Random Amplified Polymorphic DNA Analysis (RAPD). RAPD profiles clearly segregated P. brasiliensis and P. cerebriformis isolates. However, the variation on band patterns among P. cerebriformis isolates was high. Sequencing of the 28S rDNA gene showed nucleotide conservancy among P. cerebriformis isolates, providing basis for taxonomical grouping, and disclosing high divergence to P. brasiliensis supporting that they are in fact two distinct species. Moreover, DNA sequence suggests that P. cerebriformis belongs in fact to the Aspergillus genus.

Single line for assembly just-in-sequence multiple models

Programa Doutoral em Líderes para as Indústrias Tecnológicas

Molecular Epidemiology of Human Polyomavirus JC in the Biaka Pygmies and Bantu of Central Africa

Polyomavirus JC (JCV) is ubiquitous in humans and causes a chronic demyelinating disease of the central nervous system , progressive multifocal leukoencephalopathy which is common in AIDS. JCV is excreted in urine of 30-70% of adults worldwide. Based on sequence analysis of JCV complete genomes or fragments thereof, JCV can be classified into geographically derived genotypes. Types 1 and 2 are of European and Asian origin respectively while Types 3 and 6 are African in origin. Type 4, a possible recombinant of European and African genotypes (1 and 3) is common in the USA. To delineate the JCV genotypes in an aboriginal African population, random urine samples were collected from the Biaka Pygmies and Bantu from the Central African Republic. There were 43 males and 25 females aged 4-55 years, with an average age of 26 years. After PCR amplification of JCV in urine, products were directly cycle sequenced. Five of 23 Pygmy adults (22%) and four of 20 Bantu adults (20%) were positive for JC viruria. DNA sequence analysis revealed JCV Type 3 (two), Type 6 (two) and one Type 1 variant in Biaka Pygmies. All the Bantu strains were Type 6. Type 3 and 6 strains of JCV are the predominant strains in central Africa. The presence of multiple subtypes of JCV in Biaka Pygmies may be a result of extensive interactions of Pygmies with their African tribal neighbors during their itinerant movements in the equatorial forest.

Restriction fragment length polymorphism of 195 bp repeated satellite dna of Trypanosoma cruzi supports the existence of two phylogenetic groups

The restriction fragment length polymorphism of the 195 bp repeated DNA sequence of Trypanosoma cruzi was analyzed among 23 T. cruzi stocks giving a reliable picture of the whole phylogenetic variability of the species. The profiles observed with the enzymes Hinf I and Hae III were linked together and supported the existence of two groups. Group 1 shows a 195 bp repeated unit (Hinf I) and high molecular weight DNA (Hae III), while group 2 presents a ladder profile for each enzyme, which is a characteristic of tandemly repeated DNA. The two groups, respectively, clustered stocks pertaining to the two principal lineages evidenced by isoenzyme and RAPD markers. The congruence among these three independent genomic markers corroborates the existence of two real phylogenetic lineages in T. cruzi. The specific monomorphic profiles for each major phylogenetic lineage suggest the existence of ancient sexuality and cryptic biological speciation.

Mitochondrial DNA variation of Triatoma infestans populations and its implication on the specific status of T. melanosoma

DNA sequence comparison of 412 base-pairs fragments of the mitochondrial cytochrome B gene was used to infer the genetic structure of nine geographical Triatoma infestans populations and their phylogenetic relationship with T. melanosoma and T. brasiliensis. T. infestans and T. melanosoma were compared by morphometry, allozyme and cytogenetic analyses, as well as subjected to reciprocal crosses, in order to clarify the taxonomic status of the latter. No differences were found to distinguish the two species and the crosses between them yielded progeny. T. infestans populations presented four haplotypes that could be separated in two clusters: one formed by the samples from Bolivia (Andes and Chaco) and the other formed by samples from Argentina and Brazil. Silvatic and domestic T. infestans populations from Bolivia (Andes) were genetically identical.

Cloning and characterization of SmZF1, a gene encoding a Schistosoma mansoni zinc finger protein

The zinc finger motifs (Cys2His2) are found in several proteins playing a role in the regulation of transcripton. SmZF1, a Schistosoma mansoni gene encoding a zinc finger protein was initially isolated from an adult worm cDNA library, as a partial cDNA. The full sequence of the gene was obtained by subcloning and sequencing cDNA and genomic fragments. The collated gene sequence is 2181 nt and the complete cDNA sequence is 705 bp containing the full open reading frame of the gene. Analysis of the genome sequence revealed the presence of three introns interrupting the coding region. The open reading frame theoretically encodes a protein of 164 amino acids, with a calculated molecular mass of 18,667Da. The predicted protein contains three zinc finger motifs, usually present in transcription regulatory proteins. PCR amplification with specific primers for the gene allowed for the detection of the target in egg, cercariae, schistosomulum and adult worm cDNA libraries indicating the expression of the mRNA in these life cycle stages of S. mansoni. This pattern of expression suggests the gene plays a role in vital functions of different life cycle stages of the parasite. Future research will be directed to elucidate the functional role of SmZF1.

Characterization of new Schistosoma mansoni microsatellite loci in sequences obtained from public DNA databases and microsatellite enriched genomic libraries

In the last decade microsatellites have become one of the most useful genetic markers used in a large number of organisms due to their abundance and high level of polymorphism. Microsatellites have been used for individual identification, paternity tests, forensic studies and population genetics. Data on microsatellite abundance comes preferentially from microsatellite enriched libraries and DNA sequence databases. We have conducted a search in GenBank of more than 16,000 Schistosoma mansoni ESTs and 42,000 BAC sequences. In addition, we obtained 300 sequences from CA and AT microsatellite enriched genomic libraries. The sequences were searched for simple repeats using the RepeatMasker software. Of 16,022 ESTs, we detected 481 (3%) sequences that contained 622 microsatellites (434 perfect, 164 imperfect and 24 compounds). Of the 481 ESTs, 194 were grouped in 63 clusters containing 2 to 15 ESTs per cluster. Polymorphisms were observed in 16 clusters. The 287 remaining ESTs were orphan sequences. Of the 42,017 BAC end sequences, 1,598 (3.8%) contained microsatellites (2,335 perfect, 287 imperfect and 79 compounds). The 1,598 BAC end sequences 80 were grouped into 17 clusters containing 3 to 17 BAC end sequences per cluster. Microsatellites were present in 67 out of 300 sequences from microsatellite enriched libraries (55 perfect, 38 imperfect and 15 compounds). From all of the observed loci 55 were selected for having the longest perfect repeats and flanking regions that allowed the design of primers for PCR amplification. Additionally we describe two new polymorphic microsatellite loci.

Antimicrobial susceptibility determined by the E test, Löwenstein-Jensen proportion, and DNA sequencing methods among Mycobacterium tuberculosis isolates discrepancies, preliminary results

Mycobacterium tuberculosis strains resistant to streptomycin (SM), isoniazid (INH), and/or rifampin (RIF) as determined by the conventional Löwenstein-Jensen proportion method (LJPM) were compared with the E test, a minimum inhibitory concentration susceptibility method. Discrepant isolates were further evaluated by BACTEC and by DNA sequence analyses for mutations in genes most often associated with resistance to these drugs (rpsL, katG, inhA, and rpoB). Preliminary discordant E test results were seen in 75% of isolates resistant to SM and in 11% to INH. Discordance improved for these two drugs (63%) for SM and none for INH when isolates were re-tested but worsened for RIF (30%). Despite good agreement between phenotypic results and sequencing analyses, wild type profiles were detected on resistant strains mainly for SM and INH. It should be aware that susceptible isolates according to molecular methods might contain other mechanisms of resistance. Although reproducibility of the LJPM susceptibility method has been established, variable E test results for some M. tuberculosis isolates poses questions regarding its reproducibility particularly the impact of E test performance which may vary among laboratories despite adherence to recommended protocols. Further studies must be done to enlarge the evaluated samples and looked possible mutations outside of the hot spot sequenced gene among discrepant strains.

Detection of Lymnaea columella infection by Fasciola hepatica through Multiplex-PCR

From complete mitochondrial DNA sequence of Fasciola hepatica available in Genbank, specific primers were designed for a conserved and repetitive region of this trematode. A pair of primers was used for diagnosis of infected Lymnaea columella by F. hepatica during the pre-patent period simultaneously with another pair of primers which amplified the internal transcribed spacer (ITS) region of rDNA from L. columella in a single Multiplex-PCR. The amplification generated a ladder band profile specific for F. hepatica. This profile was observed in positive molluscs at different times of infection, including adult worms from the trematode. The Multiplex-PCR technique showed to be a fast and safe tool for fascioliasis diagnosis, enabling the detection of F. hepatica miracidia in L. columella during the pre-patent period and identification of transmission areas.

Determination of DNA persistence length by cryo-electron microscopy. Separation of the static and dynamic contributions to the apparent persistence length of DNA.

Axial deflection of DNA molecules in solution results from thermal motion and intrinsic curvature related to the DNA sequence. In order to measure directly the contribution of thermal motion we constructed intrinsically straight DNA molecules and measured their persistence length by cryo-electron microscopy. The persistence length of such intrinsically straight DNA molecules suspended in thin layers of cryo-vitrified solutions is about 80 nm. In order to test our experimental approach, we measured the apparent persistence length of DNA molecules with natural "random" sequences. The result of about 45 nm is consistent with the generally accepted value of the apparent persistence length of natural DNA sequences. By comparing the apparent persistence length to intrinsically straight DNA with that of natural DNA, it is possible to determine both the dynamic and the static contributions to the apparent persistence length.

Genetic variability of hepatitis A virus strain HAF-203 isolated in Brazil and expression of the VP1 gene in Escherichia coli

The hepatitis A virus (HAV) HAF-203 strain was isolated from an acute case of HAV infection. The primary isolation of HAF-203 in Brazil and its adaptation to the FRhK-4 cell lineage allowed the production of large amounts of viral particles enabling molecular characterization of the first HAV isolate in Brazil. The aim of our study was to determine the nucleotide sequence of the HAF-203 strain genome, compare it to other HAV genomes and highlight its genetic variability. The complete nucleotide sequence of the HAF-203 strain (7472 nucleotides) was compared to those obtained earlier by others for other HAV isolates. These analyses revealed 19 HAF-specific nucleotide sequence differences with 10 amino acid substitutions. Most of the non-conservative changes were located at VP1, 2C, and 3D genes, but the 3B region was the most variable. The availability of HAF-203 complementary DNA was useful for the production of the recombinant VP1 protein, which is a major determinant of viral infectivity. This recombinant protein was shown by enzyme-linked immunoassay and blotting, to be immunogenic and resemble the native protein, therefore suggesting its value as a reagent for incorporation into diagnostic tests.

Comparative genomics of ParaHox clusters of teleost fishes: gene cluster breakup and the retention of gene sets following whole genome duplications.

BACKGROUND: The evolutionary lineage leading to the teleost fish underwent a whole genome duplication termed FSGD or 3R in addition to two prior genome duplications that took place earlier during vertebrate evolution (termed 1R and 2R). Resulting from the FSGD, additional copies of genes are present in fish, compared to tetrapods whose lineage did not experience the 3R genome duplication. Interestingly, we find that ParaHox genes do not differ in number in extant teleost fishes despite their additional genome duplication from the genomic situation in mammals, but they are distributed over twice as many paralogous regions in fish genomes. RESULTS: We determined the DNA sequence of the entire ParaHox C1 paralogon in the East African cichlid fish Astatotilapia burtoni, and compared it to orthologous regions in other vertebrate genomes as well as to the paralogous vertebrate ParaHox D paralogons. Evolutionary relationships among genes from these four chromosomal regions were studied with several phylogenetic algorithms. We provide evidence that the genes of the ParaHox C paralogous cluster are duplicated in teleosts, just as it had been shown previously for the D paralogon genes. Overall, however, synteny and cluster integrity seems to be less conserved in ParaHox gene clusters than in Hox gene clusters. Comparative analyses of non-coding sequences uncovered conserved, possibly co-regulatory elements, which are likely to contain promoter motives of the genes belonging to the ParaHox paralogons. CONCLUSION: There seems to be strong stabilizing selection for gene order as well as gene orientation in the ParaHox C paralogon, since with a few exceptions, only the lengths of the introns and intergenic regions differ between the distantly related species examined. The high degree of evolutionary conservation of this gene cluster's architecture in particular - but possibly clusters of genes more generally - might be linked to the presence of promoter, enhancer or inhibitor motifs that serve to regulate more than just one gene. Therefore, deletions, inversions or relocations of individual genes could destroy the regulation of the clustered genes in this region. The existence of such a regulation network might explain the evolutionary conservation of gene order and orientation over the course of hundreds of millions of years of vertebrate evolution. Another possible explanation for the highly conserved gene order might be the existence of a regulator not located immediately next to its corresponding gene but further away since a relocation or inversion would possibly interrupt this interaction. Different ParaHox clusters were found to have experienced differential gene loss in teleosts. Yet the complete set of these homeobox genes was maintained, albeit distributed over almost twice the number of chromosomes. Selection due to dosage effects and/or stoichiometric disturbance might act more strongly to maintain a modal number of homeobox genes (and possibly transcription factors more generally) per genome, yet permit the accumulation of other (non regulatory) genes associated with these homeobox gene clusters.

Molecular evidence for a functional ecdysone signaling system in Brugia malayi.

BACKGROUND: Filarial nematodes, including Brugia malayi, the causative agent of lymphatic filariasis, undergo molting in both arthropod and mammalian hosts to complete their life cycles. An understanding of how these parasites cross developmental checkpoints may reveal potential targets for intervention. Pharmacological evidence suggests that ecdysteroids play a role in parasitic nematode molting and fertility although their specific function remains unknown. In insects, ecdysone triggers molting through the activation of the ecdysone receptor: a heterodimer of EcR (ecdysone receptor) and USP (Ultraspiracle). METHODS AND FINDINGS: We report the cloning and characterization of a B. malayi EcR homologue (Bma-EcR). Bma-EcR dimerizes with insect and nematode USP/RXRs and binds to DNA encoding a canonical ecdysone response element (EcRE). In support of the existence of an active ecdysone receptor in Brugia we also cloned a Brugia rxr (retinoid X receptor) homolog (Bma-RXR) and demonstrate that Bma-EcR and Bma-RXR interact to form an active heterodimer using a mammalian two-hybrid activation assay. The Bma-EcR ligand-binding domain (LBD) exhibits ligand-dependent transactivation via a GAL4 fusion protein combined with a chimeric RXR in mammalian cells treated with Ponasterone-A or a synthetic ecdysone agonist. Furthermore, we demonstrate specific up-regulation of reporter gene activity in transgenic B. malayi embryos transfected with a luciferase construct controlled by an EcRE engineered in a B. malayi promoter, in the presence of 20-hydroxy-ecdysone. CONCLUSIONS: Our study identifies and characterizes the two components (Bma-EcR and Bma-RXR) necessary for constituting a functional ecdysteroid receptor in B. malayi. Importantly, the ligand binding domain of BmaEcR is shown to be capable of responding to ecdysteroid ligands, and conversely, ecdysteroids can activate transcription of genes downstream of an EcRE in live B. malayi embryos. These results together confirm that an ecdysone signaling system operates in B. malayi and strongly suggest that Bma-EcR plays a central role in it. Furthermore, our study proposes that existing compounds targeting the insect ecdysone signaling pathway should be considered as potential pharmacological agents against filarial parasites.

Next-generation sequencing of experimental mouse strains.

Since the turn of the century the complete genome sequence of just one mouse strain, C57BL/6J, has been available. Knowing the sequence of this strain has enabled large-scale forward genetic screens to be performed, the creation of an almost complete set of embryonic stem (ES) cell lines with targeted alleles for protein-coding genes, and the generation of a rich catalog of mouse genomic variation. However, many experiments that use other common laboratory mouse strains have been hindered by a lack of whole-genome sequence data for these strains. The last 5 years has witnessed a revolution in DNA sequencing technologies. Recently, these technologies have been used to expand the repertoire of fully sequenced mouse genomes. In this article we review the main findings of these studies and discuss how the sequence of mouse genomes is helping pave the way from sequence to phenotype. Finally, we discuss the prospects for using de novo assembly techniques to obtain high-quality assembled genome sequences of these laboratory mouse strains, and what advances in sequencing technologies may be required to achieve this goal.