Spatiotemporal Gene Expression Research Articles

Live-cell, temporal gene expression analysis of osteogenic differentiation in adipose-derived stem cells.

Adipose-derived stem cells (ASCs) are a widely investigated type of mesenchymal stem cells with great potential for musculoskeletal regeneration. However, the use of ASCs is complicated by their cellular heterogeneity, which exists at both the population and single-cell levels. This study demonstrates a live-cell assay to investigate gene expression in ASCs undergoing osteogenesis using fluorescently tagged DNA hybridization probes called molecular beacons. A molecular beacon was designed to target the mRNA sequence for alkaline phosphatase (ALPL), a gene characteristically expressed during early osteogenesis. The percentage of cells expressing this gene in a population was monitored daily to quantify the uniformity of the differentiation process. Differentiating ASC populations were repeatedly measured in a nondestructive fashion over a 10-day period to obtain temporal gene expression data. Results showed consistent expression patterns for the investigated osteogenic genes in response to induction medium. Peak signal level, indicating when the most cells expressed ALPL at once, was observed on days 3-5. The differentiation response of sample populations was generally uniform when assessed on a well-by-well basis over time. The expression of alkaline phosphatase is consistent with previous studies of osteogenic differentiation, suggesting that molecular beacons are a viable means of monitoring the spatiotemporal gene expression of live, differentiating ASCs.

A de novo 163 kb interstitial 1q44 microdeletion in a boy with thin corpus callosum, psychomotor delay and seizures

The 1q44 deletion syndrome has shown to be a recognizable phenotype with developmental delay, short stature and corpus callosum abnormalities as relatively consistent features. However, the disorder is still clinically heterogeneous and a genotype–phenotype correlation has been challenging to establish. In particular, a delineation of a critical region for the corpus callosum development has turned out to be difficult, and many candidate genes have been proposed. We present here a patient boy with a clinical picture of the 1q44 deletion syndrome, including a thin corpus callosum, and a small de novo 1q44 deletion. The deletion spans a maximum of 163 kb, a region which only contains the two genes FAM36A and HNRNPU. This finding supports the previously suggested hypothesis that the HNRNPU is an essential gene to the development of corpus callosum. However, as patients with deletions outside this interval also have been reported to have corpus callosum abnormalities, other mechanisms are probably also involved. We also identified two conserved non-coding regions in the deleted region of the patient, and speculate that also other elements interfere with the complex interplay and spatiotemporal gene expression during embryonic development.

Transcriptional regulation of phosphate acquisition by higher plants.

Phosphorus (P), an essential macronutrient required for plant growth and development, is often limiting in natural and agro-climatic environments. To cope with heterogeneous or low phosphate (Pi) availability, plants have evolved an array of adaptive responses facilitating optimal acquisition and distribution of Pi. The root system plays a pivotal role in Pi-deficiency-mediated adaptive responses that are regulated by a complex interplay of systemic and local Pi sensing. Cross-talk with sugar, phytohormones, and other nutrient signaling pathways further highlight the intricacies involved in maintaining Pi homeostasis. Transcriptional regulation of Pi-starvation responses is particularly intriguing and involves a host of transcription factors (TFs). Although PHR1 of Arabidopsis is an extensively studied MYB TF regulating subset of Pi-starvation responses, it is not induced during Pi deprivation. Genome-wide analyses of Arabidopsis have shown that low Pi stress triggers spatiotemporal expression of several genes encoding different TFs. Functional characterization of some of these TFs reveals their diverse roles in regulating root system architecture, and acquisition and utilization of Pi. Some of the TFs are also involved in phytohormone-mediated root responses to Pi starvation. The biological roles of these TFs in transcriptional regulation of Pi homeostasis in model plants Arabidopsis thaliana and Oryza sativa are presented in this review.

Evolution of bract development and B‐class MADS box gene expression in petaloid bracts of Cornus s. l. (Cornaceae)

