Predicting A/B compartments from histone modifications using deep learning

The three-dimensional organization of genomes plays a crucial role in essential biological processes. The segregation of chromatin into A and B compartments highlights regions of activity and inactivity, providing a window into the genomic activities specific to each cell type. Yet, the steep costs associated with acquiring Hi-C data, necessary for studying this compartmentalization across various cell types, pose a significant barrier in studying cell type specific genome organization. To address this, we present a prediction tool called compartment prediction using recurrent neural networks (CoRNN), which predicts compartmentalization of 3D genome using histone modification enrichment. CoRNN demonstrates robust cross-cell-type prediction of A/B compartments with an average AuROC of 90.9%. Cell-type-specific predictions align well with known functional elements, with H3K27ac and H3K36me3 identified as highly predictive histone marks. We further investigate our mispredictions and found that they are located in regions with ambiguous compartmental status. Furthermore, our model's generalizability is validated by predicting compartments in independent tissue samples, which underscores its broad applicability.

Novel spirocyclic dimer, SpiD3, targets chronic lymphocytic leukemia survival pathways with potent preclinical effects

Chronic lymphocytic leukemia (CLL) cell survival and growth is fueled by the induction of B-cell receptor (BCR) signaling within the tumor microenvironment (TME) driving activation of NF-κB signaling and the unfolded protein response (UPR). Malignant cells have higher basal levels of UPR posing a unique therapeutic window to combat CLL cell growth using pharmacological agents that induce accumulation of misfolded proteins. Frontline CLL therapeutics that directly target BCR signaling such as Bruton-tyrosine kinase (BTK) inhibitors (e.g., ibrutinib) have enhanced patient survival. However, resistance mechanisms wherein tumor cells bypass BTK inhibition through acquired BTK mutations, and/or activation of alternative survival mechanisms have rendered ibrutinib ineffective, imposing the need for novel therapeutics. We evaluated SpiD3, a novel spirocyclic dimer, in CLL cell lines, patient-derived CLL samples, ibrutinib-resistant CLL cells, and in the Eµ-TCL1 mouse model. Our integrated multi-omics and functional analyses revealed BCR signaling, NF-κB signaling, and endoplasmic reticulum stress among the top pathways modulated by SpiD3. This was accompanied by marked upregulation of the UPR and inhibition of global protein synthesis in CLL cell lines and patient-derived CLL cells. In ibrutinib-resistant CLL cells, SpiD3 retained its anti-leukemic effects, mirrored in reduced activation of key proliferative pathways (e.g., PRAS, ERK, MYC). Translationally, we observed reduced tumor burden in SpiD3-treated Eµ-TCL1 mice. Our findings reveal that SpiD3 exploits critical vulnerabilities in CLL cells including NF-κB signaling and the UPR, culminating in profound anti-tumor properties independent of TME stimuli.

Mechanistic drivers of chromatin organization into compartments

The human genome is not just a simple string of DNA, it is a complex and dynamic entity intricately folded within the cell's nucleus. This three-dimensional organization of chromatin, the combination of DNA and proteins in the nucleus, is crucial for many biological processes and has been prominently studied for its intricate relationship to gene expression. Indeed, the transcriptional machinery does not operate in isolation but interacts intimately with the folded chromatin structure. Techniques for chromatin conformation capture, including genome-wide sequencing approaches, have revealed key organizational features of chromatin, such as the formation of loops by CCCTC-binding factor (CTCF) and the division of loci into chromatin compartments. While much of the recent research and reviews have focused on CTCF loops, we discuss several new revelations that have emerged concerning chromatin compartments, with a particular focus on what is known about mechanistic drivers of compartmentalization. These insights challenge the traditional views of chromatin organization and reveal the complexity behind the formation and maintenance of chromatin compartments.

An environmentally relevant mixture of Perfluoroalkyl Substances (PFAS) impacts proliferation, steroid hormone synthesis, and gene transcription in primary human granulosa cells

