Big genomes facilitate the comparative identification of regulatory elements

Genome sequencing and analysis of black flounder (Paralichthys orbignyanus) reveals new insights into Pleuronectiformes genomic size and structure

Black flounder (Paralichthys orbignyanus, Pleuronectiformes) is a commercially significant marine fish with promising aquaculture potential in Argentina. Despite extensive studies on Black flounder aquaculture, its limited genetic information available hampers the crucial role genetics plays in the development of this activity. In this study, we first employed Illumina sequencing technology to sequence the entire genome of Black flounder. Utilizing two independent libraries-one from a female and another from a male-with 150 bp paired-end reads, a mean insert length of 350 bp, and over 35 X-fold coverage, we achieved assemblies resulting in a genome size of ~ 538 Mbp. Analysis of the assemblies revealed that more than 98% of the core genes were present, with more than 78% of them having more than 50% coverage. This indicates a somehow complete and accurate genome at the coding sequence level. This genome contains 25,231 protein-coding genes, 445 tRNAs, 3 rRNAs, and more than 1,500 non-coding RNAs of other types. Black flounder, along with pufferfishes, seahorses, pipefishes, and anabantid fish, displays a smaller genome compared to most other teleost groups. In vertebrates, the number of transposable elements (TEs) is often correlated with genome size. However, it remains unclear whether the sizes of introns and exons also play a role in determining genome size. Hence, to elucidate the potential factors contributing to this reduced genome size, we conducted a comparative genomic analysis between Black flounder and other teleost orders to determine if the small genomic size could be explained by repetitive elements or gene features, including the whole genome genes and introns sizes. We show that the smaller genome size of flounders can be attributed to several factors, including changes in the number of repetitive elements, and decreased gene size, particularly due to lower amount of very large and small introns. Thus, these components appear to be involved in the genome reduction in Black flounder. Despite these insights, the full implications and potential benefits of genome reduction in Black flounder for reproduction and aquaculture remain incompletely understood, necessitating further research.

The Density of Regulatory Information Is a Major Determinant of Evolutionary Constraint on Noncoding DNA in Drosophila

Evolutionary analyses have estimated that ∼60% of nucleotides in intergenic regions of the Drosophila melanogaster genome are functionally relevant, suggesting that regulatory information may be encoded more densely in intergenic regions than has been revealed by most functional dissections of regulatory DNA. Here, we approached this issue through a functional dissection of the regulatory region of the gene shavenbaby (svb). Most of the ∼90 kb of this large regulatory region is highly conserved in the genus Drosophila, though characterized enhancers occupy a small fraction of this region. By analyzing the regulation of svb in different contexts of Drosophila development, we found that the regulatory information that drives svb expression in the abdominal pupal epidermis is organized in a different way than the elements that drive svb expression in the embryonic epidermis. While in the embryonic epidermis svb is activated by compact enhancers separated by large inactive DNA regions, svb expression in the pupal epidermis is driven by regulatory information distributed over broader regions of svb cis-regulatory DNA. In the same vein, we observed that other developmental genes also display a dense distribution of putative regulatory elements in their regulatory regions. Furthermore, we found that a large percentage of conserved noncoding DNA of the Drosophila genome is contained within regions of open chromatin. These results suggest that part of the evolutionary constraint on noncoding DNA of Drosophila is explained by the density of regulatory information, which may be greater than previously appreciated.

Anterior-posterior patterning of segments in Anopheles stephensi offers insights into the transition from sequential to simultaneous segmentation in holometabolous insects

The gene regulatory network for segmentation in arthropods offers valuable insights into how networks evolve owing to the breadth of species examined and the extremely detailed knowledge gained in the model organism Drosophila melanogaster. These studies have shown that Drosophila's network represents a derived state that acquired changes to accelerate segment patterning, whereas most insects specify segments gradually as the embryo elongates. Such heterochronic shifts in segmentation have potentially emerged multiple times within holometabolous insects, resulting in many mechanistic variants and difficulties in isolating underlying commonalities that permit such shifts. Recent studies identified regulatory genes that work as timing factors, coordinating gene expression transitions during segmentation. These studies predict that changes in timing factor deployment explain shifts in segment patterning relative to other developmental events. Here, we test this hypothesis by characterizing the temporal and spatial expression of the pair-rule patterning genes in the malaria vector mosquito, Anopheles stephensi. This insect is a Dipteran (fly), like Drosophila, but represents an ancient divergence within this clade, offering a useful counterpart for evo-devo studies. In mosquito embryos, we observe anterior to posterior sequential addition of stripes for many pair-rule genes and a wave of broad timer gene expression across this axis. Segment polarity gene stripes are added sequentially in the wake of the timer gene wave and the full pattern is not complete until the embryo is fully elongated. This "progressive segmentation" mode in Anopheles displays commonalities with both Drosophila's rapid segmentation mechanism and sequential modes used by more distantly related insects.

