Enabling precision medicine via standard communication of HTS provenance, analysis, and results.

A personalized approach based on a patient's or pathogen's unique genomic sequence is the foundation of precision medicine. Genomic findings must be robust and reproducible, and experimental data capture should adhere to findable, accessible, interoperable, and reusable (FAIR) guiding principles. Moreover, effective precision medicine requires standardized reporting that extends beyond wet-lab procedures to computational methods. The BioCompute framework (https://w3id.org/biocompute/1.3.0) enables standardized reporting of genomic sequence data provenance, including provenance domain, usability domain, execution domain, verification kit, and error domain. This framework facilitates communication and promotes interoperability. Bioinformatics computation instances that employ the BioCompute framework are easily relayed, repeated if needed, and compared by scientists, regulators, test developers, and clinicians. Easing the burden of performing the aforementioned tasks greatly extends the range of practical application. Large clinical trials, precision medicine, and regulatory submissions require a set of agreed upon standards that ensures efficient communication and documentation of genomic analyses. The BioCompute paradigm and the resulting BioCompute Objects (BCOs) offer that standard and are freely accessible as a GitHub organization (https://github.com/biocompute-objects) following the "Open-Stand.org principles for collaborative open standards development." With high-throughput sequencing (HTS) studies communicated using a BCO, regulatory agencies (e.g., Food and Drug Administration [FDA]), diagnostic test developers, researchers, and clinicians can expand collaboration to drive innovation in precision medicine, potentially decreasing the time and cost associated with next-generation sequencing workflow exchange, reporting, and regulatory reviews.

Multiplex detection and identification of viral, bacterial, and protozoan pathogens in human blood and plasma using a high-density resequencing pathogen microarray platform.

BACKGROUND: The implementation of nucleic acid-based tests for blood donor screening has improved the safety of the blood supply; however, the increasing number of emerging pathogen tests is burdensome. Development of multiplex testing platforms that allow simultaneous screening for different pathogens is a potential solution. STUDY DESIGN AND METHODS: The TessArray resequencing microarray is a platform that allows multiplex detection and identification of 97 different blood-borne pathogens in one single test. The objective was to evaluate the lowest concentration detected in blood or plasma, species discrimination, and applicability of the TessArray microarray platform for testing blood donors. Human blood or plasma spiked with selected pathogens (10,000, 1000, or 100 cells or copies/mL), including three viral, four bacterial, and four protozoan pathogens were each tested on this platform. The nucleic acids were extracted, amplified using multiplexed sets of pooled specific primers, fragmented, labeled, and hybridized to a microarray. Finally, the detected sequences were identified using an automated genomic database alignment algorithm. RESULTS: The performance of this platform demonstrated detection for spiked bacterial and protozoan pathogens of 100 cells/mL and viral pathogens as low as 100 copies/mL. Coded specimens, including spiked and negative controls, were identified correctly for blood specimens (31/32, 97%) and for plasma specimens (20/22, 91%) demonstrating the effectiveness of the platform. CONCLUSION: These results indicated that the TessArray microarray platform could be employed for multiplex detection and identification, with a high level of discriminatory power for numerous blood-borne pathogen targets with potential for use in blood safety.

Multiplex detection and identification of viral, bacterial, and protozoan pathogens in human blood and plasma using an expanded high-density resequencing microarray platform.

Introduction: Nucleic acid tests for blood donor screening have improved the safety of the blood supply; however, increasing numbers of emerging pathogen tests are burdensome. Multiplex testing platforms are a potential solution. Methods: The Blood Borne Pathogen Resequencing Microarray Expanded (BBP-RMAv.2) can perform multiplex detection and identification of 80 viruses, bacteria and parasites. This study evaluated pathogen detection in human blood or plasma. Samples spiked with selected pathogens, each with one of 6 viruses, 2 bacteria and 5 protozoans were tested on this platform. The nucleic acids were extracted, amplified using multiplexed sets of primers, and hybridized to a microarray. The reported sequences were aligned to a database to identify the pathogen. To directly compare the microarray to an emerging molecular approach, the amplified nucleic acids were also submitted to nanopore next generation sequencing (NGS). Results: The BBP-RMAv.2 detected viral pathogens at a concentration as low as 100 copies/ml and a range of concentrations from 1,000 to 100,000 copies/ml for all the spiked pathogens. Coded specimens were identified correctly demonstrating the effectiveness of the platform. The nanopore sequencing correctly identified most samples and the results of the two platforms were compared. Discussion: These results indicated that the BBP-RMAv.2 could be employed for multiplex detection with potential for use in blood safety or disease diagnosis. The NGS was nearly as effective at identifying pathogens in blood and performed better than BBP-RMAv.2 at identifying pathogen-negative samples.

Correction: Tracking ebolavirus genomic drift with a resequencing microarray.

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

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.

Genomic and bioinformatics analyses of HAdV-14p, reference strain of a re-emerging respiratory pathogen and analysis of B1/B2.

Unlike other human adenovirus (HAdV) species, B is divided into subspecies B1 and B2. Originally this was partly based on restriction enzyme (RE) analysis. B1 members, except HAdV-50, are commonly associated with respiratory diseases while B2 members are rarely associated with reported respiratory diseases. Recently two members of B2 have been identified in outbreaks of acute respiratory disease (ARD). One, HAdV-14, has re-emerged after an apparent 52-year absence. Genomic analysis and bioinformatics data are reported for HAdV-14 prototype for use as a reference and to understand and counter its re-emergence. The data complement and extend the original criteria for subspecies designation, unique amongst the adenoviruses, and highlight differences between B1 and B2, representing the first comprehensive analysis of this division. These data also provide finer granularity into the pathoepidemiology of the HAdVs. Whole genome analysis uncovers heterogeneous identity structures of the hexon and fiber genes amongst the HAdV-14 and the B1/B2 subspecies, which may be important in prescient vaccine development. Analysis of cell surface proteins provides insight into HAdV-14 tropism, accounting for its role as a respiratory pathogen. This HAdV-14 prototype genome is also a reference for applications of B2 adenoviruses as vectors for vaccine development and gene therapy.

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: http://binf.gmu.edu/genometools.html and http://pasteur.atcc.org:8050/GeneOrder3.0.

Development and performance of a comprehensive targeted sequencing assay for pan-ethnic screening of carrier status.

Identifying individuals as carriers of severe disease traits enables informed decision making about reproductive options. Although carrier screening has traditionally been based on ethnicity, the increasing ethnic admixture in the general population argues for the need for pan-ethnic carrier screening assays. Highly multiplexed mutation panels allow for rapid and efficient testing of hundreds of mutations concurrently. We report the development of the Pan-Ethnic Carrier Screening assay, a targeted sequencing assay for routine screening that simultaneously detects 461 common mutations in 91 different genes underlying severe, early-onset monogenic disorders. Mutation selection was aided by the use of an extensive mutation database from a clinical laboratory with expertise in newborn screening and lysosomal storage disease testing. The assay is based on the Affymetrix GeneChip microarray platform but generates genomic DNA sequence as the output. Analytical sensitivity and specificity, using genomic DNA from archived control cultures and from clinical specimens, was found to be >99% for all mutation types. This targeted sequencing assay has advantages over multiplex PCR and next-generation sequencing assays, including accuracy of mutation detection over a range of mutation types and ease of analysis and reporting of results.

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.

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, pRFdelta, 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.

Co-infections of adenovirus species in previously vaccinated patients.

Despite the success of the adenovirus vaccine administered to US military trainees, acute respiratory disease (ARD) surveillance still detected breakthrough infections (respiratory illnesses associated with the adenovirus serotypes specifically targeted by the vaccine). To explore the role of adenoviral co-infection (simultaneous infection by multiple pathogenic adenovirus species) in breakthrough disease, we examined specimens from patients with ARD by using 3 methods to detect multiple adenoviral species: a DNA microarray, a polymerase chain reaction (PCR)-enzyme-linked immunosorbent assay, and a multiplex PCR assay. Analysis of 52 samples (21 vaccinated, 31 unvaccinated) collected from 1996 to 2000 showed that all vaccinated samples had co-infections. Most of these co-infections were community-acquired serotypes of species B1 and E. Unvaccinated samples primarily contained only 1 species (species E) associated with adult respiratory illness. This study highlights the rarely reported phenomenon of adenoviral co-infections in a clinically relevant environment suitable for the generation of new recombinational variants.

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.

Human 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. ). 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.

