Genomic and bioinformatics analysis of HAdV-7, a human adenovirus of species B1 that causes acute respiratory disease: implications for vector development in human gene therapy

AbstractHuman adenovirus serotype 7 (HAdV-7) is a reemerging pathogen identified in acute respiratory disease (ARD), particularly in epidemics affecting basic military trainee populations of otherwise healthy young adults. The genome has been sequenced and annotated (GenBank accession no. AY594255). Comparative genomics and bioinformatics analyses of the HAdV-7 genome sequence provide insight into its natural history and phylogenetic relationships. A putative origin of HAdV-7 from a chimpanzee host is observed. This has implications within the current biotechnological interest of using chimpanzee adenoviruses as vectors for human gene therapy and DNA vaccine delivery. Rapid genome sequencing and analyses of this species B1 member provide an example of exploiting accurate low-pass DNA sequencing technology in pathogen characterization and epidemic outbreak surveillance through the identification, validation, and application of unique pathogen genome signatures

GeneOrder3.0: Software for comparing the order of genes in pairs of small bacterial genomes

BACKGROUND: An increasing number of whole viral and bacterial genomes are being sequenced and deposited in public databases. In parallel to the mounting interest in whole genomes, the number of whole genome analyses software tools is also increasing. GeneOrder was originally developed to provide an analysis of genes between two genomes, allowing visualization of gene order and synteny comparisons of any small genomes. It was originally developed for comparing virus, mitochondrion and chloroplast genomes. This is now extended to small bacterial genomes of sizes less than 2 Mb. RESULTS: GeneOrder3.0 has been developed and validated successfully on several small bacterial genomes (ca. 580 kb to 1.83 Mb) archived in the NCBI GenBank database. It is an updated web-based "on-the-fly" computational tool allowing gene order and synteny comparisons of any two small bacterial genomes. Analyses of several bacterial genomes show that a large amount of gene and genome re-arrangement occurs, as seen with earlier DNA software tools. This can be displayed at the protein level using GeneOrder3.0. Whole genome alignments of genes are presented in both a table and a dot plot. This allows the detection of evolutionary more distant relationships since protein sequences are more conserved than DNA sequences. CONCLUSIONS: GeneOrder3.0 allows researchers to perform comparative analysis of gene order and synteny in genomes of sizes up to 2 Mb "on-the-fly." Availability: and

A Versatile Computational Pipeline for Bacterial Genome Annotation Improvement and Comparative Analysis, with \u3cem\u3eBrucella\u3c/em\u3e as a Use Case

We present a bacterial genome computational analysis pipeline, called GenVar. The pipeline, based on the program GeneWise, is designed to analyze an annotated genome and automatically identify missed gene calls and sequence variants such as genes with disrupted reading frames (split genes) and those with insertions and deletions (indels). For a given genome to be analyzed, GenVar relies on a database containing closely related genomes (such as other species or strains) as well as a few additional reference genomes. GenVar also helps identify gene disruptions probably caused by sequencing errors. We exemplify GenVar’s capabilities by presenting results from the analysis of four Brucella genomes. Brucella is an important human pathogen and zoonotic agent. The analysis revealed hundreds of missed gene calls, new split genes and indels, several of which are species specific and hence provide valuable clues to the understanding of the genome basis of Brucella pathogenicity and host specificity

Co-infections of Adenovirus Species in Previously Vaccinated Patients

Adenoviral infections associated with respiratory illness in military trainees involve multiple co-infecting species and serotypes

Plasmids and Rickettsial Evolution: Insight from Rickettsia felis

