A characterization of cis- and trans-heritability of RNA-Seq-based gene expression.

Insights into individual differences in gene expression and its heritability (h2) can help in understanding pathways from DNA to phenotype. We estimated the heritability of gene expression of 52,844 genes measured in whole blood in the largest twin RNA-Seq sample to date (1497 individuals including 459 monozygotic twin pairs and 150 dizygotic twin pairs) from classical twin modeling and identity-by-state-based approaches. We estimated for each gene h2total, composed of cis-heritability (h2cis, the variance explained by single nucleotide polymorphisms in the cis-window of the gene), and trans-heritability (h2res, the residual variance explained by all other genome-wide variants). Mean h2total was 0.26, which was significantly higher than heritability estimates earlier found in a microarray-based study using largely overlapping (>60%) RNA samples (mean h2 = 0.14, p = 6.15 × 10-258). Mean h2cis was 0.06 and strongly correlated with beta of the top cis expression quantitative loci (eQTL, ρ = 0.76, p < 10-308) and with estimates from earlier RNA-Seq-based studies. Mean h2res was 0.20 and correlated with the beta of the corresponding trans-eQTL (ρ = 0.04, p < 1.89 × 10-3) and was significantly higher for genes involved in cytokine-cytokine interactions (p = 4.22 × 10-15), many other immune system pathways, and genes identified in genome-wide association studies for various traits including behavioral disorders and cancer. This study provides a thorough characterization of cis- and trans-h2 estimates of gene expression, which is of value for interpretation of GWAS and gene expression studies.

Bacterial communities in the small intestine respond differently to those in the caecum and colon in mice fed low- and high-fat diets.

Bacterial communities in the mouse caecum and faeces are known to be altered by changes in dietary fat. The microbiota of the mouse small intestine, by contrast, has not been extensively profiled and it is unclear whether small intestinal bacterial communities shift with dietary fat levels. We compared the microbiota in the small intestine, caecum and colon in mice fed a low-fat (LF) or high-fat (HF) diet using 16S rRNA gene sequencing. The relative abundance of major phyla in the small intestine, Bacteriodetes, Firmicutes and Proteobacteria, was similar to that in the caecum and colon; the relative abundance of Verrucomicrobia was significantly reduced in the small intestine compared to the large intestine. Several genera were uniquely detected in the small intestine and included the aerotolerant anaerobe, Lactobacillus spp. The most abundant genera in the small intestine were accounted for by anaerobic bacteria and were identical to those identified in the large intestine. An HF diet was associated with significant weight gain and adiposity and with changes in the bacterial communities throughout the intestine, with changes in the small intestine differing from those in the caecum and colon. Prominent Gram-negative bacteria including genera of the phylum Bacteroidetes and a genus of Proteobacteria significantly changed in the large intestine. The mechanistic links between these changes and the development of obesity, perhaps involving metabolic endotoxemia, remain to be determined.

Novel genetic locus implicated for HIV-1 acquisition with putative regulatory links to HIV replication and infectivity: a genome-wide association study.

Fifty percent of variability in HIV-1 susceptibility is attributable to host genetics. Thus identifying genetic associations is essential to understanding pathogenesis of HIV-1 and important for targeting drug development. To date, however, CCR5 remains the only gene conclusively associated with HIV acquisition. To identify novel host genetic determinants of HIV-1 acquisition, we conducted a genome-wide association study among a high-risk sample of 3,136 injection drug users (IDUs) from the Urban Health Study (UHS). In addition to being IDUs, HIV-controls were frequency-matched to cases on environmental exposures to enhance detection of genetic effects. We tested independent replication in the Women's Interagency HIV Study (N=2,533). We also examined publicly available gene expression data to link SNPs associated with HIV acquisition to known mechanisms affecting HIV replication/infectivity. Analysis of the UHS nominated eight genetic regions for replication testing. SNP rs4878712 in FRMPD1 met multiple testing correction for independent replication (P=1.38x10(-4)), although the UHS-WIHS meta-analysis p-value did not reach genome-wide significance (P=4.47x10(-7) vs. P<5.0x10(-8)) Gene expression analyses provided promising biological support for the protective G allele at rs4878712 lowering risk of HIV: (1) the G allele was associated with reduced expression of FBXO10 (r=-0.49, P=6.9x10(-5)); (2) FBXO10 is a component of the Skp1-Cul1-F-box protein E3 ubiquitin ligase complex that targets Bcl-2 protein for degradation; (3) lower FBXO10 expression was associated with higher BCL2 expression (r=-0.49, P=8x10(-5)); (4) higher basal levels of Bcl-2 are known to reduce HIV replication and infectivity in human and animal in vitro studies. These results suggest new potential biological pathways by which host genetics affect susceptibility to HIV upon exposure for follow-up in subsequent studies.

