Architectural Proteins and Pluripotency Factors Cooperate to Orchestrate the Transcriptional Response of hESCs to Temperature Stress

Cells respond to temperature stress via up- and downregulation of hundreds of genes. This process is thought to be regulated by the heat shock factor HSF1, which controls the release of RNAPII from promoter-proximal pausing. Here, we analyze the events taking place in hESCs upstream of RNAPII release. We find that temperature stress results in the activation or decommissioning of thousands of enhancers. This process involves alterations in the occupancy of transcription factors HSF1, AP-1, NANOG, KLF4, and OCT4 accompanied by nucleosome remodeling by BRG1 and changes in H3K27ac. Furthermore, redistribution of RAD21 and CTCF results in the formation and disassembly of interactions mediated by these two proteins. These alterations tether and untether enhancers to their cognate promoters or refashion insulated neighborhoods, thus transforming the landscape of enhancer-promoter interactions. Details of the 3D interactome remodeling process support loop extrusion initiating at random sites as a mechanism for the establishment of CTCF/cohesin loops.

Analysis of Chromatin Interactions Mediated by Specific Architectural Proteins in Drosophila Cells

Chromosome conformation capture assays have been established, modified, and enhanced for over a decade with the purpose of studying nuclear organization. A recently published method uses in situ Hi-C followed by chromatin immunoprecipitation (HiChIP) to enrich the overall yield of significant genome-wide interactions mediated by a specific protein. Here we applied a modified version of the HiChIP protocol to retrieve the significant contacts mediated by architectural protein CP190 in D. melanogaster cells.

Different enhancer classes in Drosophila bind distinct architectural proteins and mediate unique chromatin interactions and 3D architecture

Eukaryotic gene expression is regulated by enhancer–promoter interactions but the molecular mechanisms that govern specificity have remained elusive. Genome-wide studies utilizing STARR-seq identified two enhancer classes in Drosophila that interact with different core promoters: housekeeping enhancers (hkCP) and developmental enhancers (dCP). We hypothesized that the two enhancer classes are occupied by distinct architectural proteins, affecting their enhancer–promoter contacts. By evaluating ChIP-seq occupancy of architectural proteins, typical enhancer-associated proteins, and histone modifications, we determine that both enhancer classes are enriched for RNA Polymerase II, CBP, and architectural proteins but there are also distinctions. hkCP enhancers contain H3K4me3 and exclusively bind Cap-H2, Chromator, DREF and Z4, whereas dCP enhancers contain H3K4me1 and are more enriched for Rad21 and Fs(1)h-L. Additionally, we map the interactions of each enhancer class utilizing a Hi-C dataset with <1 kb resolution. Results suggest that hkCP enhancers are more likely to form multi-TSS interaction networks and be associated with topologically associating domain (TAD) borders, while dCP enhancers are more often bound to one or two TSSs and are enriched at chromatin loop anchors. The data support a model suggesting that the unique architectural protein occupancy within enhancers is one contributor to enhancer–promoter interaction specificity.

Minute-Made Data Analysis: Tools for Rapid Interrogation of Hi-C Contacts

Juicer and Juicebox, described by Durand et al. (2016a, 2016b), are two new tools for fast and reliable processing of Hi-C data, providing approaches for read processing, multiple normalization schemes, feature annotation, and dynamic browsing of chromatin contacts, thus reducing arduous Hi-C analysis into an easy yet flexible pipeline.

Evolutionarily Conserved Principles Predict 3D Chromatin Organization

Topologically associating domains (TADs), CTCF loop domains, and A/B compartments have been identified as important structural and functional components of 3D chromatin organization, yet the relationship between these features is not well understood. Using high-resolution Hi-C and HiChIP, we show that Drosophila chromatin is organized into domains we term compartmental domains that correspond precisely with A/B compartments at high resolution. We find that transcriptional state is a major predictor of Hi-C contact maps in several eukaryotes tested, including C. elegans and A. thaliana. Architectural proteins insulate compartmental domains by reducing interaction frequencies between neighboring regions in Drosophila, but CTCF loops do not play a distinct role in this organism. In mammals, compartmental domains exist alongside CTCF loop domains to form topological domains. The results suggest that compartmental domains are responsible for domain structure in all eukaryotes, with CTCF playing an important role in domain formation in mammals.

