Differences In Transcriptional Responses Research Articles

Hepatic Transcriptomic Responses to Ethinylestradiol in Embryonic Japanese Quail and Double-Crested Cormorant.

Understanding species differences in sensitivity to toxicants is a critical issue in ecotoxicology. We recently established that double-crested cormorant (DCCO) embryos are more sensitive than Japanese quail (JQ) to the developmental effects of ethinylestradiol (EE2). We explored how this difference in sensitivity between species is reflected at a transcriptomic level. The EE2 was dissolved in dimethyl sulfoxide and injected into the air cell of eggs prior to incubation at nominal concentrations of 0, 3.33, and 33.3 µg/g egg weight. At midincubation (JQ 9 days; DCCO 16 days), livers were collected from five embryos/treatment group for RNA sequencing. Data were processed and analyzed using EcoOmicsAnalyst and ExpressAnalyst. The EE2 exposure dysregulated 238 and 1,987 genes in JQ and DCCO, respectively, with 78 genes in common between the two species. These included classic biomarkers of estrogen exposure such as vitellogenin and apovitellenin. We also report DCCO-specific dysregulation of Phase I/II enzyme-coding genes and species-specific transcriptional ontogeny of vitellogenin-2. Twelve Kyoto Encyclopedia of Genes and Genomes pathways and two EcoToxModules were dysregulated in common in both species including the peroxisome proliferator-activated receptor (PPAR) signaling pathway and fatty acid metabolism. Similar to previously reported differences at the organismal level, DCCO were more responsive to EE2 exposure than JQ at the gene expression level. Our description of differences in transcriptional responses to EE2 in early life stage birds may contribute to a better understanding of the molecular basis for species differences. Environ Toxicol Chem 2024;43:772-783. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Aedes aegypti gut transcriptomes respond differently to microbiome transplants from field-caught or laboratory-reared mosquitoes.

The mosquito microbiome is critical for host development and plays a major role in many aspects of mosquito biology. While the microbiome is commonly dominated by a small number of genera, there is considerable variation in composition among mosquito species, life stages, and geography. How the host controls and is affected by this variation is unclear. Using microbiome transplant experiments, we asked whether there were differences in transcriptional responses when mosquitoes of different species were used as microbiome donors. We used microbiomes from four different donor species spanning the phylogenetic breadth of the Culicidae, collected either from the laboratory or the field. We found that when recipients received a microbiome from a donor reared in the laboratory, the response was remarkably similar regardless of donor species. However, when the donor had been collected from the field, many more genes were differentially expressed. We also found that while the transplant procedure did have some effect on the host transcriptome, this is likely to have had a limited effect on mosquito fitness. Overall, our results highlight the possibility that variation in mosquito microbiome communities is associated with variability in host-microbiome interactions and further demonstrate the utility of the microbiome transplantation technique for investigating host-microbe interactions in mosquitoes.

BHLH heterodimer complex variations regulate cell proliferation activity in the meristems of Arabidopsis thaliana

Root, shoot, and lateral meristems are the main regions of cell proliferation in plants. It has been proposed that meristems might have evolved dedicated transcriptional networks to balance cell proliferation. Here, we show that basic helix-loop-helix (bHLH) transcription factor heterodimers formed by members of the TARGET OF MONOPTEROS5 (TMO5) and LONESOME HIGHWAY (LHW) subclades are general regulators of cell proliferation in all meristems. Yet, genetics and expression analyses suggest specific functions of these transcription factors in distinct meristems, possibly due to their expression domains determining heterodimer complex variations within meristems, and to a certain extent to the absence of some of them in a given meristem. Target gene specificity analysis for heterodimer complexes focusing on the LONELY GUY gene targets further suggests differences in transcriptional responses through heterodimer diversification that could allow a common bHLH heterodimer complex module to contribute to cell proliferation control in multiple meristems.

Advances in Vitamin D Receptor Function and Evolution Based on the 3D Structure of the Lamprey Ligand-Binding Domain.