De novo insertions and deletions of predominantly paternal origin are associated with autism spectrum disorder.

Whole-exome sequencing (WES) studies have demonstrated the contribution of de novo loss-of-function single-nucleotide variants (SNVs) to autism spectrum disorder (ASD). However, challenges in the reliable detection of de novo insertions and deletions (indels) have limited inclusion of these variants in prior analyses. By applying a robust indel detection method to WES data from 787 ASD families (2,963 individuals), we demonstrate that de novo frameshift indels contribute to ASD risk (OR = 1.6; 95% CI = 1.0-2.7; p = 0.03), are more common in female probands (p = 0.02), are enriched among genes encoding FMRP targets (p = 6 × 10(-9)), and arise predominantly on the paternal chromosome (p < 0.001). On the basis of mutation rates in probands versus unaffected siblings, we conclude that de novo frameshift indels contribute to risk in approximately 3% of individuals with ASD. Finally, by observing clustering of mutations in unrelated probands, we uncover two ASD-associated genes: KMT2E (MLL5), a chromatin regulator, and RIMS1, a regulator of synaptic vesicle release.

Ethanol treatment of lymphoblastoid cell lines from alcoholics and non-alcoholics causes many subtle changes in gene expression.

To elucidate the effects of a controlled exposure to ethanol on gene expression, we studied lymphoblastoid cell lines (LCLs) from 21 alcoholics and 21 controls. We cultured each cell line for 24 h with and without 75 mM ethanol and measured gene expression using microarrays. Differences in expression between LCLs from alcoholics and controls included 13 genes previously identified as associated with alcoholism or related traits, including KCNA3, DICER1, ZNF415, CAT, SLC9A9, and PPARGC1B. The paired design allowed us to detect very small changes due to ethanol treatment: ethanol altered the expression of 37% of the probe sets (51% of the unique named genes) expressed in these LCLs, most by modest amounts. Ninety-nine percent of the named genes expressed in the LCLs were also expressed in brain. Key pathways affected by ethanol include cytokine, TNF, and NFκB signaling. Among the genes affected by ethanol were ANK3, EPHB1, SLC1A1, SLC9A9, NRD1, and SH3BP5, which were reported to be associated with alcoholism or related phenotypes in 2 genome-wide association studies. Genes that either differed in expression between alcoholics and controls or were affected by ethanol exposure are candidates for further study.

Heritability and genomics of gene expression in peripheral blood.

We assessed gene expression profiles in 2,752 twins, using a classic twin design to quantify expression heritability and quantitative trait loci (eQTLs) in peripheral blood. The most highly heritable genes (∼777) were grouped into distinct expression clusters, enriched in gene-poor regions, associated with specific gene function or ontology classes, and strongly associated with disease designation. The design enabled a comparison of twin-based heritability to estimates based on dizygotic identity-by-descent sharing and distant genetic relatedness. Consideration of sampling variation suggests that previous heritability estimates have been upwardly biased. Genotyping of 2,494 twins enabled powerful identification of eQTLs, which we further examined in a replication set of 1,895 unrelated subjects. A large number of non-redundant local eQTLs (6,756) met replication criteria, whereas a relatively small number of distant eQTLs (165) met quality control and replication standards. Our results provide a new resource toward understanding the genetic control of transcription.

Multiple distinct CHRNB3-CHRNA6 variants are genetic risk factors for nicotine dependence in African Americans and European Americans.

