GENETIC ANALYSIS OF THE FUNCTION OF THE DROSOPHILA DOUBLESEX-RELATED FACTOR dmrt93B

DMRT (Doublesex and Mab-3 related transcription factor) proteins generally associated with sexual differentiation in many organisms share a common DNA binding domain and are often expressed in reproductive tissues. Aside from doublesex, which is a central factor in the regulation of sex determination, Drosophila possesses three different dmrt genes that are of unknown function. Because the association with sexual differentiation and reproduction is not universal and some DMRT proteins have been found to play other developmental roles we chose to further characterize one of these Drosophila genes. We carried out genetic analysis of dmrt93B, which was previously found to be expressed sex-specifically in the developing somatic gonad and to affect testis morphogenesis in RNAi knockdowns. In order to disrupt this gene, the GAL4 yeast transcriptional activator followed by a polyadenylation signal was inserted after the dmrt93B start codon and introduced into the genome by homologous recombination. Analysis of the knock-in mutation as well as a small deletion removing all dmrt93B sequence demonstrate that loss of function causes partial lethality at the late pupal stage. Surprisingly, these mutations have no significant effect on gonad formation or male fertility. Analysis of GAL4-driven GFP reporter expression indicates that the dmrt93B promoter activity is highly specific to neurons in the suboesophageal and proventricular ganglion in larva and adult of both sexes suggesting a possible role in digestive tract function. Using the Capillary Feeder (CAFÉ) assay to measure daily food intake we find that reduction in this gene’s function leads to an increase in food consumption. These results suggest dmrt93 plays an important role in the formation or maintenance of neurons that affect feeding and support the idea that dmrt genes may not be restricted to roles in sexual differentiation.

Genetic analysis of factors regulating mesenchymal cell fates

Formation of cartilage and bone involves sequential processes in which undifferentiated mesenchyme aggregates into primordial condensations which subsequently grow and differentiate, resulting in morphogenesis of the adult skeleton. While much has been learned about the structural molecules which comprise cartilage and bone, little is known about the nuclear factors which regulate chondrogenesis and osteogenesis. MHox is a homeobox-containing gene which is expressed in the mesenchyme of facial, limb, and vertebral skeletal precursors during mouse embryogenesis. MHox expression has been shown to require epithelial-derived signals, suggesting that MHox may regulate the epithelial-mesenchymal interactions required for skeletal organogenesis. To determine the functions of MHox, we generated a loss-of-function mutation in the MHox gene. Mice homozygous for a mutant MHox allele exhibit defects of skeletogenesis, involving the loss or malformation of craniofacial, limb and vertebral skeletal structures. The affected skeletal elements are derived from the cranial neural crest, as well as somitic and lateral mesoderm. Analysis of the mutant phenotype during ontogeny demonstrated a defect in the formation or growth of chondrogenic and osteogenic precursors. These findings provide evidence that MHox regulates the formation of preskeletal condensations from undifferentiated mesenchyme. In addition, generation of mice doubly mutant for the MHox and S8 homeobox genes reveal that these two genes interact to control formation of the limb and craniofacial skeleton. Mice carrying mutant alleles for S8 and MHox exhibit an exaggeration of the craniofacial and limb phenotypes observed in the MHox mutant mouse. Thus, MHox and S8 are components of a combinatorial genetic code controlling generation of the skeleton of the skull and limbs. ^

Genetic analysis of tumor progression susceptibility in the mouse skin model

