Sequence-specific reconstruction from fragmentary databases using seed sequences: implementation and validation on SAGE, proteome and generic sequencing data

Motivation: DNA assembly programs classically perform an all-against-all comparison of reads to identify overlaps, followed by a multiple sequence alignment and generation of a consensus sequence. If the aim is to assemble a particular segment, instead of a whole genome or transcriptome, a target-specific assembly is a more sensible approach. GenSeed is a Perl program that implements a seed-driven recursive assembly consisting of cycles comprising a similarity search, read selection and assembly. The iterative process results in a progressive extension of the original seed sequence. GenSeed was tested and validated on many applications, including the reconstruction of nuclear genes or segments, full-length transcripts, and extrachromosomal genomes. The robustness of the method was confirmed through the use of a variety of DNA and protein seeds, including short sequences derived from SAGE and proteome projects.

Intraspecific variation in 16S rRNA gene of Mycoplasma synoviae determined by DNA sequencing

Mycoplasma synoviae (MS) is an important avian pathogen may cause both respiratory disease and joint inflammation synovitis in poultry, causing economic losses to the Brazilian poultry industry. The genotypic variation in 16S rRNA gene is unknown. Partial sequences of 16S rRNA gene of 19 strains of M. synoviae were sequenced and analyzed in order to obtain molecular characterization and evaluation of the genetic variability of strains from distinct Brazilian areas of poultry production. Different polymorphic patterns were observed. The number of polymorphic alterations in the studied strains ranged from 0 to 6. The nucleotide variations, including deletion, insertion and substitutions, ranged from 3 to 5. The genotypic diversity observed in this study may be explained by spontaneous mutations that may occur when a lineage remains in the same flock for long periods. The culling and reposition in poultry flocks may be responsible for the entry of new strains in different areas. (C) 2008 Elsevier Ltd. All rights reserved.

Network genealogy of 195-bp satellite DNA supports the superimposed hybridization hypothesis of Trypanosoma cruzi evolutionary pattern

Trypanosoma cruzi is highly diverse genetically and has been partitioned into six discrete typing units (DTUs), recently re-named T. cruzi I-VI. Although T. cruzi reproduces predominantly by binary division, accumulating evidence indicates that particular DTUs are the result of hybridization events. Two major scenarios for the origin of the hybrid lineages have been proposed. It is accepted widely that the most heterozygous TcV and TcVI DTUs are the result of genetic exchange between TcII and TcIII strains. On the other hand, the participation of a TcI parental in the current genome structure of these hybrid strains is a matter of debate. Here, sequences of the T. cruzi-specific 195-bp satellite DNA of TcI, TcII, Tat, TcV, and TcVI strains have been used for inferring network genealogies. The resulting genealogy showed a high degree of reticulation, which is consistent with more than one event of hybridization between the Tc DTUs. The data also strongly suggest that Tat is a hybrid with two distinct sets of satellite sequences, and that genetic exchange between TcI and TcII parentals occurred within the pedigree of the TcV and TcVI DTUs. Although satellite DNAs belong to the fast-evolving portion of eukaryotic genomes, in >100 satellite units of nine T. cruzi strains we found regions that display 100% identity. No DTU-specific consensus motifs were identified, inferring species-wide conservation. (C) 2010 Elsevier B.V. All rights reserved.

Genomic organization and transcription analysis of the 195-bp satellite DNA in Trypanosoma cruzi

The 195-bp satellite DNA is the most abundant Trypanosoma cruzi repetitive sequence. Here we show by RNA blotting and RT-PCR that 195 SAT is intensely transcribed. We observed a positive correlation between the level of satellite RNA and the abundance of the satellite copies in the genome of T cruzi strains and that the satellite expression is not developmentally regulated. By analyzing CL Brener individual reads, we estimated that 195 SAT corresponds to approximately 5% of the CL Brener genome. 195 SAT elements were found in only 37 annotated contigs, indicating that a large number of satellite copies were not incorporated into the assembled data. The assembled satellite units are distributed in non-syntenic regions with Trypanosoma brucei and Leishmania major genomes, enriched with surface proteins, retroelements, RHS and hypothetical proteins. Satellite repeats were not observed in annotated subtelomeric regions. We report that 12 satellite sequences are truncated by the retroelement VIPER. (C) 2008 Elsevier B.V. All rights reserved.

Plant or fungal sequences? An alternative optimized PCR protocol to avoid ITS (nrDNA) misamplification

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Low variation in ribosomal DNA and internal transcribed spacers of the symbiotic fungi of leaf-cutting ants (Attini: Formicidae)

