Analysis of gene expression in the developing mouse retina

Single-cell transcriptional logic of cell-fate specification and axon guidance in early-born retinal neurons

Retinal ganglion cells (RGCs), cone photoreceptors (cones), horizontal cells and amacrine cells are the first classes of neurons produced in the retina. However, an important question is how this diversity of cell states is transcriptionally produced. Here, we profiled 6067 single retinal cells to provide a comprehensive transcriptomic atlas showing the diversity of the early developing mouse retina. RNA velocities unveiled the dynamics of cell cycle coordination of early retinogenesis and define the transcriptional sequences at work during the hierarchical production of early cell-fate specification. We show that RGC maturation follows six waves of gene expression, with older-generated RGCs transcribing increasing amounts of guidance cues for young peripheral RGC axons that express the matching receptors. Spatial transcriptionally deduced features in subpopulations of RGCs allowed us to define novel molecular markers that are spatially restricted. Finally, the isolation of such a spatially restricted population, ipsilateral RGCs, allowed us to identify their molecular identity at the time they execute axon guidance decisions. Together, these data represent a valuable resource shedding light on transcription factor sequences and guidance cue dynamics during mouse retinal development.

Expression of CRYM in different rat organs during development and its decreased expression in degenerating pyramidal tracts in amyotrophic lateral sclerosis

The protein μ-crystallin (CRYM) is a novel component of the marsupial lens that has two functions: it is a key regulator of thyroid hormone transportation and a reductase of sulfur-containing cyclic ketimines. In this study, we examined changes of the expression pattern of CRYM in different rat organs during development using immunohistochemistry and immunoblotting. As CRYM is reportedly expressed in the corticospinal tract, we also investigated CRYM expression in human cases of amyotrophic lateral sclerosis (ALS) using immunohistochemistry. In the rat brain, CRYM was expressed in the cerebral cortex, basal ganglia, hippocampus and corticospinal tract in the early postnatal period. As postnatal development progressed, CRYM expression was restricted to large pyramidal neurons in layers V and VI of the cerebral cortex and pyramidal cells in the deep layer of CA1 in the hippocampus. Even within the same regions, CRYM-positive and negative neurons were distributed in a mosaic pattern. In the kidney, CRYM was expressed in epithelial cells of the proximal tubule and mesenchymal cells of the medulla in the early postnatal period; however, CRYM expression in the medulla was lost as mesenchymal cell numbers decreased with the rapid growth of the medulla. In human ALS brains, we observed marked loss of CRYM in the corticospinal tract, especially distally. Our results suggest that CRYM may play roles in development of cortical and hippocampal pyramidal cells in the early postnatal period, and in the later period, performs cell-specific functions in selected neuronal populations. In the kidney, CRYM may play roles in maturation of renal function. The expression patterns of CRYM may reflect significance of its interactions with T3 or ketimines in these cells and organs. The results also indicate that CRYM may be used as a marker of axonal degeneration in the corticospinal tract.

Microarray and morphological analysis of early postnatal CRB2 mutant retinas on a pure C57BL/6J genetic background

In humans, the Crumbs homologue-1 (CRB1) gene is mutated in progressive types of autosomal recessive retinitis pigmentosa and Leber congenital amaurosis. The severity of the phenotype due to human CRB1 or mouse Crb1 mutations is dependent on the genetic background. Mice on C57BL/6J background with Crb1 mutations show late onset of retinal spotting phenotype or no phenotype. Recently, we showed that conditional deletion of mouse Crb2 in the retina results in early retinal disorganization leading to severe and progressive retinal degeneration with concomitant visual loss that mimics retinitis pigmentosa due to mutations in the CRB1 gene. Recent studies in the fruit fly and zebrafish suggest roles of the Crumbs (CRB) complex members in the regulation of cellular signalling pathways including the Notch1, mechanistic target of rapamycin complex 1 (mTORC1) and the Hippo pathway. Here, we demonstrate that mice backcrossed to C57BL/6J background with loss of CRB2 in the retina show a progressive disorganization and degeneration phenotype during late retinal development. We used microarray gene profiling to study the transcriptome of retinas lacking CRB2 during late retinal development. Unexpectedly, the retinas of newborn mice lacking CRB2 showed no changes in the transcriptome during retinal development. These findings suggest that loss of CRB2 in the developing retina results in retinal disorganization and subsequent degeneration without major changes in the transcriptome of the retina. These mice might be an interesting model to study the onset of retinal degeneration upon loss of CRB proteins.

Transcriptome analysis using next generation sequencing reveals molecular signatures of diabetic retinopathy and efficacy of candidate drugs

PURPOSE

To define gene expression changes associated with diabetic retinopathy in a mouse model using next generation sequencing, and to utilize transcriptome signatures to assess molecular pathways by which pharmacological agents inhibit diabetic retinopathy.

METHODS

We applied a high throughput RNA sequencing (RNA-seq) strategy using Illumina GAIIx to characterize the entire retinal transcriptome from nondiabetic and from streptozotocin-treated mice 32 weeks after induction of diabetes. Some of the diabetic mice were treated with inhibitors of receptor for advanced glycation endproducts (RAGE) and p38 mitogen activated protein (MAP) kinase, which have previously been shown to inhibit diabetic retinopathy in rodent models. The transcripts and alternatively spliced variants were determined in all experimental groups.

RESULTS

Next generation sequencing-based RNA-seq profiles provided comprehensive signatures of transcripts that are altered in early stages of diabetic retinopathy. These transcripts encoded proteins involved in distinct yet physiologically relevant disease-associated pathways such as inflammation, microvasculature formation, apoptosis, glucose metabolism, Wnt signaling, xenobiotic metabolism, and photoreceptor biology. Significant upregulation of crystallin transcripts was observed in diabetic animals, and the diabetes-induced upregulation of these transcripts was inhibited in diabetic animals treated with inhibitors of either RAGE or p38 MAP kinase. These two therapies also showed dissimilar regulation of some subsets of transcripts that included alternatively spliced versions of arrestin, neutral sphingomyelinase activation associated factor (Nsmaf), SH3-domain GRB2-like interacting protein 1 (Sgip1), and axin.

CONCLUSIONS

