Analysis of twelve genomes of the bacterium Kerstersia gyiorum from brown-throated sloths (Bradypus variegatus), the first from a non-human host

Kerstersia gyiorum is a Gram-negative bacterium found in various animals, including humans, where it has been associated with various infections. Knowledge of the basic biology of K. gyiorum is essential to understand the evolutionary strategies of niche adaptation and how this organism contributes to infectious diseases; however, genomic data about K. gyiorum is very limited, especially from non-human hosts. In this work, we sequenced 12 K. gyiorum genomes isolated from healthy free-living brown-throated sloths (Bradypus variegatus) in the Parque Estadual das Fontes do Ipiranga (São Paulo, Brazil), and compared them with genomes from isolates of human origin, in order to gain insights into genomic diversity, phylogeny, and host specialization of this species. Phylogenetic analysis revealed that these K. gyiorum strains are structured according to host. Despite the fact that sloth isolates were sampled from a single geographic location, the intra-sloth K. gyiorum diversity was divided into three clusters, with differences of more than 1,000 single nucleotide polymorphisms between them, suggesting the circulation of various K. gyiorum lineages in sloths. Genes involved in mobilome and defense mechanisms against mobile genetic elements were the main source of gene content variation between isolates from different hosts. Sloth-specific K. gyiorum genome features include an IncN2 plasmid, a phage sequence, and a CRISPR-Cas system. The broad diversity of defense elements in K. gyiorum (14 systems) may prevent further mobile element flow and explain the low amount of mobile genetic elements in K. gyiorum genomes. Gene content variation may be important for the adaptation of K. gyiorum to different host niches. This study furthers our understanding of diversity, host adaptation, and evolution of K. gyiorum, by presenting and analyzing the first genomes of non-human isolates.

The Exometabolome of Xylella fastidiosa in Contact with Paraburkholderia phytofirmans Supernatant Reveals Changes in Nicotinamide, Amino Acids, Biotin, and Plant Hormones

Microbial competition within plant tissues affects invading pathogens' fitness. Metabolomics is a great tool for studying their biochemical interactions by identifying accumulated metabolites. Xylella fastidiosa, a Gram-negative bacterium causing Pierce's disease (PD) in grapevines, secretes various virulence factors including cell wall-degrading enzymes, adhesion proteins, and quorum-sensing molecules. These factors, along with outer membrane vesicles, contribute to its pathogenicity. Previous studies demonstrated that co-inoculating X. fastidiosa with the Paraburkholderia phytofirmans strain PsJN suppressed PD symptoms. Here, we further investigated the interaction between the phytopathogen and the endophyte by analyzing the exometabolome of wild-type X. fastidiosa and a diffusible signaling factor (DSF) mutant lacking quorum sensing, cultivated with 20% P. phytofirmans spent media. Liquid chromatography-mass spectrometry (LC-MS) and the Method for Metabolite Annotation and Gene Integration (MAGI) were used to detect and map metabolites to genomes, revealing a total of 121 metabolites, of which 25 were further investigated. These metabolites potentially relate to host adaptation, virulence, and pathogenicity. Notably, this study presents the first comprehensive profile of X. fastidiosa in the presence of a P. phytofirmans spent media. The results highlight that P. phytofirmans and the absence of functional quorum sensing affect the ratios of glutamine to glutamate (Gln:Glu) in X. fastidiosa. Additionally, two compounds with plant metabolism and growth properties, 2-aminoisobutyric acid and gibberellic acid, were downregulated when X. fastidiosa interacted with P. phytofirmans. These findings suggest that P. phytofirmans-mediated disease suppression involves modulation of the exometabolome of X. fastidiosa, impacting plant immunity.

The XadA Trimeric Autotransporter Adhesins in Xylella fastidiosa Differentially Contribute to Cell Aggregation, Biofilm Formation, Insect Transmission and Virulence to Plants