Perfluoroalkyl substances (PFAS) are a group of synthetic chemicals that are resistant to biodegradation and are environmentally persistent. PFAS are found in many consumer products and are a major source of water and soil contamination. This study investigated the effects of an environmentally relevant PFAS mixture [perfluorooctanoic acid (PFOA), perfluorooctanesulfonic acid (PFOS), perfluorohexanesulfonic acid (PFHxS)] on the transcriptome and function of human granulosa cells (hGCs). Primary hGCs were harvested from follicular aspirates of healthy, reproductive-age women who were undergoing oocyte retrieval for in vitro fertilization. LC/MS-MS was performed to identify PFAS compounds in pure follicular fluid. Cells were cultured with vehicle control or a PFAS mixture (2 nM PFHxS, 7 nM PFOA, 10 nM PFOS) for 96h. Analyses of cell proliferation/apoptosis, steroidogenesis, and gene expression were measured via MTT assays/immunofluorescence, ELISA/western blotting, and RNA sequencing/bioinformatics, respectively. PFOA, PFOS, and PFHxS were detected in 100% of follicle fluid samples. Increased cell proliferation was observed in hGCs treated with the PFAS mixture with no impacts on cellular apoptosis. The PFAS mixture also altered steroid hormone synthesis, increasing both FSH-stimulated and basal progesterone secretion and concomitant upregulation of STAR protein. RNA sequencing revealed inherent differences in transcriptomic profiles in hGCs after PFAS exposure. This study demonstrates functional and transcriptomic changes in hGCs after exposure to a PFAS mixture, improving our knowledge about the impacts of PFAS exposures and female reproductive health. These findings suggest that PFAS compounds can disrupt normal granulosa cell function with possible long-term consequences on overall reproductive health.

The zinc-finger protein Z4 cooperates with condensin II to regulate somatic chromosome pairing and 3D chromatin organization

Chromosome pairing constitutes an important level of genome organization, yet the mechanisms that regulate pairing in somatic cells and the impact on 3D chromatin organization are still poorly understood. Here, we address these questions in Drosophila, an organism with robust somatic pairing. In Drosophila, pairing preferentially occurs at loci consisting of numerous architectural protein binding sites (APBSs), suggesting a role of architectural proteins (APs) in pairing regulation. Amongst these, the anti-pairing function of the condensin II subunit CAP-H2 is well established. However, the factors that regulate CAP-H2 localization and action at APBSs remain largely unknown. Here, we identify two factors that control CAP-H2 occupancy at APBSs and, therefore, regulate pairing. We show that Z4, interacts with CAP-H2 and is required for its localization at APBSs. We also show that hyperosmotic cellular stress induces fast and reversible unpairing in a Z4/CAP-H2 dependent manner. Moreover, by combining the opposite effects of Z4 depletion and osmostress, we show that pairing correlates with the strength of intrachromosomal 3D interactions, such as active (A) compartment interactions, intragenic gene-loops, and polycomb (Pc)-mediated chromatin loops. Altogether, our results reveal new players in CAP-H2-mediated pairing regulation and the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions.

HiCrayon reveals distinct layers of multi-state 3D chromatin organization