Favorable genomic environments for cis-regulatory evolution: A novel theoretical framework

Cis-regulatory changes are arguably the primary evolutionary source of animal morphological diversity. With the recent explosion of genome-wide comparisons of the cis-regulatory content in different animal species is now possible to infer general principles underlying enhancer evolution. However, these studies have also revealed numerous discrepancies and paradoxes, suggesting that the mechanistic causes and modes of cis-regulatory evolution are still not well understood and are probably much more complex than generally appreciated. Here, we argue that the mutational mechanisms and genomic regions generating new regulatory activities must comply with the constraints imposed by the molecular properties of cis-regulatory elements (CREs) and the organizational features of long-range chromatin interactions. Accordingly, we propose a new integrative evolutionary framework for cis-regulatory evolution based on two major premises for the origin of novel enhancer activity: (i) an accessible chromatin environment and (ii) compatibility with the 3D structure and interactions of pre-existing CREs. Mechanisms and DNA sequences not fulfilling these premises, will be less likely to have a measurable impact on gene expression and as such, will have a minor contribution to the evolution of gene regulation. Finally, we discuss current comparative cis-regulatory data under the light of this new evolutionary model, and propose that the two most prominent mechanisms for the evolution of cis-regulatory changes are the overprinting of ancestral CREs and the exaptation of transposable elements.

The achaete-scute complex in Diptera: patterns of noncoding sequence evolution

The achaete‐scute complex (AS‐C) has been a useful paradigm for the study of pattern formation and its evolution. achaete‐scute genes have duplicated and evolved distinct expression patterns during the evolution of cyclorraphous Diptera. Are the expression patterns in different species driven by conserved regulatory elements? If so, when did such regulatory elements arise? Here, we have sequenced most of the AS‐C of the fly Calliphora vicina (including the genes achaete, scute and lethal of scute) to compare noncoding sequences with known cis‐regulatory sequences in Drosophila. The organization of the complex is conserved with respect to Drosophila species. There are numerous small stretches of conserved noncoding sequence that, in spite of high sequence turnover, display binding sites for known transcription factors. Synteny of the blocks of conserved noncoding sequences is maintained suggesting not only conservation of the position of regulatory elements but also an origin prior to the divergence between these two species. We propose that some of these enhancers originated by duplication with their target genes.

Evidence for deep regulatory similarities in early developmental programs across highly diverged insects

Many genes familiar from Drosophila development, such as the so-called gap, pair-rule, and segment polarity genes, play important roles in the development of other insects and in many cases appear to be deployed in a similar fashion, despite the fact that Drosophila-like "long germband" development is highly derived and confined to a subset of insect families. Whether or not these similarities extend to the regulatory level is unknown. Identification of regulatory regions beyond the well-studied Drosophila has been challenging as even within the Diptera (flies, including mosquitoes) regulatory sequences have diverged past the point of recognition by standard alignment methods. Here, we demonstrate that methods we previously developed for computational cis-regulatory module (CRM) discovery in Drosophila can be used effectively in highly diverged (250-350 Myr) insect species including Anopheles gambiae, Tribolium castaneum, Apis mellifera, and Nasonia vitripennis. In Drosophila, we have successfully used small sets of known CRMs as "training data" to guide the search for other CRMs with related function. We show here that although species-specific CRM training data do not exist, training sets from Drosophila can facilitate CRM discovery in diverged insects. We validate in vivo over a dozen new CRMs, roughly doubling the number of known CRMs in the four non-Drosophila species. Given the growing wealth of Drosophila CRM annotation, these results suggest that extensive regulatory sequence annotation will be possible in newly sequenced insects without recourse to costly and labor-intensive genome-scale experiments. We develop a new method, Regulus, which computes a probabilistic score of similarity based on binding site composition (despite the absence of nucleotide-level sequence alignment), and demonstrate similarity between functionally related CRMs from orthologous loci. Our work represents an important step toward being able to trace the evolutionary history of gene regulatory networks and defining the mechanisms underlying insect evolution.

Shadow enhancers enable Hunchback bifunctionality in the Drosophila embryo

Hunchback (Hb) is a bifunctional transcription factor that activates and represses distinct enhancers. Here, we investigate the hypothesis that Hb can activate and repress the same enhancer. Computational models predicted that Hb bifunctionally regulates the even-skipped (eve) stripe 3+7 enhancer (eve3+7) in Drosophila blastoderm embryos. We measured and modeled eve expression at cellular resolution under multiple genetic perturbations and found that the eve3+7 enhancer could not explain endogenous eve stripe 7 behavior. Instead, we found that eve stripe 7 is controlled by two enhancers: the canonical eve3+7 and a sequence encompassing the minimal eve stripe 2 enhancer (eve2+7). Hb bifunctionally regulates eve stripe 7, but it executes these two activities on different pieces of regulatory DNA--it activates the eve2+7 enhancer and represses the eve3+7 enhancer. These two "shadow enhancers" use different regulatory logic to create the same pattern.