Diabetes alters many transcripts in the retina, and two therapies that inhibit the vascular pathology similarly inhibit a portion of these changes, pointing to possible molecular mechanisms for their beneficial effects. These therapies also changed the abundance of various alternatively spliced versions of signaling transcripts, suggesting a possible role of alternative splicing in disease etiology. Our studies clearly demonstrate RNA-seq as a comprehensive strategy for identifying disease-specific transcripts, and for determining comparative profiles of molecular changes mediated by candidate drugs.

Step-wise specification of retinal stem cells during normal embryogenesis

The specification of embryonic cells to produce the retina begins at early embryonic stages as a multi-step process that gradually restricts fate potentials. First, a subset of embryonic cells becomes competent to form retina by their lack of expression of endo-mesoderm-specifying genes. From these cells, a more restricted subset is biased to form retina by virtue of their close proximity to sources of bone morphogenetic protein antagonists during neural induction. During gastrulation, the definitive RSCs (retinal stem cells) are specified as the eye field by interactions with underlying mesoderm and the expression of a network of retina-specifying genes. As the eye field is transformed into the optic vesicle and optic cup, a heterogeneous population of RPCs (retinal progenitor cells) forms to give rise to the different domains of the retina: the optic stalk, retinal pigmented epithelium and neural retina. Further diversity of RPCs appears to occur under the influences of cell-cell interactions, cytokines and combinations of regulatory genes, leading to the differentiation of a multitude of different retinal cell types. This review examines what is known about each sequential step in retinal specification during normal vertebrate development, and how that knowledge will be important to understand how RSCs might be manipulated for regenerative therapies to treat retinal diseases.

Genetic variations strongly influence phenotypic outcome in the mouse retina

Variation in genetic background can significantly influence the phenotypic outcome of both disease and non-disease associated traits. Additionally, differences in temporal and strain specific gene expression can also contribute to phenotypes in the mammalian retina. This is the first report of microarray based cross-strain analysis of gene expression in the retina investigating genetic background effects. Microarray analyses were performed on retinas from the following mouse strains: C57BL6/J, AKR/J, CAST/EiJ, and NOD.NON-H2^-nb1 at embryonic day 18.5 (E18.5) and postnatal day 30.5 (P30.5). Over 3000 differentially expressed genes were identified between strains and developmental stages. Differential gene expression was confirmed by qRT-PCR, Western blot, and immunohistochemistry. Three major gene networks were identified that function to regulate retinal or photoreceptor development, visual perception, cellular transport, and signal transduction. Many of the genes in these networks are implicated in retinal diseases such as bradyopsia, night-blindness, and cone-rod dystrophy. Our analysis revealed strain specific variations in cone photoreceptor cell patterning and retinal function. This study highlights the substantial impact of genetic background on both development and function of the retina and the level of gene expression differences tolerated for normal retinal function. These strain specific genetic variations may also be present in other tissues. In addition, this study will provide valuable insight for the development of more accurate models for human retinal diseases.

The genetics of pain: implications for evaluation and treatment of spinal disease

BACKGROUND CONTEXT

Variability in human pain experience appears to be at least partially determined by genetic inheritance. To the extent that awareness of individual pain sensitivity and the tendency to develop chronic pain after injury or surgery would be informative for clinical decision making, development and use of genetic testing for specific pain markers could contribute to improved outcomes in management of spinal disease.

PURPOSE

To review important and illustrative results from both classical and modern pain genetics studies and to introduce readers to critical definitions and concepts necessary to interpret the growing body of genetics literature relevant to spinal disease.

STUDY DESIGN/SETTING

Literature review and commentary.

METHODS

A review was performed of published English language studies in which genetic techniques were used to analyze the molecular basis of nociceptive signaling or processing with a particular emphasis on studies addressing genetic determinants of interindividual variability in pain sensitivity or predisposition to chronic pain.

RESULTS

There is compelling evidence indicating that interindividual differences in pain sensitivity and the risk of developing chronic pain syndromes are genetically determined. Despite a growing list of putative "pain genes," genetic association studies remain plagued with difficulty replicating initial findings in different cohorts.

CONCLUSIONS

Genome-wide association studies are potentially powerful means of identifying clinically relevant genetic markers predicting disease susceptibility, severity, and treatment response. However, accurate results require rigorous study design with use of large homogeneous populations and precise phenotypes.

Transcriptional profile analysis of RPGRORF15 frameshift mutation identifies novel genes associated with retinal degeneration

PURPOSE

To identify genes and molecular mechanisms associated with photoreceptor degeneration in a canine model of XLRP caused by an RPGR exon ORF15 microdeletion. Methods. Expression profiles of mutant and normal retinas were compared by using canine retinal custom cDNA microarrays. qRT-PCR, Western blot analysis, and immunohistochemistry (IHC) were applied to selected genes, to confirm and expand the microarray results.

RESULTS

At 7 and 16 weeks, respectively, 56 and 18 transcripts were downregulated in the mutant retinas, but none were differentially expressed (DE) at both ages, suggesting the involvement of temporally distinct pathways. Downregulated genes included the known retina-relevant genes PAX6, CHML, and RDH11 at 7 weeks and CRX and SAG at 16 weeks. Genes directly or indirectly active in apoptotic processes were altered at 7 weeks (CAMK2G, NTRK2, PRKCB, RALA, RBBP6, RNF41, SMYD3, SPP1, and TUBB2C) and 16 weeks (SLC25A5 and NKAP). Furthermore, the DE genes at 7 weeks (ELOVL6, GLOD4, NDUFS4, and REEP1) and 16 weeks (SLC25A5 and TARS2) are related to mitochondrial functions. qRT-PCR of 18 genes confirmed the microarray results and showed DE of additional genes not on the array. Only GFAP was DE at 3 weeks of age. Western blot and IHC analyses also confirmed the high reliability of the transcriptomic data.

CONCLUSIONS

Several DE genes were identified in mutant retinas. At 7 weeks, a combination of nonclassic anti- and proapoptosis genes appear to be involved in photoreceptor degeneration, whereas at both 7 and 16 weeks, the expression of mitochondria-related genes indicates that they may play a relevant role in the disease process.

Temporal expression of angiogenesis-related genes in developing neonatal rodent retina: a novel in vivo model to study cerebral vascular development

BACKGROUND

Experimental models to study cerebrovascular malformations are limited therefore we used the neonatal rodent retina as a model to study cerebral angiogenesis.

OBJECTIVE

We performed a gene expression analysis to define temporal changes in the expression of 96 angiogenesis-related genes during retinal vascularization.

METHODS