1α,25-dihydroxyvitamin D3 (1,25D3) regulates many physiological processes in vertebrates by binding to the vitamin D receptor (VDR). Phylogenetic analysis indicates that jawless fishes are the most basal vertebrates exhibiting a VDR gene. To elucidate the mechanism driving VDR activation during evolution, we determined the crystal structure of the VDR ligand-binding domain (LBD) complex from the basal vertebratePetromyzon marinus, sea lamprey (lVDR). Comparison of three-dimensional crystal structures of the lVDR-1,25D3 complex with higher vertebrate VDR-1,25D3 structures suggests that 1,25D3 binds to lVDR similarly to human VDR, but with unique features for lVDR around linker regions between H11 and H12 and between H9 and H10. These structural differences may contribute to the marked species differences in transcriptional responses. Furthermore, residue co-evolution analysis of VDR across vertebrates identifies amino acid positions in H9 and the large insertion domain VDR LBD specific as correlated.

Comparative Transcriptional Analysis of Patient Responders Versus Non-Responders to Glucocorticoid Treatment for Bronchopulmonary Dysplasia

Bronchopulmonary Dysplasia (BPD) is a common heterogeneous lung disease that can result from preterm birth at less than 28-weeks gestation, prenatal and postnatal inflammatory insults, ventilator associated lung injury, and oxygen-related injury. Synthetic glucocorticoids (sGCs) are commonly used pre- and postnatally to treat inflammation and improve lung physiology. Clinical responses to sGC therapy for BPD vary in patients. We hypothesize that genetic background differences in transcriptional response to sGC therapy dictate the efficacy in infants with BPD. Identifying pathways and genes that mediate these differences will allow prospective determination of which infants would respond to sGC treatment.26 preterm infants that received sGC treatment for BPD were identified. Respiratory Severity Score (RSS), an indication of BPD severity, was measured at diagnosis, 4 days, and 7 days post-sGC treatment. Patients were stratified into Responders versus Non-Responders by improvement in respiratory function after treatment. Changes in RSS were used to discriminate Responders (R >3 decrease in RSS) to treatment from Non-Responders (NR <3 decrease). 13 Responders and 13 Non-Responders were selected. They included 7 females and 19 males with an average gestational age of 24.3 weeks, and were 46% Caucasian, 31% African American, 19% Hispanic, and 4% other. 100µL of blood was collected before and after seven days of a dexamethasone treatment course.To examine differences in transcription response between Responders (n = 13) and Non-Responders (n= 13), RNA was isolated and analyzed using the Clariom S Human Transcriptome Affymetrix array. 21,500 expressed genes were profiled. Results: were imported into the Transcriptome Analysis Console (TAC) software, and genes with a significant difference (fold change >1.48 or < -1.48 and p-value <0.05) in Responders and Non-Responders were identified. Of those, 133 genes were upregulated and 74 downregulated. Ingenuity Pathway Analysis (IPA) was used to identify signaling pathways and disease processes that were uniquely altered in Responders versus Non-Responders. Non-Responders showed significant activation of neuroinflammatory signaling pathways, degranulation pathways, and lymphocyte activation disease pathways. Target genes in the top dysregulated pathways were evaluated using quantitative Polymerase Chain Reaction (qPCR). Expression changes in Matrix Metalloproteinase-25, Interleukin-12 Receptor beta, and Microsomal Glutathione Transferase-1, key mediators of inflammation, were validated in independent studies using qPCR. While response to systemic glucocorticoids in neonates with BPD is variable, these studies identified pathways that are altered in Responders versus Non-Responders and are a step towards developing pre-screening tools to stratify infants for response to sGC BPD therapy.

Methionine-chelated Zn promotes anabolism by integrating mTOR signal and autophagy pathway in juvenile yellow catfish

