Quadruped Gait Transitions from a Neural Pattern Generator with Arousal Modulated Interactions

A four-channel neural pattern generator is described in which both the frequency and the relative phase of oscillations are controlled by a scalar arousal or GO signal. The generator is used to simulate quadruped gaits; in particular, rapid transitions are simulated in the order - walk, trot, pace, and gallop - that occurs in the cat. Precise switching control is achieved by using an arousal dependent modulation of the model's inhibitory interactions. This modulation generates a different functional connectivity in a single network at different arousal levels.

Novel insights into the role of the GABAB receptor in depression and anxiety

The GABAB receptor is a functional heterodimer comprising the GABAB1 and GABAB2 subunits, with the GABAB1 subunit displaying two major isoforms, GABAB(1a) and GABAB(1b). Preclinical findings have strongly implicated the GABAB receptor in stress-related psychiatric disorders, however, the precise contribution of the GABAB receptor in depression and anxiety disorders remains unknown. Emerging data suggest that the interaction between adverse environmental conditions, such as early life stress, and a specific genetic composition can increase the risk to develop psychiatric disorders in adulthood. This thesis investigated the role of the GABAB receptor alone or in combination with early-life stress (maternal separation), in modulating antidepressant like and anxiety-related behaviours. Pharmacological blockade of the GABAB receptor with CGP52432 had antidepressant-like behavioural effects. Moreover, mice lacking the GABAB(1b) receptor subunit isoform exhibited antidepressant-like behaviours in adulthood but anxiety-like behaviour in early-life. In response to maternal separation, GABAB(1a)-/- mice exhibited early-life stress-induced anhedonia, a core symptom of depression, while GABAB(1b)-/- mice exhibited a more resilient phenotype. Moreover, when compared with wildtype or GABAB(1a)-/- mice, GABAB(1b)-/- mice that underwent maternal separation exhibited enhanced stressinduced neuronal activation in the hippocampus and in the nucleus accumbens (NAcc), a critical area for anhedonia thus suggesting that enhanced stress-induced neuronal activation in the hippocampus and NAcc in GABAB(1b)-/- mice may be important for their antidepressant-like phenotype and their resilience to stress-induced anhedonia. Pharmacological blockade of GABAB receptor and GABAB(1b) receptor subunit isoform loss of function increased adult hippocampal cell proliferation, thus suggesting that increased hippocampal neurogenesis could be a potential mechanism for the antidepressant-like effects of GABAB receptor antagonists and GABAB(1b) receptor subunit isoform disruption. Finally, this thesis investigated whether the expression of several genes involved in hippocampal neurogenesis or the antidepressant response were altered in the mouse hippocampus following chronic treatment with a GABAB receptor antagonist.

Copper-catalysed asymmetric sulfide oxidation

The focus of this thesis is the preparation of enantiopure sulfoxides by means of copper-catalysed asymmetric sulfoxidation, with particular emphasis on the synthesis of aryl benzyl and aryl alkyl sulfoxides. Chapter 1 contains a review of the methods employed for the asymmetric synthesis of sulfoxides, compounds with many applications in stereoselective synthesis and in some cases with pharmaceutical application. Chapter 1 describes asymmetric oxidation, including metal-catalysed, non metal-catalysed and enzyme-catalysed, in addition to synthetic approaches via nucleophilic substitution of appropriately substituted precursors. Kinetic resolution in oxidation of sulfoxides to the analogous sulfones is also discussed; in certain cases, access to enantioenriched sulfoxides can be achieved via a combination of asymmetric sulfoxidation and complementary kinetic resolution. The design and synthesis of a series of sulfides to enable exploration of the substituent effects of the copper-mediated oxidation was undertaken, and oxidation to the racemic sulfoxides and sulfones to provide reference samples was conducted. Oxidation of the sulfides using copper-Schiff base catalysis was undertaken leading to enantioenriched sulfoxides. The procedure employed is clean, inexpensive, not air-sensitive and utilises aqueous hydrogen peroxide as oxidant. Extensive investigation of the influence of the reaction conditions such as solvent, temperature, copper salt and ligand was undertaken to lead to the optimised conditions. While the direct attachment of one aryl substituent to the sulfide is essential for efficient enantiocontrol, in the case of the second substituent the enantiocontol is dependent on the steric rather than electronic features of the substituent. Significantly, use of naphthyl-substituted sulfides results in excellent enantiocontrol; notably 97% ee, obtained in the oxidation of 2-naphthyl benzyl sulfide, represents the highest enantioselectivity reported to date for a copper-mediated sulfur oxidation. Some insight into the mechanistic features of the copper-mediated sulfur oxidation has been developed based on this work, although further investigation is required to establish the precise nature of the catalytic species responsible for asymmetric sulfur oxidation. Full experimental details, describing the synthesis and structural characterisation, and determination of enantiopurity are included in chapter 3.

