Publisher Correction: Genomic insights into the 2016-2017 cholera epidemic in Yemen.

In the HTML version of this Letter, the affiliations for authors Andrew S. Azman, Dhirendra Kumar and Thandavarayan Ramamurthy were inverted (the PDF and print versions of the Letter were correct); the affiliations have been corrected online.

Genomic insights into the 2016-2017 cholera epidemic in Yemen.

Yemen is currently experiencing, to our knowledge, the largest cholera epidemic in recent history. The first cases were declared in September 2016, and over 1.1 million cases and 2,300 deaths have since been reported1. Here we investigate the phylogenetic relationships, pathogenesis and determinants of antimicrobial resistance by sequencing the genomes of Vibrio cholerae isolates from the epidemic in Yemen and recent isolates from neighbouring regions. These 116 genomic sequences were placed within the phylogenetic context of a global collection of 1,087 isolates of the seventh pandemic V. cholerae serogroups O1 and O139 biotype El Tor2-4. We show that the isolates from Yemen that were collected during the two epidemiological waves of the epidemic1-the first between 28 September 2016 and 23 April 2017 (25,839 suspected cases) and the second beginning on 24 April 2017 (more than 1 million suspected cases)-are V. cholerae serotype Ogawa isolates from a single sublineage of the seventh pandemic V. cholerae O1 El Tor (7PET) lineage. Using genomic approaches, we link the epidemic in Yemen to global radiations of pandemic V. cholerae and show that this sublineage originated from South Asia and that it caused outbreaks in East Africa before appearing in Yemen. Furthermore, we show that the isolates from Yemen are susceptible to several antibiotics that are commonly used to treat cholera and to polymyxin B, resistance to which is used as a marker of the El Tor biotype.

Intragenic Enhancers Attenuate Host Gene Expression.

Eukaryotic gene transcription is regulated at many steps, including RNA polymerase II (Pol II) recruitment, transcription initiation, promoter-proximal Pol II pause release, and transcription termination; however, mechanisms regulating transcription during productive elongation remain poorly understood. Enhancers, which activate gene transcription, themselves undergo Pol II-mediated transcription, but our understanding of enhancer transcription and enhancer RNAs (eRNAs) remains incomplete. Here we show that transcription at intragenic enhancers interferes with and attenuates host gene transcription during productive elongation. While the extent of attenuation correlates positively with nascent eRNA expression, the act of intragenic enhancer transcription alone, but not eRNAs, explains the attenuation. Through CRISPR/Cas9-mediated deletions, we demonstrate a physiological role for intragenic enhancer-mediated transcription attenuation in cell fate determination. We propose that intragenic enhancers not only enhance transcription of one or more genes from a distance but also fine-tune transcription of their host gene through transcription interference, facilitating differential utilization of the same regulatory element for disparate functions.

Decoding the function of bivalent chromatin in development and cancer.

Bivalent chromatin is characterized by the simultaneous presence of H3K4me3 and H3K27me3, histone modifications generally associated with transcriptionally active and repressed chromatin, respectively. Prevalent in embryonic stem cells (ESCs), bivalency is postulated to poise/prime lineage-controlling developmental genes for rapid activation during embryogenesis while maintaining a transcriptionally repressed state in the absence of activation cues; however, this hypothesis remains to be directly tested. Most gene promoters DNA hypermethylated in adult human cancers are bivalently marked in ESCs, and it was speculated that bivalency predisposes them for aberrant de novo DNA methylation and irreversible silencing in cancer, but evidence supporting this model is largely lacking. Here, we show that bivalent chromatin does not poise genes for rapid activation but protects promoters from de novo DNA methylation. Genome-wide studies in differentiating ESCs reveal that activation of bivalent genes is no more rapid than that of other transcriptionally silent genes, challenging the premise that H3K4me3 is instructive for transcription. H3K4me3 at bivalent promoters-a product of the underlying DNA sequence-persists in nearly all cell types irrespective of gene expression and confers protection from de novo DNA methylation. Bivalent genes in ESCs that are frequent targets of aberrant hypermethylation in cancer are particularly strongly associated with loss of H3K4me3/bivalency in cancer. Altogether, our findings suggest that bivalency protects reversibly repressed genes from irreversible silencing and that loss of H3K4me3 may make them more susceptible to aberrant DNA methylation in diseases such as cancer. Bivalency may thus represent a distinct regulatory mechanism for maintaining epigenetic plasticity.