Despite increasing interest in the molecular mechanisms of floral diversity, few studies have investigated the developmental and genetic bases of petaloid bracts. This study examined morphological patterns of bract initiation and expression patterns of B-class MADS-box genes in bracts of several Cornus species. We suggest that petaloid bracts in this genus may not share a single evolutionary origin. Developmental pathways of bracts and spatiotemporal expression of B-class genes in bracts and flowers were examined for four closely related dogwood species. Divergent morphological progressions and gene expression patterns were found in the two sister lineages with petaloid bracts, represented by Cornus florida and Cornus canadensis. Phylogeny-based analysis identified developmental and gene expression changes that are correlated with the evolution of petaloid bracts in C. florida and C. canadensis. Our data support the existence of independent evolutionary origins of petaloid bracts in C. canadensis and C. florida. Additionally, we suggest that functional transference within B-class gene families may have contributed to the origin of bract petaloidy in C. florida. However, the underlying mechanisms of petaloid bract development likely differ between C. florida and C. canadensis. In the future this hypothesis can be tested by functional analyses of Cornus B-class genes.

The role of microRNAs in mammalian oocytes and embryos

Advanced genomic analysis has revealed an enormous inventory of non-coding RNAs (ncRNAs), which are functionally important at transcriptional and post-transcriptional level for different cellular processes. Among the ncRNAs, microRNAs (miRNAs) have recently been highlighted extensively for their pivotal role in disease, fertility and development through post-transcriptional regulation of gene expression. The presence and spatio-temporal expression of miRNAs and miRNA processing machinery genes in oocytes and preimplantation embryos has evidenced the involvement of miRNAs for growth and maturation of mammalian oocytes, early embryonic development, stem cell lineage differentiation and implantation. Therefore, this article aims to highlight primary evidences on the importance of miRNAs and their mediated translational reprogramming in the physiology and development of mammalian oocytes and embryos.

Characterization and importance of microRNAs in mammalian gonadal functions

Recent progress in high throughput sequencing and bioinformatic analysis and other biochemical methods have fuelled our appreciation for the important role of microRNAs (miRNAs) in disease, fertility and development. These tiny RNAs were found to be potentially involved in various aspects of cellular processes of reproductive tissues by posttranscriptional regulation of protein coding genes. Mammalian gonads which exhibit strictly regulated spatiotemporal gene expression patterns are also known to express unique sets of miRNAs and genes involved in the miRNA biogenetic pathway. Studies on miRNAs and their associated processing enzymes have evidenced the contribution of these small regulatory RNAs to germ cell differentiation, post-meiotic male germ cell function and growth, and development and maturation of oocytes through pertaining tightly regulated gene expression. The existence, preferential and temporal expression of miRNAs and their processing machinery genes in different stages of testicular and ovarian cellular development have evidenced the potential role of miRNAs in testicular and ovarian physiology. MiRNAs are also found to be associated with functional regulation of gonadal somatic cells, namely Leydig cells and Sertoli cells in testis and granulosa cells/cumulus cells in the ovary in steroid synthesis. Here, we review the recent works on the involvement and diverse roles of miRNAs in the development and physiology of gonadal cells in mammalian reproduction.

Postnatal alteration of collapsin response mediator protein 4 mRNA expression in the mouse brain

Collapsin response mediator protein 4 (CRMP4) is a molecular marker for immature neurons but only limited information is available on the spatiotemporal gene expression changes of Crmp4 in the developing rodent. In the present study, the variation of CRMP4 mRNA expression in the mouse brain was investigated from postnatal day (PD) 0 (the day of birth) to adulthood by in situ hybridization. The hybridization signals were broadly detected on PD0 and regional changes in expression during development were noted. Expression patterns of CRMP4 mRNA were classified into the following three types: (i) signals that were strongest on PD0 or PD7, weak or undetectable on PD14, and absent in adulthood: this pattern was observed in most brain areas; (ii) signals that were first detected on PD0 or PD7 and persisted into adulthood: this pattern was seen in the dentate gyrus and subventricular zone of the olfactory bulb (OB); and (iii) signals that were strongest on PD0 and decreased gradually with age but were still detectable in adulthood: this pattern was identified for the first time in the mitral cell layer of the OB. Analysis using quantitative real-time RT-PCR confirmed higher expression of CRMP4 mRNA in the OB than in other adult brain regions. The persistence of CRMP4 mRNA in the adult OB, including the mitral cell layer, suggests the possibility of both neurogenetic and non-neurogenetic functional roles of CRMP4 in this region.

Learning Sparse Representations for Fruit-Fly Gene Expression Pattern Image Annotation and Retrieval

