Characterization of the pufferfish Otx2 cis-regulators reveals evolutionarily conserved genetic mechanisms for vertebrate head specification

Modification of gene expression and soluble factor secretion in the lateral ventricle choroid plexus: Analysis of the impacts on the neocortical development

The choroid plexus (ChP) is the center of soluble factor secretion into the cerebrospinal fluid in the central nervous system. It is known that various signaling factors secreted from the ChP are involved in the regulation of brain development and homeostasis. Intriguingly, the size of the ChP was prominently expanded in the brains of primates, including humans, suggesting that the expansion of the ChP contributed to mammalian brain evolution, leading to the acquisition of higher intelligence and cognitive functions. To address this hypothesis, we established transgenic (Tg) systems using regulatory elements that direct expression of candidate genes in the ChP. Overexpression of sonic hedgehog (Shh) in the developing ChP led to the expansion of the ChP with greater arborization. Shh produced in the ChP caused an increase in neural stem cells (NSCs) in the neocortical region, leading to the expansion of ventricles, ventricular zone, neocortical surface area, and neocortical surface folding. These findings suggest that the activation of Shh signaling via its enhanced secretion from the developing ChP contributed to the evolution of the neocortex. Furthermore, we found that Shh produced in the ChP enhanced NSC proliferation in the postnatal Tg brain, demonstrating that our Tg system can be used to estimate the effects of candidate factors secreted from the ChP on various aspects of brain morphogenesis and functions.

What Goes Wrong during Early Development of Artificially Reproduced European Eel Anguilla anguilla? Clues from the Larval Transcriptome and Gene Expression Patterns

Simple Summary

Closing the life cycle of the European eel in captivity is urgently needed to gain perspective for the commercial production of juvenile glass eels. Larvae are produced weekly at our facilities, but large variations in larval mortality are observed during the first week after hatching. Although much effort has been devoted to investigating ways to prevent early larval mortality, it remains unclear what the causes are. The aim of this study was to perform a transcriptomic study on European eel larvae in order to identify genes and physiological pathways that are differentially regulated in the comparison of larvae from batches that did not survive for longer than three days vs. larvae from batches that survived for at least a week up to 22 days after hatching (non-viable vs. viable larvae). In contrast to earlier published studies on European eel, we conclude that larvae exhibit immune competency. Non-viable larvae initiated an inflammatory and host protection immune response and tried to maintain osmoregulatory homeostasis. As a perspective, microbial control and salinity reduction might benefit eel larvae in terms of lower mortality and improved development by lowering the costs of immune functioning and osmoregulation.

Abstract

In eels, large variations in larval mortality exist, which would impede the viable production of juvenile glass eels in captivity. The transcriptome of European eel larvae was investigated to identify physiological pathways and genes that show differential regulation between non-viable vs. viable larvae. Expression of genes involved in inflammation and host protection was higher, suggesting that non-viable larvae suffered from microbial infection. Expression of genes involved in osmoregulation was also higher, implying that non-viable larvae tried to maintain homeostasis by strong osmoregulatory adaptation. Expression of genes involved in myogenesis, neural, and sensory development was reduced in the non-viable larvae. Expression of the major histocompatibility complex class-I (mhc1) gene, M-protein (myom2), the dopamine 2B receptor (d2br), the melatonin receptor (mtr1), and heat-shock protein beta-1 (hspb1) showed strong differential regulation and was therefore studied in 1, 8, and 15 days post-hatch (dph) larvae by RT-PCR to comprehend the roles of these genes during ontogeny. Expression patterning of these genes indicated the start of active swimming (8 dph) and feed searching behavior (15 dph) and confirmed immunocompetence immediately after hatching. This study revealed useful insights for improving larval survival by microbial control and salinity reduction.

A novel measure of non-coding genome conservation identifies genomic regulatory blocks within primates

MOTIVATION

Clusters of extremely conserved non-coding elements (CNEs) mark genomic regions devoted to cis-regulation of key developmental genes in Metazoa. We have recently shown that their span coincides with that of topologically associating domains (TADs), making them useful for estimating conserved TAD boundaries in the absence of Hi-C data. The standard approach-detecting CNEs in genome alignments and then establishing the boundaries of their clusters-requires tuning of several parameters and breaks down when comparing closely related genomes.

RESULTS

We present a novel, kurtosis-based measure of pairwise non-coding conservation that requires no pre-set thresholds for conservation level and length of CNEs. We show that it performs robustly across a large span of evolutionary distances, including across the closely related genomes of primates for which standard approaches fail. The method is straightforward to implement and enables detection and comparison of clusters of CNEs and estimation of underlying TADs across a vastly increased range of Metazoan genomes.

AVAILABILITY AND IMPLEMENTATION

The data generated for this study, and the scripts used to generate the data, can be found at https://github.com/alexander-nash/kurtosis_conservation.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Super-enhancers in transcriptional regulation and genome organization

Gene expression is precisely controlled in a stage and cell-type-specific manner, largely through the interaction between cis-regulatory elements and their associated trans-acting factors. Where these components aggregate in promoters and enhancers, they are able to cooperate to modulate chromatin structure and support the engagement in long-range 3D superstructures that shape the dynamics of a cell's genomic architecture. Recently, the term 'super-enhancer' has been introduced to describe a hyper-active regulatory domain comprising a complex array of sequence elements that work together to control the key gene networks involved in cell identity. Here, we survey the unique characteristics of super-enhancers compared to other enhancer types and summarize the recent advances in our understanding of their biological role in gene regulation. In particular, we discuss their capacity to attract the formation of phase-separated condensates, and capacity to generate three-dimensional genome structures that precisely activate their target genes. We also propose a multi-stage transition model to explain the evolutionary pressure driving the development of super-enhancers in complex organisms, and highlight the potential for involvement in tumorigenesis. Finally, we discuss more broadly the role of super-enhancers in human health disorders and related potential in therapeutic interventions.

Regulatory variants of FOXG1 in the context of its topological domain organisation

FOXG1 syndrome is caused by FOXG1 intragenic point mutations, or by long-range position effects (LRPE) of intergenic structural variants. However, the size of the FOXG1 regulatory landscape is uncertain, because the associated topologically associating domain (TAD) in fibroblasts is split into two domains in embryonic stem cells (hESC). Indeed, it has been suggested that the pathogenetic mechanism of deletions that remove the stem-cell-specific TAD boundary may be enhancer adoption due to ectopic activity of enhancer(s) located in the distal hESC-TAD. Herein we map three de novo translocation breakpoints to the proximal regulatory domain of FOXG1. The classical FOXG1 syndrome in these and in other translocation patients, and in a patient with an intergenic deletion that removes the hESC-specific TAD boundary, do not support the hypothesised enhancer adoption as a main contributor to the FOXG1 syndrome. Also, virtual 4 C and HiC-interaction data suggest that the hESC-specific TAD boundary may not be critical for FOXG1 regulation in a majority of human cells and tissues, including brain tissues and a neuronal progenitor cell line. Our data support the importance of a critical regulatory region (SRO) proximal to the hESC-specific TAD boundary. We further narrow this critical region by a deletion distal to the hESC-specific boundary, associated with a milder clinical phenotype. The distance from FOXG1 to the SRO ( > 500 kb) highlight a limitation of ENCODE DNase hypersensitivity data for functional prediction of LRPE. Moreover, the SRO has little overlap with a cluster of frequently associating regions (FIREs) located in the proximal hESC-TAD.

