Widespread occurrence of non-canonical transcription termination by human RNA polymerase III.

Human RNA polymerase (Pol) III-transcribed genes are thought to share a simple termination signal constituted by four or more consecutive thymidine residues in the coding DNA strand, just downstream of the RNA 3'-end sequence. We found that a large set of human tRNA genes (tDNAs) do not display any T(≥4) stretch within 50 bp of 3'-flanking region. In vitro analysis of tDNAs with a distanced T(≥4) revealed the existence of non-canonical terminators resembling degenerate T(≥5) elements, which ensure significant termination but at the same time allow for the production of Pol III read-through pre-tRNAs with unusually long 3' trailers. A panel of such non-canonical signals was found to direct transcription termination of unusual Pol III-synthesized viral pre-miRNA transcripts in gammaherpesvirus 68-infected cells. Genome-wide location analysis revealed that human Pol III tends to trespass into the 3'-flanking regions of tDNAs, as expected from extensive terminator read-through. The widespread occurrence of partial termination suggests that the Pol III primary transcriptome in mammals is unexpectedly enriched in 3'-trailer sequences with the potential to contribute novel functional ncRNAs.

Antibiotic-dependent correlation between drug-induced killing and loss of luminescence in Streptococcus gordonii expressing luciferase.

Measuring antibiotic-induced killing relies on time-consuming biological tests. The firefly luciferase gene (luc) was successfully used as a reporter gene to assess antibiotic efficacy rapidly in slow-growing Mycobacterium tuberculosis. We tested whether luc expression could also provide a rapid evaluation of bactericidal drugs in Streptococcus gordonii. The suicide vectors pFW5luc and a modified version of pJDC9 carrying a promoterless luc gene were used to construct transcriptional-fusion mutants. One mutant susceptible to penicillin-induced killing (LMI2) and three penicillin-tolerant derivatives (LMI103, LMI104, and LMI105) producing luciferase under independent streptococcal promoters were tested. The correlation between antibiotic-induced killing and luminescence was determined with mechanistically unrelated drugs. Chloramphenicol (20 times the MIC) inhibited bacterial growth. In parallel, luciferase stopped increasing and remained stable, as determined by luminescence and Western blots. Ciprofloxacin (200 times the MIC) rapidly killed 1.5 log10 CFU/ml in 2-4 hr. Luminescence decreased simultaneously by 10-fold. In contrast, penicillin (200 times the MIC) gave discordant results. Although killing was slow (< or = 0.5 log10 CFU/ml in 2 hr), luminescence dropped abruptly by 50-100-times in the same time. Inactivating penicillin with penicillinase restored luminescence, irrespective of viable counts. This was not due to altered luciferase expression or stability, suggesting some kind of post-translational modification. Luciferase shares homology with aminoacyl-tRNA synthetase and acyl-CoA ligase, which might be regulated by macromolecule synthesis and hence affected in penicillin-inhibited cells. Because of resemblance, luciferase might be down-regulated simultaneously. Luminescence cannot be universally used to predict antibiotic-induced killing. Thus, introducing reporter enzymes sharing mechanistic similarities with normal metabolic reactions might reveal other effects than those expected.

Application of PCR ribotyping and tDNA-PCR for Klebsiella pneumoniae identification

PCR analysis of 16S-23S internal transcribed spacer (PCR ribotyping) and tRNA intergenic spacer (tDNA-PCR) were evaluated for their effectiveness in identification of clinical strains of Klebsiella pneumoniae and differentiation with related species. For this purpose both methods were applied to forty-three clinical isolates biochemically identified as K. pneumoniae subsp. pneumoniae isolated from patients clinical specimens attended at five hospitals in three Brazilian cities. References strains of K. pneumoniae subsp. pneumoniae, K. pneumoniae subsp. ozaenae, K. oxytoca, K. planticola and Enterobacter aerogenes were also analyzed. Both PCR methods showed specific patterns for each species. A conserved PCR ribotype pattern was observed for all clinical K. pneumoniae isolates, while differing from other related analyzed species. tDNA-PCR revealed five distinct patterns among the K. pneumoniae clinical isolates studied, demonstrating a predominant group with 90,6% of isolates presenting the same pattern of K. pneumoniae type strain. Both PCR-based methods were not able to differentiate K. pneumoniae subspecies. On the basis of the results obtained, both methods were efficient to differentiate the Klebsiella species analyzed, as well as E. aerogenes. Meanwhile tDNA-PCR revealed different tRNA arrangements in K. pneumoniae, suggesting intra-species heterogeneity of their genome organization, the polymorphism of the intergenic spacers between 16S and 23S rRNA genes appears to be highly conserved whithin K. pneumoniae clinical isolates, showing that PCR ribotyping can be an useful tool for identification of K. pneumoniae isolates.

