Search

Exposure to environmental mutagens increases the risk of cancer and genetic disorders. We used Duplex Sequencing (DS), a high-accuracy error-corrected sequencing technology, to analyze mutation induction across twenty 2.4 kb intergenic and genic targets in the bone marrow of MutaMouse males exposed to benzo(a)pyrene (BaP), a widespread environmental pollutant. DS revealed a linear dose-related induction of mutations across all targets with low intra-group variability. Heterochromatic and intergenic regions exhibited the highest mutation frequencies (MF). C:G > A:T transversions at CCA, CCC and GCC trinucleotides were enriched in BaP-exposed mice consistent with the known etiology of BaP mutagenesis. However, GC-content had no effect on mutation susceptibility. A positive correlation was observed between DS and the “gold-standard” transgenic rodent gene mutation assay. Overall, we demonstrate that DS is a promising approach to study in vivo mutagenesis and yields critical insight into the genomic features governing mutation susceptibility, spectrum, and variability across the genome.

Background The transgenic rodent mutation reporter assay provides an efficient approach to identify mutagenic agents in vivo. A major advantage of this assay is that mutant reporter transgenes can be sequenced to provide information on the mode of action of a mutagen and to identify clonally expanded mutations. However, conventional DNA sequence analysis is laborious and expensive for long transgenes, such as lacZ (3096 bp), and is not normally implemented in routine screening. Methods We developed a high-throughput next-generation sequencing (NGS) approach to simultaneously sequence large numbers of barcoded mutant lacZ transgenes from different animals. We collected 3872 mutants derived from the bone marrow DNA of six Muta™Mouse males exposed to the well-established mutagen benzo[a]pyrene (BaP) and six solvent-exposed controls. Mutants within animal samples were pooled, barcoded, and then sequenced using NGS. Results We identified 1652 mutant sequences from 1006 independent mutations that underwent clonal expansion. This deep sequencing analysis of mutation spectrum demonstrated that BaP causes primarily guanine transversions (e.g. G:C → T:A), which is highly consistent with previous studies employing Sanger sequencing. Furthermore, we identified novel mutational hotspots in the lacZ transgene that were previously uncharacterized by Sanger sequencing. Deep sequencing also allowed for an unprecedented ability to correct for clonal expansion events, improving the sensitivity of the mutation reporter assay by 50 %. Conclusion These results demonstrate that the high-throughput nature and reduced costs offered by NGS provide a sensitive and fast approach for elucidating and comparing mutagenic mechanisms of various agents among tissues and enabling improved evaluation of genotoxins. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2004-4) contains supplementary material, which is available to authorized users.

Background Microarray normalizations typically apply methods that assume absence of global transcript shifts, or absence of changes in internal control features such as housekeeping genes. These normalization approaches are not appropriate for focused arrays with small sets of genes where a large portion may be expected to change. Furthermore, many microarrays lack control features that can be used for quality assurance (QA). Here, we describe a novel external control series integrated with a design feature that addresses the above issues. Results An EC dilution series that involves spike-in of a single concentration of the A. thaliana chlorophyll synthase gene to hybridize against spotted dilutions (0.000015 to 100 μM) of a single complimentary oligonucleotide representing the gene was developed. The EC series is printed in duplicate within each subgrid of the microarray and covers the full range of signal intensities from background to saturation. The design and placement of the series allows for QA examination of frequently encountered problems in hybridization (e.g., uneven hybridizations) and printing (e.g., cross-spot contamination). Additionally, we demonstrate that the series can be integrated with a LOWESS normalization to improve the detection of differential gene expression (improved sensitivity and predictivity) over LOWESS normalization on its own. Conclusion The quality of microarray experiments and the normalization methods used affect the ability to measure accurate changes in gene expression. Novel methods are required for normalization of small focused microarrays, and for incorporating measures of performance and quality. We demonstrate that dilution of oligonucleotides on the microarray itself provides an innovative approach allowing the full dynamic range of the scanner to be covered with a single gene spike-in. The dilution series can be used in a composite normalization to improve detection of differential gene expression and to provide quality control measures.

