Unraveling the transcriptional Cis-regulatory code

It is nowadays accepted that eukaryotic complexity is not dictated by the number of protein-coding genes of the genome, but rather achieved through the combinatorics of gene expression programs. Distinct aspects of the expression pattern of a gene are mediated by discrete regulatory sequences, known as cis-regulatory elements. The work described in this thesis was aimed at developing computational and statistical methods to guide the search and characterization of novel cis-regulatory elements

Endemicity analysis, parsimony and biotic elements: A formal comparison using hypothetical distributions

There is as yet no general agreement regarding the most appropriate solution to the problem of identifying areas of endemism, not even in particular cases. In this study, we compared Endemicity Analysis (EA), Parsimony Analysis of Endemicity (PAE), and Biotic Elements Analysis (BE) based on their ability to identify hypothetical predefined patterns that represent nested, overlapping, and disjoint areas of endemism supported by species with different degrees of sympatry. We found that PAE performs poorly when applied to patterns that either overlap with each other or are supported by species with imperfect sympatry. BE exhibits a counterintuitive sensitivity to the degree of congruence among the distributions of endemic species, being unable to recognize areas of endemism supported by perfectly sympatric species. In contrast, in all cases examined we found that EA results in a high proportion of correctly identified distributional patterns. In addition to highlighting the strengths and limitations of these approaches, our results show how different methods can lead to seemingly conflicting conclusions and caution about the possibility of identifying distributional patterns that are merely methodological artefacts.Fil: Casagranda, Maria Dolores. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tucumán; Argentina. Universidad Nacional de Tucumán. Facultad de Ciencias Naturales e Instituto Miguel Lillo. Instituto Superior de Entomología; ArgentinaFil: Taher, Leila Maria. Universidad Nacional de Tucumán. Facultad de Ciencias Naturales e Instituto Miguel Lillo. Instituto Superior de Entomología; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tucumán; ArgentinaFil: Szumik, Claudia Adriana. Universidad Nacional de Tucumán. Facultad de Ciencias Naturales e Instituto Miguel Lillo. Instituto Superior de Entomología; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tucumán; Argentin

Cancer Is Associated with Alterations in the Three-Dimensional Organization of the Genome

The human genome is organized into topologically associating domains (TADs), which represent contiguous regions with a higher frequency of intra-interactions as opposed to inter-interactions. TADs contribute to gene expression regulation by restricting the interactions between their regulatory elements, and TAD disruption has been associated with cancer. Here, we provide a proof of principle that mutations within TADs can be used to predict the survival of cancer patients. Specifically, we constructed a set of 1467 consensus TADs representing the three-dimensional organization of the human genome and used Cox regression analysis to identify a total of 35 prognostic TADs in different cancer types. Interestingly, only 46% of the 35 prognostic TADs comprised genes with known clinical relevance. Moreover, in the vast majority of such cases, the prognostic value of the TAD was not directly related to the presence/absence of mutations in the gene(s), emphasizing the importance of regulatory mutations. In addition, we found that 34% of the prognostic TADs show strong structural perturbations in the cancer genome, consistent with the widespread, global epigenetic dysregulation often observed in cancer patients. In summary, this study elucidates the mechanisms through which non-coding variants may influence cancer progression and opens new avenues for personalized medicine

Sost and its paralog Sostdc1 coordinate digit number in a Gli3-dependent manner.

WNT signaling is critical in most aspects of skeletal development and homeostasis, and antagonists of WNT signaling are emerging as key regulatory proteins with great promise as therapeutic agents for bone disorders. Here we show that Sost and its paralog Sostdc1 emerged through ancestral genome duplication and their expression patterns have diverged to delineate non-overlapping domains in most organ systems including musculoskeletal, cardiovascular, nervous, digestive, reproductive and respiratory. In the developing limb, Sost and Sostdc1 display dynamic expression patterns with Sost being restricted to the distal ectoderm and Sostdc1 to the proximal ectoderm and the mesenchyme. While Sostdc1(-/-) mice lack any obvious limb or skeletal defects, Sost(-/-) mice recapitulate the hand defects described for Sclerosteosis patients. However, elevated WNT signaling in Sost(-/-); Sostdc1(-/-) mice causes misregulation of SHH signaling, ectopic activation of Sox9 in the digit 1 field and preaxial polydactyly in a Gli1- and Gli3-dependent manner. In addition, we show that the syndactyly documented in Sclerosteosis is present in both Sost(-/-) and Sost(-/-); Sostdc1(-/-) mice, and is driven by misregulation of Fgf8 in the AER, a region lacking Sost and Sostdc1 expression. This study highlights the complexity of WNT signaling in skeletal biology and disease and emphasizes how redundant mechanism and non-cell autonomous effects can synergize to unveil new intricate phenotypes caused by elevated WNT signaling

Comparative computational analysis of pluripotency in human and mouse stem cells

Pluripotent cells can be subdivided into two distinct states, the naïve and the primed state, the latter being further advanced on the path of differentiation. There are substantial differences in the regulation of pluripotency between human and mouse, and in humans only stem cells that resemble the primed state in mouse are readily available. Reprogramming of human stem cells into a more naïve-like state is an important research focus. Here, we developed a pipeline to reanalyze transcriptomics data sets that describe both states, naïve and primed pluripotency, in human and mouse. The pipeline consists of identifying regulated start-ups/shut-downs in terms of molecular interactions, followed by functional annotation of the genes involved and aggregation of results across conditions, yielding sets of mechanisms that are consistently regulated in transitions towards similar states of pluripotency. Our results suggest that one published protocol for naïve human cells gave rise to human cells that indeed share putative mechanisms with the prototypical naïve mouse pluripotent cells, such as DNA damage response and histone acetylation. However, cellular response and differentiation-related mechanisms are similar between the naïve human state and the primed mouse state, so the naïve human state did not fully reflect the naïve mouse state