Genomic and bioinformatics analysis of HAdV-4, a human adenovirus causing acute respiratory disease: implications for gene therapy and vaccine vector development.

Human adenovirus serotype 4 (HAdV-4) is a reemerging viral pathogenic agent implicated in epidemic outbreaks of acute respiratory disease (ARD). This report presents a genomic and bioinformatics analysis of the prototype 35,990-nucleotide genome (GenBank accession no. AY594253). Intriguingly, the genome analysis suggests a closer phylogenetic relationship with the chimpanzee adenoviruses (simian adenoviruses) rather than with other human adenoviruses, suggesting a recent origin of HAdV-4, and therefore species E, through a zoonotic event from chimpanzees to humans. Bioinformatics analysis also suggests a pre-zoonotic recombination event, as well, between species B-like and species C-like simian adenoviruses. These observations may have implications for the current interest in using chimpanzee adenoviruses in the development of vectors for human gene therapy and for DNA-based vaccines. Also, the reemergence, surveillance, and treatment of HAdV-4 as an ARD pathogen is an opportunity to demonstrate the use of genome determination as a tool for viral infectious disease characterization and epidemic outbreak surveillance: for example, rapid and accurate low-pass sequencing and analysis of the genome. In particular, this approach allows the rapid identification and development of unique probes for the differentiation of family, species, serotype, and strain (e.g., pathogen genome signatures) for monitoring epidemic outbreaks of ARD.