The role of metabolism in the pathogenesis of adherent invasive Escherichia coli (AIEC)

Crohn’s disease (CD) is a chronic, relapsing inflammatory condition affecting the gastrointestinal tract of humans, of which there is currently no cure. The precise etiology of CD is unknown, although it has become widely accepted that it is a multifactorial disease which occurs as a result of an abnormal immune response to commensal enteric bacteria in a genetically susceptible host. Recent studies have shown that a new group of Escherichia coli, called Adherent Invasive Escherichia coli (AIEC) are present in the guts of CD patients at a higher frequency than in healthy subjects, suggesting that they may play a role in the initiation and/or maintenance of the inflammation associated with CD. Two phenotypes define an AIEC and differentiate them from other groups of E. coli. Firstly, AIEC can adhere to and invade epithelial cells; and secondly, they can replicate in macrophages. In this study, we use a strain of AIEC which has been isolated from the colonic mucosa of a CD patient, called HM605, to examine the relationship between AIEC and the macrophage. We show, using a systematic mutational approach, that while the Tricarboxylic acid (TCA) cycle, the glyoxylate pathway, the Entner-Doudoroff (ED) pathway, the Pentose Phosphate (PP) pathway and gluconeogenesis are dispensable for the intramacrophagic growth of HM605, glycolysis is an absolute requirement for the ability of this organism to replicate intracellularly. We also show that HM605 activates the inflammasome, a major driver of inflammation in mammals. Specifically, we show that macrophages infected with HM605 produce significantly higher levels of the pro-inflammatory cytokine IL-1β than macrophages infected with a non-AIEC strain, and we show by immunoblotting that this is due to cleavage of caspase-1. We also show that macrophages infected with HM605 exhibit significantly higher levels of pyroptosis, a form of inflammatory cell death, than macrophages infected with a non-AIEC strain. Therefore, AIEC strains are more pro-inflammatory than non-AIEC strains and this may have important consequences in terms of CD pathology. Moreover, we show that while inflammasome activation by HM605 is independent of intracellular bacterial replication, it is dependent on an active bacterial metabolism. Through the establishment of a genetic screen aimed at identifying mutants which activate the inflammasome to lower levels than the wild-type, we confirm our observation that bacterial metabolism is essential for successful inflammasome activation by HM605 and we also uncover new systems/structures that may be important for inflammasome activation, such as the BasS/BasR two-component system, a type VI secretion system and a K1 capsule. Finally, in this study, we also identify a putative small RNA in AIEC strain LF82, which may be involved in modulating the motility of this strain. Overall this works shows that, in the absence of specialised virulence factors, the ability of AIEC to metabolise within the host cell may be a key determinant of its pathogenesis.

Development and characterization of novel transgenic techniques for the study of neuronal circuitry

Modern neuroscience relies heavily on sophisticated tools that allow us to visualize and manipulate cells with precise spatial and temporal control. Transgenic mouse models, for example, can be used to manipulate cellular activity in order to draw conclusions about the molecular events responsible for the development, maintenance and refinement of healthy and/or diseased neuronal circuits. Although it is fairly well established that circuits respond to activity-dependent competition between neurons, we have yet to understand either the mechanisms underlying these events or the higher-order plasticity that synchronizes entire circuits. In this thesis we aimed to develop and characterize transgenic mouse models that can be used to directly address these outstanding biological questions in different ways. We present SLICK-H, a Cre-expressing mouse line that can achieve drug-inducible, widespread, neuron-specific manipulations in vivo. This model is a clear improvement over existing models because of its particularly strong, widespread, and even distribution pattern that can be tightly controlled in the absence of drug induction. We also present SLICK-V::Ptox, a mouse line that, through expression of the tetanus toxin light chain, allows long-term inhibition of neurotransmission in a small subset (<1%) of fluorescently labeled pyramidal cells. This model, which can be used to study how a silenced cell performs in a wildtype environment, greatly facilitates the in vivo study of activity-dependent competition in the mammalian brain. As an initial application we used this model to show that tetanus toxin-expressing CA1 neurons experience a 15% - 19% decrease in apical dendritic spine density. Finally, we also describe the attempt to create additional Cre-driven mouse lines that would allow conditional alteration of neuronal activity either by hyperpolarization or inhibition of neurotransmission. Overall, the models characterized in this thesis expand upon the wealth of tools available that aim to dissect neuronal circuitry by genetically manipulating neurons in vivo.

