The Addition of Transcriptomics to the Bead-Enabled Accelerated Monophasic Multi-Omics Method: A Step toward Universal Sample Preparation.

Omics-based measurements enable the study of biomolecules in a high-throughput fashion, leading to the characterization and quantification of biological systems. Multi-omics methods aim to incorporate several omics measurements for a more holistic approach, which is crucial for advancing our understanding of the diversity and redundancy of biological systems. Current multi-omics sample preparation methods have achieved proteomics, lipidomics, and metabolomics from individual samples; however, the bioinformatic tools currently available for interpreting data generated from these omics are limited. Alternately, transcriptomics has a wide arsenal of available bioinformatic tools offering intensive sample characterization but has yet to be incorporated into a unified, multi-omics sample preparation technique. Herein we describe the modified bead-enabled accelerated monophasic multi-omics (mBAMM) method, which incorporates RNA extraction for transcriptomics analysis. mBAMM was shown to enable RNA-seq without compromising the isolation of biomolecules for proteomics, lipidomics, and metabolomics. This methodology greatly improves sample characterization and represents a major innovation toward cohesive insights into biological systems.

Nonclassical autophagy activation pathways are essential for production of infectious Influenza A virus in vitro.

Autophagy is a critical mechanism deployed by eukaryotic cells in response to stress, including viral infection, to boost the innate antimicrobial responses. However, an increasing number of pathogens hijack the autophagic machinery to facilitate their own replication. Influenza A virus (IAV), responsible for several global pandemics, has an intricate dependence on autophagy for successful replication in mammalian cells. To elucidate key chokepoints in the host stress responses facilitating IAV replication, we constructed a meta-transcriptome of IAV and host gene expression dynamics during early (1-3 hpi), mid (4-6 hpi), and late (8-12 hpi) stages of the viral replication cycle at two multiplicities of infection (MOI): 1 and 5. We supplemented the global transcriptome study with phosphoproteomic analysis of stress-activated protein kinase (SAPK/JNK) signaling in lung carcinoma (predominantly used as an in vitro model of IAV replication) and normal human bronchial epithelial cells. We report significant differences in the activation profiles of autophagy regulating genes upon IAV infection at the two MOI as well as divergent dependence on ULK1 signaling within the normal and cancer cells. Regardless of the cell model, JNK-Thr187 signaling was crucial for the production of infectious viral particles.

Remedial Treatment of Corroded Iron Objects by Environmental Aeromonas Isolates.

Using bacteria to transform reactive corrosion products into stable compounds represents an alternative to traditional methods employed in iron conservation. Two environmental Aeromonas strains (CA23 and CU5) were used to transform ferric iron corrosion products (goethite and lepidocrocite) into stable ferrous iron-bearing minerals (vivianite and siderite). A genomic and transcriptomic approach was used to analyze the metabolic traits of these strains and to evaluate their pathogenic potential. Although genes involved in solid-phase iron reduction were identified, key genes present in other environmental iron-reducing species are missing from the genome of CU5. Several pathogenicity factors were identified in the genomes of both strains, but none of these was expressed under iron reduction conditions. Additional in vivo tests showed hemolytic and cytotoxic activities for strain CA23 but not for strain CU5. Both strains were easily inactivated using ethanol and heat. Nonetheless, given a lesser potential for a pathogenic lifestyle, CU5 is the most promising candidate for the development of a bio-based iron conservation method stabilizing iron corrosion. Based on all the results, a prototype treatment was established using archaeological items. On those, the conversion of reactive corrosion products and the formation of a homogenous layer of biogenic iron minerals were achieved. This study shows how naturally occurring microorganisms and their metabolic capabilities can be used to develop bio-inspired solutions to the problem of metal corrosion.IMPORTANCE Microbiology can greatly help in the quest for a sustainable solution to the problem of iron corrosion, which causes important economic losses in a wide range of fields, including the protection of cultural heritage and building materials. Using bacteria to transform reactive and unstable corrosion products into more-stable compounds represents a promising approach. The overall aim of this study was to develop a method for the conservation and restoration of corroded iron items, starting from the isolation of iron-reducing bacteria from natural environments. This resulted in the identification of a suitable candidate (Aeromonas sp. strain CU5) that mediates the formation of desirable minerals at the surfaces of the objects. This led to the proof of concept of an application method on real objects.

