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

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.

ELUCIDATING FUNCTIONAL ROLES FOR MYOGENIN IN ADULT SKELETAL MUSCLE METABOLISM, EXERCISE CAPACITY, AND REGENERATION

The four basic helix-loop-helix myogenic transcription factors, myogenin, Myf5, MRF4, and MyoD are critical for embryonic skeletal muscle development. Myogenin is necessary for the terminal differentiation of myoblasts into myofibers during embryogenesis, but little is known about the roles played by myogenin in adult skeletal muscle function and metabolism. Furthermore, while metabolism is a well-studied physiological process, how it is regulated at the transcriptional level remains poorly understood. In this study, my aim was to determine the function of myogenin in adult skeletal muscle metabolism, exercise capacity, and regeneration. To investigate this, I utilized a mouse strain harboring the Myogflox allele and a Cre recombinase transgene, enabling the efficient deletion of myogenin in the adult mouse. Myogflox/flox mice were stressed physically through involuntary treadmill running and by breeding them with a strain harboring the Duchenne’s muscular dystrophy (DMDmdx) allele. Surprisingly, Myog-deleted animals exhibited an enhanced capacity for exercise, running farther and faster than their wild-type counterparts. Increased lactate production and utilization of glucose as a fuel source indicated that Myog-deleted animals exhibited an increased glycolytic flux. Hypoglycemic Myog-deleted mice no longer possessed the ability to outrun their wild-type counterparts, implying the ability of these animals to further deplete their glucose reserves confers their enhanced exercise capacity. Moreover, Myog-deleted mice exhibited an enhanced response to long-term exercise training. The mice developed a greater proportion of type 1 oxidative muscle fibers, and displayed increased levels of succinate dehydrogenase activity, indicative of increased oxidative metabolism. Mdx:Myog-deleted mice exhibited a similar phenotype, outperforming their mdx counterparts, although lagging behind wild-type animals. The morphology of muscle tissue from mdx:Myog-deleted mice appears to mimic that of mdx animals, indicating that myogenin is dispensable for adult skeletal muscle regeneration. Through global gene expression profiling and quantitative (q)RT-PCR, I identified a unique set of putative myogenin-dependent genes involved in regulating metabolic processes. These data suggest myogenin’s functions during adulthood are distinctly different than those during embryogenesis, and myogenin acts as a high-level transcription factor regulating metabolic activity in adult skeletal muscle.

THE ROLE OF MAP KINASES ON THE FUNCTIONAL HETEROGENEITY OF HUMAN CD8+ T CELL MATURATION SUBSETS

While prior studies have focused on naïve (CD45RA+CD27+) and early stage memory (CD45RA-CD27+) CD8+ T cells, late memory CD8+ T cells (CD45RA+CD27) have received less interest because this subset of T cells is generally recognized as effectors, which produce IFNγ (but no IL-2) and perforin. However, multiple studies suggest that late memory CD8+ T cells may provide inadequate protection in infectious diseases and cancer models. To better understand the unique function of late memory CD8+ T cells, I optimized multi-color flow cytometry techniques to assess the cytokine production of each human CD8+ T cell maturation subset. I demonstrated that late memory CD8+ T cells are the predominant producer of CC chemokines (e.g. MIP-1β), but rarely produce IL-2; therefore they do not co-produce IL-2/IFNγ (polyfunctionality), which has been shown to be critical for protective immunity against chronic viral infection. These data suggest that late memory CD8+ T cells are not just cytotoxic effectors, but may have unique functional properties. Determining the molecular signature of each CD8+ T cell maturation subset will help characterize the role of late memory CD8+ T cells. Prior studies suggest that ERK1 and ERK2 play a role in cytokine production including IL-2 in T cells. Therefore, I tested whether differential expression of ERK1 and ERK2 in CD8+ T cell maturation subsets contributes to their functional signature by a novel flow cytometry technique. I found that the expression of total ERK1, but not ERK2, is significantly diminished in late memory CD8+ T cells and that ERK1 expression is strongly associated with IL-2 production and CD28 expression. I also found that IL-2 production is increased in late memory CD8+ T cells by over-expressing ERK1. Collectively, these data suggest that ERK1 is required for IL-2 production in human CD8+ T cells. In summary, this dissertation demonstrated that ERK1 is down-regulated in human late memory CD8+ T cells, leading to decreased production of IL-2. The data in this dissertation also suggested that the functional heterogeneity in human CD8+ T cell maturation subsets results from their differential ERK1 expression.

