Targeted Deletion Of Functionally Validated Enhancers Defines Their Role in Mouse Limb Development

Transcriptional enhancers are genomic DNA sequences that contain clustered transcription factor (TF) binding sites. When combinations of TFs bind to enhancer sequences they act together with basal transcriptional machinery to regulate the timing, location and quantity of gene transcription. Elucidating the genetic mechanisms responsible for differential gene expression, including the role of enhancers, during embryological and postnatal development is essential to an understanding of evolutionary processes and disease etiology. Numerous methods are in use to identify and characterize enhancers. Several high-throughput methods generate large datasets of enhancer sequences with putative roles in embryonic development. However, few enhancers have been deleted from the genome to determine their roles in the development of specific structures, such as the limb. Manipulation of enhancers at their endogenous loci, such as the deletion of such elements, leads to a better understanding of the regulatory interactions, rules and complexities that contribute to faithful and variant gene transcription – the molecular genetic substrate of evolution and disease. To understand the endogenous roles of two distinct enhancers known to be active in the mouse embryo limb bud we deleted them from the mouse genome. I hypothesized that deletion of these enhancers would lead to aberrant limb development. The enhancers were selected because of their association with p300, a protein associated with active transcription, and because the human enhancer sequences drive distinct lacZ expression patterns in limb buds of embryonic day (E) 11.5 transgenic mice. To confirm that the orthologous mouse enhancers, mouse 280 and 1442 (M280 and M1442, respectively), regulate expression in the developing limb we generated stable transgenic lines, and examined lacZ expression. In M280-lacZ mice, expression was detected in E11.5 fore- and hindlimbs in a region that corresponds to digits II-IV. M1442-lacZ mice exhibited lacZ expression in posterior and anterior margins of the fore- and hindlimbs that overlapped with digits I and V and several wrist bones. We generated mice lacking the M280 and M1442 enhancers by gene targeting. Intercrosses between M280 -/+ and M1442 -/+, respectively, generated M280 and M1442 null mice, which are born at expected Mendelian ratios and manifest no gross limb malformations. Quantitative real-time PCR of mutant E11.5 limb buds indicated that significant changes in transcriptional output of enhancer-proximal genes accompanied the deletion of both M280 and M1442. In neonatal null mice we observed that all limb bones are present in their expected positions, an observation also confirmed by histology of E18.5 distal limbs. Fine-scale measurement of E18.5 digit bone lengths found no differences between mutant and control embryos. Furthermore, when the developmental progression of cartilaginous elements was analyzed in M280 and M1442 embryos from E13.5-E15.5, transient development defects were not detected. These results demonstrate that M280 and M1442 are not required for mouse limb development. Though M280 is not required for embryonic limb development it is required for the development and/or maintenance of body size – adult M280 mice are significantly smaller than control littermates. These studies highlight the importance of experiments that manipulate enhancers in situ to understand their contribution to development.

Functional studies of a LIM -homeodomain transcription factor lmx1b in murine mid -hindbrain and limb development