Nearly finished genomes produced using gel microdroplet culturing reveal substantial intraspecies genomic diversity within the human microbiome.

The majority of microbial genomic diversity remains unexplored. This is largely due to our inability to culture most microorganisms in isolation, which is a prerequisite for traditional genome sequencing. Single-cell sequencing has allowed researchers to circumvent this limitation. DNA is amplified directly from a single cell using the whole-genome amplification technique of multiple displacement amplification (MDA). However, MDA from a single chromosome copy suffers from amplification bias and a large loss of specificity from even very small amounts of DNA contamination, which makes assembling a genome difficult and completely finishing a genome impossible except in extraordinary circumstances. Gel microdrop cultivation allows culturing of a diverse microbial community and provides hundreds to thousands of genetically identical cells as input for an MDA reaction. We demonstrate the utility of this approach by comparing sequencing results of gel microdroplets and single cells following MDA. Bias is reduced in the MDA reaction and genome sequencing, and assembly is greatly improved when using gel microdroplets. We acquired multiple near-complete genomes for two bacterial species from human oral and stool microbiome samples. A significant amount of genome diversity, including single nucleotide polymorphisms and genome recombination, is discovered. Gel microdroplets offer a powerful and high-throughput technology for assembling whole genomes from complex samples and for probing the pan-genome of naturally occurring populations.

A multiple-alignment based primer design algorithm for genetically highly variable DNA targets.

BACKGROUND: Primer design for highly variable DNA sequences is difficult, and experimental success requires attention to many interacting constraints. The advent of next-generation sequencing methods allows the investigation of rare variants otherwise hidden deep in large populations, but requires attention to population diversity and primer localization in relatively conserved regions, in addition to recognized constraints typically considered in primer design. RESULTS: Design constraints include degenerate sites to maximize population coverage, matching of melting temperatures, optimizing de novo sequence length, finding optimal bio-barcodes to allow efficient downstream analyses, and minimizing risk of dimerization. To facilitate primer design addressing these and other constraints, we created a novel computer program (PrimerDesign) that automates this complex procedure. We show its powers and limitations and give examples of successful designs for the analysis of HIV-1 populations. CONCLUSIONS: PrimerDesign is useful for researchers who want to design DNA primers and probes for analyzing highly variable DNA populations. It can be used to design primers for PCR, RT-PCR, Sanger sequencing, next-generation sequencing, and other experimental protocols targeting highly variable DNA samples.

Epitope-specific CD8+ T lymphocytes cross-recognize mutant simian immunodeficiency virus (SIV) sequences but fail to contain very early evolution and eventual fixation of epitope escape mutations during SIV infection.

Human immunodeficiency virus (HIV) and simian immunodeficiency virus (SIV) evade containment by CD8(+) T lymphocytes through focused epitope mutations. However, because of limitations in the numbers of viral sequences that can be sampled, traditional sequencing technologies have not provided a true representation of the plasticity of these viruses or the intensity of CD8(+) T lymphocyte-mediated selection pressure. Moreover, the strategy by which CD8(+) T lymphocytes contain evolving viral quasispecies has not been characterized fully. In the present study we have employed ultradeep 454 pyrosequencing of virus and simultaneous staining of CD8(+) T lymphocytes with multiple tetramers in the SIV/rhesus monkey model to explore the coevolution of virus and the cellular immune response during primary infection. We demonstrated that cytotoxic T lymphocyte (CTL)-mediated selection pressure on the infecting virus was manifested by epitope mutations as early as 21 days following infection. We also showed that CD8(+) T lymphocytes cross-recognized wild-type and mutant epitopes and that these cross-reactive cell populations were present at a time when mutant forms of virus were present at frequencies of as low as 1 in 22,000 sequenced clones. Surprisingly, these cross-reactive cells became enriched in the epitope-specific CD8(+) T lymphocyte population as viruses with mutant epitope sequences largely replaced those with epitope sequences of the transmitted virus. These studies demonstrate that mutant epitope-specific CD8(+) T lymphocytes that are present at a time when viral mutant epitope sequences are detected at extremely low frequencies fail to contain the later accumulation and fixation of the mutant epitope sequences in the viral quasispecies.

