Enhancer-promoter communication and transcriptional regulation of Igh

Design and deep learning of synthetic B-cell-specific promoters

Synthetic biology and deep learning synergistically revolutionize our ability for decoding and recoding DNA regulatory grammar. The B-cell-specific transcriptional regulation is intricate, and unlock the potential of B-cell-specific promoters as synthetic elements is important for B-cell engineering. Here, we designed and pooled synthesized 23 640 B-cell-specific promoters that exhibit larger sequence space, B-cell-specific expression, and enable diverse transcriptional patterns in B-cells. By MPRA (Massively parallel reporter assays), we deciphered the sequence features that regulate promoter transcriptional, including motifs and motif syntax (their combination and distance). Finally, we built and trained a deep learning model capable of predicting the transcriptional strength of the immunoglobulin V gene promoter directly from sequence. Prediction of thousands of promoter variants identified in the global human population shows that polymorphisms in promoters influence the transcription of immunoglobulin V genes, which may contribute to individual differences in adaptive humoral immune responses. Our work helps to decipher the transcription mechanism in immunoglobulin genes and offers thousands of non-similar promoters for B-cell engineering.

Temporal analyses reveal a pivotal role for sense and antisense enhancer RNAs in coordinate immunoglobulin lambda locus activation

Transcription enhancers are essential activators of V(D)J recombination that orchestrate non-coding transcription through complementary, unrearranged gene segments. How transcription is coordinately increased at spatially distinct promoters, however, remains poorly understood. Using the murine immunoglobulin lambda (Igλ) locus as model, we find that three enhancer-like elements in the 3' Igλ domain, Eλ3-1, HSCλ1 and HSE-1, show strikingly similar transcription factor binding dynamics and close spatial proximity, suggesting that they form an active enhancer hub. Temporal analyses show coordinate recruitment of complementary V and J gene segments to this hub, with comparable transcription factor binding dynamics to that at enhancers. We find further that E2A, p300, Mediator and Integrator bind to enhancers as early events, whereas YY1 recruitment and eRNA synthesis occur later, corresponding to transcription activation. Remarkably, the interplay between sense and antisense enhancer RNA is central to both active enhancer hub formation and coordinate Igλ transcription: Antisense Eλ3-1 eRNA represses Igλ activation whereas temporal analyses demonstrate that accumulating levels of sense eRNA boost YY1 recruitment to stabilise enhancer hub/promoter interactions and lead to coordinate transcription activation. These studies therefore demonstrate for the first time a critical role for threshold levels of sense versus antisense eRNA in locus activation.

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Genome-Wide Analysis of Proline-Rich Extensin-Like Receptor Kinases (PERKs) Gene Family Reveals Their Roles in Plant Development and Stress Conditions in Oryza sativa L

Proline-rich extensin-like receptor kinases (PERKs) play a crucial role in a wide range of biological processes in plants. In model plants like Arabidopsis, the PERK gene family has been well investigated. Conversely, no information available on the PERK gene family and their biological functions largely remained unknown in rice. This study analyzed the basic physicochemical properties, phylogeny, gene structure, cis-acting elements, Gene ontology (GO) annotation and protein-protein interaction of OsPERK gene family members using various bioinformatics tools based on the whole-genome data of O. sativa. Thus, in this work, 8 PERK genes in rice were identified, and their roles in plant development, growth, and response to various stresses were studied. A phylogenetic study revealed that OsPERKs are grouped into seven classes. Chromosomal mapping also displayed that 8 PERK genes were unevenly distributed on 12 chromosomes. Further, the prediction of subcellular localization indicated that OsPERKs were mainly located at the endomembrane system. Gene structure analysis of OsPERKs has shown a distinctive evolutionary path. In addition, synteny analysis exhibited the 40 orthologous gene pairs in Arabidopsis thaliana, Triticum aestivum, Hordeum vulgare and Medicago truncatula. Furthermore, Ka to Ks proportion shows that most OsPERK genes experienced resilient purifying selection during evolutionary processes. The OsPERK promoters contained several cis-acting regulatory, which are crucial for plant development processes, phytohormone signaling, stress, and defense response. Moreover, the expression pattern of OsPERK family members showed differential expression patterns in different tissues and various stress conditions. Taken together, these results provide clear messages for a better understanding the roles of OsPERK genes in various development stages, tissues, and multifactorial stress as well as enriched the related research of OsPERK family members in rice.

Delta.EPI: a probabilistic voting-based enhancer-promoter interaction prediction platform

Enhancer promoter interaction (EPI) involves most of gene transcriptional regulation in the high eukaryotes. Predicting the EPIs from given genomic loci or DNA sequences is not a trivial task. The benchmarking work so far for EPI predictors is more or less empirical and lacks quantitative model-based comparisons, posing challenges for molecular biologists to obtain reliable EPI predictions. Here, we present an EPI prediction platform, Delta.EPI. Based on a statistic model of the data integration, Delta.EPI is capable of comprehensively assessing the predictions from four state-of-the-art EPI predictors. Equipped with a user-friendly interface and visualization platform, Delta.EPI presents the sorted results with the confidence of EPI relevance, which may guide the molecular biologists who lack the pre-knowledge of the algorithms of EPI prediction. Last, we showcase the utility of Delta.EPI with a case study. Delta.EPI provides a powerful tool to fuel the gene regulation and 3D genome studies by ease-to-access EPI predictions. Delta.EPI can be freely access at https://ngdc.cncb.ac.cn/deltaEPI/.

Opposing Roles of FACT for Euchromatin and Heterochromatin in Yeast