In the field of chemical carcinogenesis the use of animal models has proved to be a useful tool in dissecting the multistage process of tumor formation. In this regard the outbred SENCAR mouse has been the strain of choice in the analysis of skin carcinogenesis given its high sensitivity to the chemically induced acquisition of premalignant lesions, papillomas, and the later progression of these lesions into squamous cell carcinomas (SCC).^ The derivation of an inbred strain from the SENCAR stock called SSIN, that in spite of a high sensitivity to the development of papillomas lack the ability to transform these premalignant lesions into SCC, suggested that tumor promotion and progression were under the genetic control of different sets of genes.^ In the present study the nature of susceptibility to tumor progression was investigated. Analysis of F1 hybrids between the outbred SENCAR and SSIN mice suggested that there is at least one dominant gene responsible for susceptibility to tumor progression.^ Later development of another inbred strain from the outbred SENCAR stock, that had sensitivity to both tumor promotion and progression, allowed the formulation of a more accurate genetic model. Using this newly derived line, SENCAR B/Pt. and SSIN it was determined that there is one dominant tumor progression susceptibility gene. Linkage analysis showed that this gene maps to mouse chromosome 14 and it was possible to narrow the region to a 16 cM interval.^ In order to better characterize the nature of the progression susceptibility differences between these two strains, their proliferative pattern was investigated. It was found that SENCAR B/Pt, have an enlarged proliferative compartment with overexpression of cyclin D1, p16 and p21. Further studies showed an aberrant overexpression of TGF-$\beta$ in the susceptible strain, an increase in apoptosis, p53 protein accumulation and early loss of connexin 26. These results taken together suggest that papillomas in the SENCAR B/Pt. mice have higher proliferation and may have an increase in genomic instability, these two factors would contribute to a higher sensitivity to tumor progression. ^

Biochemical and genetic analysis of two protein kinases JNKK2 and MEKK3

Mitogen-activated protein kinase (MAPK) cascades are conserved eukaryotic signaling modules consisting of a MAPK, a MAPKK and a MAP3K. MAPK cascades are involved in many cellular responses including proliferation, differentiation, apoptosis, stress and immune responses. ^ The first part of this thesis describes the cloning and biochemical analysis of JNKK2, a member of MAPKK gene family. Our results demonstrate that JNKK2 is a specific JNK activator and activates the JNK-dependent signal transduction pathway in vivo by inducing c-Jun and ATF2-mediated gene expression. We also found that JNKK2 is specifically activated by a MAP3K MEKK2 through formation of MEKK2-JNKK2-JNK1 triple complex module. JNKK2 is likely to mediate specific upstream signals to activate JNK cascade. ^ The second part of this thesis describes biochemical and gene disruption analysis of MEKK3, a member of MAP3K gene family. We showed that overexpression of MEKK3 strongly activates both JNK and p38 MAPKs but only weakly activates ERK. MEKK−/− embryos die at about embryonic day (E) 11. MEKK3−/− embryos displayed defects in blood vessel development in the yolk sacs, and in the myocardium and endocardium development at E9.5. The angiogenesis in the head, intersomitic region and placenta was also abnormal. These results demonstrate that MEKK3, a member of MAP3K MEKK/STE11 subgene family, is essential for early embryonic cardiovascular development. Furthermore, it was found that disruption of MEKK3 did not alter the expression of vascular endothelial growth factor-1 (VEGF-1), angiopoietin-1, -2 and their respective receptors Flt-1, Flk-1, Tie-1, Tie-2. Finally, MEKK3 was shown to activate myocyte-specific enhancer factor 2C (MEF2C), a crucial transcription factor for early embryonic cardiovascular development through the p38 MAPK cascade, suggesting that MEF2C is one of the key targets of the MEEKK3 signaling pathway during early embryonic cardiovascular development. ^

GENETIC ANALYSIS OF THE HIPPO PATHWAY IN MOUSE LIVER

Cancer therapy and tumor treatment remain unsolved puzzles. Genetic screening for tumor suppressor genes in Drosophila revealed the Hippo-signaling pathway as a kinase cascade consisting of five core components. Disrupting the pathway by deleting the main component genes breaks the balance of cell proliferation and apoptosis and results in epithelial tissue tumorigenesis. The pathway is therefore believed to be a tumor suppressor pathway. However, a corresponding role in mammals is yet to be determined. Our lab began to investigate the tumor suppression function of the potent mammalian Hippo pathway by putting floxed alleles into the mouse genome flanking the functional-domain-expressing exons in each component (Mst1, Mst2, Sav1, Lats1 and Lats2). These mice were then crossed with different cre-mouse lines to generate conditional knockout mice. Results indicate a ubiquitous tumor suppression function of these components, predominantly in the liver. A further liver specific analysis of the deletion mutation of these components, as well as the Yap/Taz double deletion mutation, reveals essential roles of the Hippo pathway in regulating hepatic quiescence and embryonic liver development. One of the key cellular mechanisms for the Hippo pathway’s involvement in these liver biological events is likely its cell cycle regulation function. Our work will help to develop potential therapeutic approaches for liver cancer.