A total of 72 retinas from 36 newborn C57BL/6 mice were used. Sets of neonatal mouse retinas were surgically isolated by 2-day intervals starting from postnatal day 0 to day 20 and at the 32nd day (representing adult retinas). For each of these 12 time points in the postnatal developmental period of mouse retinas, separate sets of 6 retinas from 3 mice were pooled, and their RNA was hybridized to an angiogenesis-specific gene array. Temporal expression patterns of each of the 96 angiogenesis-related genes were analyzed. For confirmation, vascular endothelial growth factor protein expression was also studied by immunohistochemistry.

RESULTS

Twenty-two of the 96 genes analyzed displayed a significantly different temporal expression profile, and the rest exhibited a static expression, as compared to the human glyceraldehyde-3-phosphate dehydrogenase gene. Among these genes, the temporal pattern of expression was variable, but peaks were seen mostly on days 8, 10, 12, and 16. This timing corresponds well to morphologic changes that occur in the retina during different stages of angiogenesis.

CONCLUSION

The neonatal rodent retina, which has a cellular architecture similar to that of the brain, has active and quantifiable angiogenic activity during the neonatal period and can be used as a simple and convenient model to study cerebral angiogenesis.

The transcriptome of Populus in elevated CO reveals increased anthocyanin biosynthesis during delayed autumnal senescence

*The delay in autumnal senescence that has occurred in recent decades has been linked to rising temperatures. Here, we suggest that increasing atmospheric CO2 may partly account for delayed autumnal senescence and for the first time, through transcriptome analysis, identify gene expression changes associated with this delay. *Using a plantation of Populus x euramericana grown in elevated [CO2] (e[CO2]) with free-air CO2 enrichment (FACE) technology, we investigated the molecular and biochemical basis of this response. A Populus cDNA microarray was used to identify genes representing multiple biochemical pathways influenced by e[CO2] during senescence. Gene expression changes were confirmed through real-time quantitative PCR, and leaf biochemical assays. *Pathways for secondary metabolism and glycolysis were significantly up-regulated by e[CO2] during senescence, in particular, those related to anthocyanin biosynthesis. Expressed sequence tags (ESTs) representing the two most significantly up-regulated transcripts in e[CO2], LDOX (leucoanthocyanidin dioxgenase) and DFR (dihydroflavonol reductase), gave (e[CO2]/ambient CO(2) (a[CO2])) expression ratios of 39.6 and 19.3, respectively. *We showed that in e[CO2] there was increased autumnal leaf sugar accumulation and up-regulation of genes determining anthocyanin biosynthesis which, we propose, prolongs leaf longevity during natural autumnal senescence.

Computational prediction of neural progenitor cell fates

Understanding how stem and progenitor cells choose between alternative cell fates is a major challenge in developmental biology. Efforts to tackle this problem have been hampered by the scarcity of markers that can be used to predict cell division outcomes. Here we present a computational method, based on algorithmic information theory, to analyze dynamic features of living cells over time. Using this method, we asked whether rat retinal progenitor cells (RPCs) display characteristic phenotypes before undergoing mitosis that could foretell their fate. We predicted whether RPCs will undergo a self-renewing or terminal division with 99% accuracy, or whether they will produce two photoreceptors or another combination of offspring with 87% accuracy. Our implementation can segment, track and generate predictions for 40 cells simultaneously on a standard computer at 5 min per frame. This method could be used to isolate cell populations with specific developmental potential, enabling previously impossible investigations.

Connexin-mediated communication controls cell proliferation and is essential in retinal histogenesis

Connexin (Cx) channels and hemichannels are involved in essential processes during nervous system development such as apoptosis, propagation of spontaneous activity and interkinetic nuclear movement. In the first part of this study, we extensively characterized Cx gene and protein expression during retinal histogenesis. We observed distinct spatio-temporal patterns among studied Cx and an overriding, ubiquitous presence of Cx45 in progenitor cells. The role of Cx-mediated communication was assessed by using broad-spectrum (carbenoxolone, CBX) and Cx36/Cx50 channel-specific (quinine) blockers. In vivo application of CBX, but not quinine, caused remarkable reduction in retinal thickness, suggesting changes in cell proliferation/apoptosis ratio. Indeed, we observed a decreased number of mitotic cells in CBX-injected retinas, with no significant changes in the expression of PCNA, a marker for cells in proliferative state. Taken together, our results pointed a pivotal role of Cx45 in the developing retina. Moreover, this study revealed that Cx-mediated communication is essential in retinal histogenesis, particularly in the control of cell proliferation.

A unique program for cell death in xylem fibers of Populus stem

Maturation of the xylem elements involves extensive deposition of secondary cell-wall material and autolytic processes resulting in cell death. We describe here a unique type of cell-death program in xylem fibers of hybrid aspen (Populus tremula x P. tremuloides) stems, including gradual degradative processes in both the nucleus and cytoplasm concurrently with the phase of active cell-wall deposition. Nuclear DNA integrity, as determined by TUNEL (terminal deoxynucleotidyl transferase-mediated dUTP nick end labeling) and Comet (single-cell gel electrophoresis) assays, was compromised early during fiber maturation. In addition, degradation of the cytoplasmic contents, as detected by electron microscopy of samples fixed by high-pressure freezing/freeze substitution (HPF-FS), was gradual and resulted in complete loss of the cytoplasmic contents well before the loss of vacuolar integrity, which is considered to be the moment of death. This type of cell death differs significantly from that seen in xylem vessels. The loss of vacuolar integrity, which is thought to initiate cell degradative processes in the xylem vessels, is one of the last processes to occur before the final autolysis of the remaining cell contents in xylem fibers. High-resolution microarray analysis in the vascular tissues of Populus stem, combined with in silico analysis of publicly available data repositories, suggests the involvement of several previously uncharacterized transcription factors, ethylene, sphingolipids and light signaling as well as autophagy in the control of fiber cell death.

Gene expression patterns in hypoxic and post-hypoxic adult rat retina with special reference to the NMDA receptor and its interactome

PURPOSE

A gene expression analysis of hypoxic rat retina was undertaken to gain a deeper understanding of the possible molecular mechanisms that underlie hypoxia-induced retinal pathologies and identify possible therapeutic targets.

METHODS