Phylogenetic analysis of Biomphalaria tenagophila (Orbigny, 1835) (Mollusca: Gastropoda)

Mitochondrial DNA of Biomphalaria tenagophila, a mollusc intermediate host of Schistosoma mansoni in Brazil, was sequenced and characterised. The genome size found for B. tenagophila was 13,722 bp and contained 13 messenger RNAs, 22 transfer RNAs (tRNA) and two ribosomal RNAs (rRNA). In addition to sequencing, the mitochondrial DNA (mtDNA) genome organization of B. tenagophila was analysed based on its content and localization of both coding and non-coding regions, regions of gene overlap and tRNA nucleotide sequences. Sequences of protein, rRNA 12S and rRNA 16S nucleotides as well as gene organization were compared between B. tenagophila and Biomphalaria glabrata, as the latter is the most important S. mansoni intermediate host in Brazil. Differences between such species were observed regarding rRNA composition. The complete sequence of the B. tenagophila mitochondrial genome was deposited in GenBank (accession EF433576). Furthermore, phylogenetic relationships were estimated among 28 mollusc species, which had their complete mitochondrial genome deposited in GenBank, using the neighbour-joining method, maximum parsimony and maximum likelihood bootstrap. B. tenagophila was positioned at a branch close to B. glabrata and Pulmonata molluscs, collectively comprising a paraphyletic group, contrary to Opistobranchia, which was positioned at a single branch and constituted a monophyletic group.

Intracellular excision and reintegration dynamics of the ICEclc genomic island of Pseudomonas knackmussii sp. strain B13.

Genomic islands are DNA elements acquired by horizontal gene transfer that are common to a large number of bacterial genomes, which can contribute specific adaptive functions, e.g. virulence, metabolic capacities or antibiotic resistances. Some genomic islands are still self-transferable and display an intricate life-style, reminiscent of both bacteriophages and conjugative plasmids. Here we studied the dynamical process of genomic island excision and intracellular reintegration using the integrative and conjugative element ICEclc from Pseudomonas knackmussii B13 as model. By using self-transfer of ICEclc from strain B13 to Pseudomonas putida and Cupriavidus necator as recipients, we show that ICEclc can target a number of different tRNA(Gly) genes in a bacterial genome, but only those which carry the GCC anticodon. Two conditional traps were designed for ICEclc based on the attR sequence, and we could show that ICEclc will insert with different frequencies in such traps producing brightly fluorescent cells. Starting from clonal primary transconjugants we demonstrate that ICEclc is excising and reintegrating at detectable frequencies, even in the absence of recipient. Recombination site analysis provided evidence to explain the characteristics of a larger number of genomic island insertions observed in a variety of strains, including Bordetella petri, Pseudomonas aeruginosa and Burkholderia.

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.

Fine-Tuning Translation Kinetics Selection as the Driving Force of Codon Usage Bias in the Hepatitis A Virus Capsid

Hepatitis A virus (HAV), the prototype of genus Hepatovirus, has several unique biological characteristics that distinguish it from other members of the Picornaviridae family. Among these, the need for an intact eIF4G factor for the initiation of translation results in an inability to shut down host protein synthesis by a mechanism similar to that of other picornaviruses. Consequently, HAV must inefficiently compete for the cellular translational machinery and this may explain its poor growth in cell culture. In this context of virus/cell competition, HAV has strategically adopted a naturally highly deoptimized codon usage with respect to that of its cellular host. With the aim to optimize its codon usage the virus was adapted to propagate in cells with impaired protein synthesis, in order to make tRNA pools more available for the virus. A significant loss of fitness was the immediate response to the adaptation process that was, however, later on recovered and more associated to a re-deoptimization rather than to an optimization of the codon usage specifically in the capsid coding region. These results exclude translation selection and instead suggest fine-tuning translation kinetics selection as the underlying mechanism of the codon usage bias in this specific genome region. Additionally, the results provide clear evidence of the Red Queen dynamics of evolution since the virus has very much evolved to re-adapt its codon usage to the environmental cellular changing conditions in order to recover the original fitness.

A review of the Tripterygion tripteronotus (Risso, 1810) complex, with a description of a new species from the Mediterranean Sea (Teleostei: Tripterygiidae)