AIMS: Studies have shown association between common variants in the α6-ß3 nicotinic receptor subunit gene cluster and nicotine dependence in European ancestry populations. We investigate whether this generalizes to African Americans, whether the association is specific to nicotine dependence and whether this region contains additional genetic contributors to nicotine dependence. DESIGN: We examined consistency of association across studies and race between the α6ß3 nicotinic receptor subunit locus and nicotine, alcohol, marijuana and cocaine dependence in three independent studies. SETTING: United States of America. PARTICIPANTS: European Americans and African Americans from three case-control studies of substance dependence. MEASUREMENTS: Subjects were evaluated using the Semi-Structured Assessment for the Genetics of Alcoholism. Nicotine dependence was determined using the Fagerström Test for Nicotine Dependence. FINDINGS: The single nucleotide polymorphism rs13273442 was associated significantly with nicotine dependence across all three studies in both ancestry groups [odds ratio (OR) = 0.75, P = 5.8 × 10(-4) European Americans; OR = 0.80, P = 0.05 African Americans]. No other substance dependence was associated consistently with this variant in either group. Another SNP in the region, rs4952, remains modestly associated with nicotine dependence in the combined data after conditioning on rs13273442. CONCLUSIONS: The common variant rs13273442 in the CHRNB3-CHNRA6 region is associated significantly with nicotine dependence in European Americans and African Americans across studies recruited for nicotine, alcohol and cocaine dependence. Although these data are modestly powered for other substances, our results provide no evidence that correlates of rs13273442 represent a general substance dependence liability. Additional variants probably account for some of the association of this region to nicotine dependence.

Sex differences in the human peripheral blood transcriptome.

BACKGROUND: Genomes of men and women differ in only a limited number of genes located on the sex chromosomes, whereas the transcriptome is far more sex-specific. Identification of sex-biased gene expression will contribute to understanding the molecular basis of sex-differences in complex traits and common diseases. RESULTS: Sex differences in the human peripheral blood transcriptome were characterized using microarrays in 5,241 subjects, accounting for menopause status and hormonal contraceptive use. Sex-specific expression was observed for 582 autosomal genes, of which 57.7% was upregulated in women (female-biased genes). Female-biased genes were enriched for several immune system GO categories, genes linked to rheumatoid arthritis (16%) and genes regulated by estrogen (18%). Male-biased genes were enriched for genes linked to renal cancer (9%). Sex-differences in gene expression were smaller in postmenopausal women, larger in women using hormonal contraceptives and not caused by sex-specific eQTLs, confirming the role of estrogen in regulating sex-biased genes. CONCLUSIONS: This study indicates that sex-bias in gene expression is extensive and may underlie sex-differences in the prevalence of common diseases.

Copy number variations in 6q14.1 and 5q13.2 are associated with alcohol dependence.

BACKGROUND: Excessive alcohol use is the third leading cause of preventable death and is highly correlated with alcohol dependence, a heritable phenotype. Many genetic factors for alcohol dependence have been found, but many remain unknown. In search of additional genetic factors, we examined the association between Diagnostic and StatisticalManual of Mental Disorders, Fourth Edition (DSM-IV) alcohol dependence and all common copy number variations (CNVs) with good reliability in the Study of Addiction: Genetics and Environment (SAGE). METHODS: All participants in SAGE were interviewed using the Semi-Structured Assessment for the Genetics of Alcoholism, as a part of 3 contributing studies. A total of 2,610 non-Hispanic European American samples were genotyped on the Illumina Human 1M array. We performed CNV calling by CNVPartition, PennCNV, and QuantiSNP, and only CNVs identified by all 3 software programs were examined. Association was conducted with the CNV (as a deletion/duplication) as well as with probes in the CNV region. Quantitative polymerase chain reaction (qPCR) was used to validate the CNVs in the laboratory. RESULTS: CNVs in 6q14.1 (p = 1.04 × 10(-6)) and 5q13.2 (p = 3.37 × 10(-4)) were significantly associated with alcohol dependence after adjusting multiple tests. On chromosome 5q13.2, there were multiple candidate genes previously associated with various neurological disorders. The region on chromosome 6q14.1 is a gene desert that has been associated with mental retardation and language delay. The CNV in 5q13.2 was validated, whereas only a component of the CNV on 6q14.1 was validated by qPCR. Thus, the CNV on 6q14.1 should be viewed with caution. CONCLUSIONS: This is the first study to show an association between DSM-IV alcohol dependence and CNVs. CNVs in regions previously associated with neurological disorders may be associated with alcohol dependence.

Significant effects of antiretroviral therapy on global gene expression in brain tissues of patients with HIV-1-associated neurocognitive disorders.