An evaluation of orogenic kinematic evolution utilizing crystalline and magnetic anisotropy in granitoids

The Silurian-Devonian Galway Granite Complex (GGC ~425-380Ma) is defined here as a suite of granitoid plutons that comprise the Main Galway Granite Batholith and the Earlier Plutons. The Main Batholith is a composite of the Carna Pluton in the west and the Kilkieran Pluton in the east and extends from Galway City ~130km to the west. The Earlier Plutons are spatially, temporally and structurally distinct, situated northwest of the Main Batholith and include the Roundstone, Omey, Inis and Letterfrack Plutons. The majority of isotopic and structural data currently available pertain to the Kilkieran Pluton, several tectonic models have already been devised for this part of the complex. These relate emplacement of the Kilkieran Pluton to extension across a large east-west Caledonian lineament, i.e. the Skird Rocks Fault, during late Caledonian transtension. No chronological data have been published that directly and accurately date the emplacement of the Carna Pluton or any of the Earlier Plutons. There is also a lack of data pertaining to the internal structure of these intrusions. Accordingly, no previous study has established the mechanisms of emplacement for the Earlier Plutons and only limited work is available for the Carna Pluton. As a consequense of this, constituents of the GGC have not previously been placed in a context relative to each other or to regional scale Silurio-Devonian kinematics. The current work focuses on the Omey, Roundstone and Carna Plutons. Here, results of detailed field and Anisotropy of Magnetic Susceptibiliy (AMS) fabric studies are presented. This work is complemented by geological mapping that focuses on fault dynamics and contact relationships. Interpretation of AMS data is aided by rock magnetic experiment data and petrographic microstructural evaluations of representative samples. A new geological map of the the Omey Pluton demonstrates that this intrusion has a defined roof and base which are gently inclined parallel to the fold hinge of the Connemara Antiform. AMS and petrographic data show the intrusion is cross cut by NNW-SSE shear zones that extend into the country rock. These pre-date and were active during magma emplacement. It is proposed that the Omey pluton was emplaced as a discordant phacolith. Pre-existing subvertical D5 faults in the host rock were reactived during emplacement, due to regional sinistral transpression, and served as centralised ascent conduits. A central portion of the Roundstone Pluton was mapped in detail for the first time. Two facies are identified, G1 forms the majority of the pluton and coeval G2 sheets cross cut G1 at the core of the pluton. NNW-SSE D5 faults mapped in the country rock extend across the pluton. These share a geometrical relationship with the distribution of submagmatic strain in the pluton and parallel the majoity of mapped subvertical G2 dykes. These data indicate that magma ascent was controlled by NNW-SSE conduits that are inherently related to those identifed in the Omey Pluton. It is proposed that the Roundstone Pluton is a punched laccolith, the symmetry and structure of which was controlled by pre-exising host rock structures and regional sinistral transpressive stress which presided during emplacement. Field relationships show the long axis of the Carna Pluton lies parallel to mulitple NNW-SSE shear zones. These are represented on a regional scale by the Clifden-Mace Fault which cross cuts the core of this intrusion. AMS and petrographic data show concentric emplacement fabrics were tectonically overprinted as magma cooled from the magmatic state due to this faulting. It is proposed that the Clifden-Mace Fault system was active during ascent and emplacement of the magma and that pluton inflation only terminated as this controlling structure went into compression due to the onset of regional transtension. U-Pb zircon laser ablation inductively coupled mass spectrometry (LA-ICP-MS) data has been compiled from four sample sites. New geochronological data from the Roundstone Pluton (RD1 = ± 3.2Ma) represent the oldest age determination obtained from any member of the GGC and demonstrates that this pluton predates the Carna Pluton by ~10Ma and probably intruded synchronously with the Omey Pluton (~422.5 ± 1.7Ma). Chronological data from the Carna Pluton (CN2 = 412.9 ± 2.5Ma; CN3 = 409.8 ± 7.2Ma; CN4 = 409.6 ± 3.6Ma) represent the first precise magma crystallisation age for this intrusion. This work shows this pluton is 10Ma older than the Kilkieran Pluton and that the supply of magma into the Carna Pluton had terminated by ~409Ma. Chronological, magnetic and field data have been utilised to evaluate the kinematic evolution of the Caledonides of western Ireland throughout the construction of the GGC. It is proposed that the GGC was constructed during four distinct episodes. The style of emplacement and the conduits used for magma transport to the site of emplacement was dependent on the orientation of local structures relative to the regional ambiant stress field. This philosophy is used to critically evaluate and progress existing hypotheses on the transition from regional transpression to regional transtension at the end of the Caledonian Orogeny.

