Identification of Stable Reference Genes for Toxicogenomic and Gene Expression Analysis.

Gene expression analysis has been becoming a popular method for studying gene function and response to different environmental stresses, including toxin/pollution exposure. Selection of a suitable reference gene is critically important for gene expression analysis due to that wrong reference genes will cause misleading and even wrong conclusion. A good reference gene should be a more stable reference gene, particularly during the toxicant exposure treatment and/or other investigation condition. In this chapter, a step-by-step protocol is present for primer design, reverse transcription PCR, primer efficiency and specificity test, qRT-PCR, and the strategy for identifying most stable reference genes for toxicogenomic and gene expression analysis. The detailed method for determining the primer gene specificity and primer efficiency are also presented in this chapter. Low primer efficiency will affect the fold changes during gene expression analysis; however, it does not affect the conclusion, up- or downregulation. Choosing a wrong reference gene may result in wrong conclusion.

Evaluating drug targets through human loss-of-function genetic variation.

Naturally occurring human genetic variants that are predicted to inactivate protein-coding genes provide an in vivo model of human gene inactivation that complements knockout studies in cells and model organisms. Here we report three key findings regarding the assessment of candidate drug targets using human loss-of-function variants. First, even essential genes, in which loss-of-function variants are not tolerated, can be highly successful as targets of inhibitory drugs. Second, in most genes, loss-of-function variants are sufficiently rare that genotype-based ascertainment of homozygous or compound heterozygous 'knockout' humans will await sample sizes that are approximately 1,000 times those presently available, unless recruitment focuses on consanguineous individuals. Third, automated variant annotation and filtering are powerful, but manual curation remains crucial for removing artefacts, and is a prerequisite for recall-by-genotype efforts. Our results provide a roadmap for human knockout studies and should guide the interpretation of loss-of-function variants in drug development.

TraDIS-Xpress: a high-resolution whole-genome assay identifies novel mechanisms of triclosan action and resistance.

Understanding the genetic basis for a phenotype is a central goal in biological research. Much has been learnt about bacterial genomes by creating large mutant libraries and looking for conditionally important genes. However, current genome-wide methods are largely unable to assay essential genes which are not amenable to disruption. To overcome this limitation, we developed a new version of "TraDIS" (transposon directed insertion-site sequencing) that we term "TraDIS-Xpress" that combines an inducible promoter into the transposon cassette. This allows controlled overexpression and repression of all genes owing to saturation of inserts adjacent to all open reading frames as well as conventional inactivation. We applied TraDIS-Xpress to identify responses to the biocide triclosan across a range of concentrations. Triclosan is endemic in modern life, but there is uncertainty about its mode of action with a concentration-dependent switch from bacteriostatic to bactericidal action unexplained. Our results show a concentration-dependent response to triclosan with different genes important in survival between static and cidal exposures. These genes include those previously reported to have a role in triclosan resistance as well as a new set of genes, including essential genes. Novel genes identified as being sensitive to triclosan exposure include those involved in barrier function, small molecule uptake, and integrity of transcription and translation. We anticipate the approach we show here, by allowing comparisons across multiple experimental conditions of TraDIS data, and including essential genes, will be a starting point for future work examining how different drug conditions impact bacterial survival mechanisms.

Agreement between two large pan-cancer CRISPR-Cas9 gene dependency data sets.

Genome-scale CRISPR-Cas9 viability screens performed in cancer cell lines provide a systematic approach to identify cancer dependencies and new therapeutic targets. As multiple large-scale screens become available, a formal assessment of the reproducibility of these experiments becomes necessary. We analyze data from recently published pan-cancer CRISPR-Cas9 screens performed at the Broad and Sanger Institutes. Despite significant differences in experimental protocols and reagents, we find that the screen results are highly concordant across multiple metrics with both common and specific dependencies jointly identified across the two studies. Furthermore, robust biomarkers of gene dependency found in one data set are recovered in the other. Through further analysis and replication experiments at each institute, we show that batch effects are driven principally by two key experimental parameters: the reagent library and the assay length. These results indicate that the Broad and Sanger CRISPR-Cas9 viability screens yield robust and reproducible findings.

The bare necessities: Uncovering essential and condition-critical genes with transposon sequencing.