Enrichment of megabase-sized DNA molecules for single-molecule optical mapping and next-generation sequencing

Next-generation sequencing (NGS) has caused a revolution, yet left a gap: long-range genetic information from native, non-amplified DNA fragments is unavailable. It might be obtained by optical mapping of megabase-sized DNA molecules. Frequently only a specific genomic region is of interest, so here we introduce a method for selection and enrichment of megabase-sized DNA molecules intended for single-molecule optical mapping: DNA from a human cell line is digested by the NotI rare-cutting enzyme and size-selected by pulsed-field gel electrophoresis. For demonstration, more than 600 sub-megabase- to megabase-sized DNA molecules were recovered from the gel and analysed by denaturation-renaturation optical mapping. Size-selected molecules from the same gel were sequenced by NGS. The optically mapped molecules and the NGS reads showed enrichment from regions defined by NotI restriction sites. We demonstrate that the unannotated genome can be characterized in a locus-specific manner via molecules partially overlapping with the annotated genome. The method is a promising tool for investigation of structural variants in enriched human genomic regions for both research and diagnostic purposes. Our enrichment method could potentially work with other genomes or target specified regions by applying other genomic editing tools, such as the CRISPR/Cas9 system.

Topologically associating domains are ancient features that coincide with Metazoan clusters of extreme noncoding conservation

Developmental genes in metazoan genomes are surrounded by dense clusters of conserved noncoding elements (CNEs). CNEs exhibit unexplained extreme levels of sequence conservation, with many acting as developmental long-range enhancers. Clusters of CNEs define the span of regulatory inputs for many important developmental regulators and have been described previously as genomic regulatory blocks (GRBs). Their function and distribution around important regulatory genes raises the question of how they relate to 3D conformation of these loci. Here, we show that clusters of CNEs strongly coincide with topological organisation, predicting the boundaries of hundreds of topologically associating domains (TADs) in human and Drosophila. The set of TADs that are associated with high levels of noncoding conservation exhibit distinct properties compared to TADs devoid of extreme noncoding conservation. The close correspondence between extreme noncoding conservation and TADs suggests that these TADs are ancient, revealing a regulatory architecture conserved over hundreds of millions of years.

Metazoan genomes contain many clusters of conserved noncoding elements. Here, the authors provide evidence that these clusters coincide with distinct topologically associating domains in humans and Drosophila, revealing a conserved regulatory genomic architecture.

Identification and functional characterization of novel transcriptional enhancers involved in regulating human GLI3 expression during early development

The zinc-finger transcription factor GLI3 acts as a primary transducer of Sonic hedgehog (Shh) signaling in a context-dependent combinatorial fashion. GLI3 participates in the patterning and growth of many organs, including the central nervous system (CNS) and limbs. Previously, we reported a subset of human intronic cis-regulators controlling many known aspects of endogenous Gli3 expression in mouse and zebrafish. Here we demonstrate in a transgenic zebrafish assay the potential of two novel tetrapod-teleost conserved non-coding elements (CNEs) docking within GLI3 intronic intervals (intron 3 and 4) to induce reporter gene expression at known sites of endogenous Gli3 transcription in embryonic domains such as the central nervous system (CNS) and limbs. Interestingly, the cell culture based assays reveal harmony with the context dependent dual nature of intra-GLI3 conserved elements. Furthermore, a transgenic zebrafish assay of previously reported limb-specific GLI3 transcriptional enhancers (previously tested in mice and chicken limb buds) induced reporter gene expression in zebrafish blood precursor cells and notochord instead of fin. These results demonstrate that the appendage-specific activity of a subset of GLI3-associated enhancers might be a tetrapod innovation. Taken together with our recent data, these results suggest that during the course of vertebrate evolution Gli3 expression control acquired a complex cis-regulatory landscape for spatiotemporal patterning of CNS and limbs. Comparative data from fish and mice suggest that the functional aspects of a subset of these cis-regulators have diverged significantly between these two lineages.

Phylum-Level Conservation of Regulatory Information in Nematodes despite Extensive Non-coding Sequence Divergence

Gene regulatory information guides development and shapes the course of evolution. To test conservation of gene regulation within the phylum Nematoda, we compared the functions of putative cis-regulatory sequences of four sets of orthologs (unc-47, unc-25, mec-3 and elt-2) from distantly-related nematode species. These species, Caenorhabditis elegans, its congeneric C. briggsae, and three parasitic species Meloidogyne hapla, Brugia malayi, and Trichinella spiralis, represent four of the five major clades in the phylum Nematoda. Despite the great phylogenetic distances sampled and the extensive sequence divergence of nematode genomes, all but one of the regulatory elements we tested are able to drive at least a subset of the expected gene expression patterns. We show that functionally conserved cis-regulatory elements have no more extended sequence similarity to their C. elegans orthologs than would be expected by chance, but they do harbor motifs that are important for proper expression of the C. elegans genes. These motifs are too short to be distinguished from the background level of sequence similarity, and while identical in sequence they are not conserved in orientation or position. Functional tests reveal that some of these motifs contribute to proper expression. Our results suggest that conserved regulatory circuitry can persist despite considerable turnover within cis elements.

To explore the phylogenetic limits of conservation of cis-regulatory elements, we used transgenesis to test the functions of enhancers of four genes from several species spanning the phylum Nematoda. While we found a striking degree of functional conservation among the examined cis elements, their DNA sequences lacked apparent conservation with the C. elegans orthologs. In fact, sequence similarity between C. elegans and the distantly related nematodes was no greater than would be expected by chance. Short motifs, similar to known regulatory sequences in C. elegans, can be detected in most of the cis elements. When tested, some of these sites appear to mediate regulatory function. However, they seem to have originated through motif turnover, rather than to have been preserved from a common ancestor. Our results suggest that gene regulatory networks are broadly conserved in the phylum Nematoda, but this conservation persists despite substantial reorganization of regulatory elements and could not be detected using naïve comparisons of sequence similarity.

Identification and analysis of a novel bmp4 enhancer in Fugu genome

Spatiotemporal expression of bone morphogenetic protein 4 (Bmp4) in epithelial and mesenchymal cells is critical for the development of many organs including teeth. Since Bmp4 has a complex and widespread regulatory area in mammals, the tissue-specific enhancers that are responsible for mesenchymal expression of Bmp4 are difficult to identify in mammals. TakiFugu rubripes (Fugu, pufferfish) has a highly compact genome size and is widely used in comparative genomics studies of gene regulatory mechanisms. In this study, we used the Fugu genome to evaluate the 15kb promoter region upstream of the Fugu bmp4 gene. By DNA segmental cloning and luciferase assay with two dental odontoblast-like cell lines, a dental ameloblast-like cell line, and a kidney fibroblast cell line, we identified a 485bp cis-regulatory enhancer between -4213 and -3728bp of the Fugu bmp4 gene. This enhancer showed strong transcriptional activity in all three dental cell lines and, to a lesser extent, also in kidney fibroblast cells. Though not located in an evolutionary conserved region, the enhancer activity for the DNA segment is intense. This is the first time a bmp4 enhancer sequence with activity in both mesenchymal and epithelial cells has been identified, which will help to decode the mechanism of tooth development in vertebrates.