Genomic and bioinformatics analyses of HAdV-4vac and HAdV-7vac, two human adenovirus (HAdV) strains that constituted original prophylaxis against HAdV-related acute respiratory disease, a reemerging epidemic disease.

Vaccine strains of human adenovirus serotypes 4 and 7 (HAdV-4vac and HAdV-7vac) have been used successfully to prevent adenovirus-related acute respiratory disease outbreaks. The genomes of these two vaccine strains have been sequenced, annotated, and compared with their prototype equivalents with the goals of understanding their genomes for molecular diagnostics applications, vaccine redevelopment, and HAdV pathoepidemiology. These reference genomes are archived in GenBank as HAdV-4vac (35,994 bp; AY594254) and HAdV-7vac (35,240 bp; AY594256). Bioinformatics and comparative whole-genome analyses with their recently reported and archived prototype genomes reveal six mismatches and four insertions-deletions (indels) between the HAdV-4 prototype and vaccine strains, in contrast to the 611 mismatches and 130 indels between the HAdV-7 prototype and vaccine strains. Annotation reveals that the HAdV-4vac and HAdV-7vac genomes contain 51 and 50 coding units, respectively. Neither vaccine strain appears to be attenuated for virulence based on bioinformatics analyses. There is evidence of genome recombination, as the inverted terminal repeat of HAdV-4vac is initially identical to that of species C whereas the prototype is identical to species B1. These vaccine reference sequences yield unique genome signatures for molecular diagnostics. As a molecular forensics application, these references identify the circulating and problematic 1950s era field strains as the original HAdV-4 prototype and the Greider prototype, from which the vaccines are derived. Thus, they are useful for genomic comparisons to current epidemic and reemerging field strains, as well as leading to an understanding of pathoepidemiology among the human adenoviruses.

Identification of a Mycobacterium tuberculosis putative classical nitroreductase gene whose expression is coregulated with that of the acr aene 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 (alpha-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.

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.

Natural variation among human adenoviruses: genome sequence and annotation of human adenovirus serotype 1.

The 36,001 base pair DNA sequence of human adenovirus serotype 1 (HAdV-1) has been determined, using a 'leveraged primer sequencing strategy' to generate high quality sequences economically. This annotated genome (GenBank AF534906) confirms anticipated similarity to closely related species C (formerly subgroup), human adenoviruses HAdV-2 and -5, and near identity with earlier reports of sequences representing parts of the HAdV-1 genome. A first round of HAdV-1 sequence data acquisition used PCR amplification and sequencing primers from sequences common to the genomes of HAdV-2 and -5. The subsequent rounds of sequencing used primers derived from the newly generated data. Corroborative re-sequencing with primers selected from this HAdV-1 dataset generated sparsely tiled arrays of high quality sequencing ladders spanning both complementary strands of the HAdV-1 genome. These strategies allow for rapid and accurate low-pass sequencing of genomes. Such rapid genome determinations facilitate the development of specific probes for differentiation of family, serotype, subtype and strain (e.g. pathogen genome signatures). These will be used to monitor epidemic outbreaks of acute respiratory disease in a defined test bed by the Epidemic Outbreak Surveillance (EOS) project.