BackgroundFruit fly embryogenesis is one of the best understood animal development systems, and the spatiotemporal gene expression dynamics in this process are captured by digital images. Analysis of these high-throughput images will provide novel insights into the functions, interactions, and networks of animal genes governing development. To facilitate comparative analysis, web-based interfaces have been developed to conduct image retrieval based on body part keywords and images. Currently, the keyword annotation of spatiotemporal gene expression patterns is conducted manually. However, this manual practice does not scale with the continuously expanding collection of images. In addition, existing image retrieval systems based on the expression patterns may be made more accurate using keywords.ResultsIn this article, we adapt advanced data mining and computer vision techniques to address the key challenges in annotating and retrieving fruit fly gene expression pattern images. To boost the performance of image annotation and retrieval, we propose representations integrating spatial information and sparse features, overcoming the limitations of prior schemes.ConclusionsWe perform systematic experimental studies to evaluate the proposed schemes in comparison with current methods. Experimental results indicate that the integration of spatial information and sparse features lead to consistent performance improvement in image annotation, while for the task of retrieval, sparse features alone yields better results.

Cloning and expression analysis of Fgf5, 6 and 7 during early chick development

FGFs with similar sequences can play different roles depending on the model organisms examined. Determining these roles requires knowledge of spatio-temporal Fgf gene expression patterns. In this study, we report the cloning of chick Fgf5, 6 and 7, and examine their gene expression patterns by whole mount in situ hybridization. We show that Fgf5’s spatio-temporally restricted expression pattern indicates a potentially novel role during inner ear development. Fgf6 and Fgf7, although belonging to different subfamilies with diverged sequences, are expressed in similar patterns within the mesoderm. Alignment of protein sequences and phylogenetic analysis demonstrate that FGF5 and FGF6 are highly conserved between chick, human, mouse and zebrafish. FGF7 is similarly conserved except for the zebrafish, which has considerably diverged.

Transcriptional Activity of Neural Retina Leucine Zipper (Nrl) Is Regulated by c-Jun N-Terminal Kinase and Tip60 during Retina Development

Neural retina leucine zipper (Nrl), a key basic motif leucine zipper (bZIP) transcription factor, modulates rod photoreceptor differentiation by activating rod-specific target genes. In searching for factors that might couple with Nrl to modulate its transcriptional activity through posttranslational modification, we observed the novel interactions of Nrl with c-Jun N-terminal kinase 1 (JNK1) and HIV Tat-interacting protein 60 (Tip60). JNK1 directly interacted with and phosphorylated Nrl at serine 50, which enhanced Nrl transcriptional activity on the rhodopsin and Ppp2r5c promoters. Use of an inactive JNK1 mutant or treatment with a JNK inhibitor (SP600125) significantly reduced JNK1-mediated phosphorylation and transcriptional activity of Nrl in cultured retinal explants. We also found that Nrl activated rhodopsin and Ppp2r5c transcription by recruiting Tip60 to promote histone H3/H4 acetylation. The binding affinity of phospho-Nrl for Tip60 was significantly greater than that of the unphosphorylated Nrl. Thus, the histone acetyltransferase-containing Tip60 behaved as a coactivator in the Nrl-dependent transcriptional regulation of the rhodopsin and Ppp2r5c genes in the developing mouse retina. A transcriptional network of interactive proteins, including Nrl, JNK1, and Tip60, may be required to precisely control spatiotemporal photoreceptor-specific gene expression during retinal development.

Krüppel-like factor 4 (KLF4) promotes the survival of natural killer cells and maintains the number of conventional dendritic cells in the spleen

The development and survival of NK cells rely on a complex, spatiotemporal gene expression pattern regulated by specific transcription factors in NK cells and tissue-specific microenvironments supported by hematopoietic cells. Here, we show that somatic deletion of the KLF4 gene, using inducible and lineage-specific cre-transgenic mice, leads to a significant reduction of NK cells (NK1.1(+) TCR-β(-)) in the blood and spleen but not in the BM, liver, or LNs. Functional and immunophenotypic analyses revealed increased apoptosis of CD27(+/-) CD11b(+) NK cells in the spleen of KLF4-deficient mice, although remaining NK cells were able to lyse tumor target cells and produce IFN-γ. A normal recovery of adoptively transferred KLF4-deficient NK cells in WT hosts suggested that the survival defect was not intrinsic of NK cells. However, BM chimeras using KLF4-deficient mice as donors indicated that reduced survival of NK cells depended on BM-derived hematopoietic cells in the spleen. The number of CD11c(hi) DCs, which are known to support NK cell survival, was reduced significantly in the spleen of KLF4-deficient mice, likely a result of a lower number of precDC progenitor cells in this tissue. Taken together, our data suggest that the pluripotency-associated gene KLF4 is required for the maintenance of DCs in the spleen and consequently, survival of differentiated NK cells in this tissue.