Quantitative system drift compensates for altered maternal inputs to the gap gene network of the scuttle fly Megaselia abdita

The segmentation gene network in insects can produce equivalent phenotypic outputs despite differences in upstream regulatory inputs between species. We investigate the mechanistic basis of this phenomenon through a systems-level analysis of the gap gene network in the scuttle fly Megaselia abdita (Phoridae). It combines quantification of gene expression at high spatio-temporal resolution with systematic knock-downs by RNA interference (RNAi). Initiation and dynamics of gap gene expression differ markedly between M. abdita and Drosophila melanogaster, while the output of the system converges to equivalent patterns at the end of the blastoderm stage. Although the qualitative structure of the gap gene network is conserved, there are differences in the strength of regulatory interactions between species. We term such network rewiring 'quantitative system drift'. It provides a mechanistic explanation for the developmental hourglass model in the dipteran lineage. Quantitative system drift is likely to be a widespread mechanism for developmental evolution.

The draft genome of the pest tephritid fruit fly Bactrocera tryoni: resources for the genomic analysis of hybridising species

Background

The tephritid fruit flies include a number of economically important pests of horticulture, with a large accumulated body of research on their biology and control. Amongst the Tephritidae, the genus Bactrocera, containing over 400 species, presents various species groups of potential utility for genetic studies of speciation, behaviour or pest control. In Australia, there exists a triad of closely-related, sympatric Bactrocera species which do not mate in the wild but which, despite distinct morphologies and behaviours, can be force-mated in the laboratory to produce fertile hybrid offspring. To exploit the opportunities offered by genomics, such as the efficient identification of genetic loci central to pest behaviour and to the earliest stages of speciation, investigators require genomic resources for future investigations.

Results

We produced a draft de novo genome assembly of Australia’s major tephritid pest species, Bactrocera tryoni. The male genome (650 -700 Mbp) includes approximately 150Mb of interspersed repetitive DNA sequences and 60Mb of satellite DNA. Assessment using conserved core eukaryotic sequences indicated 98% completeness. Over 16,000 MAKER-derived gene models showed a large degree of overlap with other Dipteran reference genomes. The sequence of the ribosomal RNA transcribed unit was also determined. Unscaffolded assemblies of B. neohumeralis and B. jarvisi were then produced; comparison with B. tryoni showed that the species are more closely related than any Drosophila species pair. The similarity of the genomes was exploited to identify 4924 potentially diagnostic indels between the species, all of which occur in non-coding regions.

Conclusions

This first draft B. tryoni genome resembles other dipteran genomes in terms of size and putative coding sequences. For all three species included in this study, we have identified a comprehensive set of non-redundant repetitive sequences, including the ribosomal RNA unit, and have quantified the major satellite DNA families. These genetic resources will facilitate the further investigations of genetic mechanisms responsible for the behavioural and morphological differences between these three species and other tephritids. We have also shown how whole genome sequence data can be used to generate simple diagnostic tests between very closely-related species where only one of the species is scaffolded.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1153) contains supplementary material, which is available to authorized users.

Assessing constraints on the path of regulatory sequence evolution

Structural and functional constraints are known to play a major role in restricting the path of evolution of protein activities. However, constraints acting on evolving transcriptional regulatory sequences, e.g. enhancers, are largely unknown. Recently, we elucidated how a novel expression pattern of the Neprilysin-1 (Nep1) gene in the optic lobe of Drosophila santomea evolved via co-option of existing enhancer activities. Drosophila santomea, which has diverged from Drosophila yakuba by approximately 400 000 years has accumulated four fixed mutations that each contribute to the full activity of this enhancer. Recreating and testing the optic lobe enhancer of the ancestor of D. santomea and D. yakuba revealed that the strong D. santomea enhancer activity evolved from a weak ancestral activity. Because each mutation on the path from the D. yakuba/santomea ancestor to modern-day D. santomea contributes to the newly derived optic lobe enhancer activity, we sought here to use this system to study the path of evolution of enhancer sequences. We inferred likely paths of evolution of this enhancer by observing the transcriptional output of all possible intermediate steps between the ancestral D. yakuba/santomea enhancer and the modern D. santomea enhancer. Many possible paths had epistatic and cooperative effects. Furthermore, we found that several paths significantly increased ectopic transcriptional activity or affected existing enhancer activities from which the novel activity was co-opted. We suggest that these attributes highlight constraints that guide the path of evolution of enhancers.