Surface adhesion strategies are widely employed by bacterial pathogens during establishment and systemic spread in their host. A variety of cell-surface appendages such as pili, fimbriae, and afimbrial adhesins are involved in these processes. The phytopathogen Xylella fastidiosa employs several of these structures for efficient colonization of its insect and plant hosts. Among the adhesins encoded in the X. fastidiosa genome, three afimbrial adhesins, XadA1, Hsf/XadA2, and XadA3, are predicted to be trimeric autotransporters with a C-terminal YadA-anchor membrane domain. We analyzed the individual contributions of XadA1, XadA2, and XadA3 to various cellular behaviors both in vitro and in vivo. Using isogenic X. fastidiosa mutants, we found that cell-cell aggregation and biofilm formation were severely impaired in the absence of XadA3. No significant reduction of cell-surface attachment was found with any mutant under flow conditions. Acquisition by insect vectors and transmission to grapevines were reduced in the XadA3 deletion mutant. While the XadA3 mutant was hypervirulent in grapevines, XadA1 or XadA2 deletion mutants conferred lower disease severity than the wild-type strain. This insight of the importance of these adhesive proteins and their individual contributions to different aspects of X. fastidiosa biology should guide new approaches to reduce pathogen transmission and disease development. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Pseudomonas-tailed lytic phages: genome mechanical analysis and putative correlation with virion morphogenesis yield

Aim: To unveil a putative correlation between phage genome flexibility and virion morphogenesis yield. Materials & methods: A deeper analysis of the mechanical properties of three Pseudomonas aeruginosa lytic phage genomes was undertaken, together with full genome cyclizability calculations. Results & conclusion: A putative correlation was established among phage genome flexibility, eclipse timeframe and virion particle morphogenesis yield, with a more flexible phage genome leading to a higher burst size and a more rigid phage genome leading to lower burst sizes. The results obtained are highly relevant to understand the influence of the phage genome plasticity on the virion morphogenesis yield inside the infected bacterial host cells and assumes particular relevance in the actual context of bacterial resistance to antibiotics.

Isolation and Molecular Characterization of a Novel Lytic Bacteriophage That Inactivates MDR Klebsiella pneumoniae Strains

The worldwide increase in serious infections caused by multidrug-resistant (MDR) K. pneumoniae emphasizes the urgent need of new therapeutic strategies for the control of this pathogen. There is growing interest in the use of bacteriophages (or phages) to treat K. pneumoniae infections, and newly isolated phages are needed. Here, we report the isolation and physical/biological/molecular characterization of a novel lytic phage and its efficacy in the control of MDR K. pneumoniae. The phage vB_KpnS_Uniso31, referred to hereafter as phage Kpn31, was isolated from hospital wastewater using K. pneumoniae CCCD-K001 as the host. Phage Kpn31 presents a siphovirus-like morphotype and was classified as Demerecviridae; Sugarlandvirus based on its complete genome sequence. The 113,444 bp Kpn31 genome does not encode known toxins or antimicrobial resistance genes, nor does it encode depolymerases related sequences. Phage Kpn31 showed an eclipse time of 15 min and a burst size of 9.12 PFU/host cell, allowing us to conclude it replicates well in K. pneumoniae CCCD-K001 with a latency period of 30 min. Phage Kpn31 was shown to be effective against at least six MDR K. pneumoniae clinical isolates in in vitro antibacterial activity assays. Based on its features, phage Kpn31 has potential for controlling infections caused by MDR K. pneumoniae.

Cyclohexanone-Based Chalcones as Alternatives for Fuel Additives

The use of small molecules, such as chalcones and their derivatives, for more efficient fuels is in increasing demand due to environmental factors. Here, three crystal structures (BH I, II, and III) of cyclohexanone-based chalcones were synthesized and described by single-crystal X-ray diffraction and Hirshfeld surface analysis. The supramolecular modeling analysis on the hyperconjugative interaction energies and QTAIM analysis at the ωB97XD/6-311++G(d,p) level of theory were carried out to analyze the intermolecular interactions in the solid-state. The structure-property relationship, frontier molecular orbital, molecular electrostatic potential, and the experimental calorific value analysis show that the three compounds are a good alternative to be used as an additive for some fuels. Our findings represent a further step forward in the development of cheaper and more efficient fuel additives and pose an opportunity for further investigation on similar analogues.

Characterization and in vitro testing of newly isolated lytic bacteriophages for the biocontrol of Pseudomonas aeruginosa

