Multiplex cis-regulatory analysis.

Recent progress in multiplex cis-regulatory analysis has increased the speed of identifying enhancers and promoters, and enabled efficient incorporation of cis-regulatory information into gene regulatory network models. Three types of barcode reporters have been developed for multiplex reporter assays in sea urchin embryos: 13-tags and 129-tags for QPCR, 130 Nanotags for NanoString, and 100 million N25-tags for next-generation sequencing. In this chapter, to facilitate adoption of high-throughput cis-regulatory analysis in sea urchin embryos, I provide practical guidelines to best utilize barcode reporters that are compatible with either QPCR or next-generation sequencing. I expect that the guidelines are also applicable to other invertebrate embryos.

Single embryo-resolution quantitative analysis of reporters permits multiplex spatial cis-regulatory analysis.

Cis-regulatory modules (CRMs) control spatiotemporal gene expression patterns in embryos. While measurement of quantitative CRM activities has become efficient, measurement of spatial CRM activities still relies on slow, one-by-one methods. To overcome this bottleneck, we have developed a high-throughput method that can simultaneously measure quantitative and spatial CRM activities. The new method builds profiles of quantitative CRM activities measured at single-embryo resolution in many mosaic embryos and uses these profiles as a 'fingerprint' of spatial patterns. We show in sea urchin embryos that the new method, Multiplex and Mosaic Observation of Spatial Information encoded in Cis-regulatory modules (MMOSAIC), can efficiently predict spatial activities of new CRMs and can detect spatial responses of CRMs to gene perturbations. We anticipate that MMOSAIC will facilitate systems-wide functional analyses of CRMs in embryos.

Barcoded DNA-tag reporters for multiplex cis-regulatory analysis.

Cis-regulatory DNA sequences causally mediate patterns of gene expression, but efficient experimental analysis of these control systems has remained challenging. Here we develop a new version of "barcoded" DNA-tag reporters, "Nanotags" that permit simultaneous quantitative analysis of up to 130 distinct cis-regulatory modules (CRMs). The activities of these reporters are measured in single experiments by the NanoString RNA counting method and other quantitative procedures. We demonstrate the efficiency of the Nanotag method by simultaneously measuring hourly temporal activities of 126 CRMs from 46 genes in the developing sea urchin embryo, otherwise a virtually impossible task. Nanotags are also used in gene perturbation experiments to reveal cis-regulatory responses of many CRMs at once. Nanotag methodology can be applied to many research areas, ranging from gene regulatory networks to functional and evolutionary genomics.

High accuracy, high-resolution prevalence measurement for the majority of locally expressed regulatory genes in early sea urchin development.

Accurate measurements of transcript abundance are a prerequisite to understand gene activity in development. Using the NanoString nCounter, an RNA counting device, we measured the prevalence of 172 transcription factors and signaling molecules in early sea urchin development. These measurements show high fidelity over more than five orders of magnitude down to a few transcripts per embryo. Most of the genes included are locally restricted in their spatial expression, and contribute to the divergent regulatory states of cells in the developing embryo. In order to obtain high-resolution expression profiles from fertilization to late gastrulation samples were collected at hourly intervals. The measured time courses agree well with, and substantially extend, prior relative abundance measurements obtained by quantitative PCR. High temporal resolution permits sequences of successively activated genes to be precisely delineated providing an ancillary tool for assembling maps of gene regulatory networks. The data are available via an interactive website for quick plotting of selected time courses.

Functional cis-regulatory genomics for systems biology.

