Introducing the Bacterial and Viral Bioinformatics Resource Center (BV-BRC): a resource combining PATRIC, IRD and ViPR.

The National Institute of Allergy and Infectious Diseases (NIAID) established the Bioinformatics Resource Center (BRC) program to assist researchers with analyzing the growing body of genome sequence and other omics-related data. In this report, we describe the merger of the PAThosystems Resource Integration Center (PATRIC), the Influenza Research Database (IRD) and the Virus Pathogen Database and Analysis Resource (ViPR) BRCs to form the Bacterial and Viral Bioinformatics Resource Center (BV-BRC) https://www.bv-brc.org/. The combined BV-BRC leverages the functionality of the bacterial and viral resources to provide a unified data model, enhanced web-based visualization and analysis tools, bioinformatics services, and a powerful suite of command line tools that benefit the bacterial and viral research communities.

Phenotype-Based Threat Assessment.

Bacterial pathogen identification, which is critical for human health, has historically relied on culturing organisms from clinical specimens. More recently, the application of machine learning (ML) to whole-genome sequences (WGSs) has facilitated pathogen identification. However, relying solely on genetic information to identify emerging or new pathogens is fundamentally constrained, especially if novel virulence factors exist. In addition, even WGSs with ML pipelines are unable to discern phenotypes associated with cryptic genetic loci linked to virulence. Here, we set out to determine if ML using phenotypic hallmarks of pathogenesis could assess potential pathogenic threat without using any sequence-based analysis. This approach successfully classified potential pathogenetic threat associated with previously machine-observed and unobserved bacteria with 99% and 85% accuracy, respectively. This work establishes a phenotype-based pipeline for potential pathogenic threat assessment, which we term PathEngine, and offers strategies for the identification of bacterial pathogens.

The PATRIC Bioinformatics Resource Center: expanding data and analysis capabilities.

