Identification of neural genes usingXenopus DNA microarrays

Ferritin H subunit gene is specifically expressed in melanophore precursor-derived white pigment cells in which reflecting platelets are formed from stage II melanosomes in the periodic albino mutant of Xenopus laevis

“White pigment cells” are derived from melanophore precursors and contain both melanophore-specific and iridophore-specific pigment organelles. Whereas melanophores differentiate in the wild type regenerating tail, white pigment cells appear in the regenerating tail in the periodic albino mutant (a^p/a^p) of Xenopus laevis. The localization and density of white pigment cells in the mutant regenerating tail are similar to those of melanophores in the wild type regenerating tail. Here, white pigment cells in the mutant regenerating tail have been compared with melanophores in the wild type regenerating tail in the presence of phenylthiourea (PTU), which inhibits melanosome maturation in melanophores but does not affect reflecting platelet formation in white pigment cells. Ultrastructural analysis shows that reflecting platelet formation in white pigment cells is different from that in iridophores. Reflecting platelets in iridophores are formed from spherical vesicles with electron-dense material, whereas they are formed from stage II melanosomes characteristic of melanophore precursors in white pigment cells. Ultrastructural features of pigment organelles, except reflecting platelets, are similar between mutant melanophores and white pigment cells. In an attempt to identify specific genes in white pigment cells, a subtracted cDNA library enriched for mutant cDNAs has been prepared. Subtracted cDNA fragments have been cloned and selected by whole mount in situ hybridization. Among cDNA fragments examined so far, the ferritin H subunit gene is specifically expressed in white pigment cells, but not in melanophores. Pigment organellogenesis and specific gene expression in white pigment cells are also discussed.

Changes in soybean global gene expression after application of lipo-chitooligosaccharide from Bradyrhizobium japonicum under sub-optimal temperature

Lipo-chitooligosaccharides (LCOs), signal compounds produced by N(2)-fixing rhizobacteria after isoflavone induction, initiate nodule formation in host legumes. Given LCOs' structural similarity to pathogen-response-eliciting chitin oligomers, foliar application of LCOs was tested for ability to induce stress-related genes under optimal growth conditions. In order to study the effects of LCO foliar spray under stressed conditions, soybean (Glycine max) seedlings grown at optimal temperature were transferred to sub-optimal temperature. After a 5-day acclimation period, the first trifoliate leaves were sprayed with 10(-7) M LCO (NodBj-V (C(18:1), MeFuc)) purified from genistein-induced Bradyrhizobium japonicum culture, and harvested at 0 and 48 h following treatment. Microarray analysis was performed using Affymetrix GeneChip® Soybean Genome Arrays. Compared to the control at 48 h after LCO treatment, a total of 147 genes were differentially expressed as a result of LCO treatment, including a number of stress-related genes and transcription factors. In addition, during the 48 h time period following foliar spray application, over a thousand genes exhibited differential expression, including hundreds of those specific to the LCO-treated plants. Our results indicated that the dynamic soybean foliar transcriptome was highly responsive to LCO treatment. Quantitative real-time PCR (qPCR) validated the microarray data.

xGit2 and xRhoGAP 11A regulate convergent extension and tissue separation in Xenopus gastrulation

In a microarray-based screen for genes that are involved in tissue separation downstream of Paraxial Protocadherin (PAPC) and Frizzled-7 (Fz7)-mediated signaling we identified xGit2 and xRhoGAP 11A, two GTPase-activating proteins (GAP) for small GTPases. xGit2 and xRhoGAP 11A are expressed in the dorsal ectoderm, and their transcription is downregulated in the involuting dorsal mesoderm by PAPC and Fz7. Overexpression of xGit2 and xRhoGAP 11A inhibits Rho activity and impairs convergent extension movements as well as tissue separation behaviour. We propose that Rho activity in the involuting mesoderm is enhanced through inhibition of xGit2 and xRhoGAP 11A transcription by PAPC and Fz7. By this mechanism xRhoGAP 11A and xGit2 are restricted to the dorsal ectoderm, while Rho signaling is inhibited.

Identification of embryonic pancreatic genes using Xenopus DNA microarrays

