An Integrative Developmental Genomics and Systems Biology Approach to Identify an In Vivo Sox Trio-Mediated Gene Regulatory Network in Murine Embryos.

Embryogenesis is an intricate process involving multiple genes and pathways. Some of the key transcription factors controlling specific cell types are the Sox trio, namely, Sox5, Sox6, and Sox9, which play crucial roles in organogenesis working in a concerted manner. Much however still needs to be learned about their combinatorial roles during this process. A developmental genomics and systems biology approach offers to complement the reductionist methodology of current developmental biology and provide a more comprehensive and integrated view of the interrelationships of complex regulatory networks that occur during organogenesis. By combining cell type-specific transcriptome analysis and in vivo ChIP-Seq of the Sox trio using mouse embryos, we provide evidence for the direct control of Sox5 and Sox6 by the transcriptional trio in the murine model and by Morpholino knockdown in zebrafish and demonstrate the novel role of Tgfb2, Fbxl18, and Tle3 in formation of Sox5, Sox6, and Sox9 dependent tissues. Concurrently, a complete embryonic gene regulatory network has been generated, identifying a wide repertoire of genes involved and controlled by the Sox trio in the intricate process of normal embryogenesis.

Genome wide binding (ChIP-Seq) of murine Bapx1 and Sox9 proteins in vivo and in vitro.

This work pertains to GEO submission GSE36672, in vivo and in vitro genome wide binding (ChIP-Seq) of Bapx1/Nkx3.2 and Sox9 proteins. We have previously shown that data from a genome wide binding assay combined with transcriptional profiling is an insightful means to divulge the mechanisms directing cell type specification and the generation of tissues and subsequent organs [1]. Our earlier work identified the role of the DNA-binding homeodomain containing protein Bapx1/Nkx3.2 in midgestation murine embryos. Microarray analysis of EGFP-tagged cells (both wildtype and null) was integrated using ChIP-Seq analysis of Bapx1/Nkx3.2 and Sox9 DNA-binding proteins in living tissue.

Gene expression profiles of Bapx1 expressing FACS sorted cells from wildtype and Bapx1-EGFP null mouse embryos.

The data described in this article refers to Chatterjee et al. (2015) "In vivo genome-wide analysis of multiple tissues identifies gene regulatory networks, novel functions and downstream regulatory genes for Bapx1 and its co-regulation with Sox9 in the mammalian vertebral column" (GEO GSE35649) [1]. Transcriptional profiling combined with genome wide binding data is a powerful tool to elucidate the molecular mechanism behind vertebrate organogenesis. It also helps to uncover multiple roles of a single gene in different organs. In the above mentioned report we reveal the function of the homeobox gene Bapx1 during the embryogenesis of five distinct organs (vertebral column, spleen, gut, forelimb and hindlimb) at a relevant developmental stage (E12.5), microarray analysis of isolated wildtype and mutant cells in is compared in conjunction with ChIP-Seq analysis. We also analyzed the development of the vertebral column by comparing microarray and ChIP-Seq data for Bapx1 with similarly generated data sets for Sox9 to generate a gene regulatory network controlling various facets of the organogenesis.

In vivo genome-wide analysis of multiple tissues identifies gene regulatory networks, novel functions and downstream regulatory genes for Bapx1 and its co-regulation with Sox9 in the mammalian vertebral column.

BACKGROUND: Vertebrate organogenesis is a highly complex process involving sequential cascades of transcription factor activation or repression. Interestingly a single developmental control gene can occasionally be essential for the morphogenesis and differentiation of tissues and organs arising from vastly disparate embryological lineages. RESULTS: Here we elucidated the role of the mammalian homeobox gene Bapx1 during the embryogenesis of five distinct organs at E12.5 - vertebral column, spleen, gut, forelimb and hindlimb - using expression profiling of sorted wildtype and mutant cells combined with genome wide binding site analysis. Furthermore we analyzed the development of the vertebral column at the molecular level by combining transcriptional profiling and genome wide binding data for Bapx1 with similarly generated data sets for Sox9 to assemble a detailed gene regulatory network revealing genes previously not reported to be controlled by either of these two transcription factors. CONCLUSIONS: The gene regulatory network appears to control cell fate decisions and morphogenesis in the vertebral column along with the prevention of premature chondrocyte differentiation thus providing a detailed molecular view of vertebral column development.

A conditional mouse line for lineage tracing of Sox9 loss-of-function cells using enhanced green fluorescent protein.

