Deciphering Bacterial and Archaeal Transcriptional Dark Matter and Its Architectural Complexity.

Transcripts are potential therapeutic targets, yet bacterial transcripts remain biological dark matter with uncharacterized biodiversity. We developed and applied an algorithm to predict transcripts for Escherichia coli K12 and E2348/69 strains (Bacteria:gamma-Proteobacteria) with newly generated ONT direct RNA sequencing data while predicting transcripts for Listeria monocytogenes strains Scott A and RO15 (Bacteria:Firmicute), Pseudomonas aeruginosa strains SG17M and NN2 strains (Bacteria:gamma-Proteobacteria), and Haloferax volcanii (Archaea:Halobacteria) using publicly available data. From >5 million E. coli K12 ONT direct RNA sequencing reads, 2,484 mRNAs are predicted and contain more than half of the predicted E. coli proteins. While the number of predicted transcripts varied by strain based on the amount of sequence data used for the predictions, across all strains examined, the average size of the predicted mRNAs is 1.6-1.7 kbp while the median size of the predicted bacterial 5'- and 3'- UTRs are 30-90 bp. Given the lack of bacterial and archaeal transcript annotation, most predictions are of novel transcripts, but we also predicted many previously characterized mRNAs and ncRNAs, including post-transcriptionally generated transcripts and small RNAs associated with pathogenesis in the E. coli E2348/69 LEE pathogenicity islands. We predicted small transcripts in the 100-200 bp range as well as >10 kbp transcripts for all strains, with the longest transcript for two of the seven strains being the nuo operon transcript, and for another two strains it was a phage/prophage transcript. This quick, easy, inexpensive, and reproducible method will facilitate the presentation of operons, transcripts, and UTR predictions alongside CDS and protein predictions in bacterial genome annotation as important resources for the research community.

Common Analysis of Direct RNA SequencinG CUrrently Leads to Misidentification of 5-Methylcytosine Modifications at GCU Motifs.

RNA modifications, such as méthylation, can be detected with Oxford Nanopore Technologies direct RNA sequencing. One commonly used tool for detecting 5-methylcytosine (m5C) modifications is Tombo, which uses an "Alternative Model" to detect putative modifications from a single sample. We examined direct RNA sequencing data from diverse taxa including virus, bacteria, fungi, and animals. The algorithm consistently identified a 5-methylcytosine at the central position of a GCU motif. However, it also identified a 5-methylcytosine in the same motif in fully unmodified in vitro transcribed RNA, suggesting that this a frequent false prediction. In the absence of further validation, several published predictions of 5-methylcytosine in human coronavirus and human cerebral organoid RNA in a GCU context should be reconsidered.

Pim kinase inhibitor co-treatment decreases alternative non-homologous end-joining DNA repair and genomic instability induced by topoisomerase 2 inhibitors in cells with FLT3 internal tandem duplication.

Acute myeloid leukemia (AML) with fms-like tyrosine kinase 3 internal tandem duplication (FLT3-ITD) relapses with new chromosome abnormalities following chemotherapy, implicating genomic instability. Error-prone alternative non-homologous end-joining (Alt-NHEJ) DNA double-strand break (DSB) repair is upregulated in FLT3-ITD-expresssing cells, driven by c-Myc. The serine/threonine kinase Pim-1 is upregulated downstream of FLT3-ITD, and inhibiting Pim increases topoisomerase 2 (TOP2) inhibitor chemotherapy drug induction of DNA DSBs and apoptosis. We hypothesized that Pim inhibition increases DNA DSBs by downregulating Alt-NHEJ, also decreasing genomic instability. Alt-NHEJ activity, measured with a green fluorescent reporter construct, increased in FLT3-ITD-transfected Ba/F3-ITD cells treated with TOP2 inhibitors, and this increase was abrogated by Pim kinase inhibitor AZD1208 co-treatment. TOP2 inhibitor and AZD1208 co-treatment downregulated cellular and nuclear expression of c-Myc and Alt-NHEJ repair pathway proteins DNA polymerase θ, DNA ligase 3 and XRCC1 in FLT3-ITD cell lines and AML patient blasts. ALT-NHEJ protein downregulation was preceded by c-Myc downregulation, inhibited by c-Myc overexpression and induced by c-Myc knockdown or inhibition. TOP2 inhibitor treatment increased chromosome breaks in metaphase spreads in FLT3-ITD-expressing cells, and AZD1208 co-treatment abrogated these increases. Thus Pim kinase inhibitor co-treatment both enhances TOP2 inhibitor cytotoxicity and decreases TOP2 inhibitor-induced genomic instability in cells with FLT3-ITD.