Aim: Two lytic phages were isolated using P. aeruginosa DSM19880 as host and fully characterized. Materials & methods: Phages were characterized physicochemically, biologically and genomically. Results & conclusion: Host range analysis revealed that the phages also infect some multidrug-resistant (MDR) P. aeruginosa clinical isolates. Increasing MOI from 1 to 1000 significantly increased phage efficiency and retarded bacteria regrowth, but phage ph0034 (reduction of 7.5 log CFU/ml) was more effective than phage ph0031 (reduction of 5.1 log CFU/ml) after 24 h. Both phages belong to Myoviridae family. Genome sequencing of phages ph0031 and ph0034 showed that they do not carry toxin, virulence, antibiotic resistance and integrase genes. The results obtained are highly relevant in the actual context of bacterial resistance to antibiotics.

Proteomic and Metabolomic Analyses of Xylella fastidiosa OMV-Enriched Fractions Reveal Association with Virulence Factors and Signaling Molecules of the DSF Family

Xylella fastidiosa releases outer membrane vesicles (OMVs) known to play a role in the systemic dissemination of this pathogen. OMVs inhibit bacterial attachment to xylem wall and traffic lipases/esterases that act on the degradation of plant cell wall. Here, we extended the characterization of X. fastidiosa OMVs by identifying proteins and metabolites potentially associated with OMVs produced by Temecula1, a Pierce's disease strain, and by 9a5c and Fb7, two citrus variegated chlorosis strains. These results strengthen that one of the OMVs multiple functions is to carry determinants of virulence, such as lipases/esterases, adhesins, proteases, porins, and a pectin lyase-like protein. For the first time, we show that the two citrus variegated chlorosis strains produce X. fastidiosa diffusible signaling factor 2 (DSF2) and citrus variegated chlorosis-DSF (likewise, Temecula1) and most importantly, that these compounds of the DSF (X. fastidiosa DSF) family are associated with OMV-enriched fractions. Altogether, our findings widen the potential functions of X. fastidiosa OMVs in intercellular signaling and host-pathogen interactions.

MARVEL, a Tool for Prediction of Bacteriophage Sequences in Metagenomic Bins

Here we present MARVEL, a tool for prediction of double-stranded DNA bacteriophage sequences in metagenomic bins. MARVEL uses a random forest machine learning approach. We trained the program on a dataset with 1,247 phage and 1,029 bacterial genomes, and tested it on a dataset with 335 bacterial and 177 phage genomes. We show that three simple genomic features extracted from contig sequences were sufficient to achieve a good performance in separating bacterial from phage sequences: gene density, strand shifts, and fraction of significant hits to a viral protein database. We compared the performance of MARVEL to that of VirSorter and VirFinder, two popular programs for predicting viral sequences. Our results show that all three programs have comparable specificity, but MARVEL achieves much better performance on the recall (sensitivity) measure. This means that MARVEL should be able to identify many more phage sequences in metagenomic bins than heretofore has been possible. In a simple test with real data, containing mostly bacterial sequences, MARVEL classified 58 out of 209 bins as phage genomes; other evidence suggests that 57 of these 58 bins are novel phage sequences. MARVEL is freely available at https://github.com/LaboratorioBioinformatica/MARVEL.

MARVEL, a Tool for Prediction of Bacteriophage Sequences in Metagenomic Bins

Here we present MARVEL, a tool for prediction of double-stranded DNA bacteriophage sequences in metagenomic bins. MARVEL uses a random forest machine learning approach. We trained the program on a dataset with 1,247 phage and 1,029 bacterial genomes, and tested it on a dataset with 335 bacterial and 177 phage genomes. We show that three simple genomic features extracted from contig sequences were sufficient to achieve a good performance in separating bacterial from phage sequences: gene density, strand shifts, and fraction of significant hits to a viral protein database. We compared the performance of MARVEL to that of VirSorter and VirFinder, two popular programs for predicting viral sequences. Our results show that all three programs have comparable specificity, but MARVEL achieves much better performance on the recall (sensitivity) measure. This means that MARVEL should be able to identify many more phage sequences in metagenomic bins than heretofore has been possible. In a simple test with real data, containing mostly bacterial sequences, MARVEL classified 58 out of 209 bins as phage genomes; other evidence suggests that 57 of these 58 bins are novel phage sequences. MARVEL is freely available at https://github.com/LaboratorioBioinformatica/MARVEL.

Genome-Centric Analysis of a Thermophilic and Cellulolytic Bacterial Consortium Derived from Composting