Rats were made severely hypoxic (6%-7% O(2)) for 3 h. Some rats were sacrificed at this time, and others were allowed to recover for 24 h under normoxic conditions. A focused oligonucleotide microarray of 1,178 genes, qRT-PCR of selected transcripts, and western analysis of hypoxia inducible factor-1alpha (HIF-1alpha) were used to compare retinas from the hypoxic and recovery groups to control animals that were not made hypoxic. SAM analysis was used to identify statistically significant changes in microarray data, and the bioinformatics programs GoMiner, Gene Set Enrichment Analysis (GSEA), and HiMAP were used to identify significant ontological categories and analyze the N-methyl-D-aspartate (NMDA) receptor interactome.

RESULTS

HIF-1alpha protein, but not mRNA, was elevated up to 15-fold during hypoxia, beginning at 0.5 h, the shortest duration examined. Of the total of 1,178 genes examined by microarray, 119 were significantly upregulated following hypoxia. Of these, 86 were still significantly upregulated following recovery. However, 24 genes were significantly downregulated following hypoxia, with 12 still significantly downregulated after recovery. Of the 1035 genes that did not change with hypoxia, the expression of 36 genes was significantly changed after recovery. Ontological analyses showed significant upregulation of a large number of genes in the glutamate receptor family, including 3 of the 5 NMDA subunits. qRT-PCR analysis further corroborated these findings. Genes known to directly interact specifically with the NR1 subunit of the NMDA receptor were identified using HiMAP databases. GSEA analysis revealed that these genes were not affected by either hypoxia or altered after recovery.

CONCLUSIONS

The identification of gene expression alterations as a function of hypoxia and recovery from hypoxia is important to understand the molecular mechanisms underlying retinal dysfunction associated with a variety of diseases. Gene changes were identified in hypoxic retina that could be linked to specific networks. Retinas recovering from hypoxia also showed distinct patterns of gene expression that were different from both normoxic control retinas and hypoxic retinas, indicating that hypoxia initiates a complex pattern of gene expression. Diseases of which hypoxia is a component may exhibit the several changes found here. Several potential therapeutic targets have been identified by our approach, including modulation of NMDA receptor expression and signaling, which until now have only been shown to play a role in responding to ischemia.

Global expression profiling in leaves of free-growing aspen

BACKGROUND

Genomic studies are routinely performed on young plants in controlled environments which is very different from natural conditions. In reality plants in temperate countries are exposed to large fluctuations in environmental conditions, in the case of perennials over several years. We have studied gene expression in leaves of a free-growing aspen (Populus tremula) throughout multiple growing seasons

RESULTS

We show that gene expression during the first month of leaf development was largely determined by a developmental program although leaf expansion, chlorophyll accumulation and the speed of progression through this program was regulated by the temperature. We were also able to define "transcriptional signatures" for four different substages of leaf development. In mature leaves, weather factors were important for gene regulation.

CONCLUSION

This study shows that multivariate methods together with high throughput transcriptional methods in the field can provide additional, novel information as to plant status under changing environmental conditions that is impossible to mimic in laboratory conditions. We have generated a dataset that could be used to e.g. identify marker genes for certain developmental stages or treatments, as well as to assess natural variation in gene expression.

Gene expression profiling of zebrafish embryonic retina

A methodology for microdissecting intact retinas from zebrafish embryos at early developmental stages for expression profiling was developed in this study. Total RNA was extracted consistently and reproducibly from the dissected retinas using a customized extraction protocol. The results from microarray experiments indicated that the purified RNA samples faithfully represented the biological differences among different types of samples. Genes that were differentially expressed in a particular neuronal layer or region of the retina were detectable by microarray experiments. In conclusion, this methodology makes it possible to obtain retinal-specific total RNA for genomics research on retinal development in zebrafish.

Microarray-based gene expression analysis during retinal maturation of albino rats

BACKGROUND

In recent years, the rat has become a commonly-used animal model for the study of retinal diseases. Similar to other tissues, the retina undergoes significant functional changes during maturation. Aiming to gain knowledge on additional aspects of retinal maturation, we performed gene expression and histological analyses of the rat retina during maturation.

METHODS

Rat retinas were dissected at three time points. Histological examination of the samples was performed, and the expression levels of retinal genes were evaluated using the rat whole-genome microarray system. Quantitative real-time PCR analysis was used to validate selected expression patterns. Various statistical and bioinformatic tools were used to identify differentially expressed genes.

RESULTS

The microarray analysis revealed a relatively high number of highly expressed non-annotated genes. We identified 603 differentially expressed genes, which were grouped into six clusters based on changes in expression levels during the first 20 weeks of life. A bioinformatic analysis of these clusters revealed sets of genes encoding proteins with functions that are likely to be relevant to retinal maturation (potassium, sodium, calcium, and chloride channels, synaptic vesicle transport, and axonogenesis). The histological analysis revealed a significant reduction of outer nuclear layer thickness and retinal ganglion cell number during maturation.

CONCLUSIONS

These data, taken together with our previously reported electrophysiological data, contribute to our understanding of the retinal maturation processes of this widely-used animal model.

Using a seed-network to query multiple large-scale gene expression datasets from the developing retina in order to identify and prioritize experimental targets

Understanding the gene networks that orchestrate the differentiation of retinal progenitors into photoreceptors in the developing retina is important not only due to its therapeutic applications in treating retinal degeneration but also because the developing retina provides an excellent model for studying CNS development. Although several studies have profiled changes in gene expression during normal retinal development, these studies offer at best only a starting point for functional studies focused on a smaller subset of genes. The large number of genes profiled at comparatively few time points makes it extremely difficult to reliably infer gene networks from a gene expression dataset. We describe a novel approach to identify and prioritize from multiple gene expression datasets, a small subset of the genes that are likely to be good candidates for further experimental investigation. We report progress on addressing this problem using a novel approach to querying multiple large-scale expression datasets using a 'seed network' consisting of a small set of genes that are implicated by published studies in rod photoreceptor differentiation. We use the seed network to identify and sort a list of genes whose expression levels are highly correlated with those of multiple seed network genes in at least two of the five gene expression datasets. The fact that several of the genes in this list have been demonstrated, through experimental studies reported in the literature, to be important in rod photoreceptor function provides support for the utility of this approach in prioritizing experimental targets for further experimental investigation. Based on Gene Ontology and KEGG pathway annotations for the list of genes obtained in the context of other information available in the literature, we identified seven genes or groups of genes for possible inclusion in the gene network involved in differentiation of retinal progenitor cells into rod photoreceptors. Our approach to querying multiple gene expression datasets using a seed network constructed from known interactions between specific genes of interest provides a promising strategy for focusing hypothesis-driven experiments using large-scale 'omics' data.