BackgroundAmino acid-chelated zinc (Zn) can increase anabolism of animals. However, the underlying mechanisms are unclear. We aimed to examine how autophagy impact anabolism following a diet containing methionine-chelated Zn (ZnMet) compared with inorganic Zn (ZnSO4). MethodsYellow catfish (weight: 4.02 ± 0.08 g) were fed two diets containing ZnSO4 or ZnMet for 8 wk. The differences in transcriptional responses and corresponding biological profiles were compared between the livers of fish fed the two Zn sources of diets. Hepatocytes of yellow catfish were incubated for 48 h in medium containing ZnSO4 (10 μM ZnSO4) or ZnMet (10 μM ZnMet) after 2 h pretreated with or without pathway inhibitors. Intracellular Zn, TG and protein contents, lipid droplet and autophagic vesicles were detected. Ultrastructural observation, enzymatic activities, qPCR assays, western blot and immunofluorescence analysis were conducted. ResultsZnMet up-regulated the expression of genes associated with anabolism and autophagy. The differentially expressed genes (DEG) analysis indicated that both mTOR and autophagy pathways were activated. ZnMet-induced activation of autophagy was mTOR-independent. In this process, forkhead box class O was deacetylated and activated, and induced autophagy, which provided substrates for energy generation. ConclusionsZnMet increased anabolism through integrating mTOR signal and autophagy pathway in yellow catfish. The present study unravels a novel mechanism for amino acid-chelated minerals improving anabolism.

Integrating LCM-Based Spatio-Temporal Transcriptomics Uncovers Conceptus and Endometrial Luminal Epithelium Communication that Coordinates the Conceptus Attachment in Pigs.

Attachment of conceptus to the endometrial luminal epithelium (LE) is a critical event for early placentation in Eutheria. Since the attachment occurs at a particular site within the uterus, a coordinated communication between three spatially distinct compartments (conceptus and endometrial LE from two anatomical regions of the uterus to which conceptus attaches and does not attach) is essential but remains to be fully characterized. Using the laser capture microdissection (LCM) technique, we firstly developed an approach that can allow us to pair the pig conceptus sample with its nearby endometrial epithelium sample without losing the native spatial information. Then, a comprehensive spatio-temporal transcriptomic profile without losing the original conceptus-endometrium coordinates was constructed. The analysis shows that an apparent difference in transcriptional responses to the conceptus exists between the endometrial LE from the two anatomically distinct regions in the uterus. In addition, we identified the communication pathways that link the conceptus and endometrial LE and found that these pathways have important roles in conceptus attachment. Furthermore, a number of genes whose expression is spatially restricted in the two different anatomical regions within the uterus were characterized for the first time and two of them (SULT2A1 and MEP1B) may cooperatively contribute to establish conceptus attachment in pigs. The results from our study have implications in understanding of conceptus/embryo attachment in pigs and other large polytocous species.

Different transcriptional responses of haploid and diploid S. cerevisiae strains to changes in cofactor preference of XR

BackgroundXylitol accumulation is a major barrier for efficient ethanol production through heterologous xylose reductase-xylitol dehydrogenase (XR-XDH) pathway in recombinant Saccharomyces cerevisiae. Mutated NADH-preferring XR is usually employed to alleviate xylitol accumulation. However, it remains unclear how mutated XR affects the metabolic network for xylose metabolism. In this study, haploid and diploid strains were employed to investigate the transcriptional responses to changes in cofactor preference of XR through RNA-seq analysis during xylose fermentation.ResultsFor the haploid strains, genes involved in xylose-assimilation (XYL1, XYL2, XKS1), glycolysis, and alcohol fermentation had higher transcript levels in response to mutated XR, which was consistent with the improved xylose consumption rate and ethanol yield. For the diploid strains, genes related to protein biosynthesis were upregulated while genes involved in glyoxylate shunt were downregulated in response to mutated XR, which might contribute to the improved yields of biomass and ethanol. When comparing the diploids with the haploids, genes involved in glycolysis and MAPK signaling pathway were significantly downregulated, while oxidative stress related transcription factors (TFs) were significantly upregulated, irrespective of the cofactor preference of XR.ConclusionsOur results not only revealed the differences in transcriptional responses of the diploid and haploid strains to mutated XR, but also provided underlying basis for better understanding the differences in xylose metabolism between the diploid and haploid strains.

