Identification, purification, sequencing and characterization of human lung fibroblast motility -stimulating factor for human soft tissue sarcoma cells

Paracrine motogenic factors, including motility cytokines and extracellular matrix molecules secreted by normal cells, can stimulate metastatic cell invasion. For extracellular matrix molecules, both the intact molecules and the degradative products may exhibit these activities, which in some cases are not shared by the intact molecules. We found that human peritumoral and lung fibroblasts secrete motility-stimulating activity for several recently established human sarcoma cell strains. The motility of lung metastasis-derived human SYN-1 sarcoma cells was preferentially stimulated by human lung and peritumoral fibroblast motility-stimulating factors (FMSFs). FMSFs were nondialyzable, susceptible to trypsin, and sensitive to dithiothreitol. Cycloheximide inhibited accumulation of FMSF activity in conditioned medium; however, addition of cycloheximide to the migration assay did not significantly affect motility-stimulating activity. Purified hepatocyte growth factor/scatter factor (HGF/SF), rabbit anti-hHGF, and RT-PCR analysis of peritumoral and lung fibroblast HGF/SF mRNA expression indicated that FMSF activity was unrelated to HGF/SF. Partial purification of FMSF by gel exclusion chromatography revealed several peaks of activity, suggesting multiple FMSF molecules or complexes.^ We purified the fibroblast motility-stimulating factor from human lung fibroblast-conditioned medium to apparent homogeneity by sequential heparin affinity chromatography and DEAE anion exchange chromatography. Lysylendopeptidase C digestion of FMSF and sequencing of peptides purified by reverse phase HPLC after digestion identified it as an N-terminal fragment of human fibronectin. Purified FMSF stimulated predominantly chemotaxis but chemokinesis as well of SYN-1 sarcoma cells and was chemotactic for a variety of human sarcoma cells, including fibrosarcoma, leiomyosarcoma, liposarcoma, synovial sarcoma and neurofibrosarcoma cells. The motility-stimulating activity present in HLF-CM was completely eliminated by either neutralization or immunodepletion with a rabbit anti-human-fibronectin antibody, thus further confirming that the fibronectin fragment was the FMSF responsible for the motility stimulation of human soft tissue sarcoma cells. Since human soft tissue sarcomas have a distinctive hematogenous metastatic pattern (predominantly lung), FMSF may play a role in this process. ^

DEVELOPMENTAL AND CELLULAR FUNCTIONS OF DELTA-CATENIN

