Mitochondrial genomes of parasitic flatworms

Complete or near-complete mitochondrial genomes are now available for 11 species or strains of parasitic flatworms belonging to the Trematoda and the Cestoda. The organization of these genomes is not strikingly different from those of other eumetazoans, although one gene (atp8) commonly found in other phyla is absent from flatworms. The gene order in most flatworms has similarities to those seen in higher protostomes such as annelids. However, the gene order has been drastically altered in Schistosoma mansoni, which obscures this possible relationship. Among the sequenced taxa, base composition varies considerably, creating potential difficulties for phylogeny reconstruction. Long non-coding regions are present in all taxa, but these vary in length from only a few hundred to similar to10 000 nucleotides. Among Schistosoma spp., the long non-coding regions are rich in repeats and length variation among individuals is known. Data from mitochondrial genomes are valuable for studies on species identification, phylogenies and biogeography.

Complete mitochondrial genomes confirm the distinctiveness of the horse-dog and sheep-dog strains of Echinococcus granulosus

Unlike other members of the genus, Echinococcus granulosus is known to exhibit considerable levels of variation in biology, physiology and molecular genetics. Indeed, some of the taxa regarded as 'genotypes' within E. granulosus might be sufficiently distinct as to merit specific status. Here, complete mitochondrial genomes are presented of 2 genotypes of E. granulosus (G1-sheep-dog strain: G4-horse-dog strain) and of another taeniid cestode, Taenia crassiceps. These genomes are characterized and compared with those of Echinococcus multilocularis and Hymenolepis diminuta. Genomes of all the species are very similar in structure, length and base-composition. Pairwise comparisons of concatenated protein-coding genes indicate that the G1 and G4 genotypes of E. granulosus are almost as distant from each other as each is from a distinct species, E. multilocularis. Sequences for the variable genes atp6 and nad3 were obtained from additional genotypes of E. granulosus, from E. vogeli and E. oligarthrus. Again, pairwise comparisons showed the distinctiveness of the G1 and G4 genotypes. Phylogenetic analyses of concatenated atp6, nad1 (partial) and cox1 (partial) genes from E. multilocularis, E. vogeli, E. oligarthrus, 5 genotypes of E. granulosus, and using T. crassiceps as an outgroup, yielded the same results. We conclude that the sheep-dog and horse-dog strains of E. granulosus should be regarded as distinct at the specific level.

New Strategies on Molecular Biology Applied to Microbial Systematics

Systematics is the study of diversity of the organisms and their relationships comprising classification, nomenclature and identification. The term classification or taxonomy means the arrangement of the organisms in groups (rate) and the nomenclature is the attribution of correct international scientific names to organisms and identification is the inclusion of unknown strains in groups derived from classification. Therefore, classification for a stable nomenclature and a perfect identification are required previously. The beginning of the new bacterial systematics era can be remembered by the introduction and application of new taxonomic concepts and techniques, from the 50’s and 60’s. Important progress were achieved using numerical taxonomy and molecular taxonomy. Molecular taxonomy, brought into effect after the emergence of the Molecular Biology resources, provided knowledge that comprises systematics of bacteria, in which occurs great evolutionary interest, or where is observed the necessity of eliminating any environmental interference. When you study the composition and disposition of nucleotides in certain portions of the genetic material, you study searching their genome, much less susceptible to environmental alterations than proteins, codified based on it. In the molecular taxonomy, you can research both DNA and RNA, and the main techniques that have been used in the systematics comprise the build of restriction maps, DNA-DNA hybridization, DNA-RNA hybridization, sequencing of DNA sequencing of sub-units 16S and 23S of rRNA, RAPD, RFLP, PFGE etc. Techniques such as base sequencing, though they are extremely sensible and greatly precise, are relatively onerous and impracticable to the great majority of the bacterial taxonomy laboratories. Several specialized techniques have been applied to taxonomic studies of microorganisms. In the last years, these have included preliminary electrophoretic analysis of soluble proteins and isoenzymes, and subsequently determination of deoxyribonucleic acid base composition and assessment of base sequence homology by means of DNA-RNA hybrid experiments beside others. These various techniques, as expected, have generally indicated a lack of taxonomic information in microbial systematics. There are numberless techniques and methodologies that make bacteria identification and classification study possible, part of them described here, allowing establish different degrees of subspecific and interspecific similarity through phenetic-genetic polymorphism analysis. However, was pointed out the necessity of using more than one technique for better establish similarity degrees within microorganisms. Obtaining data resulting from application of a sole technique isolatedly may not provide significant information from Bacterial Systematics viewpoint