Physical and functional interaction between the dopamine transporter and the synaptic vesicle protein synaptogyrin-3.

Uptake through the dopamine transporter (DAT) represents the primary mechanism used to terminate dopaminergic transmission in brain. Although it is well known that dopamine (DA) taken up by the transporter is used to replenish synaptic vesicle stores for subsequent release, the molecular details of this mechanism are not completely understood. Here, we identified the synaptic vesicle protein synaptogyrin-3 as a DAT interacting protein using the split ubiquitin system. This interaction was confirmed through coimmunoprecipitation experiments using heterologous cell lines and mouse brain. DAT and synaptogyrin-3 colocalized at presynaptic terminals from mouse striatum. Using fluorescence resonance energy transfer microscopy, we show that both proteins interact in live neurons. Pull-down assays with GST (glutathione S-transferase) proteins revealed that the cytoplasmic N termini of both DAT and synaptogyrin-3 are sufficient for this interaction. Furthermore, the N terminus of DAT is capable of binding purified synaptic vesicles from brain tissue. Functional assays revealed that synaptogyrin-3 expression correlated with DAT activity in PC12 and MN9D cells, but not in the non-neuronal HEK-293 cells. These changes were not attributed to changes in transporter cell surface levels or to direct effect of the protein-protein interaction. Instead, the synaptogyrin-3 effect on DAT activity was abolished in the presence of the vesicular monoamine transporter-2 (VMAT2) inhibitor reserpine, suggesting a dependence on the vesicular DA storage system. Finally, we provide evidence for a biochemical complex involving DAT, synaptogyrin-3, and VMAT2. Collectively, our data identify a novel interaction between DAT and synaptogyrin-3 and suggest a physical and functional link between DAT and the vesicular DA system.

Functional analysis of the amine substrate specificity domain of pepper tyramine and serotonin N-hydroxycinnamoyltransferases.

Pepper (Capsicum annuum) serotonin N-hydroxycinnamoyltransferase (SHT) catalyzes the synthesis of N-hydroxycinnamic acid amides of serotonin, including feruloylserotonin and p-coumaroylserotonin. To elucidate the domain or the key amino acid that determines the amine substrate specificity, we isolated a tyramine N-hydroxycinnamoyltransferase (THT) gene from pepper. Purified recombinant THT protein catalyzed the synthesis of N-hydroxycinnamic acid amides of tyramine, including feruloyltyramine and p-coumaroyltyramine, but did not accept serotonin as a substrate. Both the SHT and THT mRNAs were found to be expressed constitutively in all pepper organs. Pepper SHT and THT, which have primary sequences that are 78% identical, were used as models to investigate the structural determinants responsible for their distinct substrate specificities and other enzymatic properties. A series of chimeric genes was constructed by reciprocal exchange of DNA segments between the SHT and THT cDNAs. Functional characterization of the recombinant chimeric proteins revealed that the amino acid residues 129 to 165 of SHT and the corresponding residues 125 to 160 in THT are critical structural determinants for amine substrate specificity. Several amino acids are strongly implicated in the determination of amine substrate specificity, in which glycine-158 is involved in catalysis and amine substrate binding and tyrosine-149 plays a pivotal role in controlling amine substrate specificity between serotonin and tyramine in SHT. Furthermore, the indisputable role of tyrosine is corroborated by the THT-F145Y mutant that uses serotonin as the acyl acceptor. The results from the chimeras and the kinetic measurements will direct the creation of additional novel N-hydroxycinnamoyltransferases from the various N-hydroxycinnamoyltransferases found in nature.

Computational identification and functional validation of regulatory motifs in cartilage-expressed genes.

Chondrocyte gene regulation is important for the generation and maintenance of cartilage tissues. Several regulatory factors have been identified that play a role in chondrogenesis, including the positive transacting factors of the SOX family such as SOX9, SOX5, and SOX6, as well as negative transacting factors such as C/EBP and delta EF1. However, a complete understanding of the intricate regulatory network that governs the tissue-specific expression of cartilage genes is not yet available. We have taken a computational approach to identify cis-regulatory, transcription factor (TF) binding motifs in a set of cartilage characteristic genes to better define the transcriptional regulatory networks that regulate chondrogenesis. Our computational methods have identified several TFs, whose binding profiles are available in the TRANSFAC database, as important to chondrogenesis. In addition, a cartilage-specific SOX-binding profile was constructed and used to identify both known, and novel, functional paired SOX-binding motifs in chondrocyte genes. Using DNA pattern-recognition algorithms, we have also identified cis-regulatory elements for unknown TFs. We have validated our computational predictions through mutational analyses in cell transfection experiments. One novel regulatory motif, N1, found at high frequency in the COL2A1 promoter, was found to bind to chondrocyte nuclear proteins. Mutational analyses suggest that this motif binds a repressive factor that regulates basal levels of the COL2A1 promoter.