Epigenetic regulation of Plasmodium falciparum clonally variant gene expression during development in Anopheles gambiae

P. falciparum phenotypic plasticity is linked to the variant expression of clonal multigene families such as the var genes. We have examined changes in transcription and histone modifications that occur during sporogonic development of P. falciparum in the mosquito host. All var genes are silenced or transcribed at low levels in blood stages (gametocyte/ring) of the parasite in the human host. After infection of mosquitoes, a single var gene is selected for expression in the oocyst, and transcription of this gene increases dramatically in the sporozoite. The same PF3D7_1255200 var gene was activated in 4 different experimental infections. Transcription of this var gene during parasite development in the mosquito correlates with the presence of low levels of H3K9me3 at the binding site for the PF3D7_1466400 AP2 transcription factor. This chromatin state in the sporozoite also correlates with the expression of an antisense long non-coding RNA (lncRNA) that has previously been shown to promote var gene transcription during the intraerythrocytic cycle in vitro. Expression of both the sense protein-coding transcript and the antisense lncRNA increase dramatically in sporozoites. The findings suggest a complex process for the activation of a single particular var gene that involves AP2 transcription factors and lncRNAs.

Capturing native interactions: intrinsic methods to study chromatin conformation

The 3D organization of chromatin controls gene expression through spatial interactions between genomic loci. FISH and 3C‐based methods that are commonly used to study chromatin organization utilize chemical crosslinking, a step that may introduce biases in detectable chromatin interactions. In their recent study, Papantonis and colleagues (Brant et al, 2016) developed alternative new methods of detecting chromatin contacts without the use of chemical crosslinking agents. These tools increase the resolution and confidence at which interactions can be identified, and may be informative for chromatin interaction dynamics.

Long non-coding RNA produced by RNA polymerase V determines boundaries of heterochromatin

RNA-mediated transcriptional gene silencing is a conserved process where small RNAs target transposons and other sequences for repression by establishing chromatin modifications. A central element of this process are long non-coding RNAs (lncRNA), which in Arabidopsis thaliana are produced by a specialized RNA polymerase known as Pol V. Here we show that non-coding transcription by Pol V is controlled by preexisting chromatin modifications located within the transcribed regions. Most Pol V transcripts are associated with AGO4 but are not sliced by AGO4. Pol V-dependent DNA methylation is established on both strands of DNA and is tightly restricted to Pol V-transcribed regions. This indicates that chromatin modifications are established in close proximity to Pol V. Finally, Pol V transcription is preferentially enriched on edges of silenced transposable elements, where Pol V transcribes into TEs. We propose that Pol V may play an important role in the determination of heterochromatin boundaries.

DOI:http://dx.doi.org/10.7554/eLife.19092.001

The three-dimensional genome: principles and roles of long-distance interactions

The linear sequence of eukaryotic genomes is arranged in a specific manner within the three-dimensional nuclear space. Interactions between distant sites partition the genome into domains of highly associating chromatin. Interaction domains are found in many organisms, but their properties and the principles governing their establishment vary between different species. Topologically associating domains (TADs) extending over large genomic regions are found in mammals and Drosophila melanogaster, whereas other types of contact domains exist in lower eukaryotes. Here we review recent studies that explore the mechanisms by which long distance chromatin interactions determine the 3D organization of the genome and the relationship between this organization and the establishment of specific patterns of gene expression.

A Dicer-Independent Route for Biogenesis of siRNAs that Direct DNA Methylation in Arabidopsis