Microbial consortia selected from complex lignocellulolytic microbial communities are promising alternatives to deconstruct plant waste, since synergistic action of different enzymes is required for full degradation of plant biomass in biorefining applications. Culture enrichment also facilitates the study of interactions among consortium members, and can be a good source of novel microbial species. Here, we used a sample from a plant waste composting operation in the São Paulo Zoo (Brazil) as inoculum to obtain a thermophilic aerobic consortium enriched through multiple passages at 60°C in carboxymethylcellulose as sole carbon source. The microbial community composition of this consortium was investigated by shotgun metagenomics and genome-centric analysis. Six near-complete (over 90%) genomes were reconstructed. Similarity and phylogenetic analyses show that four of these six genomes are novel, with the following hypothesized identifications: a new Thermobacillus species; the first Bacillus thermozeamaize genome (for which currently only 16S sequences are available) or else the first representative of a new family in the Bacillales order; the first representative of a new genus in the Paenibacillaceae family; and the first representative of a new deep-branching family in the Clostridia class. The reconstructed genomes from known species were identified as Geobacillus thermoglucosidasius and Caldibacillus debilis. The metabolic potential of these recovered genomes based on COG and CAZy analyses show that these genomes encode several glycoside hydrolases (GHs) as well as other genes related to lignocellulose breakdown. The new Thermobacillus species stands out for being the richest in diversity and abundance of GHs, possessing the greatest potential for biomass degradation among the six recovered genomes. We also investigated the presence and activity of the organisms corresponding to these genomes in the composting operation from which the consortium was built, using compost metagenome and metatranscriptome datasets generated in a previous study. We obtained strong evidence that five of the six recovered genomes are indeed present and active in that composting process. We have thus discovered three (perhaps four) new thermophillic bacterial species that add to the increasing repertoire of known lignocellulose degraders, whose biotechnological potential can now be investigated in further studies.

Phenotype overlap in Xylella fastidiosa is controlled by the cyclic di-GMP phosphodiesterase Eal in response to antibiotic exposure and diffusible signal factor-mediated cell-cell signaling

Eal is an EAL domain protein in Xylella fastidiosa homologous to one involved in resistance to tobramycin in Pseudomonas aeruginosa. EAL and HD-GYP domain proteins are implicated in the hydrolysis of the secondary messenger bis-(3'-5')-cyclic dimeric GMP (cyclic di-GMP). Cell density-dependent communication mediated by a Diffusible Signal Factor (DSF) also modulates cyclic di-GMP levels in X. fastidiosa, thereby controlling the expression of virulence genes and genes involved in insect transmission. The possible linkage of Eal to both extrinsic factors such as antibiotics and intrinsic factors such as quorum sensing, and whether both affect virulence, was thus addressed. Expression of eal was induced by subinhibitory concentrations of tobramycin, and an eal deletion mutant was more susceptible to this antibiotic than the wild-type strain and exhibited phenotypes similar to those of an rpfF deletion mutant blocked in DSF production, such as hypermotility, reduced biofilm formation, and hypervirulence to grape. Consistent with that, the rpfF mutant was more susceptible than the wild-type strain to tobramycin. Therefore, we propose that cell-cell communication and antibiotic stress can apparently lead to similar modulations of cyclic di-GMP in X. fastidiosa, resulting in similar phenotypes. However, the effect of cell density is dominant compared to that of antibiotic stress, since eal is suppressed by RpfF, which may prevent inappropriate behavioral changes in response to antibiotic stress when DSF accumulates.

Effects of the antimicrobial peptide gomesin on the global gene expression profile, virulence and biofilm formation of Xylella fastidiosa

In the xylem vessels of susceptible hosts, such as citrus trees, Xylella fastidiosa forms biofilm-like colonies that can block water transport, which appears to correlate to disease symptoms. Besides aiding host colonization, bacterial biofilms play an important role in resistance against antimicrobial agents, for instance antimicrobial peptides (AMPs). Here, we show that gomesin, a potent AMP from a tarantula spider, modulates X. fastidiosa gene expression profile upon 60 min of treatment with a sublethal concentration. DNA microarray hybridizations revealed that among the upregulated coding sequences, some are related to biofilm production. In addition, we show that the biofilm formed by gomesin-treated bacteria is thicker than that formed by nontreated cells or cells exposed to streptomycin. We have also observed that the treatment of X. fastidiosa with a sublethal concentration of gomesin before inoculation in tobacco plants correlates with a reduction in foliar symptoms, an effect possibly due to the trapping of bacterial cells to fewer xylem vessels, given the enhancement in biofilm production. These results warrant further investigation of how X. fastidiosa would respond to the AMPs produced by citrus endophytes and by the insect vector, leading to a better understanding of the mechanism of action of these molecules on bacterial virulence.