Identification and computational analysis of gene regulatory elements

Over the last two decades, advances in experimental and computational technologies have greatly facilitated genomic research. Next-generation sequencing technologies have made de novo sequencing of large genomes affordable, and powerful computational approaches have enabled accurate annotations of genomic DNA sequences. Charting functional regions in genomes must account for not only the coding sequences, but also noncoding RNAs, repetitive elements, chromatin states, epigenetic modifications, and gene regulatory elements. A mix of comparative genomics, high-throughput biological experiments, and machine learning approaches has played a major role in this truly global effort. Here we describe some of these approaches and provide an account of our current understanding of the complex landscape of the human genome. We also present overviews of different publicly available, large-scale experimental data sets and computational tools, which we hope will prove beneficial for researchers working with large and complex genomes.

Chromatin accessibility mapping identifies mediators of basal transcription and retinoid-induced repression of OTX2 in medulloblastoma

Despite an emerging understanding of the genetic alterations giving rise to various tumors, the mechanisms whereby most oncogenes are overexpressed remain unclear. Here we have utilized an integrated approach of genomewide regulatory element mapping via DNase-seq followed by conventional reporter assays and transcription factor binding site discovery to characterize the transcriptional regulation of the medulloblastoma oncogene Orthodenticle Homeobox 2 (OTX2). Through these studies we have revealed that OTX2 is differentially regulated in medulloblastoma at the level of chromatin accessibility, which is in part mediated by DNA methylation. In cell lines exhibiting chromatin accessibility of OTX2 regulatory regions, we found that autoregulation maintains OTX2 expression. Comparison of medulloblastoma regulatory elements with those of the developing brain reveals that these tumors engage a developmental regulatory program to drive OTX2 transcription. Finally, we have identified a transcriptional regulatory element mediating retinoid-induced OTX2 repression in these tumors. This work characterizes for the first time the mechanisms of OTX2 overexpression in medulloblastoma. Furthermore, this study establishes proof of principle for applying ENCODE datasets towards the characterization of upstream trans-acting factors mediating expression of individual genes.

Identification of duplication downstream of BMP2 in a Chinese family with brachydactyly type A2 (BDA2)

Brachydactyly type A2 (BDA2, MIM 112600) is characterized by the deviation and shortening of the middle phalange of the index finger and the second toe. Using genome-wide linkage analysis in a Chinese BDA2 family, we mapped the maximum candidate interval of BDA2 to a ∼1.5 Mb region between D20S194 and D20S115 within chromosome 20p12.3 and found that the pairwise logarithm of the odds score was highest for marker D20S156 (Zmax = 6.09 at θ = 0). Based on functional and positional perspectives, the bone morphogenetic protein 2 (BMP2) gene was identified as the causal gene for BDA2 in this region, even though no point mutation was detected in BMP2. Through further investigation, we identified a 4,671 bp (Chr20: 6,809,218-6,813,888) genomic duplication downstream of BMP2. This duplication was located within the linked region, co-segregated with the BDA2 phenotype in this family, and was not found in the unaffected family members and the unrelated control individuals. Compared with the previously reported duplications, the duplication in this family has a different breakpoint flanked by the microhomologous sequence GATCA and a slightly different length. Some other microhomologous nucleotides were also found in the duplicated region. In summary, our findings support the conclusions that BMP2 is the causing gene for BDA2, that the genomic location corresponding to the duplication region is prone to structural changes associated with malformation of the digits, and that this tendency is probably caused by the abundance of microhomologous sequences in the region.

Enhancers

Transcription of eukaryotic genes is an exceedingly sophisticated and complicated process, orchestrated by layers of control mechanisms involving a myriad of transcription factors and DNA control sequences, with both groups subject to multiple modifications. The availability of various recent genomic approaches has provided previously unforeseen opportunities to examine the cis-regulatory landscape of the entire genome, resulting in the identification of a potentially overwhelming number of enhancers and novel enhancer functions. In this review, we focus on the activities of enhancers in metazoans and discuss how they serve to regulate gene expression during early development.

A lineage specific enhancer drives Otx2 expression in teleost organizer tissues