Expression of VjbR under Nutrient Limitation Conditions Is Regulated at the Post-Transcriptional Level by Specific Acidic pH Values and Urocanic Acid

VjbR is a LuxR homolog that regulates transcription of many genes including important virulence determinants of the facultative intracellular pathogen Brucella abortus. This transcription factor belongs to a family of regulators that participate in a cell-cell communication process called quorum sensing, which enables bacteria to respond to changes in cell population density by monitoring concentration of self produced autoinducer molecules. Unlike almost all other LuxR-type proteins, VjbR binds to DNA and activates transcription in the absence of any autoinducer signal. To investigate the mechanisms by which Brucella induces VjbR-mediated transcriptional activation, and to determine how inappropriate spatio-temporal expression of the VjbR target genes is prevented, we focused on the study of expression of vjbR itself. By assaying different parameters related to the intracellular lifestyle of Brucella, we identified a restricted set of conditions that triggers VjbR protein expression. Such conditions required the convergence of two signals of different nature: a specific pH value of 5.5 and the presence of urocanic acid, a metabolite involved in the connection between virulence and metabolism of Brucella. In addition, we also observed an urocanic acid, pH-dependent expression of RibH2 and VirB7, two additional intracellular survival-related proteins of Brucella. Analysis of promoter activities and determination of mRNA levels demonstrated that the urocanic acid-dependent mechanisms that induced expression of VjbR, RibH2, and VirB7 act at the post-transcriptional level. Taken together, our findings support a model whereby Brucella induces VjbR-mediated transcription by modulating expression of VjbR in response to specific signals related to the changing environment encountered within the host.

A flexible integrative approach based on random forest improves prediction of transcription factor binding sites

Transcription factor binding sites (TFBSs) are DNA sequences of 6–15 base pairs. Interaction of these TFBSs with transcription factors (TFs) is largely responsible for most spatiotemporal gene expression patterns. Here, we evaluate to what extent sequence-based prediction of TFBSs can be improved by taking into account the positional dependencies of nucleotides (NPDs) and the nucleotide sequence-dependent structure of DNA. We make use of the random forest algorithm to flexibly exploit both types of information. Results in this study show that both the structural method and the NPD method can be valuable for the prediction of TFBSs. Moreover, their predictive values seem to be complementary, even to the widely used position weight matrix (PWM) method. This led us to combine all three methods. Results obtained for five eukaryotic TFs with different DNA-binding domains show that our method improves classification accuracy for all five eukaryotic TFs compared with other approaches. Additionally, we contrast the results of seven smaller prokaryotic sets with high-quality data and show that with the use of high-quality data we can significantly improve prediction performance. Models developed in this study can be of great use for gaining insight into the mechanisms of TF binding.

A Versatile Protocol for mRNA Electroporation of Xenopus laevis Embryos

Xenopus laevis is an ideal model system for investigating the mechanisms of pattern formation. The ability to express exogenous mRNA or introduce morpholinos into cleavage-stage Xenopus embryos has allowed gain- and loss-of-function experiments that reveal molecular-genetic control of development and regeneration. However, injection of mRNAs into cleavage-stage embryos provides limited spatio-temporal control: It is difficult to limit targeting to small regions (e.g., inducing foci of expression) and the fate map does not facilitate targeting some tissues, such as those of the tail. Likewise, early injection can result in unwanted developmental defects because mRNA can be translated long before the desired time point. These are especially important limitations when studying developmental and regenerative processes during the gastrula to tailbud stages. Although transgenic techniques allow precise control over spatio-temporal expression of genes when the appropriate promoter is available, the process of creating stable transgenic animals is time-consuming. Electroporation provides an alternative method for delivering mRNA and other nucleic acids, enabling the targeting of single cells or groups of cells at any stage of development. This protocol describes detailed electroporation parameters for the transfection of mRNA into a wide range of tissues in embryos at gastrula to tailbud stages, with high efficiency and expression as early as 4 h post electroporation.