Antiretroviral therapy (ART) has reduced morbidity and mortality in HIV-1 infection; however HIV-1-associated neurocognitive disorders (HAND) persist despite treatment. The reasons for the limited efficacy of ART in the brain are unknown. Here we used functional genomics to determine ART effectiveness in the brain and to identify molecular signatures of HAND under ART. We performed genome-wide microarray analysis using Affymetrix U133 Plus 2.0 Arrays, real-time PCR, and immunohistochemistry in brain tissues from seven treated and eight untreated HAND patients and six uninfected controls. We also determined brain virus burdens by real-time PCR. Treated and untreated HAND brains had distinct gene expression profiles with ART transcriptomes clustering with HIV-1-negative controls. The molecular disease profile of untreated HAND showed dysregulated expression of 1470 genes at p<0.05, with activation of antiviral and immune responses and suppression of synaptic transmission and neurogenesis. The overall brain transcriptome changes in these patients were independent of histological manifestation of HIV-1 encephalitis and brain virus burdens. Depending on treatment compliance, brain transcriptomes from patients on ART had 83% to 93% fewer dysregulated genes and significantly lower dysregulation of biological pathways compared to untreated patients, with particular improvement indicated for nervous system functions. However a core of about 100 genes remained similarly dysregulated in both treated and untreated patient brain tissues. These genes participate in adaptive immune responses, and in interferon, cell cycle, and myelin pathways. Fluctuations of cellular gene expression in the brain correlated in Pearson's formula analysis with plasma but not brain virus burden. Our results define for the first time an aberrant genome-wide brain transcriptome of untreated HAND and they suggest that antiretroviral treatment can be broadly effective in reducing pathophysiological changes in the brain associated with HAND. Aberrantly expressed transcripts common to untreated and treated HAND may contribute to neurocognitive changes defying ART.

Copy number variation accuracy in genome-wide association studies.

BACKGROUND/AIM: Copy number variations (CNVs) are a major source of alterations among individuals and are a potential risk factor in many diseases. Numerous diseases have been linked to deletions and duplications of these chromosomal segments. Data from genome-wide association studies and other microarrays may be used to identify CNVs by several different computer programs, but the reliability of the results has been questioned. METHODS: To help researchers reduce the number of false-positive CNVs that need to be followed up with laboratory testing, we evaluated the relative performance of CNVPartition, PennCNV and QuantiSNP, and developed a statistical method for estimating sensitivity and positive predictive values of CNV calls and tested it on 96 duplicate samples in our dataset. RESULTS: We found that the positive predictive rate increases with the number of probes in the CNV and the size of the CNV, with the highest positive predicted rates in CNVs of at least 500 kb and at least 100 probes. Our analysis also indicates that identifying CNVs reported by multiple programs can greatly improve the reproducibility rate and the positive predicted rate. CONCLUSION: Our methods can be used by investigators to identify CNVs in genome-wide data with greater reliability.

Gene expression profiles of HIV-1-infected glia and brain: toward better understanding of the role of astrocytes in HIV-1-associated neurocognitive disorders.

Astrocytes are the major cellular component of the central nervous system (CNS), and they play multiple roles in brain development, normal brain function, and CNS responses to pathogens and injury. The functional versatility of astrocytes is linked to their ability to respond to a wide array of biological stimuli through finely orchestrated changes in cellular gene expression. Dysregulation of gene expression programs, generally by chronic exposure to pathogenic stimuli, may lead to dysfunction of astrocytes and contribute to neuropathogenesis. Here, we review studies that employ functional genomics to characterize the effects of HIV-1 and viral pathogenic proteins on cellular gene expression in astrocytes in vitro. We also present the first microarray analysis of primary mouse astrocytes exposed to HIV-1 in culture. In spite of different experimental conditions and microarray platforms used, comparison of the astrocyte array data sets reveals several common gene-regulatory changes that may underlie responses of these cells to HIV-1 and its proteins. We also compared the transcriptional profiles of astrocytes with those obtained in analyses of brain tissues of patients with HIV-1 dementia and macaques infected with simian immunodeficiency virus (SIV). Notably, many of the gene characteristics of responses to HIV-1 in cultured astrocytes were also altered in HIV-1 or SIV-infected brains. Functional genomics, in conjunction with other approaches, may help clarify the role of astrocytes in HIV-1 neuropathogenesis.

Effect of the methoxychlor metabolite HPTE on the rat ovarian granulosa cell transcriptome in vitro.