Molecular and genetic analysis of Myxococcus xanthus early *development

To identify more mutations that can affect the early development of Myxococcus xanthus, the synthetic transposon TnT41 was designed and constructed. By virtue of its special features, it can greatly facilitate the processes of mutation screening/selection, mapping, cloning and DNA sequencing. In addition, it allows for the systematic discovery of genes in regulatory hierarchies using their target promoters. In this study, the minimal regulatory region of the early developmentally regulated gene 4521 was used as a reporter in the TnT41 mutagenesis. Both positive (P) mutations and negative (N) mutations were isolated based on their effects on 4521 expression.^ Four of these mutations, i.e. N1, N2, P52 and P54 were analyzed in detail. Mutations N1 and N2 are insertion mutations in a gene designated sasB. The sasB gene is also identified in this study by genetic and molecular analysis of five UV-generated 4521 suppressor mutations. The sasB gene encodes a protein without meaningful homology in the databases. The sasB gene negatively regulates 4521 expression possibly through the SasS-SasR two component system. A wild-type sasB gene is required for normal M. xanthus fruiting body formation and sporulation.^ Cloning and sequencing analysis of the P52 mutation led to the identification of an operon that encodes the M. xanthus high-affinity branched-chain amino acid transporter system. This liv operon consists of five genes designated livK, livH, livM, livC, and livF, respectively. The Liv proteins are highly similar to their counterparts from other bacteria in both amino acid sequences, functional motifs and predicted secondary structures. This system is required for development since liv null mutations cause abnormality in fruiting body formation and a 100-fold decrease in sporulation efficiency.^ Mutation P54 is a TnT41 insertion in the sscM gene of the ssc chemotaxis system, which has been independently identified by Dr. Shi's lab. The sscM gene encodes a MCP (methyl-accepting chemotaxis protein) homologue. The SscM protein is predicted to contain two transmembrane domains, a signaling domain and at least one putative methylation site. Null mutations of this gene abolish the aggregation of starving cells at a very early stage, though the sporulation levels of the mutant can reach 10% that of wild-type cells. ^

Genetic analysis of Pitx2 in left -right asymmetry and cardiovascular development

Pitx2, a paired-related homeobox gene that is mutated in human Rieger Syndrome, plays a key role in transferring the early asymmetric signals to individual organs. Pitx2 encodes three isoforms, Pitx2a, Pitx2b and Pitx2c. I found that Pitx2c was the Pitx2 isoform for regulating left-right asymmetry in heart, lung and the predominant isoform in guts. Previous studies suggested that the generation of left-right asymmetry within individual organs is an all or none, random event. Phenotypic analysis of various Pitx2 allelic combinations, that encode graded levels of Pitx2c, reveals an organ-intrinsic mechanism for regulating left-right asymmetric morphogenesis based on differential response to Pitx2c levels. The heart needs low Pitx2c levels, while the lungs and duodenum require higher doses of Pitx2c. In addition, the duodenal rotation is under strict control of Pitx2c activity. Left-right asymmetry development for aortic arch arteries involves complex vascular remodeling. Left-sided expression of Pitx2c in these developing vessels implied its potential function in this process. In order to determine if Pitx2c also can regulate the left-right asymmetry of the aortic arch arteries, a Pitx2c-specific loss of function mutation is generated. Although in wild type mice, the direction of the aortic arch is always oriented toward the left side, the directions of the aortic arches in the mutants were randomized, showing that Pitx2c also determined the left-right asymmetry of these vessels. I have further showed that the cardiac neural crest wasn't involved in this vascular remodeling process. In addition, all mutant embryos had Double Outlet Right Ventricle (DORV), a common congenital heart disease. This study provided insight into the mechanism of Pitx2c-mediated late stages of left-right asymmetry development and identified the roles of Pitx2c in regulation of aortic arch remodeling and heart development. ^