Experimental dissection and characterization of the functional domains in cardiac troponin C

Contraction of vertebrate cardiac muscle is regulated by the binding of Ca$\sp{2+}$ to the troponin C (cTnC) subunit of the troponin complex. In this study, we have used site-directed mutagenesis and a variety of assay techniques to explore the functional roles of regions in cTnC, including Ca$\sp{2+}$/Mg$\sp{2+}$-binding sites III and IV, the functionally inactive site I, the N-terminal helix, the N-terminal hydrophobic pocket and the two cysteine residues with regard to their ability to form disulfide bonds. Conversion of the first Ca$\sp{2+}$ ligand from Asp to Ala inactivated sites III and IV and decreased the apparent affinity of cTnC for the thin filament. Conversion of the second ligand from Asn to Ala also inactivated these sites in the free protein but Ca$\sp{2+}$-binding was recovered upon association with troponin I and troponin T. The Ca$\sp{2+}$-concentrations required for tight thin filament-binding by proteins containing second-ligand mutations were significantly greater than that required for the wild-type protein. Mutation of site I such that the primary sequence was that of an active site with the first Ca$\sp{2+}$ ligand changed from Asp to Ala resulted in a 70% decrease in maximal Ca$\sp{2\sp+}$ dependent ATPase activity in both cardiac and fast skeletal myofibrils. Thus, the primary sequence of the inactive site I in cTnC is functionally important. Major changes in the sequence of the N-terminus had little effect on the ability of cTnC to recover maximal activity but deletion of the first nine residues resulted in a 60 to 80% decrease in maximal activity with only a minor decrease in the pCa$\sb{50}$ of activation, suggesting that the N-terminal helix must be present but that a specific sequence is not required. The formation of an inter- or intramolecular disulfide bonds caused the exposure of hydrophobic surfaces on cTnC and rendered the protein Ca$\sp{2+}$ independent. Finally, elution patterns from a hydrophobic interactions column suggest that cTnC undergoes a significant change in hydrophobicity upon Ca$\sp{2+}$ binding, the majority of which is caused by site II. These latter data show an interesting correlation between exposure of hydrophobic surfaces on and activation of cTnC. Overall, these results represent significant progress toward the elucidation of the functional roles of a variety of structural regions in cTnC. ^

A cladistic analysis of apolipoprotein B gene variation and plasma lipid and apolipoprotein B levels, using familial data

Any functionally important mutation is embedded in an evolutionary matrix of other mutations. Cladistic analysis, based on this, is a method of investigating gene effects using a haplotype phylogeny to define a set of tests which localize causal mutations to branches of the phylogeny. Previous implementations of cladistic analysis have not addressed the issue of analyzing data from related individuals, though in human studies, family data are usually needed to obtain unambiguous haplotypes. In this study, a method of cladistic analysis is described in which haplotype effects are parameterized in a linear model which accounts for familial correlations. The method was used to study the effect of apolipoprotein (Apo) B gene variation on total-, LDL-, and HDL-cholesterol, triglyceride, and Apo B levels in 121 French families. Five polymorphisms defined Apo B haplotypes: the signal peptide Insertion/deletion, Bsp 1286I, XbaI, MspI, and EcoRI. Eleven haplotypes were found, and a haplotype phylogeny was constructed and used to define a set of tests of haplotype effects on lipid and apo B levels.^ This new method of cladistic analysis, the parametric method, found significant effects for single haplotypes for all variables. For HDL-cholesterol, 3 clusters of evolutionarily-related haplotypes affecting levels were found. Haplotype effects accounted for about 10% of the genetic variance of triglyceride and HDL-cholesterol levels. The results of the parametric method were compared to those of a method of cladistic analysis based on permutational testing. The permutational method detected fewer haplotype effects, even when modified to account for correlations within families. Simulation studies exploring these differences found evidence of systematic errors in the permutational method due to the process by which haplotype groups were selected for testing.^ The applicability of cladistic analysis to human data was shown. The parametric method is suggested as an improvement over the permutational method. This study has identified candidate haplotypes for sequence comparisons in order to locate the functional mutations in the Apo B gene which may influence plasma lipid levels. ^