The pancreas is both an exocrine and endocrine endodermal organ involved in digestion and glucose homeostasis. During embryogenesis, the anlagen of the pancreas arise from dorsal and ventral evaginations of the foregut that later fuse to form a single organ. To better understand the molecular genetics of early pancreas development, we sought to isolate markers that are uniquely expressed in this tissue. Microarray analysis was performed comparing dissected pancreatic buds, liver buds, and the stomach region of tadpole stage Xenopus embryos. A total of 912 genes were found to be differentially expressed between these organs during early stages of organogenesis. K-means clustering analysis predicted 120 of these genes to be specifically enriched in the pancreas. Of these, we report on the novel expression patterns of 24 genes. Our analyses implicate the involvement of previously unsuspected signaling pathways during early pancreas development. Developmental Dynamics 238:1455-1466, 2009. (c) 2009 Wiley-Liss, Inc.

Identification of novel transcripts with differential dorso-ventral expression in Xenopus gastrula using serial analysis of gene expression

Comparison of dorsal and ventral transcriptomes of Xenopus tropicalis gastrulae using serial analysis of gene expression provides at least 86 novel differentially expressed transcripts.

Background

Recent evidence from global studies of gene expression indicates that transcriptomes are more complex than expected. Xenopus has been typically used as a model organism to study early embryonic development, particularly dorso-ventral patterning. In order to identify novel transcripts involved in dorso-ventral patterning, we compared dorsal and ventral transcriptomes of Xenopus tropicalis at the gastrula stage using serial analysis of gene expression (SAGE).

Results

Of the experimental tags, 54.5% were confidently mapped to transcripts and 125 showed a significant difference in their frequency of occurrence between dorsal and ventral libraries. We selected 20 differentially expressed tags and assigned them to specific transcripts using bioinformatics and reverse SAGE. Five mapped to transcripts with known dorso-ventral expression and the frequency of appearance for these tags in each library is in agreement with the expression described by other methods. The other 15 tags mapped to transcripts with no previously described asymmetric expression along the dorso-ventral axis. The differential expression of ten of these novel transcripts was validated by in situ hybridization and/or RT-PCR. We can estimate that this SAGE experiment provides a list of at least 86 novel transcripts with differential expression along the dorso-ventral axis. Interestingly, the expression of some novel transcripts was independent of β-catenin.

Conclusions

Our SAGE analysis provides a list of novel transcripts with differential expression in the dorso-ventral axis and a large number of orphan tags that can be used to identify novel transcripts and to improve the current annotation of the X. tropicalis genome.

Identification of asymmetrically localized transcripts along the animal-vegetal axis of the Xenopus egg

In many organisms, proper embryo development depends on the asymmetrical distribution of mRNA in the cytoplasm of the egg. Here we report comprehensive screening of RNA localized in the animal or vegetal hemisphere of the Xenopus egg. Macroarrays including over 40,000 independent embryonic cDNA clones, representing at least 17,000 unigenes, were differentially hybridized with labeled probes synthesized from the mRNA of animal or vegetal blastomeres. After two rounds of screening, we identified 33 clones of transcripts that may be preferentially distributed in the vegetal region of the early stage embryo, but transcripts localized in the animal region were not found. To assess the array results, we performed northern blot and quantitative real-time reverse transcription-polymerase chain reaction analysis. As a result, 21 transcripts of the 33 were confirmed to be localized in the vegetal region of the early stage embryo. Whole-mount in situ hybridization analysis revealed that 11 transcripts, including 7 previously reported genes, were localized in the vegetal hemisphere of the egg. These 11 transcripts were categorized into three groups according to their expression patterns in the egg. The first group, which contained four transcripts, showed uniform expression in the vegetal hemisphere, similar to VegT. The second group, which contained three transcripts, showed gradual expression from the vegetal pole to the equator, similar to Vg1. The last group, which contained three transcripts, was expressed at the germ plasm, similar to Xdazl. One transcript, Xwnt11, showed both the second and the third expression patterns.

TGF-beta signaling-mediated morphogenesis: modulation of cell adhesion via cadherin endocytosis