BACKGROUND: The genome sequence of Rickettsia felis revealed a number of rickettsial genetic anomalies that likely contribute not only to a large genome size relative to other rickettsiae, but also to phenotypic oddities that have confounded the categorization of R. felis as either typhus group (TG) or spotted fever group (SFG) rickettsiae. Most intriguing was the first report from rickettsiae of a conjugative plasmid (pRF) that contains 68 putative open reading frames, several of which are predicted to encode proteins with high similarity to conjugative machinery in other plasmid-containing bacteria. METHODOLOGY/PRINCIPAL FINDINGS: Using phylogeny estimation, we determined the mode of inheritance of pRF genes relative to conserved rickettsial chromosomal genes. Phylogenies of chromosomal genes were in agreement with other published rickettsial trees. However, phylogenies including pRF genes yielded different topologies and suggest a close relationship between pRF and ancestral group (AG) rickettsiae, including the recently completed genome of R. bellii str. RML369-C. This relatedness is further supported by the distribution of pRF genes across other rickettsiae, as 10 pRF genes (or inactive derivatives) also occur in AG (but not SFG) rickettsiae, with five of these genes characteristic of typical plasmids. Detailed characterization of pRF genes resulted in two novel findings: the identification of oriV and replication termination regions, and the likelihood that a second proposed plasmid, pRFδ, is an artifact of the original genome assembly. CONCLUSION/SIGNIFICANCE: Altogether, we propose a new rickettsial classification scheme with the addition of a fourth lineage, transitional group (TRG) rickettsiae, that is unique from TG and SFG rickettsiae and harbors genes from possible exchanges with AG rickettsiae via conjugation. We offer insight into the evolution of a plastic plasmid system in rickettsiae, including the role plasmids may have played in the acquirement of virulence traits in pathogenic strains, and the likely origin of plasmids within the rickettsial tree

Rickettsia Phylogenomics: Unwinding the Intricacies of Obligate Intracellular Life

BACKGROUND: Completed genome sequences are rapidly increasing for Rickettsia, obligate intracellular alpha-proteobacteria responsible for various human diseases, including epidemic typhus and Rocky Mountain spotted fever. In light of phylogeny, the establishment of orthologous groups (OGs) of open reading frames (ORFs) will distinguish the core rickettsial genes and other group specific genes (class 1 OGs or C1OGs) from those distributed indiscriminately throughout the rickettsial tree (class 2 OG or C2OGs). METHODOLOGY/PRINCIPAL FINDINGS: We present 1823 representative (no gene duplications) and 259 non-representative (at least one gene duplication) rickettsial OGs. While the highly reductive (approximately 1.2 MB) Rickettsia genomes range in predicted ORFs from 872 to 1512, a core of 752 OGs was identified, depicting the essential Rickettsia genes. Unsurprisingly, this core lacks many metabolic genes, reflecting the dependence on host resources for growth and survival. Additionally, we bolster our recent reclassification of Rickettsia by identifying OGs that define the AG (ancestral group), TG (typhus group), TRG (transitional group), and SFG (spotted fever group) rickettsiae. OGs for insect-associated species, tick-associated species and species that harbor plasmids were also predicted. Through superimposition of all OGs over robust phylogeny estimation, we discern between C1OGs and C2OGs, the latter depicting genes either decaying from the conserved C1OGs or acquired laterally. Finally, scrutiny of non-representative OGs revealed high levels of split genes versus gene duplications, with both phenomena confounding gene orthology assignment. Interestingly, non-representative OGs, as well as OGs comprised of several gene families typically involved in microbial pathogenicity and/or the acquisition of virulence factors, fall predominantly within C2OG distributions. CONCLUSION/SIGNIFICANCE: Collectively, we determined the relative conservation and distribution of 14354 predicted ORFs from 10 rickettsial genomes across robust phylogeny estimation. The data, available at PATRIC (PathoSystems Resource Integration Center), provide novel information for unwinding the intricacies associated with Rickettsia pathogenesis, expanding the range of potential diagnostic, vaccine and therapeutic targets

Identification of a Mycobacterium tuberculosis Putative Classical Nitroreductase Gene Whose Expression Is Coregulated with That of the acr Gene within Macrophages, in Standing versus Shaking Cultures, and under Low Oxygen Conditions