Functional studies on DNA double-strand break repair proteins (MmRad51, Ku80) in mice

RecA in Escherichia coli and it's homologue, ScRad51 in Saccharomyces cerevisiae, play important roles in recombinational repair. ScRad51 homologues have been discovered in a wide range of organisms including Schizosaccharomyces pombe, lily, chicken, mouse and human. To date there is no direct evidence to describe that mouse Rad51(MmRad51) is involved in DNA double-strand break repair. In order to elucidate the role of MmRad51 in vivo, it was mutated by the embryonic stem (ES) cell/gene targeting technology in mice. The mutant embryos arrested in development shortly after implantation. There was a decrease in cell proliferation followed by programmed cell death, and trophectoderm-derived cells were sensitive to $\gamma$-radiation. Severe chromosome loss was observed in most mitotically dividing cells. The mutant embryos lived longer and developed further in a p53 mutant background; however, double-mutant embryonic fibroblasts failed to proliferate in tissue culture, reflecting the embryos limited life span. Based on these data, MmRad51 repairs DNA damage induced by $\gamma$-radiation, is needed to maintain euplody, and plays an important role in proliferating cells.^ Ku is a heterodimer of 70 and 80 kDs subunit, which binds to DNA ends and other altered DNA structures such as hairpins, nicks, and gaps. In addition, Ku is required for DNA-PK activity through a direct association. Although the biochemical properties of Ku and DNA-PKcs have been characterized in cells, their physiological functions are not clear. In order to understand the function of Ku in vivo, we generated mice homozygous for a mutation of the Ku80 gene. Ku80-deficient mice, like scid mice, showed severe immunodeficiency due to a impairment of V(D)J recombination. Mutant mice were semiviable and runted, cells derived from mutant embryos displayed hypersensitivity to $\gamma$-radiation, a decreased growth rate, a slow entry into S phase, altered colony size distributions, and a short life span. Based on these results, mutant cells and mice appeared to prematurely age. ^

Functional analysis of MRF4 during mouse embryogenesis

MRF4 is one of four skeletal muscle specific regulatory genes, (the other three genes being MyoD, myf5, and myogenin), each of which has the unique ability to orchestrate an entire program of muscle-specific transcription when introduced into diverse cell types. These findings have led to the notion that these factors function as master regulators of muscle cell fate. Analysis of mice lacking MyoD, myf5, and myogenin have further defined their roles in the commitment and differentiation of myotomal progenitor cells. Current data strongly supports the model that MyoD and myf5 share functional redundancy in determining the muscle cell lineage, while myogenin acts downstream of MyoD and myf5, to initiate myoblast differentiation. Unlike other myogenic bHLH genes, MRF4 is expressed predominantly in the adult, suggesting that it may function to regulate adult muscle maturation and maintenance. To test this hypothesis and to eventually incorporate MRF4 into a general model for muscle specification, differentiation, maturation and maintenance, I deleted the MRF4 gene. MRF4-null mice are viable and fertile, however, they show mild rib anomalies. In addition, the expression of myogenin is dramatically upregulated only in the adult, suggesting that myogenin may compensate for the loss of MRF4 in the adult, and MRF4 may normally suppress the expression of myogenin after birth. MRF4 is also required during muscle regeneration after injury.^ To determine the degree of genetic redundancy between MRF4-myogenin; and MRF4-MyoD, I crossed the MRF4-null mice with MyoD- and myogenin-null mice respectively. There are no additional muscle phenotypes in double-null progeny from a MRF4 and myogenin cross, suggesting that the existence of residual fibers in myogenin-null mice is not due to the presence of MRF4. MRF4 expression also cannot account for the ability of myogenin-null myoblasts to differentiate in vitro. However, the combination of the MRF4-null mutation with the myogenin-null mutation results in a novel rib phenotype. This result suggests that MRF4 modifies the myogenin-null rib phenotype, and MRF4 and myogenin play redundant roles in rib development.^ MRF4 also shares dosage effects with MyoD during mouse development. (MyoD+/$-$;MRF4$-$/$-$)mice are fertile and viable, while (MyoD$-$/$-$;MRF4+/$-$) mice die between birth and two weeks after birth, and have a small skeletal structure. The double homozygous mice for MRF4 and MyoD mutations are embryonic lethal and die at around E10.5. These results suggest that MRF4 and MyoD share overlapping functions during mouse embryogenesis. ^