Effects of set cathode potentials on microbial electrosynthesis system performance and biocathode methanogen function at a metatranscriptional level

Microbial electrosynthesis exploits the catalytic activity of microorganisms to utilize a cathode as an electron donor for reducing waste CO2 to valuable fuels and chemicals. Electromethanogenesis is the process of CO2 reduction to CH4 catalyzed by methanogens using the cathode directly as a source of electrons or indirectly via H2. Understanding the effects of different set cathode potentials on the functional dynamics of electromethanogenic communities is crucial for the rational design of cathode materials. Replicate enriched electromethanogenic communities were subjected to different potentials (− 1.0 V and − 0.7 V vs. Ag/AgCl) and the potential-induced changes were analyzed using a metagenomic and metatranscriptomic approach. The most abundant and transcriptionally active organism on the biocathodes was a novel species of Methanobacterium sp. strain 34x. The cathode potential-induced changes limited electron donor availability and negatively affected the overall performance of the reactors in terms of CH4 production. Although high expression of key genes within the methane and carbon metabolism pathways was evident, there was no significant difference in transcriptional response to the different set potentials. The acetyl-CoA decarbonylase/synthase (ACDS) complex were the most highly expressed genes, highlighting the significance of carbon assimilation under limited electron donor conditions and its link to the methanogenesis pathway.

Being African American and Rural: A Double Jeopardy From COVID-19.

Being African American and Rural: A Double Jeopardy From COVID-19.

Effects of 7°C environmental temperature acclimation during a 3-week training period.

Cold environmental temperatures during exercise and recovery alter the acute response to cellular signaling and training adaptations. Approximately 3 wk is required for cold temperature acclimation to occur. To determine the impact of cold environmental temperature on training adaptations, fitness measurements, and aerobic performance, two groups of 12 untrained male subjects completed 1 h of cycling in 16 temperature acclimation sessions in either a 7°C or 20°C environmental temperature. Fitness assessments before and after acclimation occurred at standard room temperature. Muscle biopsies were taken from the vastus lateralis muscle before and after training to assess molecular markers related to mitochondrial development. Peroxisome proliferator-activated receptor-γ coactivator 1α (PGC-1α) mRNA was higher in 7°C than in 20°C in response to acute exercise before training (P = 0.012) but not after training (P = 0.813). PGC-1α mRNA was lower after training (P < 0.001). BNIP3 was lower after training in the 7°C than in the 20°C group (P = 0.017) but not before training (P = 0.549). No other differences occurred between temperature groups in VEGF, ERRα, NRF1, NRF2, TFAM, PINK1, Parkin, or BNIP3L mRNAs (P > 0.05). PGC-1α protein and mtDNA were not different before training, after training, or between temperatures (P > 0.05). Cycling power increased during the daily training (P < 0.001) but was not different between temperatures (P = 0.169). V̇o2peak increased with training (P < 0.001) but was not different between temperature groups (P = 0.460). These data indicate that a 3-wk period of acclimation/training in cold environmental temperatures alters PGC-1α gene expression acutely but this difference is not manifested in a greater increase in V̇o2peak and is dissipated as acclimation takes place.NEW & NOTEWORTHY This study examines the adaptive response of cellular signaling during exercise in cold environmental temperatures. We demonstrate that peroxisome proliferator-activated receptor-γ coactivator 1α mRNA is different between cold and room temperature environments before training but after training this difference no longer exists. This initial difference in transcriptional response between temperatures does not lead to differences in performance measures or increases in protein or mitochondria.

Identifying qualitative differences in PPARα signaling networks in human and rat hepatocytes and their significance for next generation chemical risk assessment methods

