Modeling and simulation of a synthetic genetic circuit that implements a multicelled behavior in a growing microcolony of E. coli

Synthetic Biology is a relatively new discipline, born at the beginning of the New Millennium, that brings the typical engineering approach (abstraction, modularity and standardization) to biotechnology. These principles aim to tame the extreme complexity of the various components and aid the construction of artificial biological systems with specific functions, usually by means of synthetic genetic circuits implemented in bacteria or simple eukaryotes like yeast. The cell becomes a programmable machine and its low-level programming language is made of strings of DNA. This work was performed in collaboration with researchers of the Department of Electrical Engineering of the University of Washington in Seattle and also with a student of the Corso di Laurea Magistrale in Ingegneria Biomedica at the University of Bologna: Marilisa Cortesi. During the collaboration I contributed to a Synthetic Biology project already started in the Klavins Laboratory. In particular, I modeled and subsequently simulated a synthetic genetic circuit that was ideated for the implementation of a multicelled behavior in a growing bacterial microcolony. In the first chapter the foundations of molecular biology are introduced: structure of the nucleic acids, transcription, translation and methods to regulate gene expression. An introduction to Synthetic Biology completes the section. In the second chapter is described the synthetic genetic circuit that was conceived to make spontaneously emerge, from an isogenic microcolony of bacteria, two different groups of cells, termed leaders and followers. The circuit exploits the intrinsic stochasticity of gene expression and intercellular communication via small molecules to break the symmetry in the phenotype of the microcolony. The four modules of the circuit (coin flipper, sender, receiver and follower) and their interactions are then illustrated. In the third chapter is derived the mathematical representation of the various components of the circuit and the several simplifying assumptions are made explicit. Transcription and translation are modeled as a single step and gene expression is function of the intracellular concentration of the various transcription factors that act on the different promoters of the circuit. A list of the various parameters and a justification for their value closes the chapter. In the fourth chapter are described the main characteristics of the gro simulation environment, developed by the Self Organizing Systems Laboratory of the University of Washington. Then, a sensitivity analysis performed to pinpoint the desirable characteristics of the various genetic components is detailed. The sensitivity analysis makes use of a cost function that is based on the fraction of cells in each one of the different possible states at the end of the simulation and the wanted outcome. Thanks to a particular kind of scatter plot, the parameters are ranked. Starting from an initial condition in which all the parameters assume their nominal value, the ranking suggest which parameter to tune in order to reach the goal. Obtaining a microcolony in which almost all the cells are in the follower state and only a few in the leader state seems to be the most difficult task. A small number of leader cells struggle to produce enough signal to turn the rest of the microcolony in the follower state. It is possible to obtain a microcolony in which the majority of cells are followers by increasing as much as possible the production of signal. Reaching the goal of a microcolony that is split in half between leaders and followers is comparatively easy. The best strategy seems to be increasing slightly the production of the enzyme. To end up with a majority of leaders, instead, it is advisable to increase the basal expression of the coin flipper module. At the end of the chapter, a possible future application of the leader election circuit, the spontaneous formation of spatial patterns in a microcolony, is modeled with the finite state machine formalism. The gro simulations provide insights into the genetic components that are needed to implement the behavior. In particular, since both the examples of pattern formation rely on a local version of Leader Election, a short-range communication system is essential. Moreover, new synthetic components that allow to reliably downregulate the growth rate in specific cells without side effects need to be developed. In the appendix are listed the gro code utilized to simulate the model of the circuit, a script in the Python programming language that was used to split the simulations on a Linux cluster and the Matlab code developed to analyze the data.

Molecular strategies for genetic diversity analysis and development of markers linked to resistance traits in apple

The goal of many plant scientists’ research is to explain natural phenotypic variation in term of simple changes in DNA sequence. DNA-based molecular markers are extensively used for the construction of genome-wide molecular maps and to perform genetic analysis for simple and complex traits. The PhD thesis was divided into two main research lines according to the different approaches adopted. The first research line is to analyze the genetic diversity in an Italian apple germplasm collection for the identification of markers tightly linked to targeted genes by an association genetic method. This made it possible to identify synomym and homonym accessions and triploids. The fruit red skin color trait has been used to test the reliability of the genetic approaches in this species. The second line is related to the development of molecular markers closely linked to the Rvi13 and Rvi5 scab resistance genes, previously mapped on apple’s chromosome 10 and 17 respectively by using the traditional linkage mapping method. Both region have been fine-mapped with various type of markers that could be used for marker-assisted selection in future breeding programs and to isolate the two resistance genes.

Spatial Disparity and Dynamic Trajectory of Convergence in China: Construction and Application of a Composite Index of Regional Development