Fundamental studies on the application of enzymes when brewing with unmalted oats and sorghum

The use of unmalted oats or sorghum in brewing has great potential for creating new beer types/flavors and saving costs. However, the substitution of barley malt with oat or sorghum adjunct is not only innovative but also challenging due to their specific grain characteristics. The overall objectives of this Ph.D. project were: 1) to investigate the impact of various types and levels of oats or sorghum on the quality/processability of mashes, worts, and beers; 2) to provide solutions as regards the application of industrial enzymes to overcome potential brewing problems. For these purposes, a highly precise rheological method using a controlled stress rheometer was developed and successfully applied as a tool for optimizing enzyme additions and process parameters. Further, eight different oat cultivars were compared in terms of their suitability as brewing adjuncts and two very promising types identified. In another study, the limitations of barley malt enzymes and the benefits of the application of industrial enzymes in high-gravity brewing with oats were determined. It is recommended to add enzymes to high-gravity mashes when substituting 30% or more barley malt with oats in order to prevent filtration and fermentation problems. Pilot-scale brewing trials using 10–40% unmalted oats revealed that the sensory quality of oat beers improved with increasing adjunct level. In addition, commercially available oat and sorghum flours were implemented into brewing. The use of up to 70% oat flour and 50% sorghum flour, respectively, is not only technically feasible but also economically beneficial. In a further study on sorghum was demonstrated that the optimization of industrial mashing enzymes has great potential for reducing beer production costs. A comparison of the brewing performance of red Italian and white Nigerian sorghum clearly showed that European grown sorghum is suitable for brewing purposes; 40% red sorghum beers were even found to be very low in gluten.

Langmuir-Blodgett assembly of colloidal crystals combined with controlled infiltration of conducting metal oxides

Colloidal photonic crystals (PhCs) possess a periodic dielectric structure which gives rise to a photonic band gap (PBG) and offer great potential in the ability to modify or control light at visible wavelengths. Although the refractive index contrast between the void or infill and the matrix material is paramount for photonics applications, integration into real optoelectronics devices will require a range of added functionalities such as conductivity. As such, colloidal PhCs can be used as templates to direct infiltration of other functional materials using a range of deposition strategies. The work in this thesis seeks to address two challenges; first to develop a reproducible strategy based on Langmuir-Blodgett (LB) deposition to assemble high quality colloidal PhCs based on silica with precise film thickness as most other assembly methods suffer from a lack of reproducibility thickness control. The second is to investigate the use of LBdeposited colloidal PhCs as templates for infiltration with conducting metal oxide materials using vapor phase deposition techniques. Part of this work describes the synthesis and assembly of colloidal silica spheres with different surface chemical functionalities at the air-water interface in preparation for LB deposition. Modification of surface funtionality conferred varying levels of hydrophobicity upon the particles. The behaviour of silica monolayer films at the air-water interface was characterised by Brewster Angle Microscopy and surface pressure isotherms with a view to optimising the parameters for LB deposition of multilayer colloidal PhC films. Optical characterisation of LB-fabricated colloidal PhCs indicated high quality photonic behaviour, exhibiting a pseudo PBG with a sharp Bragg diffraction peak in the visible region and reflectance intensities greater than 60%. Finally the atomic layer deposition (ALD) of nominally undoped ZnO and aluminium “doped” ZnO (Al-doped ZnO) inside the pores of a colloidal PhC assembled by the LB technique was carried out. ALD growth in this study was performed using trimethyl aluminium (TMA) and water as precursors for the alumina and diethyl zinc (DEZn) and water for the ZnO. The ZnO:Al films were grown in a laminate mode, where DEZn pulses were substituted for TMA pulses in the sequences with a Zn:Al ratio 19:1. The ALD growth of ZnO and ZnO:Al in colloidal PhCs was shown to be highly conformal, tuneable and reproducible whilst maintaining excellent photonic character. Furthermore, at high levels of infiltration the opal composite films demonstrated significant conductivity.

