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.

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.

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.

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.

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.

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.

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.