This thesis is centered on applying molecular genetics to study pattern formation during animal development. More specifically, this thesis describes the functional studies of a LIM-homeodomain gene called lmx1b during murine embryogenesis. Lmx1b expression is restricted to the mid-hindbrain junction as well as to the dorsal mesenchyme of the limb, suggesting important functions during mid-hindbrain and limb development. To test these possibilities, lmx1b homozygous mutant mice were generated and their limb and CNS phenotypes examined. Lmx1b homozygous mutant mice exhibit a large reduction of mid-hindbrain structures, and that their limbs are symmetrical along the dorsal-ventral axis as the result of a dorsal to ventral transformation. Taken together, these studies define essential functions for lmx1b in mid-hindbrain patteming and in dorsal limb cell fate determination. However, the molecular mechanisms which accounts for these phenotypes are unknown, and whether lmx1b has same or distinctive functions during the mid-hindbrain and limb development is also unclear. ^ Recently, insight into molecular mechanisms of mid-hindbrain patterning and limb development has resulted from the identification of several factors with restricted expression patterns within these regions. These include the secreted factors wnt-1, fgf-8, wnt-7a and the transcription factors pax-2, and en-1. Targeted disruption of any of these genes in mice suggests that these genes might be involved in similar regulatory pathways. Analysis of the expression of these genes in lmx1b mutants demonstrates that lmxlb is not required for the initiation, but is required to maintain their expression at the mid-hindbrain junction. Thus, lmxlb is not required for specifying mid-hindbrain cell fates, rather, it functions to ensure the establishment or maintenance of a proper organizing center at the mid-hindbrain junction. Interestingly, lmxlb functions cell non-autonomously in chimera analysis, which indicates that lmx1b might regulate the expression of secreted factors such as wnt-1 and/or fgf-8 in the organizing center. In contrast, lmx1b functions cell autonomously in the dorsal limb to govern dorsal ventral limb development and its expression is regulated by with wnt-7a and en-1. However, single and double mutant analysis suggest that all three genes have partially overlapping functions as well as independent functions. The results point toward a complicated network of cross-talks among all three limb axes. ^

ΔNp63 REGULATES A COMPLEX NETWORK OF TARGET GENES IN LIMB AND EPIDERMAL DEVELOPMENT

The skin is composed of two major compartments, the dermis and epidermis. The epidermis forms a barrier to protect the body. The stratified epithelium has self-renewing capacity throughout life, and continuous turnover is mediated by stem cells in the basal layer. p63 is structurally and functionally related to p53. In spite of their structural similarities, p63 is critical for the development and maintenance of stratified epithelial tissues, unlike p53. p63 is highly expressed in the epidermis and previously has been shown to play a critical role in the development and maintenance of the epidermis. The study of p63 has been complicated due to the existence of multiple isoforms: those with a transactivation domain (TAp63) and those lacking this domain (ΔNp63). Mice lacking p63 cannot form skin, have craniofacial and skeletal defects and die within hours after birth. These defects are due to the ability of p63 to regulate multiple processes in skin development including epithelial stem cell proliferation, differentiation, and adherence programs. To determine the roles of these isoforms in skin development and maintenance, isoform specific p63 conditional knock out mice were generated by our lab. TAp63-/- mice age prematurely, develop blisters, and display wound-healing defects that result from hyperproliferation of dermal stem cells. That results in premature depletion of these cells, which are necessary for wound repair, that indicates TAp63 plays a role in dermal/epidermal maintenance. To study the role of ΔNp63, I generated a ΔNp63-/- mouse and analyzed the skin by performing immunofluorescence for markers of epithelial differentiation. The ΔNp63-/- mice developed a thin, disorganized epithelium but differentiation markers were expressed. Interestingly, the epidermis from ΔNp63-/- mice co-expressed K14 and K10 in the same cell suggesting defects in epidermal differentiation and stratification. This phenotype is reminiscent of the DGCR8fl/fl;K14Cre and Dicerfl/fl;K14Cre mice skin. Importantly, DGCR8-/- embryonic stem cells (ESCs) display a hyperproliferation defect by failure to silence pluripotency genes. Furthermore, I have observed that epidermal cells lacking ΔNp63 display a phenotype reminiscent of embryonic stem cells instead of keratinocytes. Thus, I hypothesize that genes involved in maintaining pluripotency, like Oct4, may be upregulated in the absence of ΔNp63. To test this, q-RT PCR was performed for Oct4 mRNA with wild type and ΔNp63-/- 18.5dpc embryo skin. I found that the level of Oct4 was dramatically increased in the absence of ΔNp63-/-. Based on these results, I hypothesized that ΔNp63 induces differentiation by silencing pluripotency regulators, Oct4, Sox2 and Nanog directly through the regulation of DGCR8. I found that DGCR8 restoration resulted in repression of Oct4, Sox2 and Nanog in ΔNp63-/- epidermal cells and rescue differentiation defects. Loss of ΔNp63 resulted in pluripotency that caused defect in proper differentiation and stem cell like phenotype. This led me to culture the ΔNp63-/- epidermal cells in neuronal cell culture media in order to address whether restoration of DGCR8 can transform epidermal cells to neuronal cells. I found that DGCR8 restoration resulted in a change in cell fate. I also found that miR470 and miR145 play a role in the induction of pluripotency by repressing Oct4, Sox2 and Nanog. This indicates that ΔNp63 induces terminal differentiation through the regulation of DGCR8.