Querying gene function in bacteria has been greatly accelerated by the advent of transposon sequencing (Tn-seq) technologies (related Tn-seq strategies are known as TraDIS, INSeq, RB-TnSeq, and HITS). Pooled populations of transposon mutants are cultured in an environment and next-generation sequencing tools are used to determine areas of the genome that are important for bacterial fitness. In this review we provide an overview of Tn-seq methodologies and discuss how Tn-seq has been applied, or could be applied, to the study of oral microbiology. These applications include studying the essential genome as a means to rationally design therapeutic agents. Tn-seq has also contributed to our understanding of well-studied biological processes in oral bacteria. Other important applications include in vivo pathogenesis studies and use of Tn-seq to probe the molecular basis of microbial interactions. We also highlight recent advancements in techniques that act in synergy with Tn-seq such as clustered regularly interspaced short palindromic repeats (CRISPR) interference and microfluidic chip platforms.

6mer seed toxicity in tumor suppressive microRNAs.

Many small-interfering (si)RNAs are toxic to cancer cells through a 6mer seed sequence (positions 2-7 of the guide strand). Here we performed an siRNA screen with all 4096 6mer seeds revealing a preference for guanine in positions 1 and 2 and a high overall G or C content in the seed of the most toxic siRNAs for four tested human and mouse cell lines. Toxicity of these siRNAs stems from targeting survival genes with C-rich 3'UTRs. The master tumor suppressor miRNA miR-34a-5p is toxic through such a G-rich 6mer seed and is upregulated in cells subjected to genotoxic stress. An analysis of all mature miRNAs suggests that during evolution most miRNAs evolved to avoid guanine at the 5' end of the 6mer seed sequence of the guide strand. In contrast, for certain tumor-suppressive miRNAs the guide strand contains a G-rich toxic 6mer seed, presumably to eliminate cancer cells.

Stress effects of cyanotoxin ß-methylamino-L-alanine (BMAA) on cyanobacterial heterocyst formation and functionality.

Various species of cyanobacteria, diatoms and dinoflagellates are capable of synthesizing the non-proteinogenic neurotoxic amino acid ß-N-methylamino-L-alanine (BMAA), which is known to be a causative agent of human neurodegeneration. Similar to most cyanotoxins, the biological and ecological functions of BMAA in cyanobacteria are unknown. In this study, we show for the first time that BMAA, in micromolar amounts, inhibits the formation of heterocysts (specialized nitrogen-fixing cells) in heterocystous, diazotrophic cyanobacteria [Anabaena sp. PCC 7120, Nostoc punctiforme PCC 73102 (ATCC 29133), Nostoc sp. strain 8963] under conditions of nitrogen starvation. The inhibitory effect of BMAA is abolished by the addition of glutamate. To understand the genetic reason for the observed phenomenon, we used qPCR to study the expression of key genes involved in cell differentiation and nitrogen metabolism in the model cyanobacterium Anabaena sp. PCC 7120. We observed that in the presence of BMAA, Anabaena sp. PCC 7120 does not express two essential genes associated with heterocyst differentiation, namely, hetR and hepA. We also found that addition of BMAA to cyanobacterial cultures with mature heterocysts inhibits nifH gene expression and nitrogenase activity.

Probability of phenotypically detectable protein damage by ENU-induced mutations in the Mutagenetix database.

Computational inference of mutation effects is necessary for genetic studies in which many mutations must be considered as etiologic candidates. Programs such as PolyPhen-2 predict the relative severity of damage caused by missense mutations, but not the actual probability that a mutation will reduce/eliminate protein function. Based on genotype and phenotype data for 116,330 ENU-induced mutations in the Mutagenetix database, we calculate that putative null mutations, and PolyPhen-2-classified "probably damaging", "possibly damaging", or "probably benign" mutations have, respectively, 61%, 17%, 9.8%, and 4.5% probabilities of causing phenotypically detectable damage in the homozygous state. We use these probabilities in the estimation of genome saturation and the probability that individual proteins have been adequately tested for function in specific genetic screens. We estimate the proportion of essential autosomal genes in Mus musculus (C57BL/6J) and show that viable mutations in essential genes are more likely to induce phenotype than mutations in non-essential genes.

A Tet-Off gene expression system for validation of antifungal drug targets in a murine invasive pulmonary aspergillosis model.

Aspergillus fumigatus is one of the major causes of invasive pulmonary aspergillosis in immunocompromised patients. Novel antifungal therapy is in urgent need due to emerging resistance and adverse toxicity of current antifungal drugs. Gene products that are essential for Aspergillus viability during infection are attractive drug targets. To characterize these genes in vivo we developed a Tet-Off gene expression system in A. fumigatus, whereby the administration of doxycycline resulted in down regulation of the gene whose expression is under the control of the Tet-Off promoter. We tested the system on two potential drug targets, inosine 5'-monophosphate dehydrogenase (IMPDH) and L-ornithine N5-oxygenase (sidA) in a murine invasive pulmonary aspergillosis model. We show that depletion of IMPDH attenuated but did not completely abolish virulence in vivo whereas turning off the expression of sidA, which is required for iron acquisition, resulted in avirulence. We also investigated whether sidA expression could be controlled in a time-dependent manner in mice. Our results demonstrated that timing of doxycycline administration dramatically affects survival rate, suggesting that this genetic system can be used for testing whether an antifungal drug target is critical for fungal growth post-infection.