DNA methylation directed by 24-nucleotide (nt) small interfering RNAs (siRNAs) plays critical roles in gene regulation and transposon silencing in Arabidopsis. 24-nt siRNAs are known to be processed from double-stranded RNAs by Dicer-like 3 (DCL3) and loaded into the effector Argonaute 4 (AGO4). Here we report a distinct class of siRNAs independent of DCLs (sidRNAs). sidRNAs are present as ladders of ∼ 20-60 nt in length, often having the same 5' ends but differing in 3' ends by 1-nt steps. We further show that sidRNAs are associated with AGO4 and capable of directing DNA methylation. Finally we show that sidRNA production depends on distributive 3'-5' exonucleases. Our findings suggest an alternative route for siRNA biogenesis. Precursor transcripts are bound by AGO4 and subsequently subjected to 3'-5' exonucleolytic trimming for maturation. We propose that sidRNAs generated through this route are the initial triggers of de novo DNA methylation.

RNA-directed DNA methylation requires stepwise binding of silencing factors to long non-coding RNA

Ribonucleic acid-mediated transcriptional gene silencing (known as RNA-directed DNA methylation, or RdDM, in Arabidopsis thaliana) is important for influencing gene expression and the inhibition of transposons by the deposition of repressive chromatin marks such as histone modifications and DNA methylation. A key event in de novo methylation of DNA by RdDM is the production of long non-coding RNA (lncRNA) by RNA polymerase V (Pol V). Little is known about the events that connect Pol V transcription to the establishment of repressive chromatin modifications. Using RNA immunoprecipitation, we elucidated the order of events downstream of lncRNA production and discovered interdependency between lncRNA-associated proteins. We found that the effector protein ARGONAUTE4 (AGO4) binds lncRNA independent of the RNA-binding protein INVOLVED IN DE NOVO2 (IDN2). In contrast, IDN2 binds lncRNA in an AGO4-dependent manner. We further found that the de novo DNA methyltransferase DOMAINS REARRANGED METHYLTRANSFERASE2 (DRM2) also associates with lncRNA produced by Pol V and that this event depends on AGO4 and IDN2. We propose a model where the silencing proteins AGO4, IDN2 and DRM2 bind to lncRNA in a stepwise manner, resulting in DNA methylation of RdDM target loci.

Analysis of long non-coding RNAs produced by a specialized RNA polymerase in Arabidopsis thaliana

Long non-coding RNAs (lncRNAs) play important roles in several processes including control of gene expression. In Arabidopsis thaliana, a class of lncRNAs is produced by a specialized RNA Polymerase V (Pol V), which is involved in controlling genome activity by transcriptional gene silencing. lncRNAs produced by Pol V have been proposed to serve as scaffolds for binding of several silencing factors which further mediate the establishment of repressive chromatin modifications. We present methods for discovery and characterization of lncRNAs produced by Pol V. Chromatin Immunoprecipitation coupled with deep sequencing (ChIP-seq) allows discovery of genomic regions bound by proteins in a manner dependent on either Pol V or transcripts produced by Pol V. RNA Immunoprecipitation (RIP) allows testing lncRNA-protein interactions at identified loci. Finally, real-time RT-PCR allows detection of low abundance Pol V transcripts from total RNA. These methods may be more broadly applied to discovery and characterization of RNAs produced by distinct RNA Polymerases.

RNA polymerase V targets transcriptional silencing components to promoters of protein-coding genes

Transcriptional gene silencing controls transposons and other repetitive elements through RNA-directed DNA methylation (RdDM) and heterochromatin formation. A key component of the Arabidopsis RdDM pathway is ARGONAUTE4 (AGO4), which associates with siRNAs to mediate DNA methylation. Here, we show that AGO4 preferentially targets transposable elements embedded within promoters of protein-coding genes. This pattern of AGO4 binding cannot be simply explained by the sequences of AGO4-bound siRNAs; instead, AGO4 binding to specific gene promoters is also mediated by long non-coding RNAs (lncRNAs) produced by RNA polymerase V. lncRNA-mediated AGO4 binding to gene promoters directs asymmetric DNA methylation to these genomic regions and is involved in regulating the expression of targeted genes. Finally, AGO4 binding overlaps sites of DNA methylation affected by the biotic stress response. Based on these findings, we propose that the targets of AGO4-directed RdDM are regulatory units responsible for controlling gene expression under specific environmental conditions.