Coordination of murine limb muscle, skeleton and connective tissue development by Lmx1b

The underlying genetic defects of a congenital disease Nail-Patella Syndrome are loss-of-function mutations in the LMX1B gene. Lmx1b encodes a LIM-homeodomain transcription factor that is expressed specifically in the dorsal limb bud mesenchyme. Gain- and loss-of-function experiments suggest that Lmx1b is both necessary and sufficient to specify dorsal limb patterning. However, how Lmx1b coordinates patterning of the dorsal tissues in the limb, including muscle, skeleton and connective tissues, remains unknown. One possibility is that each tissue specifies its own pattern cell-autonomously, i.e., Lmx1b is expressed in tissues in which it functions and different tissues do not communicate with each other. Another possibility is that tissues that express Lmx1b interact with adjacent tissues and provide patterning information thereby directing the development of tissues non-cell-autonomously. Previous results showed that Lmx1b is expressed in limb connective tissue and skeleton, but is not expressed in muscle tissue. Moreover, muscles and muscle connective tissue are closely associated during development. Therefore, we hypothesize that Lmx1b controls limb muscle dorsal-ventral (DV) patterning through muscle connective tissue, but regulates skeleton and tendon/ligament development cell-autonomously. ^ To test this hypothesis, we first examined when and where the limb dorsal-ventral asymmetry is established during development. Subsequently, conditional knockout and overexpression experiments were performed to delete or activate Lmx1b in different tissues within the limb. Our results show that deletion of Lmx1b from whole limb mesenchyme results in all dorsal tissues, including muscle, tendon/ligament and skeleton, transforming into ventral structures. Skeleton-specific knockout of Lmx1b led to the dorsal duplication of distal sesamoid and metacarpal bones, but did not affect the pattern formation of other tissues, suggesting that Lmx1b controls skeleton development cell-autonomously. In addition, this skeleton-specific pattern alteration only occurs in distal limb tissues, not proximal limb tissues, indicating different regulatory mechanisms operate along the limb proximal-distal axis. Moreover, skeleton-specific ectopic expression of Lmx1b reveals a complementary skeletal-specific dorsalized phenotype. This result supports a cell-autonomous role for Lmx1b in dorsal-ventral skeletal patterning. This study enriched our understanding of limb development, and the insights from this research may also be applicable for the development of other organs. ^

Comparative studies of Hox genes in mammalian limb and vertebral pattern formation

Divergence of anterior-posterior (AP) limb pattern and differences in vertebral column morphology are the two main examples of mammalian evolution. The Hox genes (homeobox containing gene) have been implicated in driving evolution of these structures. However, regarding Hox genes, how they contribute to the generation of mammalian morphological diversities, is still unclear. Implementing comparative gene expression and phenotypic rescue studies for different mammalian Hox genes could aid in unraveling this mystery. In the first part of this thesis, the expression pattern of Hoxd13 gene, a key Hox gene in the establishment of the limb AP pattern, was examined in developing limbs of bats and mice. Bat forelimbs exhibit a pronounced asymmetric AP pattern and offer a good model to study the molecular mechanisms that contribute to the variety of mammalian limbs. The data showed that the expression domain of bat Hoxd13 was shifted prior to the asymmetric limb plate expansion, whereas its domain in mice was much more symmetric. This finding reveals a correlation between the divergence of Hoxd13 expression and the AP patterning difference in limb development. The second part of this thesis details a phenotypic rescue approach by human HOXB1-9 transgenes in mice with Hoxb1-9 deletion, The mouse mutants displayed homeosis in cervical and anterior thoracic vertebrae. The human transgenes entirely rescued the mouse mutants, suggesting that these human HOX genes have similar functions to their mouse orthologues in anterior axial skeletal patterning. The anterior expressing human HOXB transgenes such as HOXB1-3 were expressed in the mouse embryonic trunk in a similar manner as their murine orthologues. However, the anterior boundary of human HOXB9 expression domain was more posterior than that of the mouse Hoxb9 by 2-3 somites. These data provide the molecular support for the hypothesis that Hox genes are responsible for maintaining similar anterior axial skeletal architectures cervical and anterior thoracic regions, but different architectures in lumbar and posterior thoracic regions between humans and mice. ^