The molecular mechanisms governing the cell behaviors underlying morphogenesis remain a major focus of research in both developmental biology and cancer biology. TGF-beta ligands control cell fate specification via Smad-mediated signaling. However, their ability to guide cellular morphogenesis in a variety of biological contexts is poorly understood. We report on the discovery of a novel TGF-beta signaling-mediated cellular morphogenesis occurring during vertebrate gastrulation. Activin/nodal members of the TGF-beta superfamily induce the expression of two genes regulating cell adhesion during gastrulation: Fibronectin Leucine-rich Repeat Transmembrane 3 (FLRT3), a type I transmembrane protein containing extracellular leucine-rich repeats, and the small GTPase Rnd1. FLRT3 and Rnd1 interact physically and modulate cell adhesion during embryogenesis by controlling cell surface levels of cadherin through a dynamin-dependent endocytosis pathway. Our model suggests that cell adhesion can be dynamically regulated by sequestering cadherin through internalization, and subsequent redeploying internalized cadherin to the cell surface as needed. As numerous studies have linked aberrant expression of small GTPases, adhesion molecules such as cadherins, and TGF-beta signaling to oncogenesis and metastasis, it is tempting to speculate that this FLRT3/Rnd1/cadherin pathway might also control cell behavior and morphogenesis in adult tissue homeostasis.

Xenopus cDNA microarray identification of genes with endodermal organ expression

The endoderm is classically defined as the innermost layer of three Metazoan germ layers. During organogenesis, the endoderm gives rise to the digestive and respiratory tracts as well as associated organs such as the liver, pancreas, and lung. At present, however, how the endoderm forms the variety of cell types of digestive and respiratory tracts as well as the budding organs is not well understood. In order to investigate the molecular basis and mechanism of organogenesis and to identify the endodermal organ-related marker genes, we carried out microarray analysis using Xenopus cDNA chips. To achieve this goal, we isolated the Xenopus gut endoderm from three different stages of Xenopus organogenesis, and separated each stage of gut endoderm into anterior and posterior regions. Competitive hybridization of cDNA between the anterior and posterior endoderm regions, to screen genes that specifically expressed in the major organs, revealed 915 candidates. We then selected 104 clones for in situ hybridization analysis. Here, we report the identification and expression patterns of the 104 Xenopus endodermal genes, which would serve as useful markers for studying endodermal organ development.

The Xfeb gene is directly upregulated by Zic1 during early neural development

The transcription factor Zic1 plays important roles in patterning the neural plate in early vertebrate development. However, few genes that are regulated by Zic1 are known. We have identified a new direct downstream target gene of Zic1 that we have named Xfeb. Xfeb is a member of the pathogenesis-related (PR) protein superfamily and contains five tandem SCP domains. The sequence of Xfeb suggests that it may possess serine protease activity. Xfeb is expressed in the presumptive hindbrain region during neurula stages and in somite tissues later in development. Xfeb represses the hindbrain gene hoxB1 and the anterior neural gene otx2, suggesting that Xfeb is involved in regionalizing the neural plate, possibly by ensuring a posterior expression limit for otx2.

Frizzled7 mediates canonical Wnt signaling in neural crest induction

The neural crest is a multipotent cell population that migrates from the dorsal edge of the neural tube to various parts of the embryo where it differentiates into a remarkable variety of different cell types. Initial induction of neural crest is mediated by a combination of BMP, Wnt, FGF, Retinoic acid and Notch/Delta signaling. The two-signal model for neural crest induction suggests that BMP signaling induces the competence to become neural crest. The second signal involves Wnt acting through the canonical pathway and leads to expression of neural crest markers such as slug. Wnt signals from the neural plate, non-neural ectoderm and paraxial mesoderm have all been suggested to play a role in neural crest induction. We show that Xenopus frizzled7 (Xfz7) is expressed in the dorsal ectoderm including early neural crest progenitors and is a key mediator of the Wnt inductive signal. We demonstrate that Xfz7 expression is induced in response to a BMP antagonist, noggin, and that Xfz7 can induce neural crest specific genes in noggin-treated ectodermal explants (animal caps). Morpholino-mediated or dominant negative inhibition of Xfz7 inhibits Wnt induced Xslug expression in the animal cap assay and in the whole embryo leading to a loss of neural crest derived pigment cells. Full-length Xfz7 rescues the morpholino-induced phenotype, as does activated beta-catenin, suggesting that Xfz7 is signaling through the canonical pathway. We therefore demonstrate that Xfz7 is regulated by BMP antagonism and is required for neural crest induction by Wnt in the developing vertebrate embryo.

Genomic analysis of Xenopus organizer function

Background