Spatial and functional relationships among Pol V-associated loci, Pol IV-dependent siRNAs, and cytosine methylation in the Arabidopsis epigenome

Multisubunit RNA polymerases IV and V (Pols IV and V) mediate RNA-directed DNA methylation and transcriptional silencing of retrotransposons and heterochromatic repeats in plants. We identified genomic sites of Pol V occupancy in parallel with siRNA deep sequencing and methylcytosine mapping, comparing wild-type plants with mutants defective for Pol IV, Pol V, or both Pols IV and V. Approximately 60% of Pol V-associated regions encompass regions of 24-nucleotide (nt) siRNA complementarity and cytosine methylation, consistent with cytosine methylation being guided by base-pairing of Pol IV-dependent siRNAs with Pol V transcripts. However, 27% of Pol V peaks do not overlap sites of 24-nt siRNA biogenesis or cytosine methylation, indicating that Pol V alone does not specify sites of cytosine methylation. Surprisingly, the number of methylated CHH motifs, a hallmark of RNA-directed de novo methylation, is similar in wild-type plants and Pol IV or Pol V mutants. In the mutants, methylation is lost at 50%-60% of the CHH sites that are methylated in the wild type but is gained at new CHH positions, primarily in pericentromeric regions. These results indicate that Pol IV and Pol V are not required for cytosine methyltransferase activity but shape the epigenome by guiding CHH methylation to specific genomic sites.

Clinical, laboratory and genetic markers associated with erosions and remission in patients with early inflammatory arthritis: a prospective cohort study

We investigated the relationship between clinical, laboratory and genetic markers and outcome measures in 159 patients with recent onset of inflammatory arthritis (IA). The majority of patients were managed in community-based rheumatology practice. Median duration of arthritis at baseline was 3 months with median follow-up of 4.0 years (range 0-10). Markers of disease activity and 1987 ACR criteria for rheumatoid arthritis (RA) were estimated every 6 months for the first 2 years and annually thereafter. Presence of shared epitopes (SE) was established by PCR-based method. Main outcome variables were attainment of remission and presence of erosions on X-rays of hands and feet at 3 years. Remission was seen in 34.3% of patients and was independently related to age 60 and older (odds ratio (OR) 3.2; 95% confidence interval (CI), 1.2-8.7) and inversely to the presence of rheumatoid factor (RF) (OR 8.3; 95% CI, 3.2-21.3 for persistent arthritis). Patients with two SE were likely to have persistent arthritis (P=0.006), but this was not significant when corrected for RF. Independent predictors for erosions at 3 years were RF (OR 7.5; 95% CI, 1.9-29.5) and area under the curve for number of swollen joints (OR 1.08; 95% CI, 1.02-1.16). SE status was not predictive of erosions at 3 years (OR 1.6; 95% CI, 0.7-3.7). In univariate analysis, patients possessing DERAA motif on DRB1 were less likely to have erosive disease than without this motif at 4 years (OR 0.21; 95% CI, 0.0-0.9, P=0.037) but this finding was partly explained by adjusting for RF (adjusted OR 0.24; 95% CI 0.04-1.37). In this study of recent onset IA, active disease and RF were associated with poor outcome. Whilst SE did not predict erosive disease, patients with DERAA motif may be protected against erosions whilst the presence of two SE alleles suggests persistence of arthritis.

Autoantibodies to the islet cell antigen SOX-13 are associated with duration but not type of diabetes

AIMS

The autoantigen SOX-13 of the SRY-related high mobility group box is a low-frequency reactant in sera from patients with Type 1 diabetes. We further investigated the potential diagnostic role of anti-SOX-13, and in particular its ability to distinguish Type 1 from Type 2 diabetes, in two large, well-characterized cohorts.

METHODS