Ovarian granulosa cells play a central role in steroidogenesis, which is critical for female reproduction. Follicle-stimulating hormone (FSH) promotes cyclic adenosine monophosphate (cAMP)-mediated signaling to regulate granulosa cell steroidogenesis. We have shown previously that 2,2-bis-(p-hydroxyphenyl)-1,1,1-trichloroethane (HPTE) inhibits FSH- and dibutyryl cAMP-stimulated steroidogenesis and affects the messenger RNA levels of steroidogenic pathway enzymes in rat granulosa cells. However, HPTE showed a differential effect in FSH- and cAMP-stimulated cells in that HPTE more completely blocked FSH- when compared to cAMP-driven steroidogenesis. The objective of this study was to analyze the effects of HPTE on global gene expression profiles in untreated granulosa cells and those challenged with FSH or cAMP. Granulosa cells from immature rats were cultured with 0, 1, 5, or 10 microM HPTE in the presence or absence of either 3 ng FSH/ml or 1mM cAMP for 48 h. Total RNA was isolated for real-time quantitative PCR and microarray analysis using the GeneChip Rat Genome 230 2.0 and ArrayAssist Microarray Suite. An investigation of changes in gene expression across all HPTE treatments showed that HPTE altered more genes in FSH- (approximately 670 genes) than in cAMP-stimulated cells (approximately 366 genes). Analysis confirmed that HPTE more effectively inhibited FSH- than cAMP-induced steroid pathway gene expression and steroidogenesis. Furthermore, expression patterns of novel genes regulating signal transduction, transport, cell cycle, adhesion, differentiation, motility and growth, apoptosis, development, and metabolism were all altered by HPTE. This study further established that HPTE exerts differential effects within the granulosa cell steroidogenic pathway and revealed that these effects include broader changes in gene expression.

Inhibition of autoregulated TGFbeta signaling simultaneously enhances proliferation and differentiation of kidney epithelium and promotes repair following renal ischemia.

We studied autocrine transforming growth factor (TGF)beta signaling in kidney epithelium. Cultured proximal tubule cells showed regulated signaling that was high during log-phase growth, low during contact-inhibited differentiation, and rapidly increased during regeneration of wounded epithelium. Autoregulation of signaling correlated with TGFbeta receptor and Smad7 levels, but not with active TGFbeta, which was barely measurable in the growth medium. Confluent differentiated cells with low receptor and high Smad7 levels exhibited blunted responses to saturating concentrations of exogenously provided active TGFbeta, suggesting that TGFbeta signaling homeostasis was achieved by cell density-dependent modulation of signaling intermediates. Antagonism of Alk5 kinase, the TGFbeta type I receptor, dramatically accelerated the induction of differentiation in sparse, proliferating cultures and permitted better retention of differentiated features in regenerating cells of wounded, confluent cultures. Alk5 antagonism accelerated the differentiation of cells in proximal tubule primary cultures while simultaneously increasing their proliferation. Consequently, Alk5-inhibited primary cultures formed confluent, differentiated monolayers faster than untreated cultures. Furthermore, treatment with an Alk5 antagonist promoted kidney repair reflected by increased tubule differentiation and decreased tubulo-interstitial pathology during the recovery phase following ischemic injury in vivo. Our results show that autocrine TGFbeta signaling in proliferating proximal tubule cells exceeds the levels that are necessary for physiological regeneration. To that end, TGFbeta signaling is redundant and maladaptive during tubule repair by epithelial regeneration.

Oligonucleotide microarrays reveal regulated genes related to inward arterial remodeling induced by urokinase plasminogen activator.

Accumulating evidence suggests that urokinase plasminogen activator (uPA) is involved in vascular remodeling and lumen stenosis after angioplasty and stenting. We have shown previously that increased uPA expression greatly promotes neointima formation and inward arterial remodeling after balloon injury. To evaluate the role of inflammation in early mechanisms responsible for inward arterial remodeling induced by uPA and elucidate the mechanisms of remodeling, we characterized changes in the expression profiles of 8,799 genes in injured rat carotid arteries 1 and 4 days after recombinant uPA treatment compared to vehicle. We used a standard model of the balloon catheter injury of the rat carotid followed by periadventitial application to the injured vessel of either uPA dissolved in Pluronic gel, or plain gel. Vessels were harvested and analyzed by immunohistochemistry, morphometry, microarray gene expression profiling and quantitative RT-PCR. Periadventitial application of uPA significantly reduced lumen size and vessel area encompassed by the external elastic lamina at both 1 and 4 days after treatment. Inflammatory cells accumulated in the arterial adventitia at both 1 and 4 days after uPA treatment. On the 4th day, increases in the areas and arterial cell numbers of all arterial layers were found. Among 79 differentially expressed known genes 1 day after uPA application, 12 proinflammatory genes, including TNF-alpha and TACE, and 15 genes related to mitochondrial metabolism and oxidative stress regulation were identified. Four days after injury in uPA-treated arteries, 3 proinflammatory and 2 oxidation-related genes were differentially expressed. We conclude that uPA likely promotes inward arterial remodeling by regulating oxidative stress and inflammation after arterial injury.