Studies of the Xenopus organizer have laid the foundation for our understanding of the conserved signaling pathways that pattern vertebrate embryos during gastrulation. The two primary activities of the organizer, BMP and Wnt inhibition, can regulate a spectrum of genes that pattern essentially all aspects of the embryo during gastrulation. As our knowledge of organizer signaling grows, it is imperative that we begin knitting together our gene-level knowledge into genome-level signaling models. The goal of this paper was to identify complete lists of genes regulated by different aspects of organizer signaling, thereby providing a deeper understanding of the genomic mechanisms that underlie these complex and fundamental signaling events.

Results

To this end, we ectopically overexpress Noggin and Dkk-1, inhibitors of the BMP and Wnt pathways, respectively, within ventral tissues. After isolating embryonic ventral halves at early and late gastrulation, we analyze the transcriptional response to these molecules within the generated ectopic organizers using oligonucleotide microarrays. An efficient statistical analysis scheme, combined with a new Gene Ontology biological process annotation of the Xenopus genome, allows reliable and faithful clustering of molecules based upon their roles during gastrulation. From this data, we identify new organizer-related expression patterns for 19 genes. Moreover, our data sub-divides organizer genes into separate head and trunk organizing groups, which each show distinct responses to Noggin and Dkk-1 activity during gastrulation.

Conclusion

Our data provides a genomic view of the cohorts of genes that respond to Noggin and Dkk-1 activity, allowing us to separate the role of each in organizer function. These patterns demonstrate a model where BMP inhibition plays a largely inductive role during early developmental stages, thereby initiating the suites of genes needed to pattern dorsal tissues. Meanwhile, Wnt inhibition acts later during gastrulation, and is essential for maintenance of organizer gene expression throughout gastrulation, a role which may depend on its ability to block the expression of a host of ventral, posterior, and lateral fate-specifying factors.

Negative regulation of Hedgehog signaling by the cholesterogenic enzyme 7-dehydrocholesterol reductase

Cholesterol regulates Hedgehog (Hh) signaling during early vertebrate development. Smith-Lemli-Opitz syndrome (SLOS) is caused by defects in 7-dehydrocholesterol reductase (DHCR7), an enzyme catalyzing the final step of cholesterol biosynthesis. Many developmental malformations attributed to SLOS occur in tissues and organs where Hh signaling is required for development, but the precise role of DHCR7 deficiency in this disease remains murky. We report that DHCR7 and Sonic Hedgehog (Shh) are co-expressed during midline development in Xenopus embryos. DHCR7 has previously been implicated to function as a positive regulator of Hh signaling that acts to regulate the cholesterol adduction of Hh ligand or to affect Hh signaling in the responding cell. We present gain- and loss-of-function analyses suggesting that DHCR7 functions as a negative regulator of Hh signaling at the level or downstream of Smoothened (Smo) and affects intracellular Hh signaling. Our analysis also raises the possibility that the human condition SLOS is caused not only by disruption of the enzymatic role of DHCR7 as a reductase in cholesterol biosynthesis, but may also involve defects in DHCR7 resulting in derepression of Shh signaling.

Grainyhead-like 3, a transcription factor identified in a microarray screen, promotes the specification of the superficial layer of the embryonic epidermis

The Xenopus ectoderm consists of two populations of cells, superficial polarised epithelial cells and deep, non-epithelial cells. These two cell types differ in their developmental fate. In the neural ectoderm, primary neurons are derived only from the deep cells. In the epidermal ectoderm, superficial cells express high levels of differentiation markers, while most of the deep cells do not differentiate until later when they produce the stratified adult epidermis. However, few molecular differences are known between the deep and superficial cells. Here, we have undertaken a systematic approach to identify genes that show layer-restricted expression by microarray analysis of deep and superficial cells at the gastrula stage, followed by wholemount in situ hybridisation. We have identified 32 differentially expressed genes, of which 26 show higher expression in the superficial layer and 6 in the deep layer and describe their expression at the gastrula and neurula stage. One of the identified genes is the transcription factor Grhl3, which we found to be expressed in the superficial layer of the gastrula ectoderm and the neurula epidermis. By using markers identified in this work, we show that Grlh3 promotes superficial gene expression in the deep layer of the epidermis. Concomitantly, deep layer specific genes are switched off, showing that Grlh3 can promote deep cells to take on a superficial cell identity in the embryonic epidermis.

