Genetic relatedness within the genus Campylobacter inferred from rpoB sequences

The genus Campylobacter comprises 17 species, some of which are important animal and human pathogens. To gain more insight into the genetic relatedness of this genus and to improve the molecular tools available for diagnosis, a universal sequencing approach was established for the gene encoding the beta-subunit of RNA polymerase (rpoB) for the genus Campylobacter. A total of 59 strains, including the type strains of currently recognized species as well as field isolates, were investigated in the study. A primer set specific for Campylobacter species enabled straightforward amplification and sequencing of a 530 bp fragment of the rpoB gene. The 16S rRNA gene sequences of all of the strains were determined in parallel. A good congruence was obtained between 16S rRNA and rpoB gene sequence-based trees within the genus Campylobacter. The branching of the rpoB tree was similar to that of the 16S rRNA gene tree, even though a few discrepancies were observed for certain species. The resolution of the rpoB gene within the genus Campylobacter was generally much higher than that of the 16S rRNA gene sequence, resulting in a clear separation of most species and even some subspecies. The universally applicable amplification and sequencing approach for partial rpoB gene sequence determination provides a powerful tool for DNA sequence-based discrimination of Campylobacter species.

GENE EXPRESSION IN CHRONIC MYELOGENOUS LEUKEMIA

The goal of the present work was to identify and characterize gene sequences that are preferentially expressed in CML in an effort to better understand the molecular basis of the disease. As high abundance mRNAs generally encode proteins that are phenotypically characteristic of cells, positive-negative screening of a CML cDNA library was used to identify cDNA clones containing sequences preferentially transcribed in CML. One cDNA sequence that fulfilled this criterion, C-A3, has been characterized in some detail. It represents a small mRNA ((TURN)496 nucleotides) that is highly abundant ((TURN)2% of the poly(A('+))RNA) in cells from the chronic phase of CML. In situ hybridization to whole cells indicates the principal leukocytes that express C-A3 sequences are eosinophils, basophils and immature myelocytes. Surprisingly, CML patients with high numbers of myeloblasts do not have an abundance of C-A3 transcripts, although transcript levels remain elevated in patients with lymphoblasts. In AML, high transcript levels are only found sporadically and occasionally different sized transcripts can be detected. Sequences from the 3' end of the C-A3 message are present in 2-5 copies per haploid genome. The 3' end of C-A3 localizes to bands 8q21.1 and 8q23 by in situ chromosomal hybridization. This is a region that is often involved in hematopoietic malignancies. Restriction digests of human genomic DNA show a correlation between the presence of a 2.3 kb Hind III fragment and certain types of leukemia. All of the leukemic DNAs tested had this fragment. In comparison, only one of five normal DNAs had a band this size. Analysis of the nucleotide sequence indicates that C-A3 probably encodes a small, hydrophobic peptide which may be part of a larger protein. ^

Structures of four human pseudogenes for U7 small nuclear RNA.

Four U7 RNA-related sequences were isolated from a human genomic DNA library. None of the sequences completely match the published human U7 RNA sequence and all of them contain features typical of reverse-transcribed pseudogenes.

Metadata und statistic analysis of archael and bacterial sequences originating from sediments of the Håkon Mosby mud volcano (Z1-Z3, all habitats)

DNA extraction was carried out as described on the MICROBIS project pages (http://icomm.mbl.edu/microbis ) using a commercially available extraction kit. We amplified the hypervariable regions V4-V6 of archaeal and bacterial 16S rRNA genes using PCR and several sets of forward and reverse primers (http://vamps.mbl.edu/resources/primers.php). Massively parallel tag sequencing of the PCR products was carried out on a 454 Life Sciences GS FLX sequencer at Marine Biological Laboratory, Woods Hole, MA, following the same experimental conditions for all samples. Sequence reads were submitted to a rigorous quality control procedure based on mothur v30 (doi:10.1128/AEM.01541-09) including denoising of the flow grams using an algorithm based on PyroNoise (doi:10.1038/nmeth.1361), removal of PCR errors and a chimera check using uchime (doi:10.1093/bioinformatics/btr381). The reads were taxonomically assigned according to the SILVA taxonomy (SSURef v119, 07-2014; doi:10.1093/nar/gks1219) implemented in mothur and clustered at 98% ribosomal RNA gene V4-V6 sequence identity. V4-V6 amplicon sequence abundance tables were standardized to account for unequal sampling effort using 1000 (Archaea) and 2300 (Bacteria) randomly chosen sequences without replacement using mothur and then used to calculate inverse Simpson diversity indices and Chao1 richness (doi:10.2307/4615964). Bray-Curtis dissimilarities (doi:10.2307/1942268) between all samples were calculated and used for 2-dimensional non metric multidimensional scaling (NMDS) ordinations with 20 random starts (doi:10.1007/BF02289694). Stress values below 0.2 indicated that the multidimensional dataset was well represented by the 2D ordination. NMDS ordinations were compared and tested using Procrustes correlation analysis (doi:10.1007/BF02291478). All analyses were carried out with the R statistical environment and the packages vegan (available at: http://cran.r-project.org/package=vegan), labdsv (available at: http://cran.r-project.org/package=labdsv), as well as with custom R scripts. Operational taxonomic units at 98% sequence identity (OTU0.03) that occurred only once in the whole dataset were termed absolute single sequence OTUs (SSOabs; doi:10.1038/ismej.2011.132). OTU0.03 sequences that occurred only once in at least one sample, but may occur more often in other samples were termed relative single sequence OTUs (SSOrel). SSOrel are particularly interesting for community ecology, since they comprise rare organisms that might become abundant when conditions change.16S rRNA amplicons and metagenomic reads have been stored in the sequence read archive under SRA project accession number SRP042162.

The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species

While genome sequencing projects are advancing rapidly, EST sequencing and analysis remains a primary research tool for the identification and categorization of gene sequences in a wide variety of species and an important resource for annotation of genomic sequence. The TIGR Gene Indices (http://www.tigr.org/tdb/tgi.shtml) are a collection of species-specific databases that use a highly refined protocol to analyze EST sequences in an attempt to identify the genes represented by that data and to provide additional information regarding those genes. Gene Indices are constructed by first clustering, then assembling EST and annotated gene sequences from GenBank for the targeted species. This process produces a set of unique, high-fidelity virtual transcripts, or Tentative Consensus (TC) sequences. The TC sequences can be used to provide putative genes with functional annotation, to link the transcripts to mapping and genomic sequence data, to provide links between orthologous and paralogous genes and as a resource for comparative sequence analysis.

rrndb: the Ribosomal RNA Operon Copy Number Database

The Ribosomal RNA Operon Copy Number Database (rrndb) is an Internet-accessible database containing annotated information on rRNA operon copy number among prokaryotes. Gene redundancy is uncommon in prokaryotic genomes, yet the rRNA genes can vary from one to as many as 15 copies. Despite the widespread use of 16S rRNA gene sequences for identification of prokaryotes, information on the number and sequence of individual rRNA genes in a genome is not readily accessible. In an attempt to understand the evolutionary implications of rRNA operon redundancy, we have created a phylogenetically arranged report on rRNA gene copy number for a diverse collection of prokaryotic microorganisms. Each entry (organism) in the rrndb contains detailed information linked directly to external websites including the Ribosomal Database Project, GenBank, PubMed and several culture collections. Data contained in the rrndb will be valuable to researchers investigating microbial ecology and evolution using 16S rRNA gene sequences. The rrndb web site is directly accessible on the WWW at http://rrndb.cme.msu.edu.

ACTIVITY: a database on DNA/RNA sites activity adapted to apply sequence-activity relationships from one system to another

ACTIVITY is a database on DNA/RNA site sequences with known activity magnitudes, measurement systems, sequence-activity relationships under fixed experimental conditions and procedures to adapt these relationships from one measurement system to another. This database deposits information on DNA/RNA affinities to proteins and cell nuclear extracts, cutting efficiencies, gene transcription activity, mRNA translation efficiencies, mutability and other biological activities of natural sites occurring within promoters, mRNA leaders, and other regulatory regions in pro- and eukaryotic genomes, their mutant forms and synthetic analogues. Since activity magnitudes are heavily system-dependent, the current version of ACTIVITY is supplemented by three novel sub-databases: (i) SYSTEM, measurement systems; (ii) KNOWLEDGE, sequence-activity relationships under fixed experimental conditions; and (iii) CROSS_TEST, procedures adapting a relationship from one measurement system to another. These databases are useful in molecular biology, pharmacogenetics, metabolic engineering, drug design and biotechnology. The databases can be queried using SRS and are available through the Web, http://wwwmgs.bionet.nsc.ru/systems/Activity/.

Transfer RNA editing in land snail mitochondria.

Some mitochondrial tRNA genes of land snails show mismatches in the acceptor stems predicted from their gene sequences. The majority of these mismatches fall in regions where the tRNA genes overlap with adjacent downstream genes. We have synthesized cDNA from four circularized tRNAs and determined the sequences of the 5' and 3' parts of their acceptor stems. Three of the four tRNAs differ from their corresponding genes at a total of 13 positions, which all fall in the 3' part of the acceptor stems as well as the discriminator bases. The editing events detected involve changes from cytidine, thymidine, and guanosine to adenosine residues, which generally restore base-pairing in the stems. However, in one case an A-A mismatch is created from an A-C mismatch. It is suggested that this form of RNA editing may involve polyadenylylation of the maturing tRNAs as an intermediate.

Specific binding of RNA polymerase II to the human immunodeficiency virus trans-activating region RNA is regulated by cellular cofactors and Tat.

The regulation of human immunodeficiency virus type 1 (HIV-1) gene expression in response to Tat is dependent on an element downstream of the HIV-1 transcriptional initiation site designated the trans-activating region (TAR). TAR forms a stable stem-loop RNA structure in which a 3-nt bulge structure and a 6-nt loop structure are important for Tat activation. In the absence of Tat, the HIV-1 promoter generates so-called short or nonprocessive transcripts terminating at +60, while in the presence of Tat the synthesis of these short transcripts is markedly decreased and transcripts that extend through the 9.0-kb HIV-1 genome are synthesized. Tat effects on transcriptional elongation are likely due to alterations in the elongation properties of RNA polymerase II. In this study we demonstrated that a set of cellular cofactors that modulate the binding of the cellular protein TRP-185 to the TAR RNA loop sequences also functioned to markedly stimulate the specific binding of hypophosphorylated (IIa) and hyperphosphorylated (IIo) RNA polymerase II to TAR RNA. The concentrations of RNA polymerase II required for this interaction with TAR RNA were similar to those required to initiate in vitro transcription from the HIV-1 long terminal repeat. RNA gel retardation analysis with wild-type and mutant TAR RNAs indicated that the TAR RNA loop and bulge sequences were critical for the binding of RNA polymerase II. The addition of wild-type but not mutant Tat protein to gel retardation analysis with TAR RNA and RNA polymerase II resulted in the loss of binding of RNA polymerase II binding to TAR RNA. These results suggest that Tat may function to alter RNA polymerase II, which is paused due to its binding to HIV-1 TAR RNA with resultant stimulation of its transcriptional elongation properties.

Whole genome analysis reveals a high incidence of non-optimal codons in secretory signal sequences of Escherichia coli

Translational pausing may occur due to a number of mechanisms, including the presence of non-optimal codons, and it is thought to play a role in the folding of specific polypeptide domains during translation and in the facilitation of signal peptide recognition during see-dependent protein targeting. In this whole genome analysis of Escherichia coli we have found that non-optimal codons in the signal peptide-encoding sequences of secretory genes are overrepresented relative to the mature portions of these genes; this is in addition to their overrepresentation in the 5'-regions of genes encoding non-secretory proteins. We also find increased non-optimal codon usage at the 3' ends of most E. coli genes, in both non-secretory and secretory sequences. Whereas presumptive translational pausing at the 5' and 3' ends of E. coli messenger RNAs may clearly have a general role in translation, we suggest that it also has a specific role in sec-dependent protein export, possibly in facilitating signal peptide recognition. This finding may have important implications for our understanding of how the majority of non-cytoplasmic proteins are targeted, a process that is essential to all biological cells. (C) 2004 Elsevier Inc. All rights reserved.

Reference gene selection and RNA preservation protocol in the cat flea, Ctenocephalides felis, for gene expression studies

This work was supported by a Knowledge Transfer Network BBSRC Industrial Case (#414 BB/L502467/1) studentship in association Zoetis Inc.

Structural Basis of RNA Recognition by Mycobacterium tuberculosis MazF Homologues

The MazEF toxin-antitoxin (TA) system consists of the antitoxin MazE and the toxin MazF. MazF is a sequence-specific endoribonuclease that upon activation causes cellular growth arrest and increass the level of persisters. Moreover, MazF-induced cells are in a quasi-dormant state that cells remain metabolically active while stop dividing. The quasi-dormancy is similar to the nonreplicating state of M. tuberculosis during latent tuberculosis, thus suggesting the role of mazEF in M. tuberculosis dormancy and persistence. M. tuberculosis has nine mazEF TA modules, each with different RNA cleavage specificities and implicated in selective gene expression during stress conditions. To date only the Bacillus subtilis MazF-RNA complex structure has been determined. As M. tuberculosis MazF homologues recognize distinct RNA sequences, their molecular mechanisms of substrate specificity remain unclear. By taking advantage of X-ray crystallography, we have determined structures of two M. tuberculosis MazF-RNA complexes, MazF-mt1 (Rv2801c) and MazF-mt3 (Rv1991c) in complex with an uncleavable RNA substrate. These structures have provided the molecular basis of sequence-specific RNA recognition and cleavage by MazF toxins.

Both MazF-mt1-RNA and MazF-mt3-RNA complexes showed similar structural organization with one molecule of RNA bound to a MazF-mt1 or MazF-mt3 dimer and occupying the same pocket within the MazF dimer interface. Similar to B. subtilis MazF-RNA complex, MazF-mt1 and MazF-mt3 displayed a conserved active site architecture, where two highly conserved residues, Arg and Thr, form hydrogen bonds with the scissile phosphate group in the cleavage site of the bound RNA. The MazF-mt1-RNA complex also showed specific interactions with its three-base RNA recognition element. Compared with the B. subtilis MazF-RNA complex, our structures showed that residues involved in sequence-specific recognition of target RNA vary between the MazF homologues, therefore explaining the molecular basis for their different RNA recognition sequences. In addition, local conformational changes of the loops in the RNA binding site of MazF-mt1 appear to play a role in MazF targeting different RNA lengths and sequences. In contrast, the MazF-mt3-RNA complex is in a non-optimal RNA binding state with a symmetry-related MazF-mt3 molecule found to make interactions with the bound RNA in the crystal. The crystal-packing interactions were further examined by isothermal titration calorimetry (ITC) studies on selected MazF-mt3 mutants. Our attempts to utilize a MazF-mt3 mutant bearing mutations involved in crystal contacts all crystallized with few nucleotides, which are still found to interact with a symmetry mate. However, these different crystal forms revealed the conformational flexibility of loops in the RNA binding interface of MazF-mt3, suggesting their role in RNA binding and recognition, which will require further studies on additional MazF-mt3-RNA complex interactions.

In conclusion, the structures of the MazF-mt1-RNA and MazF-mt3-RNA complexes provide the first structural information on any M. tuberculosis MazF homologues. Supplemented with structure-guided mutational studies on MazF toxicity in vivo, this study has addressed the structural basis of different RNA cleavage specificities among MazF homologues. Our work will guide future studies on the function of other M. tuberculosis MazF and MazE-MazF homologues, and will help delineate their physiological roles in M. tuberculosis stress responses and pathogenesis.

On the role of RNA silencing in the pathogenicity and evolution of viroids and viral satellites

Viroids and most viral satellites have small, noncoding, and highly structured RNA genomes. How they cause disease symptoms without encoding proteins and why they have characteristic secondary structures are two longstanding questions. Recent studies have shown that both viroids and satellites are capable of inducing RNA silencing, suggesting a possible role of this mechanism in the pathology and evolution of these subviral RNAs. Here we show that preventing RNA silencing in tobacco, using a silencing suppressor, greatly reduces the symptoms caused by the Y satellite of cucumber mosaic virus. Furthermore, tomato plants expressing hairpin RNA, derived from potato spindle tuber viroid, developed symptoms similar to those of potato spindle tuber viroid infection. These results provide evidence suggesting that viroids and satellites cause disease symptoms by directing RNA silencing against physiologically important host genes. We also show that viroid and satellite RNAs are significantly resistant to RNA silencing-mediated degradation, suggesting that RNA silencing is an important selection pressure shaping the evolution of the secondary structures of these pathogens.

Characterisation of a glutathione reductase gene and its genetic locus from pea (Pisum sativum L.)

A cDNA encoding the chloroplast/mitochondrial form of glutathione reductase (GR:EC 1,6,4,2) from pea (Pisum sativum L.) was used to map a single GR locus, named GORI. In two domesticated genotypes of pea (cv, Birte and JI 399) it is likely that the GORI locus contains a single gene. However, in a semi-domesticated land race of pea sequences were detected but closely related sets of GR gene sequences were in JI 281 represent either a second intact gene or a partial or pseudogene copy. A GR gene was cloned from ev. Birte, sequenced and its structure analysed. No features of the transcription or structure of the gene suggested a mechanism for generating any more than one form of . From these data plus previously published biochemical evidence was suggested a second, distinct gene encoding for the cytosolic form of GR should be present in peas. The GORI-encoded GR mRNA can be detected in all main organs of the plant and no alternative spliced species was present which could perhaps account for the generation of multiple isoforms of GR. The mismatch between the number of charge-separable isoforms in pea and the proposed number suggests that different GR isoforms arise by some form of post-transnational modification.

Archaea in the Mycorrhizosphere of Boreal Forest Trees

Archaea were long thought to be a group of ancient bacteria, which mainly lived in extreme environments. Due to the development of DNA sequencing methods and molecular phylogenetic analyses, it was shown that the living organisms are in fact divided into three domains; the Archaea, Bacteria and the Eucarya. Since the beginning of the previous decade, it was shown that archaea generally inhabit moderate environments and that these non-extremophilic archaea are more ubiquitous than the extremophiles. Group 1 of non-extreme archaea affiliate with the phylum Crenarchaeota. The most commonly found soil archaea belong to the subgroup 1.1b. However, the Crenarchaeota found in the Fennoscandian boreal forest soil belong to the subgroup 1.1c. The organic top layer of the boreal forest soil, the humus, is dominated by ectomycorrhizal fungal hyphae. These colonise virtually all tree fine root tips in the humus layer and have been shown to harbour distinct bacterial populations different from those in the humus. The archaea have also been shown to colonise both boreal forest humus and the rhizospheres of plants. In this work, studies on the archaeal communities in the ectomycorrhizospheres of boreal forest trees were conducted in microcosms. Archaea belonging to the group 1.1c Crenarchaeota and Euryarchaeota of the genera Halobacterium and Methanolobus were detected. The archaea generally colonised fungal habitats, such as ectomycorrhizas and external mycelia, rather than the non-mycorrhizal fine roots of trees. The species of ectomycorrhizal fungus had a great impact on the archaeal community composition. A stable euryarchaeotal community was detected especially in the mycorrhizas, of most of the tested Scots pine colonising ectomycorrhizal fungi. The Crenarchaeota appeared more sporadically in these habitats, but had a greater diversity than the Euryarchaeota. P. involutus mycorrhizas had a higher diversity of 1.1c Crenarchaeota than the other ectomycorrhizal fungi. The detection level of archaea in the roots of boreal trees was generally low although archaea have been shown to associate with roots of different plants. However, alder showed a high diversity of 1.1c Crenarchaeota, exceeding that of any of the tested mycorrhizas. The archaeal 16S rRNA genes detected from the non-mycorrhizal roots were different from those of the P. involutus mycorrhizas. In the phylogenetic analyses, the archaeal 16S rRNA gene sequences obtained from non-mycorrhizal fine roots fell in a separate cluster within the group 1.1c Crenarchaeota than those from the mycorrhizas. When the roots of the differrent tree species were colonised by P. involutus, the diversity and frequency of the archaeal populations of the different tree species were more similar to each other. Both Cren- and Euryarchaeota were enriched in cultures to which C-1 substrates were added. The 1.1c Crenarchaeota grew anaerobically in mineral medium with CH4 and CO2 as the only available C sources, and in yeast extract media with CO2 and CH4 or H2. The crenarchaeotal diversity was higher in aerobic cultures on mineral medium with CH4 or CH3OH than in the anaerobic cultures. Ecological functions of the mycorrhizal 1.1c Crenarchaeota in both anaerobic and aerobic cycling of C-1 compounds were indicated. The phylogenetic analyses did not divide the detected Crenarchaeota into anaerobic and aerobic groups. This may suggest that the mycorrhizospheric crenarchaeotal communities consist of closely related groups of anaerobic and aerobic 1.1c Crenarchaeota, or the 1.1c Crenarchaeota may be facultatively anaerobic. Halobacteria were enriched in non-saline anaerobic yeast extract medium cultures in which CH4 was either added or produced, but were not detected in the aerobic cultures. They may potentially be involved in anaerobic CH4 cycling in ectomycorrhizas. The CH4 production of the mycorrhizal samples was over 10 times higher than for humus devoid of mycorrhizal hyphae, indicating a high CH4 production potential of the mycorrhizal metanogenic community. Autofluorescent methanogenic archaea were detected by microscopy and 16S rRNA gene sequences of the genus Methanolobus were obtained. The archaeal community depended on both tree species and the type of ectomycorrhizal fungus colonising the roots and the Cren- and Euryarchaeota may have different ecological functions in the different parts of the boreal forest tree rhizosphere and mycorrhizosphere. By employing the results of this study, it may be possible to isolate both 1.1c Crenarchaeota as well as non-halophilic halobacteria and aerotolerant methanogens from mycorrhizospheres. These archaea may be used as indicators for change in the boreal forest soil ecosystem due to different factors, such as exploitations of forests and the rise in global temperature. More information about the microbial populations with apparently low cell numbers but significant ecological impacts, such as the boreal forest soil methanogens, may be of crucial importance to counteract human impacts on such globally important ecosystems as the boreal forests.