Improved kinematic sensing for motion control applications

New compensation methods are presented that can greatly reduce the slit errors (i.e. transition location errors) and interval errors induced due to non-idealities in optical incremental encoders (square-wave). An M/T-type, constant sample-time digital tachometer (CSDT) is selected for measuring the velocity of the sensor drives. Using this data, three encoder compensation techniques (two pseudoinverse based methods and an iterative method) are presented that improve velocity measurement accuracy. The methods do not require precise knowledge of shaft velocity. During the initial learning stage of the compensation algorithm (possibly performed in-situ), slit errors/interval errors are calculated through pseudoinversebased solutions of simple approximate linear equations, which can provide fast solutions, or an iterative method that requires very little memory storage. Subsequent operation of the motion system utilizes adjusted slit positions for more accurate velocity calculation. In the theoretical analysis of the compensation of encoder errors, encoder error sources such as random electrical noise and error in estimated reference velocity are considered. Initially, the proposed learning compensation techniques are validated by implementing the algorithms in MATLAB software, showing a 95% to 99% improvement in velocity measurement. However, it is also observed that the efficiency of the algorithm decreases with the higher presence of non-repetitive random noise and/or with the errors in reference velocity calculations. The performance improvement in velocity measurement is also demonstrated experimentally using motor-drive systems, each of which includes a field-programmable gate array (FPGA) for CSDT counting/timing purposes, and a digital-signal-processor (DSP). Results from open-loop velocity measurement and closed-loop servocontrol applications, on three optical incremental square-wave encoders and two motor drives, are compiled. While implementing these algorithms experimentally on different drives (with and without a flywheel) and on encoders of different resolutions, slit error reductions of 60% to 86% are obtained (typically approximately 80%).

Déficits immunitaires primaires: diagnostic, prise en charge, quelques perspectives.

Severe primary immunodeficiencies (PID) are rare; their global incidence is comparable to that of childhood leukemia; they include more than 100 different entities. Clinical manifestations are: unusually severe or frequent infections or infections that do not respond to adequate treatment; an increased risk of certain malignancies; sometimes auto-immune manifestations. Delayed diagnosis and management of PID can lead to severe and irreversible complications or to death. PID can become manifest only in the adult; in common variable immune deficiency, the median age at diagnosis is between the 2nd and the 3rd decade of life. PID are often transmitted genetically; recent progresses in molecular biology have allowed more precise and earlier, including antenatal, diagnosis. Molecular treatment of 3 infants with a severe immunodeficiency has recently been achieved in April 2000. Those progresses were mostly based on the study of immunodeficiency databases. We present here the work of a Belgian group specialized in PID; meetings have started in June 1997. This group establishes guidelines for the diagnosis and treatment of PID, adapted to the local situation. The elaboration of a national register of PID is also underway; this has to provide all guaranties of anonymity to patients and families. Such a register already exists at the European level; it has provided the basis for new diagnostic and therapeutic possibilities. The inclusion of Belgian data in this register should allow essential progresses essential for our patients.

Leveraging nanoscale plasmonic modes to achieve reproducible enhancement of light.

The strongly enhanced and localized optical fields that occur within the gaps between metallic nanostructures can be leveraged for a wide range of functionality in nanophotonic and optical metamaterial applications. Here, we introduce a means of precise control over these nanoscale gaps through the application of a molecular spacer layer that is self-assembled onto a gold film, upon which gold nanoparticles (NPs) are deposited electrostatically. Simulations using a three-dimensional finite element model and measurements from single NPs confirm that the gaps formed by this process, between the NP and the gold film, are highly reproducible transducers of surface-enhanced resonant Raman scattering. With a spacer layer of roughly 1.6 nm, all NPs exhibit a strong Raman signal that decays rapidly as the spacer layer is increased.

Promoting developmental transcription.

