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 m5C at GCU motifs.

RNA modifications, such as methylation, 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 viruses, bacteria, fungi, and animals. The algorithm consistently identified a m5C at the central position of a GCU motif. However, it also identified a m5C in the same motif in fully unmodified in vitro transcribed RNA, suggesting that this is a frequent false prediction. In the absence of further validation, several published predictions of m5C in a GCU context should be reconsidered, including those from human coronavirus and human cerebral organoid samples.

Immunogenomic profile at baseline predicts host susceptibility to clinical malaria.

Introduction: Host gene and protein expression impact susceptibility to clinical malaria, but the balance of immune cell populations, cytokines and genes that contributes to protection, remains incompletely understood. Little is known about the determinants of host susceptibility to clinical malaria at a time when acquired immunity is developing. Methods: We analyzed peripheral blood mononuclear cells (PBMCs) collected from children who differed in susceptibility to clinical malaria, all from a small town in Mali. PBMCs were collected from children aged 4-6 years at the start, peak and end of the malaria season. We characterized the immune cell composition and cytokine secretion for a subset of 20 children per timepoint (10 children with no symptomatic malaria age-matched to 10 children with >2 symptomatic malarial illnesses), and gene expression patterns for six children (three per cohort) per timepoint. Results: We observed differences between the two groups of children in the expression of genes related to cell death and inflammation; in particular, inflammatory genes such as CXCL10 and STAT1 and apoptotic genes such as XAF1 were upregulated in susceptible children before the transmission season began. We also noted higher frequency of HLA-DR+ CD4 T cells in protected children during the peak of the malaria season and comparable levels cytokine secretion after stimulation with malaria schizonts across all three time points. Conclusion: This study highlights the importance of baseline immune signatures in determining disease outcome. Our data suggests that differences in apoptotic and inflammatory gene expression patterns can serve as predictive markers of susceptibility to clinical malaria.

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.

Accumulation of endosymbiont genomes in an insect autosome followed by endosymbiont replacement.

Eukaryotic genomes can acquire bacterial DNA via lateral gene transfer (LGT).1 A prominent source of LGT is Wolbachia,2 a widespread endosymbiont of arthropods and nematodes that is transmitted maternally through female germline cells.3,4 The DNA transfer from the Wolbachia endosymbiont wAna to Drosophila ananassae is extensive5-7 and has been localized to chromosome 4, contributing to chromosome expansion in this lineage.6 As has happened frequently with claims of bacteria-to-eukaryote LGT, the contribution of wAna transfers to the expanded size of D. ananassae chromosome 4 has been specifically contested8 owing to an assembly where Wolbachia sequences were classified as contaminants and removed.9 Here, long-read sequencing with DNA from a Wolbachia-cured line enabled assembly of 4.9 Mbp of nuclear Wolbachia transfers (nuwts) in D. ananassae and a 24-kbp nuclear mitochondrial transfer. The nuwts are <8,000 years old in at least two locations in chromosome 4 with at least one whole-genome integration followed by rapid extensive duplication of most of the genome with regions that have up to 10 copies. The genes in nuwts are accumulating small indels and mobile element insertions. Among the highly duplicated genes are cifA and cifB, two genes associated with Wolbachia-mediated Drosophila cytoplasmic incompatibility. The wAna strain that was the source of nuwts was subsequently replaced by a different wAna endosymbiont. Direct RNA Nanopore sequencing of Wolbachia-cured lines identified nuwt transcripts, including spliced transcripts, but functionality, if any, remains elusive.

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.

Nearly Complete Genome Sequence of Brugia pahangi FR3.

Brugia pahangi is a zoonotic parasite that is closely related to human-infecting filarial nematodes. Here, we report the nearly complete genome of Brugia pahangi, including assemblies of four autosomes and an X chromosome, with only seven gaps. The Y chromosome is still not completely assembled.

Complete Genome Sequence of wBp, the Wolbachia Endosymbiont of Brugia pahangi FR3.

Lymphatic filariasis is a devastating disease caused by filarial nematode roundworms, which contain obligate Wolbachia endosymbionts. Here, we assembled the genome of wBp, the Wolbachia endosymbiont of the filarial nematode Brugia pahangi, from Illumina, Pacific Biosciences, and Oxford Nanopore data. The complete, circular genome is 1,072,967 bp.

Complete Mitochondrial Genome Sequence of Mansonella perstans.

The 13,647-bp complete mitochondrial genome of Mansonella perstans was sequenced and is syntenic to the mitochondrial genome of Mansonella ozzardi Phylogenetic analysis of the mitochondrial genome is consistent with the known phylogeny of ONC5 group filarial nematodes.

Drug Repurposing of Bromodomain Inhibitors as Potential Novel Therapeutic Leads for Lymphatic Filariasis Guided by Multispecies Transcriptomics.