DNA is stored in the nucleus of a cell in a folded state; however, only the necessary genetic information is extracted from the required group of genes. The key to extracting genetic information is chromatin ambivalence. Depending on the chromosomal region, chromatin is characterized into low-density "euchromatin" and high-density "heterochromatin", with various factors being involved in its regulation. Here, we focus on chromatin regulation and gene expression by the yeast FACT complex, which functions in both euchromatin and heterochromatin. FACT is known as a histone H2A/H2B chaperone and was initially reported as an elongation factor associated with RNA polymerase II. In budding yeast, FACT activates promoter chromatin by interacting with the transcriptional activators SBF/MBF via the regulation of G1/S cell cycle genes. In fission yeast, FACT plays an important role in the formation of higher-order chromatin structures and transcriptional repression by binding to Swi6, an HP1 family protein, at heterochromatin. This FACT property, which refers to the alternate chromatin-regulation depending on the binding partner, is an interesting phenomenon. Further analysis of nucleosome regulation within heterochromatin is expected in future studies.

Evolutive emergence and divergence of an Ig regulatory node: An environmental sensor getting cues from the aryl hydrocarbon receptor?

One gene, the immunoglobulin heavy chain (IgH) gene, is responsible for the expression of all the different antibody isotypes. Transcriptional regulation of the IgH gene is complex and involves several regulatory elements including a large element at the 3' end of the IgH gene locus (3'RR). Animal models have demonstrated an essential role of the 3'RR in the ability of B cells to express high affinity antibodies and to express different antibody classes. Additionally, environmental chemicals such as aryl hydrocarbon receptor (AhR) ligands modulate mouse 3'RR activity that mirrors the effects of these chemicals on antibody production and immunocompetence in mouse models. Although first discovered as a mediator of the toxicity induced by the high affinity ligand 2,3,7,8-tetracholordibenzo-p-dioxin (dioxin), understanding of the AhR has expanded to a physiological role in preserving homeostasis and maintaining immunocompetence. We posit that the AhR also plays a role in human antibody production and that the 3'RR is not only an IgH regulatory node but also an environmental sensor receiving signals through intrinsic and extrinsic pathways, including the AhR. This review will 1) highlight the emerging role of the AhR as a key transducer between environmental signals and altered immune function; 2) examine the current state of knowledge regarding IgH gene regulation and the role of the AhR in modulation of Ig production; 3) describe the evolution of the IgH gene that resulted in species and population differences; and 4) explore the evidence supporting the environmental sensing capacity of the 3'RR and the AhR as a transducer of these cues. This review will also underscore the need for studies focused on human models due to the premise that understanding genetic differences in the human population and the signaling pathways that converge at the 3'RR will provide valuable insight into individual sensitivities to environmental factors and antibody-mediated disease conditions, including emerging infections such as SARS-CoV-2.

Low input capture Hi-C (liCHi-C) identifies promoter-enhancer interactions at high-resolution

Long-range interactions between regulatory elements and promoters are key in gene transcriptional control; however, their study requires large amounts of starting material, which is not compatible with clinical scenarios nor the study of rare cell populations. Here we introduce low input capture Hi-C (liCHi-C) as a cost-effective, flexible method to map and robustly compare promoter interactomes at high resolution. As proof of its broad applicability, we implement liCHi-C to study normal and malignant human hematopoietic hierarchy in clinical samples. We demonstrate that the dynamic promoter architecture identifies developmental trajectories and orchestrates transcriptional transitions during cell-state commitment. Moreover, liCHi-C enables the identification of disease-relevant cell types, genes and pathways potentially deregulated by non-coding alterations at distal regulatory elements. Finally, we show that liCHi-C can be harnessed to uncover genome-wide structural variants, resolve their breakpoints and infer their pathogenic effects. Collectively, our optimized liCHi-C method expands the study of 3D chromatin organization to unique, low-abundance cell populations, and offers an opportunity to uncover factors and regulatory networks involved in disease pathogenesis.

Genome-Wide Identification and Expression Profiling of Aconitase Gene Family Members Reveals Their Roles in Plant Development and Adaptation to Diverse Stress in Triticum aestivum L

Global warming is a serious threat to food security and severely affects plant growth, developmental processes, and, eventually, crop productivity. Respiratory metabolism plays a critical role in the adaptation of diverse stress in plants. Aconitase (ACO) is the main enzyme, which catalyzes the revocable isomerization of citrate to isocitrate in the Krebs cycle. The function of ACO gene family members has been extensively studied in model plants, for instance Arabidopsis. However, their role in plant developmental processes and various stress conditions largely remained unknown in other plant species. Thus, we identified 15 ACO genes in wheat to elucidate their function in plant developmental processes and different stress environments. The phylogenetic tree revealed that TaACO genes were classified into six groups. Further, gene structure analysis of TaACOs has shown a distinctive evolutionary path. Synteny analysis showed the 84 orthologous gene pairs in Brachypodium distachyon, Aegilops tauschii, Triticum dicoccoides, Oryza sativa, and Arabidopsis thaliana. Furthermore, Ka/Ks ratio revealed that most TaACO genes experienced strong purifying selection during evolution. Numerous cis-acting regulatory elements were detected in the TaACO promoters, which play a crucial role in plant development processes, phytohormone signaling, and are related to defense and stress. To understand the function of TaACO genes, the expression profiling of TaACO genes were investigated in different tissues, developmental stages, and stress conditions. The transcript per million values of TaACOs genes were retrieved from the Wheat Expression Browser Database. We noticed the differential expression of the TaACO genes in different tissues and various stress conditions. Moreover, gene ontology analysis has shown enrichment in the tricarboxylic acid metabolic process (GO:0072350), citrate metabolic process (GO:0006101), isocitrate metabolic process GO:0006102, carbohydrate metabolic (GO:0005975), and glyoxylate metabolic process (GO:0046487). Therefore, this study provided valuable insight into the ACO gene family in wheat and contributed to the further functional characterization of TaACO during different plant development processes and various stress conditions.

Molecular characterization of the genome-wide BOR transporter family and their responses to boron conditions in common wheat (Triticum aestivum L.)