Leaf-cutting ants of the genera Atta and Acromyrmex (tribe Attini) are symbiotic with basidiomycete fungi of the genus Leucoagaricus (tribe Leucocoprineae), which they cultivate on vegetable matter inside their nests. We determined the variation of the 28S, 18S, and 5.8S ribosomal DNA (rDNA) gene loci and the rapidly evolving internal transcribed spacers 1 and 2 (ITS1 and ITS2) of 15 sympatric and allopatric fungi associated with colonies of 11 species of leafcutter ants living up to 2,600 km apart in Brazil. We found that the fungal rDNA and ITS sequences from different species of ants were identical (or nearly identical) to each other, whereas 10 GenBank Leucoagaricus species showed higher ITS variation. Our findings suggest that Atta and Acromyrmex leafcutters living in geographic sites that are very distant from each other cultivate a single fungal species made up of closely related lineages of Leucoagaricus gongylophorus. We discuss the strikingly high similarity in the ITS1 and ITS2 regions of the Atta and Acromyrmex symbiotic L. gongylophorus studied by us, in contrast to the lower similarity displayed by their non-symbiotic counterparts. We suggest that the similarity of our L. gongylophorus isolates is an indication of the recent association of the fungus with these ants, and propose that both the intense lateral transmission of fungal material within leafcutter nests and the selection of more adapted fungal strains are involved in the homogenization of the symbiotic fungal stock.

Distinctively variable sequence-based nuclear DNA markers for multilocus phylogeography of the soybean- and rice-infecting fungal pathogen Rhizoctonia solani AG-1 IA

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A influência do esqueleto açúcar-fostato no transporte da molécula de DNA

This dissertation analyses the influence of sugar-phosphate structure in the electronic transport in the double stretch DNA molecule, with the sequence of the base pairs modeled by two types of quasi-periodic sequences: Rudin-Shapiro and Fibonacci. For the sequences, the density of state was calculated and it was compared with the density of state of a piece of human DNA Ch22. After, the electronic transmittance was investigated. In both situations, the Hamiltonians are different. On the analysis of density of state, it was employed the Dyson equation. On the transmittance, the time independent Schrödinger equation was used. In both cases, the tight-binding model was applied. The density of states obtained through Rudin-Shapiro sequence reveal to be similar to the density of state for the Ch22. And for transmittance only until the fifth generation of the Fibonacci sequence was acquired. We have considered long range correlations in both transport mechanism

Conducão eletrônica e propriedades termodinâmicas da molécula de DNA

In this thesis, we study the thermo-electronic properties of the DNA molecule. For this purpose, we used three types of models with the DNA, all assuming a at geometry (2D), each built by a sequence of quasiperiodic (Fibonacci and / or Rudin-Shapiro) and a sequence of natural DNA, part of the human chromosome Ch22. The first two models have two types of components that are the nitrogenous bases (guanine G, cytosine C, adenine A and thymine T) and a cluster sugar-phosphate (SP), while the third has only the nitrogenous bases. In the first model we calculate the density of states using the formalism of Dyson and transmittance for the time independent Schr odinger equation . In the second model we used the renormalizationprocedure for the profile of the transmittance and consequently the I (current) versus V (voltage). In the third model we calculate the density of states formalism by Dean and used the results together with the Fermi-Dirac statistics for the chemical potential and the quantum specific heat. Finally, we compare the physical properties found for the quasi-periodic sequences and those that use a portion of the genomic DNA sequence (Ch22).

Comparison of the sequences of the internal transcribed Spacer regions and PbGP43 genes of Paracoccidioides brasiliensis from patients and armadillos (Dasypus novemcinctus)

Paracoccidioides brasiliensis isolates from 10 nine-banded armadillos (Dasypus novemcinctus) were comparable with 19 clinical isolates by sequence analysis of the PbGP43 gene and ribosomal internal transcribed spacer 1 (ITS1) and ITS2 and by random amplified polymorphic DNA. In this original ITS study, eight isolates differed by one or three sites among five total substitution sites.

EBNA-1 gene sequences in Brazilian and American patients with Hodgkin's disease

We examined the types of Epstein-Barr virus-associated nuclear antigen-1 (EBNA-1) gene carboxy (C)-terminal mutations occurring in Hodgkin's disease (HD) and reactive tissues from two different geographic regions. Previously reported EBNA-1 C-terminal region amino acid sequence variants, based on the amino acid at codon 487, include Prototype (P)-ala, which is found in the B95.8-derived prototype virus, P-thr, Variant (V)-leu, V-val, and V-pro. Using polymerase chain reaction (PCR) to amplify portions of the EBNA-1 gene, followed by DNA sequencing, we found a single EBNA-1 gene sequence variant in each tissue, whether reactive or neoplastic and whether from Brazil or the United States. Variant EBNA-1 gene sequences were more common in both neoplastic and non-neoplastic tissues from different geographic areas than the so-called prototype sequence. In the 17 Brazilian HD cases, 4 cases had P-thr variants and 13 had V-leu variants. In the six reactive tissues from Brazil, one had a P-ala variant, two had P-thr variants, and three had V-leu variants. In the 12 American HD cases, 2 had P-ala variants, 6 had P-thr variants, and 4 had V-leu variants. The 11 American reactive tissues included 2 P ala variants, 5 P-thr variants, and 4 V-leu variants. In both countries, there were similar variant EBNA-1 sequences present in normal tissues and HD cases. Compared with the P ala and P-thr cases, the V-leu cases were more likely to have the 30-bp latent membrane protein 1 (LMP1) gene deletion (P = 0.0075). In addition, cases of HD with the V-leu were statistically associated with a substitution of asparagine for glutamine at codon 322 of the C-terminal portion of the LMP1 gene. Our results suggest that any variation in EBNA-1 gene sequence is caused by a polymorphism present in pre-existing viral strains in the underlying population, and not a mutation occurring during oncogenesis. (C) 1999 by the American Society of Hematology.