China is a large country characterized by remarkable growth and distinct regional diversity. Spatial disparity has always been a hot issue since China has been struggling to follow a balanced growth path but still confronting with unprecedented pressures and challenges. To better understand the inequality level benchmarking spatial distributions of Chinese provinces and municipalities and estimate dynamic trajectory of sustainable development in China, I constructed the Composite Index of Regional Development (CIRD) with five sub pillars/dimensions involving Macroeconomic Index (MEI), Science and Innovation Index (SCI), Environmental Sustainability Index (ESI), Human Capital Index (HCI) and Public Facilities Index (PFI), endeavoring to cover various fields of regional socioeconomic development. Ranking reports on the five sub dimensions and aggregated CIRD were provided in order to better measure the developmental degrees of 31 or 30 Chinese provinces and municipalities over 13 years from 1998 to 2010 as the time interval of three “Five-year Plans”. Further empirical applications of this CIRD focused on clustering and convergence estimation, attempting to fill up the gap in quantifying the developmental levels of regional comprehensive socioeconomics and estimating the dynamic convergence trajectory of regional sustainable development in a long run. Four clusters were benchmarked geographically-oriented in the map on the basis of cluster analysis, and club-convergence was observed in the Chinese provinces and municipalities based on stochastic kernel density estimation.

Genetic evidence supporting the association of protease and protease inhibitor genes with inflammatory bowel disease: a systematic review