Gene expression is controlled by interactions between trans-regulatory factors and cis-regulatory DNA sequences, and these interactions constitute the essential functional linkages of gene regulatory networks (GRNs). Validation of GRN models requires experimental cis-regulatory tests of predicted linkages to authenticate their identities and proposed functions. However, cis-regulatory analysis is, at present, at a severe bottleneck in genomic system biology because of the demanding experimental methodologies currently in use for discovering cis-regulatory modules (CRMs), in the genome, and for measuring their activities. Here we demonstrate a high-throughput approach to both discovery and quantitative characterization of CRMs. The unique aspect is use of DNA sequence tags to "barcode" CRM expression constructs, which can then be mixed, injected together into sea urchin eggs, and subsequently deconvolved. This method has increased the rate of cis-regulatory analysis by >100-fold compared with conventional one-by-one reporter assays. The utility of the DNA-tag reporters was demonstrated by the rapid discovery of 81 active CRMs from 37 previously unexplored sea urchin genes. We then obtained simultaneous high-resolution temporal characterization of the regulatory activities of more than 80 CRMs. On average 2-3 CRMs were discovered per gene. Comparison of endogenous gene expression profiles with those of the CRMs recovered from each gene showed that, for most cases, at least one CRM is active in each phase of endogenous expression, suggesting that CRM recovery was comprehensive. This approach will qualitatively alter the practice of GRN construction as well as validation, and will impact many additional areas of regulatory system biology.

Exclusive developmental functions of gatae cis-regulatory modules in the Strongylocentrorus purpuratus embryo.

The gatae gene of Strongylocentrotus purpuratus is orthologous to vertebrate gata-4,5,6 genes. This gene is expressed in the endomesoderm in the blastula and later the gut of the embryo, and is required for normal development. A gatae BAC containing a GFP reporter knocked into exon one of the gene was able to reproduce all aspects of endogenous gatae expression in the embryo. To identify putative gatae cis-regulatory modules we carried out an interspecific sequence conservation analysis with respect to a Lytechinus variegatus gatae BAC, which revealed 25 conserved non-coding sequence patches. These were individually tested in gene transfer experiments, and two modules capable of driving localized reporter expression in the embryo were identified. Module 10 produces early expression in mesoderm and endoderm cells up to the early gastrula stage, while module 24 generates late endodermal expression at gastrula and pluteus stages. Module 10 was then deleted from the gatae BAC by reciprocal recombination, resulting in total loss of reporter expression in the time frame in which it is normally active. Similar deletion of module 24 led to ubiquitous GFP expression in the gastrula and pluteus. These results show that Module 10 is uniquely necessary and sufficient to account for the early phase of gatae expression during endomesoderm specification. In addition, they imply a functional cis-regulatory module exclusion, whereby only a single module can associate with the basal promoter and drive gene expression at any given time.

Cis-regulatory control of the nodal gene, initiator of the sea urchin oral ectoderm gene network.

Expression of the nodal gene initiates the gene regulatory network which establishes the transcriptional specification of the oral ectoderm in the sea urchin embryo. This gene encodes a TGFbeta ligand, and in Strongylocentrotus purpuratus its transcription is activated in the presumptive oral ectoderm at about the 30-cell stage. Thereafter Nodal signaling occurs among all cells of the oral ectoderm territory, and nodal expression is required for expression of oral ectoderm regulatory genes. The cis-regulatory system of the nodal gene transduces anisotropically distributed cytoplasmic cues that distinguish the future oral and aboral domains of the early embryo. Here we establish the genomic basis for the initiation and maintenance of nodal gene expression in the oral ectoderm. Functional cis-regulatory control modules of the nodal gene were identified by interspecific sequence conservation. A 5' cis-regulatory module functions both to initiate expression of the nodal gene and to maintain its expression by means of feedback input from the Nodal signal transduction system. These functions are mediated respectively by target sites for bZIP transcription factors, and by SMAD target sites. At least one SMAD site is also needed for the initiation of expression. An intron module also contains SMAD sites which respond to Nodal feedback, and in addition acts to repress vegetal expression. These observations explain the main features of nodal expression in the oral ectoderm: since the activity of bZIP factors is redox sensitive, and the initial polarization of oral vs. aboral fate is manifested in a redox differential, the bZIP sites account for the activation of nodal on the oral side; and since the immediate early signal transduction response factors for Nodal are SMAD factors, the SMAD sites account for the feedback maintenance of nodal gene expression.

