Evolution of lineage-specific functions in ancient cis-regulatory modules

Morphological evolution is driven both by coding sequence variation and by changes in regulatory sequences. However, how cis-regulatory modules (CRMs) evolve to generate entirely novel expression domains is largely unknown. Here, we reconstruct the evolutionary history of a lens enhancer located within a CRM that not only predates the lens, a vertebrate innovation, but bilaterian animals in general. Alignments of orthologous sequences from different deuterostomes sub-divide the CRM into a deeply conserved core and a more divergent flanking region. We demonstrate that all deuterostome flanking regions, including invertebrate sequences, activate gene expression in the zebrafish lens through the same ancient cluster of activator sites. However, levels of gene expression vary between species due to the presence of repressor motifs in flanking region and core. These repressor motifs are responsible for the relatively weak enhancer activity of tetrapod flanking regions. Ray-finned fish, however, have gained two additional lineage-specific activator motifs which in combination with the ancient cluster of activators and the core constitute a potent lens enhancer. The exploitation and modification of existing regulatory potential in flanking regions but not in the highly conserved core might represent a more general model for the emergence of novel regulatory functions in complex CRMs.     

The structure and evolution of cis-regulatory regions: the shavenbaby story

In this paper, we provide a historical account of the contribution of a single line of research to our current understanding of the structure of cis-regulatory regions and the genetic basis for morphological evolution. We revisit the experiments that shed light on the evolution of larval cuticular patterns within the genus Drosophila and the evolution and structure of the shavenbaby gene. We describe the experiments that led to the discovery that multiple genetic changes in the cis-regulatory region of shavenbaby caused the loss of dorsal cuticular hairs (quaternary trichomes) in first instar larvae of Drosophila sechellia. We also discuss the experiments that showed that the convergent loss of quaternary trichomes in D. sechellia and Drosophila ezoana was generated by parallel genetic changes in orthologous enhancers of shavenbaby. We discuss the observation that multiple shavenbaby enhancers drive overlapping patterns of expression in the embryo and that these apparently redundant enhancers ensure robust shavenbaby expression and trichome morphogenesis under stressful conditions. All together, these data, collected over 13 years, provide a fundamental case study in the fields of gene regulation and morphological evolution, and highlight the importance of prolonged, detailed studies of single genes.     

Developmental cis-regulatory analysis of the cyclin D gene in the sea urchin Strongylocentrotus purpuratus

Cyclin D genes regulate the cell cycle, growth and differentiation in response to intercellular signaling. While the promoters of vertebrate cyclin D genes have been analyzed, the cis-regulatory sequences across an entire cyclin D locus have not. Doing so would increase understanding of how cyclin D genes respond to the regulatory states established by developmental gene regulatory networks, linking cell cycle and growth control to the ontogenetic program. Therefore, we conducted a cis-regulatory analysis on the cyclin D gene, SpcycD, of the sea urchin, Strongylocentrotus purpuratus, during embryogenesis, identifying upstream and intronic sequences, located within six defined regions bearing one or more cis-regulatory modules each.     

R11: a cis-regulatory node of the sea urchin embryo gene network that controls early expression of SpDelta in micromeres

A gene regulatory network (GRN) controls the process by which the endomesoderm of the sea urchin embryo is specified. In this GRN, the program of gene expression unique to the skeletogenic micromere lineage is set in train by activation of the pmar1 gene. Through a double repression system, this gene is responsible for localization of expression of downstream regulatory and signaling genes to cells of this lineage. One of these genes, delta, encodes a Notch ligand, and its expression in the right place and time is crucial to the specification of the endomesoderm. Here we report a cis-regulatory element R11 that is responsible for localizing the expression of delta by means of its response to the pmar1 repression system. R11 was identified as an evolutionarily conserved genomic sequence located about 13 kb downstream of the last exon of the delta gene. We demonstrate here that this cis-regulatory element is able to drive the expression of a reporter gene in the same cells and at the same time that the endogenous delta gene is expressed, and that temporally, spatially, and quantitatively it responds to the pmar1 repression system just as predicted for the delta gene in the endomesoderm GRN. This work illustrates the application of cis-regulatory analysis to the validation of predictions of the GRN model. In addition, we introduce new methodological tools for quantitative measurement of the output of expression constructs that promise to be of general value for cis-regulatory analysis in sea urchin embryos.     

Effects of gene regulatory reprogramming on gene expression in human and mouse developing hearts