Traditionally, conditional knockout studies in mouse have utilized the Cre or Flpe technology to activate the expression of reporter genes such as lacZ or PLAP. Employing these reporter genes, however, requires tissue fixation. To make way for downstream in vivo or in vitro applications, we have inserted enhanced green fluorescent protein (EGFP) into the endogenous Sox9 locus and generated a novel conditional Sox9 null allele, by flanking the entire Sox9 coding region with loxP sites and inserting an EGFP reporter gene into the 3'-UTR allowing for EGFP to be expressed upon Sox9 loss of function yet under the control of the endogenous Sox9 promoter. Mating this new allele to any Cre-expressing line, the fate of Sox9 null cells can be traced in the cell type of interest in vivo or in vitro after fluorescence-activated cell sorting.

Mouse strain specific gene expression differences for illumina microarray expression profiling in embryos.

BACKGROUND: In the field of mouse genetics the advent of technologies like microarray based expression profiling dramatically increased data availability and sensitivity, yet these advanced methods are often vulnerable to the unavoidable heterogeneity of in vivo material and might therefore reflect differentially expressed genes between mouse strains of no relevance to a targeted experiment. The aim of this study was not to elaborate on the usefulness of microarray analysis in general, but to expand our knowledge regarding this potential "background noise" for the widely used Illumina microarray platform surpassing existing data which focused primarily on the adult sensory and nervous system, by analyzing patterns of gene expression at different embryonic stages using wild type strains and modern transgenic models of often non-isogenic backgrounds. RESULTS: Wild type embryos of 11 mouse strains commonly used in transgenic and molecular genetic studies at three developmental time points were subjected to Illumina microarray expression profiling in a strain-by-strain comparison. Our data robustly reflects known gene expression patterns during mid-gestation development. Decreasing diversity of the input tissue and/or increasing strain diversity raised the sensitivity of the array towards the genetic background. Consistent strain sensitivity of some probes was attributed to genetic polymorphisms or probe design related artifacts. CONCLUSION: Our study provides an extensive reference list of gene expression profiling background noise of value to anyone in the field of developmental biology and transgenic research performing microarray expression profiling with the widely used Illumina microarray platform. Probes identified as strain specific background noise further allow for microarray expression profiling on its own to be a valuable tool for establishing genealogies of mouse inbred strains.

New Bapx1(Cre-EGFP) mouse lines for lineage tracing and conditional knockout studies.

To gain insight into the roles of various genes in development and to circumvent embryonic lethality that hinders genetic studies, lineage tracing and conditional knockout techniques have been widely performed on mice using the increasing numbers of gene-targeted Cre mouse lines. Employing the internal ribosome entry site (IRES) and the 2A peptide multicistronic expression strategies, we report two new Bapx1 mouse lines with functional Bapx1 whereby Cre and enhanced green fluorescence protein (EGFP) are expressed discretely under the control of the Bapx1 promoter. These mouse lines, when mated with the Rosa26R-lacZ reporter line, can be used to trace the lineage of Bapx1-expressing cells whereas stage-specific, spatial expression of Bapx1 can be visualized by the EGFP fluorescence. In addition, both of our Bapx1(Cre-EGFP) mouse lines can be used to enrich for Bapx1-specific cells and also serve as effective conditional knockout tools to investigate gene functions in the skeleton and/or visceral organs.

Generation of mice with a novel conditional null allele of the Sox9 gene.

Sox9 is expressed in multiple tissues during mouse development and adulthood. Mutations in the Sox9 gene or changes in expression levels can be attributed to many congenital diseases. Heterozygous loss-of-function mutations in the human SOX9 gene cause Campomelic dysplasia, a semi-lethal skeletal malformation syndrome. Disruption of Sox9 by conventional gene targeting leads to perinatal lethality in heterozygous mice, hence hampering the feasibility to obtain the homozygous Sox9 null mice for in vivo functional studies. In this study, we generated a conditional allele of Sox9 (Sox9 ( tm4.Tlu )) by flanking exon 1 with loxP sites. Homozygous mice for the Sox9 ( tm4.Tlu ) allele (Sox9 ( flox/flox )) are viable, fertile and indistinguishable from wildtype (WT) mice, indicating that the Sox9 ( tm4.Tlu ) allele is a fully functional Sox9 allele. Furthermore, we demonstrated that Cre-mediated recombination using a Col2a1-Cre line resulted in specific ablation of Sox9 activity in cartilage tissues.

Comparison of IRES and F2A-based locus-specific multicistronic expression in stable mouse lines.