Characterization of molecular damage induced by UV photons and carbon ions on biomimetic heterostructures

The study of the effect of radiation on living tissues is a rather complex task to address mainly because they are made of a set of complex functional biological structures and interfaces. Particularly if one is looking for where damage is taking place in a first stage and what are the underlying reaction mechanisms. In this work a new approach is addressed to study the effect of radiation by making use of well identified molecular hetero-structures samples which mimic the biological environment. These were obtained by assembling onto a solid support deoxyribonucleic acid (DNA) and phospholipids together with a soft water-containing polyelectrolyte precursor in layered structures and by producing lipid layers at liquid/air interface with DNA as subphase. The effects of both ultraviolet (UV) radiation and carbon ions beams were systematically investigated in these heterostructures, namely damage on DNA by means vacuum ultraviolet (VUV), infrared (IR), X-Ray Photoelectron (XPS) and impedance spectroscopy. Experimental results revealed that UV affects furanose, PO2-, thymines, cytosines and adenines groups. The XPS spectrometry carried out on the samples allowed validate the VUV and IR results and to conclude that ionized phosphate groups, surrounded by the sodium counterions, congregate hydration water molecules which play a role of UV protection. The ac electrical conductivity measurements revealed that the DNA electrical conduction is arising from DNA chain electron hopping between base-pairs and phosphate groups, with the hopping distance equal to the distance between DNA base-pairs and is strongly dependent on UV radiation exposure, due loss of phosphate groups. Characterization of DNA samples exposed to a 4 keV C3+ ions beam revealed also carbon-oxygen bonds break, phosphate groups damage and formation of new species. Results from radiation induced damage carried out on biomimetic heterostructures having different compositions revealed that damage is dependent on sample composition, with respect to functional targeted groups and extent of damage. Conversely, LbL films of 1,2-dipalmitoyl-sn-Glycero-3-[Phospho-rac-(1-glycerol)] (Sodium Salt) (DPPG) liposomes, alternated with poly(allylamine hydrochloride) (PAH) revealed to be unaffected, even by prolonged UV irradiation exposure, in the absence of water molecules. However, DPPG molecules were damaged by the UV radiation in presence of water with cleavage of C-O, C=O and –PO2- bonds. Finally, the study of DNA interaction with the ionic lipids at liquid/air interfaces revealed that electrical charge of the lipid influences the interaction of phospholipid with DNA. In the presence of DNA in the subphase, the effects from UV irrladiation were seen to be smaller, which means that ionic products from biomolecules degradation stabilize the intact DPPG molecules. This mechanism may explain why UV irradiation does not cause immediate cell collapse, thus providing time for the cellular machinery to repair elements damaged by UV.

Burkholderia sartisoli sp. nov., isolated from a polycyclic aromatic hydrocarbon-contaminated soil.

A Gram-negative, rod-shaped, aerobic bacterium, designated strain RP007(T), was isolated from a polycyclic aromatic hydrocarbon-contaminated soil in New Zealand. Two additional strains were recovered from a compost heap in Belgium (LMG 18808) and from the rhizosphere of maize in the Netherlands (LMG 24204). The three strains had virtually identical 16S rRNA gene sequences and whole-cell protein profiles, and they were identified as members of the genus Burkholderia, with Burkholderia phenazinium as their closest relative. Strain RP007(T) had a DNA G+C content of 63.5 mol% and could be distinguished from B. phenazinium based on a range of biochemical characteristics. Strain RP007(T) showed levels of DNA-DNA relatedness towards the type strain of B. phenazinium and those of other recognized Burkholderia species of less than 30 %. The results of 16S rRNA gene sequence analysis, DNA-DNA hybridization experiments and physiological and biochemical tests allowed the differentiation of strain RP007(T) from all recognized species of the genus Burkholderia. Strains RP007(T), LMG 18808 and LMG 24204 are therefore considered to represent a single novel species of the genus Burkholderia, for which the name Burkholderia sartisoli sp. nov. is proposed. The type strain is RP007(T) (=LMG 24000(T) =CCUG 53604(T) =ICMP 13529(T)).

Genometric analyses of the organization of circular chromosomes: a universal pressure determines the direction of ribosomal RNA genes transcription relative to chromosome replication.