Boron (B) deficiency is an agricultural problem that causes significant yield losses in many countries. B transporters (BORs) are responsible for B uptake and distribution and play important roles in yield formation. A comprehensive analysis of the BOR family members in common wheat is still lacking. In the present study, to clarify the molecular characterization and response to B status, genome-wide TaBOR genes and expression patterns were investigated. Fourteen TaBOR genes were identified in common wheat by a homology search. The corresponding phylogenetic tree indicated that 14 TaBOR genes were separately classified into subfamilies of TaBOR1, TaBOR3, and TaBOR4. All TaBOR genes had 12-14 extrons and 11-13 introns. Most TaBOR proteins contained 10 conserved motifs, and motifs 1, 2, 3, 4, and 6 constituted the conserved bicarbonate (HCO₃ ^-) domain. Fourteen TaBOR genes were mapped on 13 chromosomes mainly distributed in the first, third, fifth, and seventh homologous groups. The promoters of TaBOR genes consisted of phytohormones, light responses, and stress-related cis-elements. GO analysis indicated that TaBOR genes were enriched in terms of transmembrane transport and ion homeostasis. TaBOR genes showed diverse expression profiles in different tissues. The members of the TaBOR1 subfamily showed high expression in grains, leaves, roots, stems, and spikes, but members of the TaBOR4 subfamily were highly expressed only in spikes and grains. RT-qPCR indicated that TaBOR1-5A, TaBOR1-5B, and TaBOR1-5D were induced by low B concentrations and had much higher expression in roots than in shoots. TaBOR3-3A, TaBOR3-3B, TaBOR3-3D, TaBOR4-1A, TaBOR4-1B, TaBOR4-1D, and TaBOR3-4B were induced by low and high B concentrations and had high expression in roots and shoots. TaBOR3-4D and TaBOR3-7B were upregulated by low and high B concentrations, respectively, but had expression only in roots. Our results provide basic information on the TaBOR family, which is beneficial for elucidating the functions of TaBOR genes to overcome the problem of B deficiency.

Genome-Wide Analysis and Characterization of the Proline-Rich Extensin-like Receptor Kinases (PERKs) Gene Family Reveals Their Role in Different Developmental Stages and Stress Conditions in Wheat (Triticum aestivum L.)

Proline-rich extensin-like receptor kinases (PERKs) are a class of receptor kinases implicated in multiple cellular processes in plants. However, there is a lack of information on the PERK gene family in wheat. Therefore, we identified 37 PERK genes in wheat to understand their role in various developmental processes and stress conditions. Phylogenetic analysis of PERK genes from Arabidopsis thaliana, Oryza sativa, Glycine max, and T. aestivum grouped them into eight well-defined classes. Furthermore, synteny analysis revealed 275 orthologous gene pairs in B. distachyon, Ae. tauschii, T. dicoccoides, O. sativa and A. thaliana. Ka/Ks values showed that most TaPERK genes, except TaPERK1, TaPERK2, TaPERK17, and TaPERK26, underwent strong purifying selection during evolutionary processes. Several cis-acting regulatory elements, essential for plant growth and development and the response to light, phytohormones, and diverse biotic and abiotic stresses, were predicted in the promoter regions of TaPERK genes. In addition, the expression profile of the TaPERK gene family revealed differential expression of TaPERK genes in various tissues and developmental stages. Furthermore, TaPERK gene expression was induced by various biotic and abiotic stresses. The RT-qPCR analysis also revealed similar results with slight variation. Therefore, this study's outcome provides valuable information for elucidating the precise functions of TaPERK in developmental processes and diverse stress conditions in wheat.

May the Odds Be Ever in Your Favor: Non-deterministic Mechanisms Diversifying Cell Surface Molecule Expression

How does the information in the genome program the functions of the wide variety of cells in the body? While the development of biological organisms appears to follow an explicit set of genomic instructions to generate the same outcome each time, many biological mechanisms harness molecular noise to produce variable outcomes. Non-deterministic variation is frequently observed in the diversification of cell surface molecules that give cells their functional properties, and is observed across eukaryotic clades, from single-celled protozoans to mammals. This is particularly evident in immune systems, where random recombination produces millions of antibodies from only a few genes; in nervous systems, where stochastic mechanisms vary the sensory receptors and synaptic matching molecules produced by different neurons; and in microbial antigenic variation. These systems employ overlapping molecular strategies including allelic exclusion, gene silencing by constitutive heterochromatin, targeted double-strand breaks, and competition for limiting enhancers. Here, we describe and compare five stochastic molecular mechanisms that produce variety in pathogen coat proteins and in the cell surface receptors of animal immune and neuronal cells, with an emphasis on the utility of non-deterministic variation.

Genome-Wide Identification and Characterization of the Brassinazole-resistant (BZR) Gene Family and Its Expression in the Various Developmental Stage and Stress Conditions in Wheat (Triticum aestivum L.)

Brassinosteroids (BRs) play crucial roles in various biological processes, including plant developmental processes and response to diverse biotic and abiotic stresses. However, no information is currently available about this gene family in wheat (Triticum aestivum L.). In the present investigation, we identified the BZR gene family in wheat to understand the evolution and their role in diverse developmental processes and under different stress conditions. In this study, we performed the genome-wide analysis of the BZR gene family in the bread wheat and identified 20 TaBZR genes through a homology search and further characterized them to understand their structure, function, and distribution across various tissues. Phylogenetic analyses lead to the classification of TaBZR genes into five different groups or subfamilies, providing evidence of evolutionary relationship with Arabidopsis thaliana, Zea mays, Glycine max, and Oryza sativa. A gene exon/intron structure analysis showed a distinct evolutionary path and predicted the possible gene duplication events. Further, the physical and biochemical properties, conserved motifs, chromosomal, subcellular localization, and cis-acting regulatory elements were also examined using various computational approaches. In addition, an analysis of public RNA-seq data also shows that TaBZR genes may be involved in diverse developmental processes and stress tolerance mechanisms. Moreover, qRT-PCR results also showed similar expression with slight variation. Collectively, these results suggest that TaBZR genes might play an important role in plant developmental processes and various stress conditions. Therefore, this work provides valuable information for further elucidate the precise role of BZR family members in wheat.