False discovery rate: the Achilles' heel of proteogenomics.

Proteogenomics refers to the integrated analysis of the genome and proteome that leverages mass-spectrometry (MS)-based proteomics data to improve genome annotations, understand gene expression control through proteoforms and find sequence variants to develop novel insights for disease classification and therapeutic strategies. However, proteogenomic studies often suffer from reduced sensitivity and specificity due to inflated database size. To control the error rates, proteogenomics depends on the target-decoy search strategy, the de-facto method for false discovery rate (FDR) estimation in proteomics. The proteogenomic databases constructed from three- or six-frame nucleotide database translation not only increase the search space and compute-time but also violate the equivalence of target and decoy databases. These searches result in poorer separation between target and decoy scores, leading to stringent FDR thresholds. Understanding these factors and applying modified strategies such as two-pass database search or peptide-class-specific FDR can result in a better interpretation of MS data without introducing additional statistical biases. Based on these considerations, a user can interpret the proteogenomics results appropriately and control false positives and negatives in a more informed manner. In this review, first, we briefly discuss the proteogenomic workflows and limitations in database construction, followed by various considerations that can influence potential novel discoveries in a proteogenomic study. We conclude with suggestions to counter these challenges for better proteogenomic data interpretation.

Polyketide Quinones Are Alternate Intermediate Electron Carriers during Mycobacterial Respiration in Oxygen-Deficient Niches.

Mycobacterium tuberculosis (Mtb) adaptation to hypoxia is considered crucial to its prolonged latent persistence in humans. Mtb lesions are known to contain physiologically heterogeneous microenvironments that bring about differential responses from bacteria. Here we exploit metabolic variability within biofilm cells to identify alternate respiratory polyketide quinones (PkQs) from both Mycobacterium smegmatis (Msmeg) and Mtb. PkQs are specifically expressed in biofilms and other oxygen-deficient niches to maintain cellular bioenergetics. Under such conditions, these metabolites function as mobile electron carriers in the respiratory electron transport chain. In the absence of PkQs, mycobacteria escape from the hypoxic core of biofilms and prefer oxygen-rich conditions. Unlike the ubiquitous isoprenoid pathway for the biosynthesis of respiratory quinones, PkQs are produced by type III polyketide synthases using fatty acyl-CoA precursors. The biosynthetic pathway is conserved in several other bacterial genomes, and our study reveals a redox-balancing chemicocellular process in microbial physiology.

A comparative study on effect of heat stress on physiological and cellular responses of crossbred cattle and riverine buffalo in subtropical region of India.

The study was carried out to compare the in vitro and in vivo heat shock responses of cattle and buffaloes. The expression of heat responsive genes (HSP70 and HSF family) were studied in vitro in peripheral blood mononuclear cells (PBMCs) of cattle and buffalo. In vivo observations on animals were carried out to investigate the physiological responses of cattle and buffalo at different THI over a period of 14 months. The study indicated that onset and severity of heat stress at different THI varied significantly between cattle and buffalo. Rectal temperature (RT) showed a significant (p < 0.05) increase at THI 67 in buffaloes and at THI 68 in cattle. Significant (p < 0.01) differences in RT between the species were observed at THI 71, 72, and 73. Respiration rate (RR) significantly (p < 0.05) increased at THI 70 in both the species and significant (p < 0.05) differences in RR were observed between the species at THI 65, 68, 69, and 74. THI had significant (p < 0.05) effect on blood glucose and blood electrolytes of the species with increased levels at higher THI. Serum AST and ALT levels showed less pronounced changes over increasing THI. Heat stress-associated expressions of HSP 70 genes followed temporal changes with incremental THI. The expression of HSPA8 was consistent at lower THI whereas upregulation of HSPA1A and HSPA1L was evident at higher THI. The study concludes that changes in physiological parameters such as RT and RR occur in a phasic pattern in both species and onset of heat stress was early in buffalo as compared to cattle.