The evolution of small insertions and deletions in the coding genes of Drosophila melanogaster

Studies of protein evolution have focused on amino acid substitutions with much less systematic analysis on insertion and deletions (indels) in protein coding genes. We hence surveyed 7,500 genes between Drosophila melanogaster and D. simulans, using D. yakuba as an outgroup for this purpose. The evolutionary rate of coding indels is indeed low, at only 3% of that of nonsynonymous substitutions. As coding indels follow a geometric distribution in size and tend to fall in low-complexity regions of proteins, it is unclear whether selection or mutation underlies this low rate. To resolve the issue, we collected genomic sequences from an isogenic African line of D. melanogaster (ZS30) at a high coverage of 70× and analyzed indel polymorphism between ZS30 and the reference genome. In comparing polymorphism and divergence, we found that the divergence to polymorphism ratio (i.e., fixation index) for smaller indels (size ≤ 10 bp) is very similar to that for synonymous changes, suggesting that most of the within-species polymorphism and between-species divergence for indels are selectively neutral. Interestingly, deletions of larger sizes (size ≥ 11 bp and ≤ 30 bp) have a much higher fixation index than synonymous mutations and 44.4% of fixed middle-sized deletions are estimated to be adaptive. To our surprise, this pattern is not found for insertions. Protein indel evolution appear to be in a dynamic flux of neutrally driven expansion (insertions) together with adaptive-driven contraction (deletions), and these observations provide important insights for understanding the fitness of new mutations as well as the evolutionary driving forces for genomic evolution in Drosophila species.

Comparative transcriptomics of early dipteran development

Background

Modern sequencing technologies have massively increased the amount of data available for comparative genomics. Whole-transcriptome shotgun sequencing (RNA-seq) provides a powerful basis for comparative studies. In particular, this approach holds great promise for emerging model species in fields such as evolutionary developmental biology (evo-devo).

Results

We have sequenced early embryonic transcriptomes of two non-drosophilid dipteran species: the moth midge Clogmia albipunctata, and the scuttle fly Megaselia abdita. Our analysis includes a third, published, transcriptome for the hoverfly Episyrphus balteatus. These emerging models for comparative developmental studies close an important phylogenetic gap between Drosophila melanogaster and other insect model systems. In this paper, we provide a comparative analysis of early embryonic transcriptomes across species, and use our data for a phylogenomic re-evaluation of dipteran phylogenetic relationships.

Conclusions

We show how comparative transcriptomics can be used to create useful resources for evo-devo, and to investigate phylogenetic relationships. Our results demonstrate that de novo assembly of short (Illumina) reads yields high-quality, high-coverage transcriptomic data sets. We use these data to investigate deep dipteran phylogenetic relationships. Our results, based on a concatenation of 160 orthologous genes, provide support for the traditional view of Clogmia being the sister group of Brachycera (Megaselia, Episyrphus, Drosophila), rather than that of Culicomorpha (which includes mosquitoes and blackflies).

An EST database of the Caribbean fruit fly, Anastrepha suspensa (Diptera: Tephritidae)

Invasive tephritid fruit flies are a great threat to agriculture worldwide and warrant serious pest control measures. Molecular strategies that promote embryonic lethality in these agricultural pests are limited by the small amount of nucleotide sequence data available for tephritids. To increase the dataset for sequence mining, we generated an EST database by 454 sequencing of the caribfly, Anastrepha suspensa, a model tephritid pest. This database yielded 95,803 assembled sequences with 24% identified as independent transcripts. The percentage of caribfly sequences with hits to the closely related tephritid, Rhagoletis pomonella, transcriptome was higher (28%) than to Drosophila proteins/genes (18%) in NCBI. The database contained genes specifically expressed in embryos, genes involved in the cell death, sex-determination, and RNAi pathways, and transposable elements and microsatellites. This study significantly expands the nucleotide data available for caribflies and will be a valuable resource for gene isolation and genomic studies in tephritid insects.

Molecular Techniques for the Detection and Differentiation of Host and Parasitoid Species and the Implications for Fruit Fly Management

Parasitoid detection and identification is a necessary step in the development and implementation of fruit fly biological control strategies employing parasitoid augmentive release. In recent years, DNA-based methods have been used to identify natural enemies of pest species where morphological differentiation is problematic. Molecular techniques also offer a considerable advantage over traditional morphological methods of fruit fly and parasitoid discrimination as well as within-host parasitoid identification, which currently relies on dissection of immature parasitoids from the host, or lengthy and labour-intensive rearing methods. Here we review recent research focusing on the use of molecular strategies for fruit fly and parasitoid detection and differentiation and discuss the implications of these studies on fruit fly management.

Computational strategies for the genome-wide identification of cis-regulatory elements and transcriptional targets

