Dual species transcriptomics reveals conserved metabolic and immunologic processes in interactions between human neutrophils and Neisseria gonorrhoeae.

Neisseria gonorrhoeae (the gonococcus, Gc) causes the sexually transmitted infection gonorrhea. Gc is a prominent threat to human health by causing severe lifelong sequelae, including infertility and chronic pelvic pain, which is amplified by the emergence of "superbug" strains resistant to all current antibiotics. Gc is highly adapted to colonize human mucosal surfaces, where it survives despite initiating a robust inflammatory response and influx of polymorphonuclear leukocytes (PMNs, neutrophils) that typically clear bacteria. Here, dual-species RNA-sequencing was used to define Gc and PMN transcriptional profiles alone and after infection. Core host and bacterial responses were assessed for two strains of Gc and three human donors' PMNs. Comparative analysis of Gc transcripts revealed overlap between Gc responses to PMNs, iron, and hydrogen peroxide; 98 transcripts were differentially expressed across both Gc strains in response to PMN co-culture, including iron-responsive and oxidative stress response genes. We experimentally determined that the iron-dependent TbpB is suppressed by PMN co-culture, and iron-limited Gc have a survival advantage when cultured with PMNs. Analysis of PMN transcripts modulated by Gc infection revealed differential expression of genes driving cell adhesion, migration, inflammatory responses, and inflammation resolution pathways. Production of pro-inflammatory cytokines, including IL1B and IL8, the adhesion factor ICAM1, and prostaglandin PGE2 were induced in PMNs in response to Gc. Together, this study represents a comprehensive and experimentally validated dual-species transcriptomic analysis of two isolates of Gc and primary human PMNs that gives insight into how this bacterium survives innate immune onslaught to cause disease.

A comparative analysis of library prep approaches for sequencing low input translatome samples.

BACKGROUND: Cell type-specific ribosome-pulldown has become an increasingly popular method for analysis of gene expression. It allows for expression analysis from intact tissues and monitoring of protein synthesis in vivo. However, while its utility has been assessed, technical aspects related to sequencing of these samples, often starting with a smaller amount of RNA, have not been reported. In this study, we evaluated the performance of five library prep protocols for ribosome-associated mRNAs when only 250 pg-4 ng of total RNA are used. RESULTS: We obtained total and RiboTag-IP RNA, in three biological replicates. We compared 5 methods of library preparation for Illumina Next Generation sequencing: NuGEN Ovation RNA-Seq system V2 Kit, TaKaRa SMARTer Stranded Total RNA-Seq Kit, TaKaRa SMART-Seq v4 Ultra Low Input RNA Kit, Illumina TruSeq RNA Library Prep Kit v2 and NEBNext® Ultra™ Directional RNA Library Prep Kit using slightly modified protocols each with 4 ng of total RNA. An additional set of samples was processed using the TruSeq kit with 70 ng, as a 'gold standard' control and the SMART-Seq v4 with 250 pg of total RNA. TruSeq-processed samples had the best metrics overall, with similar results for the 4 ng and 70 ng samples. The results of the SMART-Seq v4 processed samples were similar to TruSeq (Spearman correlation > 0.8) despite using lower amount of input RNA. All RiboTag-IP samples had an increase in the intronic reads compared with the corresponding whole tissue, suggesting that the IP captures some immature mRNAs. The SMARTer-processed samples had a higher representation of ribosomal and non-coding RNAs leading to lower representation of protein coding mRNA. The enrichment or depletion of IP samples compared to corresponding input RNA was similar across all kits except for SMARTer kit. CONCLUSION: RiboTag-seq can be performed successfully with as little as 250 pg of total RNA when using the SMART-Seq v4 kit and 4 ng when using the modified protocols of other library preparation kits. The SMART-Seq v4 and TruSeq kits resulted in the highest quality libraries. RiboTag IP RNA contains some immature transcripts.

Single molecule sequencing and genome assembly of a clinical specimen of Loa loa, the causative agent of loiasis.