The co-visualization of chromatin conformation with 1D 'omics data is key to the multi-omics driven data analysis of 3D genome organization. Chromatin contact maps are often shown as 2D heatmaps and visually compared to 1D genomic data by simple juxtaposition. While common, this strategy is imprecise, placing the onus on the reader to align features with each other. To remedy this, we developed HiCrayon, an interactive tool that facilitates the integration of 3D chromatin organization maps and 1D datasets. This visualization method integrates data from genomic assays directly into the chromatin contact map by coloring interactions according to 1D signal. HiCrayon is implemented using R shiny and python to create a graphical user interface (GUI) application, available in both web or containerized format to promote accessibility. HiCrayon is implemented in R, and includes a graphical user interface (GUI), as well as a slimmed-down web-based version that lets users quickly produce publication-ready images. We demonstrate the utility of HiCrayon in visualizing the effectiveness of compartment calling and the relationship between ChIP-seq and various features of chromatin organization. We also demonstrate the improved visualization of other 3D genomic phenomena, such as differences between loops associated with CTCF/cohesin vs. those associated with H3K27ac. We then demonstrate HiCrayon's visualization of organizational changes that occur during differentiation and use HiCrayon to detect compartment patterns that cannot be assigned to either A or B compartments, revealing a distinct 3rd chromatin compartment. Overall, we demonstrate the utility of co-visualizing 2D chromatin conformation with 1D genomic signals within the same matrix to reveal fundamental aspects of genome organization. Local version: https://github.com/JRowleyLab/HiCrayon Web version: https://jrowleylab.com/HiCrayon.

Defining a cohort of anemia-activated cis elements reveals a mechanism promoting erythroid precursor function

Acute anemia elicits broad transcriptional changes in erythroid progenitors and precursors. We previously discovered a cis-regulatory transcriptional enhancer at the sterile alpha motif domain-14 enhancer locus (S14E), defined by a CANNTG-spacer-AGATAA composite motif and occupied by GATA1 and TAL1 transcription factors, is required for survival in severe anemia. However, S14E is only 1 of dozens of anemia-activated genes containing similar motifs. In a mouse model of acute anemia, we identified populations of expanding erythroid precursors, which increased expression of genes that contain S14E-like cis elements. We reveal that several S14E-like cis elements provide important transcriptional control of newly identified anemia-inducing genes, including the Ssx-2 interacting protein (Ssx2ip). Ssx2ip expression was determined to play an important role in erythroid progenitor/precursor cell activities, cell cycle regulation, and cell proliferation. Over a weeklong course of acute anemia recovery, we observed that erythroid gene activation mediated by S14E-like cis elements occurs during a phase coincident with low hematocrit and high progenitor activities, with distinct transcriptional programs activated at earlier and later time points. Our results define a genome-wide mechanism in which S14E-like enhancers control transcriptional responses during erythroid regeneration. These findings provide a framework to understand anemia-specific transcriptional mechanisms, ineffective erythropoiesis, anemia recovery, and phenotypic variability within human populations.

Regulation of CTCF loop formation during pancreatic cell differentiation

Transcription reprogramming during cell differentiation involves targeting enhancers to genes responsible for establishment of cell fates. To understand the contribution of CTCF-mediated chromatin organization to cell lineage commitment, we analyzed 3D chromatin architecture during the differentiation of human embryonic stem cells into pancreatic islet organoids. We find that CTCF loops are formed and disassembled at different stages of the differentiation process by either recruitment of CTCF to new anchor sites or use of pre-existing sites not previously involved in loop formation. Recruitment of CTCF to new sites in the genome involves demethylation of H3K9me3 to H3K9me2, demethylation of DNA, recruitment of pioneer factors, and positioning of nucleosomes flanking the new CTCF sites. Existing CTCF sites not involved in loop formation become functional loop anchors via the establishment of new cohesin loading sites containing NIPBL and YY1 at sites between the new anchors. In both cases, formation of new CTCF loops leads to strengthening of enhancer promoter interactions and increased transcription of genes adjacent to loop anchors. These results suggest an important role for CTCF and cohesin in controlling gene expression during cell differentiation.

Characterization of m6A modifiers and RNA modifications in uterine fibroids