Lineage-specific regulatory elements underlie adaptation of species and play a role in disease susceptibility. We compared functionally conserved and lineage-specific enhancers by cross-mapping 5042 human and 6564 mouse heart enhancers. Of these, 79 per cent are lineage-specific, lacking a functional orthologue. Heart enhancers tend to cluster and, commonly, there are multiple heart enhancers in a heart locus providing a regulatory stability to the locus. We observed little cross-clustering, however, between lineage-specific and functionally conserved heart enhancers suggesting regulatory function acquisition and development in loci previously lacking heart activity. We also identified 862 human-specific heart enhancers: 417 featuring sequence conservation with mouse (class II) and 445 with neither sequence nor function conservation (class III). Ninety-eight per cent of class III enhancers were deleted from the mouse genome, and we estimated a similar-sized enhancer gain in the human lineage. Human-specific enhancers display no detectable decrease in the negative selection pressure and are strongly associated with genes partaking in the heart regulatory programmes. The loss of a heart enhancer could be compensated by activity of a redundant heart enhancer; however, we observed redundancy in only 15 per cent of class II and III enhancer loci indicating a large-scale reprogramming of the heart regulatory programme in mammals.     

A complex ensemble of cis-regulatory elements controls the expression of a Vicia faba non-storage seed protein gene

We have identified cis-regulatory elements within the 5-upstream region of a Vicia faba non-storage seed protein gene, called usp, by studying the expression of usp-promoter deletion fragments fused to reporter genes in transgenic tobacco seeds. 0.4 kb of usp upstream sequence contain at least six, but probably more, distinct cis-regulatory elements which are responsible for seemingly all quantitative, spatial and temporal aspects of expression. Expression-increasing and-decreasing elements are interspersed and include an AT-rich sequence, a G-box element and a CATGCATG motif. The latter acts as a negative element in contrast to what has been found for the same motif in legumin-and vicilin-type seed storage protein gene promoters. Seed specificity of expression is mainly determined by the –68/+51 region which confers, however, only very low levels of expression. The data support the combinatiorial model of promoter function.     

Global spread and genetic variants of the two CYP9M10 haplotype forms associated with insecticide resistance in Culex quinquefasciatus Say

Insecticide resistance develops as a genetic factor (allele) conferring lower susceptibility to insecticides proliferates within a target insect population under strong positive selection. Intriguingly, a resistance allele pre-existing in a population often bears a series of further adaptive allelic variants through new mutations. This phenomenon occasionally results in replacement of the predominating resistance allele by fitter new derivatives, and consequently, development of greater resistance at the population level. The overexpression of the cytochrome P450 gene CYP9M10 is associated with pyrethroid resistance in the southern house mosquito Culex quinquefasciatus. Previously, we have found two genealogically related overexpressing CYP9M10 haplotypes, which differ in gene copy number (duplicated and non-duplicated). The duplicated haplotype was derived from the non-duplicated overproducer probably recently. In the present study, we investigated allelic series of CYP9M10 involved in three C. quinquefasciatus laboratory colonies recently collected from three different localities. Duplicated and non-duplicated overproducing haplotypes coexisted in African and Asian colonies indicating a global distribution of both haplotype lineages. The duplicated haplotypes both in the Asian and African colonies were associated with higher expression levels and stronger resistance than non-duplicated overproducing haplotypes. There were slight variation in expression level among the non-duplicated overproducing haplotypes. The nucleotide sequences in coding and upstream regions among members of this group also showed a little diversity. Non-duplicated overproducing haplotypes with relatively higher expression were genealogically closer to the duplicated haplotypes than the other non-duplicated overproducing haplotypes, suggesting multiple cis-acting mutations before duplication.     

Confocal quantification of cis-regulatory reporter gene expression in living sea urchin

Quantification of GFP reporter gene expression at single cell level in living sea urchin embryos can now be accomplished by a new method of confocal laser scanning microscopy (CLSM). Eggs injected with a tissue-specific GFP reporter DNA construct were grown to gastrula stage and their fluorescence recorded as a series of contiguous Z-section slices that spanned the entire embryo. To measure the depth-dependent signal decay seen in the successive slices of an image stack, the eggs were coinjected with a freely diffusible internal fluorescent standard, rhodamine dextran. The measured rhodamine fluorescence was used to generate a computational correction for the depth-dependent loss of GFP fluorescence per slice. The intensity of GFP fluorescence was converted to the number of GFP molecules using a conversion constant derived from CLSM imaging of eggs injected with a measured quantity of GFP protein. The outcome is a validated method for accurately counting GFP molecules in given cells in reporter gene transfer experiments, as we demonstrate by use of an expression construct expressed exclusively in skeletogenic cells.     

The genetic basis of evolutionary change in gene expression levels

The regulation of gene expression is an important determinant of organismal phenotype and evolution. However, the widespread recognition of this fact occurred long after the synthesis of evolution and genetics. Here, we give a brief sketch of thoughts regarding gene regulation in the history of evolution and genetics. We then review the development of genome-wide studies of gene regulatory variation in the context of the location and mode of action of the causative genetic changes. In particular, we review mapping of the genetic basis of expression variation through expression quantitative trait locus studies and measuring the cis/trans component of expression variation in allele-specific expression studies. We conclude by proposing a systematic integration of ideas that combines global mapping studies, cis/trans tests and modern population genetics methodologies, in order to directly estimate the forces acting on regulatory variation within and between species.     

