Hugh D. Mitchell

Challenges in biomarker discovery: combining expert insights with statistical analysis of complex omics data

INTRODUCTION: The advent of high throughput technologies capable of comprehensive analysis of genes, transcripts, proteins and other significant biological molecules has provided an unprecedented opportunity for the identification of molecular markers of disease processes. However, it has simultaneously complicated the problem of extracting meaningful molecular signatures of biological processes from these complex datasets. The process of biomarker discovery and characterization provides opportunities for more sophisticated approaches to integrating purely statistical and expert knowledge-based approaches. AREAS COVERED: In this review we will present examples of current practices for biomarker discovery from complex omic datasets and the challenges that have been encountered in deriving valid and useful signatures of disease. We will then present a high-level review of data-driven (statistical) and knowledge-based methods applied to biomarker discovery, highlighting some current efforts to combine the two distinct approaches. EXPERT OPINION: Effective, reproducible and objective tools for combining data-driven and knowledge-based approaches to identify predictive signatures of disease are key to future success in the biomarker field. We will describe our recommendations for possible approaches to this problem including metrics for the evaluation of biomarkers.

Salmonella modulates metabolism during growth under conditions that induce expression of virulence genes

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a facultative pathogen that uses complex mechanisms to invade and proliferate within mammalian host cells. To investigate possible contributions of metabolic processes to virulence in S. Typhimurium grown under conditions known to induce expression of virulence genes, we used a metabolomics-driven systems biology approach coupled with genome-scale modeling. First, we identified distinct metabolite profiles associated with bacteria grown in either rich or virulence-inducing media and report the most comprehensive coverage of the S. Typhimurium metabolome to date. Second, we applied an omics-informed genome-scale modeling analysis of the functional consequences of adaptive alterations in S. Typhimurium metabolism during growth under our conditions. Modeling efforts highlighted a decreased cellular capability to both produce and utilize intracellular amino acids during stationary phase culture in virulence conditions, despite significant abundance increases for these molecules as observed by our metabolomics measurements. Furthermore, analyses of omics data in the context of the metabolic model indicated rewiring of the metabolic network to support pathways associated with virulence. For example, cellular concentrations of polyamines were perturbed, as well as the predicted capacity for secretion and uptake.

A multi-omic systems approach to elucidating Yersinia virulence mechanisms

The underlying mechanisms that lead to dramatic differences between closely related pathogens are not always readily apparent. For example, the genomes of Yersinia pestis (YP) the causative agent of plague with a high mortality rate and Yersinia pseudotuberculosis (YPT) an enteric pathogen with a modest mortality rate are highly similar with some species specific differences; however the molecular causes of their distinct clinical outcomes remain poorly understood. In this study, a temporal multi-omic analysis of YP and YPT at physiologically relevant temperatures was performed to gain insights into how an acute and highly lethal bacterial pathogen, YP, differs from its less virulent progenitor, YPT. This analysis revealed higher gene and protein expression levels of conserved major virulence factors in YP relative to YPT, including the Yop virulon and the pH6 antigen. This suggests that adaptation in the regulatory architecture, in addition to the presence of unique genetic material, may contribute to the increased pathogenecity of YP relative to YPT. Additionally, global transcriptome and proteome responses of YP and YPT revealed conserved post-transcriptional control of metabolism and the translational machinery including the modulation of glutamate levels in Yersiniae. Finally, the omics data was coupled with a computational network analysis, allowing an efficient prediction of novel Yersinia virulence factors based on gene and protein expression patterns.

The highly conserved MraZ protein is a transcriptional regulator in Escherichia coli