Using bacterial genomes and essential genes for the development of new antibiotics.

The shrinking antibiotic development pipeline together with the global increase in antibiotic resistant infections requires that new molecules with antimicrobial activity are developed. Traditional empirical screening approaches of natural and non-natural compounds have identified the majority of antibiotics that are currently available, however this approach has produced relatively few new antibiotics over the last few decades. The vast amount of bacterial genome sequence information that has become available since the sequencing of the first bacterial genome more than 20years ago holds potential for contributing to the discovery of novel antimicrobial compounds. Comparative genomic approaches can identify genes that are highly conserved within and between bacterial species, and thus may represent genes that participate in key bacterial processes. Whole genome mutagenesis studies can also identify genes necessary for bacterial growth and survival under different environmental conditions, making them attractive targets for the development of novel inhibitory compounds. In addition, transcriptomic and proteomic approaches can be used to characterize RNA and protein levels on a cellular scale, providing information on bacterial physiology that can be applied to antibiotic target identification. Finally, bacterial genomes can be mined to identify biosynthetic pathways that produce many intrinsic antimicrobial compounds and peptides. In this review, we provide an overview of past and current efforts aimed at using bacterial genomic data in the discovery and development of novel antibiotics.

Antibacterial drug discovery in the resistance era.

The looming antibiotic-resistance crisis has penetrated the consciousness of clinicians, researchers, policymakers, politicians and the public at large. The evolution and widespread distribution of antibiotic-resistance elements in bacterial pathogens has made diseases that were once easily treatable deadly again. Unfortunately, accompanying the rise in global resistance is a failure in antibacterial drug discovery. Lessons from the history of antibiotic discovery and fresh understanding of antibiotic action and the cell biology of microorganisms have the potential to deliver twenty-first century medicines that are able to control infection in the resistance era.

Selection and Evaluation of Reference Genes for Reverse Transcription-Quantitative PCR Expression Studies in a Thermophilic Bacterium Grown under Different Culture Conditions.

The phylum Deinococcus-Thermus is a deeply-branching lineage of bacteria widely recognized as one of the most extremophilic. Members of the Thermus genus are of major interest due to both their bioremediation and biotechnology potentials. However, the molecular mechanisms associated with these key metabolic pathways remain unknown. Reverse-transcription quantitative PCR (RT-qPCR) is a high-throughput means of studying the expression of a large suite of genes over time and under different conditions. The selection of a stably-expressed reference gene is critical when using relative quantification methods, as target gene expression is normalized to expression of the reference gene. However, little information exists as to reference gene selection in extremophiles. This study evaluated 11 candidate reference genes for use with the thermophile Thermus scotoductus when grown under different culture conditions. Based on the combined stability values from BestKeeper and NormFinder software packages, the following are the most appropriate reference genes when comparing: (1) aerobic and anaerobic growth: TSC_c19900, polA2, gyrA, gyrB; (2) anaerobic growth with varied electron acceptors: TSC_c19900, infA, pfk, gyrA, gyrB; (3) aerobic growth with different heating methods: gyrA, gap, gyrB; (4) all conditions mentioned above: gap, gyrA, gyrB. The commonly-employed rpoC does not serve as a reliable reference gene in thermophiles, due to its expression instability across all culture conditions tested here. As extremophiles exhibit a tendency for polyploidy, absolute quantification was employed to determine the ratio of transcript to gene copy number in a subset of the genes. A strong negative correlation was found to exist between ratio and threshold cycle (CT) values, demonstrating that CT changes reflect transcript copy number, and not gene copy number, fluctuations. Even with the potential for polyploidy in extremophiles, the results obtained via absolute quantification indicate that relative quantification is appropriate for RT-qPCR studies with this thermophile.

Validation of potential reference genes for qPCR in maize across abiotic stresses, hormone treatments, and tissue types.