Big results from small samples: evaluation of amplification protocols for gene expression profiling.

Microarrays have revolutionized many areas of biology due to our technical ability to quantify tens of thousands of transcripts within a single experiment. However, there are still many areas that cannot benefit from this technology due to the amount of biological material needed for microarray analysis. In response to this demand, chemistries have been developed that boast the capability of generating targets from nanogram amounts of total RnA, reflecting minimal amounts of biological material, on the order of several hundred or thousand cells. Herein, we describe the evaluation of four chemistries for RnA amplification in terms of reproducibility, sensitivity, accuracy, and comparability to results from a single round of T7 amplification. No evidence for false-positive measurements of differential expression was observed. In contrast, clear differences between chemistries in sensitivity and accuracy were detected. PCR validation showed an interaction of probe sequence on the array and target labeling chemistry, resulting in a chemistry-dependent probe set sensitivity varying over an order of magnitude.

Transcriptional profiling of C57 and DBA strains of mice in the absence and presence of morphine.

BACKGROUND: The mouse C57BL/6 (C57) and DBA/2J (DBA) inbred strains differ substantially in many aspects of their response to drugs of abuse. The development of microarray analyses represents a genome-wide method for measuring differences across strains, focusing on expression differences. In the current study, we carried out microarray analysis in C57 and DBA mice in the nucleus accumbens of drug-naïve and morphine-treated animals. RESULTS: We identified mRNAs with altered expression between the two strains. We validated the mRNA expression changes of several such mRNAs, including Gnb1, which has been observed to be regulated by several drugs of abuse. In addition, we validated alterations in the enzyme activity of one mRNA product, catechol-O-methyltransferase (Comt). Data mining of expression and behavioral data indicates that both Gnb1 and Comt expression correlate with aspects of drug response in C57/DBA recombinant inbred strains. Pathway analysis was carried out to identify pathways showing significant alterations as a result of treatment and/or due to strain differences. These analyses identified axon guidance genes, particularly the semaphorins, as showing altered expression in the presence of morphine, and plasticity genes as showing altered expression across strains. Pathway analysis of genes showing strain by treatment interaction suggest that the phosphatidylinositol signaling pathway may represent an important difference between the strains as related to morphine exposure. CONCLUSION: mRNAs with differing expression between the two strains could potentially contribute to strain-specific responses to drugs of abuse. One such mRNA is Comt and we hypothesize that altered expression of Comt may represent a potential mechanism for regulating the effect of, and response to, multiple substances of abuse. Similarly, a role for Gnb1 in responses to multiple drugs of abuse is supported by expression data from our study and from other studies. Finally, the data support a role for semaphorin signaling in morphine effects, and indicate that altered expression of genes involved in phosphatidylinositol signaling and plasticity might also affect the altered drug responses in the two strains.

TCDD-induced alterations in gene expression profiles of the developing mouse paw do not influence morphological differentiation of this potential target tissue.

The aryl-hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that mediates the toxicity of certain halogenated aromatic hydrocarbons including 2,3,7,8-tetra-chlorodibenzo-p-dioxin (TCDD). These compounds are potent developmental toxicants that can alter gene expression and disrupt processes of proliferation and differentiation. It has not yet been determined which tissues during development are most sensitive to these compounds, nor which genes are directly associated with the toxicities. We developed a transgenic (TG) mouse model to delineate the temporal and spatial context of transcriptionally active AhR by utilizing a dioxin responsive element-linked LacZ reporter system. The present study focuses on the pattern of TCDD-induced transgene expression localized to the footpad and digits of the paws between gestational days (GD) 13 and 18. Paw morphology was evaluated at several developmental stages following TCDD exposure. Gene expression profiles acquired by microarray technology were evaluated in the paws of fetuses exposed at GD 14.5. The results showed that TCDD exposure in utero induced LacZ expression in the developing paws. This expression appeared to be localized to the mesenchymal cell layer. Gross morphological changes were not observed in the paws prior to or after birth following TCDD exposure in utero. However, significant alterations in gene expression profiles in the developing paws were observed at 24 h following TCDD exposure in utero. These results indicate that the developing paw is a target tissue of TCDD in terms of altered gene expression, further validating the use of this AhR responsive reporter gene TG mouse model in studying AhR ligand-mediated responsiveness. However, the linkage of these changes to detectable biological outcomes in the paw remains unclear.