Selective pressures related to gene function and chromosomal architecture are acting on genome sequences and can be revealed, for instance, by appropriate genometric methods. Cumulative nucleotide skew analyses, i.e., GC, TA, and ORF orientation skews, predict the location of the origin of DNA replication for 88 out of 100 completely sequenced bacterial chromosomes. These methods appear fully reliable for proteobacteria, Gram-positives, and spirochetes as well as for euryarchaeotes. Based on this genome architecture information, coorientation analyses reveal that in prokaryotes, ribosomal RNA (rRNA) genes encoding the small and large ribosomal subunits are all transcribed in the same direction as DNA replication; that is, they are located along the leading strand. This result offers a simple and reliable method for circumscribing the region containing the origin of the DNA replication and reveals a strong selective pressure acting on the orientation of rRNA genes similar to the weaker one acting on the orientation of ORFs. Rate of coorientation of transfer RNA (tRNA) genes with DNA replication appears to be taxon-specific. Analyzing nucleotide biases such as GC and TA skews of genes and plotting one against the other reveals a taxonomic clusterization of species. All ribosomal RNA genes are enriched in Gs and depleted in Cs, the only so far known exception being the rRNA genes of deuterostomian mitochondria. However, this exception can be explained by the fact that in the chromosome of the human mitochondrion, the model of the deuterostomian organelle genome, DNA replication, and rRNA transcription proceed in opposite directions. A general rule is deduced from prokaryotic and mitochondrial genomes: ribosomal RNA genes that are transcribed in the same direction as the DNA replication are enriched in Gs, and those transcribed in the opposite direction are depleted in Gs.

Vitellogenin B2 gene in Xenopus laevis: isolation, in vitro transcription and relation to other vitellogenin genes.

The isolation of the four Xenopus laevis vitellogenin genes has been completed by the purification from a DNA library of the B2 gene together with its flanking sequences. The overlapping DNA fragments analyzed cover 34 kilobases. The B2 gene which has a length of 17.5 kilobases was characterized by heteroduplex and R-loop mapping in the electron microscope and by in vitro transcription in a HeLa whole-cell extract. Its structural organization is compared with that of the closely related B1 gene. The mRNA-coding sequence of about 6 kilobases is interrupted 34 times in the B1 gene and 33 times in the B2 gene. Sequence homology between the two genes was not only found in exons. In addition, 54% of the intron sequences as well as 63% and 48.5% respectively of the 5' and 3' flanking sequences, show enough homology to form stable duplexes. These findings are compared with earlier results obtained with the two other closely related members of the vitellogenin gene family, the A1 and the A2 genes.

Monophyly of beta-tubulin and H+-ATPase gene variants in Glomus intraradices: consequences for molecular evolutionary studies of AM fungal genes.

Arbuscular mycorrhizal fungi (AMF) are an ecologically important group of fungi. Previous studies showed the presence of divergent copies of beta-tubulin and V-type vacuolar H+-ATPase genes in AMF genomes and suggested horizontal gene transfer from host plants or mycoparasites to AMF. We sequenced these genes from DNA isolated from an in vitro cultured isolate of Glomus intraradices that was free of any obvious contaminants. We found two highly variable beta-tubulin sequences and variable H+-ATPase sequences. Despite this high variation, comparison of the sequences with those in gene banks supported a glomeromycotan origin of G. intraradices beta-tubulin and H+-ATPase sequences. Thus, our results are in sharp contrast with the previously reported polyphyletic origin of those genes. We present evidence that some highly divergent sequences of beta-tubulin and H+-ATPase deposited in the databases are likely to be contaminants. We therefore reject the prediction of horizontal transfer to AMF genomes. High differences in GC content between glomeromycotan sequences and sequences grouping in other lineages are shown and we suggest they can be used as an indicator to detect such contaminants. H+-ATPase phylogeny gave unexpected results and failed to resolve fungi as a natural group. beta-Tubulin phylogeny supported Glomeromeromycota as sister group of the Chytridiomycota. Contrasts between our results and trees previously generated using rDNA sequences are discussed.

M-protein and other intrinsic virulence factors of Streptococcus pyogenes are encoded on an ancient pathogenicity island.