Genome-Wide Identification and Characterization of PIN-FORMED (PIN) Gene Family Reveals Role in Developmental and Various Stress Conditions in Triticum aestivum L

PIN-FORMED (PIN) genes play a crucial role in regulating polar auxin distribution in diverse developmental processes, including tropic responses, embryogenesis, tissue differentiation, and organogenesis. However, the role of PIN-mediated auxin transport in various plant species is poorly understood. Currently, no information is available about this gene family in wheat (Triticum aestivum L.). In the present investigation, we identified the PIN gene family in wheat to understand the evolution of PIN-mediated auxin transport and its role in various developmental processes and under different biotic and abiotic stress conditions. In this study, we performed genome-wide analysis of the PIN gene family in common wheat and identified 44 TaPIN genes through a homology search, further characterizing them to understand their structure, function, and distribution across various tissues. Phylogenetic analyses led to the classification of TaPIN genes into seven different groups, providing evidence of an evolutionary relationship with Arabidopsis thaliana and Oryza sativa. A gene exon/intron structure analysis showed a distinct evolutionary path and predicted the possible gene duplication events. Further, the physical and biochemical properties, conserved motifs, chromosomal, subcellular localization, transmembrane domains, and three-dimensional (3D) structure were also examined using various computational approaches. Cis-elements analysis of TaPIN genes showed that TaPIN promoters consist of phytohormone, plant growth and development, and stress-related cis-elements. In addition, expression profile analysis also revealed that the expression patterns of the TaPIN genes were different in different tissues and developmental stages. Several members of the TaPIN family were induced during biotic and abiotic stress. Moreover, the expression patterns of TaPIN genes were verified by qRT-PCR. The qRT-PCR results also show a similar expression with slight variation. Therefore, the outcome of this study provides basic genomic information on the expression of the TaPIN gene family and will pave the way for dissecting the precise role of TaPINs in plant developmental processes and different stress conditions.

Immunoglobulin Expression in Cancer Cells and Its Critical Roles in Tumorigenesis

Traditionally, immunoglobulin (Ig) was believed to be produced by only B-lineage cells. However, increasing evidence has revealed a high level of Ig expression in cancer cells, and this Ig is named cancer-derived Ig. Further studies have shown that cancer-derived Ig shares identical basic structures with B cell-derived Ig but exhibits several distinct characteristics, including restricted variable region sequences and aberrant glycosylation. In contrast to B cell-derived Ig, which functions as an antibody in the humoral immune response, cancer-derived Ig exerts profound protumorigenic effects via multiple mechanisms, including promoting the malignant behaviors of cancer cells, mediating tumor immune escape, inducing inflammation, and activating the aggregation of platelets. Importantly, cancer-derived Ig shows promising potential for application as a diagnostic and therapeutic target in cancer patients. In this review, we summarize progress in the research area of cancer-derived Ig and discuss the perspectives of applying this novel target for the management of cancer patients.

SLy2-overexpression impairs B-cell development in the bone marrow and the IgG response towards pneumococcal conjugate-vaccine

Background

Infections with Streptococcus pneumoniae can cause severe diseases in humans including pneumonia. Although guidelines for vaccination have been established, S. pneumoniae is still responsible for a serious burden of disease around the globe. Currently, two pneumococcal immunizations are available, namely the pure polysaccharide vaccine Pneumovax23 (P23) and the conjugate‐vaccine Prevenar13 (PCV13). We recently reported impaired thymus‐independent antibody responses towards P23 in mice overexpressing the immunoinhibitory adapter SLy2. The purpose of this study was to evaluate adaptive B‐cell responses towards the thymus‐dependent vaccine PCV13 in SLy2‐overexpressing mice and to study their survival rate during pneumococcal lung infection. Moreover, we investigated B‐cell developmental stages within the bone marrow (BM) in the context of excessive SLy2‐expression.

Methods

B‐cell subsets and their surface immune globulins were investigated by flow cytometry. For class‐switch assays, isolated splenic B cells were stimulated in vitro with lipopolysaccharide and interleukin‐4 and antibody secretion was quantified via LEGENDplex. To study PCV13‐specific responses, mice were immunized and serum antibody titers (immunoglobulin M, immunoglobulins IgG₁, IgG₂, and IgG₃) were examined by enzyme‐linked immunosorbent assay. Survival rates of mice were assessed within 7 days upon intranasal challenge with S. pneumoniae.

Results

Our data demonstrate impaired IgG₁ and IgG₃ antibody responses towards the pneumococcal conjugate‐vaccine PCV13 in SLy2‐overexpressing mice. This was accompanied by reduced frequencies and numbers of BM‐resident plasmablasts. In addition, we found drastically reduced counts of B‐cell precursors in the BM of SLy2‐Tg mice. The survival rate upon intranasal challenge with S. pneumoniae was mostly comparable between the genotypes.

Conclusion

Our findings demonstrate an important role of the adapter protein SLy2 in the context of adaptive antibody responses against pneumococcal conjugate‐vaccine. Interestingly, deficits in humoral immunity seemed to be compensated by cellular immune effectors upon bacterial challenge. Our study further shows a novel relevance of SLy2 for plasmablasts and B‐cell progenitors in the BM.

Mechanism and regulation of class switch recombination by IgH transcriptional control elements

Class switch recombination (CSR) plays an important role in humoral immunity by generating antibodies with different effector functions. CSR to a particular antibody isotype is induced by external stimuli, and occurs between highly repetitive switch (S) sequences. CSR requires transcription across S regions, which generates long non-coding RNAs and secondary structures that promote accessibility of S sequences to activation-induced cytidine deaminase (AID). AID initiates DNA double-strand breaks (DSBs) intermediates that are repaired by general DNA repair pathways. Switch transcription is controlled by various regulatory elements, including enhancers and insulators. The current paradigm posits that transcriptional control of CSR involves long-range chromatin interactions between regulatory elements and chromatin loops-stabilizing factors, which promote alignment of partner S regions in a CSR centre (CSRC) and initiation of CSR. In this review, we focus on the role of IgH transcriptional control elements in CSR and the chromatin-based mechanisms underlying this control.