Transcription factors (TFs) are key proteins that decode the information in our genome to express a precise and unique set of proteins and RNA molecules in each cell type in our body. These factors play a pivotal role in all biological processes, including the determination of a cell's fate during development and the maintenance of a cell's physiological function. To achieve this, a TF binds to specific DNA sequences in the noncoding part of the genome, recruits chromatin modifiers and cofactors, and directs the transcription initiation rate of its "target genes." Therefore, a key challenge in deciphering a transcriptional switch is to identify the direct target genes of the master regulators that control the switch, the cis-regulatory elements implementing (auto-)regulatory loops, and the target genes of all the TFs in the downstream regulatory network. A better knowledge of a TF's targetome during specification and differentiation of a particular cell type will generate mechanistic insight into its developmental program. Here, I review computational strategies and methods to predict transcriptional targets by genome-wide searches for TF binding sites using position weight matrices, motif clusters, phylogenetic footprinting, chromatin binding and accessibility data, enhancer classification, motif enrichment, and gene expression signatures.

Robustness and evolvability in the B-system of flower development

BACKGROUND

Gene duplication has often been invoked as a key mechanism responsible for evolution of new morphologies. The floral homeotic B-group gene family has undergone a number of gene duplication events, and yet the functions of these genes appear to be largely conserved. However, detailed comparative analysis has indicated that such duplicate genes have considerable cryptic variability in their functions. In the Solanaceae, two duplicate B-group gene lineages have been retained in three subfamilies. Comparisons of orthologous genes across members of the Solanaceae have demonstrated that the combined function of all four B-gene members is to establish petal and stamen identity, but that this function was partitioned differently in each species. These observations emphasize both the robustness and the evolvability of the B-system.

SCOPE

We provide an overview of how the B-function genes can robustly specify petal and stamen identity and at the same time evolve through changes in protein-protein interaction, gene expression patterns, copy number variation or alterations in the downstream genes they control. By using mathematical models we explore regulatory differences between species and how these impose constraints on downstream gene regulation.

CONCLUSIONS

Evolvability of the B-genes can be understood through the multiple ways in which the B-system can be robust. Quantitative approaches should allow for the incorporation of more biological realism in the representations of these regulatory systems and this should contribute to understanding the constraints under which different B-systems can function and evolve. This, in turn, can provide a better understanding of the ways in which B-genes have contributed to flower diversity.

Changes in selective effects over time facilitate turnover of enhancer sequences

Correct gene expression is often critical and consequently stabilizing selection on expression is widespread. Yet few genes possess highly conserved regulatory DNA, and for the few enhancers that have been carefully characterized, substantial functional reorganization has often occurred. Given that natural selection removes mutations of even very small deleterious effect, how can transcription factor binding evolve so readily when it underlies a conserved phenotype? As a first step toward addressing this question, I combine a computational model for regulatory function that incorporates many aspects of our present biological knowledge with a model for the fitness effects of misexpression. I then use this model to study the evolution of enhancers. Several robust behaviors emerge: First, the selective effects of mutations at a site change dramatically over time due to substitutions elsewhere in the enhancer, and even the overall degree of constraint across the enhancer can change considerably. Second, many of the substitutions responsible for changes in binding occur at sites where previously the mutation would have been strongly deleterious, suggesting that fluctuations in selective effects at a site are important for functional turnover. Third, most substitutions contributing to the repatterning of binding and constraint are effectively neutral, highlighting the importance of genetic drift-even for enhancers underlying conserved phenotypes. These findings have important implications for phylogenetic inference of function and for interpretations of selection coefficients estimated for regulatory DNA.

Hox gene Ultrabithorax regulates distinct sets of target genes at successive stages of Drosophila haltere morphogenesis

Hox genes encode highly conserved transcription factors that regionalize the animal body axis by controlling complex developmental processes. Although they are known to operate in multiple cell types and at different stages, we are still missing the batteries of genes targeted by any one Hox gene over the course of a single developmental process to achieve a particular cell and organ morphology. The transformation of wings into halteres by the Hox gene Ultrabithorax (Ubx) in Drosophila melanogaster presents an excellent model system to study the Hox control of transcriptional networks during successive stages of appendage morphogenesis and cell differentiation. We have used an inducible misexpression system to switch on Ubx in the wing epithelium at successive stages during metamorphosis--in the larva, prepupa, and pupa. We have then used extensive microarray expression profiling and quantitative RT-PCR to identify the primary transcriptional responses to Ubx. We find that Ubx targets range from regulatory genes like transcription factors and signaling components to terminal differentiation genes affecting a broad repertoire of cell behaviors and metabolic reactions. Ubx up- and down-regulates hundreds of downstream genes at each stage, mostly in a subtle manner. Strikingly, our analysis reveals that Ubx target genes are largely distinct at different stages of appendage morphogenesis, suggesting extensive interactions between Hox genes and hormone-controlled regulatory networks to orchestrate complex genetic programs during metamorphosis.