{\it Hox\/} genes and specification of regional identity in the axial skeleton: Analysis of the individual and combined mutant phenotypes of the paralogous group 4 murine {\it Hox\/} genes; {\it hoxa\/}-4, {\it hoxb\/}-4 and {\it hoxd\/}-4

The Hox gene products are transcription factors involved in specifying regional identity along the anteroposterior body axis. In Drosophila, where these genes are known as HOM-C (Homeotic-complex) genes and where they have been most extensively studied, they are expressed in restricted domains along the anteroposterior axis with different anterior limits. Genetic analysis of a large number of gain- and loss-of-function alleles of these genes has revealed that these genes are important in specifying segmental identity at their anterior limits of expression. Furthermore, there is a functional dominance of posterior genes over anterior genes, such that posterior genes can dominantly specify their developmental programs in spite of the expression of more anterior genes in the same segment. In the mouse, there are four clusters of HOM-C genes, called Hox genes. Thus, there may be up to four genes, called paralogs, that are more highly homologous to each other and to their Drosophila homolog than they are to the other mouse Hox genes. The single mutants for two paralogous genes, hoxa-4 and hoxd-4, presented in this dissertation, are similar to several other mouse Hox mutants in that they show partial, incompletely penetrant homeotic transformations of vertebrae at their anterior limit of expression. These mutants were then bred with hoxb-4 mutants (Ramirez-Solis, et al. 1993) to generate the three possible double mutant combinations as well as the triple mutant. The skeletal phenotypes of these group 4 Hox compound mutants displayed clear alterations in regional identity, such that a nearly complete transformation towards the morphology of the first cervical vertebra occurs. These results suggest a certain degree of functional redundancy among paralogous genes in specifying regional identity. Furthermore, there was a remarkable dose-dependent increase in the number of vertebrae transformed to a first cervical vertebra identity, including the second through the fifth cervical vertebrae in the triple mutant. Thus, these genes are required in a larger anteroposterior domain than is revealed by the single mutant phenotypes alone, such that multiple mutations in these genes result in transformations of vertebrae that are not at their anterior limit of expression. ^

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

Genetic epidemiology of gallbladder disease in Mexican Americans and cholesterol 7$\alpha$-hydroxylase gene sequence variation

Among Mexican Americans, the second largest minority group in the United States, the prevalence of gallbladder disease is markedly elevated. Previous data from both genetic admixture and family studies indicate that there is a genetic component to the occurrence of gallbladder disease in Mexican Americans. However, prior to this thesis no formal genetic analysis of gallbladder disease had been carried out nor had any contributing genes been identified.^ The results of complex segregation analysis in a sample of 232 Mexican American pedigrees documented the existence of a major gene having two alleles with age- and gender-specific effects influencing the occurrence of gallbladder disease. The estimated frequency of the allele increasing susceptibility was 0.39. The lifetime probabilities that an individual will be affected by gallbladder disease were 1.0, 0.54, and 0.00 for females of genotypes "AA", "Aa", and "aa", respectively, and 0.68, 0.30, and 0.00 for males, respectively. This analysis provided the first conclusive evidence for the existence of a common single gene having a large effect on the occurrence of gallbladder disease.^ Human cholesterol 7$\alpha$-hydroxylase is the rate-limiting enzyme in bile acid synthesis. The results of an association study in both a random sample and a matched case/control sample showed that there is a significant association between cholesterol 7$\alpha$-hydroxylase gene variation and the occurrence of gallbladder disease in Mexican Americans males but not in females. These data have implicated a specific gene, 7$\alpha$-hydroxylase, in the etiology of gallbladder disease in this population.^ Finally, I asked whether the inferred major gene from complex segregation analysis is genetically linked to the cholesterol 7$\alpha$-hydroxylase gene. Three pedigrees predicted to be informative for linkage analysis by virtue of supporting the major gene hypothesis and having parents with informative genotypes and multiple offspring were selected for this linkage analysis. In each of these pedigrees, the recombination fractions maximized at 0 with a positive, albeit low, LOD score. The results of this linkage analysis provide preliminary and suggestive evidence that the cholesterol 7$\alpha$-hydroxylase gene and the inferred gallbladder disease susceptibility gene are genetically linked. ^