Curated genome annotation of Oryza sativa ssp. japonica and comparative genome analysis with Arabidopsis thaliana.

We present here the annotation of the complete genome of rice Oryza sativa L. ssp. japonica cultivar Nipponbare. All functional annotations for proteins and non-protein-coding RNA (npRNA) candidates were manually curated. Functions were identified or inferred in 19,969 (70%) of the proteins, and 131 possible npRNAs (including 58 antisense transcripts) were found. Almost 5000 annotated protein-coding genes were found to be disrupted in insertional mutant lines, which will accelerate future experimental validation of the annotations. The rice loci were determined by using cDNA sequences obtained from rice and other representative cereals. Our conservative estimate based on these loci and an extrapolation suggested that the gene number of rice is approximately 32,000, which is smaller than previous estimates. We conducted comparative analyses between rice and Arabidopsis thaliana and found that both genomes possessed several lineage-specific genes, which might account for the observed differences between these species, while they had similar sets of predicted functional domains among the protein sequences. A system to control translational efficiency seems to be conserved across large evolutionary distances. Moreover, the evolutionary process of protein-coding genes was examined. Our results suggest that natural selection may have played a role for duplicated genes in both species, so that duplication was suppressed or favored in a manner that depended on the function of a gene.

Evolutionary change of the numbers of homeobox genes in bilateral animals.

It has been known that the conservation or diversity of homeobox genes is responsible for the similarity and variability of some of the morphological or physiological characters among different organisms. To gain some insights into the evolutionary pattern of homeobox genes in bilateral animals, we studied the change of the numbers of these genes during the evolution of bilateral animals. We analyzed 2,031 homeodomain sequences compiled from 11 species of bilateral animals ranging from Caenorhabditis elegans to humans. Our phylogenetic analysis using a modified reconciled-tree method suggested that there were at least about 88 homeobox genes in the common ancestor of bilateral animals. About 50-60 genes of them have left at least one descendant gene in each of the 11 species studied, suggesting that about 30-40 genes were lost in a lineage-specific manner. Although similar numbers of ancestral genes have survived in each species, vertebrate lineages gained many more genes by duplication than invertebrate lineages, resulting in more than 200 homeobox genes in vertebrates and about 100 in invertebrates. After these gene duplications, a substantial number of old duplicate genes have also been lost in each lineage. Because many old duplicate genes were lost, it is likely that lost genes had already been differentiated from other groups of genes at the time of gene loss. We conclude that both gain and loss of homeobox genes were important for the evolutionary change of phenotypic characters in bilateral animals.

A simple method for predicting the functional differentiation of duplicate genes and its application to MIKC-type MADS-box genes.

A simple statistical method for predicting the functional differentiation of duplicate genes was developed. This method is based on the premise that the extent of functional differentiation between duplicate genes is reflected in the difference in evolutionary rate because the functional change of genes is often caused by relaxation or intensification of functional constraints. With this idea in mind, we developed a window analysis of protein sequences to identify the protein regions in which the significant rate difference exists. We applied this method to MIKC-type MADS-box proteins that control flower development in plants. We examined 23 pairs of sequences of floral MADS-box proteins from petunia and found that the rate differences for 14 pairs are significant. The significant rate differences were observed mostly in the K domain, which is important for dimerization between MADS-box proteins. These results indicate that our statistical method may be useful for predicting protein regions that are likely to be functionally differentiated. These regions may be chosen for further experimental studies.

Type I MADS-box genes have experienced faster birth-and-death evolution than type II MADS-box genes in angiosperms.