Animal growth and development depend on the precise control of gene expression at the level of transcription. A central role in the regulation of developmental transcription is attributed to transcription factors that bind DNA enhancer elements, which are often located far from gene transcription start sites. Here, we review recent studies that have uncovered significant regulatory functions in developmental transcription for the TFIID basal transcription factors and for the DNA core promoter elements that are located close to transcription start sites.

BMP signaling in the development of the mouse esophagus and forestomach.

The stratification and differentiation of the epidermis are known to involve the precise control of multiple signaling pathways. By contrast, little is known about the development of the mouse esophagus and forestomach, which are composed of a stratified squamous epithelium. Based on prior work in the skin, we hypothesized that bone morphogenetic protein (BMP) signaling is a central player. To test this hypothesis, we first used a BMP reporter mouse line harboring a BRE-lacZ allele, along with in situ hybridization to localize transcripts for BMP signaling components, including various antagonists. We then exploited a Shh-Cre allele that drives recombination in the embryonic foregut epithelium to generate gain- or loss-of-function models for the Bmpr1a (Alk3) receptor. In gain-of-function (Shh-Cre;Rosa26(CAG-loxpstoploxp-caBmprIa)) embryos, high levels of ectopic BMP signaling stall the transition from simple columnar to multilayered undifferentiated epithelium in the esophagus and forestomach. In loss-of-function experiments, conditional deletion of the BMP receptor in Shh-Cre;Bmpr1a(flox/flox) embryos allows the formation of a multilayered squamous epithelium but this fails to differentiate, as shown by the absence of expression of the suprabasal markers loricrin and involucrin. Together, these findings suggest multiple roles for BMP signaling in the developing esophagus and forestomach.

Expression in aneuploid Drosophila S2 cells.

Extensive departures from balanced gene dose in aneuploids are highly deleterious. However, we know very little about the relationship between gene copy number and expression in aneuploid cells. We determined copy number and transcript abundance (expression) genome-wide in Drosophila S2 cells by DNA-Seq and RNA-Seq. We found that S2 cells are aneuploid for >43 Mb of the genome, primarily in the range of one to five copies, and show a male genotype ( approximately two X chromosomes and four sets of autosomes, or 2X;4A). Both X chromosomes and autosomes showed expression dosage compensation. X chromosome expression was elevated in a fixed-fold manner regardless of actual gene dose. In engineering terms, the system "anticipates" the perturbation caused by X dose, rather than responding to an error caused by the perturbation. This feed-forward regulation resulted in precise dosage compensation only when X dose was half of the autosome dose. Insufficient compensation occurred at lower X chromosome dose and excessive expression occurred at higher doses. RNAi knockdown of the Male Specific Lethal complex abolished feed-forward regulation. Both autosome and X chromosome genes show Male Specific Lethal-independent compensation that fits a first order dose-response curve. Our data indicate that expression dosage compensation dampens the effect of altered DNA copy number genome-wide. For the X chromosome, compensation includes fixed and dose-dependent components.

Modeling the evolution of regulatory elements by simultaneous detection and alignment with phylogenetic pair HMMs.

The computational detection of regulatory elements in DNA is a difficult but important problem impacting our progress in understanding the complex nature of eukaryotic gene regulation. Attempts to utilize cross-species conservation for this task have been hampered both by evolutionary changes of functional sites and poor performance of general-purpose alignment programs when applied to non-coding sequence. We describe a new and flexible framework for modeling binding site evolution in multiple related genomes, based on phylogenetic pair hidden Markov models which explicitly model the gain and loss of binding sites along a phylogeny. We demonstrate the value of this framework for both the alignment of regulatory regions and the inference of precise binding-site locations within those regions. As the underlying formalism is a stochastic, generative model, it can also be used to simulate the evolution of regulatory elements. Our implementation is scalable in terms of numbers of species and sequence lengths and can produce alignments and binding-site predictions with accuracy rivaling or exceeding current systems that specialize in only alignment or only binding-site prediction. We demonstrate the validity and power of various model components on extensive simulations of realistic sequence data and apply a specific model to study Drosophila enhancers in as many as ten related genomes and in the presence of gain and loss of binding sites. Different models and modeling assumptions can be easily specified, thus providing an invaluable tool for the exploration of biological hypotheses that can drive improvements in our understanding of the mechanisms and evolution of gene regulation.