A Meta-Analysis of Wolbachia Transcriptomics Reveals a Stage-Specific Wolbachia Transcriptional Response Shared Across Different Hosts.

Wolbachia is a genus containing obligate, intracellular endosymbionts with arthropod and nematode hosts. Numerous studies have identified differentially expressed transcripts in Wolbachia endosymbionts that potentially inform the biological interplay between these endosymbionts and their hosts, albeit with discordant results. Here, we re-analyze previously published Wolbachia RNA-Seq transcriptomics data sets using a single workflow consisting of the most up-to-date algorithms and techniques, with the aim of identifying trends or patterns in the pan-Wolbachia transcriptional response. We find that data from one of the early studies in filarial nematodes did not allow for robust conclusions about Wolbachia differential expression with these methods, suggesting the original interpretations should be reconsidered. Across datasets analyzed with this unified workflow, there is a general lack of global gene regulation with the exception of a weak transcriptional response resulting in the upregulation of ribosomal proteins in early larval stages. This weak response is observed across diverse Wolbachia strains from both nematode and insect hosts suggesting a potential pan-Wolbachia transcriptional response during host development that diverged more than 700 million years ago.

The genome of Onchocerca volvulus, agent of river blindness.

Human onchocerciasis is a serious neglected tropical disease caused by the filarial nematode Onchocerca volvulus that can lead to blindness and chronic disability. Control of the disease relies largely on mass administration of a single drug, and the development of new drugs and vaccines depends on a better knowledge of parasite biology. Here, we describe the chromosomes of O. volvulus and its Wolbachia endosymbiont. We provide the highest-quality sequence assembly for any parasitic nematode to date, giving a glimpse into the evolution of filarial parasite chromosomes and proteomes. This resource was used to investigate gene families with key functions that could be potentially exploited as targets for future drugs. Using metabolic reconstruction of the nematode and its endosymbiont, we identified enzymes that are likely to be essential for O. volvulus viability. In addition, we have generated a list of proteins that could be targeted by Federal-Drug-Agency-approved but repurposed drugs, providing starting points for anti-onchocerciasis drug development.

Phenotypic, genomic, and transcriptional characterization of Streptococcus pneumoniae interacting with human pharyngeal cells.

BACKGROUND: Streptococcus pneumoniae is a leading cause of childhood morbidity and mortality worldwide, despite the availability of effective pneumococcal vaccines. Understanding the molecular interactions between the bacterium and the host will contribute to the control and prevention of pneumococcal disease. RESULTS: We used a combination of adherence assays, mutagenesis and functional genomics to identify novel factors involved in adherence. By contrasting these processes in two pneumococcal strains, TIGR4 and G54, we showed that adherence and invasion capacities vary markedly by strain. Electron microscopy showed more adherent bacteria in association with membranous pseudopodia in the TIGR4 strain. Operons for cell wall phosphorylcholine incorporation (lic), manganese transport (psa) and phosphate utilization (phn) were up-regulated in both strains on exposure to epithelial cells. Pneumolysin, pili, stress protection genes (adhC-czcD) and genes of the type II fatty acid synthesis pathway were highly expressed in the naturally more invasive strain, TIGR4. Deletion mutagenesis of five gene regions identified as regulated in this study revealed attenuation in adherence. Most strikingly, ∆SP_1922 which was predicted to contain a B-cell epitope and revealed significant attenuation in adherence, appeared to be expressed as a part of an operon that includes the gene encoding the cytoplasmic pore-forming toxin and vaccine candidate, pneumolysin. CONCLUSION: This work identifies a list of novel potential pneumococcal adherence determinants.

The two-component response regulator LiaR regulates cell wall stress responses, pili expression and virulence in group B Streptococcus.