In this paper, we evaluate the PPARα signaling network in rats, examining transcriptional responses in primary hepatocytes exposed to a PPARα specific ligand, GW7647. These transcriptomic studies were complemented with ChIP-seq studies of PPARα binding and transcription binding motif identification for PPARα responsive genes. We also conducted a limited study of GW7647 dosing the in intact rat to examine differences in transcriptional responses for primary hepatocytes in vitro and in the intact liver. The rat network has a much larger number of down-regulated genes and pathways than we had found in the human and the PPARα binding motifs in rat differed for upregulated and down regulated genes. Based on these results and comparison with our previous work with the human PPARα signaling network, we identified qualitative differences in the transcriptional networks controlled by PPARα activation in the two species that provide an explanation of the interspecies differences in the responses of humans and rodents to GW7647 and likely to other PPARα agonists. These studies also allow some observations on the manner in which in vitro, fit-for-purpose assays in human hepatocytes could form the basis for risk assessment without recourse to in-life studies in rodents or other test species.

Unique Transcriptional Architecture in Airway Epithelial Cells and Macrophages Shapes Distinct Responses following Influenza Virus Infection Ex Vivo.

Airway epithelial cells and macrophages differ markedly in their responses to influenza A virus (IAV) infection. To investigate transcriptional responses underlying these differences, purified subsets of type II airway epithelial cells (ATII) and alveolar macrophages (AM) recovered from the lungs of mock- or IAV-infected mice at 9 h postinfection were subjected to RNA sequencing. This time point was chosen to allow for characterization of cell types first infected with the virus inoculum, prior to multicycle virus replication and the infiltration of inflammatory cells into the airways. In the absence of infection, AM predominantly expressed genes related to immunity, whereas ATII expressed genes consistent with their physiological roles in the lung. Following IAV infection, AM almost exclusively activated cell-intrinsic antiviral pathways that were dependent on interferon (IFN) regulatory factor 3/7 (IRF3/7) and/or type I IFN signaling. In contrast, IAV-infected ATII activated a broader range of physiological responses, including cell-intrinsic antiviral pathways, which were both independent of and dependent on IRF3/7 and/or type I IFN. These data suggest that transcriptional profiles hardwired during development are a major determinant underlying the different responses of ATII and AM to IAV infection.IMPORTANCE Airway epithelial cells (AEC) and airway macrophages (AM) represent major targets of influenza A virus (IAV) infection in the lung, yet the two cell types respond very differently to IAV infection. We have used RNA sequencing to define the host transcriptional responses in each cell type under steady-state conditions as well as following IAV infection. To do this, different cell subsets isolated from the lungs of mock- and IAV-infected mice were subjected to RNA sequencing. Under steady-state conditions, AM and AEC express distinct transcriptional activities, consistent with distinct physiological roles in the airways. Not surprisingly, these cells also exhibited major differences in transcriptional responses following IAV infection. These studies shed light on how the different transcriptional architectures of airway cells from two different lineages drive transcriptional responses to IAV infection.

A comparative transcriptome analysis of the novel obligate methanotroph Methylomonas sp. DH-1 reveals key differences in transcriptional responses in C1 and secondary metabolite pathways during growth on methane and methanol

BackgroundMethanotrophs play an important role in biotechnological applications, with their ability to utilize single carbon (C1) feedstock such as methane and methanol to produce a range of high-value compounds. A newly isolated obligate methanotroph strain, Methylomonas sp. DH-1, became a platform strain for biotechnological applications because it has proven capable of producing chemicals, fuels, and secondary metabolites from methane and methanol. In this study, transcriptome analysis with RNA-seq was used to investigate the transcriptional change of Methylomonas sp. DH-1 on methane and methanol. This was done to improve knowledge about C1 assimilation and secondary metabolite pathways in this promising, but under-characterized, methane-bioconversion strain.ResultsWe integrated genomic and transcriptomic analysis of the newly isolated Methylomonas sp. DH-1 grown on methane and methanol. Detailed transcriptomic analysis indicated that (i) Methylomonas sp. DH-1 possesses the ribulose monophosphate (RuMP) cycle and the Embden–Meyerhof–Parnas (EMP) pathway, which can serve as main pathways for C1 assimilation, (ii) the existence and the expression of a complete serine cycle and a complete tricarboxylic acid (TCA) cycle might contribute to methane conversion and energy production, and (iii) the highly active endogenous plasmid pDH1 may code for essential metabolic processes. Comparative transcriptomic analysis on methane and methanol as a sole carbon source revealed different transcriptional responses of Methylomonas sp. DH-1, especially in C1 assimilation, secondary metabolite pathways, and oxidative stress. Especially, these results suggest a shift of central metabolism when substrate changed from methane to methanol in which formaldehyde oxidation pathway and serine cycle carried more flux to produce acetyl-coA and NADH. Meanwhile, downregulation of TCA cycle when grown on methanol may suggest a shift of its main function is to provide de novo biosynthesis, but not produce NADH.ConclusionsThis study provides insights into the transcriptomic profile of Methylomonas sp. DH-1 grown on major carbon sources for C1 assimilation, providing in-depth knowledge on the metabolic pathways of this strain. These observations and analyses can contribute to future metabolic engineering with the newly isolated, yet under-characterized, Methylomonas sp. DH-1 to enhance its biochemical application in relevant industries.