Catenins have diverse and powerful roles in embryogenesis, homeostasis or disease progression, as best exemplified by the well-known beta-catenin. The less studied delta-catenin likewise contains a central Armadillo-domain. In common with other p120 sub-class members, it acts in a variety of intracellular compartments and modulates cadherin stability, small GTPase activities and gene transcription. In mammals, delta-catenin exhibits neural specific expression, with its knock-out in mice correspondingly producing cognitive defects and synaptic dysfunctions. My work instead employed the amphibian, Xenopus laevis, to explore delta-catenin’s physiological functions in a distinct vertebrate system. Initial isolation and characterization indicated delta-catenin’s expression in Xenopus. Unlike the pattern observed for mammals, delta-catenin was detected in most adult Xenopus tissues, although enriched in embryonic structures of neural fate as visualized using RNA in-situ hybridization. To determine delta-catenin’s requirement in amphibian development, I employed anti-sense morpholinos to knock-down gene products, finding that delta-catenin depletion results in developmental defects in gastrulation, neural crest migration and kidney tubulogenesis, phenotypes that were specific based upon rescue experiments. In biochemical and cellular assays, delta-catenin knock-down reduced cadherin levels and cell adhesion, and impaired activation of RhoA and Rac1, small GTPases that regulate actin dynamics and morphogenetic movements. Indeed, exogenous C-cadherin, or dominant-negative RhoA or dominant-active Rac1, significantly rescued delta-catenin depletion. Thus, my results indicate delta-catenin’s essential roles in Xenopus development, with contributing functional links to cadherins and Rho family small G proteins. In examining delta-catenin’s nuclear roles, I identified delta-catenin as an interacting partner and substrate of the caspase-3 protease, which plays critical roles in apoptotic as well as non-apoptotic processes. Delta-catenin’s interaction with and sensitivity to caspase-3 was confirmed using assays involving its cleavage in vitro, as well as within Xenopus apoptotic extracts or mammalian cell lines. The cleavage site, a highly conserved caspase consensus motif (DELD) within Armadillo-repeat 6 of delta-catenin, was identified through peptide sequencing. Cleavage thus generates an amino- (1-816) and carboxyl-terminal (817-1314) fragment each containing about half of the central Armadillo-domain. I found that cleavage of delta-catenin both abolishes its association with cadherins, and impairs its ability to modulate small GTPases. Interestingly, the carboxyl-terminal fragment (817-1314) possesses a conserved putative nuclear localization signal that I found is needed to facilitate delta-catenin’s nuclear targeting. To probe for novel nuclear roles of delta-catenin, I performed yeast two-hybrid screening of a mouse brain cDNA library, resolving and then validating its interaction with an uncharacterized KRAB family zinc finger protein I named ZIFCAT. My results indicate that ZIFCAT is nuclear, and suggest that it may associate with DNA as a transcriptional repressor. I further determined that other p120 sub-class catenins are similarly cleaved by caspase-3, and likewise bind ZIFCAT. These findings potentially reveal a simple yet novel signaling pathway based upon caspase-3 cleavage of p120 sub-family members, facilitating the coordinate modulation of cadherins, small GTPases and nuclear functions. Together, my work suggested delta-catenin’s essential roles in Xenopus development, and has revealed its novel contributions to cell junctions (via cadherins), cytoskeleton (via small G proteins), and nucleus (via ZIFCAT). Future questions include the larger role and gene targets of delta-catenin in nucleus, and identification of upstream signaling events controlling delta-catenin’s activities in development or disease progression.

The C-terminal domain of MinC inhibits assembly of the Z ring in Escherichia coli.

In Escherichia coli, the Min system, consisting of three proteins, MinC, MinD, and MinE, negatively regulates FtsZ assembly at the cell poles, helping to ensure that the Z ring will assemble only at midcell. Of the three Min proteins, MinC is sufficient to inhibit Z-ring assembly. By binding to MinD, which is mostly localized at the membrane near the cell poles, MinC is sequestered away from the cell midpoint, increasing the probability of Z-ring assembly there. Previously, it has been shown that the two halves of MinC have two distinct functions. The N-terminal half is sufficient for inhibition of FtsZ assembly, whereas the C-terminal half of the protein is required for binding to MinD as well as to a component of the division septum. In this study, we discovered that overproduction of the C-terminal half of MinC (MinC(122-231)) could also inhibit cell division and that this inhibition was at the level of Z-ring disassembly and dependent on MinD. We also found that fusing green fluorescent protein to either the N-terminal end of MinC(122-231), the C terminus of full-length MinC, or the C terminus of MinC(122-231) perturbed MinC function, which may explain why cell division inhibition by MinC(122-231) was not detected previously. These results suggest that the C-terminal half of MinC has an additional function in the regulation of Z-ring assembly.

Sequencing of TGF-beta pathway genes in familial cases of intracranial aneurysm.