The gap gene network

Gap genes are involved in segment determination during the early development of the fruit fly Drosophila melanogaster as well as in other insects. This review attempts to synthesize the current knowledge of the gap gene network through a comprehensive survey of the experimental literature. I focus on genetic and molecular evidence, which provides us with an almost-complete picture of the regulatory interactions responsible for trunk gap gene expression. I discuss the regulatory mechanisms involved, and highlight the remaining ambiguities and gaps in the evidence. This is followed by a brief discussion of molecular regulatory mechanisms for transcriptional regulation, as well as precision and size-regulation provided by the system. Finally, I discuss evidence on the evolution of gap gene expression from species other than Drosophila. My survey concludes that studies of the gap gene system continue to reveal interesting and important new insights into the role of gene regulatory networks in development and evolution.

Long-range regulatory synergy is required to allow control of the TAC1 locus by MEK/ERK signalling in sensory neurones

Changes in the expression of the neuropeptide substance P (SP) in different populations of sensory neurones are associated with the progression of chronic inflammatory disease. Thus, understanding the genomic and cellular mechanisms driving the expression of the TAC1 gene, which encodes SP, in sensory neurones is essential to understanding its role in inflammatory disease. We used a novel combination of computational genomics, primary-cell culture and mouse transgenics to determine the genomic and cellular mechanisms that control the expression of TAC1 in sensory neurones. Intriguingly, we demonstrated that the promoter of the TAC1 gene must act in synergy with a remote enhancer, identified using comparative genomics, to respond to MAPK signalling that modulates the expression of TAC1 in sensory neurones. We also reveal that noxious stimulation of sensory neurones triggers this synergy in larger diameter sensory neurones – an expression of SP associated with hyperalgesia. This noxious stimulation of TAC1 enhancer-promotor synergy could be strongly blocked by antagonism of the MEK pathway. This study provides a unique insight into the role of long-range enhancer-promoter synergy and selectivity in the tissue-specific response of promoters to specific signal transduction pathways and suggests a possible new avenue for the development of novel anti-inflammatory therapies.

Assessing computational methods of cis-regulatory module prediction

Computational methods attempting to identify instances of cis-regulatory modules (CRMs) in the genome face a challenging problem of searching for potentially interacting transcription factor binding sites while knowledge of the specific interactions involved remains limited. Without a comprehensive comparison of their performance, the reliability and accuracy of these tools remains unclear. Faced with a large number of different tools that address this problem, we summarized and categorized them based on search strategy and input data requirements. Twelve representative methods were chosen and applied to predict CRMs from the Drosophila CRM database REDfly, and across the human ENCODE regions. Our results show that the optimal choice of method varies depending on species and composition of the sequences in question. When discriminating CRMs from non-coding regions, those methods considering evolutionary conservation have a stronger predictive power than methods designed to be run on a single genome. Different CRM representations and search strategies rely on different CRM properties, and different methods can complement one another. For example, some favour homotypical clusters of binding sites, while others perform best on short CRMs. Furthermore, most methods appear to be sensitive to the composition and structure of the genome to which they are applied. We analyze the principal features that distinguish the methods that performed well, identify weaknesses leading to poor performance, and provide a guide for users. We also propose key considerations for the development and evaluation of future CRM-prediction methods.

Transcriptional regulation involves multiple transcription factors binding to DNA sequences. A limited repertoire of transcription factors performs this complex regulatory step through various spatial and temporal interactions between themselves and their binding sites. These transcription factor binding interactions are clustered as distinct modules: cis-regulatory modules (CRMs). Computational methods attempting to identify instances of CRMs in the genome face a challenging problem because a majority of these interactions between transcription factors remain unknown. To investigate the reliability and accuracy of these methods, we chose twelve representative methods and applied them to predict CRMs on both the fly and human genomes. Our results show that the optimal choice of method varies depending on species and composition of the sequences in question. Different CRM representations and search strategies rely on different CRM properties, and different methods can complement one another. We provide a guide for users and key considerations for developers. We also expect that, along with new technology generating new types of genomic data, future CRM prediction methods will be able to reveal transcription binding interactions in three-dimensional space.

Massive turnover of functional sequence in human and other mammalian genomes