Analysis of Variation in Clubfoot Candidate Genes

Isolated clubfoot, a common birth defect occurring in more than 135,000 livebirths worldwide each year, is associated with significant health care and financial burdens. Clubfoot is defined by forefoot adduction, hindfoot varus, midfoot cavus and hindfoot equinus. Isolated clubfoot, which is the focus of these studies, is distinct from syndromic clubfoot because there are no other associated malformations. Population, family, twin and segregation analysis studies provide evidence that genetic and environmental factors play an etiologic role in isolated clubfoot. The studies described in this thesis were performed to define the role of genetic variation in isolated clubfoot. Interrogation of a deletion region associated with syndromic clubfoot, suggested that CASP8 and CASP10, two apoptotic genes, play a role in isolated clubfoot. To explore the role of apoptotic genes in clubfoot, SNPs spanning genes involved in the apoptotic pathway in the six chromosomal deletion regions, and limb patterning genes, HOXD and HOXA, were interrogated. SNPs in mitochondrial mediated apoptotic genes and several SNPs in HOXA and HOXD genes were modestly associated with clubfoot with the most significant SNP, rs3801776, located in the basal promoter of HOXA9. Several significant associations were found with SNPs in NFAT2 and TNIP2. Significant gene interactions were detected between SNPs in HOX and apoptotic genes. These findings suggest a model for clubfoot in which variation in one gene is not sufficient to cause the malformation but requires variation several genes to perturb protein expression sufficiently to alter muscle and foot development. These results significantly impact our knowledge base by delineating underlying mechanisms causing clubfoot.

Xenobiotic Metabolism Genes and Clubfoot

Idiopathic or isolated clubfoot is a common orthopedic birth defect that affects approximately 135,000 children worldwide. It is characterized by equinus, varus and adductus deformities of the ankle and foot. Correction of clubfoot involves months of serial manipulations, castings and bracing, with surgical correction needed in forty percent of cases. Multifactorial etiology has been suggested in numerous studies with both environmental and genetic factors playing an etiologic role. Maternal smoking during pregnancy is the only common environmental factor that has consistently been shown to increase the risk for clubfoot. Moreover, a positive family history of clubfoot and maternal smoking increases the risk of clubfoot twenty fold. These findings suggest that genetic variation in smoking metabolism genes may increase susceptibility to clubfoot. Based on this reasoning, we interrogated eight candidate genes, chosen based on their involvement in phase 1 and 2 cigarette smoke metabolism. Twenty-two SNPs and two null alleles in eight genes (CYP1A1, CYP1A2, CYP1B1, CYP2A6, EPHX1, NAT2, GSTM1 and GSTT1) were genotyped in a dataset composed of nonHispanic white and Hispanic multiplex and simplex families. Only one SNP in CYP1A1, rs1048943, had significantly altered transmission in the aggregate and multiplex NHW datasets (p=0.003 and p=0.009). Perturbation of CYP1A1 by rs1048943 polymorphism causes an increase in the amount of harmful, adduct forming metabolic intermediates. A significant gene interaction between EPHX1 and NAT2 was also found (p=0.007). This interaction may affect the metabolism of harmful metabolic intermediates. Additionally, marginal interactions were found for other xenobiotic genes and these interactions may play a contributory role in clubfoot. Importantly, for CYP1A2, significant maternal (p=0.03; RR=1.24; 95% CI: 1.04-1.44) and fetal (p=0.01; RR=1.33; 95% CI: 1.13-1.54) genotypic effects were identified suggesting that both maternal and fetal genotypes impact normal limb development. No association was found for maternal smoking status and tobacco metabolism genes. Together, these results suggest that xenobiotic metabolism genes may play a contributory role in the etiology of clubfoot regardless of maternal smoking status and may impact foot development through perturbation of tobacco metabolic pathways.