BACKGROUND: The increasing number of completely sequenced bacterial genomes allows comparing their architecture and genetic makeup. Such new information highlights the crucial role of lateral genetic exchanges in bacterial evolution and speciation. RESULTS: Here we analyzed the twelve sequenced genomes of Streptococcus pyogenes by a naïve approach that examines the preferential nucleotide usage along the chromosome, namely the usage of G versus C (GC-skew) and T versus A (TA-skew). The cumulative GC-skew plot presented an inverted V-shape composed of two symmetrical linear segments, where the minimum and maximum corresponded to the origin and terminus of DNA replication. In contrast, the cumulative TA-skew presented a V-shape, which segments were interrupted by several steep slopes regions (SSRs), indicative of a different nucleotide composition bias. Each S. pyogenes genome contained up to nine individual SSRs, encompassing all described strain-specific prophages. In addition, each genome contained a similar unique non-phage SSR, the core of which consisted of 31 highly homologous genes. This core includes the M-protein, other mga-related factors and other virulence genes, totaling ten intrinsic virulence genes. In addition to a high content in virulence-related genes and to a peculiar nucleotide bias, this SSR, which is 47 kb-long in a M1GAS strain, harbors direct repeats and a tRNA gene, suggesting a mobile element. Moreover, its complete absence in a M-protein negative group A Streptococcus natural isolate demonstrates that it could be spontaneously lost, but in vitro deletion experiments indicates that its excision occurred at very low rate. The stability of this SSR, combined to its presence in all sequenced S. pyogenes sequenced genome, suggests that it results from an ancient acquisition. CONCLUSION: Thus, this non-phagic SSR is compatible with a pathogenicity island, acquired before S. pyogenes speciation. Its potential excision might bear relevance for vaccine development, because vaccines targeting M-protein might select for M-protein-negative variants that still carry other virulence determinants.

Biosynthesis of pyochelin and dihydroaeruginoic acid requires the iron-regulated pchDCBA operon in Pseudomonas aeruginosa.

The high-affinity siderophore salicylate is an intermediate in the biosynthetic pathway of pyochelin, another siderophore and chelator of transition metal ions, in Pseudomonas aeruginosa. The 2.5-kb region upstream of the salicylate biosynthetic genes pchBA was sequenced and found to contain two additional, contiguous genes, pchD and pchC, having the same orientation. The deduced amino acid sequence of the 60-kDa PchD protein was similar to those of the EntE protein (2,3-dihydroxybenzoate-AMP ligase) of Escherichia coli and other adenylate-forming enzymes, suggesting that salicylate might be adenylated at the carboxyl group by PchD. The 28-kDa PchC protein showed similarities to thioesterases of prokaryotic and eukaryotic origin and might participate in the release of the product(s) formed from activated salicylate. One potential product, dihydroaeruginoate (Dha), was identified in culture supernatants of iron-limited P. aeruginosa cells. The antifungal antibiotic Dha is thought to arise from the reaction of salicylate with cysteine, followed by cyclization of cysteine. Inactivation of the chromosomal pchD gene by insertion of the transcription and translation stop element omega Sm/Sp abolished the production of Dha and pyochelin, implying that PchD-mediated activation of salicylate may be a common first step in the synthesis of both metabolites. Furthermore, the pchD::omega Sm/Sp mutation had a strong polar effect on the expression of the pchBA genes, i.e., on salicylate synthesis, indicating that the pchDCBA genes constitute a transcriptional unit. A full-length pchDCBA transcript of ca. 4.4 kb could be detected in iron-deprived, growing cells of P. aeruginosa. Transcription of pchD started at tandemly arranged promoters, which overlapped with two Fur boxes (binding sites for the ferric uptake regulator) and the promoter of the divergently transcribed pchR gene encoding an activator of pyochelin biosynthesis. This promoter arrangement allows tight iron-mediated repression of the pchDCBA operon.

Transcriptional regulatory patterns of the myelin basic protein and malic enzyme genes by the thyroid hormone receptors alpha1 and beta1.

While there is evidence that the two ubiquitously expressed thyroid hormone (T3) receptors, TRalpha1 and TRbeta1, have distinct functional specificities, the mechanism by which they discriminate potential target genes remains largely unexplained. In this study, we demonstrate that the thyroid hormone response elements (TRE) from the malic enzyme and myelin basic protein genes (METRE and MBPTRE) respectively, are not functionally equivalent. The METRE, which is a direct repeat motif with a 4-base pair gap between the two half-site hexamers binds thyroid hormone receptor as a heterodimer with 9-cis-retinoic acid receptor (RXR) and mediates a high T3-dependent activation in response to TRalpha1 or TRbeta1 in NIH3T3 cells. In contrast, the MBPTRE, which consists of an inverted palindrome formed by two hexamers spaced by 6 base pairs, confers an efficient transactivation by TRbeta1 but a poor transactivation by TRalpha1. While both receptors form heterodimers with RXR on MBPTRE, the poor transactivation by TRalpha1 correlates also with its ability to bind efficiently as a monomer. This monomer, which is only observed with TRalpha1 bound to MBPTRE, interacts neither with N-CoR nor with SRC-1, explaining its functional inefficacy. However, in Xenopus oocytes, in which RXR proteins are not detectable, the transactivation mediated by TRalpha1 and TRbeta1 is equivalent and independent of a RXR supply, raising the question of the identity of the thyroid hormone receptor partner in these cells. Thus, in mammalian cells, the binding characteristics of TRalpha1 to MBPTRE (i.e. high monomer binding efficiency and low transactivation activity) might explain the particular pattern of T3 responsiveness of MBP gene expression during central nervous system development.