BACKGROUND AND PURPOSE: Familial aggregation of intracranial aneurysms (IA) strongly suggests a genetic contribution to pathogenesis. However, genetic risk factors have yet to be defined. For families affected by aortic aneurysms, specific gene variants have been identified, many affecting the receptors to transforming growth factor-beta (TGF-beta). In recent work, we found that aortic and intracranial aneurysms may share a common genetic basis in some families. We hypothesized, therefore, that mutations in TGF-beta receptors might also play a role in IA pathogenesis. METHODS: To identify genetic variants in TGF-beta and its receptors, TGFB1, TGFBR1, TGFBR2, ACVR1, TGFBR3, and ENG were directly sequenced in 44 unrelated patients with familial IA. Novel variants were confirmed by restriction digestion analyses, and allele frequencies were analyzed in cases versus individuals without known intracranial disease. Similarly, allele frequencies of a subset of known SNPs in each gene were also analyzed for association with IA. RESULTS: No mutations were found in TGFB1, TGFBR1, TGFBR2, or ACVR1. Novel variants identified in ENG (p.A60E) and TGFBR3 (p.W112R) were not detected in at least 892 reference chromosomes. ENG p.A60E showed significant association with familial IA in case-control studies (P=0.0080). No association with IA could be found for any of the known polymorphisms tested. CONCLUSIONS: Mutations in TGF-beta receptor genes are not a major cause of IA. However, we identified rare variants in ENG and TGFBR3 that may be important for IA pathogenesis in a subset of families.

Studies on the {\it in vivo} function of a phosphatidylinositol transfer protein from {\it Saccharomyces cerevisiae}