An obligatory caravanserai stop on the silk road to neural induction: Inhibition of BMP/GDF signaling

Work in Xenopus laevis produced the first molecular explanation for neural specification, the default model, where inactivation of the BMP pathway in ectodermal cells changes fates from epidermal to neural. This review covers the present status of our understanding of neural specification, with emphasis on Xenopus, but including relevant facts in other model systems. While recent experiments have increased the complexity of the molecular picture, they have also provided additional support for the default model and the central position of the BMP pathway. We conclude that synergy between accumulated knowledge and technical progress will maintain Xenopus at the forefront of research in neural development.

Genetic and genomic prospects for Xenopus tropicalis research

Research using Xenopus laevis has made enormous contributions to our understanding of vertebrate development, control of the eukaryotic cell cycle and the cytoskeleton. One limitation, however, has been the lack of systematic genetic studies in Xenopus to complement molecular and cell biological investigations. Work with the closely related diploid frog Xenopus tropicalis is beginning to address this limitation. Here, we review the resources that will make genetic studies using X. tropicalis a reality.

Macroarray-based analysis of tail regeneration in Xenopus laevis larvae

Xenopus larvae possess a remarkable ability to regenerate their tails after they have been severed. To gain an understanding of the molecular mechanisms underlying tail regeneration, we performed a cDNA macroarray-based analysis of gene expression. A Xenopus cDNA macroarray representing 42,240 independent clones was differentially hybridized with probes synthesized from the total RNA of normal and regenerating tails. Temporal expression analysis revealed that the up-regulated genes could be grouped into early or late responding genes. A comparative expression analysis revealed that most genes showed similar expression patterns between tail development and regeneration. However, some genes showed regeneration-specific expression. Finally, we identified 48 up-regulated genes that fell into several categories based on their putative functions. These categories reflect the various processes that take place during regeneration, such as inflammation response, wound healing, cell proliferation, cell differentiation, and control of cell structure. Thus, we have identified a panel of genes that appear to be involved in the process of regeneration.

Xenomics

Xenopus genomics, or Xenomics for short, is coming of age. Indeed, biological insight into processes such as growth factor signaling and patterning of the early embryo is now being gained by combining the value of Xenopus as a model organism for cell and developmental biology with genomic approaches. In this review I address these recent advances and explore future possibilities gained from combining this powerful experimental system with genomic approaches, as well as how our quest to understand basic biological principles will be greatly facilitated though the marriage of Xenopus and genomics.

Regulatory targets for transcription factor AP2 in Xenopus embryos

The transcription factor AP2 (TFAP2) has an important role in regulating gene expression in both epidermis and neural crest cells. In order to further characterize these functions we have used a hormone inducible TFAP2alpha fusion protein in a Xenopus animal cap assay to identify downstream targets of this factor. The most common pattern comprised genes predominantly expressed in the epidermis. A second group was expressed at high levels in the neural crest, but all of these were also expressed in the epidermis as well as in other tissues in which TFAP2alpha has not been detected, suggesting modular control involving both TFAP2-dependent and TFAP2-independent components. In addition, a few strongly induced genes did not overlap at all in expression pattern with that of TFAP2alpha in the early embryo, and were also activated precociously in the experimentally manipulated ectoderm, and thus likely represent inappropriate regulatory interactions. A final group was identified that were repressed by TFAP2alpha and were expressed in the neural plate. These results provide further support for the importance of TFAP2alpha in ectoderm development, and also highlight the molecular linkage between the epidermis and neural crest in the Xenopus embryo.

XenDB: full length cDNA prediction and cross species mapping in Xenopus laevis

Background

Research using the model system Xenopus laevis has provided critical insights into the mechanisms of early vertebrate development and cell biology. Large scale sequencing efforts have provided an increasingly important resource for researchers. To provide full advantage of the available sequence, we have analyzed 350,468 Xenopus laevis Expressed Sequence Tags (ESTs) both to identify full length protein encoding sequences and to develop a unique database system to support comparative approaches between X. laevis and other model systems.

Description