Comparing variability in diagnosis of upper respiratory tract infections in patients using syndromic, next generation sequencing, and PCR-based methods.

Early and accurate diagnosis of respiratory pathogens and associated outbreaks can allow for the control of spread, epidemiological modeling, targeted treatment, and decision making-as is evident with the current COVID-19 pandemic. Many respiratory infections share common symptoms, making them difficult to diagnose using only syndromic presentation. Yet, with delays in getting reference laboratory tests and limited availability and poor sensitivity of point-of-care tests, syndromic diagnosis is the most-relied upon method in clinical practice today. Here, we examine the variability in diagnostic identification of respiratory infections during the annual infection cycle in northern New Mexico, by comparing syndromic diagnostics with polymerase chain reaction (PCR) and sequencing-based methods, with the goal of assessing gaps in our current ability to identify respiratory pathogens. Of 97 individuals that presented with symptoms of respiratory infection, only 23 were positive for at least one RNA virus, as confirmed by sequencing. Whereas influenza virus (n = 7) was expected during this infection cycle, we also observed coronavirus (n = 7), respiratory syncytial virus (n = 8), parainfluenza virus (n = 4), and human metapneumovirus (n = 1) in individuals with respiratory infection symptoms. Four patients were coinfected with two viruses. In 21 individuals that tested positive using PCR, RNA sequencing completely matched in only 12 (57%) of these individuals. Few individuals (37.1%) were diagnosed to have an upper respiratory tract infection or viral syndrome by syndromic diagnostics, and the type of virus could only be distinguished in one patient. Thus, current syndromic diagnostic approaches fail to accurately identify respiratory pathogens associated with infection and are not suited to capture emerging threats in an accurate fashion. We conclude there is a critical and urgent need for layered agnostic diagnostics to track known and unknown pathogens at the point of care to control future outbreaks.

Correction: California condor microbiomes: Bacterial variety and functional properties in captive-bred individuals.

[This corrects the article DOI: 10.1371/journal.pone.0225858.].

A Gene Cluster That Encodes Histone Deacetylase Inhibitors Contributes to Bacterial Persistence and Antibiotic Tolerance in Burkholderia thailandensis.