TADsplimer reveals splits and mergers of topologically associating domains for epigenetic regulation of transcription

We present TADsplimer, the first computational tool to systematically detect topologically associating domain (TAD) splits and mergers across the genome between Hi-C samples. TADsplimer recaptures splits and mergers of TADs with high accuracy in simulation analyses and defines hundreds of TAD splits and mergers between pairs of different cell types, such as endothelial cells and fibroblasts. Our work reveals a key role for TAD remodeling in epigenetic regulation of transcription and delivers the first tool for the community to perform dynamic analysis of TAD splits and mergers in numerous biological and disease models.

Aryl hydrocarbon receptor-induced activation of the human IGH hs1.2 enhancer: Mutational analysis of putative regulatory binding motifs

The human hs1.2 enhancer within the Ig heavy chain gene (IGH) is polymorphic and associated with a number of autoimmune diseases. The polymorphic region is characterized by tandem repeats of an ∼53-bp invariant sequence containing possible binding sites for several transcription factors. Our previous studies suggest the human hs1.2 enhancer is sensitive to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), an environmental toxicant and high affinity ligand of the aryl hydrocarbon receptor (AhR). TCDD induced hs1.2 enhancer activity in an AhR-dependent manner and the number of invariant sequences influenced the magnitude of activity. To better understand the regulation of human hs1.2 enhancer activity, the objective of the current study was to utilize mutational analysis and luciferase reporter constructs to evaluate the contribution of putative transcription factor binding sites to overall hs1.2 enhancer activity and modulation by TCDD. Basal and LPS-induced activity of the hs1.2 enhancer appeared to be most affected by mutation of sites outside of the invariant sequence or deletion of the entire invariant sequence; whereas sites influencing the effect of TCDD were dependent on the cellular activation state (i.e. unstimulated vs. LPS stimulation) and relatively independent of the putative AhR binding site within the invariant sequence. These results suggest that AhR activation affects human hs1.2 activity through an as yet undetermined non-canonical pathway. A better understanding regarding the role of the hs1.2 enhancer in human Ig expression and how AhR ligands modulate its activity may lead to insights into overall Ig regulation and mechanisms of dysfunction.

Charting the cis-regulome of activated B cells by coupling structural and functional genomics

Cis-regulomes underlying immune-cell-specific genomic states have been extensively analyzed by structure-based chromatin profiling. By coupling such approaches with a high-throughput enhancer screen (self-transcribing active regulatory region sequencing (STARR-seq)), we assembled a functional cis-regulome for lipopolysaccharide-activated B cells. Functional enhancers, in contrast with accessible chromatin regions that lack enhancer activity, were enriched for enhancer RNAs (eRNAs) and preferentially interacted in vivo with B cell lineage-determining transcription factors. Interestingly, preferential combinatorial binding by these transcription factors was not associated with differential enrichment of their sites. Instead, active enhancers were resolved by principal component analysis (PCA) from all accessible regions by co-varying transcription factor motif scores involving a distinct set of signaling-induced transcription factors. High-resolution chromosome conformation capture (Hi-C) analysis revealed multiplex, activated enhancer-promoter configurations encompassing numerous multi-enhancer genes and multi-genic enhancers engaged in the control of divergent molecular pathways. Motif analysis of pathway-specific enhancers provides a catalog of diverse transcription factor codes for biological processes encompassing B cell activation, cycling and differentiation.

Comparison of computational methods for 3D genome analysis at single-cell Hi-C level

Hi-C is a high-throughput chromosome conformation capture technology that is becoming routine in the literature. Although the price of sequencing has been dropping dramatically, high-resolution Hi-C data are not always an option for many studies, such as in single cells. However, the performance of current computational methods based on Hi-C at the ultra-sparse data condition has yet to be fully assessed. Therefore, in this paper, after briefly surveying the primary computational methods for Hi-C data analysis, we assess the performance of representative methods on data normalization, identification of compartments, Topologically Associating Domains (TADs) and chromatin loops under the condition of ultra-low resolution. We showed that most state-of-the-art methods do not work properly for that condition. Then, we applied the three best-performing methods on real single-cell Hi-C data, and their performance indicates that compartments may be a statistical feature emerging from the cell population, while TADs and chromatin loops may dynamically exist in single cells.

Long-read sequencing unveils IGH-DUX4 translocation into the silenced IGH allele in B-cell acute lymphoblastic leukemia

IGH@ proto-oncogene translocation is a common oncogenic event in lymphoid lineage cancers such as B-ALL, lymphoma and multiple myeloma. Here, to investigate the interplay between IGH@ proto-oncogene translocation and IGH allelic exclusion, we perform long-read whole-genome and transcriptome sequencing along with epigenetic and 3D genome profiling of Nalm6, an IGH-DUX4 positive B-ALL cell line. We detect significant allelic imbalance on the wild-type over the IGH-DUX4 haplotype in expression and epigenetic data, showing IGH-DUX4 translocation occurs on the silenced IGH allele. In vitro, this reduces the oncogenic stress of DUX4 high-level expression. Moreover, patient samples of IGH-DUX4 B-ALL have similar expression profile and IGH breakpoints as Nalm6, suggesting a common mechanism to allow optimal dosage of non-toxic DUX4 expression.

The IGH@ proto-oncogene translocation is a known genomic driver in several blood cancers. Here, the authors show that IGH-DUX4 translocation occurs on the silenced IGH allele avoiding toxic high-level expression of DUX4 in B-ALL.

Association of GTF2IRD1-GTF2I polymorphisms with neuromyelitis optica spectrum disorders in Han Chinese patients