Group B Streptococcus (GBS) remains the leading cause of early onset sepsis among term infants. Evasion of innate immune defences is critical to neonatal GBS disease pathogenesis. Effectors of innate immunity, as well as numerous antibiotics, frequently target the peptidoglycan layer of the Gram-positive bacterial cell wall. The intramembrane-sensing histidine kinase (IM-HK) class of two-component regulatory systems has been identified as important to the Gram-positive response to cell wall stress. We have characterized the GBS homologue of LiaR, the response regulator component of the Lia system, to determine its role in GBS pathogenesis. LiaR is expressed as part of a three-gene operon (liaFSR) with a promoter located upstream of liaF. A LiaR deletion mutant is more susceptible to cell wall-active antibiotics (vancomycin and bacitracin) as well as antimicrobial peptides (polymixin B, colistin, and nisin) compared to isogenic wild-type GBS. LiaR mutant GBS are significantly attenuated in mouse models of both GBS sepsis and pneumonia. Transcriptional profiling with DNA microarray and Northern blot demonstrated that LiaR regulates expression of genes involved in microbial defence against host antimicrobial systems including genes functioning in cell wall synthesis, pili formation and cell membrane modification. We conclude that the LiaFSR system, the first member of the IM-HK regulatory systems to be studied in GBS, is involved in sensing perturbations in the integrity of the cell wall and activates a transcriptional response that is important to the pathogenesis of GBS infection.

Extensive genomic diversity of closely related Wolbachia strains.

Using microarray-based comparative genome hybridization (mCGH), the genomic content of Wolbachia pipientis wMel from Drosophila melanogaster was compared to the closely related Wolbachia from D. innubila (wInn), D. santomea (wSan), and three strains from D. simulans (wAu, wRi, wSim). A large number of auxiliary genes are identified in these five strains, with most absent/divergent genes being unique to a given strain. Each strain caused an average of approximately 60 genes to be removed from the core genome. As such, these organisms do not appear to have the streamlined genomes expected of obligate intracellular bacteria. Prophage, hypothetical and ankyrin repeat genes are over-represented in the absent/divergent genes, with 21-87% of absent/divergent genes coming from prophage regions. The only wMel region absent/divergent in all five query strains is that containing WD_0509 to WD_0511, including a DNA mismatch repair protein MutL-2, a degenerate RNase, and a conserved hypothetical protein. A region flanked by the two portions of the WO-B prophage in wMel is found in four of the five Wolbachia strains as well as on a plasmid of a rickettsial endosymbiont of Ixodes scapularis, suggesting lateral gene transfer between these two obligate intracellular species. Overall, these insect-associated Wolbachia have highly mosaic genomes, with lateral gene transfer playing an important role in their diversity and evolution.

Comparative genomics of Neisseria meningitidis: core genome, islands of horizontal transfer and pathogen-specific genes.

To better understand Neisseria meningitidis genomes and virulence, microarray comparative genome hybridization (mCGH) data were collected from one Neisseria cinerea, two Neisseria lactamica, two Neisseria gonorrhoeae and 48 Neisseria meningitidis isolates. For N. meningitidis, these isolates are from diverse clonal complexes, invasive and carriage strains, and all major serogroups. The microarray platform represented N. meningitidis strains MC58, Z2491 and FAM18, and N. gonorrhoeae FA1090. By comparing hybridization data to genome sequences, the core N. meningitidis genome and insertions/deletions (e.g. capsule locus, type I secretion system) related to pathogenicity were identified, including further characterization of the capsule locus, bioinformatics analysis of a type I secretion system, and identification of some metabolic pathways associated with intracellular survival in pathogens. Hybridization data clustered meningococcal isolates from similar clonal complexes that were distinguished by the differential presence of six distinct islands of horizontal transfer. Several of these islands contained prophage or other mobile elements, including a novel prophage and a transposon carrying portions of a type I secretion system. Acquisition of some genetic islands appears to have occurred in multiple lineages, including transfer between N. lactamica and N. meningitidis. However, island acquisition occurs infrequently, such that the genomic-level relationship is not obscured within clonal complexes. The N. meningitidis genome is characterized by the horizontal acquisition of multiple genetic islands; the study of these islands reveals important sets of genes varying between isolates and likely to be related to pathogenicity.