Persister cells are genetically identical variants in a bacterial population that have phenotypically modified their physiology to survive environmental stress. In bacterial pathogens, persisters are able to survive antibiotic treatment and reinfect patients in a frustrating cycle of chronic infection. To better define core persistence mechanisms for therapeutics development, we performed transcriptomics analyses of Burkholderia thailandensis populations enriched for persisters via three methods: flow sorting for low proton motive force, meropenem treatment, and culture aging. Although the three persister-enriched populations generally displayed divergent gene expression profiles that reflect the multimechanistic nature of stress adaptations, there were several common gene pathways activated in two or all three populations. These include polyketide and nonribosomal peptide synthesis, Clp proteases, mobile elements, enzymes involved in lipid metabolism, and ATP-binding cassette (ABC) transporter systems. In particular, identification of genes that encode polyketide synthases (PKSs) and fatty acid catabolism factors indicates that generation of secondary metabolites, natural products, and complex lipids could be part of the metabolic program that governs the persistence state. We also found that loss-of-function mutations in the PKS-encoding gene locus BTH_I2366, which plays a role in biosynthesis of histone deacetylase (HDAC) inhibitors, resulted in increased sensitivity to antibiotics targeting DNA replication. Furthermore, treatment of multiple bacterial pathogens with a fatty acid synthesis inhibitor, CP-640186, potentiated the efficacy of meropenem against the persister populations. Altogether, our results suggest that bacterial persisters may exhibit an outwardly dormant physiology but maintain active metabolic processes that are required to maintain persistence.IMPORTANCE The discovery of antibiotics such as penicillin and streptomycin marked a historic milestone in the 1940s and heralded a new era of antimicrobial therapy as the modern standard for medical treatment. Yet, even in those early days of discovery, it was noted that a small subset of cells (∼1 in 105) survived antibiotic treatment and continued to persist, leading to recurrence of chronic infection. These persisters are phenotypic variants that have modified their physiology to survive environmental stress. In this study, we have performed three transcriptomic screens to identify persistence genes that are common between three different stressor conditions. In particular, we identified genes that function in the synthesis of secondary metabolites, small molecules, and complex lipids, which are likely required to maintain the persistence state. Targeting universal persistence genes can lead to the development of clinically relevant antipersistence therapeutics for infectious disease management.

Genome Sequence of a Staphylococcus xylosus Clinical Isolate, Strain SMA0341-04 (UGA5), from Siaya County Referral Hospital in Siaya, Kenya.

We report here the genome sequence of a Staphylococcus xylosus clinical isolate, strain SMA0341-04 (UGA5), which contains one chromosome and at least one plasmid. Notably, strain SMA0341-04 (UGA5) contains the tetracycline efflux major facilitator superfamily (MFS) transporter (tetK) gene.

Genome Sequence of Staphylococcus pettenkoferi Strain SMA0010-04 (UGA20), a Clinical Isolate from Siaya County Referral Hospital in Siaya, Kenya.

Here, we report the sequence of a Staphylococcus pettenkoferi clinical isolate, strain SMA0010-04 (UGA20), which contains the PC1 beta-lactamase (blaZ) gene.

Genome Sequences of a Staphylococcus aureus Clinical Isolate, Strain SMA0034-04 (UGA22), from Siaya County Referral Hospital in Siaya, Kenya.

Here, we report the genome sequences of a Staphylococcus aureus clinical isolate, strain SMA0034-04 (UGA22), which contains one chromosome and one plasmid. We also reveal that isolate SMA0034-04 (UGA22) contains loci in the genome that encode multiple exotoxins.

Draft Genome Sequences of Two Staphylococcus warneri Clinical Isolates, Strains SMA0023-04 (UGA3) and SMA0670-05 (UGA28), from Siaya County Referral Hospital, Siaya, Kenya.

We report the complete draft genome sequences of two Staphylococcus warneri clinical isolates, strains SMA0023-04 (UGA3) and SMA0670-05 (UGA28), each of which contains one chromosome and at least one plasmid. Isolate SMA0023-04 (UGA3) contains tetracycline efflux major facilitator superfamily (MFS) transporter (tetK), macrolide resistance (msrC and mphC), and beta-lactamase (blaZ) genes on its plasmids.

California condor microbiomes: Bacterial variety and functional properties in captive-bred individuals.