Genomic plasticity associated with antimicrobial resistance in Vibrio cholerae.

The Bay of Bengal is known as the epicenter for seeding several devastating cholera outbreaks across the globe. Vibrio cholerae, the etiological agent of cholera, has extraordinary competency to acquire exogenous DNA by horizontal gene transfer (HGT) and adapt them into its genome for structuring metabolic processes, developing drug resistance, and colonizing the human intestine. Antimicrobial resistance (AMR) in V. cholerae has become a global concern. However, little is known about the identity of the resistance traits, source of AMR genes, acquisition process, and stability of the genetic elements linked with resistance genes in V. cholerae Here we present details of AMR profiles of 443 V. cholerae strains isolated from the stool samples of diarrheal patients from two regions of India. We sequenced the whole genome of multidrug-resistant (MDR) and extensively drug-resistant (XDR) V. cholerae to identify AMR genes and genomic elements that harbor the resistance traits. Our genomic findings were further confirmed by proteome analysis. We also engineered the genome of V. cholerae to monitor the importance of the autonomously replicating plasmid and core genome in the resistance profile. Our findings provided insights into the genomes of recent cholera isolates and identified several acquired traits including plasmids, transposons, integrative conjugative elements (ICEs), pathogenicity islands (PIs), prophages, and gene cassettes that confer fitness to the pathogen. The knowledge generated from this study would help in better understanding of V. cholerae evolution and management of cholera disease by providing clinical guidance on preferred treatment regimens.

A facile and efficient multicomponent ultrasound-assisted "on water" synthesis of benzodiazepine ring.

A facile and efficient multicomponent synthesis of benzodiazepine ring in water under ultrasound irradiation is reported first time. The current procedure escapes traditional chromatography and purification process and provided the product in excellent yields of 95% as compared to conventional methods. The approach was also validated on gram-scale synthesis.

GLIS3 Transcriptionally Activates WNT Genes to Promote Differentiation of Human Embryonic Stem Cells into Posterior Neural Progenitors.

Anterior-posterior (A-P) specification of the neural tube involves initial acquisition of anterior fate followed by the induction of posterior characteristics in the primitive anterior neuroectoderm. Several morphogens have been implicated in the regulation of A-P neural patterning; however, our understanding of the upstream regulators of these morphogens remains incomplete. Here, we show that the Krüppel-like zinc finger transcription factor GLI-Similar 3 (GLIS3) can direct differentiation of human embryonic stem cells (hESCs) into posterior neural progenitor cells in lieu of the default anterior pathway. Transcriptomic analyses reveal that this switch in cell fate is due to rapid activation of Wingless/Integrated (WNT) signaling pathway. Mechanistically, through genome-wide RNA-Seq, ChIP-Seq, and functional analyses, we show that GLIS3 binds to and directly regulates the transcription of several WNT genes, including the strong posteriorizing factor WNT3A, and that inhibition of WNT signaling is sufficient to abrogate GLIS3-induced posterior specification. Our findings suggest a potential role for GLIS3 in the regulation of A-P specification through direct transcriptional activation of WNT genes. Stem Cells 2018 Stem Cells 2019;37:202-215.

An alkali tolerant α-l-rhamnosidase from Fusarium moniliforme MTCC-2088 used in de-rhamnosylation of natural glycosides.

Analkali tolerant α-l-rhamnosidase has been purified to homogeneity from the culture filtrate of a new fungal strain, Fusarium moniliforme MTCC-2088, using concentration by ultrafiltration and cation exchange chromatography on CM cellulose column. The molecular mass of the purified enzyme has been found to be 36.0 kDa using SDS-PAGE analysis. The Km value using p-nitrophenyl-α-l-rhamnopyranoside as the variable substrate in 0.2 M sodium phosphate buffer pH10.5 at50 °C was 0.50 mM. The catalytic rate constant was15.6 s-1giving the values of kcat/Km is 3.12 × 104M-1 s-1. The pH and temperature optima of the enzyme were 10.5 and 50 °C, respectively. The purified enzyme had better stability at 10 °C in basic pH medium. The enzyme derhamnosylated natural glycosides like naringin to prunin, rutin to isoquercitrin and hesperidin to hesperetin glucoside. The purified α-l-rhamnosidase has potential for enhancement of wine aroma.