Integration of an abdominal Hox complex with Pax2 yields cell-specific EGF secretion from Drosophila sensory precursor cells

Cis-regulatory modules (CRMs) ensure specific developmental outcomes by mediating both proper spatiotemporal gene expression patterns and appropriate transcriptional levels. In Drosophila, the precise transcriptional control of the serine protease rhomboid regulates EGF signaling to specify distinct cell types. Recently, we identified a CRM that activates rhomboid expression and thereby EGF secretion from a subset of abdominal sensory organ precursor cells (SOPs) to induce an appropriate number of lipid-processing cells called oenocytes. Here, we use scanning mutagenesis coupled with reporter assays, biochemistry and genetics to dissect the transcriptional mechanisms regulating SOP-specific rhomboid activation. Our results show that proper spatial activity of the rhomboid CRM is dependent upon direct integration of the abdomen-specific Hox factor Abdominal-A and the SOP-restricted Pax2 factor. In addition, we show that the Extradenticle and Homothorax Hox co-factors are differentially integrated on the rhomboid CRM by abdominal versus thoracic Hox proteins in the presence of Pax2. Last, we show that Abdominal-A uses both Pax2-dependent and Pax2-independent mechanisms to stimulate rhomboid CRM activity to induce proper oenocyte numbers. Thus, these data demonstrate how a CRM integrates Hox and neural transcriptional inputs to regulate the appropriate spatial pattern and levels of EGF secretion to specify an essential cell fate.

A Machine Learning Approach for Identifying Novel Cell Type–Specific Transcriptional Regulators of Myogenesis

Transcriptional enhancers integrate the contributions of multiple classes of transcription factors (TFs) to orchestrate the myriad spatio-temporal gene expression programs that occur during development. A molecular understanding of enhancers with similar activities requires the identification of both their unique and their shared sequence features. To address this problem, we combined phylogenetic profiling with a DNA–based enhancer sequence classifier that analyzes the TF binding sites (TFBSs) governing the transcription of a co-expressed gene set. We first assembled a small number of enhancers that are active in Drosophila melanogaster muscle founder cells (FCs) and other mesodermal cell types. Using phylogenetic profiling, we increased the number of enhancers by incorporating orthologous but divergent sequences from other Drosophila species. Functional assays revealed that the diverged enhancer orthologs were active in largely similar patterns as their D. melanogaster counterparts, although there was extensive evolutionary shuffling of known TFBSs. We then built and trained a classifier using this enhancer set and identified additional related enhancers based on the presence or absence of known and putative TFBSs. Predicted FC enhancers were over-represented in proximity to known FC genes; and many of the TFBSs learned by the classifier were found to be critical for enhancer activity, including POU homeodomain, Myb, Ets, Forkhead, and T-box motifs. Empirical testing also revealed that the T-box TF encoded by org-1 is a previously uncharacterized regulator of muscle cell identity. Finally, we found extensive diversity in the composition of TFBSs within known FC enhancers, suggesting that motif combinatorics plays an essential role in the cellular specificity exhibited by such enhancers. In summary, machine learning combined with evolutionary sequence analysis is useful for recognizing novel TFBSs and for facilitating the identification of cognate TFs that coordinate cell type–specific developmental gene expression patterns.

Nematostella vectensis achaete-scute homolog NvashA regulates embryonic ectodermal neurogenesis and represents an ancient component of the metazoan neural specification pathway

achaete-scute homologs (ash) regulate neural development in all bilaterian model animals indicating that they represent a component of the ancestral neurogenic pathway. We test this by investigating four ash genes during development of a basal metazoan, the cnidarian sea anemone Nematostella vectensis. Spatiotemporal expression of ash genes in the early embryo and larval stages suggests that they regulate neurogenesis. More specifically, NvashA is co-expressed with neural genes in the embryonic ectoderm. Knockdown of NvashA results in decreased expression of eight neural markers, including the six novel neural targets identified here. Conversely, overexpression of NvashA induces increased expression of all eight genes, but only within their normal axial domains. Overexpression of NvashB-D differentially increases expression of NvashA targets. The expression patterns and differential ability of ash genes to regulate neural gene expression reveals surprising molecular complexity in these 'simple' animals. These data suggest that achaete-scute homologs functioned in the ancestral metazoan neurogenic pathway and provide a foundation to investigate further the evolution of neurogenesis and the origin of complex central nervous systems.