Novel insights into the genomic basis of citrus canker based on the genome sequences of two strains of Xanthomonas fuscans subsp. aurantifolii

BACKGROUND

Citrus canker is a disease that has severe economic impact on the citrus industry worldwide. There are three types of canker, called A, B, and C. The three types have different phenotypes and affect different citrus species. The causative agent for type A is Xanthomonas citri subsp. citri, whose genome sequence was made available in 2002. Xanthomonas fuscans subsp. aurantifolii strain B causes canker B and Xanthomonas fuscans subsp. aurantifolii strain C causes canker C.

RESULTS

We have sequenced the genomes of strains B and C to draft status. We have compared their genomic content to X. citri subsp. citri and to other Xanthomonas genomes, with special emphasis on type III secreted effector repertoires. In addition to pthA, already known to be present in all three citrus canker strains, two additional effector genes, xopE3 and xopAI, are also present in all three strains and are both located on the same putative genomic island. These two effector genes, along with one other effector-like gene in the same region, are thus good candidates for being pathogenicity factors on citrus. Numerous gene content differences also exist between the three cankers strains, which can be correlated with their different virulence and host range. Particular attention was placed on the analysis of genes involved in biofilm formation and quorum sensing, type IV secretion, flagellum synthesis and motility, lipopolysacharide synthesis, and on the gene xacPNP, which codes for a natriuretic protein.

CONCLUSION

We have uncovered numerous commonalities and differences in gene content between the genomes of the pathogenic agents causing citrus canker A, B, and C and other Xanthomonas genomes. Molecular genetics can now be employed to determine the role of these genes in plant-microbe interactions. The gained knowledge will be instrumental for improving citrus canker control.

The P450 oxidoreductase, RedA, controls development beyond the mound stage in Dictyostelium discoideum

BACKGROUND

NADPH-cytochrome-P450 oxidoreductase (CPR) is a ubiquitous enzyme that belongs to a family of diflavin oxidoreductases and is required for activity of the microsomal cytochrome-P450 monooxygenase system. CPR gene-disruption experiments have demonstrated that absence of this enzyme causes developmental defects both in mouse and insect.

RESULTS

Annotation of the sequenced genome of D. discoideum revealed the presence of three genes (redA, redB and redC) that encode putative members of the diflavin oxidoreductase protein family. redA transcripts are present during growth and early development but then decline, reaching undetectable levels after the mound stage. redB transcripts are present in the same levels during growth and development while redC expression was detected only in vegetative growing cells. We isolated a mutant strain of Dictyostelium discoideum following restriction enzyme-mediated integration (REMI) mutagenesis in which redA was disrupted. This mutant develops only to the mound stage and accumulates a bright yellow pigment. The mound-arrest phenotype is cell-autonomous suggesting that the defect occurs within the cells rather than in intercellular signaling.

CONCLUSION

The developmental arrest due to disruption of redA implicates CPR in the metabolism of compounds that control cell differentiation.

Cyclic AMP and calcium interplay as second messengers in melatonin-dependent regulation of Plasmodium falciparum cell cycle

The host hormone melatonin increases cytoplasmic Ca(2+) concentration and synchronizes Plasmodium cell cycle (Hotta, C.T., M.L. Gazarini, F.H. Beraldo, F.P. Varotti, C. Lopes, R.P. Markus, T. Pozzan, and C.R. Garcia. 2000. Nat. Cell Biol. 2:466-468). Here we show that in Plasmodium falciparum melatonin induces an increase in cyclic AMP (cAMP) levels and cAMP-dependent protein kinase (PKA) activity (40 and 50%, respectively). When red blood cells infected with P. falciparum are treated with cAMP analogue adenosine 3',5'-cyclic monophosphate N6-benzoyl/PKA activator (6-Bz-cAMP) there is an alteration of the parasite cell cycle. This effect appears to depend on activation of PKA (abolished by the PKA inhibitors adenosine 3',5'-cyclic monophosphorothioate/8 Bromo Rp isomer, PKI [cell permeable peptide], and H89). An unexpected cross talk was found to exist between the cAMP and the Ca(2+)-dependent signaling pathways. The increases in cAMP by melatonin are inhibited by blocker of phospholipase C U73122, and addition of 6-Bz-cAMP increases cytosolic Ca(2+) concentration, through PKA activation. These findings suggest that in Plasmodium a highly complex interplay exists between the Ca(2+) and cAMP signaling pathways, but also that the control of the parasite cell cycle by melatonin requires the activation of both second messenger controlled pathways.