BACKGROUND: More than 20% of the world's population is at risk for infection by filarial nematodes and >180 million people worldwide are already infected. Along with infection comes significant morbidity that has a socioeconomic impact. The eight filarial nematodes that infect humans are Wuchereria bancrofti, Brugia malayi, Brugia timori, Onchocerca volvulus, Loa loa, Mansonella perstans, Mansonella streptocerca, and Mansonella ozzardi, of which three have published draft genome sequences. Since all have humans as the definitive host, standard avenues of research that rely on culturing and genetics have often not been possible. Therefore, genome sequencing provides an important window into understanding the biology of these parasites. The need for large amounts of high quality genomic DNA from homozygous, inbred lines; the availability of only short sequence reads from next-generation sequencing platforms at a reasonable expense; and the lack of random large insert libraries has limited our ability to generate high quality genome sequences for these parasites. However, the Pacific Biosciences single molecule, real-time sequencing platform holds great promise in reducing input amounts and generating sufficiently long sequences that bypass the need for large insert paired libraries. RESULTS: Here, we report on efforts to generate a more complete genome assembly for L. loa using genetically heterogeneous DNA isolated from a single clinical sample and sequenced on the Pacific Biosciences platform. To obtain the best assembly, numerous assemblers and sequencing datasets were analyzed, combined, and compared. Quiver-informed trimming of an assembly of only Pacific Biosciences reads by HGAP2 was selected as the final assembly of 96.4 Mbp in 2,250 contigs. This results in ~9% more of the genome in ~85% fewer contigs from ~80% less starting material at a fraction of the cost of previous Roche 454-based sequencing efforts. CONCLUSIONS: The result is the most complete filarial nematode assembly produced thus far and demonstrates the utility of single molecule sequencing on the Pacific Biosciences platform for genetically heterogeneous metazoan genomes.

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.

Comparison of long-read sequencing technologies in interrogating bacteria and fly genomes.

The newest generation of DNA sequencing technology is highlighted by the ability to generate sequence reads hundreds of kilobases in length. Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) have pioneered competitive long read platforms, with more recent work focused on improving sequencing throughput and per-base accuracy. We used whole-genome sequencing data produced by three PacBio protocols (Sequel II CLR, Sequel II HiFi, RS II) and two ONT protocols (Rapid Sequencing and Ligation Sequencing) to compare assemblies of the bacteria Escherichia coli and the fruit fly Drosophila ananassae. In both organisms tested, Sequel II assemblies had the highest consensus accuracy, even after accounting for differences in sequencing throughput. ONT and PacBio CLR had the longest reads sequenced compared to PacBio RS II and HiFi, and genome contiguity was highest when assembling these datasets. ONT Rapid Sequencing libraries had the fewest chimeric reads in addition to superior quantification of E. coli plasmids versus ligation-based libraries. The quality of assemblies can be enhanced by adopting hybrid approaches using Illumina libraries for bacterial genome assembly or polishing eukaryotic genome assemblies, and an ONT-Illumina hybrid approach would be more cost-effective for many users. Genome-wide DNA methylation could be detected using both technologies, however ONT libraries enabled the identification of a broader range of known E. coli methyltransferase recognition motifs in addition to undocumented D. ananassae motifs. The ideal choice of long read technology may depend on several factors including the question or hypothesis under examination. No single technology outperformed others in all metrics examined.

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.

Intratumor genetic heterogeneity in squamous cell carcinoma of the oral cavity.

BACKGROUND: We sought to evaluate intratumor heterogeneity in squamous cell carcinoma of the oral cavity (OCC) and specifically determine the effect of physical separation and histologic differentiation within the same tumor. METHODS: We performed whole exome sequencing on five biopsy sites-two from well-differentiated, two from poorly differentiated regions, and one from normal parenchyma-from five primary OCC specimens. RESULTS: We found high levels of intratumor heterogeneity and, in four primary tumors, identified only 0 to 2 identical mutations in all subsites. We found that the heterogeneity inversely correlated with physical separation and that pairs of well-differentiated samples were more similar to each other than analogous poorly differentiated specimens. Only TP53 mutations, but not other purported "driver mutations" in head and neck squamous cell carcinoma, were found in multiple biopsy sites. CONCLUSION: These data highlight the challenges to characterization of the mutational landscape of OCC with single site biopsy and have implications for personalized medicine.

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.

Complete Genome Sequences of Mycobacterium kansasii Strains Isolated from Rhesus Macaques.

Mycobacterium kansasii is a nontuberculous mycobacterium. It causes opportunistic infections with pulmonary and extrapulmonary manifestations. We report here the complete genome sequences of two M. kansasii strains isolated from rhesus macaques. We performed genome comparisons with human and environmental isolates of M. kansasii to assess the genomic diversity of this species.