A genetic characterization of {\it runt\/} activity during {\it Drosophila\/} embryogenesis

The Drosophila melanogaster gene runt encodes a novel transcriptional regulator that was originally identified on the basis of its key role in embryonic pattern formation. For my thesis I undertook a genetic analysis of runt activity to identify loci that interact with this unique transcriptional regulator. Specifically, I screened the genome with deficiencies for loci that interact with runt in a dose-dependent fashion during early embryogenesis. From this screen I discovered a vital dose-dependent interaction between runt and the achaete-scute complex (AS-C). The characterization of this interaction led to the exciting discovery of important roles for runt in sex determination and neurogenesis (Duffy and Gergen 1991, Duffy et al. 1991). I demonstrated that in sex determination runt is necessary for the normal transcriptional activation of the master sex-determining gene Sx1 and has all the properties of an X:A numerator element. I also showed that runt is required during the early stages of neurogenesis for the normal development of a subset of CNS ganglion mother cells and neurons. In addition, the screen, which focused on the identification and characterization of maternal loci that influence the activity of runt during segmentation, identified several new maternal loci, one of which affects the activity of the maternal posterior group genes on embryonic pattern formation. ^

Accessing novel developmental mechanisms in the mouse by gene trapping

Genetic analysis is a powerful method for analyzing the function of specific genes in development. I sought to identify novel genes in the mouse using a genetic analysis relying on the expression pattern and phenotype of mutated genes. To this end, I have conducted a gene trap screen using the vector $\rm SA\beta geo,$ a promoterless DNA construct that encodes a fusion protein with lacZ and neomycin resistance activities. Productive integration and expression of the $\beta$geo protein in embryonic stem (ES) cells requires integration into an active transcription unit. The endogenous regulatory elements direct reporter gene expression which reflects the expression of the endogenous gene. Of eight mouse lines generated from gene trap ES cell clones, four showed differential regulation of $\beta$geo activity during embryogenesis. These four were analyzed in more detail.^ Three of the lines RNA 1, RNA2 and RNA 3 had similar expression patterns, within subsets of cells in sites of embryonic hematopoiesis. Cloning of the trapped genes revealed that all three integrations had occurred within 45S rRNA precursor transcription units. These results imply that there exists in these cells some mechanism responsible for the efficient production of the $\beta$geo protein from an RNA polymerase I transcript that is not present in most of the cells in the embryo.^ The fourth line, GT-2, showed widespread, dynamic expression. Many of the sites of expression were important classic embryonic induction systems. Cloning of the sequences fused to the $5\sp\prime$ end of the $\beta$geo sequence revealed that the trapped gene contained significant sequence homology with a previously identified human sequence HumORF5. An open reading frame of this sequence is homologous to a group of eukaryotic proteins that are members of the RNA helicase superfamily I.^ Analysis of the gene trap lines suggests that potentially novel developmental mechanisms have been uncovered. In the case of RNA 1, 2 and 3, the differential production of ribosomal RNAs may be required for differentiation or function of the $\beta$geo positive hematopoietic cells. In the GT-2 line, a previously unsuspected temporal and spatial regulation of a putative RNA helicase implies a role for this activity during specific aspects of mouse development. ^

Genes to blood pressure: The effects of variation in genes of the renin-angiotensin system on the hormonal and renal control of blood pressure