Around the world, scavenging birds such as vultures and condors have been experiencing drastic population declines. Scavenging birds have a distinct digestive process to deal with higher amounts of bacteria in their primary diet of carcasses in varying levels of decay. These observations motivate us to present an analysis of captive and healthy California condor (Gymnogyps californianus) microbiomes to characterize a population raised together under similar conditions. Shotgun metagenomic DNA sequences were analyzed from fecal and cloacal samples of captive birds. Classification of shotgun DNA sequence data with peptide signatures using the Sequedex package provided both phylogenetic and functional profiles, as well as individually annotated reads for targeted confirmatory analysis. We observed bacterial species previously associated with birds and gut microbiomes, including both virulent and opportunistic pathogens such as Clostridium perfringens, Propionibacterium acnes, Shigella flexneri, and Fusobacterium mortiferum, common flora such as Lactobacillus johnsonii, Lactobacillus ruminus, and Bacteroides vulgatus, and mucosal microbes such as Delftia acidovorans, Stenotrophomonas maltophilia, and Corynebacterium falsnii. Classification using shotgun metagenomic reads from phylogenetic marker genes was consistent with, and more specific than, analysis based on 16S rDNA data. Classification of samples based on either phylogenetic or functional profiles of genomic fragments differentiated three types of samples: fecal, mature cloacal and immature cloacal, with immature birds having approximately 40% higher diversity of microbes.

A metagenomic viral discovery approach identifies potential zoonotic and novel mammalian viruses in Neoromicia bats within South Africa.

Species within the Neoromicia bat genus are abundant and widely distributed in Africa. It is common for these insectivorous bats to roost in anthropogenic structures in urban regions. Additionally, Neoromicia capensis have previously been identified as potential hosts for Middle East respiratory syndrome (MERS)-related coronaviruses. This study aimed to ascertain the gastrointestinal virome of these bats, as viruses excreted in fecal material or which may be replicating in rectal or intestinal tissues have the greatest opportunities of coming into contact with other hosts. Samples were collected in five regions of South Africa over eight years. Initial virome composition was determined by viral metagenomic sequencing by pooling samples and enriching for viral particles. Libraries were sequenced on the Illumina MiSeq and NextSeq500 platforms, producing a combined 37 million reads. Bioinformatics analysis of the high throughput sequencing data detected the full genome of a novel species of the Circoviridae family, and also identified sequence data from the Adenoviridae, Coronaviridae, Herpesviridae, Parvoviridae, Papillomaviridae, Phenuiviridae, and Picornaviridae families. Metagenomic sequencing data was insufficient to determine the viral diversity of certain families due to the fragmented coverage of genomes and lack of suitable sequencing depth, as some viruses were detected from the analysis of reads-data only. Follow up conventional PCR assays targeting conserved gene regions for the Adenoviridae, Coronaviridae, and Herpesviridae families were used to confirm metagenomic data and generate additional sequences to determine genetic diversity. The complete coding genome of a MERS-related coronavirus was recovered with additional amplicon sequencing on the MiSeq platform. The new genome shared 97.2% overall nucleotide identity to a previous Neoromicia-associated MERS-related virus, also from South Africa. Conventional PCR analysis detected diverse adenovirus and herpesvirus sequences that were widespread throughout Neoromicia populations in South Africa. Furthermore, similar adenovirus sequences were detected within these populations throughout several years. With the exception of the coronaviruses, the study represents the first report of sequence data from several viral families within a Southern African insectivorous bat genus; highlighting the need for continued investigations in this regard.

Molecular Dissection of Induced Platinum Resistance through Functional and Gene Expression Analysis in a Cell Culture Model of Bladder Cancer.