A cis-encoded sRNA controls the expression of fabH2 in Yersinia

YsrH is a novel cis-encoded sRNA located on the opposite strand to fabH2, which is essential for fatty acid biosynthesis in bacteria. In this study, YsrH-mediated regulation of fabH2 expression was investigated in Yersinia pseudotuberculosis. Constitutive and inducible over-expression of YsrH decreased the mRNA level of fabH2, while expression of downstream fabD and fabG remained unaffected. Polynucleotide phosphorylase (PNPase) also played an important role in this regulation process by mediating YsrH decay in the exponential phase. Thus, our data defines a cis-encoded sRNA that regulates fatty acid synthesis via a regulatory mechanism also involving PNPase.     

Conservation and diversity in the cis-regulatory networks that integrate information controlling expression of Hoxa2 in hindbrain and cranial neural crest cells in vertebrates

The Hoxa2 and Hoxb2 genes are members of paralogy group II and display segmental patterns of expression in the developing vertebrate hindbrain and cranial neural crest cells. Functional analyses have demonstrated that these genes play critical roles in regulating morphogenetic pathways that direct the regional identity and anteroposterior character of hindbrain rhombomeres and neural crest-derived structures. Transgenic regulatory studies have also begun to characterize enhancers and cis-elements for those mouse and chicken genes that direct restricted patterns of expression in the hindbrain and neural crest. In light of the conserved role of Hoxa2 in neural crest patterning in vertebrates and the similarities between paralogs, it is important to understand the extent to which common regulatory networks and elements have been preserved between species and between paralogs. To investigate this problem, we have cloned and sequenced the intergenic region between Hoxa2 and Hoxa3 in the chick HoxA complex and used it for making comparative analyses with the respective human, mouse, and horn shark regions. We have also used transgenic assays in mouse and chick embryos to test the functional activity of Hoxa2 enhancers in heterologous species. Our analysis reveals that three of the critical individual components of the Hoxa2 enhancer region from mouse necessary for hindbrain expression (Krox20, BoxA, and TCT motifs) have been partially conserved. However, their number and organization are highly varied for the same gene in different species and between paralogs within a species. Other essential mouse elements appear to have diverged or are absent in chick and shark. We find the mouse r3/r5 enhancer fails to work in chick embryos and the chick enhancer works poorly in mice. This implies that new motifs have been recruited or utilized to mediate restricted activity of the enhancer in other species. With respect to neural crest regulation, cis-components are embedded among the hindbrain control elements and are highly diverged between species. Hence, there has been no widespread conservation of sequence identity over the entire enhancer domain from shark to humans, despite the common function of these genes in head patterning. This provides insight into how apparently equivalent regulatory regions from the same gene in different species have evolved different components to potentiate their activity in combination with a selection of core components.     

Overall and allele-specific expression of the SMC1A gene in female Cornelia de Lange syndrome patients and healthy controls

Cornelia de Lange syndrome (CdLS) is a rare multisystem disorder characterized by facial dysmorphisms, limb anomalies, and growth and cognitive deficits. Mutations in genes encoding subunits (SMC1A, SMC3, RAD21) or regulators (NIPBL, HDAC8) of the cohesin complex account for approximately 65% of clinically diagnosed CdLS cases. The SMC1A gene (Xp11.22), responsible for 5% of CdLS cases, partially escapes X chromosome inactivation in humans and the allele on the inactive X chromosome is variably expressed. In this study, we evaluated overall and allele-specific SMC1A expression. Real-time PCR analysis conducted on 17 controls showed that SMC1A expression in females is 50% higher than in males. Immunoblotting experiments confirmed a 44% higher protein level in healthy females than in males, and showed no significant differences in SMC1A protein levels between controls and patients. Pyrosequencing was used to assess the reciprocal level of allelic expression in six female carriers of different SMC1A mutations and 15 controls who were heterozygous at a polymorphic transcribed SMC1A locus. The two alleles were expressed at a 1:1 ratio in the control group and at a 2:1 ratio in favor of the wild type allele in the test group. Since a dominant negative effect is considered the pathogenic mechanism in SMC1A-defective female patients, the level of allelic preferential expression might be one of the factors contributing to the wide phenotypic variability observed in these patients. An extension of this study to a larger cohort containing mild to borderline cases could enhance our understanding of the clinical spectrum of SMC1A-linked CdLS.     

HSM — a hybrid system based approach for modelling intracellular networks

The paper proposes a hybrid system based approach for modelling of intracellular networks and introduces a restricted subclass of hybrid systems – HSM – with an objective of still being able to provide sufficient power for the modelling of biological systems, while imposing some restrictions that facilitate analysis of systems described by such models.