Objective. Essential hypertension affects 25% of the US adult population and is a leading contributor to morbidity and mortality. Because BP is a multifactorial phenotype that resists simple genetic analysis, intermediate phenotypes within the complex network of BP regulatory systems may be more accessible to genetic dissection. The Renin-Angiotensin System (RAS) is known to influence intermediate and long-term blood pressure regulation through alterations in vascular tone and renal sodium and fluid resorption. This dissertation examines associations between renin (REN), angiotensinogen (AGT), angiotensin-converting enzyme (ACE) and angiotensin II type 1 receptor (AT1) gene variation and interindividual differences in plasma hormone levels, renal hemodynamics, and BP homeostasis.^ Methods. A total of 150 unrelated men and 150 unrelated women, between 20.0 and 49.9 years of age and free of acute or chronic illness except for a history of hypertension (11 men and 7 women, all off medications), were studied after one week on a controlled sodium diet. RAS plasma hormone levels, renal hemodynamics and BP were determined prior to and during angiotensin II (Ang II) infusion. Individuals were genotyped by PCR for a variable number tandem repeat (VNTR) polymorphism in REN, and for the following restriction fragment length polymorphisms (RFLP): AGT M235T, ACE I/D, and AT1 A1166C. Associations between clinical measurements and allelic variation were examined using multiple linear regression statistical models.^ Results. Women homozygous for the AT1 1166C allele demonstrated higher intracellular levels of sodium (p = 0.044). Men homozygous for the AGT T235 allele demonstrated a blunted decrement in renal plasma flow in response to Ang II infusion (p = 0.0002). There were no significant associations between RAS gene variation and interindividual variation in RAS plasma hormone levels or BP.^ Conclusions. Rather than identifying new BP controlling genes or alleles, the study paradigm employed in this thesis (i.e., measured genes, controlled environments and interventions) may provide mechanistic insight into how candidate genes affect BP homeostasis. ^

Expression cloning of a candidate prostate cancer oncogene

Genetic analysis, both karyotyping and comparative genomic hybridization, of prostate cancer cell lines and specimens have revealed multiple areas of concordant increases in DNA content. An increase of DNA in specific regions of the genome in cancer is often associated with the amplification of oncogenes. Based on these observations we have hypothesized that oncogenes are involved in the initiation or progression of prostate cancer. An expression cloning approach was utilized to identify candidate oncogenes in prostate cancer. ^ A full-length, unidirectional cDNA expression library was constructed from DU145 prostate cancer cells. The cDNA library was screened using CP12, a rat prostate epithelial cell line. In soft agarose assays, CP12 (parental or vector transfected) do not form colonies. However, upon the introduction of a number of known oncogenes CP12 becomes anchorage independent in soft agarose. Based on this in-vitro phenotypic shift, a DU145 cDNA library was stably transfected into CP12, and selected for anchorage independence. Two hundred fifty nine anchorage independent clones were isolated. Some colonies contained more than one insert, bringing the candidate oncogene pool to approximately 400. Seven inserts were sequenced at random. Using the sequences obtained, GenBank was screened, and matches were found with p53, PARG1, a mitochondrial ATPase, RNF6, and three unknown genes that mapped to Unigene clusters. As the pool of cDNA inserts appeared promising, overexpressed genes were further selected. From 259 clones, 17 clones were overexpressed more than 6-fold in DU145 compared to Normal Prostate. From the 17 clones, 12 cDNA inserts were strongly expressed in DU145 and were isolated for sequencing. ^ Two of the sequences, 1G6 and 3E9, were identical. Expression of 1G6/2G9/3E9 was tested by RT-PCR. 1G6/2G9/3E9 was not expressed in normal prostate, but was expressed in all prostate cancer cell lines tested as well as six prostate cancer samples. When retransfected into CP12, 1G6/2G9/3E9 induced the formation of foci and anchorage independent colonies. Thus, functional and expression data suggest that 1G6/2G9/3E9 may be a prostate cancer oncogene. ^