Plant MADS-box genes form a large gene family for transcription factors and are involved in various aspects of developmental processes, including flower development. They are known to be subject to birth-and-death evolution, but the detailed features of this mode of evolution remain unclear. To have a deeper insight into the evolutionary pattern of this gene family, we enumerated all available functional and nonfunctional (pseudogene) MADS-box genes from the Arabidopsis and rice genomes. Plant MADS-box genes can be classified into types I and II genes on the basis of phylogenetic analysis. Conducting extensive homology search and phylogenetic analysis, we found 64 presumed functional and 37 nonfunctional type I genes and 43 presumed functional and 4 nonfunctional type II genes in Arabidopsis. We also found 24 presumed functional and 6 nonfunctional type I genes and 47 presumed functional and 1 nonfunctional type II genes in rice. Our phylogenetic analysis indicated there were at least about four to eight type I genes and approximately 15-20 type II genes in the most recent common ancestor of Arabidopsis and rice. It has also been suggested that type I genes have experienced a higher rate of birth-and-death evolution than type II genes in angiosperms. Furthermore, the higher rate of birth-and-death evolution in type I genes appeared partly due to a higher frequency of segmental gene duplication and weaker purifying selection in type I than in type II genes.

Generation and analysis of end sequence database for T-DNA tagging lines in rice.

We analyzed 6749 lines tagged by the gene trap vector pGA2707. This resulted in the isolation of 3793 genomic sequences flanking the T-DNA. Among the insertions, 1846 T-DNAs were integrated into genic regions, and 1864 were located in intergenic regions. Frequencies were also higher at the beginning and end of the coding regions and upstream near the ATG start codon. The overall GC content at the insertion sites was close to that measured from the entire rice (Oryza sativa) genome. Functional classification of these 1846 tagged genes showed a distribution similar to that observed for all the genes in the rice chromosomes. This indicates that T-DNA insertion is not biased toward a particular class of genes. There were 764, 327, and 346 T-DNA insertions in chromosomes 1, 4 and 10, respectively. Insertions were not evenly distributed; frequencies were higher at the ends of the chromosomes and lower near the centromere. At certain sites, the frequency was higher than in the surrounding regions. This sequence database will be valuable in identifying knockout mutants for elucidating gene function in rice. This resource is available to the scientific community at http://www.postech.ac.kr/life/pfg/risd.

Antiquity and evolution of the MADS-box gene family controlling flower development in plants.

MADS-box genes in plants control various aspects of development and reproductive processes including flower formation. To obtain some insight into the roles of these genes in morphological evolution, we investigated the origin and diversification of floral MADS-box genes by conducting molecular evolutionary genetics analyses. Our results suggest that the most recent common ancestor of today's floral MADS-box genes evolved roughly 650 MYA, much earlier than the Cambrian explosion. They also suggest that the functional classes T (SVP), B (and Bs), C, F (AGL20 or TM3), A, and G (AGL6) of floral MADS-box genes diverged sequentially in this order from the class E gene lineage. The divergence between the class G and E genes apparently occurred around the time of the angiosperm/gymnosperm split. Furthermore, the ancestors of three classes of genes (class T genes, class B/Bs genes, and the common ancestor of the other classes of genes) might have existed at the time of the Cambrian explosion. We also conducted a phylogenetic analysis of MADS-domain sequences from various species of plants and animals and presented a hypothetical scenario of the evolution of MADS-box genes in plants and animals, taking into account paleontological information. Our study supports the idea that there are two main evolutionary lineages (type I and type II) of MADS-box genes in plants and animals.

Systematic reverse genetic screening of T-DNA tagged genes in rice for functional genomic analyses: MADS-box genes as a test case.

We have generated 47 DNA pools and 235 subpools from 21,049 T-DNA insertion lines of rice. DNA pools of 500-1,000 lines were adequate for screening a T-DNA insertion within a 2-kb region. To examine the efficacy of the DNA pools, we selected MADS-box genes, which play an important role in controlling various aspects of plant development. A total of 34 MIKC-type MADS-box genes have now been identified from rice sequence databases. Our PCR screening for T-DNA insertions within 12 MADS-box genes resulted in the identification of five insertions in four different genes. These DNA pools will be valuable when isolating T-DNA insertional mutants in various rice genes. The DNA pool screening service and the mutant seeds are available upon request to genean@postech.ac.kr.