Efficient and stoichiometric expression of genes concatenated by bi- or multi-cistronic vectors has become an invaluable tool not only in basic biology to track and visualize proteins in vivo, but also for vaccine development and in the clinics for gene therapy. To adequately compare, in vivo, the effectiveness of two of the currently popular co-expression strategies - the internal ribosome entry site (IRES) derived from the picornavirus and the 2A peptide from the foot-and-mouth disease virus (FDMV) (F2A), we analyzed two locus-specific knock-in mouse lines co-expressing SRY-box containing gene 9 (Sox9) and enhanced green fluorescent protein (EGFP) linked by the IRES (Sox9(IRES-EGFP)) or the F2A (Sox9(F2A-EGFP)) sequence. Both the constructs expressed Sox9 and EGFP proteins in the appropriate Sox9 expression domains, with the IRES construct expressing reduced levels of EGFP compared to that of the F2A. The latter, on the other hand, produced about 42.2% Sox9-EGFP fusion protein, reflecting an inefficient ribosome 'skipping' mechanism. To investigate if the discrepancy in the 'skipping' process was locus-dependent, we further analyzed the FLAG(3)-Bapx1(F2A-EGFP) mouse line and found similar levels of fusion protein being produced. To assess if EGFP was hindering the 'skipping' mechanism, we examined another mouse line co-expressing Bagpipe homeobox gene 1 homolog (Bapx1), Cre recombinase and EGFP (Bapx1(F2A-Cre-F2A-EGFP)). While the 'skipping' was highly efficient between Bapx1 and Cre, the 'skipping' between Cre and EGFP was highly inefficient. We have thus demonstrated in our comparison study that the efficient and close to equivalent expression of genes linked by F2A is achievable in stable mouse lines, but the EGFP reporter may cause undesirable inhibition of the 'skipping' at the F2A sequence. Hence, the use of other reporter genes should be explored when utilizing F2A peptides.

A more cost effective and rapid high percentage germ-line transmitting chimeric mouse generation procedure via microinjection of 2-cell, 4-cell, and 8-cell embryos with ES and iPS cells.

The long-standing traditional method of delivering embryonic stem (ES) cells adjacent to the inner cell mass (ICM) of blastocysts to generate chimeras improved with the advent of laser- or Piezo assisted 8-cell embryo microinjection. Building on this technology but omitting either the laser or the Piezo to penetrate the zona pellucida and making use of earlier embryonic stages (2-cell and 4-cell), we were able to significantly speed up and economize our ES cell microinjection and chimera production throughput. We demonstrate here that embryonic (ES) and induced pluripotent stem (iPS) cells can stay fully pluripotent when delivered into 2-cell- and 4-cell-stage embryos, long before they would naturally be incorporated into the ICM of a blastocyst (E3.5) and give rise to high percentage and germline transmitting chimeras.

Ultrasensitive cell-based bioassay for the measurement of global estrogenic activity of flavonoid mixtures revealing additive, restrictive, and enhanced actions in binary and higher order combinations.

Flavonoids present in food, botanicals, and body fluids occur as complex mixtures, and data on their combinatorial estrogenic effects are sparse. Human cell lines that permanently express estrogen receptor (ER) alpha and ERbeta proteins were developed for the measurement of the global estrogenicity of flavonoids in such complex mixtures. The presence of estrogenic ligands, known and unknown, in these mixtures can be detected by activation of an ER-driven luciferase reporter gene. We also examined the effect of hydroxylation on the estrogenic activities of four common flavonoids-apigenin, kaempferol, luteolin, and quercetin, alone and in combination. An inverse relationship was observed between the number of hydroxyl groups in flavonoids and ERalpha bioactivity. When submaximal doses of apigenin, luteolin, kaempferol, genistein, and estradiol were combined in binary and higher order mixtures, the experimental estrogenic effects matched those obtained by summing effects extrapolated from dose-response curves of individual compounds. The estrogenic activities of mixtures containing quercetin were observed to deviate from additivity, suggesting that it was a partial agonist/antagonist. Our assay reveals superagonistic, additive, and antagonistic ERalpha or ERbeta actions of flavonoids and adds to our understanding of the estrogenic effects of phytoestrogens in complex mixtures.

New estrogenic prenylflavone from Epimedium brevicornum inhibits the growth of breast cancer cells.

Estrogens maintain female sexual health. The hormone also drives the growth of estrogen receptor (ER) positive breast tumors, and ER modulators, like tamoxifen, are used to reduce tumor recurrence. To identify phytoestrogens with possible health benefits, we screened several Traditional Chinese Medicines and encountered an extract from the leaves of Epimedium brevicornum (EB), with strong (EC50: 1.3 microg/mL) and specific ER-stimulatory activity. It increased estrogen-responsive human breast cancer cell proliferation at low doses, but paradoxically caused profound inhibition of growth at higher doses. Using bioassay-guided fractionation, we isolated and characterized a new prenylflavone, breviflavone B, which exerted biphasic stimulatory and inhibitory effects on breast cancer cell proliferation, mimicking the effects of EB. In contrast to estradiol and genistein, high doses (> 2 microM) of breviflavone B almost eliminated ERalpha protein; a process that may be mediated through increased proteasome degradation. Pre-clinical studies are needed to explore whether these prenylflavones are of value in estrogen-deficiency states and for prophylaxis of breast cancer.