Induction of Cell Cycle and NK Cell Responses by Live-Attenuated Oral Vaccines against Typhoid Fever.

The mechanisms by which oral, live-attenuated vaccines protect against typhoid fever are poorly understood. Here, we analyze transcriptional responses after vaccination with Ty21a or vaccine candidate, M01ZH09. Alterations in response profiles were related to vaccine-induced immune responses and subsequent outcome after wild-type Salmonella Typhi challenge. Despite broad genetic similarity, we detected differences in transcriptional responses to each vaccine. Seven days after M01ZH09 vaccination, marked cell cycle activation was identified and associated with humoral immunogenicity. By contrast, vaccination with Ty21a was associated with NK cell activity and validated in peripheral blood mononuclear cell stimulation assays confirming superior induction of an NK cell response. Moreover, transcriptional signatures of amino acid metabolism in Ty21a recipients were associated with protection against infection, including increased incubation time and decreased severity. Our data provide detailed insight into molecular immune responses to typhoid vaccines, which could aid the rational design of improved oral, live-attenuated vaccines against enteric pathogens.

Colonic transcriptional response to 1α,25(OH)2 vitamin D3 in African- and European-Americans

Colorectal cancer (CRC) is a significant health burden especially among African Americans (AA). Epidemiological studies have correlated low serum vitamin D with CRC risk, and, while hypovitaminosis D is more common and more severe in AA, the mechanisms by which vitamin D modulates CRC risk and how these differ by race are not well understood. Active vitamin D (1α,25(OH)2D3) has chemoprotective effects primarily through transcriptional regulation of target genes in the colon. We hypothesized that transcriptional response to 1α,25(OH)2D3 differs between AA and European Americans (EA) irrespective of serum vitamin D and that regulatory variants could impact transcriptional response. We treated ex vivo colon cultures from 34 healthy subjects (16 AA and 18 EA) with 0.1μM 1α,25(OH)2D3 or vehicle control for 6h and performed genome-wide transcriptional profiling. We found 8 genes with significant differences in transcriptional response to 1α,25(OH)2D3 between AA and EA with definitive replication of inter-ethnic differences for uridine phosphorylase 1 (UPP1) and zinc finger-SWIM containing 4 (ZSWIM4). We performed expression quantitative trait loci (eQTL) mapping and identified response cis-eQTLs for ZSWIM4 as well as for histone deacetylase 3 (HDAC3), the latter of which showed a trend toward significant inter-ethnic differences in transcriptional response. Allele frequency differences of eQTLs for ZSWIM4 and HDAC3 accounted for observed transcriptional differences between populations. Taken together, our results demonstrate that transcriptional response to 1α,25(OH)2D3 differs between AA and EA independent of serum 25(OH)D levels. We provide evidence in support of a genetic regulatory mechanism underlying transcriptional differences between populations for ZSWIM4 and HDAC3. Further work is needed to elucidate how response eQTLs modify vitamin D response and whether genotype and/or transcriptional response correlate with chemopreventive effects. Relevant biomarkers, such as tissue-specific 1α,25(OH)2D3 transcriptional response, could identify individuals likely to benefit from vitamin D for CRC prevention as well as elucidate basic mechanisms underlying CRC disparities.