We compared specimens of Tripterygion tripteronotus from 52 localities of the Mediterranean Sea and adjacent waters, using four gene sequences (12S rRNA, tRNA-valine, 16S rRNA and COI) and morphological characters. Two well-differentiated clades with a mean genetic divergence of 6.89±0.73% were found with molecular data, indicating the existence of two different species. These two species have disjunctive geographic distribution areas without any molecular hybrid populations. Subtle but diagnostic morphological differences were also present between the two species. T. tripteronotus is restricted to the northern Mediterranean basin, from the NE coast of Spain to Greece and Turkey, including the islands of Malta and Cyprus. T. tartessicum n. sp. is geographically distributed along the southern coast of Spain, from Cape of La Nao to the Gulf of Cadiz, the Balearic Islands and northern Africa, from Morocco to Tunisia. According to molecular data, these two species could have diverged during the Pliocene glaciations 2.7-3.6 Mya.

Fine-Tuning Translation Kinetics Selection as the Driving Force of Codon Usage Bias in the Hepatitis A Virus Capsid

Hepatitis A virus (HAV), the prototype of genus Hepatovirus, has several unique biological characteristics that distinguish it from other members of the Picornaviridae family. Among these, the need for an intact eIF4G factor for the initiation of translation results in an inability to shut down host protein synthesis by a mechanism similar to that of other picornaviruses. Consequently, HAV must inefficiently compete for the cellular translational machinery and this may explain its poor growth in cell culture. In this context of virus/cell competition, HAV has strategically adopted a naturally highly deoptimized codon usage with respect to that of its cellular host. With the aim to optimize its codon usage the virus was adapted to propagate in cells with impaired protein synthesis, in order to make tRNA pools more available for the virus. A significant loss of fitness was the immediate response to the adaptation process that was, however, later on recovered and more associated to a re-deoptimization rather than to an optimization of the codon usage specifically in the capsid coding region. These results exclude translation selection and instead suggest fine-tuning translation kinetics selection as the underlying mechanism of the codon usage bias in this specific genome region. Additionally, the results provide clear evidence of the Red Queen dynamics of evolution since the virus has very much evolved to re-adapt its codon usage to the environmental cellular changing conditions in order to recover the original fitness.

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver.

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes.

BRF1 mutations alter RNA polymerase III-dependent transcription and cause neurodevelopmental anomalies.

RNA polymerase III (Pol III) synthesizes tRNAs and other small noncoding RNAs to regulate protein synthesis. Dysregulation of Pol III transcription has been linked to cancer, and germline mutations in genes encoding Pol III subunits or tRNA processing factors cause neurogenetic disorders in humans, such as hypomyelinating leukodystrophies and pontocerebellar hypoplasia. Here we describe an autosomal recessive disorder characterized by cerebellar hypoplasia and intellectual disability, as well as facial dysmorphic features, short stature, microcephaly, and dental anomalies. Whole-exome sequencing revealed biallelic missense alterations of BRF1 in three families. In support of the pathogenic potential of the discovered alleles, suppression or CRISPR-mediated deletion of brf1 in zebrafish embryos recapitulated key neurodevelopmental phenotypes; in vivo complementation showed all four candidate mutations to be pathogenic in an apparent isoform-specific context. BRF1 associates with BDP1 and TBP to form the transcription factor IIIB (TFIIIB), which recruits Pol III to target genes. We show that disease-causing mutations reduce Brf1 occupancy at tRNA target genes in Saccharomyces cerevisiae and impair cell growth. Moreover, BRF1 mutations reduce Pol III-related transcription activity in vitro. Taken together, our data show that BRF1 mutations that reduce protein activity cause neurodevelopmental anomalies, suggesting that BRF1-mediated Pol III transcription is required for normal cerebellar and cognitive development.

A genomic island present along the bacterial chromosome of the Parachlamydiaceae UWE25, an obligate amoebal endosymbiont, encodes a potentially functional F-like conjugative DNA transfer system.