The reverse transcription quantitative polymerase chain reaction (RT-qPCR) is a powerful and widely used technique for the measurement of gene expression. Reference genes, which serve as endogenous controls ensure that the results are accurate and reproducible, are vital for data normalization. To bolster the literature on reference gene selection in maize, ten candidate reference genes, including eight traditionally used internal control genes and two potential candidate genes from our microarray datasets, were evaluated for expression level in maize across abiotic stresses (cold, heat, salinity, and PEG), phytohormone treatments (abscisic acid, salicylic acid, jasmonic acid, ethylene, and gibberellins), and different tissue types. Three analytical software packages, geNorm, NormFinder, and Bestkeeper, were used to assess the stability of reference gene expression. The results revealed that elongation factor 1 alpha (EF1α), tubulin beta (ß-TUB), cyclophilin (CYP), and eukaryotic initiation factor 4A (EIF4A) were the most reliable reference genes for overall gene expression normalization in maize, while GRP (Glycine-rich RNA-binding protein), GLU1(beta-glucosidase), and UBQ9 (ubiquitin 9) were the least stable and most unsuitable genes. In addition, the suitability of EF1α, ß-TUB, and their combination as reference genes was confirmed by validating the expression of WRKY50 in various samples. The current study indicates the appropriate reference genes for the urgent requirement of gene expression normalization in maize across certain abiotic stresses, hormones, and tissue types.

Mood-stabilizers differentially affect housekeeping gene expression in human cells.

Recent studies have revealed that antidepressants affect the expression of constitutively expressed "housekeeping genes" commonly used as normalizing reference genes in quantitative polymerase chain reaction (qPCR) experiments. There has yet to be an investigation however on the effects of mood-stabilizers on housekeeping gene stability. The current study utilized lymphoblastoid cell lines (LCLs) derived from patients with mood disorders to investigate the effects of a range of doses of lithium (0, 1, 2 and 5 mM) and sodium valproate (0, 0.06, 0.03 and 0.6 mM) on the stability of 12 housekeeping genes. RNA was extracted from LCLs and qPCR was used to generate cycle threshold (Ct ) values which were input into RefFinder analyses. The study revealed drug-specific effects on housekeeping gene stability. The most stable housekeeping genes in LCLs treated: acutely with sodium valproate were ACTB and RPL13A; acutely with lithium were GAPDH and ATP5B; chronically with lithium were ATP5B and CYC1. The stability of GAPDH and B2M were particularly affected by duration of lithium treatment. The study adds to a growing literature that the selection of appropriate housekeeping genes is important for the accurate normalization of target gene expression in experiments investigating the molecular effects of mood disorder pharmacotherapies.

Assessing the essentiality of the decaprenyl-phospho-d-arabinofuranose pathway in Mycobacterium tuberculosis using conditional mutants.

In Mycobacterium tuberculosis the decaprenyl-phospho-d-arabinofuranose (DPA) pathway is a validated target for the drugs ethambutol and benzothiazinones. To identify other potential drug targets in the pathway, we generated conditional knock-down mutants of each gene involved using the TET-PIP OFF system. dprE1, dprE2, ubiA, prsA, rv2361c, tkt and rpiB were confirmed to be essential under non-permissive conditions, whereas rv3807c was not required for survival. In the most vulnerable group, DprE1-depleted cells died faster in vitro and intracellularly than those lacking UbiA and PrsA. Downregulation of DprE1 and UbiA resulted in similar phenotypes, namely swelling of the bacteria, cell wall damage and lysis as observed at the single cell level, by real time microscopy and electron microscopy. By contrast, depletion of PrsA led to cell elongation and implosion, which was suggestive of a more pleiotropic effect. Drug sensitivity assays with known DPA-inhibitors supported the use of conditional knock-down strains for target-based whole-cell screens. Together, our work provides strong evidence for the vulnerability of all but one of the enzymes in the DPA pathway and generates valuable tools for the identification of lead compounds targeting the different biosynthetic steps. PrsA, phosphoribosyl-pyrophosphate synthetase, appears to be a particularly attractive new target for drug discovery.

Validating reference genes within a mouse model system of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) toxicity.

BACKGROUND: Quantitative real-time PCR (qPCR) is the "gold-standard" technique for measuring mRNA abundances. To correctly compare samples and generate biologically valid results, qPCR data usually require comprehensive normalization to account for sample content variation between reactions. The most common normalization approaches use one or more endogenous controls (reference or house-keeping genes) to adjust the measured levels of experimental genes appropriately. Ideal reference genes are those that display minimal variation across experimental conditions, and thus can vary widely across different biological systems. In particular, toxicogenomic studies of transcriptionally-disruptive toxins, like 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), require careful consideration of reference genes. RESULTS: We examined seven candidate reference genes in 199 mice varying in genotype and time/dose of TCDD exposure. We assessed gene-stability in four ways: (1) the variance of the raw Cq values across biological replicates, (2) the fold-change from basal mRNA levels following treatment, (3) the inter- and intra-group stability evaluated using the NormFinder algorithm, (4) the comparative ΔCq method for each candidate gene. Univariate analyses showed Hprt and Eef1a1 are the two most stable individual reference genes. It has been suggested that using multiple genes would produce a more consistent normalization factor; multivariate analysis was performed using NormFinder. In general, stability increased with the number of genes used, but specific gene-combinations synergized. CONCLUSIONS: We have validated seven reference genes for use in analyzing mRNA abundances in mouse models of TCDD toxicity. The use of multiple reference genes increases stability, providing more consistent normalization and more reliable results. The number of reference genes used should be maximized, based on experimental capabilities (platform, sample availability, etc.). Our results show the benefit of validating reference genes using multiple methods prior to generating large biological datasets.