Variants at the GTF2I repeat domain containing 1 (GTF2IRD1)–GTF2I locus are associated with primary Sjögren’s syndrome, systemic lupus erythematosus, and rheumatoid arthritis. Numerous studies have indicated that this susceptibility locus is shared by multiple autoimmune diseases. However, until now there were no studies of the correlation between GTF2IRD1–GTF2I polymorphisms and neuromyelitis optica spectrum disorders (NMOSD). This case control study assessed this association by recruiting 305 participants with neuromyelitis optica spectrum disorders and 487 healthy controls at the Department of Neurology, from September 2014 to April 2017. Peripheral blood was collected, DNA extracteds and the genetic association between GTF2IRD1–GTF2I polymorphisms and neuromyelitis optica spectrum disorders in the Chinese Han population was analyzed by genotyping. We found that the T allele of rs117026326 was associated with an increased risk of neuromyelitis optica spectrum disorders (odds ratio (OR) = 1.364, 95% confidence interval (CI) 1.019–1.828; P = 0.037). This association persisted after stratification analysis for aquaporin-4 immunoglobulin G antibodies (AQP4-IgG) positivity (OR = 1.397, 95% CI 1.021–1.912; P = 0.036) and stratification according to coexisting autoimmune diseases (OR = 1.446, 95% CI 1.072–1.952; P = 0.015). Furthermore, the CC genotype of rs73366469 was frequent in AQP4-IgG-seropositive patients (OR = 3.15, 95% CI 1.183–8.393, P = 0.022). In conclusion, the T allele of rs117026326 was associated with susceptibility to neuromyelitis optica spectrum disorders, and the CC genotype of rs73366469 conferred susceptibility to AQP4-IgG-seropositivity in Han Chinese patients. The protocol was approved by the Ethics Committee of West China Hospital of Sichuan University, China (approval number: 2016-31) on March 2, 2016.

Discovery platform for inhibitors of IgH gene enhancer activity

Immunoglobulin heavy chain (IgH) translocations are common and early oncogenic events in B cell and plasma cell malignancies including B cell non-Hodgkin's lymphoma (NHL) and multiple myeloma (MM). IgH translocations bring oncogenes into close proximity with potent enhancer elements within the IgH locus, leading to oncogene up-regulation. As IgH enhancer activity is tightly controlled by B cell lineage-specific signaling and transcriptional networks, we hypothesized that IgH enhancers are potentially druggable targets/elements. To test this, we developed a molecular imaging-based high-throughput screening platform for discovering inhibitors of IgH enhancer-driven transcriptional activity. As proof of concept, we identified a low micromolar potency molecule (compound 30666) that inhibited immunoglobulin production by MM cells and blocked expression of an array of IgH translocation-induced oncogenes (CCND1, FGFR3/MMSET, and MYC) in MM and NHL cell lines. Prolonged exposure to 30666 significantly reduced the viability of IgH translocation-positive NHL and MM cells, but was less effective against cells lacking IgH translocations. Compound 30666 exhibited suitable pharmacological properties, including metabolic stability in liver microsomes and oral bioavailability in mice, and demonstrated preclinical anti-MM activity in a plasmacytoma mouse model. Our work suggests that IgH enhancers are attractive and potentially druggable targets for IgH translocation driven malignancies.

Decoding topologically associating domains with ultra-low resolution Hi-C data by graph structural entropy

Submegabase-size topologically associating domains (TAD) have been observed in high-throughput chromatin interaction data (Hi-C). However, accurate detection of TADs depends on ultra-deep sequencing and sophisticated normalization procedures. Here we propose a fast and normalization-free method to decode the domains of chromosomes (deDoc) that utilizes structural information theory. By treating Hi-C contact matrix as a representation of a graph, deDoc partitions the graph into segments with minimal structural entropy. We show that structural entropy can also be used to determine the proper bin size of the Hi-C data. By applying deDoc to pooled Hi-C data from 10 single cells, we detect megabase-size TAD-like domains. This result implies that the modular structure of the genome spatial organization may be fundamental to even a small cohort of single cells. Our algorithms may facilitate systematic investigations of chromosomal domains on a larger scale than hitherto have been possible.

Characteristic arrangement of nucleosomes is predictive of chromatin interactions at kilobase resolution

High-throughput chromosome conformation capture (3C) technologies, such as Hi-C, have made it possible to survey 3D genome structure. However, obtaining 3D profiles at kilobase resolution at low cost remains a major challenge. Therefore, we herein present an algorithm for precise identification of chromatin interaction sites at kilobase resolution from MNase-seq data, termed chromatin interaction site detector (CISD), and a CISD-based chromatin loop predictor (CISD_loop) that predicts chromatin–chromatin interactions (CCIs) from low-resolution Hi-C data. We show that the predictions of CISD and CISD_loop overlap closely with chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) anchors and loops, respectively. The validity of CISD/CISD_loop was further supported by a 3C assay at about 5 kb resolution. Finally, we demonstrate that only modest amounts of MNase-seq and Hi-C data are sufficient to achieve ultrahigh resolution CCI maps. Our results suggest that CCIs may result in characteristic nucleosomes arrangement patterns flanking the interaction sites, and our algorithms may facilitate precise and systematic investigations of CCIs on a larger scale than hitherto have been possible.

3CDB: a manually curated database of chromosome conformation capture data