To better understand the transcriptomic interplay of organisms associated with lymphatic filariasis, we conducted multispecies transcriptome sequencing (RNA-Seq) on the filarial nematode Brugia malayi, its Wolbachia endosymbiont wBm, and its laboratory vector Aedes aegypti across the entire B. malayi life cycle. In wBm, transcription of the noncoding 6S RNA suggests that it may be a regulator of bacterial cell growth, as its transcript levels correlate with bacterial replication rates. For A. aegypti, the transcriptional response reflects the stress that B. malayi infection exerts on the mosquito with indicators of increased energy demand. In B. malayi, expression modules associated with adult female samples consistently contained an overrepresentation of genes involved in chromatin remodeling, such as the bromodomain-containing proteins. All bromodomain-containing proteins encoded by B. malayi were observed to be upregulated in the adult female, embryo, and microfilaria life stages, including 2 members of the bromodomain and extraterminal (BET) protein family. The BET inhibitor JQ1(+), originally developed as a cancer therapeutic, caused lethality of adult worms in vitro, suggesting that it may be a potential therapeutic that can be repurposed for treating lymphatic filariasis.IMPORTANCE The current treatment regimen for lymphatic filariasis is mostly microfilaricidal. In an effort to identify new drug candidates for lymphatic filariasis, we conducted a three-way transcriptomics/systems biology study of one of the causative agents of lymphatic filariasis, Brugia malayi, its Wolbachia endosymbiont wBm, and its vector host Aedes aegypti at 16 distinct B. malayi life stages. B. malayi upregulates the expression of bromodomain-containing proteins in the adult female, embryo, and microfilaria stages. In vitro, we find that the existing cancer therapeutic JQ1(+), which is a bromodomain and extraterminal protein inhibitor, has adulticidal activity in B. malayi.

Complete Genome Sequence of wAna, the Wolbachia Endosymbiont of Drosophila ananassae.

Here, we present the complete genome sequence of the Wolbachia endosymbiont wAna, isolated from Drosophila ananassae and derived from Oxford Nanopore and Illumina sequencing. We anticipate that this will aid in Wolbachia comparative genomics and the assembly of D. ananassae specifically in regions containing extensive lateral gene transfer events.

TwinBLAST: When Two Is Better than One.

Analysis of sequence read pairs can be essential for characterizing structural variation, including junction-spanning pairs of reads (JSPRs) suggesting recent lateral/horizontal gene transfer. TwinBLAST can be used to facilitate this analysis of JSPRs by enabling the visualization and curation of two BLAST reports side by side in a single interface.

Multispecies Transcriptomics Data Set of Brugia malayi, Its Wolbachia Endosymbiont wBm, and Aedes aegypti across the B. malayi Life Cycle.

Here, we present a comprehensive transcriptomics data set of Brugia malayi, its Wolbachia endosymbiont wBm, and its vector host. This study samples from 16 stages across the entire B. malayi life cycle, including stage 1 through 4 larvae, adult males and females, embryos, immature microfilariae, and mature microfilariae.

Targeted enrichment outperforms other enrichment techniques and enables more multi-species RNA-Seq analyses.

Enrichment methodologies enable the analysis of minor members in multi-species transcriptomic data. We compared the standard enrichment of bacterial and eukaryotic mRNA to a targeted enrichment using an Agilent SureSelect (AgSS) capture for Brugia malayi, Aspergillus fumigatus, and the Wolbachia endosymbiont of B. malayi (wBm). Without introducing significant systematic bias, the AgSS quantitatively enriched samples, resulting in more reads mapping to the target organism. The AgSS-enriched libraries consistently had a positive linear correlation with their unenriched counterparts (r2 = 0.559-0.867). Up to a 2,242-fold enrichment of RNA from the target organism was obtained following a power law (r2 = 0.90), with the greatest fold enrichment achieved in samples with the largest ratio difference between the major and minor members. While using a single total library for prokaryote and eukaryote enrichment from a single RNA sample could be beneficial for samples where RNA is limiting, we observed a decrease in reads mapping to protein coding genes and an increase in multi-mapping reads to rRNAs in AgSS enrichments from eukaryotic total RNA libraries compared to eukaryotic poly(A)-enriched libraries. Our results support a recommendation of using AgSS targeted enrichment on poly(A)-enriched libraries for eukaryotic captures, and total RNA libraries for prokaryotic captures, to increase the robustness of multi-species transcriptomic studies.

Lateral gene transfer between prokaryotes and eukaryotes.

Lateral gene transfer (LGT) is an all-encompassing term for the movement of DNA between diverse organisms. LGT is synonymous with horizontal gene transfer, and the terms are used interchangeably throughout the scientific literature. While LGT has been recognized within the bacteria domain of life for decades, inter-domain LGTs are being increasingly described. LGTs between bacteria and complex multicellular organisms are of interest because they challenge the long-held dogma that such transfers could only occur in closely-related, single-celled organisms. Scientists will continue to challenge our understanding of LGT as we sequence more, diverse organisms, as we sequence more endosymbiont-colonized arthropods, and as we continue to appreciate LGT events, both young and old.