Analysis of gene expression during neurite outgrowth and regeneration

Background

The ability of a neuron to regenerate functional connections after injury is influenced by both its intrinsic state and also by extrinsic cues in its surroundings. Investigations of the transcriptional changes undergone by neurons during in vivo models of injury and regeneration have revealed many transcripts associated with these processes. Because of the complex milieu of interactions in vivo, these results include not only expression changes directly related to regenerative outgrowth and but also unrelated responses to surrounding cells and signals. In vitro models of neurite outgrowth provide a means to study the intrinsic transcriptional patterns of neurite outgrowth in the absence of extensive extrinsic cues from nearby cells and tissues.

Results

We have undertaken a genome-wide study of transcriptional activity in embryonic superior cervical ganglia (SCG) and dorsal root ganglia (DRG) during a time course of neurite outgrowth in vitro. Gene expression observed in these models likely includes both developmental gene expression patterns and regenerative responses to axotomy, which occurs as the result of tissue dissection. Comparison across both models revealed many genes with similar gene expression patterns during neurite outgrowth. These patterns were minimally affected by exposure to the potent inhibitory cue Semaphorin3A, indicating that this extrinsic cue does not exert major effects at the level of nuclear transcription. We also compared our data to several published studies of DRG and SCG gene expression in animal models of regeneration, and found the expression of a large number of genes in common between neurite outgrowth in vitro and regeneration in vivo.

Conclusion

Many gene expression changes undergone by SCG and DRG during in vitro outgrowth are shared between these two tissue types and in common with in vivo regeneration models. This suggests that the genes identified in this in vitro study may represent new candidates worthy of further study for potential roles in the therapeutic regrowth of neuronal connections.

Complement induction in spinal cord microglia results in anaphylatoxin C5a-mediated pain hypersensitivity

Microarray expression profiles reveal substantial changes in gene expression in the ipsilateral dorsal horn of the spinal cord in response to three peripheral nerve injury models of neuropathic pain. However, only 54 of the 612 regulated genes are commonly expressed across all the neuropathic pain models. Many of the commonly regulated transcripts are immune related and include the complement components C1q, C3, and C4, which we find are expressed only by microglia. C1q and C4 are, moreover, the most strongly regulated of all 612 regulated genes. In addition, we find that the terminal complement component C5 and the C5a receptor (C5aR) are upregulated in spinal microglia after peripheral nerve injury. Mice null for C5 had reduced neuropathic pain sensitivity, excluding C3a as a pain effector. C6-deficient rats, which cannot form the membrane attack complex, have a normal neuropathic pain phenotype. However, C5a applied intrathecally produces a dose-dependent, slow-onset cold pain in naive animals. Furthermore, a C5aR peptide antagonist reduces cold allodynia in neuropathic pain models. We conclude that induction of the complement cascade in spinal cord microglia after peripheral nerve injury contributes to neuropathic pain through the release and action of the C5a anaphylatoxin peptide.

Environmental and hormonal regulation of the activity-dormancy cycle in the cambial meristem involves stage-specific modulation of transcriptional and metabolic networks

We have performed transcript and metabolite profiling of isolated cambial meristem cells of the model tree aspen during the course of their activity-dormancy cycle to better understand the environmental and hormonal regulation of this process in perennial plants. Considerable modulation of cambial transcriptome and metabolome occurs throughout the activity-dormancy cycle. However, in addition to transcription, post-transcriptional control is also an important regulatory mechanism as exemplified by the regulation of cell-cycle genes during the reactivation of cambial cell division in the spring. Genes related to cold hardiness display temporally distinct induction patterns in the autumn which could explain the step-wise development of cold hardiness. Factors other than low temperature regulate the induction of early cold hardiness-related genes whereas abscisic acid (ABA) could potentially regulate the induction of late cold hardiness-related genes in the autumn. Starch breakdown in the autumn appears to be regulated by the 'short day' signal and plays a key role in providing substrates for the production of energy, fatty acids and cryoprotectants. Catabolism of sucrose and fats provides energy during the early stages of reactivation in the spring, whereas the reducing equivalents are generated through activation of the pentose phosphate shunt. Modulation of gibberellin (GA) signaling and biosynthesis could play a key role in the regulation of cambial activity during the activity-dormancy cycle as suggested by the induction of PttRGA which encodes a negative regulator of growth in the autumn and that of a GA-20 oxidase, a key gibberellin biosynthesis gene during reactivation in spring. In summary, our data reveal the dynamics of transcriptional and metabolic networks and identify potential targets of environmental and hormonal signals in the regulation of the activity-dormancy cycle in cambial meristem.

Candidate regulators of mammary branching morphogenesis identified by genome-wide transcript analysis

The mammary gland develops in a process known as branching morphogenesis, whereby a distal epithelial bud extends and bifurcates to form an extensive ductal network. Compared with other branched organs, such as the lung and kidney, little is known about the molecular basis of branching in the mammary gland. Here we report a microarray profiling strategy to identify novel genes that may regulate mammary branching. We microdissected terminal end bud (TEB) and mature duct microenvironments from beta-actin-green fluorescent protein reporter mice and compared their RNA expression profiles with epithelium-free mammary stroma by means of microarray. We identified 1,074 genes enriched in the TEB microenvironment, 222 genes enriched in the mature duct microenvironment, and 385 genes enriched in both TEB and mature duct microenvironments. The microarray correctly predicted the expression of genes known to be enriched in the epithelium (Ets-5) and stroma (MMP-14) of TEBs and in the mature duct microenvironment (MMP-3). The microarray also correctly predicted the localization of previously uncharacterized genes, such as the TEB-enriched SPRR-1a, the duct-enriched casein-gamma, and the general epithelial marker pleiotrophin. Analysis of genes enriched in TEBs revealed several genes in the Wnt (Wnt-2, Wnt-5a, Wnt-7b, Dsh-3, Frizzled-1, Frizzled-2), hedgehog (Dhh), ephrin (Ephrin-B1, Eph-A2), and transcription factor (Twist-1, Twist-2, Snail) families. In situ hybridization verified that these genes were enriched in the TEB epithelium (Wnt-5a, Wnt-7b, Dhh, Eph-A2) or TEB stroma (Wnt-2, Frizzled-1, Ephrin-B1). We discuss the potential roles of these genes in mammary branching morphogenesis.