The mraZ and mraW genes are highly conserved in bacteria, both in sequence and in their position at the head of the division and cell wall (dcw) gene cluster. Located directly upstream of the mraZ gene, the Pmra promoter drives the transcription of mraZ and mraW, as well as many essential cell division and cell wall genes, but no regulator of Pmra has been found to date. Although MraZ has structural similarity to the AbrB transition state regulator and the MazE antitoxin and MraW is known to methylate the 16S rRNA, mraZ and mraW null mutants have no detectable phenotypes. Here we show that overproduction of Escherichia coli MraZ inhibited cell division and was lethal in rich medium at high induction levels and in minimal medium at low induction levels. Co-overproduction of MraW suppressed MraZ toxicity, and loss of MraW enhanced MraZ toxicity, suggesting that MraZ and MraW have antagonistic functions. MraZ-green fluorescent protein localized to the nucleoid, suggesting that it binds DNA. Consistent with this idea, purified MraZ directly bound a region of DNA containing three direct repeats between Pmra and the mraZ gene. Excess MraZ reduced the expression of an mraZ-lacZ reporter, suggesting that MraZ acts as a repressor of Pmra, whereas a DNA-binding mutant form of MraZ failed to repress expression. Transcriptome sequencing (RNA-seq) analysis suggested that MraZ also regulates the expression of genes outside the dcw cluster. In support of this, purified MraZ could directly bind to a putative operator site upstream of mioC, one of the repressed genes identified by RNA-seq.

A Network Integration Approach to Predict Conserved Regulators Related to Pathogenicity of Influenza and SARS-CoV Respiratory Viruses

Respiratory infections stemming from influenza viruses and the Severe Acute Respiratory Syndrome corona virus (SARS-CoV) represent a serious public health threat as emerging pandemics. Despite efforts to identify the critical interactions of these viruses with host machinery, the key regulatory events that lead to disease pathology remain poorly targeted with therapeutics. Here we implement an integrated network interrogation approach, in which proteome and transcriptome datasets from infection of both viruses in human lung epithelial cells are utilized to predict regulatory genes involved in the host response. We take advantage of a novel “crowd-based” approach to identify and combine ranking metrics that isolate genes/proteins likely related to the pathogenicity of SARS-CoV and influenza virus. Subsequently, a multivariate regression model is used to compare predicted lung epithelial regulatory influences with data derived from other respiratory virus infection models. We predicted a small set of regulatory factors with conserved behavior for consideration as important components of viral pathogenesis that might also serve as therapeutic targets for intervention. Our results demonstrate the utility of integrating diverse ‘omic datasets to predict and prioritize regulatory features conserved across multiple pathogen infection models.

Transcriptomic and proteomic dynamics in the metabolism of a diazotrophic cyanobacterium, Cyanothece sp. PCC 7822 during a diurnal light–dark cycle

BackgroundCyanothece sp. PCC 7822 is an excellent cyanobacterial model organism with great potential to be applied as a biocatalyst for the production of high value compounds. Like other unicellular diazotrophic cyanobacterial species, it has a tightly regulated metabolism synchronized to the light–dark cycle. Utilizing transcriptomic and proteomic methods, we quantified the relationships between transcription and translation underlying central and secondary metabolism in response to nitrogen free, 12 hour light and 12 hour dark conditions.ResultsBy combining mass-spectrometry based proteomics and RNA-sequencing transcriptomics, we quantitatively measured a total of 6766 mRNAs and 1322 proteins at four time points across a 24 hour light–dark cycle. Photosynthesis, nitrogen fixation, and carbon storage relevant genes were expressed during the preceding light or dark period, concurrent with measured nitrogenase activity in the late light period. We describe many instances of disparity in peak mRNA and protein abundances, and strong correlation of light dependent expression of both antisense and CRISPR-related gene expression. The proteins for nitrogenase and the pentose phosphate pathway were highest in the dark, whereas those for glycolysis and the TCA cycle were more prominent in the light. Interestingly, one copy of the psbA gene encoding the photosystem II (PSII) reaction center protein D1 (psbA4) was highly upregulated only in the dark. This protein likely cannot catalyze O2 evolution and so may be used by the cell to keep PSII intact during N2 fixation. The CRISPR elements were found exclusively at the ends of the large plasmid and we speculate that their presence is crucial to the maintenance of this plasmid.ConclusionsThis investigation of parallel transcriptional and translational activity within Cyanothece sp. PCC 7822 provided quantitative information on expression levels of metabolic pathways relevant to engineering efforts. The identification of expression patterns for both mRNA and protein affords a basis for improving biofuel production in this strain and for further genetic manipulations. Expression analysis of the genes encoded on the 6 plasmids provided insight into the possible acquisition and maintenance of some of these extra-chromosomal elements.