Identification and domain mapping of Dictyostelium discoideum type-1 protein phosphatase inhibitor-2

The protein phosphatase type-1 catalytic subunit (PP1c) does not exist freely in the cell and its activity must be very strictly controlled. Several protein inhibitors of PP1c have been described including the classical mammalian inhibitor-1 (I-1) and inhibitor-2 (I-2). Association of these inhibitors with PP1c appears to involve multiple contacts and in the case of I-2 no less than five I-2 interaction subdomains have been proposed. In this report, we provide both in vitro and in vivo evidence that the Dictyostelium discoideum genome encodes a protein (DdI-2) that is an ortholog of mammalian I-2, being the first PP1c interacting protein characterized in this social amoeba. Despite the low overall sequence similarity of DdI-2 with other I-2 sequences and its long N-terminal extension, the five PP1c interaction motifs proposed for mammalian I-2 are reasonably conserved in the Dictyostelium ortholog. We demonstrate that DdI-2 interacts with and inhibits D. discoideum PP1c (DdPP1c), which we have previously characterized. Moreover, using yeast two-hybrid assays we show that a stable interaction of DdI-2 with DdPP1c requires multiple contacts.

Androgen responsive intronic non-coding RNAs

BACKGROUND

Transcription of large numbers of non-coding RNAs originating from intronic regions of human genes has been recently reported, but mechanisms governing their biosynthesis and biological functions are largely unknown. In this work, we evaluated the existence of a common mechanism of transcription regulation shared by protein-coding mRNAs and intronic RNAs by measuring the effect of androgen on the transcriptional profile of a prostate cancer cell line.

RESULTS

Using a custom-built cDNA microarray enriched in intronic transcribed sequences, we found 39 intronic non-coding RNAs for which levels were significantly regulated by androgen exposure. Orientation-specific reverse transcription-PCR indicated that 10 of the 13 were transcribed in the antisense direction. These transcripts are long (0.5-5 kb), unspliced and apparently do not code for proteins. Interestingly, we found that the relative levels of androgen-regulated intronic transcripts could be correlated with the levels of the corresponding protein-coding gene (asGAS6 and asDNAJC3) or with the alternative usage of exons (asKDELR2 and asITGA6) in the corresponding protein-coding transcripts. Binding of the androgen receptor to a putative regulatory region upstream from asMYO5A, an androgen-regulated antisense intronic transcript, was confirmed by chromatin immunoprecipitation.

CONCLUSION

Altogether, these results indicate that at least a fraction of naturally transcribed intronic non-coding RNAs may be regulated by common physiological signals such as hormones, and further corroborate the notion that the intronic complement of the transcriptome play functional roles in the human gene-expression program.

Genome mapping and expression analyses of human intronic noncoding RNAs reveal tissue-specific patterns and enrichment in genes related to regulation of transcription

BACKGROUND

RNAs transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression. However, the complement of human genes in which introns are transcribed, and the number of intronic transcriptional units and their tissue expression patterns are not known.

RESULTS

A survey of mRNA and EST public databases revealed more than 55,000 totally intronic noncoding (TIN) RNAs transcribed from the introns of 74% of all unique RefSeq genes. Guided by this information, we designed an oligoarray platform containing sense and antisense probes for each of 7,135 randomly selected TIN transcripts plus the corresponding protein-coding genes. We identified exonic and intronic tissue-specific expression signatures for human liver, prostate and kidney. The most highly expressed antisense TIN RNAs were transcribed from introns of protein-coding genes significantly enriched (p = 0.002 to 0.022) in the 'Regulation of transcription' Gene Ontology category. RNA polymerase II inhibition resulted in increased expression of a fraction of intronic RNAs in cell cultures, suggesting that other RNA polymerases may be involved in their biosynthesis. Members of a subset of intronic and protein-coding signatures transcribed from the same genomic loci have correlated expression patterns, suggesting that intronic RNAs regulate the abundance or the pattern of exon usage in protein-coding messages.

