A FUNCTIONAL VARIANT OF HLA-DR1 DELINEATED BY T-LYMPHOCYTE CLONE REACTIVITY, PROTEIN CONFORMATION, AND GENOMIC RESTRICTION SITES

This research characterized a serologically indistinguishable form of HLA-DR1 that: (1) cannot stimulate some DR1-restricted or specific T-lymphocyte clones; (2) displays an unusual electrophoretic pattern on two dimensional gels; and (3) is marked by a polymorphic restriction site of the alpha gene. Inefficient stimulation of some DR1-restricted clones was a property of DR1$\sp{+}$ cells that shared HLA-B14 on the same haplotype and/or were carriers of 21-hydroxylase (21-OH) deficiency. Nonclassical 21-OH deficiency frequently demonstrates genetic linkage with HLA-B14;DR1 haplotypes and associates with duplications of C4B and one 21-OH gene. Cells having both stimulatory (DR1$\sb{\rm n}$) and nonstimulatory (DR1$\sb{\rm x}$) parental haplotypes did not mediate proliferation of these clones. However, heterozygous DR1$\sb{\rm x}$, 2 and DR1$\sb{\rm x}$, 7 cells were efficient stimulators of DR2 and DR7 specific clones, respectively, suggesting that a trans acting factor may modify DR1 alleles or products to yield a dominant DR1$\sb{\rm x}$ phenotype. Incompetent stimulator populations did not secrete an intercellular soluble or contact dependent suppressor factor nor did they express interleukin-2 receptors competing for T-cell growth factors. Two dimensional gel analysis of anti-DR immunoprecipitates revealed, in addition to normal DR$\alpha$ and DR$\beta$ chains, a 50kD species from DR1$\sb{\rm x}$ but not from the majority of DR1$\sb{\rm n}$ or non-DR1 cells. The 50kD structure was stable under reducing conditions in SDS and urea, had antigenic homology with DR, and dissociated after boiling into 34kD and 28kD peptide chains apparently identical with DR$\alpha$ and DR$\beta$ as shown by limited digest peptide maps. N-linked glycosylation and sialation of DRgp50 appeared to be unchanged from normal DR$\alpha$ and DR$\beta$. Bg1II digestion and $DR\alpha$ probing of DR1$\sb{\rm x}$ genomic DNA revealed a 4.5kb fragment while DR1$\sb{\rm n}$ DNA yielded 3.8 and 0.76kb fragments; all restriction sites mapped to the 3$\sp\prime$ untranslated region of $DR\alpha$. Collectively, these data suggest that DRgp50 represents a novel combinatorial association between constitutive chains of DR that may interfere with or compete for normal T cell receptor recognition of DR1 as both an alloantigen and restricting element. Furthermore, extensive chromosomal abnormalities previously mapped to the class III region of B14;DR1 haplotypes may extend into the adjacent class II region with consequent intrusion on immune function. ^

NETWORK TOPOLOGY IN HUMAN PROTEIN INTERACTION DATA PREDICTS FUNCTIONAL ASSOCIATION

High-throughput assays, such as yeast two-hybrid system, have generated a huge amount of protein-protein interaction (PPI) data in the past decade. This tremendously increases the need for developing reliable methods to systematically and automatically suggest protein functions and relationships between them. With the available PPI data, it is now possible to study the functions and relationships in the context of a large-scale network. To data, several network-based schemes have been provided to effectively annotate protein functions on a large scale. However, due to those inherent noises in high-throughput data generation, new methods and algorithms should be developed to increase the reliability of functional annotations. Previous work in a yeast PPI network (Samanta and Liang, 2003) has shown that the local connection topology, particularly for two proteins sharing an unusually large number of neighbors, can predict functional associations between proteins, and hence suggest their functions. One advantage of the work is that their algorithm is not sensitive to noises (false positives) in high-throughput PPI data. In this study, we improved their prediction scheme by developing a new algorithm and new methods which we applied on a human PPI network to make a genome-wide functional inference. We used the new algorithm to measure and reduce the influence of hub proteins on detecting functionally associated proteins. We used the annotations of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) as independent and unbiased benchmarks to evaluate our algorithms and methods within the human PPI network. We showed that, compared with the previous work from Samanta and Liang, our algorithm and methods developed in this study improved the overall quality of functional inferences for human proteins. By applying the algorithms to the human PPI network, we obtained 4,233 significant functional associations among 1,754 proteins. Further comparisons of their KEGG and GO annotations allowed us to assign 466 KEGG pathway annotations to 274 proteins and 123 GO annotations to 114 proteins with estimated false discovery rates of <21% for KEGG and <30% for GO. We clustered 1,729 proteins by their functional associations and made pathway analysis to identify several subclusters that are highly enriched in certain signaling pathways. Particularly, we performed a detailed analysis on a subcluster enriched in the transforming growth factor β signaling pathway (P<10-50) which is important in cell proliferation and tumorigenesis. Analysis of another four subclusters also suggested potential new players in six signaling pathways worthy of further experimental investigations. Our study gives clear insight into the common neighbor-based prediction scheme and provides a reliable method for large-scale functional annotations in this post-genomic era.