Gene expression profiling of zebrafish embryonic retinal pigment epithelium in vivo

PURPOSE

Eye development and photoreceptor maintenance is dependent on the retinal pigment epithelium (RPE), a thin layer of cells that underlies the neural retina. Despite its importance, development of RPE has not been studied by a genomic approach. In this study, a microarray expression-profiling methodology was established for studying RPE development.

METHODS

The intact retina with RPE attached was dissected from developing embryos, and differentially expressed genes in RPE were inferred by comparing the dissected tissues with retinas without RPE, in microarray and statistical analyses.

RESULTS

Of the probesets used, 8810 were significantly expressed in RPE at 52 hours postfertilization (hpf), of which 1443 may have biologically meaningful expression levels. Further, 78 and 988 probesets were found to be significantly over- or underexpressed in RPE, respectively, compared with retina. Also, 79.2% (38/48) of the known overexpressed probesets were independently validated as RPE-related transcripts.

CONCLUSIONS

The results strongly suggest that this methodology can obtain in vivo RPE-specific gene expression from the zebrafish embryos and identify novel RPE markers.

DNA repair in mammalian embryos

Mammalian cells have developed complex mechanisms to identify DNA damage and activate the required response to maintain genome integrity. Those mechanisms include DNA damage detection, DNA repair, cell cycle arrest and apoptosis which operate together to protect the conceptus from DNA damage originating either in parental gametes or in the embryo's somatic cells. DNA repair in the newly fertilized preimplantation embryo is believed to rely entirely on the oocyte's machinery (mRNAs and proteins deposited and stored prior to ovulation). DNA repair genes have been shown to be expressed in the early stages of mammalian development. The survival of the embryo necessitates that the oocyte be sufficiently equipped with maternal stored products and that embryonic gene expression commences at the correct time. A Medline based literature search was performed using the keywords 'DNA repair' and 'embryo development' or 'gametogenesis' (publication dates between 1995 and 2006). Mammalian studies which investigated gene expression were selected. Further articles were acquired from the citations in the articles obtained from the preliminary Medline search. This paper reviews mammalian DNA repair from gametogenesis to preimplantation embryos to late gestational stages.

Retroviral misexpression of cVax disturbs retinal ganglion cell axon fasciculation and intraretinal pathfinding in vivo and guidance of nasal ganglion cell axons in vivo

The transcription factor cVax (Vax2) is expressed in the ventral neural retina and restricted expression is a prerequisite for at least three prominent aspects of retinal dorsal-ventral patterning: polarized expression of EphB/B-ephrin molecules, the retinotectal projection and the distribution of rod photoreceptors across the retina. In the chick retina, the fasciculation pattern of ganglion cell axons also differs between the dorsal and ventral eye. To investigate the molecular mechanisms involved, the nerve fiber layer was analyzed after retroviral misexpression of several factors known to regulate the positional specification of retinal ganglion cells. Forced cVax expression ventralized the fasciculation pattern and caused axon pathfinding errors near the optic disc. Ectopic expression of different ephrin molecules indicated that axon fasciculation is, at least in part, mediated by the EphB system. Finally, we report that retroviral misexpression of cVax increased the pool of EphA4 receptors phosphorylated on tyrosine residues and altered the guidance preference of nasal axons in vitro. These results identify novel functions for cVax in intraretinal axon fasciculation and pathfinding as well as suggest a mechanism to explain how restricted cVax expression may influence map formation along the dorso-ventral and antero-posterior axes of the optic tectum.

A biphasic pattern of gene expression during mouse retina development

BACKGROUND

Between embryonic day 12 and postnatal day 21, six major neuronal and one glia cell type are generated from multipotential progenitors in a characteristic sequence during mouse retina development. We investigated expression patterns of retina transcripts during the major embryonic and postnatal developmental stages to provide a systematic view of normal mouse retina development,

RESULTS

A tissue-specific cDNA microarray was generated using a set of sequence non-redundant EST clones collected from mouse retina. Eleven stages of mouse retina, from embryonic day 12.5 (El2.5) to postnatal day 21 (PN21), were collected for RNA isolation. Non-amplified RNAs were labeled for microarray experiments and three sets of data were analyzed for significance, hierarchical relationships, and functional clustering. Six individual gene expression clusters were identified based on expression patterns of transcripts through retina development. Two developmental phases were clearly divided with postnatal day 5 (PN5) as a separate cluster. Among 4,180 transcripts that changed significantly during development, approximately 2/3 of the genes were expressed at high levels up until PN5 and then declined whereas the other 1/3 of the genes increased expression from PN5 and remained at the higher levels until at least PN21. Less than 1% of the genes observed showed a peak of expression between the two phases. Among the later increased population, only about 40% genes are correlated with rod photoreceptors, indicating that multiple cell types contributed to gene expression in this phase. Within the same functional classes, however, different gene populations were expressed in distinct developmental phases. A correlation coefficient analysis of gene expression during retina development between previous SAGE studies and this study was also carried out.

CONCLUSION

This study provides a complementary genome-wide view of common gene dynamics and a broad molecular classification of mouse retina development. Different genes in the same functional clusters are expressed in the different developmental stages, suggesting that cells might change gene expression profiles from differentiation to maturation stages. We propose that large-scale changes in gene regulation during development are necessary for the final maturation and function of the retina.

Nocturnal changes in leaf growth of Populus deltoides are controlled by cytoplasmic growth

Growing leaves do not expand at a constant rate but exhibit pronounced diel growth rhythms. However, the mechanisms giving rise to distinct diel growth dynamics in different species are still largely unknown. As a first step towards identifying genes controlling rate and timing of leaf growth, we analysed the transcriptomes of rapidly expanding and fully expanded leaves of Populus deltoides Bartr. ex. Marsh at points of high and low expansion at night. Tissues with well defined temporal growth rates were harvested using an online growth-monitoring system based on a digital image sequence processing method developed for quantitative mapping of dicot leaf growth. Unlike plants studied previously, leaf growth in P. deltoides was characterised by lack of a base-tip gradient across the lamina, and by maximal and minimal growth at dusk and dawn, respectively. Microarray analysis revealed that the nocturnal decline in growth coincided with a concerted down-regulation of ribosomal protein genes, indicating deceleration of cytoplasmic growth. In a subsequent time-course experiment, Northern blotting and real-time RT-PCR confirmed that the ribosomal protein gene RPL12 and a cell-cycle gene H2B were down-regulated after midnight following a decrease in cellular carbohydrate concentrations. Thus, we propose that the spatio-temporal growth pattern in leaves of P. deltoides primarily arises from cytoplasmic growth whose activity increases from afternoon to midnight and thereafter decreases in this species.

