Pathway analysis of genome-wide association study data in primary biliary cirrhosisa

Genome-wide association studies (GWAS) have successfully identified several genetic loci associated with inherited predisposition to primary biliary cirrhosis (PBC), the most common autoimmune disease of the liver. Pathway-based tests constitute a novel paradigm for GWAS analysis. By evaluating genetic variation across a biological pathway (gene set), these tests have the potential to determine the collective impact of variants with subtle effects that are individually too weak to be detected in traditional single variant GWAS analysis. To identify biological pathways associated with the risk of development of PBC, GWAS of PBC from Italy (449 cases and 940 controls) and Canada (530 cases and 398 controls) were independently analyzed. The linear combination test (LCT), a recently developed pathway-level statistical method was used for this analysis. For additional validation, pathways that were replicated at the P <0.05 level of significance in both GWAS on LCT analysis were also tested for association with PBC in each dataset using two complementary GWAS pathway approaches. The complementary approaches included a modification of the gene set enrichment analysis algorithm (i-GSEA4GWAS) and Fisher's exact test for pathway enrichment ratios. Twenty-five pathways were associated with PBC risk on LCT analysis in the Italian dataset at P<0.05, of which eight had an FDR<0.25. The top pathway in the Italian dataset was the TNF/stress related signaling pathway (p=7.38×10 -4, FDR=0.18). Twenty-six pathways were associated with PBC at the P<0.05 level using the LCT in the Canadian dataset with the regulation and function of ChREBP in liver pathway (p=5.68×10-4, FDR=0.285) emerging as the most significant pathway. Two pathways, phosphatidylinositol signaling system (Italian: p=0.016, FDR=0.436; Canadian: p=0.034, FDR=0.693) and hedgehog signaling (Italian: p=0.044, FDR=0.636; Canadian: p=0.041, FDR=0.693), were replicated at LCT P<0.05 in both datasets. Statistically significant association of both pathways with PBC genetic susceptibility was confirmed in the Italian dataset on i-GSEA4GWAS. Results for the phosphatidylinositol signaling system were also significant in both datasets on applying Fisher's exact test for pathway enrichment ratios. This study identified a combination of known and novel pathway-level associations with PBC risk. If functionally validated, the findings may yield fresh insights into the etiology of this complex autoimmune disease with possible preventive and therapeutic application.^

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.

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 (

New methods for quantification and analysis of quantitative real-time polymerase chain reaction data

Quantitative real-time polymerase chain reaction (qPCR) is a sensitive gene quantitation method that has been widely used in the biological and biomedical fields. The currently used methods for PCR data analysis, including the threshold cycle (CT) method, linear and non-linear model fitting methods, all require subtracting background fluorescence. However, the removal of background fluorescence is usually inaccurate, and therefore can distort results. Here, we propose a new method, the taking-difference linear regression method, to overcome this limitation. Briefly, for each two consecutive PCR cycles, we subtracted the fluorescence in the former cycle from that in the later cycle, transforming the n cycle raw data into n-1 cycle data. Then linear regression was applied to the natural logarithm of the transformed data. Finally, amplification efficiencies and the initial DNA molecular numbers were calculated for each PCR run. To evaluate this new method, we compared it in terms of accuracy and precision with the original linear regression method with three background corrections, being the mean of cycles 1-3, the mean of cycles 3-7, and the minimum. Three criteria, including threshold identification, max R2, and max slope, were employed to search for target data points. Considering that PCR data are time series data, we also applied linear mixed models. Collectively, when the threshold identification criterion was applied and when the linear mixed model was adopted, the taking-difference linear regression method was superior as it gave an accurate estimation of initial DNA amount and a reasonable estimation of PCR amplification efficiencies. When the criteria of max R2 and max slope were used, the original linear regression method gave an accurate estimation of initial DNA amount. Overall, the taking-difference linear regression method avoids the error in subtracting an unknown background and thus it is theoretically more accurate and reliable. This method is easy to perform and the taking-difference strategy can be extended to all current methods for qPCR data analysis.^

Functional analysis of the sea urchin {\it orthodenticle\/} -related gene, {\it SPOTX\/}, in aboral ectoderm-specific gene activation and aboral ectoderm cell fate specification

An important question in developmental biology is how embryonic cell types are derived from a fertilized egg. To address this question, this thesis investigates the mechanisms by which the aboral ectoderm-specific Spec2a gene is spatially and temporally regulated during sea urchin embryogenesis. The Spec2a gene of the sea urchin Strongylocentratus purpuratus has served as a valuable maker to understand the basis of lineage-specific gene activation and the role of transcription factors in cell fate specification. The hypothesis is that transcription factors responsible for cell type-specific gene activation are key components in the initial cell specification step. The Spec2a gene, which encodes a small cytosolic calcium-binding protein, is expressed exclusively in aboral ectoderm cell lineages. The 1516-bp control region of the Spec2a gene contains a 188-bp enhancer element required for temporal activation and aboral ectoderm/mesenchyme cell expression, while an unidentified element upstream of the enhancer represses expression in mesenchyme cells. Using an enhancer activation assay, combined with site-directed mutagenesis, I showed that three TAATCC/T sites within the enhancer are responsible for enhancer activity. Mutagenizing these sites and a fourth one just upstream abolished all activity from the Spec2a control region. A 77-bp DNA fragment from the Spec2a enhancer containing two of the TAATCC/T sites is sufficient for aboral ectoderm/mesenchyme cell expression. A cDNA encoding SpOtx, an orthodenticle-related protein, was cloned from S. purpuratus and shown to bind with high affinity to the TAATCC/T sequences within the Spec2a control region. SpOtx transcripts were found initially in all cells of the cleaving embryo, but they gradually became restricted to oral ectoderm and endoderm cells, suggesting that SpOtx might play a role in the initial temporal activation of the Spec2a gene and most likely has additional functions in the developing embryo. To reveal the broader biological functions of SpOtx, I injected SpOtx mRNA into living sea urchin eggs to determine what effects overexpressing the SpOtx protein might have on embryo development. SpOtx mRNA-injected embryos displayed dramatic alterations in development. Instead of developing into pluteus larvae with 15 different cell types, uniform epithelia balls were formed. These balls consisted of a thin layer of squamous cells with short cilia highly reminiscent of aboral ectoderm. Immunohistochemical staining and RT-PCR demonstrated that the SpOtx-injected embryoids expressed aboral ectoderm markers uniformly, but showed very weak or no expression of markers for non-aboral ectoderm cell types. These data strongly suggested that overexpression of SpOtx redirected the normal fate of non-aboral ectoderm cells to that of aboral ectoderm. These results show that SpOtx is involved in aboral ectoderm differentiation by activating aboral ectoderm-specific genes and that modulating its expression can lead to changes in cell fate. ^

Biological effects of PEA3 expression in {\it HER-2\/}/{\it neu\/}-overexpressing human cancer cells and the molecular mechanisms of PEA3-mediated transcriptional repression on {\it HER-2\/}/{\it neu\/}

HER-2/neu is a receptor tyrosine kinase highly homologous with epidermal growth factor receptor. Overexpression and/or amplification of HER-2/neu has been implicated in the genesis of a number of human cancers, especially breast and ovarian cancers. Transcriptional upregulation has been shown to contribute significantly to the overexpression of this gene. Studies on the transcriptional regulation of HER-2/neu gene are important for understanding the mechanism of cell transformation and developing the therapeutic strategies to block HER-2/neu-mediated cancers. PEA3 is a DNA binding transcriptional factor and its consensus sequence exists on the HER-2/neu promoter. To examine the role of PEA3 in HER-2/neu expression and cell transformation, we transfected PEA3 into the human breast and ovarian cancer cells that overexpress HER-2/neu and showed that PEA3 dramatically represses HER-2/neu transcription. PEA3 suppresses the oncogenic neu-mediated transformation in mouse fibroblast NIH 3T3 cells. Expression of PEA3 selectively blocks the growth of human cancer cells that overexpress HER-2/neu and inhibits their colony formation. It does not occur in the cancer cells expressing basal level of HER-2/neu. Further studies in the orthotopic ovarian cancer model demonstrated that expression of PEA3 preferentially inhibits growth and tumor development of human cancer cells that overexpress HER-2/neu, the tumor-bearing mice survived significantly longer if treated by injection of the PEA3-liposome complex intraperitoneally. Immunoblotting and immunohistochemical analysis of the tumor tissues indicated that PEA3 mediates the tumor suppression activity through targeting HER-2/neu-p185. Thus, PEA3 is a negative regulator of HER-2/neu gene expression and functions as a tumor suppressor gene in the HER-2/neu-overexpressing human cancer cells.^ The molecular mechanisms of PEA3 mediated transcriptional repression were investigated. PEA3 binds specifically at the PEA3 site on HER-2/neu promoter and this promoter-binding is required for the PEA3 mediated transcriptional repression. Mutation of the PEA3 binding site on HER-2/neu promoter causes decreased transcriptional activity, indicating that the PEA3 binding site is an enhancer-like element in the HER-2/neu-overexpressing cells. We therefore hypothesized that in the HER-2/neu-overexpressing cells, PEA3 competes with a transactivator for binding to the PEA3 site, preventing the putative factor from activating the transcription of HER-2/neu. This hypothesis was supported by the data which demonstrate that PEA3 competes with another nuclear protein for binding to the HER-2/neu promoter in vitro, and expression of a truncated protein which encodes the DNA binding domain of PEA3 is sufficient to repress HER-2/neu transcription in the HER-2/neu-overexpressing human cancer cells. ^

The cDNA cloning of murine HIP /L29 and the functional analysis of HIP/L29 domains

Human heparin/heparan sulfate interacting protein/L29 (HIP/L29) is a heparin/heparan sulfate (Hp/HS) binding protein found in many adult human tissues. Potential functions of this protein are promotion of embryo adhesion, modulation of blood coagulation, and control of cell growth. While these activities are diverse, the ability of human HIP/L29 to interact with Hp/HS at the cell surface may be a unifying mechanism of action since Hp/HS influences all of these processes. A murine ortholog has been identified that has 78.8% homology over the entire sequence and identity over the N-terminal 64 amino acids when compared to human HIP/L29. Northern, Western, and immunohistochemical analysis shows that murine HIP/L29 mRNA and protein are expressed in a tissue specific manner. Murine HIP/L29 is enriched in the membrane fraction of NmuMG cells where it is eluted with high salt, suggesting that it is a peripheral membrane protein. The ability of murine HIP/L29 to bind Hp is verified by studies using native and recombinant forms of murine HIP/L29. A synthetic peptide (HIP peptide-2) derived from the identical N-terminal region of HIP/L29 proteins was tested for the ability to bind Hp and support cell adhesion. This peptide was chosen because it conforms to a proposed consensus sequence for Hp/HS binding peptides. HIP peptide-2 binds Hp in a dose-dependent, saturable, and selective manner and supports Hp-dependent cell adhesion. However, a scrambled form of this peptide displayed similar activities indicating a lack of peptide sequence specificity required for activity. Lastly, an unbiased approach was used to identify sequences within human and mouse HIP/L29 proteins necessary for Hp/HS binding. A panel of recombinant proteins was made that collectively are deficient in every human HIP/L29 domain. The activities of these deletion mutants and recombinant murine HIP/L29 were compared to the activity of recombinant human HIP/L29 in a number of assays designed to look at differences in the ability to bind Hp/HS. These studies suggest that each domain within human HIP/L29 is important for binding to Hp/HS and divergences in the C-terminus of human and mouse HIP/L29 account for a decrease in murine HIP/L29 affinity for Hp/HS. It is apparent that multiple domains within human and mouse HIP/L29 contribute to the function of Hp/HS binding. The interaction of multiple HIP/L29 domains with Hp/HS will influence the biological activity of HIP/L29 proteins. ^

Interactions and structural analysis of the zinc -binding motif in decorin, a small leucine rich proteoglycan in extracellular matrix

Decorin, a dermatan/chondroitin sulfate proteoglycan, is ubiquitously distributed in the extracellular matrix (ECM) of mammals. Decorin belongs to the small leucine rich proteoglycan (SLRP) family, a proteoglycan family characterized by a core protein dominated by Leucine Rich Repeat motifs. The decorin core protein appears to mediate the binding of decorin to ECM molecules, such as collagens and fibronectin. It is believed that the interactions of decorin with these ECM molecules contribute to the regulation of ECM assembly, cell adhesions, and cell proliferation. These basic biological processes play critical roles during embryonic development and wound healing and are altered in pathological conditions such as fibrosis and tumorgenesis. ^ In this dissertation, we discover that decorin core protein can bind to Zn2+ ions with high affinity. Zinc is an essential trace element in mammals. Zn2+ ions play a catalytic role in the activation of many enzymes and a structural role in the stabilization of protein conformation. By examining purified recombinant decorin and its core protein fragments for Zn2+ binding activity using Zn2+-chelating column chromatography and Zn2+-equilibrium dialysis approaches, we have located the Zn2+ binding domain to the N-terminal sequence of the decorin core protein. The decorin N-terminal domain appears to contain two Zn2+ binding sites with similar high binding affinity. The sequence of the decorin N-terminal domain does not resemble any other reported zinc-binding motifs and, therefore, represents a novel Zn 2+ binding motif. By investigating the influence of Zn2+ ions on decorin binding interactions, we found a novel Zn2+ dependent interaction with fibrinogen, the major plasma protein in blood clots. Furthermore, a recombinant peptide (MD4) consisting of a 41 amino acid sequence of mouse decorin N-terminal domain can prolong thrombin induced fibrinogen/fibrin clot formation. This suggests that in the presence of Zn2+ the decorin N-terminal domain has an anticoagulation activity. The changed Zn2+-binding activities of the truncated MD4 peptides and site-directed mutagenesis generated mutant peptides revealed that the functional MD4 peptide might contain both a structural zinc-binding site in the cysteine cluster region and a catalytic zinc site that could be created by the flanking sequences of the cysteine cluster region. A model of a loop-like structure for MD4 peptide is proposed. ^