Construction of the control system of target molecule expression in Escherichia coli: application to a validation platform for bactericidal and bacteriostatic profiles due to suppression of a target molecule.

Validation of bactericidal profiles owing to a deficiency of target bacterial molecule provides opportunities to discover antimicrobial drug candidates. In this study, we constructed genetic-engineered Escherichia coli strains, in which the target gene expression is conditionally regulated by a tryptophan promoter, while the target protein expression is regulated by N-end rule-based proteolysis. Among 10 genes, whose correspondent proteins are target candidates of antibiotics for community acquired respiratory tract infection, it was clearly demonstrated that the suppression of DnaB, GlmU, or DnaX results in a bactericidal profile, while the suppression of FabB, PyrG, DnaG, Der, PyrH, Era, or IspA leads to a bacteriostatic profile. This study is the first to predict the antibacterial inhibition profiles of Der, DnaG, DnaX, Era, GlmU, IspA, PyrG, and PyrH, and confirms previous findings for DnaB and FabB. The results suggested that the system constructed in this study is a novel and useful tool to validate whether the target bacterial molecule has appropriate properties as a target of antimicrobial agents.

Reference genes for qPCR assays in toxic metal and salinity stress in two flatworm model organisms.

The flatworm species Schmidtea mediterranea and Macrostomum lignano have become new and innovative model organisms in stem cell, regeneration and tissue homeostasis research. Because of their unique stem cell system, (lab) technical advantages and their phylogenetic position within the Metazoa, they are also ideal candidate model organisms for toxicity assays. As stress and biomarker screenings are often performed at the transcriptional level, the aim of this study was to establish a set of reference genes for qPCR experiments for these two model organisms in different stress situations. We examined the transcriptional stability of nine potential reference genes (actb, tubb, ck2, cox4, cys, rpl13, gapdh, gm2ap, plscr1) to assess those that are most stable during altered stress conditions (exposure to carcinogenic metals and salinity stress). The gene expression stability was evaluated by means of geNorm and NormFinder algorithms. Sets of best reference genes in these analyses varied between different stress situations, although gm2ap and actb were stably transcribed during all tested combinations. In order to demonstrate the impact of bad normalisation, the stress-specific gene hsp90 was normalised to different sets of reference genes. In contrast to the normalisation according to GeNorm and NormFinder, normalisation of hsp90 in Macrostomum lignano during cadmium stress did not show a significant difference when normalised to only gapdh. On the other hand an increase of variability was noticed when normalised to all nine tested reference genes together. Testing appropriate reference genes is therefore strongly advisable in every new experimental condition.

Oral toxicity of Photorhabdus culture media on gene expression of the adult sweetpotato whitefly, Bemisia tabaci.

The oral toxicity of culture media of the symbiotic bacteria, Photorhabdus temperata, mutually associated with entomopathogenic nematode Heterorhabditis megidis and Photorhabdus luminescens ssp. laumondii (TT01) mutually associated with Heterorhabditis bacteriophora, were investigated in the adults of Bemisia tabaci. The oral ingestion of sucrose diet solutions (20%) containing bacteria-free supernatant of the culture media from symbiotic bacteria gradually increased mortalities and was completely lethal at 60 h after the treatments, whereas the mortalities of the controls, sucrose solutions with or without media that uncultured with bacteria, were less than 17% up to 84 h of incubation. The effects of oral ingestion of symbiont culture media were demonstrated on the expression rates of several genes of B. tabaci using quantitative real-time RT-PCR analysis. Genes associated with immunity (knottin) and nervous system (acetylcholine receptor, acetylcholine esterase and sodium channel) were up-regulated while genes involved in metabolism (cytochromep450 and carboxylesterase) were down-regulated, but genes involved in development (ecdysone receptor), reproduction (vitellogenin) and stress (hsp70, hsp90 and shsp) did not change transcription rates. Our results provide information for the understanding of the mechanism of symbiont pathogenic factors for the manipulation of host physiology at the transcription level.