SOX-13 autoantibody status was ascertained using a radioimmunoprecipitation assay in (i) a random sample of 546 participants in an Australian community-based study (the Fremantle Diabetes Study; FDS) of whom 119 had Type 1 and 427 Type 2 diabetes, and (ii) a sample of 333 subjects with Type 2 diabetes from the United Kingdom Prospective Diabetes Study (UKPDS) stratified by age, anti-glutamic acid decarboxylase (GAD) and islet cell antibody (ICA) status, and requirement for insulin therapy within 6 years of diagnosis.

RESULTS

The frequencies of anti-SOX-13 in the FDS subjects were 16.0% and 14.8% for Type 1 and Type 2 patients, respectively, and levels were similar. In the UKPDS subjects, the frequency was 4.5%. In a logistic regression model involving demographic, anthropometric and metabolic variables, only diabetes duration was significantly associated with anti-SOX-13 positivity, especially for duration > 5 years (P < 0.002). When the coexistence of autoantibodies was assessed in the two study samples, there were no significant associations between anti-SOX-13 and ICA, anti-GAD or ICA512/IA-2.

CONCLUSIONS

Whilst the frequency of anti-SOX-13 may be increased in some populations of diabetic patients, this reactivity does not usefully distinguish Type 1 from Type 2 diabetes. However, the association with diabetes duration suggests that anti-SOX-13 may be a non-specific marker of tissue damage associated with chronic hyperglycaemia.

Antimitochondrial autoantibodies in saliva and sera from patients with primary biliary cirrhosis

BACKGROUND AND AIMS

Primary biliary cirrhosis (PBC) is a cholestatic autoimmune liver disease characterized by antimitochondrial autoantibodies (AMA) in serum, for which the reactants are E2 subunits of the three 2-oxoacid dehydrogenase (2-OAD) enzymes, particularly pyruvate dehydrogenase complex (PDC-E2). Some 70% of patients with PBC have a coexisting autoimmune disease including Sjögren's syndrome. We aimed to ascertain the frequency and isotype of AMA in saliva in PBC.

METHODS

Serum and saliva from 12 patients with PBC were tested for AMA by immunoblotting on bovine heart mitochondria, and by an automated microassay based on inhibition of the enzymatic activity of PDC.

RESULTS

Autoantibodies of the immunoglobulin (Ig)G, IgM, and IgA immunoglobulin isotypes against the E2 subunits of 2-OAD enzymes were demonstrable in PBC in serum (12 of 12 cases) and saliva (nine of 12 cases). Salivary autoantibodies, like serum autoantibodies, were predominantly reactive with PDC and of the IgG isotype. Results for serum and saliva corresponded closely with regard to reactivity with individual enzymes of the 2-OAD enzyme family, and to the autoantibody isotype that was predominantly expressed, and also in the capacity to inhibit the enzymatic activity of PDC.

CONCLUSIONS

The presence of AMA in saliva to 2-OAD enzymes indicates that salivary glands could participate in the pathogenetic process of PBC. The detection of salivary AMA by a semi-automated enzyme inhibition assay offers possibilities for rapid population screening for detection of preclinical PBC among at-risk individuals, middle-aged to older women.

Antibodies to SOX13 (ICA12) are associated with type 1 diabetes

SOX13 is an islet cell autoantigen (ICA12), identified by antibody screening of an islet cDNA library, using sera from patients with Type 1 diabetes. We ascertained the frequency of antibody reactivity to SOX13 and compared it with other Type 1 diabetes autoantibody reactivities. Antibodies were measured by radioimmunoprecipitation (RIP) using (35) S labelled SOX13 expressed in rabbit reticulocyte lysate. Sera from 109 subjects with Type 1 diabetes, 29 with Type 2 diabetes, 144 with other autoimmune diseases and from 201 controls were tested for anti-SOX13, and results were compared with the frequency of antibodies to glutamic acid decarboxylase (anti-GAD), islet cell antigen 512 (anti-ICA512) and islet cell cytoplasm (ICA). Anti-SOX13 were detected in 20 (18.3%) of 109 subjects with Type 1 diabetes, and more frequently in adults than in children (29% vs 10%). Anti-SOX13 usually occurred with anti-GAD but rarely with anti-ICA512. Seven sera positive for anti-SOX13 did not react with either GAD, ICA512 or islet cell cytoplasm indicating that anti-SOX13 represented a distinct population of antibodies. Reactivity to SOX13 represents a further autoantibody response in adults with Type 1 diabetes and may provide a useful disease marker in subjects in whom other autoantibody tests are negative.