Evolutionary conserved networks of human height identify multiple Mendelian causes of short stature

Height is a heritable and highly heterogeneous trait. Short stature affects 3% of the population and in most cases is genetic in origin. After excluding known causes, 67% of affected individuals remain without diagnosis. To identify novel candidate genes for short stature, we performed exome sequencing in 254 unrelated families with short stature of unknown cause and identified variants in 63 candidate genes in 92 (36%) independent families. Based on systematic characterization of variants and functional analysis including expression in chondrocytes, we classified 13 genes as strong candidates. Whereas variants in at least two families were detected for all 13 candidates, two genes had variants in 6 (UBR4) and 8 (LAMA5) families, respectively. To facilitate their characterization, we established a clustered network of 1025 known growth and short stature genes, which yielded 29 significantly enriched clusters, including skeletal system development, appendage development, metabolic processes, and ciliopathy. Eleven of the candidate genes mapped to 21 of these clusters, including CPZ, EDEM3, FBRS, IFT81, KCND1, PLXNA3, RASA3, SLC7A8, UBR4, USP45, and ZFHX3. Fifty additional growth-related candidates we identified await confirmation in other affected families. Our study identifies Mendelian forms of growth retardation as an important component of idiopathic short stature

CRISPR/Cas9-Mediated Knock-Out of KrasG12D Mutated Pancreatic Cancer Cell Lines

In 90% of pancreatic ductal adenocarcinoma cases, genetic alteration of the proto-oncogene Kras has occurred, leading to uncontrolled proliferation of cancerous cells. Targeting Kras has proven to be difficult and the battle against pancreatic cancer is ongoing. A promising approach to combat cancer was the discovery of the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated (Cas) system, which can be used to genetically modify cells. To assess the potential of a CRISPR/CRISPR-associated protein 9 (Cas9) method to eliminate Kras mutations in cells, we aimed to knock-out the c.35G>A (p.G12D) Kras mutation. Therefore, three cell lines with a heterozygous Kras mutation (the human cell lines SUIT-2 and Panc-1 and the cell line TB32047 from a KPC mouse model) were used. After transfection, puromycin selection and single-cell cloning, proteins from two negative controls and five to seven clones were isolated to verify the knock-out and to analyze changes in key signal transduction proteins. Western blots showed a specific knock-out in the KrasG12D protein, but wildtype Kras was expressed by all of the cells. Signal transduction analysis (for Erk, Akt, Stat3, AMPKα, and c-myc) revealed expression levels similar to the wildtype. The results described herein indicate that knocking-out the KrasG12D mutation by CRISPR/Cas9 is possible. Additionally, under regular growth conditions, the knock-out clones resembled wildtype cells

How enhancers regulate wavelike gene expression patterns

A key problem in development is to understand how genes turn on or off at the right place and right time during embryogenesis. Such decisions are made by non-coding sequences called ‘enhancers.’ Much of our models of how enhancers work rely on the assumption that genes are activated de novo as stable domains across embryonic tissues. Such a view has been strengthened by the intensive landmark studies of the early patterning of the anterior-posterior (AP) axis of the Drosophila embryo, where indeed gene expression domains seem to arise more or less stably. However, careful analysis of gene expression patterns in other model systems (including the AP patterning in vertebrates and short-germ insects like the beetle Tribolium castaneum) painted a different, very dynamic view of gene regulation, where genes are oftentimes expressed in a wavelike fashion. How such gene expression waves are mediated at the enhancer level is so far unclear. Here, we establish the AP patterning of the short-germ beetle Tribolium as a model system to study dynamic and temporal pattern formation at the enhancer level. To that end, we established an enhancer prediction system in Tribolium based on time- and tissue-specific ATAC-seq and an enhancer live reporter system based on MS2 tagging. Using this experimental framework, we discovered several Tribolium enhancers, and assessed the spatiotemporal activities of some of them in live embryos. We found our data consistent with a model in which the timing of gene expression during embryonic pattern formation is mediated by a balancing act between enhancers that induce rapid changes in gene expression patterns (that we call ‘dynamic enhancers’) and enhancers that stabilize gene expression patterns (that we call ‘static enhancers’). However, more data is needed for a strong support for this or any other alternative models

Transcriptional Regulation of Enhancers Active in Protodomains of the Developing Cerebral Cortex

SummaryElucidating the genetic control of cerebral cortical (pallial) development is essential for understanding function, evolution, and disorders of the brain. Transcription factors (TFs) that embryonically regulate pallial regionalization are expressed in gradients, raising the question of how discrete domains are generated. We provide evidence that small enhancer elements active in protodomains integrate broad transcriptional information. CreERT2 and GFP expression from 14 different enhancer elements in stable transgenic mice allowed us to define a comprehensive regional fate map of the pallium. We explored transcriptional mechanisms that control the activity of the enhancers using informatics, in vivo occupancy by TFs that regulate cortical patterning (CoupTFI, Pax6, and Pbx1), and analysis of enhancer activity in Pax6 mutants. Overall, the results provide insights into how broadly expressed patterning TFs regulate the activity of small enhancer elements that drive gene expression in pallial protodomains that fate map to distinct cortical regions