CONCLUSION

We have identified diverse intronic RNA expression patterns, pointing to distinct regulatory roles. This gene-oriented approach, using a combined intron-exon oligoarray, should permit further comparative analysis of intronic transcription under various physiological and pathological conditions, thus advancing current knowledge about the biological functions of these noncoding RNAs.

Large-scale Transcriptome Analyses Reveal New Genetic Marker Candidates of Head, Neck, and Thyroid Cancer

A detailed genome mapping analysis of 213,636 expressed sequence tags (EST) derived from nontumor and tumor tissues of the oral cavity, larynx, pharynx, and thyroid was done. Transcripts matching known human genes were identified; potential new splice variants were flagged and subjected to manual curation, pointing to 788 putatively new alternative splicing isoforms, the majority (75%) being insertion events. A subset of 34 new splicing isoforms (5% of 788 events) was selected and 23 (68%) were confirmed by reverse transcription-PCR and DNA sequencing. Putative new genes were revealed, including six transcripts mapped to well-studied chromosomes such as 22, as well as transcripts that mapped to 253 intergenic regions. In addition, 2,251 noncoding intronic RNAs, eventually involved in transcriptional regulation, were found. A set of 250 candidate markers for loss of heterozygosis or gene amplification was selected by identifying transcripts that mapped to genomic regions previously known to be frequently amplified or deleted in head, neck, and thyroid tumors. Three of these markers were evaluated by quantitative reverse transcription-PCR in an independent set of individual samples. Along with detailed clinical data about tumor origin, the information reported here is now publicly available on a dedicated Web site as a resource for further biological investigation. This first in silico reconstruction of the head, neck, and thyroid transcriptomes points to a wealth of new candidate markers that can be used for future studies on the molecular basis of these tumors. Similar analysis is warranted for a number of other tumors for which large EST data sets are available.

Antisense intronic non-coding RNA levels correlate to the degree of tumor differentiation in prostate cancer

A large fraction of transcripts are expressed antisense to introns of known genes in the human genome. Here we show the construction and use of a cDNA microarray platform enriched in intronic transcripts to assess their biological relevance in pathological conditions. To validate the approach, prostate cancer was used as a model, and 27 patient tumor samples with Gleason scores ranging from 5 to 10 were analyzed. We find that a considerably higher fraction (6.6%, [23/346]) of intronic transcripts are significantly correlated (P< or =0.001) to the degree of prostate tumor differentiation (Gleason score) when compared to transcripts from unannotated genomic regions (1%, [6/539]) or from exons of known genes (2%, [27/1369]). Among the top twelve transcripts most correlated to tumor differentiation, six are antisense intronic messages as shown by orientation-specific RT-PCR or Northern blot analysis with strand-specific riboprobe. Orientation-specific real-time RT-PCR with six tumor samples, confirmed the correlation (P=0.024) between the low/high degrees of tumor differentiation and antisense intronic RASSF1 transcript levels. The need to use intron arrays to reveal the transcriptome profile of antisense intronic RNA in cancer has clearly emerged.

DNA microarray-based genome comparison of a pathogenic and a nonpathogenic strain of Xylella fastidiosa delineates genes important for bacterial virulence

Xylella fastidiosa is a phytopathogenic bacterium that causes serious diseases in a wide range of economically important crops. Despite extensive comparative analyses of genome sequences of Xylella pathogenic strains from different plant hosts, nonpathogenic strains have not been studied. In this report, we show that X. fastidiosa strain J1a12, associated with citrus variegated chlorosis (CVC), is nonpathogenic when injected into citrus and tobacco plants. Furthermore, a DNA microarray-based comparison of J1a12 with 9a5c, a CVC strain that is highly pathogenic and had its genome completely sequenced, revealed that 14 coding sequences of strain 9a5c are absent or highly divergent in strain J1a12. Among them, we found an arginase and a fimbrial adhesin precursor of type III pilus, which were confirmed to be absent in the nonpathogenic strain by PCR and DNA sequencing. The absence of arginase can be correlated to the inability of J1a12 to multiply in host plants. This enzyme has been recently shown to act as a bacterial survival mechanism by down-regulating host nitric oxide production. The lack of the adhesin precursor gene is in accordance with the less aggregated phenotype observed for J1a12 cells growing in vitro. Thus, the absence of both genes can be associated with the failure of the J1a12 strain to establish and spread in citrus and tobacco plants. These results provide the first detailed comparison between a nonpathogenic strain and a pathogenic strain of X. fastidiosa, constituting an important step towards understanding the molecular basis of the disease.