Phosphatidylinositol transfer proteins (PI-TP's) catalyze the transfer of phosphatidylinositol and phosphatidylcholine between membranes in vitro. However the in vivo function of these proteins is unknown. In this thesis we have used a combined biochemical and genetic approach to determine the importance of PI-TP in vivo. An oligonucleotide based on the amino terminal sequence of the PI-TP from Saccharomyces cerevisiae, was used to screen a yeast genomic library for the gene encoding PI-TP (PIT1 gene). Yeast strains transformed with the positive clones showed overproduction of transfer activities and transfer protein in the 100,000 x g supernatants. The 5$\sp\prime$ terminus of the PIT1 gene correlates with the predicted codons for residues 3-30 of the determined protein sequence. Tetrad analysis of a heterozygous diploid (PIT1/pit1::LEU2) revealed that the PIT1 gene is essential for cell growth. Non-viable spores could be rescued by transformation of the above diploid prior to sporulation, with a plasmid borne copy of the wild type gene. Sequencing of the entire PIT1 gene has revealed that the PIT1 gene is identical to the SEC14 gene. The sec14 ts mutant which exhibits conditional defects at the Golgi stage of protein secretion, is also temperature sensitive for PI-TP activity in vitro. These findings represent the first instance in which a physiological function has been assigned to any phospholipid transfer protein. ^

Fibrillins: Sequence, phylogeny and carboxy terminal domain proteolytic processing

Fibrillin-1 and -2 are large secreted glycoproteins that are known to be components of extracellular matrix microfibrils located in the vasculature, basement membrane and various connective tissues. These microfibrils are often associated with a superstructure known as the elastic fiber. During the development of elastic tissues, fibrillin microfibrils precede the appearance of elastin and may provide a scaffolding for the deposition and crosslinking of elastin. Using RT/PCR, we cloned and sequenced 3.85Kbp of the FBN2 gene. Five differences were found between our contig sequence and that published by Zhang et al. (1995). Like many extracellular matrix proteins, the fibrillins are modular proteins. We compared analogous domains of the two fibrillins and also members of the latent TGF-$\beta$ binding protein (LTBP) family to determine their phylogenetic relationship. We found that the two families are homologous. LTBP-2 is the most similar to the fibrillin family while FBN-1 is the most similar to the LTBP family. The fibrillin-1 carboxy terminal domain is proteolytically processed. Two eukaryotic protein expression systems, baculoviral and CHO-K1, were developed to examine the proteolytic processing of the carboxy terminal domain of the fibrillin-1 protein. Both expression systems successfully processed the domain and both processed a mutant less efficiently. In the CHO-K1 cells, processing occurred intracellularly. ^

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.^

A MOLECULAR APPROACH FOR DETECTING GENETIC DAMAGE USING THE ECO R1 FAMILY OF HUMAN DNA (SEQUENCING)

A clone of the primary Eco R1 family of human DNA sequences has been used as an indicator sequence for detecting alterations induced by a toxic agent. Specific clones of this family have been examined and compared to the consensus sequence to determine the normal variability of this family. Though variations were observed, data indicated that such clones can be used to study induced DNA modifications. This DNA was exposed to the toxic agent dimethyl sulfate under various conditions and a distinct pattern of aberrations was shown to occur. It is suggested that this approach be used to characterize patterns of damage induced by various agents in the ultimate development of a system capable of monitoring human genotoxic exposure. ^

DEVELOPMENT OF NOVEL METHODS TO MINIMIZE THE IMPACT OF SEQUENCING ERRORS IN THE NEXT-GENERATION SEQUENCING DATA ANALYSIS

Next-generation sequencing (NGS) technology has become a prominent tool in biological and biomedical research. However, NGS data analysis, such as de novo assembly, mapping and variants detection is far from maturity, and the high sequencing error-rate is one of the major problems. . To minimize the impact of sequencing errors, we developed a highly robust and efficient method, MTM, to correct the errors in NGS reads. We demonstrated the effectiveness of MTM on both single-cell data with highly non-uniform coverage and normal data with uniformly high coverage, reflecting that MTM’s performance does not rely on the coverage of the sequencing reads. MTM was also compared with Hammer and Quake, the best methods for correcting non-uniform and uniform data respectively. For non-uniform data, MTM outperformed both Hammer and Quake. For uniform data, MTM showed better performance than Quake and comparable results to Hammer. By making better error correction with MTM, the quality of downstream analysis, such as mapping and SNP detection, was improved. SNP calling is a major application of NGS technologies. However, the existence of sequencing errors complicates this process, especially for the low coverage (

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

The role of the conserved N -terminal domain of STAT3 in STAT3 dephosphorylation and nuclear export

Rapid redistribution of STAT subcellular localization is an essential feature of cytokine signaling. To elucidate the molecular basis of STAT3 function, which plays a critical role in controlling innate immune responses in vivo, we initiated studies to determine the mechanisms controlling STAT3 nuclear trafficking. We found that STAT3 is transported to the nucleus in the absence of cytokine treatment, as judged by indirect immunofluorescence studies in the presence of leptomycin B, an inhibitor of CRM1-dependent nuclear export, suggesting that the non-phosphorylated STAT3 protein contains a functional nuclear import signal. An isoform lacking the STAT3 N-terminal domain (Δ133STAT3) retains the ability to undergo constitutive nuclear localization, indicating that this region is not essential for cytokine-independent nuclear import. Δ133STAT3 is also transported to the nucleus following stimulation with interleukin-6 (IL-6). Interestingly, IL-6-dependent tyrosine phosphorylation of Δ133STAT3 appears to be prolonged and the nuclear export of the protein delayed in cells expressing endogenous STAT3, consistent with defective Δ133STAT3 dephosphorylation. Endogenous STAT3 does not promote the nuclear export of Δ133STAT3, although dimerization between endogenous Stat3 and Δ133STAT3 is detected readily. Thus, the STAT3 N-terminal domain is not required for dimerization with full-length STAT3, yet appears to play a role in proper export of Stat3 from the nucleus following cytokine stimulation. STAT3-deficient cells reconstituted with Δ133STAT3 show enhanced and prolonged Stat1 signaling in response to IL-6, suggesting that induction of the STAT3-dependent negative regulator SOCS3 is impaired. In fact, Δ133STAT3 fails to induce SOCS3 mRNA efficiently. These studies collectively indicate that the STAT3 N-terminal region may be important for IL-6-dependent target gene activation and nuclear dephosphorylation, while dispensable for nuclear import. STAT3 is an oncogene. STAT3 is constitutively activated in primary tumors of many types. Thus far, research in the design of STAT3 protein inhibitors has focused on the SH2 and DNA-binding domains of STAT3. Interference with these domains eliminates all signaling through STAT3. If the N-terminal domain is involved in tetramerization on a subset of target genes, inhibition of this region may lead to a more selective inhibition of some STAT3 functions while leaving others intact. ^