We report herein the development, functional and molecular characterization of an isogenic, paired bladder cancer cell culture model system for studying platinum drug resistance. The 5637 human bladder cancer cell line was cultured over ten months with stepwise increases in oxaliplatin concentration to generate a drug resistant 5637R sub cell line. The MTT assay was used to measure the cytotoxicity of several bladder cancer drugs. Liquid scintillation counting allowed quantification of cellular drug uptake and efflux of radiolabeled oxaliplatin and carboplatin. The impact of intracellular drug inactivation was assessed by chemical modulation of glutathione levels. Oxaliplatin- and carboplatin-DNA adduct formation and repair was measured using accelerator mass spectrometry. Resistance factors including apoptosis, growth factor signaling and others were assessed with RNAseq of both cell lines and included confirmation of selected transcripts by RT-PCR. Oxaliplatin, carboplatin, cisplatin and gemcitabine were significantly less cytotoxic to 5637R cells compared to the 5637 cells. In contrast, doxorubicin, methotrexate and vinblastine had no cell line dependent difference in cytotoxicity. Upon exposure to therapeutically relevant doses of oxaliplatin, 5637R cells had lower drug-DNA adduct levels than 5637 cells. This difference was partially accounted for by pre-DNA damage mechanisms such as drug uptake and intracellular inactivation by glutathione, as well as faster oxaliplatin-DNA adduct repair. In contrast, both cell lines had no significant differences in carboplatin cell uptake, efflux and drug-DNA adduct formation and repair, suggesting distinct resistance mechanisms for these two closely related drugs. The functional studies were augmented by RNAseq analysis, which demonstrated a significant change in expression of 83 transcripts, including 50 known genes and 22 novel transcripts. Most of the transcripts were not previously associated with bladder cancer chemoresistance. This model system and the associated phenotypic and genotypic data has the potential to identify some novel details of resistance mechanisms of clinical importance to bladder cancer.

Forest floor community metatranscriptomes identify fungal and bacterial responses to N deposition in two maple forests.

Anthropogenic N deposition alters patterns of C and N cycling in temperate forests, where forest floor litter decomposition is a key process mediated by a diverse community of bacteria and fungi. To track forest floor decomposer activity we generated metatranscriptomes that simultaneously surveyed the actively expressed bacterial and eukaryote genes in the forest floor, to compare the impact of N deposition on the decomposers in two natural maple forests in Michigan, USA, where replicate field plots had been amended with N for 16 years. Site and N amendment responses were compared using about 74,000 carbohydrate active enzyme transcript sequences (CAZymes) in each metatranscriptome. Parallel ribosomal RNA (rRNA) surveys of bacterial and fungal biomass and taxonomic composition showed no significant differences in either biomass or OTU richness between the two sites or in response to N. Site and N amendment were not significant variables defining bacterial taxonomic composition, but they were significant for fungal community composition, explaining 17 and 14% of the variability, respectively. The relative abundance of expressed bacterial and fungal CAZymes changed significantly with N amendment in one of the forests, and N-response trends were also identified in the second forest. Although the two ambient forests were similar in community biomass, taxonomic structure and active CAZyme profile, the shifts in active CAZyme profiles in response to N-amendment differed between the sites. One site responded with an over-expression of bacterial CAZymes, and the other site responded with an over-expression of both fungal and different bacterial CAZymes. Both sites showed reduced representation of fungal lignocellulose degrading enzymes in N-amendment plots. The metatranscriptome approach provided a holistic assessment of eukaryote and bacterial gene expression and is applicable to other systems where eukaryotes and bacteria interact.

Genome Sequence of a Chromium-Reducing Strain, Bacillus cereus S612.

We report here the genome sequence of an effective chromium-reducing bacterium, Bacillus cereus strain S612. The size of the draft genome sequence is approximately 5.4 Mb, with a G+C content of 35%, and it is predicted to contain 5,450 protein-coding genes.

Genome Sequence of Aeribacillus pallidus Strain GS3372, an Endospore-Forming Bacterium Isolated in a Deep Geothermal Reservoir.

The genome of strain GS3372 is the first publicly available strain of Aeribacillus pallidus. This endospore-forming thermophilic strain was isolated from a deep geothermal reservoir. The availability of this genome can contribute to the clarification of the taxonomy of the closely related Anoxybacillus, Geobacillus, and Aeribacillus genera.

Genome Sequence of Bacillus alveayuensis Strain 24KAM51, a Halotolerant Thermophile Isolated from a Hydrothermal Vent.

Bacillus alveayuensis strain 24KAM51 was isolated from a marine hydrothermal vent in Milos, Greece. Its genome depicts interesting features of halotolerance and resistance to heavy metals.