Altered transmission of HOX and apoptotic SNPs identify a potential common pathway for clubfoot.

Clubfoot is a common birth defect that affects 135,000 newborns each year worldwide. It is characterized by equinus deformity of one or both feet and hypoplastic calf muscles. Despite numerous study approaches, the cause(s) remains poorly understood although a multifactorial etiology is generally accepted. We considered the HOXA and HOXD gene clusters and insulin-like growth factor binding protein 3 (IGFBP3) as candidate genes because of their important roles in limb and muscle morphogenesis. Twenty SNPs from the HOXA and HOXD gene clusters and 12 SNPs in IGFBP3 were genotyped in a sample composed of non-Hispanic white and Hispanic multiplex and simplex families (discovery samples) and a second sample of non-Hispanic white simplex trios (validation sample). Four SNPs (rs6668, rs2428431, rs3801776, and rs3779456) in the HOXA cluster demonstrated altered transmission in the discovery sample, but only rs3801776, located in the HOXA basal promoter region, showed altered transmission in both the discovery and validation samples (P = 0.004 and 0.028). Interestingly, HOXA9 is expressed in muscle during development. An SNP in IGFBP3, rs13223993, also showed altered transmission (P = 0.003) in the discovery sample. Gene-gene interactions were identified between variants in HOXA, HOXD, and IGFBP3 and with previously associated SNPs in mitochondrial-mediated apoptotic genes. The most significant interactions were found between CASP3 SNPS and variants in HOXA, HOXD, and IGFBP3. These results suggest a biologic model for clubfoot in which perturbation of HOX and apoptotic genes together affect muscle and limb development, which may cause the downstream failure of limb rotation into a plantar grade position.

Functional studies of the homeobox gene {\it Goosecoid\/} during mouse embryogenesis

A fundamental question in developmental biology is to understand the mechanisms that govern the development of an adult individual from a single cell. Goosecoid (Gsc) is an evolutionarily conserved homeobox gene that has been cloned in vertebrates and in Drosophila. In mice, Gsc is first expressed during gastrulation stages where it marks anterior structures of the embryo, this pattern of expression is conserved among vertebrates. Later, expression is observed during organogenesis of the head, limbs and the trunk. The conserved pattern of expression of Gsc during gastrulation and gain of function experiments in Xenopus suggested a function for Gsc in the development of anterior structures in vertebrates. Also, its expression pattern in mouse suggested a role in morphogenesis of the head, limbs and trunk. To determine the functional requirement of Gsc in mice a loss of function mutation was generated by homologous recombination in embryonic stem cells and mice mutant for Gsc were generated.^ Gsc-null mice survived to birth but died hours after delivery. Phenotypic analysis revealed craniofacial and rib cage abnormalities that correlated with the second phase of Gsc expression in the head and trunk but no anomalies were found that correlated with its pattern of expression during gastrulation or limb development.^ To determine the mode of action of Gsc during craniofacial development aggregation chimeras were generated between Gsc-null and wild-type embryos. Chimeras were generated by the aggregation of cleavage stage embryos, taking advantage of two different Gsc-null alleles generated during gene targeting. Chimeras demonstrated a cell-autonomous function for Gsc during craniofacial development and a requirement for Gsc function in cartilage and mesenchymal tissues.^ Thus, during embryogenesis in mice, Gsc is not an essential component of gastrulation as had been suggested in previous experiments. Gsc is required for craniofacial development where it acts cell autonomously in cartilage and mesenchymal tissues. Gsc is also required for proper development of the rib cage but it is dispensable for limb development in mice. ^