Molecular targets of elevated [CO2] in leaves and stems of Populus deltoides: implications for future tree growth and carbon sequestration

We report the first comprehensive analysis of the effects of elevated [CO₂] on gene expression in source leaf and stem sink tissues in woody plants. We have taken advantage of coppiced Populus deltoides (Bartr.) stands grown for 3 years under three different and constant elevated [CO₂] in the agriforest mesocosms of Biosphere 2. Leaf area per tree was doubled by elevated [CO₂] but although growth at 800 v. 400 µmol mol^-1 CO₂ resulted in a significant increase in stem biomass, growth was not stimulated at 1200 µmol mol^-1 CO₂. Growth under elevated [CO₂] also resulted in significant increases in stem wood density. Analysis of expression data for the 13 490 clones present on POP1 microarrays revealed 95 and 277 [CO₂]-responsive clones in leaves and stems respectively, with the response being stronger at 1200 µmol mol^-1. When these [CO₂]-responsive genes were assigned to functional categories, metabolism-related genes were the most responsive to elevated [CO₂]. However within this category, expression of genes relating to bioenergetic processes was unchanged in leaves whereas the expression of genes for storage proteins and of those involved in control of wall expansion was enhanced. In contrast to leaves, the genes up-regulated in stems under elevated [CO₂] were primarily enzymes responsible for lignin formation and polymerisation or ethylene response factors, also known to induce lignin biosynthesis. Concomitant with this enhancement of lignin biosynthesis in stems, there was a pronounced repression of genes related to cell wall formation and cell growth. These changes in gene expression have clear consequences for long-term carbon sequestration, reducing the carbon-fertilisation effect, and the potential for increased lignification may negatively impact on future wood quality for timber and paper production.

Up-regulation and increased phosphorylation of protein kinase C (PKC) delta, mu and theta in the degenerating rd1 mouse retina

The rd1 mouse serves as a model for inherited photoreceptor degeneration: retinitis pigmentosa. Microarray techniques were employed to compare the transcriptomes of rd1 and congenic wild-type retinas at postnatal day 11, when degenerative processes have started but most photoreceptors are still present. Of the several genes that were differentially expressed, focus was put on those associated with the protein kinase C (PKC) signaling pathway, in particular PKCdelta, mu and theta. Microarray identified these as being up-regulated in the rd1 retina, which was confirmed by QRT-PCR. Western blotting and immunostaining, using antibodies against either total or phosphorylated variants of the PKC isoforms, revealed increased expression and phosphorylation of PKCdelta, mu and theta in the rd1 retina at the protein level as well. Our results suggest that these PKC isoforms are involved in rd1 degeneration.

Molecular gradients and development of retinotopic maps

Gradients of axon guidance molecules have long been postulated to control the development of the organization of neural connections into topographic maps. We review progress in identifying molecules required for mapping and the mechanisms by which they act, focusing on the visual system, the predominant model for map development. The Eph family of receptor tyrosine kinases and their ligands, the ephrins, remain the only molecules that meet all criteria for graded topographic guidance molecules, although others fulfill some criteria. Recent reports further define their modes of action and new roles for them, including EphB/ephrin-B control of dorsal-ventral mapping, bidirectional signaling of EphAs/ephrin-As, bifunctional action of ephrins as attractants or repellents in a context-dependent manner, and complex interactions between multiple guidance molecules. In addition, spontaneous patterned neural activity has recently been shown to be required for map refinement during a brief critical period. We speculate on additional activities required for map development and suggest a synthesis of molecular and cellular mechanisms within the context of the complexities of map development.

The transcriptome of Populus in elevated CO2

The consequences of increasing atmospheric carbon dioxide for long-term adaptation of forest ecosystems remain uncertain, with virtually no studies undertaken at the genetic level. A global analysis using cDNA microarrays was conducted following 6 yr exposure of Populus x euramericana (clone I-214) to elevated [CO(2)] in a FACE (free-air CO(2) enrichment) experiment. Gene expression was sensitive to elevated [CO(2)] but the response depended on the developmental age of the leaves, and < 50 transcripts differed significantly between different CO(2) environments. For young leaves most differentially expressed genes were upregulated in elevated [CO(2)], while in semimature leaves most were downregulated in elevated [CO(2)]. For transcripts related only to the small subunit of Rubisco, upregulation in LPI 3 and downregulation in LPI 6 leaves in elevated CO(2) was confirmed by anova. Similar patterns of gene expression for young leaves were also confirmed independently across year 3 and year 6 microarray data, and using real-time RT-PCR. This study provides the first clues to the long-term genetic expression changes that may occur during long-term plant response to elevated CO(2).

Identification of regulatory targets of tissue-specific transcription factors: application to retina-specific gene regulation

Identification of tissue-specific gene regulatory networks can yield insights into the molecular basis of a tissue's development, function and pathology. Here, we present a computational approach designed to identify potential regulatory target genes of photoreceptor cell-specific transcription factors (TFs). The approach is based on the hypothesis that genes related to the retina in terms of expression, disease and/or function are more likely to be the targets of retina-specific TFs than other genes. A list of genes that are preferentially expressed in retina was obtained by integrating expressed sequence tag, SAGE and microarray datasets. The regulatory targets of retina-specific TFs are enriched in this set of retina-related genes. A Bayesian approach was employed to integrate information about binding site location relative to a gene's transcription start site. Our method was applied to three retina-specific TFs, CRX, NRL and NR2E3, and a number of potential targets were predicted. To experimentally assess the validity of the bioinformatic predictions, mobility shift, transient transfection and chromatin immunoprecipitation assays were performed with five predicted CRX targets, and the results were suggestive of CRX regulation in 5/5, 3/5 and 4/5 cases, respectively. Together, these experiments strongly suggest that RP1, GUCY2D, ABCA4 are novel targets of CRX.

Transcriptional factors involved in photoreceptor differentiation