Despite the availability of dozens of animal genome sequences, two key questions remain unanswered: First, what fraction of any species' genome confers biological function, and second, are apparent differences in organismal complexity reflected in an objective measure of genomic complexity? Here, we address both questions by applying, across the mammalian phylogeny, an evolutionary model that estimates the amount of functional DNA that is shared between two species' genomes. Our main findings are, first, that as the divergence between mammalian species increases, the predicted amount of pairwise shared functional sequence drops off dramatically. We show by simulations that this is not an artifact of the method, but rather indicates that functional (and mostly noncoding) sequence is turning over at a very high rate. We estimate that between 200 and 300 Mb (∼6.5%-10%) of the human genome is under functional constraint, which includes five to eight times as many constrained noncoding bases than bases that code for protein. In contrast, in D. melanogaster we estimate only 56-66 Mb to be constrained, implying a ratio of noncoding to coding constrained bases of about 2. This suggests that, rather than genome size or protein-coding gene complement, it is the number of functional bases that might best mirror our naïve preconceptions of organismal complexity.

High level of structural polymorphism driven by mobile elements in the Hox genomic region of the Chaetognath Spadella cephaloptera

Little is known about the relationships between genome polymorphism, mobile element dynamics, and population size among animal populations. The chaetognath species Spadella cephaloptera offers a unique perspective to examine this issue because they display a high level of genetic polymorphism at the population level. Here, we have investigated in detail the extent of nucleotide and structural polymorphism in a region harboring Hox1 and several coding genes and presumptive functional elements. Sequencing of several bacterial artificial chromosome inserts representative of this nuclear region uncovered a high level of structural heterogeneity, which is mainly caused by the polymorphic insertion of a diversity of genetic mobile elements. By anchoring this variation through individual genotyping, we demonstrated that sequence diversity could be attributed to the allelic pool of a single population, which was confirmed by detection of extensive recombination within the genomic region studied. The high average level of nucleotide heterozygosity provides clues of selection in both coding and noncoding domains. This pattern stresses how selective processes remarkably cope with intense sequence turnover due to substitutions, mobile element insertions, and recombination to preserve the integrity of functional landscape. These findings suggest that genome polymorphism could provide pivotal information for future functional annotation of genomes.

Sex and the single embryo: early deveiopment in the Mediterranean fruit fly, Ceratitis capitata

Background

In embryos the maternal-to-zygotic transition (MTZ) integrates post-transcriptional regulation of maternal transcripts with transcriptional activation of the zygotic genome. Although the molecular mechanisms underlying this event are being clarified in Drosophila melanogaster, little is know about the embryogenic processes in other insect species. The recent publication of expressed sequence tags (ESTs) from embryos of the global pest species Ceratitis capitata (medfly) has enabled the investigation of embryogenesis in this species and has allowed a comparison of the embryogenic processes in these two related dipteran species, C. capitata and D. melanogaster, that shared a common ancestor 80-100 mya.

Results

Using a novel PCR-based sexing method, which takes advantage of a putative LTR retrotransposon MITE insertion on the medfly Y chromosome, the transcriptomes of individual early male and female embryos were analysed using RT-PCR. This study is focused on two crucial aspects of the onset of embryonic development: sex determination and cellular blastoderm formation. Together with the three known medfly genes (Cctransformer, Cctransformer2 and Ccdoublesex), the expression patterns of other medfly genes that are similar to the D. melanogaster sex-determination genes (sisterlessA, groucho, deadpan, Sex-lethal, female lethal d, sans fille and intersex) and four cellular blastoderm formation genes (Rho1, spaghetti squash, slow-as-molasses and serendipity-α) were analyzed, allowing us to sketch a preliminary outline of the embryonic process in the medfly. Furthermore, a putative homologue of the Zelda gene has been considered, which in D. melanogaster encodes a DNA-binding factor responsible for the maternal-to-zygotic transition.

Conclusions

Our novel sexing method facilitates the study of i) when the MTZ transition occurs in males and females of C. capitata, ii) when and how the maternal information of "female-development" is reprogrammed in the embryos and iii) similarities and differences in the regulation of gene expression in C. capitata and D. melanogaster. We suggest a new model for the onset of the sex determination cascade in the medfly: the maternally inherited Cctra transcripts in the female embryos are insufficient to produce enough active protein to inhibit the male mode of Cctra splicing. The slow rate of development and the inefficiency of the splicing mechanism in the pre-cellular blastoderm facilitates the male-determining factor (M) activity, which probably acts by inhibiting CcTRA protein activity.

Dissecting the regulatory switches of development: lessons from enhancer evolution in Drosophila

Cis-regulatory modules are non-protein-coding regions of DNA essential for the control of gene expression. One class of regulatory modules is embryonic enhancers, which drive gene expression during development as a result of transcription factor protein binding at the enhancer sequences. Recent comparative studies have begun to investigate the evolution of the sequence architecture within enhancers. These analyses are illuminating the way that developmental biologists think about enhancers by revealing their molecular mechanism of function.

Comparative genomics of gene regulation-conservation and divergence of cis-regulatory information