Functional studies of homeobox gene {\it Cart1\/} during mouse development

Cart1 is a paired-class homeobox-containing gene that is expressed in head mesenchyme, branchial arches, limb buds, and various cartilages during embryogenesis. To understand the role of Cart1 during mammalian development, I generated Cart1-mutant mice by gene targeting in mouse embryonic stem cells. Cart1-homozygous mutants were born alive but all died soon after birth. Most had acrania (absence of the cranial vault) and meroanencephaly (absence of part of the brain). In situ hybridization studies showed that Cart1 is expressed specifically in forebrain mesenchyme but not in midbrain or hindbrain mesenchyme nor in the neural tube. Developmental studies revealed a transient deficiency of forebrain mesenchyme cells due to apoptosis associated with a delay in neural tube closure in that region. Subsequently, the forebrain region became filled with mesenchyme and closed, however, the midbrain neural tube region never initiated closure and remained open. These results suggest that Cart1 is required for the survival of forebrain mesenchyme and that its absence disrupts cranial neural tube morphogenesis by blocking the initiation of closure in the midbrain region, and this ultimately leads to the generation of lethal craniofacial defects. Prenatal treatment of Cart1 homozygous mutants with folic acid suppressed the development of the acrania/meroanencephaly phenotype. Thus, Cart1 mutant mice provide a novel animal model for understanding the cellular, molecular, and genetic etiology of neural tube defects and for the development of prenatal therapeutic protocols using folic acid. ^

Generation of an hTERT-restricted, CMV-augmented oncolytic adenovirus capable of tumor-specific transgene expression: Development, characterization, and optimization

Current shortcomings in cancer therapy require the generation of new, broadly applicable, potent, targeted treatments. Here, an adenovirus is engineered to replicate specifically in cells with active human telomerase promotion using a modified hTERT promoter, fused to a CMV promoter element. The virus was also modified to contain a visible reporter transgene, GFP. The virus, Ad/hTC-GFP-E1 was characterized in vitro and demonstrated tumor specific activity both by dose and over time course experiments in a variety of cell lines. In vivo, Ad/hTC-GFP-E1 was affected at suppressing tumor growth and providing a survival benefit without causing any measurable toxicity. To increase the host range of the vector, the fiber region was modified to contain an RGD-motif. The vector, AdRGD/hTC-GFP-E1, was recharacterized in vitro, revealing heightened levels of infectivity and toxicity however maintaining a therapeutic window between cancer and normal cell toxicity. AdRGD/hTC-GFP-E1 was administered in vivo by limb perfusion and was observed to be tumor specific both in expression and replication. To further enhance the efficacy of viral vectors in lung delivery, asthma medications were investigated for their abilities to enhance transgene delivery and expression. A combination of bronchodilators, mast cell inhibitors, and mucolytic agents was devised which demonstrated fold increases in expression in immunocompetent mouse lungs as single agents and more homogenous, intense levels of expression when done in combination of all agents. To characterize the methods in which some cancers are resistant or may become resistant to oncolytic treatments, several small molecule inhibitors of metabolic pathways were applied in combination with oncolytic infection in vitro. SP600125 and PD 98059, respective JNK and ERK inhibitors, successfully suppressed oncolytic toxicity, however did not affect infectivity or transgene expression of Ad/hTC-GFP-E1. JNK and ERK inhibition did significantly suppress viral replication, however, as analyzed by lysate transfer and titration assays. In contrast, SB 203580, an inhibitor for p38, did not demonstrate any protective effects with infected cells. Flow cytometric analysis indicated a possible correlation with G1 arrest and suppressed viral production, however more compounds must be investigated to clarify this observation. ^