RASL11A, member of a novel small monomeric GTPase gene family, is down-regulated in prostate tumors

We performed a genome-wide search for novel loci encoding for Ras-related proteins based on the genome mapping coordinates of the cancer-derived EST dataset at GenBank. Partial sequences from two novel human genes were identified and subsequently used for full length transcript cloning. RASL11A and ARL9 belong to two novel subfamilies coding for small GTPases that we found to be highly conserved among eukaryotes. The Arl9/Arl10 subfamily displays a conserved interswitch toggle that places it evolutionarily closer to the Arf family. Rasl11 proteins are more closely related to the Ras branch of GTPases. All orthologues newly identified here exhibit an Asn residue in place of the highly conserved Thr35 of the G domain, suggesting that the universal switch mechanism of small GTPases may be structurally different in this subfamily. We determined by Northern blot that RASL11A is transcribed in several human tissues and that it is down-regulated in prostate tumors as measured by quantitative real-time PCR. These results highlight a previously uncharacterized subfamily of Ras-related genes that may have a tumor suppressor role in prostate cancer.

Analysis and functional annotation of an expressed sequence tag collection for tropical crop sugarcane

To contribute to our understanding of the genome complexity of sugarcane, we undertook a large-scale expressed sequence tag (EST) program. More than 260,000 cDNA clones were partially sequenced from 26 standard cDNA libraries generated from different sugarcane tissues. After the processing of the sequences, 237,954 high-quality ESTs were identified. These ESTs were assembled into 43,141 putative transcripts. Of the assembled sequences, 35.6% presented no matches with existing sequences in public databases. A global analysis of the whole SUCEST data set indicated that 14,409 assembled sequences (33% of the total) contained at least one cDNA clone with a full-length insert. Annotation of the 43,141 assembled sequences associated almost 50% of the putative identified sugarcane genes with protein metabolism, cellular communication/signal transduction, bioenergetics, and stress responses. Inspection of the translated assembled sequences for conserved protein domains revealed 40,821 amino acid sequences with 1415 Pfam domains. Reassembling the consensus sequences of the 43,141 transcripts revealed a 22% redundancy in the first assembling. This indicated that possibly 33,620 unique genes had been identified and indicated that >90% of the sugarcane expressed genes were tagged.

The generation and utilization of a cancer-oriented representation of the human transcriptome by using expressed sequence tags

Whereas genome sequencing defines the genetic potential of an organism, transcript sequencing defines the utilization of this potential and links the genome with most areas of biology. To exploit the information within the human genome in the fight against cancer, we have deposited some two million expressed sequence tags (ESTs) from human tumors and their corresponding normal tissues in the public databases. The data currently define approximately 23,500 genes, of which only approximately 1,250 are still represented only by ESTs. Examination of the EST coverage of known cancer-related (CR) genes reveals that <1% do not have corresponding ESTs, indicating that the representation of genes associated with commonly studied tumors is high. The careful recording of the origin of all ESTs we have produced has enabled detailed definition of where the genes they represent are expressed in the human body. More than 100,000 ESTs are available for seven tissues, indicating a surprising variability of gene usage that has led to the discovery of a significant number of genes with restricted expression, and that may thus be therapeutically useful. The ESTs also reveal novel nonsynonymous germline variants (although the one-pass nature of the data necessitates careful validation) and many alternatively spliced transcripts. Although widely exploited by the scientific community, vindicating our totally open source policy, the EST data generated still provide extensive information that remains to be systematically explored, and that may further facilitate progress toward both the understanding and treatment of human cancers.

ESTWeb: bioinformatics services for EST sequencing projects

ESTWeb is an internet based software package designed for uniform data processing and storage for large-scale EST sequencing projects. The package provides for: (a) reception of sequencing chromatograms; (b) sequence processing such as base-calling, vector screening, comparison with public databases; (c) storage of data and analysis in a relational database, (d) generation of a graphical report of individual sequence quality; and (e) issuing of reports with statistics of productivity and redundancy. The software facilitates real-time monitoring and evaluation of EST sequence acquisition progress along an EST sequencing project.