Uterine leiomyoma or fibroids are the most common prevalent noncancerous tumors of the uterine muscle layer. Common symptoms associated with fibroids include pelvic pain, heavy menstrual bleeding, anemia, and pelvic pressure. These tumors are a leading cause of gynecological care but lack long-term therapy as the origin and development of fibroids are not well understood. Several next-generation sequencing technologies have been performed to identify the underlying genetic and epigenetic basis of fibroids. However, there remains a systemic gap in our understanding of molecular and biological process that define uterine fibroids. Recent epitranscriptomics studies have unraveled RNA modifications that are associated with all forms of RNA and are thought to influence both normal physiological functions and the progression of diseases. We quantified RNA expression profiles by analyzing publicly available RNA-seq data for 15 known epigenetic mediators to identify their expression profile in uterine fibroids compared to myometrium. To validate our findings, we performed RT-qPCR on a separate cohort of uterine fibroids targeting these modifiers confirming our RNA-seq data. We then examined protein profiles of key m6A modifiers in fibroids and their matched myometrium. In concordance with our RNA expression profiles, no significant differences were observed in these proteins in uterine fibroids compared to myometrium. To determine abundance of RNA modifications, mRNA and small RNA from fibroids and matched myometrium were analyzed by UHPLC MS/MS. In addition to the prevalent N6-methyladenosine (m6A), we identified 11 other known modifiers but did not identify any aberrant expression in fibroids. We then mined a previously published dataset and identified differential expression of m6A modifiers that were specific to fibroid genetic sub-type. Our analysis also identified m6A consensus motifs on genes previously identified to be dysregulated in uterine fibroids. Overall, using state-of-the-art mass spectrometry, RNA expression and protein profiles, we characterized and identified differentially expressed m6A modifiers in relation to driver mutations. Despite the use of several different approaches, we identified limited differential expression of RNA modifiers and associated modifications in uterine fibroids. However, considering the highly heterogenous genomic and cellular nature of fibroids, and the possible contribution of single molecule m6A modifications to fibroid pathology, there is a need for greater in-depth characterization of m6A marks and modifiers in a larger and varied patient cohort.

EHD1-dependent traffic of IGF-1 receptor to the cell surface is essential for Ewing sarcoma tumorigenesis and metastasis

Overexpression of the EPS15 Homology Domain containing 1 (EHD1) protein has been linked to tumorigenesis but whether its core function as a regulator of intracellular traffic of cell surface receptors plays a role in oncogenesis remains unknown. We establish that EHD1 is overexpressed in Ewing sarcoma (EWS), with high EHD1 mRNA expression specifying shorter patient survival. ShRNA-knockdown and CRISPR-knockout with mouse Ehd1 rescue established a requirement of EHD1 for tumorigenesis and metastasis. RTK antibody arrays identified IGF-1R as a target of EHD1 regulation in EWS. Mechanistically, we demonstrate a requirement of EHD1 for endocytic recycling and Golgi to plasma membrane traffic of IGF-1R to maintain its surface expression and downstream signaling. Conversely, EHD1 overexpression-dependent exaggerated oncogenic traits require IGF-1R expression and kinase activity. Our findings define the RTK traffic regulation as a proximal mechanism of EHD1 overexpression-dependent oncogenesis that impinges on IGF-1R in EWS, supporting the potential of IGF-1R and EHD1 co-targeting.

Chromatin alternates between A and B compartments at kilobase scale for subgenic organization

Nuclear compartments are prominent features of 3D chromatin organization, but sequencing depth limitations have impeded investigation at ultra fine-scale. CTCF loops are generally studied at a finer scale, but the impact of looping on proximal interactions remains enigmatic. Here, we critically examine nuclear compartments and CTCF loop-proximal interactions using a combination of in situ Hi-C at unparalleled depth, algorithm development, and biophysical modeling. Producing a large Hi-C map with 33 billion contacts in conjunction with an algorithm for performing principal component analysis on sparse, super massive matrices (POSSUMM), we resolve compartments to 500 bp. Our results demonstrate that essentially all active promoters and distal enhancers localize in the A compartment, even when flanking sequences do not. Furthermore, we find that the TSS and TTS of paused genes are often segregated into separate compartments. We then identify diffuse interactions that radiate from CTCF loop anchors, which correlate with strong enhancer-promoter interactions and proximal transcription. We also find that these diffuse interactions depend on CTCF's RNA binding domains. In this work, we demonstrate features of fine-scale chromatin organization consistent with a revised model in which compartments are more precise than commonly thought while CTCF loops are more protracted.