Isolation, characterization and functional analysis of the {\it xnf7\/} gene from {\it Xenopus laevis\/}

A fundamental problem in developmental biology concerns the mechanisms involved in the establishment of the embryonic axis. We are studying Xenopus nuclear factor 7 (xnf7) which we believe to be involved in dorsal-ventral patterning in Xenopus laevis. Xnf7 is a maternal gene product that is retained in the cytoplasm during early embryogenesis until the mid-blastula transition (MBT) when it reenters the nuclei. It is a member of a novel zinc finger proteins, the B-box family, consisting mainly of transcription factors and protooncogenes.^ The xnf7 gene is reexpressed during embryogenesis at the gastrula-neurula stage of development, with its zygotic expression limited to the central nervous system (CNS). In this study we showed that there are two different cDNAs coding for xnf7, xnf7-O and xnf7-B. They differ by 39 amino acid changes scattered throughout the cDNA. The expression of both forms of xnf7 is limited primarily to the central nervous system (CNS) and dorsal axial structures during later stages of embryogenesis.^ In order to study the spatial and temporal regulation of the gene, we screened a Xenopus genomic library using part of xnf7 cDNA as a probe. A genomic clone corresponding to the xnf7-O type was isolated, its 5$\sp\prime$ putative regulatory region sequenced, and its transcriptional initiation site mapped. The putative promoter region contained binding sites for Sp1, E2F, USF, a Pu box and AP1. CAT/xnf7 fusion genes were constructed containing various 5$\sp\prime$ deleted regions of the xnf7 promoter linked to a CAT (Chloramphenicol Acetyl Transferase) reporter vector. These constructs were injected into Xenopus oocytes and embryos to study the regions of the xnf7 promoter responsible for basal, temporal and spatial regulation of the gene. The activity of the fusion genes was measured by the conversion of chloramphenicol to its acetylated forms, and the spatial distribution of the transcripts by whole mount in situ hybridization. We showed that the elements involved in basal regulation of xnf7 lie within 121 basepairs upstream of the transcriptional inititiation site. A DNase I footprint analysis performed using oocyte extract showed that a E2F and 2 Sp1 sites were protected. During development, the fusion genes were expressed following the MBT, in accordance with the timing of the endogenous xnf7 gene. Spatially, the expression of the fusion gene containing 421 basepairs of the promoter was localized to the dorsal region of the embryo in a pattern that was almost identical to that detected with the endogenous transcripts. Therefore, the elements involved in spatial and temporal regulation of the xnf7 gene during development were contained within 421 basepairs upstream of the transcriptional initiation site. Future work will further define the elements involved in the spatial and temporal regulation and the trans-factors that interact with them. ^

Physical and functional mapping of a novel tumor suppressor locus for prostate cancer within chromosome 10p12.31-q11

Prostate cancer remains the second leading cause of male cancer deaths in the United States, yet the molecular mechanisms underlying this disease remain largely unknown. Cytogenetic and molecular analyses of prostate tumors suggest a consistent association with the loss of chromosome 10. Previously, we have defined a novel tumor suppressor locus PAC-1 within chromosome 10pter-q11. Introduction of the short arm of chromosome 10 into a prostatic adenocarcinoma cell line PC-3H resulted in dramatic tumor suppression and restoration of a programmed cell death pathway. Using a combined approach of comparative genomic hybridization and microsatellite analysis of PC-3H, I have identified a region of hemizygosity within 10p12-p15. This region has been shown to be involved in frequent loss of heterozygosity in gliomas and melanoma. To functionally dissect the region within chromosome 10p containing PAC-1, we developed a strategy of serial microcell fusion, a technique that allows the transfer of defined fragments of chromosome 10p into PC-3H. Serial microcell fusion was used to transfer defined 10p fragments into a mouse A9 fibrosarcoma cell line. Once characterized by FISH and microsatellite analyses, the 10p fragments were subsequently transferred into PC-3H to generate a panel of microcell hybrid clones containing overlapping deletions of chromosome 10p. In vivo and microsatellite analyses of these PC hybrids identified a small chromosome 10p fragment (an estimated 31 Mb in size inclusive of the centromere) that when transferred into the PC-3H background, resulted in significant tumor suppression and limited a region of functional tumor suppressor activity to chromosome 10p12.31-q11. This region coincides with a region of LOH demonstrated in prostate cancer. These studies demonstrate the utility of this approach as a powerful tool to limit regions of functional tumor suppressor activity. Furthermore, these data used in conjunction with data generated by the Human Genome Project lent a focused approach to identify candidate tumor suppressor genes involved in prostate cancer. ^

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