Phylogenomics Controlling for Base Compositional Bias Reveals a Single Origin of Eusociality in Corbiculate Bees.

As increasingly large molecular data sets are collected for phylogenomics, the conflicting phylogenetic signal among gene trees poses challenges to resolve some difficult nodes of the Tree of Life. Among these nodes, the phylogenetic position of the honey bees (Apini) within the corbiculate bee group remains controversial, despite its considerable importance for understanding the emergence and maintenance of eusociality. Here, we show that this controversy stems in part from pervasive phylogenetic conflicts among GC-rich gene trees. GC-rich genes typically have a high nucleotidic heterogeneity among species, which can induce topological conflicts among gene trees. When retaining only the most GC-homogeneous genes or using a nonhomogeneous model of sequence evolution, our analyses reveal a monophyletic group of the three lineages with a eusocial lifestyle (honey bees, bumble bees, and stingless bees). These phylogenetic relationships strongly suggest a single origin of eusociality in the corbiculate bees, with no reversal to solitary living in this group. To accurately reconstruct other important evolutionary steps across the Tree of Life, we suggest removing GC-rich and GC-heterogeneous genes from large phylogenomic data sets. Interpreted as a consequence of genome-wide variations in recombination rates, this GC effect can affect all taxa featuring GC-biased gene conversion, which is common in eukaryotes.

Intracellular antioxidant and DNA repair enzymes as correlates of stress resistance and longevity in vertebrates

In animals, both stress resistance and longevity appear to be influenced by the insulin/insulin-like growth factor-l signaling (lIS) pathway, the basic organization of which is highly conserved from invertebrates to vertebrates. Reduced lIS or genetic disruption of the lIS pathway leads to the activation of forkhead box transcription factors, which is thought to upregulate the expression of genes involved in enhancing stress resistance, including perhaps key antioxidant enzymes as well as DNA repair enzymes. Enhanced antioxidant and DNA repair capacities may underlie the enhanced cellular stress resistance observed in long-lived animals, however little data is available that directly supports this idea. I used three. experimental approaches to test the association of intracellular antioxidant and DNA base excision repair (BER) capacities with stress resistance and longevity: (1) a comparison of multiple vertebrate endotherm species of varying body masses and longevities; (2) a comparison of long-lived Snell dwarf mice and their normallittermates; and (3) a comparison of hypometabolic animals undergoing hibernation or estivation with their active counterparts. The activities of the five major intracellular antioxidant enzymes as well as the two rate-limiting enzymes in the BER pathway, apurininc/apyrimidinic (AP) endonuclease and polymerase ~, were measured. These measurements were performed in one or more of the following: (1) cultured dermal fibroblasts; (2) brain tissue; (3) heart tissue; (4) liver tissue. My results indicate that antioxidant enzymes are not universally upregulated in association with enhanced stress resistance and longevity. I also did not find that BER enzyme activity was positively correlated with longevity, in an inter-species context, though there was evidence for enhanced BER in long-lived Snell dwarf mice. Thus, while there were instances in which enhanced antioxidant and BER enzyme activities were associated with increased stress resistance and/or longevity, this was not universally the case, indicating that other mechanisms must be involved. These results suggest the need to re-examine existing 'oxidative stress' hypotheses of longevity and probe further into the molecular physiology of longevity to discover its mechanistic basis.

Improving Short DNA Sequence Alignment with Parallel Computing

Variations in different types of genomes have been found to be responsible for a large degree of physical diversity such as appearance and susceptibility to disease. Identification of genomic variations is difficult and can be facilitated through computational analysis of DNA sequences. Newly available technologies are able to sequence billions of DNA base pairs relatively quickly. These sequences can be used to identify variations within their specific genome but must be mapped to a reference sequence first. In order to align these sequences to a reference sequence, we require mapping algorithms that make use of approximate string matching and string indexing methods. To date, few mapping algorithms have been tailored to handle the massive amounts of output generated by newly available sequencing technologies. In otrder to handle this large amount of data, we modified the popular mapping software BWA to run in parallel using OpenMPI. Parallel BWA matches the efficiency of multithreaded BWA functions while providing efficient parallelism for BWA functions that do not currently support multithreading. Parallel BWA shows significant wall time speedup in comparison to multithreaded BWA on high-performance computing clusters, and will thus facilitate the analysis of genome sequencing data.