Considerations and caveats for analyzing chromatin compartments

Genomes are organized into nuclear compartments, separating active from inactive chromatin. Chromatin compartments are readily visible in a large number of species by experiments that map chromatin conformation genome-wide. When analyzing these maps, a common step is the identification of genomic intervals that interact within A (active) and B (inactive) compartments. It has also become increasingly common to identify and analyze subcompartments. We review different strategies to identify A/B and subcompartment intervals, including a discussion of various machine-learning approaches to predict these features. We then discuss the strengths and limitations of current strategies and examine how these aspects of analysis may have impacted our understanding of chromatin compartments.

Somatic chromosome pairing has a determinant impact on 3D chromatin organization

In the nucleus, chromatin is intricately structured into multiple layers of 3D organization important for genome activity. How distinct layers influence each other is not well understood. In particular, the contribution of chromosome pairing to 3D chromatin organization has been largely neglected. Here, we address this question in Drosophila, an organism that shows robust chromosome pairing in interphasic somatic cells. The extent of chromosome pairing depends on the balance between pairing and anti-pairing factors, with the anti-pairing activity of the CAP-H2 condensin II subunit being the best documented. Here, we identify the zinc-finger protein Z4 as a strong anti-pairer that interacts with and mediates the chromatin binding of CAP-H2. We also report that hyperosmotic cellular stress induces fast and reversible chromosome unpairing that depends on Z4/CAP-H2. And, most important, by combining Z4 depletion and osmostress, we show that chromosome pairing reinforces intrachromosomal 3D interactions. On the one hand, pairing facilitates RNAPII occupancy that correlates with enhanced intragenic gene-loop interactions. In addition, acting at a distance, pairing reinforces chromatin-loop interactions mediated by Polycomb (Pc). In contrast, chromosome pairing does not affect which genomic intervals segregate to active (A) and inactive (B) compartments, with only minimal effects on the strength of A-A compartmental interactions. Altogether, our results unveil the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions, unraveling the interwoven relationship between different layers of chromatin organization and the essential contribution of chromosome pairing.

Elevating SOX2 Downregulates MYC through a SOX2:MYC Signaling Axis and Induces a Slowly Cycling Proliferative State in Human Tumor Cells

Slowly cycling/infrequently proliferating tumor cells present a clinical challenge due to their ability to evade treatment. Previous studies established that high levels of SOX2 in both fetal and tumor cells restrict cell proliferation and induce a slowly cycling state. However, the mechanisms through which elevated SOX2 levels inhibit tumor cell proliferation have not been identified. To identify common mechanisms through which SOX2 elevation restricts tumor cell proliferation, we initially performed RNA-seq using two diverse tumor cell types. SOX2 elevation in both cell types downregulated MYC target genes. Consistent with these findings, elevating SOX2 in five cell lines representing three different human cancer types decreased MYC expression. Importantly, the expression of a dominant-negative MYC variant, omomyc, recapitulated many of the effects of SOX2 on proliferation, cell cycle, gene expression, and biosynthetic activity. We also demonstrated that rescuing MYC activity in the context of elevated SOX2 induces cell death, indicating that the downregulation of MYC is a critical mechanistic step necessary to maintain survival in the slowly cycling state induced by elevated SOX2. Altogether, our findings uncover a novel SOX2:MYC signaling axis and provide important insights into the molecular mechanisms through which SOX2 elevation induces a slowly cycling proliferative state.

Implications of Dosage Deficiencies in CTCF and Cohesin on Genome Organization, Gene Expression, and Human Neurodevelopment