Discriminative identification of transcriptional responses of promoters and enhancers after stimulus.

Promoters and enhancers regulate the initiation of gene expression and maintenance of expression levels in spatial and temporal manner. Recent findings stemming from the Cap Analysis of Gene Expression (CAGE) demonstrate that promoters and enhancers, based on their expression profiles after stimulus, belong to different transcription response subclasses. One of the most promising biological features that might explain the difference in transcriptional response between subclasses is the local chromatin environment. We introduce a novel computational framework, PEDAL, for distinguishing effectively transcriptional profiles of promoters and enhancers using solely histone modification marks, chromatin accessibility and binding sites of transcription factors and co-activators. A case study on data from MCF-7 cell-line reveals that PEDAL can identify successfully the transcription response subclasses of promoters and enhancers from two different stimulations. Moreover, we report subsets of input markers that discriminate with minimized classification error MCF-7 promoter and enhancer transcription response subclasses. Our work provides a general computational approach for identifying effectively cell-specific and stimulation-specific promoter and enhancer transcriptional profiles, and thus, contributes to improve our understanding of transcriptional activation in human.

Patterns of Transcriptional Response to 1,25-Dihydroxyvitamin D3 and Bacterial Lipopolysaccharide in Primary Human Monocytes.

The active form of vitamin D, 1,25-dihydroxyvitamin D3 (1,25D), plays an important immunomodulatory role, regulating transcription of genes in the innate and adaptive immune system. The present study examines patterns of transcriptome-wide response to 1,25D, and the bacterial lipopolysaccharide (LPS) in primary human monocytes, to elucidate pathways underlying the effects of 1,25D on the immune system. Monocytes obtained from healthy individuals of African-American and European-American ancestry were treated with 1,25D, LPS, or both, simultaneously. The addition of 1,25D during stimulation with LPS induced significant upregulation of genes in the antimicrobial and autophagy pathways, and downregulation of proinflammatory response genes compared to LPS treatment alone. A joint Bayesian analysis enabled clustering of genes into patterns of shared transcriptional response across treatments. The biological pathways enriched within these expression patterns highlighted several mechanisms through which 1,25D could exert its immunomodulatory role. Pathways such as mTOR signaling, EIF2 signaling, IL-8 signaling, and Tec Kinase signaling were enriched among genes with opposite transcriptional responses to 1,25D and LPS, respectively, highlighting the important roles of these pathways in mediating the immunomodulatory activity of 1,25D. Furthermore, a subset of genes with evidence of interethnic differences in transcriptional response was also identified, suggesting that in addition to the well-established interethnic variation in circulating levels of vitamin D, the intensity of transcriptional response to 1,25D and LPS also varies between ethnic groups. We propose that dysregulation of the pathways identified in this study could contribute to immune-mediated disease risk.

Abstract B19: Colonic transcriptional response to 1,25(OH)2 vitamin D in African and European Americans