A Wnt-bmp feedback circuit controls intertissue signaling dynamics in tooth organogenesis.

Many vertebrate organs form through the sequential and reciprocal exchange of signaling molecules between juxtaposed epithelial and mesenchymal tissues. We undertook a systems biology approach that combined the generation and analysis of large-scale spatiotemporal gene expression data with mouse genetic experiments to gain insight into the mechanisms that control epithelial-mesenchymal signaling interactions in the developing mouse molar tooth. We showed that the shift in instructive signaling potential from dental epithelium to dental mesenchyme was accompanied by temporally coordinated genome-wide changes in gene expression in both compartments. To identify the mechanism responsible, we developed a probabilistic technique that integrates regulatory evidence from gene expression data and from the literature to reconstruct a gene regulatory network for the epithelial and mesenchymal compartments in early tooth development. By integrating these epithelial and mesenchymal gene regulatory networks through the action of diffusible extracellular signaling molecules, we identified a key epithelial-mesenchymal intertissue Wnt-Bmp (bone morphogenetic protein) feedback circuit. We then validated this circuit in vivo with compound genetic mutations in mice that disrupted this circuit. Moreover, mathematical modeling demonstrated that the structure of the circuit accounted for the observed reciprocal signaling dynamics. Thus, we have identified a critical signaling circuit that controls the coordinated genome-wide expression changes and reciprocal signaling molecule dynamics that occur in interacting epithelial and mesenchymal compartments during organogenesis.

Gene order and chromosome dynamics coordinate spatiotemporal gene expression during the bacterial growth cycle

In Escherichia coli crosstalk between DNA supercoiling, nucleoid-associated proteins and major RNA polymerase σ initiation factors regulates growth phase-dependent gene transcription. We show that the highly conserved spatial ordering of relevant genes along the chromosomal replichores largely corresponds both to their temporal expression patterns during growth and to an inferred gradient of DNA superhelical density from the origin to the terminus. Genes implicated in similar functions are related mainly in trans across the chromosomal replichores, whereas DNA-binding transcriptional regulators interact predominantly with targets in cis along the replichores. We also demonstrate that macrodomains (the individual structural partitions of the chromosome) are regulated differently. We infer that spatial and temporal variation of DNA superhelicity during the growth cycle coordinates oxygen and nutrient availability with global chromosome structure, thus providing a mechanistic insight into how the organization of a complete bacterial chromosome encodes a spatiotemporal program integrating DNA replication and global gene expression.

Comparative assessment of fungal cellobiohydrolase I richness and composition in cDNA generated using oligo(dT) primers or random hexamers

Understanding soil fungal distribution and activities, particularly at the level of gene expression, is important in unveiling mechanisms regulating their activities in situ. Recent identification of fungal genes involved in carbon cycling has provided the foundation for developing reverse-transcriptase PCR assays to monitor spatiotemporal gene expression patterns in soils and other complex microbial systems. The polyadenylated 3′ ends of eukaryotic mRNA transcripts enables the use of oligo(dT) primers for cDNA synthesis, but this can result in the overrepresentation of the 3′ end of transcripts in cDNA pools. In an effort to increase the uniformity of transcripts represented in cDNA pools, random hexamers have been used. The use of both priming methods is abundant in the literature, but we do not know how these methods perform relative to each other. We performed comparative richness and compositional analyses of the fungal glycosyl hydrolase family 7 cellobiohydrolase I gene cbhI amplified from soil cDNAs that had been generated using either oligo(dT) primers or random hexamers. Our results demonstrate that similar cbhI richness and composition were recovered using both approaches. Richness estimates and compositional profiles of cbhI sequence libraries generated from random hexamer-primed cDNA were more variable than from libraries generated from oligo(dT) primed cDNA. However, our overall results indicate that, on average, comparable richness and composition were recovered from soil cDNAs when either priming method was used.