A comprehensive collection of systems biology data characterizing the host response to viral infection

The Systems Biology for Infectious Diseases Research program was established by the U.S. National Institute of Allergy and Infectious Diseases to investigate host-pathogen interactions at a systems level. This program generated 47 transcriptomic and proteomic datasets from 30 studies that investigate in vivo and in vitro host responses to viral infections. Human pathogens in the Orthomyxoviridae and Coronaviridae families, especially pandemic H1N1 and avian H5N1 influenza A viruses and severe acute respiratory syndrome coronavirus (SARS-CoV), were investigated. Study validation was demonstrated via experimental quality control measures and meta-analysis of independent experiments performed under similar conditions. Primary assay results are archived at the GEO and PeptideAtlas public repositories, while processed statistical results together with standardized metadata are publically available at the Influenza Research Database (www.fludb.org) and the Virus Pathogen Resource (www.viprbrc.org). By comparing data from mutant versus wild-type virus and host strains, RNA versus protein differential expression, and infection with genetically similar strains, these data can be used to further investigate genetic and physiological determinants of host responses to viral infection.

MERS-CoV Accessory ORFs Play Key Role for Infection and Pathogenesis.

ABSTRACT While dispensable for viral replication, coronavirus (CoV) accessory open reading frame (ORF) proteins often play critical roles during infection and pathogenesis. Utilizing a previously generated mutant, we demonstrate that the absence of all four Middle East respiratory syndrome CoV (MERS-CoV) accessory ORFs (deletion of ORF3, -4a, -4b, and -5 [dORF3-5]) has major implications for viral replication and pathogenesis. Importantly, attenuation of the dORF3-5 mutant is primarily driven by dysregulated host responses, including disrupted cell processes, augmented interferon (IFN) pathway activation, and robust inflammation. In vitro replication attenuation also extends to in vivo models, allowing use of dORF3-5 as a live attenuated vaccine platform. Finally, examination of ORF5 implicates a partial role in modulation of NF-κB-mediated inflammation. Together, the results demonstrate the importance of MERS-CoV accessory ORFs for pathogenesis and highlight them as potential targets for surveillance and therapeutic treatments moving forward. IMPORTANCE The initial emergence and periodic outbreaks of MERS-CoV highlight a continuing threat posed by zoonotic pathogens to global public health. In these studies, mutant virus generation demonstrates the necessity of accessory ORFs in regard to MERS-CoV infection and pathogenesis. With this in mind, accessory ORF functions can be targeted for both therapeutic and vaccine treatments in response to MERS-CoV and related group 2C coronaviruses. In addition, disruption of accessory ORFs in parallel may offer a rapid response platform to attenuation of future emergent strains based on both SARS- and MERS-CoV accessory ORF mutants. IMPORTANCE The initial emergence and periodic outbreaks of MERS-CoV highlight a continuing threat posed by zoonotic pathogens to global public health. In these studies, mutant virus generation demonstrates the necessity of accessory ORFs in regard to MERS-CoV infection and pathogenesis. With this in mind, accessory ORF functions can be targeted for both therapeutic and vaccine treatments in response to MERS-CoV and related group 2C coronaviruses. In addition, disruption of accessory ORFs in parallel may offer a rapid response platform to attenuation of future emergent strains based on both SARS- and MERS-CoV accessory ORF mutants.

CpG preconditioning regulates miRNA expression that modulates genomic reprogramming associated with neuroprotection against ischemic injury