Chromosome conformation capture (3C) is a biochemical technology to analyse contact frequencies between selected genomic sites in a cell population. Its recent genomic variants, e.g. Hi-C/ chromatin interaction analysis by paired-end tag (ChIA-PET), have enabled the study of nuclear organization at an unprecedented level. However, due to the inherent low resolution and ultrahigh cost of Hi-C/ChIA-PET, 3C is still the gold standard for determining interactions between given regulatory DNA elements, such as enhancers and promoters. Therefore, we developed a database of 3C determined functional chromatin interactions (3CDB;http://3cdb.big.ac.cn). To construct 3CDB, we searched PubMed and Google Scholar with carefully designed keyword combinations and retrieved more than 5000 articles from which we manually extracted 3319 interactions in 17 species. Moreover, we proposed a systematic evaluation scheme for data reliability and classified the interactions into four categories. Contact frequencies are not directly comparable as a result of various modified 3C protocols employed among laboratories. Our evaluation scheme provides a plausible solution to this long-standing problem in the field. A user-friendly web interface was designed to assist quick searches in 3CDB. We believe that 3CDB will provide fundamental information for experimental design and phylogenetic analysis, as well as bridge the gap between molecular and systems biologists who must now contend with noisy high-throughput data.Database URL:http://3cdb.big.ac.cn.

Related Mechanisms of Antibody Somatic Hypermutation and Class Switch Recombination

The primary antibody repertoire is generated by mechanisms involving the assembly of the exons that encode the antigen-binding variable regions of immunoglobulin heavy (IgH) and light (IgL) chains during the early development of B lymphocytes. After antigen-dependent activation, mature B lymphocytes can further alter their IgH and IgL variable region exons by the process of somatic hypermutation (SHM), which allows the selection of B cells in which SHMs resulted in the production of antibodies with increased antigen affinity. In addition, during antigen-dependent activation, B cells can also change the constant region of their IgH chain through a DNA double-strand-break (DSB) dependent process referred to as IgH class switch recombination (CSR), which generates B cell progeny that produce antibodies with different IgH constant region effector functions that are best suited for a elimination of a particular pathogen or in a particular setting. Both the mutations that underlie SHM and the DSBs that underlie CSR are initiated in target genes by activation-induced cytidine deaminase (AID). This review describes in depth the processes of SHM and CSR with a focus on mechanisms that direct AID cytidine deamination in activated B cells and mechanisms that promote the differential outcomes of such cytidine deamination.

IGHV1-69 polymorphism modulates anti-influenza antibody repertoires, correlates with IGHV utilization shifts and varies by ethnicity

IGHV polymorphism provides a rich source of humoral immune system diversity. One important example is the IGHV1-69 germline gene where the biased use of alleles that encode the critical CDR-H2 Phe54 (F-alleles) to make broadly neutralizing antibodies (HV1-69-sBnAb) to the influenza A hemagglutinin stem domain has been clearly established. However, whether IGHV1-69 polymorphism can also modulate B cell function and Ab repertoire expression through promoter and copy number (CN) variations has not been reported, nor has whether IGHV1-69 allelic distribution is impacted by ethnicity. Here we studied a cohort of NIH H5N1 vaccinees and demonstrate for the first time the influence of IGHV1-69 polymorphism on V-segment usage, somatic hypermutation and B cell expansion that elucidates the dominance of F-alleles in HV1-69-sBnAbs. We provide evidence that Phe54/Leu54 (F/L) polymorphism correlates with shifted repertoire usage of other IGHV germline genes. In addition, we analyzed ethnically diverse individuals within the 1000 genomes project and discovered marked variations in F- and L- genotypes and CN among the various ethnic groups that may impact HV1-69-sBnAb responses. These results have immediate implications for understanding HV1-69-sBnAb responses at the individual and population level and for the design and implementation of "universal" influenza vaccine.

Long-Range Control of V(D)J Recombination & Allelic Exclusion: Modeling Views

Allelic exclusion of immunoglobulin (Ig) and T-cell receptor (TCR) genes ensures the development of B and T lymphocytes operating under the mode of clonal selection. This phenomenon associates asynchronous V(D)J recombination events at Ig or TCR alleles and inhibitory feedback control. Despite years of intense research, however, the mechanisms that sustain asymmetric choice in random Ig/TCR dual allele usage and the production of Ig/TCR monoallelic expressing B and T lymphocytes remain unclear and open for debate. In this chapter, we first recapitulate the biological evidence that almost from the start appeared to link V(D)J recombination and allelic exclusion. We review the theoretical models previously proposed to explain this connection. Finally, we introduce our own mathematical modeling views based on how the developmental dynamics of individual lymphoid cells combine to sustain allelic exclusion.

Domain-Specific and Stage-Intrinsic Changes in Tcrb Conformation during Thymocyte Development

Considerable cross-talk exists between mechanisms controlling genome architecture and gene expression. AgR loci are excellent models for these processes because they are regulated at both conformational and transcriptional levels to facilitate their assembly by V(D)J recombination. Upon commitment to the double-negative stage of T cell development, Tcrb adopts a compact conformation that promotes long-range recombination between Vβ gene segments (Trbvs) and their DβJβ targets. Formation of a functional VβDβJβ join signals for robust proliferation of double-negative thymocytes and their differentiation into double-positive (DP) cells, where Trbv recombination is squelched (allelic exclusion). DP differentiation also is accompanied by decontraction of Tcrb, which has been thought to separate the entire Trbv cluster from DβJβ segments (spatial segregation-based model for allelic exclusion). However, DP cells also repress transcription of unrearranged Trbvs, which may contribute to allelic exclusion. We performed a more detailed study of developmental changes in Tcrb topology and found that only the most distal portion of the Trbv cluster separates from DβJβ segments in DP thymocytes, leaving most Trbvs spatially available for rearrangement. Preferential dissociation of distal Trbvs is independent of robust proliferation or changes in transcription, chromatin, or architectural factors, which are coordinately regulated across the entire Trbv cluster. Segregation of distal Trbvs also occurs on alleles harboring a functional VβDβJβ join, suggesting that this process is independent of rearrangement status and is DP intrinsic. Our finding that most Trbvs remain associated with DβJβ targets in DP cells revises allelic exclusion models from their current conformation-dominant to a transcription-dominant formulation.

The immunoglobulin heavy chain VH6-1 promoter regulates Ig transcription in non-B cells

Background

Non-B cell immunoglobulins (Igs) are widely expressed in epithelial cancer cells. The past 20 years of research have demonstrated that non-B cell Igs are associated with cancer cell proliferation, the cellular cytoskeleton and cancer stem cells. In this study we explored the transcriptional mechanism of IgM production in non-B cells.

Methods

The promoter region of a V-segment of the heavy mu chain gene (VH6-1) was cloned from a colon cancer cell line HT-29. Next, the promoter activities in non-B cells and B-cells were detected using the dual-luciferase reporter assay. Then the transcription factor binding to the promoter regions was evaluated by electrophoretic mobility shift assays (EMSAs) and gel supershift experiments.

Results

Our data showed that the sequence 1200 bp upstream of VH6-1 exhibited promoter activity in both B and non-B cells. No new regulatory elements were identified within the region 1200 bp to 300 bp upstream of VH6-1. In addition, Oct-1 was found to bind to the octamer element of the Ig gene promoter in cancer cells, in contrast to B cells, which utilize the transcriptional factor Oct-2.

Conclusion

The regulatory mechanisms among different cell types controlling the production of IgM heavy chains are worth discussing.

The transcription-splicing protein NonO/p54nrb and three NonO-interacting proteins bind to distal enhancer region and augment rhodopsin expression

Phototransduction machinery in vertebrate photoreceptors is contained within the membrane discs of outer segments. Daily renewal of 10% of photoreceptor outer segments requires stringent control of gene expression. Rhodopsin constitutes over 90% of the protein in rod discs, and its altered expression or transport is associated with photoreceptor dysfunction and/or death. Two cis-regulatory sequences have been identified upstream of the rhodopsin transcription start site. While the proximal promoter binds to specific transcription factors, including NRL and CRX, the rhodopsin enhancer region (RER) reportedly contributes to precise and high-level expression of rhodopsin in vivo. Here, we report the identification of RER-bound proteins by mass spectrometry. We validate the binding of NonO (p54(nrb)), a protein implicated in coupling transcription to splicing, and three NonO-interacting proteins-hnRNP M, Ywhaz and Ppp1ca. NonO and its interactors can activate rhodopsin promoter in HEK293 cells and function synergistically with NRL and CRX. DNA-binding domain of NonO is critical for rhodopsin promoter activation. Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) analysis demonstrates high occupancy of NonO at rhodopsin and a subset of phototransduction genes. Furthermore, shRNA knockdown of NonO in mouse retina leads to loss of rhodopsin expression and rod cell death, which can be partially rescued by a C-terminal NonO construct. RNA-seq analysis of the NonO shRNA-treated retina revealed splicing defects and altered expression of genes, specifically those associated with phototransduction. Our studies identify an important contribution of NonO and its interacting modulator proteins in enhancing rod-specific gene expression and controlling rod homeostasis.