Reassessment of temperature-humidity index for measuring heat stress in crossbred dairy cattle of a sub-tropical region.

The temperature-humidity index (THI) has been extensively applied for assessing heat stress in moderate to hot conditions in dairy cattle. However, there exist wide variation between researchers in defining an appropriate range of THI values for denoting different levels of stress. The present study was aimed to reassess previously described heat stress indicators of dairy cattle of sub-tropical region of India. From comparative evaluation of meteorological data over previous four years (2014-2017) the period of year when high THI prevailed in the region was determined. Accordingly, the time period of sample collection and observation on animals was decided, so that a THI range of 68-86 could be covered. After analyzing physiological, biochemical parameters and expression profile of heat shock response (HSR) genes of animals in response to different THI, it was evident from the study that animal undergoes few or little changes at THI 72, but major physiological changes occurred after THI reached 74. At THI range 74-79, no drastic change in these parameters occurred suggesting animals undergo transient acclimatization in this range to maintain homeostasis. Once THI reached and crossed 80, this homeostasis was perturbed and animals experienced major physiological changes again. Overall, the study suggests that THI values indicating level of heat stress are dependent on the geographic location, as well as type of animal and therefore, existing THI should be recalibrated for different climatic region for accurate assessment of the heat stress.

Integrated Transcriptomic-Proteomic Analysis Using a Proteogenomic Workflow Refines Rat Genome Annotation.

Proteogenomic re-annotation and mRNA splicing information can lead to the discovery of various protein forms for eukaryotic model organisms like rat. However, detection of novel proteoforms using mass spectrometry proteomics data remains a formidable challenge. We developed EuGenoSuite, an open source multiple algorithmic proteomic search tool and utilized it in our in-house integrated transcriptomic-proteomic pipeline to facilitate automated proteogenomic analysis. Using four proteogenomic pipelines (integrated transcriptomic-proteomic, Peppy, Enosi, and ProteoAnnotator) on publicly available RNA-sequence and MS proteomics data, we discovered 363 novel peptides in rat brain microglia representing novel proteoforms for 249 gene loci in the rat genome. These novel peptides aided in the discovery of novel exons, translation of annotated untranslated regions, pseudogenes, and splice variants for various loci; many of which have known disease associations, including neurological disorders like schizophrenia, amyotrophic lateral sclerosis, etc. Novel isoforms were also discovered for genes implicated in cardiovascular diseases and breast cancer for which rats are considered model organisms. Our integrative multi-omics data analysis not only enables the discovery of new proteoforms but also generates an improved reference for human disease studies in the rat model.

Proteogenomics of rare taxonomic phyla: A prospective treasure trove of protein coding genes.

Sustainable innovations in sequencing technologies have resulted in a torrent of microbial genome sequencing projects. However, the prokaryotic genomes sequenced so far are unequally distributed along their phylogenetic tree; few phyla contain the majority, the rest only a few representatives. Accurate genome annotation lags far behind genome sequencing. While automated computational prediction, aided by comparative genomics, remains a popular choice for genome annotation, substantial fraction of these annotations are erroneous. Proteogenomics utilizes protein level experimental observations to annotate protein coding genes on a genome wide scale. Benefits of proteogenomics include discovery and correction of gene annotations regardless of their phylogenetic conservation. This not only allows detection of common, conserved proteins but also the discovery of protein products of rare genes that may be horizontally transferred or taxonomy specific. Chances of encountering such genes are more in rare phyla that comprise a small number of complete genome sequences. We collated all bacterial and archaeal proteogenomic studies carried out to date and reviewed them in the context of genome sequencing projects. Here, we present a comprehensive list of microbial proteogenomic studies, their taxonomic distribution, and also urge for targeted proteogenomics of underexplored taxa to build an extensive reference of protein coding genes.