Tuberculosis remains a leading killer worldwide, and new approaches for its treatment and prevention are urgently needed. This effort will benefit greatly from a better understanding of gene regulation in Mycobacterium tuberculosis, particularly with respect to this pathogen's response to its host environment. We examined the behavior of two promoters from the divergently transcribed M. tuberculosis genes acr/hspX/Rv2031c (α-crystallin homolog) and Rv2032/acg (acr-coregulated gene) by using a promoter-GFP fusion assay in Mycobacterium bovis BCG. We found that Rv2032 is a novel macrophage-induced gene whose expression is coregulated with that of acr. Relative levels of intracellular induction for both promoters were significantly affected by shallow standing versus shaking bacterial culture conditions prior to macrophage infection, and both promoters were strongly induced under low oxygen conditions. Deletion analyses showed that DNA sequences within a 43-bp region were required for expression of these promoters under all conditions. Multiple sequence alignment and database searches performed with PROBE indicated that Rv2032 is one of eight M. tuberculosis genes of previously unknown function that belong to an unusual superfamily of classical nitroreductases, which may have a role for bacteria within the host environment. These findings show that mycobacterial culture conditions can greatly influence the results and interpretation of subsequent gene regulation experiments. We propose that these differences might be exploited for dissection of the regulatory factors that affect mycobacterial gene expression within the host

Tracking ebolavirus genomic drift with a resequencing microarray.

Filoviruses are emerging pathogens that cause acute fever with high fatality rate and present a global public health threat. During the 2013-2016 Ebola virus outbreak, genome sequencing allowed the study of virus evolution, mutations affecting pathogenicity and infectivity, and tracing the viral spread. In 2018, early sequence identification of the Ebolavirus as EBOV in the Democratic Republic of the Congo supported the use of an Ebola virus vaccine. However, field-deployable sequencing methods are needed to enable a rapid public health response. Resequencing microarrays (RMA) are a targeted method to obtain genomic sequence on clinical specimens rapidly, and sensitively, overcoming the need for extensive bioinformatic analysis. This study presents the design and initial evaluation of an ebolavirus resequencing microarray (Ebolavirus-RMA) system for sequencing the major genomic regions of four Ebolaviruses that cause disease in humans. The design of the Ebolavirus-RMA system is described and evaluated by sequencing repository samples of three Ebolaviruses and two EBOV variants. The ability of the system to identify genetic drift in a replicating virus was achieved by sequencing the ebolavirus glycoprotein gene in a recombinant virus cultured under pressure from a neutralizing antibody. Comparison of the Ebolavirus-RMA results to the Genbank database sequence file with the accession number given for the source RNA and Ebolavirus-RMA results compared to Next Generation Sequence results of the same RNA samples showed up to 99% agreement

A Family of acr-Coregulated Mycobacterium tuberculosis Genes Shares a Common DNA Motif and Requires Rv3133c (dosR or devR) for Expression

Previous work has shown that the divergently transcribed Mycobacterium tuberculosis genes acr (hspX, Rv2031c) and acg (Rv2032) are induced under conditions of shallow standing culture and low oxygen and intracellularly within macrophages. We used a combination of computational and experimental methods to identify promoters for eight additional genes that are regulated in a similar manner and that comprise an acr-coregulated promoter (ACP) family. Transcriptional regulation of these ACP family members was evaluated by using a plasmid-based promoter-green fluorescent protein fusion system and flow cytometry. All promoters showed increased expression in shallow standing versus shaking cultures, in low- versus high-oxygen conditions, and intracellularly within macrophages versus extracellularly in tissue culture medium. However, there were quantitative differences in expression among promoters and among conditions for each promoter. A conserved 18-bp palindromic sequence motif was identified in all ACPs by Gibbs sampling-based computational analyses. Two such motifs overlap regions in the acr and acg promoters that were previously shown to be required for their expression. In addition, we found that 5% carbon dioxide was required for growth of Mycobacterium bovis BCG under microaerophilic (1.3% O(2)) culture conditions and fully prevented the growth cessation typically associated with rapid removal of oxygen. These findings are likely to be relevant to the in vivo environment and will contribute to our understanding of the pathogenesis of tuberculosis infection