A comparative study of antibody expressions in primary biliary cirrhosis and autoimmune cholangitis using phage display

Primary biliary cirrhosis (PBC) and autoimmune cholangitis (AIC) are serologic expressions of an autoimmune liver disease affecting biliary ductular cells. Previously we screened a phage-displayed random peptide library with polyclonal IgG from 2 Australian patients with PBC and derived peptides that identified a single conformational (discontinuous) epitope in the inner lipoyl domain of the E2 subunit of the pyruvate dehydrogenase complex (PDC-E2), the characteristic autoantigen in PBC. Here we have used phage display to investigate the reactivity of PBC sera from 2 ethnically and geographically distinct populations, Japanese and Australian, and the 2 serologic expressions, PBC and AIC. Random 7-mer and 12-mer peptide libraries were biopanned with IgG from 3 Japanese patients with PBC and 3 with AIC who did not have anti-PDC-E2. The phage clones (phagotopes) obtained were tested by capture enzyme-linked immunosorbent assay (ELISA) for reactivity with affinity-purified anti-PDC-E2, and compared with those obtained from Australian patients with PBC. Peptide sequences of the derived phagotopes and sequences derived by biopanning with irrelevant antisera were aligned to develop a guide tree based on physicochemical similarity. Both Australian and Japanese PBC-derived phagotopes were distributed in branches of the guide tree that contained the peptide sequences MH and FV previously identified as part of an immunodominant conformational epitope of PDC-E2, indicating that epitope selection was not influenced by the racial origin of the PBC sera. Biopanning with either PBC or AIC-derived IgG yielded phagotopes that reacted with anti-PDC-E2 by capture ELISA, further establishing that there is a similar autoimmune targeting in PBC and AIC.

Expression of protein tyrosine phosphatase-like molecule ICA512/IA-2 induces growth arrest in yeast cells and transfected mammalian cell lines

The ICA512/IA-2 molecule, a protein with similarity to receptor-type protein tyrosine phosphatases, was discovered during studies to identify autoantigens in Type 1 diabetes. The biological function of ICA512/IA-2 is unknown. We describe striking effects of ICA512/IA-2 on viability and growth of both yeast cells and cultured mammalian cells. In transformed yeast Saccharomyces cerevisiae cells, expression of ICA512/IA-2 induced growth retardation as judged by measurements of optical density and counts of colony-forming units. In contrast, expression of the intracellular domain (amino acids 600-979) of ICA512/IA-2 in yeast or mammalian cells had no such effects. In investigations on apoptosis, expression of ICA512/IA-2 in yeast cells caused loss of plasma membrane asymmetry, but not release of cytochrome c from mitochondria which did occur in a control system after expression of the pro-apoptotic molecule Bax. Possible interactions between ICA512/IA-2 and components of the cytoskeleton were not supported by studies on staining of fixed yeast cells with phalloidin-Texas Red. With transfected mammalian cell lines COS-7 and NIH3T3, expression of ICA512/IA-2 likewise induced growth arrest, with some of the morphological features of apoptosis. Thus obligatory expression of ICA512/IA-2 in eukaryotic cells causes disruption of cellular activities, with growth arrest in yeast and nuclear pycnosis/fragmentation in mammalian cells. A possible explanation is that growth inhibition reflects a part of the presently unknown function of ICA512/IA-2.

Antibodies to diabetes-associated autoantigens in Indian patients with Type 1 diabetes: prevalence of anti-ICA512/IA2 and anti-SOX13