Integrating transcriptome and proteome profiling: Strategies and applications.

Discovering the gene expression signature associated with a cellular state is one of the basic quests in majority of biological studies. For most of the clinical and cellular manifestations, these molecular differences may be exhibited across multiple layers of gene regulation like genomic variations, gene expression, protein translation and post-translational modifications. These system wide variations are dynamic in nature and their crosstalk is overwhelmingly complex, thus analyzing them separately may not be very informative. This necessitates the integrative analysis of such multiple layers of information to understand the interplay of the individual components of the biological system. Recent developments in high throughput RNA sequencing and mass spectrometric (MS) technologies to probe transcripts and proteins made these as preferred methods for understanding global gene regulation. Subsequently, improvements in "big-data" analysis techniques enable novel conclusions to be drawn from integrative transcriptomic-proteomic analysis. The unified analyses of both these data types have been rewarding for several biological objectives like improving genome annotation, predicting RNA-protein quantities, deciphering gene regulations, discovering disease markers and drug targets. There are different ways in which transcriptomics and proteomics data can be integrated; each aiming for different research objectives. Here, we review various studies, approaches and computational tools targeted for integrative analysis of these two high-throughput omics methods.

Proteogenomic Tools and Approaches to Explore Protein Coding Landscapes of Eukaryotic Genomes.

Proteogenomic strategies aim to refine genome-wide annotations of protein coding features by using actual protein level observations. Most of the currently applied proteogenomic approaches include integrative analysis of multiple types of high-throughput omics data, e.g., genomics, transcriptomics, proteomics, etc. Recent efforts towards creating a human proteome map were primarily targeted to experimentally detect at least one protein product for each gene in the genome and extensively utilized proteogenomic approaches. The 14 year long wait to get a draft human proteome map, after completion of similar efforts to sequence the genome, explains the huge complexity and technical hurdles of such efforts. Further, the integrative analysis of large-scale multi-omics datasets inherent to these studies becomes a major bottleneck to their success. However, recent developments of various analysis tools and pipelines dedicated to proteogenomics reduce both the time and complexity of such analysis. Here, we summarize notable approaches, studies, software developments and their potential applications towards eukaryotic genome annotation and clinical proteogenomics.

Probing the Missing Human Proteome: A Computational Perspective.

The missing human proteome comprises predicted protein-coding genes with no credible protein level evidence detected so far and constitutes ~18% of the human protein coding genes (neXtProt release 19/9/2014). The missing proteins may be of pharmacological interest as many of these are membrane receptors, thus requiring comprehensive characterization. In the present study, we explored various computational parameters, crucial during protein searches from tandem mass spectrometry (MS) data, for their impact on missing protein identification. Variables taken into consideration are differences in search database composition, shared peptides, semitryptic searches, post-translational modifications (PTMs), and transcriptome guided proteogenomic searches. We used a multialgorithmic approach for protein detection from publicly available mass spectra from recent studies covering diverse human tissues and cell types. Using the aforementioned approaches, we successfully detected 24 missing proteins (22-PE2, 1-PE4, and 1-PE5). Maximum of these identifications could be attributed to differences in reference proteome databases, exemplifying use of a single standard database for human protein detection from MS data. Our results suggest that search strategies with modified parameters can be rewarding alternatives for extensive profiling of missing proteins. We conclude that using complementary spectral data searches incorporating different parameters like PTMs, against a comprehensive and compact search database, might lead to discoveries of the proteins attributed so far as the missing human proteome.

Transcriptome and venom proteome of the box jellyfish Chironex fleckeri.