As part of the European research consortium IBDase, we addressed the role of proteases and protease inhibitors (P/PIs) in inflammatory bowel disease (IBD), characterized by chronic mucosal inflammation of the gastrointestinal tract, which affects 2.2 million people in Europe and 1.4 million people in North America. We systematically reviewed all published genetic studies on populations of European ancestry (67 studies on Crohn's disease [CD] and 37 studies on ulcerative colitis [UC]) to identify critical genomic regions associated with IBD. We developed a computer algorithm to map the 807 P/PI genes with exact genomic locations listed in the MEROPS database of peptidases onto these critical regions and to rank P/PI genes according to the accumulated evidence for their association with CD and UC. 82 P/PI genes (75 coding for proteases and 7 coding for protease inhibitors) were retained for CD based on the accumulated evidence. The cylindromatosis/turban tumor syndrome gene (CYLD) on chromosome 16 ranked highest, followed by acylaminoacyl-peptidase (APEH), dystroglycan (DAG1), macrophage-stimulating protein (MST1) and ubiquitin-specific peptidase 4 (USP4), all located on chromosome 3. For UC, 18 P/PI genes were retained (14 proteases and 4 protease inhibitors), with a considerably lower amount of accumulated evidence. The ranking of P/PI genes as established in this systematic review is currently used to guide validation studies of candidate P/PI genes, and their functional characterization in interdisciplinary mechanistic studies in vitro and in vivo as part of IBDase. The approach used here overcomes some of the problems encountered when subjectively selecting genes for further evaluation and could be applied to any complex disease and gene family.

Design and construction of a new Drosophila species, D.synthetica, by synthetic regulatory evolution

Here, I merge the principles of synthetic biology1,2 and regulatory evolution3-11 to create a new species12-15 with a minimal set of known elements. Using preexisting transgenes and recessive mutations of Drosophila melanogaster, a transgenic population arises with small eyes and a different venation pattern that fulfills the criteria of a new species according to Mayr's "Biological Species Concept"7,10. The genetic circuit entails the loss of a non-essential transcription factor and the introduction of cryptic enhancers. Subsequent activation of those enhancers causes hybrid lethality. The transition from "transgenic organisms" towards "synthetic species", such as Drosophila synthetica, constitutes a safety mechanism to avoid hybridization with wild type populations and preserve natural biodiversity16-18. Drosophila synthetica is the first transgenic organism that cannot hybridize with the original wild type population but remains fertile when crossed with other transgenic animals.

A second generation radiation hybrid map to aid the assembly of the bovine genome sequence

BACKGROUND: Several approaches can be used to determine the order of loci on chromosomes and hence develop maps of the genome. However, all mapping approaches are prone to errors either arising from technical deficiencies or lack of statistical support to distinguish between alternative orders of loci. The accuracy of the genome maps could be improved, in principle, if information from different sources was combined to produce integrated maps. The publicly available bovine genomic sequence assembly with 6x coverage (Btau_2.0) is based on whole genome shotgun sequence data and limited mapping data however, it is recognised that this assembly is a draft that contains errors. Correcting the sequence assembly requires extensive additional mapping information to improve the reliability of the ordering of sequence scaffolds on chromosomes. The radiation hybrid (RH) map described here has been contributed to the international sequencing project to aid this process. RESULTS: An RH map for the 30 bovine chromosomes is presented. The map was built using the Roslin 3000-rad RH panel (BovGen RH map) and contains 3966 markers including 2473 new loci in addition to 262 amplified fragment-length polymorphisms (AFLP) and 1231 markers previously published with the first generation RH map. Sequences of the mapped loci were aligned with published bovine genome maps to identify inconsistencies. In addition to differences in the order of loci, several cases were observed where the chromosomal assignment of loci differed between maps. All the chromosome maps were aligned with the current 6x bovine assembly (Btau_2.0) and 2898 loci were unambiguously located in the bovine sequence. The order of loci on the RH map for BTA 5, 7, 16, 22, 25 and 29 differed substantially from the assembled bovine sequence. From the 2898 loci unambiguously identified in the bovine sequence assembly, 131 mapped to different chromosomes in the BovGen RH map. CONCLUSION: Alignment of the BovGen RH map with other published RH and genetic maps showed higher consistency in marker order and chromosome assignment than with the current 6x sequence assembly. This suggests that the bovine sequence assembly could be significantly improved by incorporating additional independent mapping information.

Genetic analysis of syndactyly in German Holstein cattle

Congenital syndactyly with a variable number of affected feet was observed in eight black and white German Holstein calves. Analysis of the pedigree data revealed that all affected individuals could be traced back to a single founder. The pedigree was consistent with monogenic autosomal recessive inheritance and variable expressivity. Bovine syndactyly or "mulefoot" has been previously shown to map on the telomeric end of bovine chromosome 15 and we performed PCR genotyping of microsatellite markers spanning 27 cM of this chromosomal region to test the new cases for genetic linkage with the phenotype. The haplotype segregation confirmed the suggested inheritance pattern of the mulefoot mutation in this family and markers RM004, BM848 and BMS820 showed significant linkage to the phenotype. The results confirmed the chromosomal location of the mulefoot gene in this pedigree. Furthermore the study demonstrated that although marker testing has been available for nearly a decade the use of mulefoot carriers in cattle breeding remains uncontrolled. The presented family provides a resource for positional cloning of the causative mutation.

Genetic diversity of the cestode Echinococcus multilocularis in red foxes at a continental scale in Europe

BACKGROUND: Alveolar echinococcosis (AE) is a severe helminth disease affecting humans, which is caused by the fox tapeworm Echinococcus multilocularis. AE represents a serious public health issue in larger regions of China, Siberia, and other regions in Asia. In Europe, a significant increase in prevalence since the 1990s is not only affecting the historically documented endemic area north of the Alps but more recently also neighbouring regions previously not known to be endemic. The genetic diversity of the parasite population and respective distribution in Europe have now been investigated in view of generating a fine-tuned map of parasite variants occurring in Europe. This approach may serve as a model to study the parasite at a worldwide level. METHODOLOGY/PRINCIPAL FINDINGS: The genetic diversity of E. multilocularis was assessed based upon the tandemly repeated microsatellite marker EmsB in association with matching fox host geographical positions. Our study demonstrated a higher genetic diversity in the endemic areas north of the Alps when compared to other areas. CONCLUSIONS/SIGNIFICANCE: The study of the spatial distribution of E. multilocularis in Europe, based on 32 genetic clusters, suggests that Europe can be considered as a unique global focus of E. multilocularis, which can be schematically drawn as a central core located in Switzerland and Jura Swabe flanked by neighbouring regions where the parasite exhibits a lower genetic diversity. The transmission of the parasite into peripheral regions is governed by a "mainland-island" system. Moreover, the presence of similar genetic profiles in both zones indicated a founder event.

The EmsB tandemly repeated multilocus microsatellite: a new tool to investigate genetic diversity of Echinococcus granulosus sensu lato

Cystic echinococcosis (CE) is a widespread and severe zoonotic disease caused by infection with the larval stage of the eucestode Echinococcus granulosus sensu lato. The polymorphism exhibited by nuclear and mitochondrial markers conventionally used for the genotyping of different parasite species and strains does not reach the level necessary for the identification of genetic variants linked to restricted geographical areas. EmsB is a tandemly repeated multilocus microsatellite that proved its usefulness for the study of genetic polymorphisms within the species E. multilocularis, the causative agent of alveolar echinococcosis. In the present study, EmsB was used to characterize E. granulosus sensu lato samples collected from different host species (sheep, cattle, dromedaries, dogs, and human patients) originating from six different countries (Algeria, Mauritania, Romania, Serbia, Brazil, and the People's Republic of China). The conventional mitochondrial cox1 and nad1 markers identified genotypes G1, G3, G5, G6, and G7, which are clustered into three groups corresponding to the species E. granulosus sensu stricto, E. ortleppi, and E. canadensis. With the same samples, EmsB provided a higher degree of genetic discrimination and identified variations that correlated with the relatively small-scale geographic origins of the samples. In addition, one of the Brazilian single hydatid cysts presented a hybrid genotypic profile that suggested genetic exchanges between E. granulosus sensu stricto and E. ortleppi. In summary, the EmsB microsatellite exhibits an interesting potential for the elaboration of a detailed map of the distribution of genetic variants and therefore for the determination and tracking of the source of CE.

Construction of improved temperature-sensitive and mobilizable vectors and their use for constructing mutations in the adhesin-encoding acm gene of poorly transformable clinical Enterococcus faecium strains.

Inactivation by allelic exchange in clinical isolates of the emerging nosocomial pathogen Enterococcus faecium has been hindered by lack of efficient tools, and, in this study, transformation of clinical isolates was found to be particularly problematic. For this reason, a vector for allelic replacement (pTEX5500ts) was constructed that includes (i) the pWV01-based gram-positive repAts replication region, which is known to confer a high degree of temperature intolerance, (ii) Escherichia coli oriR from pUC18, (iii) two extended multiple-cloning sites located upstream and downstream of one of the marker genes for efficient cloning of flanking regions for double-crossover mutagenesis, (iv) transcriptional terminator sites to terminate undesired readthrough, and (v) a synthetic extended promoter region containing the cat gene for allelic exchange and a high-level gentamicin resistance gene, aph(2'')-Id, to distinguish double-crossover recombination, both of which are functional in gram-positive and gram-negative backgrounds. To demonstrate the functionality of this vector, the vector was used to construct an acm (encoding an adhesin to collagen from E. faecium) deletion mutant of a poorly transformable multidrug-resistant E. faecium endocarditis isolate, TX0082. The acm-deleted strain, TX6051 (TX0082Deltaacm), was shown to lack Acm on its surface, which resulted in the abolishment of the collagen adherence phenotype observed in TX0082. A mobilizable derivative (pTEX5501ts) that contains oriT of Tn916 to facilitate conjugative transfer from the transformable E. faecalis strain JH2Sm::Tn916 to E. faecium was also constructed. Using this vector, the acm gene of a nonelectroporable E. faecium wound isolate was successfully interrupted. Thus, pTEX5500ts and its mobilizable derivative demonstrated their roles as important tools by helping to create the first reported allelic replacement in E. faecium; the constructed this acm deletion mutant will be useful for assessing the role of acm in E. faecium pathogenesis using animal models.

Mapping and cloning of genes from portions of the human genome using somatic cell hybrids

Human x rodent somatic cell hybrids have played an important role in human genetics research. They have been especially useful for assigning genes to chromosomes and isolating DNA markers from specific regions of the human genome.^ By employing a combination of somatic cell genetic, recombinant DNA, and cytogenetic techniques, human DNA excision repair gene ERCC4 was mapped regionally to human 16p13.13-13.2, even though the gene has not been cloned. Human x Chinese hamster ovary (CHO) cell hybrids selected for human ERCC4 activity and containing 16p13.1-p13.3 as the only human genetic material were identified. These hybrids were used to order DNA markers located in 16p13.1-p13.3. New DNA markers physically close to ERCC4 were isolated from such hybrids. Using amplified human DNA from the hybrids as probe in fluorescent in situ hybridization, the short arm breakpoint in the chromosome 16 inversion associated with acute myelomonocytic leukemia (AMML) was found to be physically close to the ERCC4 gene. The physical mapping and eventually, the cloning of the ERCC4 gene, will benefit the understanding of the DNA repair system and the study of other important biomedical problems such as tumorigenesis.^ To facilitate the cloning of ERCC4 gene and, in general, the cloning of genes from any defined regions of the human genome, a method was developed for the direct isolation of human transcribed genes ffom somatic cell hybrids. cDNA was prepared from human x rodent hybrid by using consensus 5$\sp\prime$ splice site sequences as primers. These primers were designed to select immature, unspliced messenger RNA (still retaining species specific repeat sequences) as templates. Screening of a derived cDNA library for human repeat sequences resulted in the isolation of human clones at the anticipated frequency with characteristics expected of exons of transcribed human genes. The usefulness of the splice site specific primers was analyzed and the cDNA synthesis conditions with these primers were optimized. The procedure was shown to be sensitive enough to clone weakly expressed genes. Studying the expression of the represented genes with the isolated clones was shown to be feasible. Such regional specific human gene fragments will be very valuable for many human genetic studies such as the search of inherited disease genes and the construction of a cDNA map of the human genome. ^

Genetic Control of Plant Development by Overriding a Geometric Division Rule

Formative cell divisions are critical for multicellular patterning. In the early plant embryo, such divisions follow from orienting the division plane. A major unanswered question is how division plane orientation is genetically controlled, and in particular whether this relates to cell geometry. We have generated a complete 4D map of early Arabidopsis embryogenesis and used computational analysis to demonstrate that several divisions follow a rule that uses the smallest wall area going through the center of the cell. In other cases, however, cell division clearly deviates from this rule, which invariably leads to asymmetric cell division. By analyzing mutant embryos and through targeted genetic perturbation, we show that response to the hormone auxin triggers a deviation from the ``shortest wall'' rule. Our work demonstrates that a simple default rule couples division orientation to cell geometry in the embryo and that genetic regulation can create patterns by overriding the default rule.

Comparative genetic analyses point to HCP5 as susceptibility locus for HCV-associated hepatocellular carcinoma

BACKGROUND & AIMS: Recently, genetic variations in MICA (lead single nucleotide polymorphism [SNP] rs2596542) were identified by a genome-wide association study (GWAS) to be associated with hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC) in Japanese patients. In the present study, we sought to determine whether this SNP is predictive of HCC development in the Caucasian population as well. METHODS: An extended region around rs2596542 was genotyped in 1924 HCV-infected patients from the Swiss Hepatitis C Cohort Study (SCCS). Pair-wise correlation between key SNPs was calculated both in the Japanese and European populations (HapMap3: CEU and JPT). RESULTS: To our surprise, the minor allele A of rs2596542 in proximity of MICA appeared to have a protective impact on HCC development in Caucasians, which represents an inverse association as compared to the one observed in the Japanese population. Detailed fine-mapping analyses revealed a new SNP in HCP5 (rs2244546) upstream of MICA as strong predictor of HCV-related HCC in the SCCS (univariable p=0.027; multivariable p=0.0002, odds ratio=3.96, 95% confidence interval=1.90-8.27). This newly identified SNP had a similarly directed effect on HCC in both Caucasian and Japanese populations, suggesting that rs2244546 may better tag a putative true variant than the originally identified SNPs. CONCLUSIONS: Our data confirms the MICA/HCP5 region as susceptibility locus for HCV-related HCC and identifies rs2244546 in HCP5 as a novel tagging SNP. In addition, our data exemplify the need for conducting meta-analyses of cohorts of different ethnicities in order to fine map GWAS signals.

Developmental self-construction and -configuration of functional neocortical neuronal networks.

The prenatal development of neural circuits must provide sufficient configuration to support at least a set of core postnatal behaviors. Although knowledge of various genetic and cellular aspects of development is accumulating rapidly, there is less systematic understanding of how these various processes play together in order to construct such functional networks. Here we make some steps toward such understanding by demonstrating through detailed simulations how a competitive co-operative ('winner-take-all', WTA) network architecture can arise by development from a single precursor cell. This precursor is granted a simplified gene regulatory network that directs cell mitosis, differentiation, migration, neurite outgrowth and synaptogenesis. Once initial axonal connection patterns are established, their synaptic weights undergo homeostatic unsupervised learning that is shaped by wave-like input patterns. We demonstrate how this autonomous genetically directed developmental sequence can give rise to self-calibrated WTA networks, and compare our simulation results with biological data.

Map positions of third chromosomal female sterile and lethal mutations of Drosophila melanogaster

Chromosomal mutations induced by ethyl methanesulfonate (EMS) treatment can cause female sterility or maternal-effect lethality in Drosophila. EMS is particularly useful to researchers because it creates mutations independent of position effects. However, because researchers have little control over the chromosomal site of mutation, post-mutagenic genetic mapping is required to determine the cytological location of the mutation. To make a valuable set of mutants more useful to the research community, we have mapped the uncharacterized part of the female-sterile – maternal-effect lethal Tübingen collection. We mapped 49 female-sterile – maternal-effect lethal alleles and 72 lethal alleles to individual deficiency intervals on the third chromosome. In addition, we analyzed the phenotype of ovaries resulting from female sterile mutations. The observed phenotypes range from tumorous ovaries and early blocks in oogenesis, to later blocks, slow growth, blocks in stage 10, to apparently full development of the ovary. The mapping and phenotypic characterization of these 121 mutations provide the necessary information for the researcher to consider a specific mutant as a candidate for their gene of interest.Key words: Drosophila melanogaster, oogenesis, female sterile, maternal-effect lethal, EMS-induced mutations.