Cytosine-phosphate-guanine (CpG) preconditioning reprograms the genomic response to stroke to protect the brain against ischemic injury. The mechanisms underlying genomic reprogramming are incompletely understood. MicroRNAs (miRNAs) regulate gene expression; however, their role in modulating gene responses produced by CpG preconditioning is unknown. We evaluated brain miRNA expression in response to CpG preconditioning before and after stroke using microarray. Importantly, we have data from previous gene microarrays under the same conditions, which allowed integration of miRNA and gene expression data to specifically identify regulated miRNA gene targets. CpG preconditioning did not significantly alter miRNA expression before stroke, indicating that miRNA regulation is not critical for the initiation of preconditioning-induced neuroprotection. However, after stroke, differentially regulated miRNAs between CpG- and saline-treated animals associated with the upregulation of several neuroprotective genes, implicating these miRNAs in genomic reprogramming that increases neuroprotection. Statistical analysis revealed that the miRNA targets were enriched in the gene population regulated in the setting of stroke, implying that miRNAs likely orchestrate this gene expression. These data suggest that miRNAs regulate endogenous responses to stroke and that manipulation of these miRNAs may have the potential to acutely activate novel neuroprotective processes that reduce damage.

Middle East Respiratory Syndrome Coronavirus Nonstructural Protein 16 Is Necessary for Interferon Resistance and Viral Pathogenesis

Coronavirus (CoV) emergence in both humans and livestock represents a significant threat to global public health, as evidenced by the sudden emergence of severe acute respiratory syndrome CoV (SARS-CoV), MERS-CoV, porcine epidemic diarrhea virus, and swine delta CoV in the 21st century. These studies describe an approach that effectively targets the highly conserved 2′O-MTase activity of CoVs for attenuation. With clear understanding of the IFN/IFIT (IFN-induced proteins with tetratricopeptide repeats)-based mechanism, NSP16 mutants provide a suitable target for a live attenuated vaccine platform, as well as therapeutic development for both current and future emergent CoV strains. Importantly, other approaches targeting other conserved pan-CoV functions have not yet proven effective against MERS-CoV, illustrating the broad applicability of targeting viral 2′O-MTase function across CoVs. ABSTRACT Coronaviruses (CoVs) encode a mixture of highly conserved and novel genes, as well as genetic elements necessary for infection and pathogenesis, raising the possibility of common targets for attenuation and therapeutic design. In this study, we focused on highly conserved nonstructural protein 16 (NSP16), a viral 2′O-methyltransferase (2′O-MTase) that encodes critical functions in immune modulation and infection. Using reverse genetics, we disrupted a key motif in the conserved KDKE motif of Middle East respiratory syndrome CoV (MERS-CoV) NSP16 (D130A) and evaluated the effect on viral infection and pathogenesis. While the absence of 2′O-MTase activity had only a marginal impact on propagation and replication in Vero cells, dNSP16 mutant MERS-CoV demonstrated significant attenuation relative to the control both in primary human airway cell cultures and in vivo. Further examination indicated that dNSP16 mutant MERS-CoV had a type I interferon (IFN)-based attenuation and was partially restored in the absence of molecules of IFN-induced proteins with tetratricopeptide repeats. Importantly, the robust attenuation permitted the use of dNSP16 mutant MERS-CoV as a live attenuated vaccine platform protecting from a challenge with a mouse-adapted MERS-CoV strain. These studies demonstrate the importance of the conserved 2′O-MTase activity for CoV pathogenesis and highlight NSP16 as a conserved universal target for rapid live attenuated vaccine design in an expanding CoV outbreak setting. IMPORTANCE Coronavirus (CoV) emergence in both humans and livestock represents a significant threat to global public health, as evidenced by the sudden emergence of severe acute respiratory syndrome CoV (SARS-CoV), MERS-CoV, porcine epidemic diarrhea virus, and swine delta CoV in the 21st century. These studies describe an approach that effectively targets the highly conserved 2′O-MTase activity of CoVs for attenuation. With clear understanding of the IFN/IFIT (IFN-induced proteins with tetratricopeptide repeats)-based mechanism, NSP16 mutants provide a suitable target for a live attenuated vaccine platform, as well as therapeutic development for both current and future emergent CoV strains. Importantly, other approaches targeting other conserved pan-CoV functions have not yet proven effective against MERS-CoV, illustrating the broad applicability of targeting viral 2′O-MTase function across CoVs.