We ascertained frequencies of autoantibodies to a suite of islet cell antigens including ICA512/IA2 and SOX13 in Asian Indians with Type 1 diabetes and in other forms of diabetes. Autoantibodies to ICA512/IA2 and SOX13 were tested by radioimmunoprecipitation assay, and results were amalgamated with previous data on antibodies to glutamic acid decarboxylase (GAD) and to islet cell cytoplasmic antigens (ICA). The frequency of anti-SOX13 was higher in Asian Indians than in Europids. Overall, the combined frequency for all autoantibodies to diabetes-associated antigens in Type 1 diabetes in Indians approached the frequency reported for Europids. There was an unexpectedly high frequency of autoantibody reactions to any one of the autoantigens tested (24%) in fibrocalculous pancreatic diabetes, however, individual autoantibody frequencies were relatively low. Our data indicate that, whatever the population studied, testing for multiple autoantigenic reactivities is more informative than more limited testing, and that there may be regional (presumably ethnically based) differences in levels of particular autoantibodies in cases of Type 1 diabetes.

Post-Panning Computer-Aided Analysis of Phagotope Collections Selected with Neurocysticercosis Patient Polyclonal Antibodies Separation of Disease-Relevant and Irrelevant Peptide Sequences

The homology of peptide sequences selected from a 7mer phage display library with antibodies elicited by the multicelled parasite Taenia solium in cerebrospinal fluid and serum of neurocysticercosis (NCC) patients and by antibodies of uninfected control patients with similar neurological complications of other ethiology (non-NCC) were analyzed using a PILEUP-Tudos sequence alignments program. The analysis generated dendrograms bearing two types of sequence clusters, those containing (1) only NCC patients-derived peptides and (2) both NCC- and control non-CC -- patient derivatives. By using ELISA, peptides that were selected by the antibodies were identified predominantly in the NCC-derived clusters. In repeated analysis in which sequences were added or removed, the first type of clusters maintained their structure, while the second type of clusters were split into many separate homology units dispersed throughout the guide tree. These results are interpreted as the ability of the analysis to segregate NCC-specific peptide sequences from other sequences. Altogether, this study demonstrates the high potential of the PILEUP-Tudos computer program to analyze phagotope collections recovered through biopanning with polyclonal antibodies elicited in patients by complex and as yet unknown multiple pathogenic antigens and to separate all phagotopes that are disease-relevant on the basis of the sequence homology.

Conformational epitopes on the diabetes autoantigen GAD65 identified by peptide phage display and molecular modeling

The major diabetes autoantigen, glutamic acid decarboxylase (GAD65), contains a region of sequence similarity, including six identical residues PEVKEK, to the P2C protein of coxsackie B virus, suggesting that cross-reactivity between coxsackie B virus and GAD65 can initiate autoimmune diabetes. We used the human islet cell mAbs MICA3 and MICA4 to identify the Ab epitopes of GAD65 by screening phage-displayed random peptide libraries. The identified peptide sequences could be mapped to a homology model of the pyridoxal phosphate (PLP) binding domain of GAD65. For MICA3, a surface loop containing the sequence PEVKEK and two adjacent exposed helixes were identified in the PLP binding domain as well as a region of the C terminus of GAD65 that has previously been identified as critical for MICA3 binding. To confirm that the loop containing the PEVKEK sequence contributes to the MICA3 epitope, this loop was deleted by mutagenesis. This reduced binding of MICA3 by 70%. Peptide sequences selected using MICA4 were rich in basic or hydroxyl-containing amino acids, and the surface of the GAD65 PLP-binding domain surrounding Lys358, which is known to be critical for MICA4 binding, was likewise rich in these amino acids. Also, the two phage most reactive with MICA4 encoded the motif VALxG, and the reverse of this sequence, LAV, was located in this same region. Thus, we have defined the MICA3 and MICA4 epitopes on GAD65 using the combination of phage display, molecular modeling, and mutagenesis and have provided compelling evidence for the involvement of the PEVKEK loop in the MICA3 epitope.