Properly organizing DNA within the nucleus is critical to ensure normal downstream nuclear functions. CTCF and cohesin act as major architectural proteins, working in concert to generate thousands of high-intensity chromatin loops. Due to their central role in loop formation, a massive research effort has been dedicated to investigating the mechanism by which CTCF and cohesin create these loops. Recent results lead to questioning the direct impact of CTCF loops on gene expression. Additionally, results of controlled depletion experiments in cell lines has indicated that genome architecture may be somewhat resistant to incomplete deficiencies in CTCF or cohesin. However, heterozygous human genetic deficiencies in CTCF and cohesin have illustrated the importance of their dosage in genome architecture, cellular processes, animal behavior, and disease phenotypes. Thus, the importance of considering CTCF or cohesin levels is especially made clear by these heterozygous germline variants that characterize genetic syndromes, which are increasingly recognized in clinical practice. Defined primarily by developmental delay and intellectual disability, the phenotypes of CTCF and cohesin deficiency illustrate the importance of architectural proteins particularly in neurodevelopment. We discuss the distinct roles of CTCF and cohesin in forming chromatin loops, highlight the major role that dosage of each protein plays in the amplitude of observed effects on gene expression, and contrast these results to heterozygous mutation phenotypes in murine models and clinical patients. Insights highlighted by this comparison have implications for future research into these newly emerging genetic syndromes.

Ecdysoneless Protein Regulates Viral and Cellular mRNA Splicing to Promote Cervical Oncogenesis

High-risk human papillomaviruses (HPV), exemplified by HPV16/18, are causally linked to human cancers of the anogenital tract, skin, and upper aerodigestive tract. Previously, we identified Ecdysoneless (ECD) protein, the human homolog of the Drosophila ecdysoneless gene, as a novel HPV16 E6-interacting protein. Here, we show that ECD, through its C-terminal region, selectively binds to high-risk but not to low-risk HPV E6 proteins. We demonstrate that ECD is overexpressed in cervical and head and neck squamous cell carcinoma (HNSCC) cell lines as well as in tumor tissues. Using The Cancer Genome Atlas dataset, we show that ECD mRNA overexpression predicts shorter survival in patients with cervical and HNSCC. We demonstrate that ECD knockdown in cervical cancer cell lines led to impaired oncogenic behavior, and ECD co-overexpression with E7 immortalized primary human keratinocytes. RNA-sequencing analyses of SiHa cells upon ECD knockdown showed to aberrations in E6/E7 RNA splicing, as well as RNA splicing of several HPV oncogenesis-linked cellular genes, including splicing of components of mRNA splicing machinery itself. Taken together, our results support a novel role of ECD in viral and cellular mRNA splicing to support HPV-driven oncogenesis. IMPLICATIONS: This study links ECD overexpression to poor prognosis and shorter survival in HNSCC and cervical cancers and identifies a critical role of ECD in cervical oncogenesis through regulation of viral and cellular mRNA splicing.

Evolutionary History and Activity of RNase H1-Like Proteins in Arabidopsis thaliana

RNase H1 is an endonuclease specific toward the RNA strand of RNA:DNA hybrids. Members of this protein family are present in most living organisms and are essential for removing RNA that base pairs with DNA. It prevents detrimental effects of RNA:DNA hybrids and is involved in several biological processes. Arabidopsis thaliana has been previously shown to contain three genes encoding RNase H1 proteins that localize to three distinct cellular compartments. We show that these genes originate from two gene duplication events. One occurred in the common ancestor of dicots and produced nuclear and organellar RNase H1 paralogs. Second duplication occurred in the common ancestor of Brassicaceae and produced mitochondrial- and plastid-localized proteins. These proteins have the canonical RNase H1 activity, which requires at least four ribonucleotides for endonucleolytic digestion. Analysis of mutants in the RNase H1 genes revealed that the nuclear RNH1A and mitochondrial RNH1B are dispensable for development under normal growth conditions. However, the presence of at least one organellar RNase H1 (RNH1B or RNH1C) is required for embryonic development. The plastid-localized RNH1C affects plastid DNA copy number and sensitivity to replicative stress. Our results present the evolutionary history of RNH1 proteins in A. thaliana, demonstrate their canonical RNase H1 activity and indicate their role in early embryonic development.