Transcriptional regulation and functional analysis of the human Wilms' tumor 1 gene

The Wilms' tumor 1 gene (WT1) encodes a zinc-finger transcription factor and is expressed in urogenital, hematopoietic and other tissues. It is expressed in a temporal and spatial manner in both embryonic and adult stages. To obtain a better understanding of the biological function of WT1, we studied two aspects of WT1 regulation: one is the identification of tissue-specific cis-regulatory elements that regulate its expression, the other is the downstream genes which are modulated by WT1.^ My studies indicate that in addition to the promoter, other regulatory elements are required for the tissue specific expression of this gene. A 259-bp hematopoietic specific enhancer in intron 3 of the WT1 gene increased the transcriptional activity of the WT1 promoter by 8- to 10-fold in K562 and HL60 cells. Sequence analysis revealed both GATA and c-Myb motifs in the enhancer fragment. Mutation of the GATA motif decreased the enhancer activity by 60% in K562 cells. Electrophoretic mobility shift assays showed that both GATA-1 and GATA-2 proteins in K562 nuclear extracts bind to this motif. Cotransfection of the enhancer containing reporter construct with a GATA-1 or GATA-2 expression vector showed that both GATA-1 and GATA-2 transactivated this enhancer, increasing the CAT reporter activity 10-15 fold and 5-fold respectively. Similar analysis of the c-Myb motif by cotransfection with the enhancer CAT reporter construct and a c-Myb expression vector showed that c-Myb transactivated the enhancer by 5-fold. A DNase I-hypersensitive site has been identified in the 258 bp enhancer region. These data suggest that GATA-1 and c-Myb are responsible for the activity of this enhancer in hematopoietic cells and may bind to the enhancer in vivo. In the process of searching for cis-regulatory elements in transgenic mice, we have identified a 1.0 kb fragment that is 50 kb downstream from the promoter and is required for the central nervous system expression of WT1.^ In the search for downstream target genes of WT1, we noted that the proto-oncogene N-myc is coexpressed with the tumor suppressor gene WT1 in the developing kidney and is overexpressed in many Wilms' tumors. Sequence analysis revealed eleven consensus WT1 binding sites located in the 1 kb mouse N-myc promoter. We further showed that the N-myc promoter was down-regulated by WT1 in transient transfection assays. Electrophoretic mobility shift assays showed that oligonucleotides containing the WT1 motifs could bind WT1 protein. Furthermore, a Denys-Drash syndrome mutant of WT1, R394W, that has a mutation in the DNA binding domain, failed to repress the N-myc promoter. This suggests that the repression of the N-myc promoter is mediated by DNA binding of WT1. This finding helps to elucidate the relationship of WT1 and N-myc in tumorigenesis and renal development. ^

Functional analysis of histones H3 and H4 in Tup1 repression and in GCN5 activation

Tup1 forms a complex with Ssn6 in yeast. Ssn6-Tup1 complex is recruited via direct interactions with specific DNA binding proteins to a specific promoter region and mediates repression of several sets of genes including a-cell specific genes (asg) in $\alpha$ cells. It has been shown that repression of asgs also requires histone H4 and that Tup1 can directly interact with H3 and H4 in vitro. To address whether histone H3 is required for the repression of asgs, I have examined the effect of H3 and H4 mutations on the expression of a $\alpha$2-controlled LacZ reporter. Assay of $\beta$-glactosidase shows that mutations in either H3 or H4 cause a weak derepression of the reporter gene. Some double mutations result in a stronger derepression, while others do not. The H3 N-terminal deletion also leads to a slightly decreased expression of the reporter gene in $\alpha$ cells. Our data suggest that the N-termini of both H3 and H4 are cooperatively involved in the repression of a-cell specific genes in $\alpha$ cells, possibly through their interaction with Tup1.^ GCN5 was originally identified as a transcriptional regulator required to activate a subset of genes in yeast. Recently, it has been shown that GCN5 encodes the catalytic subunit of a nuclear histone acetyltransferase, providing the first direct link between histone acetylation and gene transcription. Recombinant Gcn5p (rGcn5p) exhibits a limited substrate specificity in vitro. However, neither the specificity of this enzyme in vivo nor the importance of particular acetylated residues to transcription or cell growth are well defined. In order to define the sites of histone acetylation mediated by Gcn5p in vivo and assess the significance of histone acetylation, more than 30 yeast strains have been constructed to bear specific H3 and/or H4 mutations in the presence or absence of GCN5 function. Our genetic data suggest that Gcn5p may have additional targets in vivo that are not identified as the targets of rGcn5p by previous studies. Western analysis using antibodies specifically recognizing particular acetylated isoforms of H3 and H4 led us to conclude that Gcn5p is necessary for full acetylation of multiple sites in both H3 and H4 in vivo. Consistent with these observations, rGcn5p still acetylates histones H3 and H4 bearing mutations either in H3 K14 or H4 K8,16, sites previously identified as the targets of acetylation by rGcn5p in H3 and H4. Our data also demonstrated that Gcn5p-mediated acetylation events are important for normal progression of the cell cycle and for transcriptional activation. Furthermore, a critical overall level of acetylation is essential for cell viability. ^