Using a suffix array based clustering approach, we have identified 25,971 clusters and 40,877 singleton sequences. Generation of a consensus sequence for each cluster resulted in 31,353 tentative contig and 4,801 singleton sequences. Using both BLASTX and FASTY comparison to five model organisms and the NR protein database, more than 15,000 sequences are predicted to encode full length proteins and these have been matched to publicly available IMAGE clones when available. Each sequence has been compared to the KOG database and ~67% of the sequences have been assigned a putative functional category. Based on sequence homology to mouse and human, putative GO annotations have been determined.

Conclusion

The results of the analysis have been stored in a publicly available database XenDB . A unique capability of the database is the ability to batch upload cross species queries to identify potential Xenopus homologues and their associated full length clones. Examples are provided including mapping of microarray results and application of 'in silico' analysis. The ability to quickly translate the results of various species into 'Xenopus-centric' information should greatly enhance comparative embryological approaches.

Supplementary material can be found at .

Global analysis of RAR-responsive genes in the Xenopus neurula using cDNA microarrays

Retinoid signaling is important for patterning the vertebrate hindbrain and midaxial regions. We recently showed that signaling through retinoic acid receptors (RARs) is essential for anteroposterior patterning along the entire body axis. To further investigate the mechanisms through which RARs act, we used microarray analysis to investigate the effects of modulating RAR activity on target gene expression. We identified 334 up-regulated genes (92% of which were validated), including known RA-responsive genes, known genes that have never been proposed as RA targets and many hypothetical and unidentified genes (n = 166). Sixty-seven validated down-regulated genes were identified, including known RA-responsive genes and anterior marker genes. The expression patterns of selected up-regulated genes (n = 45) were examined at neurula stages using whole-mount in situ hybridization. We found that most of these genes were expressed in the neural tube and many were expressed in anterior tissues such as neural crest, brain, eye anlagen, and cement gland. Some were expressed in tissues such as notochord, somites, pronephros, and blood islands, where retinoic acid (RA) plays established roles in organogenesis. Members of this set of newly identified RAR target genes are likely to play important roles in neural patterning and organogenesis under the control of RAR signaling pathways, and their further characterization will expand our understanding of RA signaling during development.

A Xenopus DNA microarray approach to identify novel direct BMP target genes involved in early embryonic development

Bone morphogenetic proteins (BMPs), a subgroup of the transforming growth factor-beta (TGF-beta) superfamily, were originally isolated from bone on the basis of their ability to induce ectopic bone development. Although BMPs are involved in a wide range of developmental and physiological functions, very few vertebrate target genes in this pathway have been identified. To identify target genes regulated by the BMP growth factor family in Xenopus, large-scale microarray analyses were conducted to discover genes directly activated by this factor in dissociated animal cap tissues treated with a combination of the protein synthesis inhibitor cycloheximide and BMP2. Consequent expression patterns and behaviors of the most highly induced genes were analyzed by in situ and reverse transcriptase-polymerase chain reaction analyses. Here, we describe two sets of the most highly induced direct BMP target genes identified using microarrays prepared from two different stages of early Xenopus development. A wide variety of genes are induced by BMP2, ranging from cell cycle controllers, enzymes, signal transduction cascade components, and components of the blood and vascular system. The finding reinforces the notion that BMP signals play important roles in a variety of biological processes.

Xenopus as a model system to study transcriptional regulatory networks

Development is controlled by a complex series of events requiring sequential gene activation. Understanding the logic of gene networks during development is necessary for a complete understanding of how genes contribute to phenotype. Pioneering work initiated in the sea urchin and Drosophila has demonstrated that reasonable transcriptional regulatory network diagrams representing early development in multicellular animals can be generated through use of appropriate genomic, genetic, and biochemical tools. Establishment of similar regulatory network diagrams for vertebrate development is a necessary step. The amphibian Xenopus has long been used as a model for vertebrate early development and has contributed greatly to the elucidation of gene regulation. Because the best and most extensively studied transcriptional regulatory network in Xenopus is that underlying the formation and function of Spemann's organizer, we describe the current status of our understanding of this gene regulatory network and its relationship to mesodermal patterning. Seventy-four transcription factors currently known to be expressed in the mesoendoderm of Xenopus gastrula were characterized according to their modes of action, DNA binding consensus sequences, and target genes. Among them, nineteen transcription factors were characterized sufficiently in detail, allowing us to generate a gene regulatory network diagram. Additionally, we discuss recent amphibian work using a combined DNA microarray and bioinformatics approach that promises to accelerate regulatory network studies.