Detection of regulatory SNPs in human genome using ChIP-seq ENCODE data

A vast amount of SNPs derived from genome-wide association studies are represented by non-coding ones, therefore exacerbating the need for effective identification of regulatory SNPs (rSNPs) among them. However, this task remains challenging since the regulatory part of the human genome is annotated much poorly as opposed to coding regions. Here we describe an approach aggregating the whole set of ENCODE ChIP-seq data in order to search for rSNPs, and provide the experimental evidence of its efficiency. Its algorithm is based on the assumption that the enrichment of a genomic region with transcription factor binding loci (ChIP-seq peaks) indicates its regulatory function, and thereby SNPs located in this region are more likely to influence transcription regulation. To ensure that the approach preferably selects functionally meaningful SNPs, we performed enrichment analysis of several human SNP datasets associated with phenotypic manifestations. It was shown that all samples are significantly enriched with SNPs falling into the regions of multiple ChIP-seq peaks as compared with the randomly selected SNPs. For experimental verification, 40 SNPs falling into overlapping regions of at least 7 TF binding loci were selected from OMIM. The effect of SNPs on the binding of the DNA fragments containing them to the nuclear proteins from four human cell lines (HepG2, HeLaS3, HCT-116, and K562) has been tested by EMSA. A radical change in the binding pattern has been observed for 29 SNPs, besides, 6 more SNPs also demonstrated less pronounced changes. Taken together, the results demonstrate the effective way to search for potential rSNPs with the aid of ChIP-seq data provided by ENCODE project.

Noncoding transcription at enhancers: general principles and functional models

Mammalian genomes are extensively transcribed outside the borders of protein-coding genes. Genome-wide studies recently demonstrated that cis-regulatory genomic elements implicated in transcriptional control, such as enhancers and locus-control regions, represent major sites of extragenic noncoding transcription. Enhancer-templated transcripts provide a quantitatively small contribution to the total amount of cellular nonribosomal RNA; nevertheless, the possibility that enhancer transcription and the resulting enhancer RNAs may, in some cases, have functional roles, rather than represent mere transcriptional noise at accessible genomic regions, is supported by an increasing amount of experimental data. In this article we review the current knowledge on enhancer transcription and its functional implications.

Bidirectional crosstalk between platelets and monocytes initiated by Toll-like receptor: an important step in the early defense against fungal infections?

Monocytes are important in the defense against fungal infections due to their phagocytic and immunoregulatory functions. Platelets also contribute in such immune responses through their release of soluble mediators, including chemokines as well as several other soluble mediators. Both monocytes and platelets express several Toll-like receptors (TLRs) that can recognize fungal molecules and thus initiate intracellular signaling events. TLR ligation on monocytes and platelets may thereby be an early immunological event and function as an initiator of a local proinflammatory crosstalk between platelets and monocytes resulting in (i) monocyte-induced increase of platelet activation and (ii) platelet-associated enhancement of the monocyte activation/function. These effects may have clinical implications both for the efficiency of antifungal treatment and for the predisposition to fungal infections, for example, increased predisposition in patients with thrombocytopenia/monocytopenia due to chemotherapy- or disease-induced bone marrow failure.

Biochemistry and biology of the inducible multifunctional transcription factor TFII-I: 10 years later

Exactly twenty years ago TFII-I was discovered as a biochemical entity that was able to bind to and function via a core promoter element called the Initiator (Inr). Since then several different properties of this signal-induced multifunctional factor were discovered. Here I update these ever expanding functions of TFII-I--focusing primarily on the last ten years since the first review appeared in this journal.