Whole-Genome Sequences of Bacteremia Isolates of Bordetella holmesii.

Bordetella holmesii causes respiratory and invasive diseases in humans, but its pathogenesis remains poorly understood. We report here the genome sequences of seven bacteremia isolates of B. holmesii, including the type strain. Comparative analysis of these sequences may aid studies of B. holmesii biology and assist in the development of species-specific diagnostic strategies.

Efficient Enrichment of Bacterial mRNA from Host-Bacteria Total RNA Samples.

Despite numerous advances in genomics and bioinformatics, technological hurdles remain to examine host-microbe transcriptomics. Sometimes the transcriptome of either or both can be ascertained merely by generating more sequencing reads. However, many cases exist where bacterial mRNA needs to be enriched further to enable cost-effective sequencing of the pathogen or endosymbiont. While a suitable method is commercially available for mammalian samples of this type, development of such methods has languished for invertebrate samples. Furthermore, a common method across multiple taxa would facilitate comparisons between bacteria in invertebrate vectors and their vertebrate hosts. Here, a method is described to concurrently remove polyadenylated transcripts, prokaryotic rRNA, and eukaryotic rRNA, including those with low amounts of starting material (e.g. 100 ng). In a Wolbachia-Drosophila system, this bacterial mRNA enrichment yielded a 3-fold increase in Wolbachia mRNA abundance and a concomitant 3.3-fold increase in the percentage of transcripts detected. More specifically, 70% of the genome could be recovered by transcriptome sequencing compared to 21% in the total RNA. Sequencing of similar bacterial mRNA-enriched samples generated from Ehrlichia-infected canine cells covers 93% of the Ehrlichia genome, suggesting ubiquitous transcription across the entire Ehrlichia chaffeensis genome. This technique can potentially be used to enrich bacterial mRNA in many studies of host-microbe interactions.

Effect of salt and osmotic stress upon expression of the ethylene receptor ETR1 in Arabidopsis thaliana.

In hormone perception, varying the concentrations of hormone, receptor, or downstream signaling elements can modulate signal transduction. Previous research has demonstrated that ethylene biosynthesis in plants is regulated by abiotic factors. Here we report that exposure of Arabidopsis plants to NaCl reduced expression of the ethylene receptor ETR1. The change in gene expression was reflected at the protein level based on immunoblot analysis. Further analysis supports a general effect of osmotic stress upon the expression level of ETR1. The reduction in ETR1 levels should cause increased sensitivity of the plant to ethylene. These results suggest that plant responses to abiotic stress are modulated by changes in the expression level of ethylene receptors.

Genome sequences of rhinovirus a isolates from wisconsin pediatric respiratory studies.

Full-length or nearly full-length RNA genome sequences for 98 rhinovirus (RV) A isolates (from the Enterovirus genus of the Picornaviridae family), representing 43 different genotypes, were resolved as part of ongoing studies to define RV genetic diversity and its potential link to respiratory disease.

Genome sequences of rhinovirus B isolates from wisconsin pediatric respiratory studies.

Nearly full-length RNA genome sequences for 39 rhinovirus B isolates (RV-B), representing 13 different genotypes, were resolved as part of ongoing studies at the University of Wisconsin that attempt to link rhinovirus (RV) diversity and respiratory disease in infants.

Genome sequences of rhinovirus C isolates from wisconsin pediatric respiratory studies.

Human rhinovirus (RV) isolates from the RV-C species are recently discovered infectious agents that are closely linked to asthma and wheezing etiologies in infants. Clinical study samples collected at the University of Wisconsin-Madison describe 41 nearly complete genome sequences representing 21 RV-C genotypes.

Ligand-induced degradation of the ethylene receptor ETR2 through a proteasome-dependent pathway in Arabidopsis.