Transcription factor SEF: Purification and functional characterization

Cytokine-induced transcription of the serum amyloid A3 (SAA3) gene promoter requires a transcriptional enhancer that contains three functional elements: two C/EBP-binding sites and a third site that interacts with a constitutively expressed transcription factor, SAA3 enhancer factor (SEF). Deletion or site-specific mutations in the SEF-binding site drastically reduced SAA3 promoter activity, strongly suggesting that SEF is important in SAA3 promoter function. To further elucidate its role in the regulation of the SAA3 gene, we purified SEF from HeLa cell nuclear extracts to near homogeneity by using conventional liquid chromatography and DNA-affinity chromatography. Ultraviolet cross-linking and Southwestern experiments indicated that SEF consisted of a single polypeptide with an apparent molecular mass of 65 kDa. Protein sequencing, oligonucleotide competition and antibody supershift experiments identified SEF as transcription factor LBP-1c/CP2/LSF. Cotransfection of SEF expression plasmid with SAA3-luciferase reporter resulted in 3- to 5-fold activation of SAA3 promoter. Interestingly, when SEF-transfected cells were treated with either conditioned medium (CM) or interleukin (IL) 1, the SAA3 promoter was synergistically activated in a dose-dependent manner. Furthermore, when SEF-binding site was mutated, the response of SAA3 promoter to IL-1 or CM stimulation was abolished or drastically decreased, suggesting that SEF may functionally cooperate with an IL-1-inducible transcription factor. Indeed, our functional studies showed that NFκB is a key transcription factor that mediates the IL-1-induced expression of SAA3 gene, and that SEF can synergize with NFκBp65 to activate SAA3 promoter. By coimmunoprecipitation experiments, we found that SEF could specifically interact with NFκBp65, and that the association of these two factors was enhanced upon IL-1 and CM stimulation. This suggests that the molecular basis for the functional synergy between SEF and NFκB may be due to the ability of SEF to physically interact with NPκB. In addition to its interaction with SEF, NFκB-dependent activation also requires the weak κB site in the C element and its interaction with C/EBP. Besides its role in regulating SAA3 gene expression, we provide evidence that SEF could also bind in a sequence-specific manner to the promoters of α2-macroglobulin, Aα fibrinogen, and 6–16 genes and to an intronic enhancer of the human Wilm's tumor 1 gene, suggesting a functional role in the regulation of these genes. By coimmunoprecipitation experiments, we determined that SEF could specifically associate with both Stat3 and Stat2 upon cytokine stimulation. To examine the functional roles of such interactions, we evaluated the effects of SEF on the transcriptional regulation of two reporter genes: Aα fibrinogen and 6–16, which are IL-6- and interferon-α-responsive, respectively. Our results showed that cotransfection of SEF expression plasmid can activate the expression of Aα fibrinogen gene and 6–16 gene. Moreover, SEF can dramatically enhance the interferon-α-induced expression of 6–16 gene and IL-6-induced expression of Aα fibrinogen gene, suggesting that SEF may functionally cooperate with ISGF3 and Stat3 to mediate interferon-α and IL-6 signaling. ^ Our findings that SEF can interact with multiple cytokine-inducible transcription factors to mediate the expression of target genes open a new avenue of investigation of cooperative transcriptional regulation of gene expression, and should further our understanding of differential gene expression in response to a specific stimulus. In summary, our data provide evidence that SEF can mediate the signaling of different cytokines by interacting with various cytokine-inducible transcription factors. ^

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