In mouse Otx2 plays essential roles in anterior-posterior axis formation and head development in anterior visceral endoderm and anterior mesendoderm. The Otx2 expression in these sites is regulated by VE and CM enhancers at the 5' proximal to the translation start site, and we proposed that these enhancers would have been established in ancestral sarcoptergians after divergence from actinopterigians for the use of Otx2 as the head organizer gene (Kurokawa et al., 2010). This would make doubtful an earlier proposal of ours that a 1.1 kb fragment located at +14.4 to +15.5 kb 3' (3'En) of fugu Otx2a gene harbors enhancers phylogenetically and functionally homologous to mouse VE and CM enhancers (Kimura-Yoshida et al., 2007). In the present study, we demonstrate that fugu Otx2a is not expressed in the dorsal margin of blastoderm, shield and early anterior mesendoderm, and that the fugu Otx2a 3'En do not exhibit activities at these sites of fugu embryos. We conclude that the fugu Otx2a 3'En does not harbor an organizer enhancer, but encodes an enhancer for the expression in later anterior mesendodermal tissues. Instead, in fugu embryos Otx2b is expressed in the dorsal margin of blastoderm at blastula stage and shield at 50% epiboly, and this expression is directed by an enhancer, 5'En, located at -1000 to -800 bp, which is uniquely conserved among teleost Otx2b orthologues.

Controlled somatic and germline copy number variation in the mouse model

Changes in the number of chromosomes, but also variations in the copy number of chromosomal regions have been described in various pathological conditions, such as cancer and aneuploidy, but also in normal physiological condition. Our classical view of DNA replication and mitotic preservation of the chromosomal integrity is now challenged as new technologies allow us to observe such mosaic somatic changes in copy number affecting regions of chromosomes with various sizes. In order to go further in the understanding of copy number influence in normal condition we could take advantage of the novel strategy called Targeted Asymmetric Sister Chromatin Event of Recombination (TASCER) to induce recombination during the G2 phase so that we can generate deletions and duplications of regions of interest prior to mitosis. Using this approach in the mouse we could address the effects of copy number variation and segmental aneuploidy in daughter cells and allow us to explore somatic mosaics for large region of interest in the mouse.

Positional integratomic approach in identification of genomic candidate regions for Parkinson's disease

MOTIVATION

Recent abundance of data from studies employing high-throughput technologies to reveal alterations in human disease on genomic, transcriptomic, proteomic and other levels, offer the possibility to integrate this information into a comprehensive picture of molecular events occurring in human disease. Diversity of data originating from these studies presents a methodological obstacle in the integration process, also due to difficulties in choosing the optimal unified denominator that would allow inclusion of variables from various types of studies. We present a novel approach for integration of such multi-origin data based on positions of genetic alterations occurring in human diseases. Parkinson's disease (PD) was chosen as a model for evaluation of our methodology.

METHODS

Datasets from various types of studies in PD (linkage, genome-wide association, transcriptomic and proteomic studies) were obtained from online repositories or were extracted from available research papers. Subsequently, human genome assembly was subdivided into 10 kb regions, and significant signals from aforementioned studies were arranged into their corresponding regions according to their genomic position. For each region, rank product values were calculated and significance values were estimated by permuting the original dataset.

RESULTS

Altogether, 179 regions (representing 33 contiguous genomic regions) had significant accumulation of signals when P-value cut-off was set at 0.0001. Identified regions with significant accumulation of signals contained 29 plausible candidate genes for PD. In conclusion, we present a novel approach for identification of candidate regions and genes for various human disorders, based on the positional integration of data across various types of omic studies.

Wide-ranging functions of E2F4 in transcriptional activation and repression revealed by genome-wide analysis

The E2F family of transcription factors has important roles in cell cycle progression. E2F4 is an E2F family member that has been proposed to be primarily a repressor of transcription, but the scope of its binding activity and functions in transcriptional regulation is not fully known. We used ChIP sequencing (ChIP-seq) to identify around 16 000 E2F4 binding sites which potentially regulate 7346 downstream target genes with wide-ranging functions in DNA repair, cell cycle regulation, apoptosis, and other processes. While half of all E2F4 binding sites (56%) occurred near transcription start sites (TSSs), ∼20% of sites occurred more than 20 kb away from any annotated TSS. These distal sites showed histone modifications suggesting that E2F4 may function as a long-range regulator, which we confirmed by functional experimental assays on a subset. Overexpression of E2F4 and its transcriptional cofactors of the retinoblastoma (Rb) family and its binding partner DP-1 revealed that E2F4 acts as an activator as well as a repressor. E2F4 binding sites also occurred near regulatory elements for miRNAs such as let-7a and mir-17, suggestive of regulation of miRNAs by E2F4. Taken together, our genome-wide analysis provided evidence of versatile roles of E2F4 and insights into its functions.

When needles look like hay: how to find tissue-specific enhancers in model organism genomes

A major prerequisite for the investigation of tissue-specific processes is the identification of cis-regulatory elements. No generally applicable technique is available to distinguish them from any other type of genomic non-coding sequence. Therefore, researchers often have to identify these elements by elaborate in vivo screens, testing individual regions until the right one is found. Here, based on many examples from the literature, we summarize how functional enhancers have been isolated from other elements in the genome and how they have been characterized in transgenic animals. Covering computational and experimental studies, we provide an overview of the global properties of cis-regulatory elements, like their specific interactions with promoters and target gene distances. We describe conserved non-coding elements (CNEs) and their internal structure, nucleotide composition, binding site clustering and overlap, with a special focus on developmental enhancers. Conflicting data and unresolved questions on the nature of these elements are highlighted. Our comprehensive overview of the experimental shortcuts that have been found in the different model organism communities and the new field of high-throughput assays should help during the preparation phase of a screen for enhancers. The review is accompanied by a list of general guidelines for such a project.

Otx2 and Otx1 protect diencephalon and mesencephalon from caudalization into metencephalon during early brain regionalization

Otx2 is expressed in each step and site of head development. To dissect each Otx2 function we have identified a series of Otx2 enhancers. The Otx2 expression in the anterior neuroectoderm is regulated by the AN enhancer and the subsequent expression in forebrain and midbrain later than E8.5 by FM1 and FM2 enhancers; the Otx1 expression takes place at E8.0. In telencephalon later than E9.5 Otx1 continues to be expressed in the entire pallium, while the Otx2 expression is confined to the most medial pallium. To determine the Otx functions in forebrain and midbrain development we have generated mouse mutants that lack both FM1 and FM2 enhancers (DKO: Otx2(ΔFM1ΔFM2/ΔFM1ΔFM2)) and examined the TKO (Otx1(-/-)Otx2(ΔFM1ΔFM2/ΔFM1ΔFM2)) phenotype. The mutants develop normally until E8.0, but subsequently by E9.5 the diencephalon, including thalamic eminence and prethalamus, and the mesencephalon are caudalized into metencephalon consisting of isthmus and rhombomere 1; the caudalization does not extend to rhombomere 2 and more caudal rhombomeres. In rostral forebrain, neopallium, ganglionic eminences and hypothalamus in front of prethalamus develop; we propose that they become insensitive to the caudalization with the switch from the Otx2 expression under the AN enhancer to that under FM1 and FM2 enhancers. In contrast, the medial pallium requires Otx1 and Otx2 for its development later than E9.5, and the Otx2 expression in diencepalon and mesencephalon later than E9.5 is also directed by an enhancer other than FM1 and FM2 enhancers.

Human intronic enhancers control distinct sub-domains of Gli3 expression during mouse CNS and limb development

BACKGROUND

The zinc-finger transcription factor GLI3 is an important mediator of Sonic hedgehog signaling and crucial for patterning of many aspects of the vertebrate body plan. In vertebrates, the mechanism of SHH signal transduction and its action on target genes by means of activating or repressing forms of GLI3 have been studied most extensively during limb development and the specification of the central nervous system. From these studies it has emerged, that Gli3 expression must be subject to a tight spatiotemporal regulation. However, the genetic mechanisms and the cis-acting elements controlling the expression of Gli3 remained largely unknown.

RESULTS

Here, we demonstrate in chicken and mouse transgenic embryos that human GLI3-intronic conserved non-coding sequence elements (CNEs) autonomously control individual aspects of Gli3 expression. Their combined action shows many aspects of a Gli3-specific pattern of transcriptional activity. In the mouse limb bud, different CNEs enhance Gli3-specific expression in evolutionary ancient stylopod and zeugopod versus modern skeletal structures of the autopod. Limb bud specificity is also found in chicken but had not been detected in zebrafish embryos. Three of these elements govern central nervous system specific gene expression during mouse embryogenesis, each targeting a subset of endogenous Gli3 transcription sites. Even though fish, birds, and mammals share an ancient repertoire of gene regulatory elements within Gli3, the functions of individual enhancers from this catalog have diverged significantly. During evolution, ancient broad-range regulatory elements within Gli3 attained higher specificity, critical for patterning of more specialized structures, by abolishing the potential for redundant expression control.

CONCLUSION

These results not only demonstrate the high level of complexity in the genetic mechanisms controlling Gli3 expression, but also reveal the evolutionary significance of cis-acting regulatory networks of early developmental regulators in vertebrates.

Fine mapping of the chromosome 10q11-q21 linkage region in Alzheimer's disease cases and controls

We have previously reported strong linkage on chromosome 10q in pedigrees transmitting Alzheimer's disease through the mother, overlapping with many significant linkage reports including the largest reported study. Here, we report the most comprehensive fine mapping of this region to date. In a sample of 638 late-onset Alzheimer's disease (LOAD) cases and controls including 104 maternal LOAD cases, we genotyped 3,884 single nucleotide polymorphisms (SNPs) covering 15.2 Mb. We then used imputations and publicly available data to generate an extended dataset including 4,329 SNPs for 1,209 AD cases and 839 controls in the same region. Further, we screened eight genes in this region for rare alleles in 283 individuals by nucleotide sequencing, and we tested for possible monoallelic expression as it might underlie our maternal parent of origin linkage. We excluded the possibility of multiple rare coding risk variants for these genes and monoallelic expression when we could test for it. One SNP, rs10824310 in the PRKG1 gene, showed study-wide significant association without a parent of origin effect, but the effect size estimate is not of sufficient magnitude to explain the linkage, and no association is observed in an independent genome-wide association studies (GWAS) report. Further, no causative variants were identified though sequencing. Analysis of cases with maternal disease origin pointed to a few regions of interest that included the genes PRKG1 and PCDH15 and an intergenic interval of 200 Kb. It is likely that non-transcribed rare variants or other mechanisms involving these genomic regions underlie the observed linkage and parent of origin effect. Acquiring additional support and clarifying the mechanisms of such involvement is important for AD and other complex disorder genetics research.

Exonic remnants of whole-genome duplication reveal cis-regulatory function of coding exons

Using a comparative genomics approach to reconstruct the fate of genomic regulatory blocks (GRBs) and identify exonic remnants that have survived the disappearance of their host genes after whole-genome duplication (WGD) in teleosts, we discover a set of 38 candidate cis-regulatory coding exons (RCEs) with predicted target genes. These elements demonstrate evolutionary separation of overlapping protein-coding and regulatory information after WGD in teleosts. We present evidence that the corresponding mammalian exons are still under both coding and non-coding selection pressure, are more conserved than other protein coding exons in the host gene and several control sets, and share key characteristics with highly conserved non-coding elements in the same regions. Their dual function is corroborated by existing experimental data. Additionally, we show examples of human exon remnants stemming from the vertebrate 2R WGD. Our findings suggest that long-range cis-regulatory inputs for developmental genes are not limited to non-coding regions, but can also overlap the coding sequence of unrelated genes. Thus, exonic regulatory elements in GRBs might be functionally equivalent to those in non-coding regions, calling for a re-evaluation of the sequence space in which to look for long-range regulatory elements and experimentally test their activity.

Role of SoxB1 transcription factors in development

SoxB1 factors, which include Sox1, 2, and 3, share more than 90% amino acid identity in their DNA binding HMG box and participate in diverse developmental events. They are known to exert cell-type-specific functions in concert with other transcription factors on Sox factor-dependent regulatory enhancers. Due to the high degree of sequence similarity both within and outside the HMG box, SoxB1 members show almost identical biological activities. As a result, they exhibit strong functional redundancy in regions where SoxB1 members are coexpressed, such as neural stem/progenitor cells in the developing central nervous system.

[Conserved non-coding elements in human genome]

Study of comparative genomics has revealed that about 5% of the human genome are under purifying selection, 3.5% of which are conserved non-coding elements (CNEs). While the coding regions comprise of only a small part. In human, the CNEs are functionally important, which may be associated with the process of the establishment and maintain of chromatin architecture, transcription regulation, and pre-mRNA processing. They are also related to ontogeny of mammals and human diseases. This review outlined the identification, functional significance, evolutionary origin, and effects on human genetic defects of the CNEs.

Conservation of gene linkage in dispersed vertebrate NK homeobox clusters

Nk homeobox genes are important regulators of many different developmental processes including muscle, heart, central nervous system and sensory organ development. They are thought to have arisen as part of the ANTP megacluster, which also gave rise to Hox and ParaHox genes, and at least some NK genes remain tightly linked in all animals examined so far. The protostome-deuterostome ancestor probably contained a cluster of nine Nk genes: (Msx)-(Nk4/tinman)-(Nk3/bagpipe)-(Lbx/ladybird)-(Tlx/c15)-(Nk7)-(Nk6/hgtx)-(Nk1/slouch)-(Nk5/Hmx). Of these genes, only NKX2.6-NKX3.1, LBX1-TLX1 and LBX2-TLX2 remain tightly linked in humans. However, it is currently unclear whether this is unique to the human genome as we do not know which of these Nk genes are clustered in other vertebrates. This makes it difficult to assess whether the remaining linkages are due to selective pressures or because chance rearrangements have "missed" certain genes. In this paper, we identify all of the paralogs of these ancestrally clustered NK genes in several distinct vertebrates. We demonstrate that tight linkages of Lbx1-Tlx1, Lbx2-Tlx2 and Nkx3.1-Nkx2.6 have been widely maintained in both the ray-finned and lobe-finned fish lineages. Moreover, the recently duplicated Hmx2-Hmx3 genes are also tightly linked. Finally, we show that Lbx1-Tlx1 and Hmx2-Hmx3 are flanked by highly conserved noncoding elements, suggesting that shared regulatory regions may have resulted in evolutionary pressure to maintain these linkages. Consistent with this, these pairs of genes have overlapping expression domains. In contrast, Lbx2-Tlx2 and Nkx3.1-Nkx2.6, which do not seem to be coexpressed, are also not associated with conserved noncoding sequences, suggesting that an alternative mechanism may be responsible for the continued clustering of these genes.

Synorth: exploring the evolution of synteny and long-range regulatory interactions in vertebrate genomes

Synorth is a web resource for exploring and categorizing the syntenic relationships in gene regulatory blocks across multiple genomes.

Genomic regulatory blocks are chromosomal regions spanned by long clusters of highly conserved noncoding elements devoted to long-range regulation of developmental genes, often immobilizing other, unrelated genes into long-lasting syntenic arrangements. Synorth is a web resource for exploring and categorizing the syntenic relationships in genomic regulatory blocks across multiple genomes, tracing their evolutionary fate after teleost whole genome duplication at the level of genomic regulatory block loci, individual genes, and their phylogenetic context.

Transcriptional features of genomic regulatory blocks

CAGE tag mapping of transcription start sites across different human tissues shows that genomic regulatory blocks have unique features that are the likely cause of their ability to respond to regulatory inputs from very long distances.

Background

Genomic regulatory blocks (GRBs) are chromosomal regions spanned by highly conserved non-coding elements (HCNEs), most of which serve as regulatory inputs of one target gene in the region. The target genes are most often transcription factors involved in embryonic development and differentiation. GRBs often contain extensive gene deserts, as well as additional 'bystander' genes intertwined with HCNEs but whose expression and function are unrelated to those of the target gene. The tight regulation of target genes, complex arrangement of regulatory inputs, and the differential responsiveness of genes in the region call for the examination of fundamental rules governing transcriptional activity in GRBs. Here we use extensive CAGE tag mapping of transcription start sites across different human tissues and differentiation stages combined with expression data and a number of sequence and epigenetic features to discover these rules and patterns.

Results

We show evidence that GRB target genes have properties that set them apart from their bystanders as well as other genes in the genome: longer CpG islands, a higher number and wider spacing of alternative transcription start sites, and a distinct composition of transcription factor binding sites in their core/proximal promoters. Target gene expression correlates with the acetylation state of HCNEs in the region. Additionally, target gene promoters have a distinct combination of activating and repressing histone modifications in mouse embryonic stem cell lines.

Conclusions

GRB targets are genes with a number of unique features that are the likely cause of their ability to respond to regulatory inputs from very long distances.

Systematic human/zebrafish comparative identification of cis-regulatory activity around vertebrate developmental transcription factor genes

Pan-vertebrate developmental cis-regulatory elements are discernible as highly conserved noncoding elements (HCNEs) and are often dispersed over large areas around the pleiotropic genes whose expression they control. On the loci of two developmental transcription factor genes, SOX3 and PAX6, we demonstrate that HCNEs conserved between human and zebrafish can be systematically and reliably tested for their regulatory function in multiple stable transgenes in zebrafish, and their genomic reach estimated with confidence using synteny conservation and HCNE density along these loci. HCNEs of both human and zebrafish function as specific developmental enhancers in zebrafish. We show that human HCNEs result in expression patterns in zebrafish equivalent to those in mouse, establishing zebrafish as a suitable model for large-scale testing of human developmental enhancers. Orthologous human and zebrafish enhancers underwent functional evolution within their sequence and often directed related but non-identical expression patterns. Despite an evolutionary distance of 450 million years, one pax6 HCNE drove expression in identical areas when comparing zebrafish vs. human HCNEs. HCNEs from the same area often drive overlapping patterns, suggesting that multiple regulatory inputs are required to achieve robust and precise complex expression patterns exhibited by developmental genes.

Positional distribution of human transcription factor binding sites

We developed a method for estimating the positional distribution of transcription factor (TF) binding sites using ChIP-chip data, and applied it to recently published experiments on binding sites of nine TFs: OCT4, SOX2, NANOG, HNF1A, HNF4A, HNF6, FOXA2, USF1 and CREB1. The data were obtained from a genome-wide coverage of promoter regions from 8-kb upstream of the transcription start site (TSS) to 2-kb downstream. The number of target genes of each TF ranges from few hundred to several thousand. We found that for each of the nine TFs the estimated binding site distribution is closely approximated by a mixture of two components: a narrow peak, localized within 300-bp upstream of the TSS, and a distribution of almost uniform density within the tested region. Using Gene Ontology (GO) and Enrichment analysis, we were able to associate (for each of the TFs studied) the target genes of both types of binding with known biological processes. Most GO terms were enriched either among the proximal targets or among those with a uniform distribution of binding sites. For example, the three stemness-related TFs have several hundred target genes that belong to ‘development’ and ‘morphogenesis’ whose binding sites belong to the uniform distribution.

Searching for functional genetic variants in non-coding DNA

1. The search for DNA sequence variants for complex human polygenic conditions has been a strong focus of recent genetic research. While gene loci have been identified, few variants in the coding sequences of these genes have been found, suggesting that non-coding sequence variation may underlie many complex conditions. 2. Non-coding DNA harbours regulatory elements capable of making changes to gene expression. However, regulatory DNA sequences are currently difficult to recognize and their function is poorly understood, complicating the task of assigning potential functional significance to non-coding variation. 3. Comparative genomics, the study of evolutionary DNA conservation, has enabled the emergent field of non-coding DNA identification in human disease analysis. 4. This brief review will focus on the potential of a relatively high throughput technique based on comparative genomics, that may aid in the identification of functionally important non-coding sequence variation in complex diseases.

Possible involvement of SINEs in mammalian-specific brain formation

Retroposons, such as short interspersed elements (SINEs) and long interspersed elements (LINEs), are the major constituents of higher vertebrate genomes. Although there are many examples of retroposons' acquiring function, none has been implicated in the morphological innovations specific to a certain taxonomic group. We previously characterized a SINE family, AmnSINE1, members of which constitute a part of conserved noncoding elements (CNEs) in mammalian genomes. We proposed that this family acquired genomic functionality or was exapted after retropositioning in a mammalian ancestor. Here we identified 53 new AmnSINE1 loci and refined 124 total loci, two of which were further analyzed. Using a mouse enhancer assay, we demonstrate that one SINE locus, AS071, 178 kbp from the gene FGF8 (fibroblast growth factor 8), is an enhancer that recapitulates FGF8 expression in two regions of the developing forebrain, namely the diencephalon and the hypothalamus. Our gain-of-function analysis revealed that FGF8 expression in the diencephalon controls patterning of thalamic nuclei, which act as a relay center of the neocortex, suggesting a role for FGF8 in mammalian-specific forebrain patterning. Furthermore, we demonstrated that the locus, AS021, 392 kbp from the gene SATB2, controls gene expression in the lateral telencephalon, which is thought to be a signaling center during development. These results suggest important roles for SINEs in the development of the mammalian neuronal network, a part of which was initiated with the exaptation of AmnSINE1 in a common mammalian ancestor.

Ancora: a web resource for exploring highly conserved noncoding elements and their association with developmental regulatory genes

Ancora is a web resource that provides data and tools for exploring genomic organization of highly conserved noncoding elements for multiple genomes.

Metazoan genomes contain arrays of highly conserved noncoding elements (HCNEs) that span developmental regulatory genes and define regulatory domains. We describe Ancora , a web resource that provides data and tools for exploring genomic organization of HCNEs for multiple genomes. Ancora includes a genome browser that shows HCNE locations and features novel HCNE density plots as a powerful tool to discover developmental regulatory genes and distinguish their regulatory elements and domains.

Long-range gene control and genetic disease

The past two decades have seen great progress in the elucidation of the genetic basis of human genetic disease. Many clinical phenotypes have been linked with mutations or deletions in specific causative genes. However, it is often less recognized that in addition to the integrity of the protein-coding sequences, human health critically also depends on the spatially, temporally, and quantitatively correct expression of those genes. Genetic disease can therefore equally be caused by disruption of the regulatory mechanisms that ensure proper gene expression. The term "position effect" is used in those situations where the expression level of a gene is deleteriously affected by an alteration in its chromosomal environment, while maintaining an intact transcription unit. Here, we review recent advances in our understanding of the possible mechanisms of a number of "position effect" disease cases and discuss the findings with respect to current models for genome organization and long-range control of gene expression.

Retroviral enhancer detection insertions in zebrafish combined with comparative genomics reveal genomic regulatory blocks - a fundamental feature of vertebrate genomes

A large-scale enhancer detection screen was performed in the zebrafish using a retroviral vector carrying a basal promoter and a fluorescent protein reporter cassette. Analysis of insertional hotspots uncovered areas around developmental regulatory genes in which an insertion results in the same global expression pattern, irrespective of exact position. These areas coincide with vertebrate chromosomal segments containing identical gene order; a phenomenon known as conserved synteny and thought to be a vestige of evolution. Genomic comparative studies have found large numbers of highly conserved noncoding elements (HCNEs) spanning these and other loci. HCNEs are thought to act as transcriptional enhancers based on the finding that many of those that have been tested direct tissue specific expression in transient or transgenic assays. Although gene order in hox and other gene clusters has long been known to be conserved because of shared regulatory sequences or overlapping transcriptional units, the chromosomal areas found through insertional hotspots contain only one or a few developmental regulatory genes as well as phylogenetically unrelated genes. We have termed these regions genomic regulatory blocks (GRBs), and show that they underlie the phenomenon of conserved synteny through all sequenced vertebrate genomes. After teleost whole genome duplication, a subset of GRBs were retained in two copies, underwent degenerative changes compared with tetrapod loci that exist as single copy, and that therefore can be viewed as representing the ancestral form. We discuss these findings in light of evolution of vertebrate chromosomal architecture and the identification of human disease mutations.

Genomic regulatory blocks underlie extensive microsynteny conservation in insects

Insect genomes contain larger blocks of conserved gene order (microsynteny) than would be expected under a random breakage model of chromosome evolution. We present evidence that microsynteny has been retained to keep large arrays of highly conserved noncoding elements (HCNEs) intact. These arrays span key developmental regulatory genes, forming genomic regulatory blocks (GRBs). We recently described GRBs in vertebrates, where most HCNEs function as enhancers and HCNE arrays specify complex expression programs of their target genes. Here we present a comparison of five Drosophila genomes showing that HCNE density peaks centrally in large synteny blocks containing multiple genes. Besides developmental regulators that are likely targets of HCNE enhancers, HCNE arrays often span unrelated neighboring genes. We describe differences in core promoters between the target genes and the unrelated genes that offer an explanation for the differences in their responsiveness to enhancers. We show examples of a striking correspondence between boundaries of synteny blocks, HCNE arrays, and Polycomb binding regions, confirming that the synteny blocks correspond to regulatory domains. Although few noncoding elements are highly conserved between Drosophila and the malaria mosquito Anopheles gambiae, we find that A. gambiae regions orthologous to Drosophila GRBs contain an equivalent distribution of noncoding elements highly conserved in the yellow fever mosquito Aëdes aegypti and coincide with regions of ancient microsynteny between Drosophila and mosquitoes. The structural and functional equivalence between insect and vertebrate GRBs marks them as an ancient feature of metazoan genomes and as a key to future studies of development and gene regulation.

Developmental malformations of the eye: the role of PAX6, SOX2 and OTX2

Eye development initiates as an evagination of the early neural plate, before the closure of the neural tube. Structural malformations of the eye such as anophthalmia and microphthalmia arise very early in development. It is not surprising therefore that three of the genes currently identified to play a significant role in these developmental eye anomalies are also major players in brain development and regionalization. However, as has been emerging for a high proportion of transcriptional regulators studied, these genes have evolved to play multiple roles throughout development, and perhaps even in adult tissue maintenance. This complex spatiotemporal expression pattern requires elaborate regulatory systems which we are beginning to unravel. A major component of these complex regulatory networks is a series of cis-acting elements, highly conserved through evolution, which spread large distances from the coding region of each gene. We describe how cross regulation for PAX6, SOX2 and perhaps OTX2 has now been uncovered, pointing to the mechanisms that can fine-tune the expression of such essential developmental components. These interactions also help us understand why there is significant phenotypic overlap between mutations at these three loci.

Wide-scale analysis of human functional transcription factor binding reveals a strong bias towards the transcription start site

Background

Transcription factors (TF) regulate expression by binding to specific DNA sequences. A binding event is functional when it affects gene expression. Functionality of a binding site is reflected in conservation of the binding sequence during evolution and in over represented binding in gene groups with coherent biological functions. Functionality is governed by several parameters such as the TF-DNA binding strength, distance of the binding site from the transcription start site (TSS), DNA packing, and more. Understanding how these parameters control functionality of different TFs in different biological contexts is a must for identifying functional TF binding sites and for understanding regulation of transcription.

Methodology/Principal Findings

We introduce a novel method to screen the promoters of a set of genes with shared biological function (obtained from the functional Gene Ontology (GO) classification) against a precompiled library of motifs, and find those motifs which are statistically over-represented in the gene set. More than 8000 human (and 23,000 mouse) genes, were assigned to one of 134 GO sets. Their promoters were searched (from 200 bp downstream to 1000 bp upstream the TSS) for 414 known DNA motifs. We optimized the sequence similarity score threshold, independently for every location window, taking into account nucleotide heterogeneity along the promoters of the target genes. The method, combined with binding sequence and location conservation between human and mouse, identifies with high probability functional binding sites for groups of functionally-related genes. We found many location-sensitive functional binding events and showed that they clustered close to the TSS. Our method and findings were tested experimentally.

Conclusions/Significance

We identified reliably functional TF binding sites. This is an essential step towards constructing regulatory networks. The promoter region proximal to the TSS is of central importance for regulation of transcription in human and mouse, just as it is in bacteria and yeast.

Cis-regulation of the amphioxus engrailed gene: Insights into evolution of a muscle-specific enhancer

To gain insights into the relation between evolution of cis-regulatory DNA and evolution of gene function, we identified tissue-specific enhancers of the engrailed gene of the basal chordate amphioxus (Branchiostoma floridae) and compared their ability to direct expression in both amphioxus and its nearest chordate relative, the tunicate Ciona intestinalis. In amphioxus embryos, the native engrailed gene is expressed in three domains - the eight most anterior somites, a few cells in the central nervous system (CNS) and a few ectodermal cells. In contrast, in C. intestinalis, in which muscle development is highly divergent, engrailed expression is limited to the CNS. To characterize the tissue-specific enhancers of amphioxus engrailed, we first showed that 7.8kb of upstream DNA of amphioxus engrailed directs expression to all three domains in amphioxus that express the native gene. We then identified the amphioxus engrailed muscle-specific enhancer as the 1.2kb region of upstream DNA with the highest sequence identity to the mouse en-2 jaw muscle enhancer. This amphioxus enhancer directed expression to both the somites in amphioxus and to the larval muscles in C. intestinalis. These results show that even though expression of the native engrailed has apparently been lost in developing C. intestinalis muscles, they express the transcription factors necessary to activate transcription from the amphioxus engrailed enhancer, suggesting that gene networks may not be completely disrupted if an individual component is lost.

Crucial roles of Foxa2 in mouse anterior-posterior axis polarization via regulation of anterior visceral endoderm-specific genes

Anterior visceral endoderm (AVE) plays essential roles with respect to anterior-posterior axis development in the early mouse embryo. To assess the genetic cascade involved in AVE formation, the cis-regulatory elements directing expression of vertebrate Otx2 genes in the AVE were analyzed via generation of transgenic mice. Otx2 expression in AVE is regulated directly by the forkhead transcription factor, Foxa2. Moreover, Foxa2 is essential for expression of the Wnt antagonists, Dkk1 and Cerl, in visceral endoderm during the pre- to early streak stages; however, Foxa2 appears to be dispensable for subsequent Dkk1 expression associated with forebrain induction. Thus, we propose that Foxa2 is crucial in early anterior-posterior axis polarization in terms of regulation of expression of AVE-specific genes. These findings provide profound insights into conserved roles of Foxa2 transcription factors in anterior specification throughout the evolution of the chordate body plan.

Two or four bristles: functional evolution of an enhancer of scute in Drosophilidae

Changes in cis-regulatory sequences are proposed to underlie much of morphological evolution. Yet, little is known about how such modifications translate into phenotypic differences. To address this problem, we focus on the dorsocentral bristles of Drosophilidae. In Drosophila melanogaster, development of these bristles depends on a cis-regulatory element, the dorsocentral enhancer, to activate scute in a cluster of cells from which two bristles on the posterior scutum arise. A few species however, such as D. quadrilineata, bear anterior dorsocentral bristles as well as posterior ones, a derived feature. This correlates with an anterior expansion of the scute expression domain. Here, we show that the D. quadrilineata enhancer has evolved, and is now active in more anterior regions. When used to rescue scute expression in transgenic D. melanogaster, the D. quadrilineata enhancer is able to induce anterior bristles. Importantly, these properties are not displayed by homologous enhancers from control species bearing only two posterior bristles. We also provide evidence that upstream regulation of the enhancer, by the GATA transcription factor Pannier, has been evolutionarily conserved. This work illustrates how, in the context of a conserved trans-regulatory landscape, evolutionary tinkering of pre-existing enhancers can modify gene expression patterns and contribute to morphological diversification.

Evolutionary change in function of the dorsocentral enhancer (DCE) of scute has resulted in altered bristle formation between two species of Drosophila.

The Fugu tyrp1 promoter directs specific GFP expression in zebrafish: tools to study the RPE and the neural crest-derived melanophores

In vertebrates, pigment cells account for a small percentage of the total cell population and they intermingle with other cell types. This makes it difficult to isolate them for analyzes of their functions in the context of development. To alleviate such difficulty, we generated two stable transgenic zebrafish lines (pt101 and pt102) that express green fluorescent protein (GFP) in melanophores under the control of the 1 kb Fugu tyrp1 promoter. In pt101, GFP is expressed in both retinal pigment epithelium (RPE) cells and the neural crest-derived melanophores (NCDM), whereas in pt102, GFP is predominately expressed in the NCDM. Our results indicate that the Fugu tyrp1 promoter can direct transgene expression in a cell-type-specific manner in zebrafish. In addition, our findings provide evidence supporting differential regulations of melanin-synthesizing genes in RPE cells and the NCDM in zebrafish. Utilizing the varying GFP expression levels in these fish, we have isolated melanophores via flow cytometry and revealed the capability of sorting the NCDM from RPE cells as well. Thus, these transgenic lines are useful tools to study melanophores in zebrafish.

Comprehensive characterization of the cis-regulatory code responsible for the spatio-temporal expression of olSix3.2 in the developing medaka forebrain

BACKGROUND

Embryonic development is coordinated by sets of cis-regulatory elements that are collectively responsible for the precise spatio-temporal organization of regulatory gene networks. There is little information on how these elements, which are often associated with highly conserved noncoding sequences, are combined to generate precise gene expression patterns in vertebrates. To address this issue, we have focused on Six3, an important regulator of vertebrate forebrain development.

RESULTS

Using computational analysis and exploiting the diversity of teleost genomes, we identified a cluster of highly conserved noncoding sequences surrounding the Six3 gene. Transgenesis in medaka fish demonstrates that these sequences have enhancer, silencer, and silencer blocker activities that are differentially combined to control the entire distribution of Six3.

CONCLUSION

This report provides the first example of the precise regulatory code necessary for the expression of a vertebrate gene, and offers a unique framework for defining the interplay of trans-acting factors that control the evolutionary conserved use of Six3.

Functional single nucleotide polymorphism-based association studies

Association studies hold great promise for the elucidation of the genetic basis of diseases. Studies based on functional single nucleotide polymorphisms (SNPs) or on linkage disequilibrium (LD) represent two main types of designs. LD-based association studies can be comprehensive for common causative variants, but they perform poorly for rare alleles. Conversely, functional SNP-based studies are efficient because they focus on the SNPs with the highest a priori chance of being associated. Our poor ability to predict the functional effect of SNPs, however, hampers attempts to make these studies comprehensive. Recent progress in comparative genomics, and evidence that functional elements tend to lie in conserved regions, promises to change the landscape, permitting functional SNP association studies to be carried out that comprehensively assess common and rare alleles. SNP genotyping technologies are already sufficient for such studies, but studies will require continued genomic sequencing of multiple species, research on the functional role of conserved sequences and additional SNP discovery and validation efforts (including targeted SNP discovery to identify the rare alleles in functional regions). With these resources, we expect that comprehensive functional SNP association studies will soon be possible.

Strong bias in the location of functional promoter polymorphisms

A considerable proportion of heritable human phenotypic variation is thought to result from altered gene expression. Unfortunately, it is currently impossible to use bioinformatic analysis to discriminate between DNA sequence variants that are likely to influence gene expression and those that are not. In an attempt to define some of the characteristics of promoter polymorphisms with functional effects on gene expression, we examined 674 haplotypes representing 247 unique gene promoters using a standardized reporter gene assay system. Sequence variants that altered gene expression by 1.5-fold or more were strongly biased toward a location in the core and proximal promoter regions, 50% being within the first 100 bases 5' to the transcription start site. No bias was seen in the allele frequencies of functional and nonfunctional sequence variants. Only 33% of the functional variants were found in known consensus transcription factor binding sequences or motifs, which suggests that either there are many unknown transcription factor binding motifs or other, unknown mechanisms are involved. The genes with functional polymorphisms that are reported here for the first time include AGTRL2, CAT, CHRNA5, CTSG, CYP2D6, DLD, ERCC1, GABRA1, GABRP, HNRPH3, HIP1, IGKV1-9, KCNJ15, KCNK6, KLK1, MSMB, MYOC, NPY2R, NOTCH4, ORM2, PEDF, PTPRCAP, ST16 (IL24), SULT1A1, and TSHR.

New technologies, new findings, and new concepts in the study of vertebrate cis-regulatory sequences

All vertebrates share a similar early embryonic body plan and use the same regulatory genes for their development. The availability of numerous sequenced vertebrate genomes and significant advances in bioinformatics have resulted in the finding that the genomic regions of many of these developmental regulatory genes also contain highly conserved noncoding sequence. In silico discovery of conserved noncoding regions and of transcription factor binding sites as well as the development of methods for high throughput transgenesis in Xenopus and zebrafish are dramatically increasing the speed with which regulatory elements can be discovered, characterized, and tested in the context of whole live embryos. We review here some of the recent technological developments that will likely lead to a surge in research on how vertebrate genomes encode regulation of transcriptional activity, how regulatory sequences constrain genomic architecture, and ultimately how vertebrate form has evolved.

Ancient duplicated conserved noncoding elements in vertebrates: a genomic and functional analysis

Fish-mammal genomic comparisons have proved powerful in identifying conserved noncoding elements likely to be cis-regulatory in nature, and the majority of those tested in vivo have been shown to act as tissue-specific enhancers associated with genes involved in transcriptional regulation of development. Although most of these elements share little sequence identity to each other, a small number are remarkably similar and appear to be the product of duplication events. Here, we searched for duplicated conserved noncoding elements in the human genome, using comparisons with Fugu to select putative cis-regulatory sequences. We identified 124 families of duplicated elements, each containing between two and five members, that are highly conserved within and between vertebrate genomes. In 74% of cases, we were able to assign a specific set of paralogous genes with annotation relating to transcriptional regulation and/or development to each family, thus removing much of the ambiguity in identifying associated genes. We find that duplicate elements have the potential to up-regulate reporter gene expression in a tissue-specific manner and that expression domains often overlap, but are not necessarily identical, between family members. Over two thirds of the families are conserved in duplicate in fish and appear to predate the large-scale duplication events thought to have occurred at the origin of vertebrates. We propose a model whereby gene duplication and the evolution of cis-regulatory elements can be considered in the context of increased morphological diversity and the emergence of the modern vertebrate body plan.