BACKGROUND: The genome of Protochlamydia amoebophila UWE25, a Parachlamydia-related endosymbiont of free-living amoebae, was recently published, providing the opportunity to search for genomic islands (GIs). RESULTS: On the residual cumulative G+C content curve, a G+C-rich 19-kb region was observed. This sequence is part of a 100-kb chromosome region, containing 100 highly co-oriented ORFs, flanked by two 17-bp direct repeats. Two identical gly-tRNA genes in tandem are present at the proximal end of this genetic element. Several mobility genes encoding transposases and bacteriophage-related proteins are located within this chromosome region. Thus, this region largely fulfills the criteria of GIs. The G+C content analysis shows that several modules compose this GI. Surprisingly, one of them encodes all genes essential for F-like conjugative DNA transfer (traF, traG, traH, traN, traU, traW, and trbC), involved in sex pilus retraction and mating pair stabilization, strongly suggesting that, similarly to the other F-like operons, the parachlamydial tra unit is devoted to DNA transfer. A close relatedness of this tra unit to F-like tra operons involved in conjugative transfer is confirmed by phylogenetic analyses performed on concatenated genes and gene order conservation. These analyses and that of gly-tRNA distribution in 140 GIs suggest a proteobacterial origin of the parachlamydial tra unit. CONCLUSIONS: A GI of the UWE25 chromosome encodes a potentially functional F-like DNA conjugative system. This is the first hint of a putative conjugative system in chlamydiae. Conjugation most probably occurs within free-living amoebae, that may contain hundreds of Parachlamydia bacteria tightly packed in vacuoles. Such a conjugative system might be involved in DNA transfer between internalized bacteria. Since this system is absent from the sequenced genomes of Chlamydiaceae, we hypothesize that it was acquired after the divergence between Parachlamydiaceae and Chlamydiaceae, when the Parachlamydia-related symbiont was an intracellular bacteria. It suggests that this heterologous DNA was acquired from a phylogenetically-distant bacteria sharing an amoebal vacuole. Since Parachlamydiaceae are emerging agents of pneumonia, this GI might be involved in pathogenicity. In future, conjugative systems might be developed as genetic tools for Chlamydiales.

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.

Complete DNA sequence of Kuraishia capsulata illustrates novel genomic features among budding yeasts (Saccharomycotina)

The numerous yeast genome sequences presently available provide a rich source of information for functional as well as evolutionary genomics but unequally cover the large phylogenetic diversity of extant yeasts. We present here the complete sequence of the nuclear genome of the haploid-type strain of Kuraishia capsulata (CBS1993(T)), a nitrate-assimilating Saccharomycetales of uncertain taxonomy, isolated from tunnels of insect larvae underneath coniferous barks and characterized by its copious production of extracellular polysaccharides. The sequence is composed of seven scaffolds, one per chromosome, totaling 11.4 Mb and containing 6,029 protein-coding genes, ~13.5% of which being interrupted by introns. This GC-rich yeast genome (45.7%) appears phylogenetically related with the few other nitrate-assimilating yeasts sequenced so far, Ogataea polymorpha, O. parapolymorpha, and Dekkera bruxellensis, with which it shares a very reduced number of tRNA genes, a novel tRNA sparing strategy, and a common nitrate assimilation cluster, three specific features to this group of yeasts. Centromeres were recognized in GC-poor troughs of each scaffold. The strain bears MAT alpha genes at a single MAT locus and presents a significant degree of conservation with Saccharomyces cerevisiae genes, suggesting that it can perform sexual cycles in nature, although genes involved in meiosis were not all recognized. The complete absence of conservation of synteny between K. capsulata and any other yeast genome described so far, including the three other nitrate-assimilating species, validates the interest of this species for long-range evolutionary genomic studies among Saccharomycotina yeasts.

Loss of the RNA polymerase III repressor MAF1 confers obesity resistance.

MAF1 is a global repressor of RNA polymerase III transcription that regulates the expression of highly abundant noncoding RNAs in response to nutrient availability and cellular stress. Thus, MAF1 function is thought to be important for metabolic economy. Here we show that a whole-body knockout of Maf1 in mice confers resistance to diet-induced obesity and nonalcoholic fatty liver disease by reducing food intake and increasing metabolic inefficiency. Energy expenditure in Maf1(-/-) mice is increased by several mechanisms. Precursor tRNA synthesis was increased in multiple tissues without significant effects on mature tRNA levels, implying increased turnover in a futile tRNA cycle. Elevated futile cycling of hepatic lipids was also observed. Metabolite profiling of the liver and skeletal muscle revealed elevated levels of many amino acids and spermidine, which links the induction of autophagy in Maf1(-/-) mice with their extended life span. The increase in spermidine was accompanied by reduced levels of nicotinamide N-methyltransferase, which promotes polyamine synthesis, enables nicotinamide salvage to regenerate NAD(+), and is associated with obesity resistance. Consistent with this, NAD(+) levels were increased in muscle. The importance of MAF1 for metabolic economy reveals the potential for MAF1 modulators to protect against obesity and its harmful consequences.