Regenerative medicine constitutes a potentially promising therapy for blind people suffering from retinal degenerative diseases such as retinitis pigmentosa and age-related macular degeneration. For the realization of retinal regeneration, it is necessary to establish 1) a method to produce functional photoreceptor cells in vitro and 2) successful transplantation of the donor cells to connect their axons to the recipient secondary neurons so that they can function properly. The results of experimental transplantation of human retinal photoreceptor cells from cadaveric eyes or of fetal retinal cells into the retina of RP patients have not been satisfactory, but encouraging enough to indicate that the transplantation of developing retinal cells may have beneficial results. Recently, attempts have been made to generate photoreceptor-like cells from stem cells, but it remains to be seen whether they are in fact photoreceptor cells. It is therefore important to fully understand the mechanisms involved in the development of these cells, and to characterize them not only by transcriptome but also by functional analysis.

Increased measurement accuracy for sequence-verified microarray probes

Microarrays have been extensively used to investigate genome-wide expression patterns. Although this technology has been tremendously successful, it has suffered from suboptimal individual measurement precision. Significant improvements in this respect have been recently made. In an effort to further explore the underlying variability, we have attempted to globally assess the accuracy of individual probe sequences used to query gene expression. For mammalian Affymetrix microarrays, we identify an unexpectedly large number of probes (greater than 19% of the probes on each platform) that do not correspond to their appropriate mRNA reference sequence (RefSeq). Compared with data derived from inaccurate probes, we find that data derived from sequence-verified probes show 1) increased precision in technical replicates, 2) increased accuracy translating data from one generation microarray to another, 3) increased accuracy translating data from oligonucleotide to cDNA microarrays, and 4) improved capture of biological information in human clinical specimens. The logical conclusion of this work is that probes containing the most reliable sequence information provide the most accurate results. Our data reveal that the identification and removal of inaccurate probes can significantly improve this technology.

Valproic acid teratogenicity: a toxicogenomics approach

Embryonic development is a highly coordinated set of processes that depend on hierarchies of signaling and gene regulatory networks, and the disruption of such networks may underlie many cases of chemically induced birth defects. The antiepileptic drug valproic acid (VPA) is a potent inducer of neural tube defects (NTDs) in human and mouse embryos. As with many other developmental toxicants however, the mechanism of VPA teratogenicity is unknown. Using microarray analysis, we compared the global gene expression responses to VPA in mouse embryos during the critical stages of teratogen action in vivo with those in cultured P19 embryocarcinoma cells in vitro. Among the identified VPA-responsive genes, some have been associated previously with NTDs or VPA effects [vinculin, metallothioneins 1 and 2 (Mt1, Mt2), keratin 1-18 (Krt1-18)], whereas others provide novel putative VPA targets, some of which are associated with processes relevant to neural tube formation and closure [transgelin 2 (Tagln2), thyroid hormone receptor interacting protein 6, galectin-1 (Lgals1), inhibitor of DNA binding 1 (Idb1), fatty acid synthase (Fasn), annexins A5 and A11 (Anxa5, Anxa11)], or with VPA effects or known molecular actions of VPA (Lgals1, Mt1, Mt2, Id1, Fasn, Anxa5, Anxa11, Krt1-18). A subset of genes with a transcriptional response to VPA that is similar in embryos and the cell model can be evaluated as potential biomarkers for VPA-induced teratogenicity that could be exploited directly in P19 cell-based in vitro assays. As several of the identified genes may be activated or repressed through a pathway of histone deacetylase (HDAC) inhibition and specificity protein 1 activation, our data support a role of HDAC as an important molecular target of VPA action in vivo.

Large-scale identification and characterization of genes with asymmetric expression patterns in the developing chick retina

To understand the molecular basis of topographic retinotectal projection, an overall view of the asymmetrically expressed molecules in the developing retina is needed. We performed a large-scale screening using restriction landmark cDNA scanning (RLCS) in the embryonic day 8 (E8) chick retina. RLCS is a cDNA display system, in which a large number of cDNA species are displayed as two-dimensional spots with intensities reflecting their expression levels as mRNA. We searched for spots that gave different signal intensities between the nasal and temporal retinas or between the dorsal and ventral retinas, and detected about 200 spots that were preferential on one side in the retina. The asymmetric expression of each gene was verified by Northern blotting and in situ hybridization. By subsequent analyses using molecular cloning, DNA sequencing, and database searching, 33 asymmetric molecules along the nasotemporal (N-T) axis and 20 along the dorsoventral (D-V) axis were identified. These included transcription factors, secretory factors, transmembrane proteins, and intracellular proteins with various putative functions. Their expression profiles revealed by in situ hybridization are highly diverse and individual. Moreover, many of them begin to be expressed in the retina from the early developmental stages, suggesting that they are implicated in the establishment and maintenance of regional specificity in the developing retina. The molecular repertoire revealed by this work will provide candidates for future studies to elucidate the molecular mechanisms of topographic retinotectal map formation.

Gene expression profiling of mammary gland development reveals putative roles for death receptors and immune mediators in post-lactational regression

INTRODUCTION

In order to gain a better understanding of the molecular processes that underlie apoptosis and tissue regression in mammary gland, we undertook a large-scale analysis of transcriptional changes during the mouse mammary pregnancy cycle, with emphasis on the transition from lactation to involution.

METHOD

Affymetrix microarrays, representing 8618 genes, were used to compare mammary tissue from 12 time points (one virgin, three gestation, three lactation and five involution stages). Six animals were used for each time point. Common patterns of gene expression across all time points were identified and related to biological function.

RESULTS

The majority of significantly induced genes in involution were also differentially regulated at earlier stages in the pregnancy cycle. This included a marked increase in inflammatory mediators during involution and at parturition, which correlated with leukaemia inhibitory factor-Stat3 (signal transducer and activator of signalling-3) signalling. Before involution, expected increases in cell proliferation, biosynthesis and metabolism-related genes were observed. During involution, the first 24 hours after weaning was characterized by a transient increase in expression of components of the death receptor pathways of apoptosis, inflammatory cytokines and acute phase response genes. After 24 hours, regulators of intrinsic apoptosis were induced in conjunction with markers of phagocyte activity, matrix proteases, suppressors of neutrophils and soluble components of specific and innate immunity.

CONCLUSION

We provide a resource of mouse mammary gene expression data for download or online analysis. Here we highlight the sequential induction of distinct apoptosis pathways in involution and the stimulation of immunomodulatory signals, which probably suppress the potentially damaging effects of a cellular inflammatory response while maintaining an appropriate antimicrobial and phagocytic environment.