Abstract Background: 1,25(OH)2 vitamin D (1,25vitD) is a steroid hormone with anti-cancer effects that functions through direct transcriptional mechanisms. In colorectal cancer (CRC), 1,25vitD is thought to play a protective role; however, it is not known to what extent inter–ethnic differences in transcriptional response could underlie differences in disease susceptibility in African-Americans (AA) whose incidence of CRC is the highest in the US. To study differences in short-term transcriptional response in the colon, we have optimized an ex vivo culture system using colonic biopsies obtained during colonoscopy, thereby eliminating a number of confounders. We use this system to study transcriptional responses between ethnicities. Methods: Healthy subjects (n=69; 34 in discovery cohort &amp; 35 in replication cohort) undergoing colonoscopy were consented and colonic biopsies were obtained at the rectosigmoid junction. Two colonic biopsies each were treated with 0.1μM 1,25vitD or vehicle control (EtOH) and incubated for 6 hours at 37°C. For the discovery cohort, RNA from two treatment replicates was pooled and hybridized to Illumina Human HT-12 expression chips. Transcriptional response was analyzed using the packages LUMI and LIMMA in R adjusting for covariates and multiple testing. For the replication cohort, expression for select genes was performed using quantitative real-time PCR. In silico analysis was done to identify vitamin D receptor (VDR) binding peaks using ChIP-sequence data from LS180 cell lines. Results: In the discovery cohort, 16 AA and 18 European Americans (EA) were included. As proof of concept, the top two significantly up-regulated genes were CD14 (p=1.15x10-25) and CYP24A1 (p=7.65x10-23). There were 10 genes that showed significantly different responses to 1,25vitD between AA and EA including UPP1, UCKL1, ZSWIM4, ERCC1, MFSDA2, EPHA2, CLRN3, KIAA1324, TRIP8 and WDR55. In the replication cohort, 20 AA and 15 EA were included. Of 5 genes tested to date, 3 genes (UPP1, ERCC1 &amp; ZSWIM4) showed evidence of replication with similar direction and effect sizes as the discovery cohort. Two genes (UCKL1 &amp; UPP1) have VDR binding peaks within an intron (UCKL1) or 21kb downstream (UPP1). Within these peaks, there are SNPs with allele frequency differences between Africans and Europeans. In particular, the UCKL1 SNP rs58190754 that is located within the DR3 binding motif and is predicted to be functionally significant for VDR binding has an allele frequency of 10% in YRI &amp; AA and &amp;lt;1% in Europeans. Conclusions: Using short-term ex vivo culture of primary colonic tissue, we found significantly up- and down-regulated genes in response to 1,25vitD including a number of established 1,25vitD targets. There were 10 genes with significantly different transcriptional responses to 1,25vitD treatment between AA and EA of which 3 genes have been replicated in preliminary work. Two genes (UPP1 &amp; UCKL1) harbor VDR binding peaks with candidate SNPs that could be functionally relevant for gene expression differences between populations. These genes are involved in pyrimidine salvage and could have a role in cancer susceptibility. Future work will genotype SNPs and determine whether genotype predicts expression. In addition, functional assays of pyrimidine salvage are being tested. Finally, replication of additional genes is ongoing. Further studies of these genes may reveal important pathways underlying CRC disparities. Citation Format: Sonia Kupfer, Brandon Mapes, Shigeki Nakagome, David Witonsky, Anna Di Rienzo. Colonic transcriptional response to 1,25(OH)2 vitamin D in African and European Americans. [abstract]. In: Proceedings of the Eighth AACR Conference on The Science of Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; Nov 13-16, 2015; Atlanta, GA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2016;25(3 Suppl):Abstract nr B19.

Deeply diverged alleles in the Arabidopsis AREB1 transcription factor drive genome-wide differences in transcriptional response to the environment.

Gene regulatory variation is an important driver of the evolution of physiological and developmental responses to the environment. The abscisic acid (ABA) signaling pathway has long been studied as a key component of the cellular response to abiotic stresses in plants. We identify two haplotypes in an Arabidopsis thaliana transcription factor, AREB1, which plays a central role in ABA-mediated response to osmotic stress. These two haplotypes show the sequence signature of long-term maintenance of genetic diversity, suggesting a role for a diversifying selection process such as balancing selection. We find that the two haplotypes, distinguished by a large number of single nucleotide polymorphisms and the presence or absence of four small insertion/deletions in AREB1 intron 1 and exon 2, are at roughly equal frequencies in Arabidopsis, and show high linkage disequilibrium and deep sequence divergence. We use a transgenic approach, along with mRNA Sequencing-based assay of genome-wide expression levels, and find considerable functional divergence between alleles representing the two haplotype groups. Specifically, we find that, under benign soil-water conditions, transgenic lines containing different AREB1 alleles differ in the expression of a large number of genes associated with pathogen response. There are relatively modest gene expression differences between the two transgenic lines under restricted soil water content. Our finding of pathogen-related activity expands the known roles of AREB1 in A. thaliana and reveals the molecular basis of gene-by-environment interaction in a putatively adaptive plant regulatory protein.