We recently witnessed a tremendous increase in genomics studies on gene regulation and in entirely sequenced genomes from closely related species. This has triggered analyses that suggest a wide range of evolutionary dynamics of gene regulation, from rapid turnover of transcription-factor binding sites to conservation of enhancer function across large evolutionary distances. Many examples show that enhancers can evolve beyond recognizable sequence similarity while retaining function. However, bioinformatics approaches are increasingly able to detect conserved regulatory elements through characteristic evolutionary sequence signatures. Cis-regulatory changes are also a major source of morphological evolution, which might be facilitated by many biochemically functional elements that are selectively neutral and by the buffering function of redundant enhancers and 'shadow' enhancers.

Evolution of transcriptional control in mammals

Changes in gene expression directed by transcriptional regulators can give rise to new phenotypes. While gene expression profiles can be maintained across large evolutionary distances, transcription factor-DNA interactions diverge rapidly. The application of new genome-wide methodologies has begun refining our global understanding of when and where mammalian transcription factors interact with DNA, thereby providing new insight into the mechanisms of transcriptional evolution. The interplay between cis and trans regulation of gene expression is an increasingly active area of investigation, and recent studies suggest that mutations in cis-regulatory DNA can explain many inter-species differences in gene expression.

Functional evolution of cis-regulatory modules at a homeotic gene in Drosophila

It is a long-held belief in evolutionary biology that the rate of molecular evolution for a given DNA sequence is inversely related to the level of functional constraint. This belief holds true for the protein-coding homeotic (Hox) genes originally discovered in Drosophila melanogaster. Expression of the Hox genes in Drosophila embryos is essential for body patterning and is controlled by an extensive array of cis-regulatory modules (CRMs). How the regulatory modules functionally evolve in different species is not clear. A comparison of the CRMs for the Abdominal-B gene from different Drosophila species reveals relatively low levels of overall sequence conservation. However, embryonic enhancer CRMs from other Drosophila species direct transgenic reporter gene expression in the same spatial and temporal patterns during development as their D. melanogaster orthologs. Bioinformatic analysis reveals the presence of short conserved sequences within defined CRMs, representing gap and pair-rule transcription factor binding sites. One predicted binding site for the gap transcription factor KRUPPEL in the IAB5 CRM was found to be altered in Superabdominal (Sab) mutations. In Sab mutant flies, the third abdominal segment is transformed into a copy of the fifth abdominal segment. A model for KRUPPEL-mediated repression at this binding site is presented. These findings challenge our current understanding of the relationship between sequence evolution at the molecular level and functional activity of a CRM. While the overall sequence conservation at Drosophila CRMs is not distinctive from neighboring genomic regions, functionally critical transcription factor binding sites within embryonic enhancer CRMs are highly conserved. These results have implications for understanding mechanisms of gene expression during embryonic development, enhancer function, and the molecular evolution of eukaryotic regulatory modules.

Conserved noncoding elements and the evolution of animal body plans

The genomes of vertebrates, flies, and nematodes contain highly conserved noncoding elements (CNEs). CNEs cluster around genes that regulate development, and where tested, they can act as transcriptional enhancers. Within an animal group CNEs are the most conserved sequences but between groups they are normally diverged beyond recognition. Alternative CNEs are, however, associated with an overlapping set of genes that control development in all animals. Here, we discuss the evidence that CNEs are part of the core gene regulatory networks (GRNs) that specify alternative animal body plans. The major animal groups arose >550 million years ago. We propose that the cis-regulatory inputs identified by CNEs arose during the "re-wiring" of regulatory interactions that occurred during early animal evolution. Consequently, different animal groups, with different core GRNs, contain alternative sets of CNEs. Due to the subsequent stability of animal body plans, these core regulatory sequences have been evolving in parallel under strong purifying selection in different animal groups.

What can you do with 0.1x genome coverage? A case study based on a genome survey of the scuttle fly Megaselia scalaris (Phoridae)

BACKGROUND

The declining cost of DNA sequencing is making genome sequencing a feasible option for more organisms, including many of interest to ecologists and evolutionary biologists. While obtaining high-depth, completely assembled genome sequences for most non-model organisms remains challenging, low-coverage genome survey sequences (GSS) can provide a wealth of biologically useful information at low cost. Here, using a random pyrosequencing approach, we sequence the genome of the scuttle fly Megaselia scalaris and evaluate the utility of our low-coverage GSS approach.

RESULTS

Random pyrosequencing of the M. scalaris genome provided a depth of coverage (0.05x0.1x) much lower than typical GSS studies. We demonstrate that, even with extremely low-coverage sequencing, bioinformatics approaches can yield extensive information about functional and repetitive elements. We also use our GSS data to develop genomic resources such as a nearly complete mitochondrial genome sequence and microsatellite markers for M. scalaris.

CONCLUSION

We conclude that low-coverage genome surveys are effective at generating useful information about organisms currently lacking genomic sequence data.