Author Correction: Widespread long-range cis-regulatory elements in the maize genome

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Analysis of Hi-C data using SIP effectively identifies loops in organisms from C. elegans to mammals

Chromatin loops are a major component of 3D nuclear organization, visually apparent as intense point-to-point interactions in Hi-C maps. Identification of these loops is a critical part of most Hi-C analyses. However, current methods often miss visually evident CTCF loops in Hi-C data sets from mammals, and they completely fail to identify high intensity loops in other organisms. We present SIP, Significant Interaction Peak caller, and SIPMeta, which are platform independent programs to identify and characterize these loops in a time- and memory-efficient manner. We show that SIP is resistant to noise and sequencing depth, and can be used to detect loops that were previously missed in human cells as well as loops in other organisms. SIPMeta corrects for a common visualization artifact by accounting for Manhattan distance to create average plots of Hi-C and HiChIP data. We then demonstrate that the use of SIP and SIPMeta can lead to biological insights by characterizing the contribution of several transcription factors to CTCF loop stability in human cells. We also annotate loops associated with the SMC component of the dosage compensation complex (DCC) in Caenorhabditis elegans and demonstrate that loop anchors represent bidirectional blocks for symmetrical loop extrusion. This is in contrast to the asymmetrical extrusion until unidirectional blockage by CTCF that is presumed to occur in mammals. Using HiChIP and multiway ligation events, we then show that DCC loops form a network of strong interactions that may contribute to X Chromosome-wide condensation in C. elegans hermaphrodites.

Widespread long-range cis-regulatory elements in the maize genome

Genetic mapping studies on crops suggest that agronomic traits can be controlled by gene-distal intergenic loci. Despite the biological importance and the potential agronomic utility of these loci, they remain virtually uncharacterized in all crop species to date. Here, we provide genetic, epigenomic and functional molecular evidence to support the widespread existence of gene-distal (hereafter, distal) loci that act as long-range transcriptional cis-regulatory elements (CREs) in the maize genome. Such loci are enriched for euchromatic features that suggest their regulatory functions. Chromatin loops link together putative CREs with genes and recapitulate genetic interactions. Putative CREs also display elevated transcriptional enhancer activities, as measured by self-transcribing active regulatory region sequencing. These results provide functional support for the widespread existence of CREs that act over large genomic distances to control gene expression.

Immediate and deferred epigenomic signatures of in vivo neuronal activation in mouse hippocampus

Activity-driven transcription plays an important role in many brain processes, including those underlying memory and epilepsy. Here we combine genetic tagging of nuclei and ribosomes with RNA sequencing, chromatin immunoprecipitation with sequencing, assay for transposase-accessible chromatin using sequencing and Hi-C to investigate transcriptional and chromatin changes occurring in mouse hippocampal excitatory neurons at different time points after synchronous activation during seizure and sparse activation by novel context exploration. The transcriptional burst is associated with an increase in chromatin accessibility of activity-regulated genes and enhancers, de novo binding of activity-regulated transcription factors, augmented promoter-enhancer interactions and the formation of gene loops that bring together the transcription start site and transcription termination site of induced genes and may sustain the fast reloading of RNA polymerase complexes. Some chromatin occupancy changes and interactions, particularly those driven by AP1, remain long after neuronal activation and could underlie the changes in neuronal responsiveness and circuit connectivity observed in these neuroplasticity paradigms, perhaps thereby contributing to metaplasticity in the adult brain.

Maintenance of CTCF- and Transcription Factor-Mediated Interactions from the Gametes to the Early Mouse Embryo