Identification of human chromosome 22 transcribed sequences with ORF expressed sequence tags

Transcribed sequences in the human genome can be identified with confidence only by alignment with sequences derived from cDNAs synthesized from naturally occurring mRNAs. We constructed a set of 250,000 cDNAs that represent partial expressed gene sequences and that are biased toward the central coding regions of the resulting transcripts. They are termed ORF expressed sequence tags (ORESTES). The 250,000 ORESTEs were assembled into 81,429 contigs. of these, 1,181 (1.45%) were found to match sequences in chromosome 22 with at least one ORESTES contig for 162 (65.6%) of the 247 known genes, for 67 (44.6%) of the 150 related genes, and for 45 of the 148 (30.4%) EST-predicted genes on this chromosome. Using a set of stringent criteria to validate our sequences, we identified a further 219 previously unannotated transcribed sequences on chromosome 22. of these, 171 were in fact also defined by EST or full length cDNA sequences available in GenBank but not utilized in the initial annotation of the first human chromosome sequence. Thus despite representing less than 15% of all expressed human sequences in the public databases at the time of the present analysis, ORESTEs sequences defined 48 transcribed sequences on chromosome 22 not defined by other sequences. All of the transcribed sequences defined by ORESTEs coincided with DNA regions predicted as encoding exons by GENSCAN.

DNA gains at 8q23.2: a potential early marker in head and neck carcinomas

Gains or amplifications involving chromosome arm 8q are one of the most recurrent chromosomal alterations in head and neck tumors. To characterize previously reported gains, we performed fluorescence in situ hybridization (FISH) using the sequences BAC RP1179E1 and 8-centromere PMJ 128 as probes. Gains and/or amplifications were detected in all 19 cases evaluated by FISH. The FISH analysis, but not G-banding, revealed homogeneously staining region in three cases. We conclude that gains of one or more genes on chromosome arm 8q may be important for the early stages of head and neck carcinomas. (C) 2003 Elsevier B.V. All rights reserved.

Isolation and characterization of microsatellite DNA markers for spinner dolphin (Stenella longirostris)

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Genomic DNA microarray comparison of gene expression patterns in Paracoccidioides brasiliensis mycelia and yeasts in vitro

Paracoccidioides brasiliensis is a thermally dimorphic fungus, and causes the most prevalent systemic mycosis in Latin America. Infection is initiated by inhalation of conidia or mycelial fragments by the host, followed by further differentiation into the yeast form. Information regarding gene expression by either form has rarely been addressed with respect to multiple time points of growth in culture. Here, we report on the construction of a genomic DNA microarray, covering approximately 25% of the genome of the organism, and its utilization in identifying genes and gene expression patterns during growth in vitro. Cloned, amplified inserts from randomly sheared genomic DNA (gDNA) and known control genes were printed onto glass slides to generate a microarray of over 12 000 elements. To examine gene expression, mRNA was extracted and amplified from mycelial or yeast cultures grown in semi-defined medium for 5, 8 and 14 days. Principal components analysis and hierarchical clustering indicated that yeast gene expression profiles differed greatly from those of mycelia, especially at earlier time points, and that mycelial gene expression changed less than gene expression in yeasts over time. Genes upregulated in yeasts were found to encode proteins shown to be involved in methionine/cysteine metabolism, respiratory and metabolic processes (of sugars, amino acids, proteins and lipids), transporters (small peptides, sugars, ions and toxins), regulatory proteins and transcription factors. Mycelial genes involved in processes such as cell division, protein catabolism, nucleotide biosynthesis and toxin and sugar transport showed differential expression. Sequenced clones were compared with Histoplasma capsulatum and Coccidioides posadasii genome sequences to assess potentially common pathways across species, such as sulfur and lipid metabolism, amino acid transporters, transcription factors and genes possibly related to virulence. We also analysed gene expression with time in culture and found that while transposable elements and components of respiratory pathways tended to increase in expression with time, genes encoding ribosomal structural proteins and protein catabolism tended to sharply decrease in expression over time, particularly in yeast. These findings expand our knowledge of the different morphological forms of P. brasiliensis during growth in culture.