FANCM AND FAAP24 MAINTAIN GENOMIC STABILITY THROUGH COOPERATIVE AND UNIQUE FUNCTIONS

Fanconi anemia (FA) is a rare recessive genetic disease with an array of clinical manifestations including multiple congenital abnormalities, progressive bone marrow failure and profound cancer susceptibility. A hallmark of cells derived from FA patients is hypersensitivity to DNA interstrand crosslinking agents such as mitomycin C (MMC) and cisplatin, suggesting that FA- and FA-associated proteins play important roles in protecting cells from DNA interstrand crosslink (ICL) damage. Two genes involved in the FA pathway, FANCM and FAAP24, are of particular interest because they contain DNA interacting domains. However, there are no definitive patient mutations for these two genes, and the resulting lack of human genetic model system renders their functional studies difficult. In this study, I established isogenic human FANCM- and FAAP24-null mutants through homologous replacement-mediated gene targeting in HCT-116 cells, and systematically investigated the functions of FANCM and FAAP24 inchromosome stability, FA pathway activation, DNA damage checkpoint signaling, and ICL repair. I found that the FANCM-/-/FAAP24-/- double mutant was much more sensitive to DNA crosslinking agents than FANCM-/- and FAAP24-/- single mutants, suggesting that FANCM and FAAP24 possess epistatic as well as unique functions in response to ICL damage. I demonstrated that FANCM and FAAP24 coordinately support the activation of FA pathway by promoting chromatin localization of FA core complex and FANCD2 monoubiqutination. They also cooperatively function to suppress sister chromatid exchange and radial chromosome formation, likely by limiting crossovers in recombination repair. In addition, I defined novel non-overlapping functions of FANCM and FAAP24 in response to ICL damage. FAAP24 plays a major role in activating ICL-induced ATR-dependent checkpoint, which is independent of its interaction with FANCM. On the other hand, FANCM promotes recombination-independent ICL repair independently of FAAP24. Mechanistically, FANCM facilitates recruitment of nucleotide excision repair machinery and lesion bypass factors to ICL damage sites through its translocase activity. Collectively, my studies provide mechanistic insights into how genome integrity is both coordinately and independently protected by FANCM and FAAP24.

Identification and characterization of cellular genes differentially expressed in a subset of tumors with non -functional WT1

Wilms tumor (WT) or nephroblastoma is a genetically heterogeneous pediatric renal tumor that accounts for 6–7% of all childhood cancers in the U.S. WT1, located at 11p13, is the sole WT gene cloned to date. Additional genomic regions containing genes that play a role in the development of Wilms tumor include 11p15, 7p, 16q, 1p, 17q and 19q. This heterogeneity has made it extremely difficult to develop an understanding of the pathways involved in the development of WT, even in the 5–20% of tumors that show mutations at the WT1 locus. My research addresses this gap in our current comprehension of the development of WT. ^ I have used two complementary approaches to extend the current understanding of molecular changes involved in the development of WT. In order to minimize complexities due to genetic heterogeneity, I confined my analysis to the WT1 pathway by assessing those genetically defined tumors that carry WT1 mutations. WT1 encodes a zinc finger transcription factor, and in vitro studies have identified many genes that are potentially regulated in vivo by WT1. However, there is very little in vivo data that suggests that they are transcriptionally regulated endogenously by WT1. In one approach I assessed the role of WT1 in the in vivo regulation of PDGFA and IGF2, two genes that are strong contenders for endogenous regulation by WT1. Using primary tissue samples, I found no correlation between the level of RNA expression of WT1 with either PDGFA or IGF2, suggesting that WT1 does not play a critical role in their expression in either normal kidney or WT. ^ In a parallel strategy, using differential display analysis I compared global gene expression in a subset of tumors with known homozygous inactivating WT1 mutations (WT1-tumors) to the gene expression in a panel of appropriate control tissues (fetal kidney, normal kidney, rhabdoid tumor and pediatric renal cell carcinoma). Transcripts that are aberrantly expressed in this subset of Wilms tumors are candidates for endogenous transcriptional regulation by WT1 as well as for potentially functioning in the development of WT. By this approach I identified several differentially expressed transcripts. I further characterized two of these transcripts, identifying a candidate WT gene in the process. I then performed a detailed analysis of this WT candidate gene, which maps to 7p. Future studies will shed more light on the role of these differentially expressed genes in WT. ^