The PathoSystems Resource Integration Center (PATRIC) is the bacterial Bioinformatics Resource Center funded by the National Institute of Allergy and Infectious Diseases (https://www.patricbrc.org). PATRIC supports bioinformatic analyses of all bacteria with a special emphasis on pathogens, offering a rich comparative analysis environment that provides users with access to over 250 000 uniformly annotated and publicly available genomes with curated metadata. PATRIC offers web-based visualization and comparative analysis tools, a private workspace in which users can analyze their own data in the context of the public collections, services that streamline complex bioinformatic workflows and command-line tools for bulk data analysis. Over the past several years, as genomic and other omics-related experiments have become more cost-effective and widespread, we have observed considerable growth in the usage of and demand for easy-to-use, publicly available bioinformatic tools and services. Here we report the recent updates to the PATRIC resource, including new web-based comparative analysis tools, eight new services and the release of a command-line interface to access, query and analyze data.

The Role of luxS in Histophilus somni Virulence and Biofilm Formation.

S-Ribosylhomocysteinase (LuxS) is required for the synthesis of the autoinducer-2 (AI-2) quorum-sensing signaling molecule in many Gram-negative bacteria. The bovine (and ovine) opportunistic pathogen Histophilus somni contains luxS and forms a biofilm containing an exopolysaccharide (EPS) in the matrix. Since biofilm formation is regulated by quorum sensing in many bacteria, the roles of luxS in H. somni virulence and biofilm formation were investigated. Although culture supernatants from H. somni were ineffective at inducing bioluminescence in the Vibrio harveyi reporter strain BB170, H. somniluxS complemented the biosynthesis of AI-2 in the luxS-deficient Escherichia coli strain DH5α. H. somni strain 2336 luxS was inactivated by transposon mutagenesis. RNA expression profiles revealed that many genes were significantly differentially expressed in the luxS mutant compared to that in the wild-type, whether the bacteria were grown planktonically or in a biofilm. Furthermore, the luxS mutant had a truncated and asialylated lipooligosaccharide (LOS) and was substantially more serum sensitive than the wild-type. Not surprisingly, the luxS mutant was attenuated in a mouse model for H. somni virulence, and some of the altered phenotypes were partially restored after the mutation was complemented with a functional luxS However, no major differences were observed between the wild-type and the luxS mutant in regard to outer membrane protein profiles, biofilm formation, EPS production, or intracellular survival. These results indicate that luxS plays a role in H. somni virulence in the context of LOS biosynthesis but not biofilm formation or other phenotypic properties examined.

Phylogenomic analysis of Haemophilus parasuis and proposed reclassification to Glaesserella parasuis, gen. nov., comb. nov.

The Gram-negative bacterium Haemophilus parasuis is the etiologic agent of Glässer's disease in pigs, and causes significant economic losses to the swine industry. This bacterium has been classified as a member of the family Pasteurellaceae in the genus Haemophilus, but phylogenetic relatedness has not been adequately examined to support this genus classification. Phenotypically, all 38 strains of H. parasuis tested were positive for catalase activity, oxidase activity, V-factor requirement, and acid formation from maltose and d-galactose without gas. All strains were negative for X-factor requirement, formation of indole from tryptophan, urease, l-arabinose, and α-glucosidase activity. To determine whether H. parasuis belongs to one of the current Pasteurellaceae genera 40 H. parasuis genomes, plus those of representative Pasteurellaceae, were subjected to phylogenetic analysis of concatenated, multi-protein alignments. Sequence variation at 16S rRNA and rpoB loci allowed the 15 reference serovars of H. parasuis to be integrated into the whole-genome tree. The phylogenetic analysis showed H. parasuis to be a distinct and tight clade whose sister taxon is the genus Bibersteinia. Within H. parasuis two clades were identified with individual serovars distributed between the two. As a result, H. parasuis was confirmed as a member of the family Pasteurellaceae, but was distinct from other genera in this family. Therefore, we propose the name Glaesserella parasuis, gen. nov., comb. nov. for bacterial strains currently classified as H. parasuis. The reference strain of this species is ATCC 19417 (1374)T, NCTC 4557T, DSM 21448T, CCUG 3712T.

Comparative phylogenomics and evolution of the Brucellae reveal a path to virulence.

Brucella species include important zoonotic pathogens that have a substantial impact on both agriculture and human health throughout the world. Brucellae are thought of as "stealth pathogens" that escape recognition by the host innate immune response, modulate the acquired immune response, and evade intracellular destruction. We analyzed the genome sequences of members of the family Brucellaceae to assess its evolutionary history from likely free-living soil-based progenitors into highly successful intracellular pathogens. Phylogenetic analysis split the genus into two groups: recently identified and early-dividing "atypical" strains and a highly conserved "classical" core clade containing the major pathogenic species. Lateral gene transfer events brought unique genomic regions into Brucella that differentiated them from Ochrobactrum and allowed the stepwise acquisition of virulence factors that include a type IV secretion system, a perosamine-based O antigen, and systems for sequestering metal ions that are absent in progenitors. Subsequent radiation within the core Brucella resulted in lineages that appear to have evolved within their preferred mammalian hosts, restricting their virulence to become stealth pathogens capable of causing long-term chronic infections.

Characterization and RNA-seq analysis of underperformer, an activation-tagged potato mutant.

The potato cv. Bintje and a Bintje activation-tagged mutant, underperformer (up) were compared. Mutant up plants grown in vitro were dwarf, with abundant axillary shoot growth, greater tuber yield, altered tuber traits and early senescence compared to wild type. Under in vivo conditions, the dwarf and early senescence phenotypes of the mutant remained, but the up plants exhibited a lower tuber yield and fewer axillary shoots compared to wild type. Southern blot analyses indicated a single T-DNA insertion in the mutant, located on chromosome 10. Initial PCR-based gene expression studies indicated transcriptional activation/repression of several genes in the mutant flanking the insertion. The gene immediately flanking the right border of the T-DNA insertion, which encoded an uncharacterized Broad complex, Tramtrac, Bric-a-brac; also known as Pox virus and Zinc finger (BTB/POZ) domain-containing protein (StBTB/POZ1) containing an Armadillo repeat region, was up-regulated in the mutant. Global gene expression comparisons between Bintje and up using RNA-seq on leaves from 60 day-old plants revealed a dataset of over 1,600 differentially expressed genes. Gene expression analyses suggested a variety of biological processes and pathways were modified in the mutant, including carbohydrate and lipid metabolism, cell division and cell cycle activity, biotic and abiotic stress responses, and proteolysis.

An active ac/ds transposon system for activation tagging in tomato cultivar m82 using clonal propagation.

Tomato (Solanum lycopersicum) is a model organism for Solanaceae in both molecular and agronomic research. This project utilized Agrobacterium tumefaciens transformation and the transposon-tagging construct Activator (Ac)/Dissociator (Ds)-ATag-Bar_gosGFP to produce activation-tagged and knockout mutants in the processing tomato cultivar M82. The construct carried hygromycin resistance (hyg), green fluorescent protein (GFP), and the transposase (TPase) of maize (Zea mays) Activator major transcript X054214.1 on the stable Ac element, along with a 35S enhancer tetramer and glufosinate herbicide resistance (BAR) on the mobile Ds-ATag element. An in vitro propagation strategy was used to produce a population of 25 T0 plants from a single transformed plant regenerated in tissue culture. A T1 population of 11,000 selfed and cv M82 backcrossed progeny was produced from the functional T0 line. This population was screened using glufosinate herbicide, hygromycin leaf painting, and multiplex polymerase chain reaction (PCR). Insertion sites of transposed Ds-ATag elements were identified through thermal asymmetric interlaced PCR, and resulting product sequences were aligned to the recently published tomato genome. A population of 509 independent, Ds-only transposant lines spanning all 12 tomato chromosomes has been developed. Insertion site analysis demonstrated that more than 80% of these lines harbored Ds insertions conducive to activation tagging. The capacity of the Ds-ATag element to alter transcription was verified by quantitative real-time reverse transcription-PCR in two mutant lines. The transposon-tagged lines have been immortalized in seed stocks and can be accessed through an online database, providing a unique resource for tomato breeding and analysis of gene function in the background of a commercial tomato cultivar.

Comparative Genomic Analysis of Bacterial Data in BV-BRC: An Example Exploring Antimicrobial Resistance.

As genomic and related data continue to expand, research biologists are often hampered by the computational hurdles required to analyze their data. The National Institute of Allergy and Infectious Diseases (NIAID) established the Bioinformatics Resource Centers (BRC) to assist researchers with their analysis of genome sequence and other omics-related data. Recently, the PAThosystems Resource Integration Center (PATRIC), the Influenza Research Database (IRD), and the Virus Pathogen Database and Analysis Resource (ViPR) BRCs merged to form the Bacterial and Viral Bioinformatics Resource Center (BV-BRC) at https://www.bv-brc.org/ . The combined BV-BRC leverages the functionality of the original resources for bacterial and viral research communities with a unified data model, enhanced web-based visualization and analysis tools, and bioinformatics services. Here we demonstrate how antimicrobial resistance data can be analyzed in the new resource.

An automated system for interrogating the evolution of microbial endosymbiosis.

Inter-kingdom endosymbiotic interactions between bacteria and eukaryotic cells are critical to human health and disease. However, the molecular mechanisms that drive the emergence of endosymbiosis remain obscure. Here, we describe the development of a microfluidic system, named SEER (S̲ystem for the E̲volution of E̲ndosymbiotic R̲elationships), that automates the evolutionary selection of bacteria with enhanced intracellular survival and persistence within host cells, hallmarks of endosymbiosis. Using this system, we show that a laboratory strain of Escherichia coli that initially possessed limited abilities to survive within host cells, when subjected to SEER selection, rapidly evolved to display a 55-fold enhancement in intracellular survival. Notably, molecular dissection of the evolved strains revealed that a single-point mutation in a flexible loop of CpxR, a gene regulator that controls bacterial stress responses, substantially contributed to this intracellular survival. Taken together, these results establish SEER as the first microfluidic system for investigating the evolution of endosymbiosis, show the importance of CpxR in endosymbiosis, and set the stage for evolving bespoke inter-kingdom endosymbiotic systems with novel or emergent properties.

Transposon tagging in diploid strawberry.

Fragaria vesca was transformed with a transposon tagging construct harbouring amino terminally deleted maize transposase and EGFP (Ac element), NPTII, CaMV 35S promoter (P35S) driving transposase and mannopine synthase promoter (Pmas) driving EGFP (Ds element). Of 180 primary transgenics, 48 were potential launch pads, 72 were multiple insertions or chimaeras, and 60 exhibited somatic transposition. T1 progeny of 32 putative launch pads were screened by multiplex PCR for transposition. Evidence of germ-line transposition occurred in 13 putative launch pads; however, the transposition frequency was too low in three for efficient recovery of transposants. The transposition frequency in the remaining launch pads ranged from 16% to 40%. After self-pollination of the T0 launch pads, putative transposants in the T1 generation were identified by multiplex PCR. Sequencing of hiTAIL-PCR products derived from nested primers within the Ds end sequences (either P35S at the left border or the inverted repeat at the right border) of T1 plants revealed transposition of the Ds element to distant sites in the strawberry genome. From more than 2400 T1 plants screened, 103 unique transposants have been identified, among which 17 were somatic transpositions observed in the T0 generation. Ds insertion sites were dispersed among various gene elements [exons (15%), introns (23%), promoters (30%), 3' UTRs (17%) as well as intergenically (15%)]. Three-primer (one on either side of the Ds insertion and one within the Ds T-DNA) PCR could be used to identify homozygous T2 transposon-tagged plants. The mutant collection has been catalogued in an on-line database.

The role of uspE in virulence and biofilm formation by Histophilus somni.

UspE is a global regulator in Escherichia coli. To study the function of Histophilus somni uspE, strain 2336::TnuspE was identified from a bank of mutants generated with EZ::Tn5™ Tnp Transposome™ that were biofilm deficient. The 2336::TnuspE mutant was highly attenuated in mice, the electrophoretic profile of its lipooligosaccharide (LOS) indicated the LOS was truncated, and the mutant was significantly more serum-sensitive compared to the wildtype strain. In addition to forming a deficient biofilm, exopolysaccharide (EPS) production was also compromised, but the electrophoretic profile of outer membrane proteins was not altered. RNA sequence analysis revealed that the transcription levels of some stress response chaperones, transport proteins, and a large number of ribosomal protein genes in 2336::TnuspE were significantly differentially regulated compared to strain 2336. Therefore, uspE may differentially function in direct and indirect expression of H. somni genes, but its attenuation may be linked to poor biofilm formation and rapid clearance of the bacteria resulting from a compromised LOS structure. Our results support that uspE is a global stress regulatory gene in H. somni.

Phylogeny of gammaproteobacteria.

The phylogeny of the large bacterial class Gammaproteobacteria has been difficult to resolve. Here we apply a telescoping multiprotein approach to the problem for 104 diverse gammaproteobacterial genomes, based on a set of 356 protein families for the whole class and even larger sets for each of four cohesive subregions of the tree. Although the deepest divergences were resistant to full resolution, some surprising patterns were strongly supported. A representative of the Acidithiobacillales routinely appeared among the outgroup members, suggesting that in conflict with rRNA-based phylogenies this order does not belong to Gammaproteobacteria; instead, it (and, independently, "Mariprofundus") diverged after the establishment of the Alphaproteobacteria yet before the betaproteobacteria/gammaproteobacteria split. None of the orders Alteromonadales, Pseudomonadales, or Oceanospirillales were monophyletic; we obtained strong support for clades that contain some but exclude other members of all three orders. Extreme amino acid bias in the highly A+T-rich genome of Candidatus Carsonella prevented its reliable placement within Gammaproteobacteria, and high bias caused artifacts that limited the resolution of the relationships of other insect endosymbionts, which appear to have had multiple origins, although the unbiased genome of the endosymbiont Sodalis acted as an attractor for them. Instability was observed for the root of the Enterobacteriales, with nearly equal subsets of the protein families favoring one or the other of two alternative root positions; the nematode symbiont Photorhabdus was identified as a disruptor whose omission helped stabilize the Enterobacteriales root.

Implementing reverse genetics in Rosaceae: analysis of T-DNA flanking sequences of insertional mutant lines in the diploid strawberry, Fragaria vesca.

Reverse genetics is used for functional genomics research in model plants. To establish a model system for the systematic reverse genetics research in the Rosaceae family, we analyzed genomic DNA flanking the T-DNA insertions in 191 transgenic plants of the diploid strawberry, Fragaria vesca. One hundred and seventy-six T-DNA flanking sequences were amplified from the right border (RB) and 37 from the left border (LB) by thermal asymmetric interlaced PCR. Analysis of the T-DNA nick positions revealed that T-DNA was most frequently nicked at the cleavage sites. Analysis of 11 T-DNA integration sites indicated that T-DNA was integrated into the F. vesca genome by illegitimate recombination, as reported in other model plants: Arabidopsis, rice and tobacco. First, deletion of DNA was found at T-DNA integration target sites in all transgenic plants tested. Second, microsimilarities of a few base pairs between the left and/or right ends of the T-DNA and genomic sites were found in all transgenic plants tested. Finally, filler DNA was identified in four break-points. Out of 191 transgenic plants, T-DNA flanking sequences of 79 plants (41%) showed significant similarity to genes, elements or proteins of other plant species and 67 (35%) of the sequences are still unknown strawberry gene fragments. T-DNA flanking sequences of 126 plants (66%) showed homology to plant ESTs. This is the first report of T-DNA integration in a sizeable population of a rosaceous species. We have shown in this paper that T-DNA integration in strawberry is not random but directed by sequence microsimilarities in the host genome.

Identifying essential genes in Arabidopsis thaliana.

Eight years after publication of the Arabidopsis genome sequence and two years before completing the first phase of an international effort to characterize the function of every Arabidopsis gene, plant biologists remain unable to provide a definitive answer to the following basic question: what is the minimal gene set required for normal growth and development? The purpose of this review is to summarize different strategies employed to identify essential genes in Arabidopsis, an important component of the minimal gene set in plants, to present an overview of the datasets and specific genes identified to date, and to discuss the prospects for future saturation of this important class of genes. The long-term goal of this collaborative effort is to facilitate basic research in plant biology and complement ongoing research with other model organisms.

Polymicrobial Biofilm Interaction Between Histophilus somni and Pasteurella multocida.

Histophilus somni and Pasteurella multocida are two of multiple agents responsible for bovine respiratory disease (BRD) in cattle. Following respiratory infection of calves with H. somni, P. multocida may also be isolated from the lower respiratory tract. Because H. somni may form a biofilm during BRD, we sought to determine if P. multocida can co-exist with H. somni in a polymicrobial biofilm in vitro and in vivo. Interactions between the two species in the biofilm were characterized and quantified by fluorescence in situ hybridization (FISH). The biofilm matrix of each species was examined using fluorescently tagged lectins (FTL) specific for the exopolysaccharide (EPS) using confocal laser scanning microscopy. Bacterial interactions were determined by auto-aggregation and biofilm morphology. Pasteurella multocida and H. somni were evenly distributed in the in vitro biofilm, and both species contributed to the polymicrobial biofilm matrix. The average biomass and biofilm thickness, and the total carbohydrate and protein content of the biofilm, were greatest when both species were present. Polymicrobial bacterial suspensions auto-aggregated faster than single species suspensions, suggesting physical interactions between the two species. Almost 300 P. multocida genes were significantly differentially regulated when the bacteria were in a polymicrobial biofilm compared to a mono-species biofilm, as determined by RNA-sequencing. As expected, host genes associated with inflammation and immune response were significantly upregulated at the infection site following H. somni challenge. Encapsulated P. multocida isolates not capable of forming a substantial biofilm enhanced an in vitro polymicrobial biofilm with H. somni, indicating they contributed to the polymicrobial biofilm matrix. Indirect evidence indicated that encapsulated P. multocida also contributed to a polymicrobial biofilm in vivo. Only the EPS of H. somni could be detected by FTL staining of bovine tissues following challenge with H. somni. However, both species were isolated and an immune response to the biofilm matrix of both species was greater than the response to planktonic cells, suggesting encapsulated P. multocida may take advantage of the H. somni biofilm to persist in the host during chronic BRD. These results may have important implications for the management and prevention of BRD.

Evaluating the impact of international airline suspensions on the early global spread of COVID-19.

Global airline networks play a key role in the global importation of emerging infectious diseases. Detailed information on air traffic between international airports has been demonstrated to be useful in retrospectively validating and prospectively predicting case emergence in other countries. In this paper, we use a well-established metric known as effective distance on the global air traffic data from IATA to quantify risk of emergence for different countries as a consequence of direct importation from China, and compare it against arrival times for the first 24 countries. Using this model trained on official first reports from WHO, we estimate time of arrival (ToA) for all other countries. We then incorporate data on airline suspensions to recompute the effective distance and assess the effect of such cancellations in delaying the estimated arrival time for all other countries. Finally we use the infectious disease vulnerability indices to explain some of the estimated reporting delays.

Genome sequences of three agrobacterium biovars help elucidate the evolution of multichromosome genomes in bacteria.

The family Rhizobiaceae contains plant-associated bacteria with critical roles in ecology and agriculture. Within this family, many Rhizobium and Sinorhizobium strains are nitrogen-fixing plant mutualists, while many strains designated as Agrobacterium are plant pathogens. These contrasting lifestyles are primarily dependent on the transmissible plasmids each strain harbors. Members of the Rhizobiaceae also have diverse genome architectures that include single chromosomes, multiple chromosomes, and plasmids of various sizes. Agrobacterium strains have been divided into three biovars, based on physiological and biochemical properties. The genome of a biovar I strain, A. tumefaciens C58, has been previously sequenced. In this study, the genomes of the biovar II strain A. radiobacter K84, a commercially available biological control strain that inhibits certain pathogenic agrobacteria, and the biovar III strain A. vitis S4, a narrow-host-range strain that infects grapes and invokes a hypersensitive response on nonhost plants, were fully sequenced and annotated. Comparison with other sequenced members of the Alphaproteobacteria provides new data on the evolution of multipartite bacterial genomes. Primary chromosomes show extensive conservation of both gene content and order. In contrast, secondary chromosomes share smaller percentages of genes, and conserved gene order is restricted to short blocks. We propose that secondary chromosomes originated from an ancestral plasmid to which genes have been transferred from a progenitor primary chromosome. Similar patterns are observed in select Beta- and Gammaproteobacteria species. Together, these results define the evolution of chromosome architecture and gene content among the Rhizobiaceae and support a generalized mechanism for second-chromosome formation among bacteria.

Complete genome sequence of citrus huanglongbing bacterium, 'Candidatus Liberibacter asiaticus' obtained through metagenomics.

Citrus huanglongbing is the most destructive disease of citrus worldwide. It is spread by citrus psyllids and is associated with a low-titer, phloem-limited infection by any of three uncultured species of alpha-Proteobacteria, 'Candidatus Liberibacter asiaticus', 'Ca. L. americanus', and 'Ca. L. africanus'. A complete circular 'Ca. L. asiaticus' genome has been obtained by metagenomics, using the DNA extracted from a single 'Ca. L. asiaticus'-infected psyllid. The 1.23-Mb genome has an average 36.5% GC content. Annotation revealed a high percentage of genes involved in both cell motility (4.5%) and active transport in general (8.0%), which may contribute to its virulence. 'Ca. L. asiaticus' appears to have a limited ability for aerobic respiration and is likely auxotrophic for at least five amino acids. Consistent with its intracellular nature, 'Ca. L. asiaticus' lacks type III and type IV secretion systems as well as typical free-living or plant-colonizing extracellular degradative enzymes. 'Ca. L. asiaticus' appears to have all type I secretion system genes needed for both multidrug efflux and toxin effector secretion. Multi-protein phylogenetic analysis confirmed 'Ca. L. asiaticus' as an early-branching and highly divergent member of the family Rhizobiaceae. This is the first genome sequence of an uncultured alpha-proteobacteria that is both an intracellular plant pathogen and insect symbiont.

GeneTrees: a phylogenomics resource for prokaryotes.

The GeneTrees phylogenomics system pursues comparative genomic analyses from the perspective of gene phylogenies for individual genes. The GeneTrees project has the goal of providing detailed evolutionary models for all protein-coding gene components of the fully sequenced genomes. Currently, a database of alignments and trees for all protein sequences for 325 fully sequenced and annotated prokaryote genomes is available. The prokaryote database contains 890,000 protein sequences organized into over 100,000 alignments, each described by a phylogenetic tree. An original homology group discovery tool assembles sets of related proteins from all versus all pairwise alignments. Multiple alignments for each homology group are stored and subjected to phylogenetic tree inference. A graphical web interface provides visual exploration of the GeneTrees database. Homology groups can be queried by sequence identifiers or annotation terms. Genomes can be browsed visually on a gene map of each chromosome or plasmid. Phylogenetic trees with support values are displayed in conjunction with the associated sequence alignment. A variety of classes of information can be selected to label the tree tips to aid in visual evaluation of annotation and gene function. This web interface is available at http://genetrees.vbi.vt.edu.