Background Disruption of thyroid hormone signalling can alter growth, development and energy metabolism. Thyroid hormones exert their effects through interactions with thyroid receptors that directly bind thyroid response elements and can alter transcriptional activity of target genes. The effects of short-term thyroid hormone perturbation on hepatic mRNA transcription in juvenile mice were evaluated, with the goal of identifying genes containing active thyroid response elements. Thyroid hormone disruption was induced from postnatal day 12 to 15 by adding goitrogens to dams' drinking water (hypothyroid). A subgroup of thyroid hormone-disrupted pups received intraperitoneal injections of replacement thyroid hormones four hours prior to sacrifice (replacement). An additional group received only thyroid hormones four hours prior to sacrifice (hyperthyroid). Hepatic mRNA was extracted and hybridized to Agilent mouse microarrays. Results Transcriptional profiling enabled the identification of 28 genes that appeared to be under direct thyroid hormone-regulation. The regulatory regions of the genome adjacent to these genes were examined for half-site sequences that resemble known thyroid response elements. A bioinformatics search identified 33 thyroid response elements in the promoter regions of 13 different genes thought to be directly regulated by thyroid hormones. Thyroid response elements found in the promoter regions of Tor1a, 2310003H01Rik, Hect3d and Slc25a45 were further validated by confirming that the thyroid receptor is associated with these sequences in vivo and that it can bind directly to these sequences in vitro. Three different arrangements of thyroid response elements were identified. Some of these thyroid response elements were located far up-stream (> 7 kb) of the transcription start site of the regulated gene. Conclusions Transcriptional profiling of thyroid hormone disrupted animals coupled with a novel bioinformatics search revealed new thyroid response elements associated with genes previously unknown to be responsive to thyroid hormone. The work provides insight into thyroid response element sequence motif characteristics.

Background Although analysis of microRNAs (miRNAs) by DNA microarrays is gaining in popularity, these new technologies have not been adequately validated. We examined within and between platform reproducibility of four miRNA array technologies alongside TaqMan PCR arrays. Results Two distinct pools of reference materials were selected in order to maximize differences in miRNA content. Filtering for miRNA that yielded signal above background revealed 54 miRNA probes (matched by sequence) across all platforms. Using this probeset as well as all probes that were present on an individual platform, within-platform analyses revealed Spearman correlations of >0.9 for most platforms. Comparing between platforms, rank analysis of the log ratios of the two reference pools also revealed high correlation (range 0.663-0.949). Spearman rank correlation and concordance correlation coefficients for miRNA arrays against TaqMan qRT-PCR arrays were similar for all of the technologies. Platform performances were similar to those of previous cross-platform exercises on mRNA and miRNA microarray technologies. Conclusions These data indicate that miRNA microarray platforms generated highly reproducible data and can be recommended for the study of changes in miRNA expression.

Background Thyroid hormones play an essential role in early vertebrate development as well as other key processes. One of its modes of action is to bind to the thyroid hormone receptor (TR) which, in turn, binds to thyroid response elements (TREs) in promoter regions of target genes. The sequence motif for TREs remains largely undefined as does the precise chromosomal location of the TR binding sites. A chromatin immunoprecipitation on microarray (ChIP-chip) experiment was conducted using mouse cerebellum post natal day (PND) 4 and PND15 for the thyroid hormone receptor (TR) beta 1 to map its binding sites on over 5000 gene promoter regions. We have performed a detailed computational analysis of these data. Results By analysing a recent spike-in study, the optimal normalization and peak identification approaches were determined for our dataset. Application of these techniques led to the identification of 211 ChIP-chip peaks enriched for TR binding in cerebellum samples. ChIP-PCR validation of 25 peaks led to the identification of 16 true positive TREs. Following a detailed literature review to identify all known mouse TREs, a position weight matrix (PWM) was created representing the classic TRE sequence motif. Various classes of promoter regions were investigated for the presence of this PWM, including permuted sequences, randomly selected promoter sequences, and genes known to be regulated by TH. We found that while the occurrence of the TRE motif is strongly correlated with gene regulation by TH for some genes, other TH-regulated genes do not exhibit an increased density of TRE half-site motifs. Furthermore, we demonstrate that an increase in the rate of occurrence of the half-site motifs does not always indicate the specific location of the TRE within the promoter region. To account for the fact that TR often operates as a dimer, we introduce a novel dual-threshold PWM scanning approach for identifying TREs with a true positive rate of 0.73 and a false positive rate of 0.2. Application of this approach to ChIP-chip peak regions revealed the presence of 85 putative TREs suitable for further in vitro validation. Conclusions This study further elucidates TRβ gene regulation in mouse cerebellum, with 211 promoter regions identified to bind to TR. While we have identified 85 putative TREs within these regions, future work will study other mechanisms of action that may mediate the remaining observed TR-binding activity.