BACKGROUND: The box jellyfish, Chironex fleckeri, is the largest and most dangerous cubozoan jellyfish to humans. It produces potent and rapid-acting venom and its sting causes severe localized and systemic effects that are potentially life-threatening. In this study, a combined transcriptomic and proteomic approach was used to identify C. fleckeri proteins that elicit toxic effects in envenoming. RESULTS: More than 40,000,000 Illumina reads were used to de novo assemble ∼ 34,000 contiguous cDNA sequences and ∼ 20,000 proteins were predicted based on homology searches, protein motifs, gene ontology and biological pathway mapping. More than 170 potential toxin proteins were identified from the transcriptome on the basis of homology to known toxins in publicly available sequence databases. MS/MS analysis of C. fleckeri venom identified over 250 proteins, including a subset of the toxins predicted from analysis of the transcriptome. Potential toxins identified using MS/MS included metalloproteinases, an alpha-macroglobulin domain containing protein, two CRISP proteins and a turripeptide-like protease inhibitor. Nine novel examples of a taxonomically restricted family of potent cnidarian pore-forming toxins were also identified. Members of this toxin family are potently haemolytic and cause pain, inflammation, dermonecrosis, cardiovascular collapse and death in experimental animals, suggesting that these toxins are responsible for many of the symptoms of C. fleckeri envenomation. CONCLUSIONS: This study provides the first overview of a box jellyfish transcriptome which, coupled with venom proteomics data, enhances our current understanding of box jellyfish venom composition and the molecular structure and function of cnidarian toxins. The generated data represent a useful resource to guide future comparative studies, novel protein/peptide discovery and the development of more effective treatments for jellyfish stings in humans. (Length: 300).

Topical application of substance P promotes wound healing in streptozotocin-induced diabetic rats.

Substance P (SP) is known to stimulate angiogenesis, fibroblasts proliferation and expressions of cytokines and growth factors involved in wound healing. However, SP level reduces in dermis in diabetics and, hence, it was hypothesized that exogenously applied SP could be helpful in improving wound healing in diabetic rats. Excision skin wound was created on the back of diabetic rats and rats were divided into three groups i.e. (i) saline-, (ii) gel- and (iii) SP-treated. Normal saline, pluronic gel and SP (10(-6)M) in gel were topically applied once daily for 19days. SP treatment significantly increased the wound closure, levels of interleukin-10, and expressions of vascular endothelial growth factor, transforming growth factor-beta1, heme oxygenase-1 and endothelial nitric oxide synthase, whereas it significantly decreased the expression of tumor necrosis factor-alpha, interleukin-1beta and matrix metalloproteinases-9 in the granulation/healing tissue. The inflammatory cells were present for long time in normal saline-treated group. Histological evaluation revealed better extracellular matrix formation with marked fibroblast proliferation and collagen deposition in SP-treated group. Early epithelial layer formation, increased microvessel density and greater growth associated protein-43 positive nerve fibers were also evidenced in SP-treated group. In conclusion, SP treatment markedly accelerated cutaneous wound healing in diabetic rats.

Betulinic acid attenuates lung injury by modulation of inflammatory cytokine response in experimentally-induced polymicrobial sepsis in mice.

Sepsis commonly progresses to acute lung injury (ALI), an inflammatory lung disease with high morbidity and mortality. Septic ALI is characterized by excessive production of proinflammatory mediators. It remained refractory to present therapies and new therapies need to be developed to improve further clinical outcomes. Betulinic acid (BA), a pentacyclic lupane group triterpenoid has been shown to have anti-inflammatory activities in many studies. However, its therapeutic efficacy in polymicrobial septic ALI is yet unknown. Therefore, we investigated the effects of BA on septic ALI using cecal ligation and puncture (CLP) model in mice. Vehicle or BA (3, 10, and 30mg/kg) was administered intraperitoneally, 3 times (0, 24 and 48h) before CLP and CLP was done on 49(th)h of the study. Survival rate was observed till 120h post CLP. Lung tissues were collected for analysis by sacrificing mice 18h post CLP. BA at 10 and 30mg/kg dose significantly reduced sepsis-induced mortality and lung injury as implied by attenuated lung histopathological changes, decreased protein and neutrophils infiltration. BA also decreased lung NF-κB expression, cytokine, intercellular adhesion molecule-1, monocyte chemoattractant protein-1 and matrix metalloproteinase-9 levels. These evidences suggest that, the protective effects of BA on lungs are associated with defending action against inflammatory response and BA could be a potential modulatory agent of inflammation in sepsis-induced ALI.