The epigenetic information present in mammalian gametes and whether it is transmitted to the progeny are relatively unknown. We find that many promoters in mouse sperm are occupied by RNA polymerase II (Pol II) and Mediator. The same promoters are accessible in GV and MII oocytes and preimplantation embryos. Sperm distal ATAC-seq sites containing motifs for various transcription factors are conserved in monkeys and humans. ChIP-seq analyses confirm that Foxa1, ERα, and AR occupy distal enhancers in sperm. Accessible sperm enhancers containing H3.3 and H2A.Z are also accessible in oocytes and preimplantation embryos. Furthermore, their interactions with promoters in the gametes persist during early development. Sperm- or oocyte-specific interactions mediated by CTCF and cohesin are only present in the paternal or maternal chromosomes, respectively, in the zygote and 2-cell stages. These interactions converge in both chromosomes by the 8-cell stage. Thus, mammalian gametes contain complex patterns of 3D interactions that can be transmitted to the zygote after fertilization.

Condensin II Counteracts Cohesin and RNA Polymerase II in the Establishment of 3D Chromatin Organization

Interaction domains in Drosophila chromosomes form by segregation of active and inactive chromatin in the absence of CTCF loops, but the role of transcription versus other architectural proteins in chromatin organization is unclear. Here, we find that positioning of RNAPII via transcription elongation is essential in the formation of gene loops, which in turn interact to form compartmental domains. Inhibition of transcription elongation or depletion of cohesin decreases gene looping and formation of active compartmental domains. In contrast, depletion of condensin II, which also localizes to active chromatin, causes increased gene looping, formation of compartmental domains, and stronger intra-chromosomal compartmental interactions. Condensin II has a similar role in maintaining inter-chromosomal interactions responsible for pairing between homologous chromosomes, whereas inhibition of transcription elongation or cohesin depletion has little effect on homolog pairing. The results suggest distinct roles for cohesin and condensin II in the establishment of 3D nuclear organization in Drosophila.

Organizational principles of 3D genome architecture

Studies of 3D chromatin organization have suggested that chromosomes are hierarchically organized into large compartments composed of smaller domains called topologically associating domains (TADs). Recent evidence suggests that compartments are smaller than previously thought and that the transcriptional or chromatin state is responsible for interactions leading to the formation of small compartmental domains in all organisms. In vertebrates, CTCF forms loop domains, probably via an extrusion process involving cohesin. CTCF loops cooperate with compartmental domains to establish the 3D organization of the genome. The continuous extrusion of the chromatin fibre by cohesin may also be responsible for the establishment of enhancer-promoter interactions and stochastic aspects of the transcription process. These observations suggest that the 3D organization of the genome is an emergent property of chromatin and its components, and thus may not be only a determinant but also a consequence of its function.

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.

Characterization of epitopes for virus-neutralizing monoclonal antibodies to Ross River virus E2 using phage-displayed random peptide libraries

Ross River virus (RRV) is the predominant cause of epidemic polyarthritis in Australia, yet the antigenic determinants are not well defined. We aimed to characterize epitope(s) on RRV-E2 for a panel of monoclonal antibodies (MAbs) that recognize overlapping conformational epitopes on the E2 envelope protein of RRV and that neutralize virus infection of cells in vitro. Phage-displayed random peptide libraries were probed with the MAbs T1E7, NB3C4, and T10C9 using solution-phase and solid-phase biopanning methods. The peptides VSIFPPA and KTAISPT were selected 15 and 6 times, respectively, by all three of the MAbs using solution-phase biopanning. The peptide LRLPPAP was selected 8 times by NB3C4 using solid-phase biopanning; this peptide shares a trio of amino acids with the peptide VSIFPPA. Phage that expressed the peptides VSIFPPA and LRLPPAP were reactive with T1E7 and/or NB3C4, and phage that expressed the peptides VSIFPPA, LRLPPAP, and KTAISPT partially inhibited the reactivity of T1E7 with RRV. The selected peptides resemble regions of RRV-E2 adjacent to sites mutated in neutralization escape variants of RRV derived by culture in the presence of these MAbs (E2 210-219 and 238-245) and an additional region of E2 172-182. Together these sites represent a conformational epitope of E2 that is informative of cellular contact sites on RRV.