Large-scale intersubspecific recombination in the plant-pathogenic bacterium Xylella fastidiosa is associated with the host shift to mulberry.

Homologous recombination plays an important role in the structuring of genetic variation of many bacteria; however, its importance in adaptive evolution is not well established. We investigated the association of intersubspecific homologous recombination (IHR) with the shift to a novel host (mulberry) by the plant-pathogenic bacterium Xylella fastidiosa. Mulberry leaf scorch was identified about 25 years ago in native red mulberry in the eastern United States and has spread to introduced white mulberry in California. Comparing a sequence of 8 genes (4,706 bp) from 21 mulberry-type isolates to published data (352 isolates representing all subspecies), we confirmed previous indications that the mulberry isolates define a group distinct from the 4 subspecies, and we propose naming the taxon X. fastidiosa subsp. morus. The ancestry of its gene sequences was mixed, with 4 derived from X. fastidiosa subsp. fastidiosa (introduced from Central America), 3 from X. fastidiosa subsp. multiplex (considered native to the United States), and 1 chimeric, demonstrating that this group originated by large-scale IHR. The very low within-type genetic variation (0.08% site polymorphism), plus the apparent inability of native X. fastidiosa subsp. multiplex to infect mulberry, suggests that this host shift was achieved after strong selection acted on genetic variants created by IHR. Sequence data indicate that a single ancestral IHR event gave rise not only to X. fastidiosa subsp. morus but also to the X. fastidiosa subsp. multiplex recombinant group which infects several hosts but is the only type naturally infecting blueberry, thus implicating this IHR in the invasion of at least two novel native hosts, mulberry and blueberry.

Recent evolutionary radiation and host plant specialization in the Xylella fastidiosa subspecies native to the United States.

The bacterial pathogen, Xylella fastidiosa, infects many plant species in the Americas, making it a good model for investigating the genetics of host adaptation. We used multilocus sequence typing (MLST) to identify isolates of the native U.S. subsp. multiplex that were largely unaffected by intersubspecific homologous recombination (IHR) and to investigate how their evolutionary history influences plant host specialization. We identified 110 "non-IHR" isolates, 2 minimally recombinant "intermediate" ones (including the subspecific type), and 31 with extensive IHR. The non-IHR and intermediate isolates defined 23 sequence types (STs) which we used to identify 22 plant hosts (73% trees) characteristic of the subspecies. Except for almond, subsp. multiplex showed no host overlap with the introduced subspecies (subspecies fastidiosa and sandyi). MLST sequences revealed that subsp. multiplex underwent recent radiation (<25% of subspecies age) which included only limited intrasubspecific recombination (ρ/θ = 0.02); only one isolated lineage (ST50 from ash) was older. A total of 20 of the STs grouped into three loose phylogenetic clusters distinguished by nonoverlapping hosts (excepting purple leaf plum): "almond," "peach," and "oak" types. These host differences were not geographical, since all three types also occurred in California. ST designation was a good indicator of host specialization. ST09, widespread in the southeastern United States, only infected oak species, and all peach isolates were ST10 (from California, Florida, and Georgia). Only ST23 had a broad host range. Hosts of related genotypes were sometimes related, but often host groupings crossed plant family or even order, suggesting that phylogenetically plastic features of hosts affect bacterial pathogenicity.

Detecting genetic introgression: high levels of intersubspecific recombination found in Xylella fastidiosa in Brazil.

Documenting the role of novel mutation versus homologous recombination in bacterial evolution, and especially in the invasion of new hosts, is central to understanding the long-term dynamics of pathogenic bacteria. We used multilocus sequence typing (MLST) to study this issue in Xylella fastidiosa subsp. pauca from Brazil, a bacterium causing citrus variegated chlorosis (CVC) and coffee leaf scorch (CLS). All 55 citrus isolates typed (plus one coffee isolate) defined three similar sequence types (STs) dominated by ST11 (85%), while the remaining 22 coffee isolates defined two STs, mainly ST16 (74%). This low level of variation masked unusually large allelic differences (>1% divergence with no intermediates) at five loci (leuA, petC, malF, cysG, and holC). We developed an introgression test to detect whether these large differences were due to introgression via homologous recombination from another X. fastidiosa subspecies. Using additional sequencing around these loci, we established that the seven randomly chosen MLST targets contained seven regions of introgression totaling 2,172 bp of 4,161 bp (52%), only 409 bp (10%) of which were detected by other recombination tests. This high level of introgression suggests the hypothesis that X. fastidiosa subsp. pauca became pathogenic on citrus and coffee (crops cultivated in Brazil for several hundred years) only recently after it gained genetic variation via intersubspecific recombination, facilitating a switch from native hosts. A candidate donor is the subspecies infecting plum in the region since 1935 (possibly X. fastidiosa subsp. multiplex). This hypothesis predicts that nonrecombinant native X. fastidiosa subsp. pauca (not yet isolated) does not cause disease in citrus or coffee.