Protein degradation plays an important role in modulating ethylene signal transduction in plants. Here we show that the ethylene receptor ETR2 is one such target for degradation and that its degradation is dependent upon perception of the signaling ligand ethylene. The ETR2 protein is initially induced by ethylene treatment, consistent with an increase in transcript levels. At ethylene concentrations above 1 mul/liter, however, ETR2 protein levels subsequently decrease in a post-transcriptional fashion. Genetic and chemical approaches indicate that ethylene perception by the receptors initiates the reduction in ETR2 protein levels. The ethylene-induced decrease in ETR2 levels is not affected by cycloheximide, an inhibitor of protein biosynthesis, but is affected by proteasome inhibitors, indicating a role for the proteasome in ETR2 degradation. Ethylene-induced degradation still occurs in seedlings treated with brefeldin A, indicating that degradation of ETR2 does not require exit from its subcellular location at the endoplasmic reticulum. These data support a model in which ETR2 is degraded by a proteasome-dependent pathway in response to ethylene binding. Implications of this model for ethylene signaling are discussed.

Toward a comprehensive atlas of the physical interactome of Saccharomyces cerevisiae.

Defining protein complexes is critical to virtually all aspects of cell biology. Two recent affinity purification/mass spectrometry studies in Saccharomyces cerevisiae have vastly increased the available protein interaction data. The practical utility of such high throughput interaction sets, however, is substantially decreased by the presence of false positives. Here we created a novel probabilistic metric that takes advantage of the high density of these data, including both the presence and absence of individual associations, to provide a measure of the relative confidence of each potential protein-protein interaction. This analysis largely overcomes the noise inherent in high throughput immunoprecipitation experiments. For example, of the 12,122 binary interactions in the general repository of interaction data (BioGRID) derived from these two studies, we marked 7504 as being of substantially lower confidence. Additionally, applying our metric and a stringent cutoff we identified a set of 9074 interactions (including 4456 that were not among the 12,122 interactions) with accuracy comparable to that of conventional small scale methodologies. Finally we organized proteins into coherent multisubunit complexes using hierarchical clustering. This work thus provides a highly accurate physical interaction map of yeast in a format that is readily accessible to the biological community.

Effect of ethylene pathway mutations upon expression of the ethylene receptor ETR1 from Arabidopsis.

The ethylene receptor family of Arabidopsis consists of five members, one of these being ETR1. The effect of ethylene pathway mutations upon expression of ETR1 was examined. For this purpose, ETR1 levels were quantified in mutant backgrounds containing receptor loss-of-function mutations, ethylene-insensitive mutations, and constitutive ethylene response mutations. Ethylene-insensitive mutations of ETR1 resulted in a posttranscriptional increase in levels of the mutant receptor. Treatment of seedlings with silver, which leads to ethylene insensitivity, also resulted in an increase in levels of ETR1. Loss-of-function mutations of ETR1 resulted in both transcriptional and posttranscriptional changes in levels of the receptor. Most other ethylene pathway mutations, including a newly isolated T-DNA insertion mutation in the gene encoding the ethylene receptor ERS1, had relatively minor effects upon the expression of ETR1. Our results indicate that mutations in ETR1 can affect expression at the posttranscriptional level, and suggest that these posttranscriptional changes may contribute to the phenotypes observed in the mutants. Our results also refine the model on how mutations in ethylene receptors are able to confer dominant ethylene insensitivity upon plants.

Localization of the Raf-like kinase CTR1 to the endoplasmic reticulum of Arabidopsis through participation in ethylene receptor signaling complexes.

The plant hormone ethylene is perceived by a five-member family of receptors related to the bacterial histidine kinases. The Raf-like kinase CTR1 functions downstream of the ethylene receptors as a negative regulator of ethylene signal transduction. CTR1 is shown here to be associated with membranes of the endoplasmic reticulum in Arabidopsis as a result of its interactions with ethylene receptors. Membrane association of CTR1 is reduced by mutations that eliminate ethylene receptors and by a mutation in CTR1 that reduces its ability to bind to the ethylene receptor ETR1. Direct evidence that CTR1 is part of an ethylene receptor signaling complex was obtained by co-purification of the ethylene receptor ETR1 with a tagged version of CTR1 from an Arabidopsis membrane extract. The histidine kinase activity of ETR1 is not required for its association with CTR1, based on co-purification of tagged ETR1 mutants and CTR1 after expression in a transgenic yeast system. These data demonstrate that CTR1 is part of an ethylene receptor signaling complex in Arabidopsis and support a model in which localization of CTR1 to the endoplasmic reticulum is necessary for its function. Additional data that demonstrate a post-transcriptional effect of ethylene upon the expression of CTR1 suggest that production of ethylene receptor signaling complexes may be coordinately regulated.