Active Genomic Regions Research Articles

Convergent Degenerated Regulatory Elements Associated with Limb Loss in Limbless Amphibians and Reptiles.

Limbs are a defining characteristic of tetrapods, yet numerous taxa, primarily among amphibians and reptiles, have independently lost limbs as an adaptation to new ecological niches. To elucidate the genetic factors contributing to this convergent limb loss, we present a 12 Gb chromosome-level assembly of the Banna caecilian (Ichthyophis bannanicus), a limbless amphibian. Our comparative analysis, which includes the reconstruction of amphibian karyotype evolution, reveals constrained gene length evolution in a subset of developmental genes across 3 large genomes. Investigation of limb development genes uncovered the loss of Grem1 in caecilians and Tulp3 in snakes. Interestingly, caecilians and snakes share a significantly larger number of convergent degenerated conserved noncoding elements than limbless lizards, which have a shorter evolutionary history of limb loss. These convergent degenerated conserved noncoding elements overlap significantly with active genomic regions during mouse limb development and are conserved in limbed species, suggesting their essential role in limb patterning in the tetrapod common ancestor. While most convergent degenerated conserved noncoding elements emerged in the jawed vertebrate ancestor, coinciding with the origin of paired appendage, more recent degenerated conserved noncoding elements also contribute to limb development, as demonstrated through functional experiments. Our study provides novel insights into the regulatory elements associated with limb development and loss, offering an evolutionary perspective on the genetic basis of morphological specialization.

Recruitment of homodimeric proneural factors by conserved CAT-CAT E-boxes drives major epigenetic reconfiguration in cortical neurogenesis.

Proneural factors of the basic helix-loop-helix family coordinate neurogenesis and neurodifferentiation. Among them, NEUROG2 and NEUROD2 subsequently act to specify neurons of the glutamatergic lineage. Disruption of these factors, their target genes and binding DNA motifs has been linked to various neuropsychiatric disorders. Proneural factors bind to specific DNA motifs called E-boxes (hexanucleotides of the form CANNTG, composed of two CAN half sites on opposed strands). While corticogenesis heavily relies on E-box activity, the collaboration of proneural factors on different E-box types and their chromatin remodeling mechanisms remain largely unknown. Here, we conducted a comprehensive analysis using chromatin immunoprecipitation followed bysequencing (ChIP-seq) data for NEUROG2 and NEUROD2, along with time-matched single-cell RNA-seq, ATAC-seq and DNA methylation data from the developing mouse cortex. Our findings show that these factors are highly enriched in transiently active genomic regions during intermediate stages of neuronal differentiation. Although they primarily bind CAG-containing E-boxes, their binding in dynamic regions is notably enriched in CAT-CAT E-boxes (i.e. CATATG, denoted as 5'3' half sites for dimers), which undergo significant DNA demethylation and exhibit the highest levels of evolutionary constraint. Aided by HT-SELEX data reanalysis, structural modeling and DNA footprinting, we propose that these proneural factors exert maximal chromatin remodeling influence during intermediate stages of neurogenesis by binding as homodimers to CAT-CAT motifs. This study provides an in-depth integrative analysis of the dynamic regulation of E-boxes during neuronal development, enhancing our understanding of the mechanisms underlying the binding specificity of critical proneural factors.

Rapid degradation of DHX36 revealing its transcriptional role by interacting with G‐quadruplex

AbstractAccumulating evidence indicates that G‐quadruplexes (G4s) are involved in transcriptional regulation. Previous studies have demonstrated that DHX36 preferentially resolves G4s, suggesting its potential impact on gene transcription mediated by these structures. However, systematic validation is required to establish a link between DHX36 activity and its roles in transcriptional regulation. In this study, we investigate the role of DHX36 in transcription. First, we employ the cleavage under targets and tagmentation (CUT&Tag), an efficient method for mapping protein–DNA interactions, to identify the binding sites in the chromatin of MCF‐7 cells. Subsequently, we use the auxin‐inducible degron (AID) protein degradation system and improved nascent RNA sequencing method acrylonitrile‐mediated uridine‐to‐cytidine conversion sequencing (AMUC‐seq) to pinpoint genes directly regulated by DHX36. Our results reveal a significant enrichment of G4 structures at DHX36 target sites, predominantly located in active genomic regions. In vitro assays further demonstrate DHX36's interaction with G4 sequences from three specific oncogenes. These findings underscore the potential role of DHX36 in modulating gene transcription through G4 structures.

Involvement of the splicing factor SART1 in the BRCA1-dependent homologous recombination repair of DNA double-strand breaks

Although previous studies have reported that pre-mRNA splicing factors (SFs) are involved in the repair of DNA double-strand breaks (DSBs) via homologous recombination (HR), their exact role in promoting HR remains poorly understood. Here, we showed that SART1, an SF upregulated in several types of cancer, promotes DSB end resection, an essential first step of HR. The resection-promoting function of SART1 requires phosphorylation at threonine 430 and 695 by ATM/ATR. SART1 is recruited to DSB sites in a manner dependent on transcription and its RS domain. SART1 is epistatic with BRCA1, a major HR factor, in the promotion of resection, especially transcription-associated resection in the G2 phase. SART1 and BRCA1 accumulate at DSB sites in an interdependent manner, and epistatically counteract the resection blockade posed by 53BP1 and RIF1. Furthermore, chromosome analysis demonstrated that SART1 and BRCA1 epistatically suppressed genomic alterations caused by DSB misrepair in the G2 phase. Collectively, these results indicate that SART1 and BRCA1 cooperatively facilitate resection of DSBs arising in transcriptionally active genomic regions in the G2 phase, thereby promoting faithful repair by HR, and suppressing genome instability.

DIPG-05. CHARACTERIZATION OF THE EFFECTS OF H3K27M VARIANTS ON THE DNA-DAMAGE RESPONSE AND TUMOR EVOLUTION IN DIFFUSE MIDLINE GLIOMA

Abstract Diffuse midline gliomas (DMGs) are universally-fatal pediatric brain tumors driven by K27M alterations in histone 3 (H3K27M). Histone 3 is encoded by several variants, including the canonical H3.1, and H3.3 which replaces H3.1 at transcriptionally active and repetitive genomic regions. Both H3.3K27M and H3.1K27M are present as driver events in DMGs and recent single-cell-RNA and whole-genome sequencing experiments have demonstrated that they define distinct subsets of disease, with varying secondary genetic alterations and different gene-expression patterns. The cause of the divergence in routes of tumor evolution between these subsets remains an open question. Previous work by our group has shown that H3.3K27M DMGs show greater genomic instability than H3.1K27M tumors – with a greater number and higher complexity of genomic structural variants resulting from errors in double-stranded DNA break repair. This raises the possibility that histone variants differentially influence the generation of, or response to, DNA breaks. Here we present an early-stage project that aims to investigate whether H3K27M variants are sufficient to produce differences in genomic instability and in directing the path of tumor evolution. Using isogenic cell-line models expressing H3.1K27M or H3.3K27M, we aim to characterize whether H3K27M variants produce varying levels of replicative stress and P53 activation, leading to varying selective pressures for P53 loss. In turn, we seek to demonstrate that loss of P53 produces a greater growth advantage in H3.3K27M cells, and that P53 loss leads to increased tolerance of genomic instability in DMG models. Together, these experiments seek to provide new insight on the mechanisms of divergence between H3.1K27M and H3.3K27M-driven DMGs, and a better understanding of the sources of genomic instability in DMGs. Given that ionizing radiation remains the only approved treatment modality for DMGs, a greater understanding of DMG response to DNA-damage is critical to improving therapeutic outcomes.

The BAF chromatin remodeler synergizes with RNA polymerase II and transcription factors to evict nucleosomes

Chromatin accessibility is a hallmark of active transcription and entails ATP-dependent nucleosome remodeling, which is carried out by complexes such as Brahma-associated factor (BAF). However, the mechanistic links between transcription, nucleosome remodeling and chromatin accessibility are unclear. Here, we used a chemical–genetic approach coupled with time-resolved chromatin profiling to dissect the interplay between RNA Polymerase II (RNAPII), BAF and DNA-sequence-specific transcription factors in mouse embryonic stem cells. We show that BAF dynamically unwraps and evicts nucleosomes at accessible chromatin regions, while RNAPII promoter-proximal pausing stabilizes BAF chromatin occupancy and enhances ATP-dependent nucleosome eviction by BAF. We find that although RNAPII and BAF dynamically probe both transcriptionally active and Polycomb-repressed genomic regions, pluripotency transcription factor chromatin binding confers locus specificity for productive chromatin remodeling and nucleosome eviction by BAF. Our study suggests a paradigm for how functional synergy between dynamically acting chromatin factors regulates locus-specific nucleosome organization and chromatin accessibility.

Blimp-1 and c-Maf regulate Il10 andnegatively regulate common and unique proinflammatory gene networks in IL-12 plus IL-27-driven T helper-1 cells.

CD4 + Th1 cells producing IFN-γ are required to eradicate intracellular pathogens, however if uncontrolled these cells can cause immunopathology. The cytokine IL-10 is produced by multiple immune cells including Th1 cells during infection and regulates the immune response to minimise collateral host damage. In this study we aimed to elucidate the transcriptional network of genes controlling the expression of Il10 and proinflammatory cytokines, including Ifng in Th1 cells differentiated from mouse naive CD4 + T cells. We applied computational analysis of gene regulation derived from temporal profiling of gene expression clusters obtained from bulk RNA sequencing (RNA-seq) of flow cytometry sorted naïve CD4 + T cells from mouse spleens differentiated in vitro into Th1 effector cells with IL-12 and IL-27 to produce Ifng and Il10, compared to IL-27 alone which express Il10 only , or IL-12 alone which express Ifng and no Il10, or medium control driven-CD4 + T cells which do not express effector cytokines . Data were integrated with analysis of active genomic regions from these T cells using an assay for transposase-accessible chromatin with sequencing (ATAC)-seq, integrated with literature derived-Chromatin-immunoprecipitation (ChIP)-seq data and the RNA-seq data, to elucidate the transcriptional network of genes controlling expression of Il10 and pro-inflammatory effector genes in Th1 cells. The co-dominant role for the transcription factors, Prdm1 (encoding Blimp-1) and Maf (encoding c-Maf) , in cytokine gene regulation in Th1 cells, was confirmed using T cells obtained from mice with T-cell specific deletion of these transcription factors. We show that the transcription factors Blimp-1 and c-Maf each have unique and common effects on cytokine gene regulation and not only co-operate to induce Il10 gene expression in IL-12 plus IL-27 differentiated mouse Th1 cells, but additionally directly negatively regulate key proinflammatory cytokines including Ifng, thus providing mechanisms for reinforcement of regulated Th1 cell responses. These data show that Blimp-1 and c-Maf positively and negatively regulate a network of both unique and common anti-inflammatory and pro-inflammatory genes to reinforce a Th1 response in mice that will eradicate pathogens with minimum immunopathology.

Identifying predictors of glioma evolution from longitudinal sequencing.

Clonal evolution drives cancer progression and therapeutic resistance. Recent studies have revealed divergent longitudinal trajectories in gliomas, but early molecular features steering posttreatment cancer evolution remain unclear. Here, we collected sequencing and clinical data of initial-recurrent tumor pairs from 544 adult diffuse gliomas and performed multivariate analysis to identify early molecular predictors of tumor evolution in three diffuse glioma subtypes. We found that CDKN2A deletion at initial diagnosis preceded tumor necrosis and microvascular proliferation that occur at later stages of IDH-mutant glioma. Ki67 expression at diagnosis was positively correlated with acquiring hypermutation at recurrence in the IDH-wild-type glioma. In all glioma subtypes, MYC gain or MYC-target activation at diagnosis was associated with treatment-induced hypermutation at recurrence. To predict glioma evolution, we constructed CELLO2 (Cancer EvoLution for LOngitudinal data version 2), a machine learning model integrating features at diagnosis to forecast hypermutation and progression after treatment. CELLO2 successfully stratified patients into subgroups with distinct prognoses and identified a high-risk patient group featured by MYC gain with worse post-progression survival, from the low-grade IDH-mutant-noncodel subtype. We then performed chronic temozolomide-induction experiments in glioma cell lines and isogenic patient-derived gliomaspheres and demonstrated that MYC drives temozolomide resistance by promoting hypermutation. Mechanistically, we demonstrated that, by binding to open chromatin and transcriptionally active genomic regions, c-MYC increases the vulnerability of key mismatch repair genes to treatment-induced mutagenesis, thus triggering hypermutation. This study reveals early predictors of cancer evolution under therapy and provides a resource for precision oncology targeting cancer dynamics in diffuse gliomas.

Prolyl endopeptidase remodels macrophage function as a novel transcriptional coregulator and inhibits fibrosis

Macrophages are immune cells crucial for host defense and homeostasis maintenance, and their dysregulation is involved in multiple pathological conditions, such as liver fibrosis. The transcriptional regulation in macrophage is indispensable for fine-tuning of macrophage functions, but the details have not been fully elucidated. Prolyl endopeptidase (PREP) is a dipeptidyl peptidase with both proteolytic and non-proteolytic functions. In this study, we found that Prep knockout significantly contributed to transcriptomic alterations in quiescent and M1/M2-polarized bone marrow-derived macrophages (BMDMs), as well as aggravated fibrosis in an experimental nonalcoholic steatohepatitis (NASH) model. Mechanistically, PREP predominantly localized to the macrophage nuclei and functioned as a transcriptional coregulator. Using CUT&Tag and co-immunoprecipitation, we found that PREP was mainly distributed in active cis-regulatory genomic regions and physically interacted with the transcription factor PU.1. Among PREP-regulated downstream genes, genes encoding profibrotic cathepsin B and D were overexpressed in BMDMs and fibrotic liver tissue. Our results indicate that PREP in macrophages functions as a transcriptional coregulator that finely tunes macrophage functions, and plays a protective role against liver fibrosis pathogenesis.

Multiplex-GAM: genome-wide identification of chromatin contacts yields insights overlooked by Hi-C

Technology for measuring 3D genome topology is increasingly important for studying gene regulation, for genome assembly and for mapping of genome rearrangements. Hi-C and other ligation-based methods have become routine but have specific biases. Here, we develop multiplex-GAM, a faster and more affordable version of genome architecture mapping (GAM), a ligation-free technique that maps chromatin contacts genome-wide. We perform a detailed comparison of multiplex-GAM and Hi-C using mouse embryonic stem cells. When examining the strongest contacts detected by either method, we find that only one-third of these are shared. The strongest contacts specifically found in GAM often involve ‘active’ regions, including many transcribed genes and super-enhancers, whereas in Hi-C they more often contain ‘inactive’ regions. Our work shows that active genomic regions are involved in extensive complex contacts that are currently underestimated in ligation-based approaches, and highlights the need for orthogonal advances in genome-wide contact mapping technologies.

Epigenetic plasticity safeguards heterochromatin configuration in mammals.

Heterochromatin is a key architectural feature of eukaryotic chromosomes critical for cell type-specific gene expression and genome stability. In the mammalian nucleus, heterochromatin segregates from transcriptionally active genomic regions and exists in large, condensed, and inactive nuclear compartments. However, the mechanisms underlying the spatial organization of heterochromatin need to be better understood. Histone H3 lysine 9 trimethylation (H3K9me3) and lysine 27 trimethylation (H3K27me3) are two major epigenetic modifications that enrich constitutive and facultative heterochromatin, respectively. Mammals have at least five H3K9 methyltransferases (SUV39H1, SUV39H2, SETDB1, G9a and GLP) and two H3K27 methyltransferases (EZH1 and EZH2). In this study, we addressed the role of H3K9 and H3K27 methylation in heterochromatin organization using a combination of mutant cells for five H3K9 methyltransferases and an EZH1/2 dual inhibitor, DS3201. We showed that H3K27me3, which is normally segregated from H3K9me3, was redistributed to regions targeted by H3K9me3 after the loss of H3K9 methylation and that the loss of both H3K9 and H3K27 methylation resulted in impaired condensation and spatial organization of heterochromatin. Our data demonstrate that the H3K27me3 pathway safeguards heterochromatin organization after the loss of H3K9 methylation in mammalian cells.

Epigenetic modifications of tumor necrosis factor-alpha in joint cartilage tissue from osteoarthritis patients - CONSORT.

Background:Osteoarthritis (OA) remains one of the most common osteopathy for centuries, which can be attributed to multiple risk factors including mechanical and biochemical ones. More and more studies verified that inflammatory cytokines play important roles in the progression of OA, such as tumor necrosis factor-alpha (TNF-α). In this study, we aimed to investigate the relationship between epigenetic manifestations of TNF-? and the pathogenesis of OA.Methods:Totally, 37 OA patients’ cartilage was collected through the knee joint and 13 samples of articular cartilage as healthy control was collected through traumatic amputation. Real-time PCR, Western blot and ELISA analysis were performed to observe the expression of target genes and proteins in collected samples.Results:Compared with the healthy control group, TNF-? was over-expressing in cartilage which was collected from OA patients. DNA hypomethylation, histone hyperacetylation and histone methylation were observed in the TNF-? promoter in OA compared with normal patients, and we also studied series of enzymes associated with epigenetics. The results showed that by increasing DNA methylation and decreasing histone acetylation in the TNF-? promoter, and TNF-? over-expression in OA cartilage was suppressed, histone methylation has no significant correlation with OA.Conclusion:In conclusion, the changes of epigenetic status regulate TNF-α expression in the cells, which are pivotal to the OA disease process. These results may give us a better understanding of OA and may provide new therapeutic options.

Sites of active gene regulation in the prenatal frontal cortex and their role in neuropsychiatric disorders.

Common genetic variation appears to largely influence risk for neuropsychiatric disorders through effects on gene regulation. It is therefore possible to shed light on the biology of these conditions by testing for enrichment of associated genetic variation within regulatory genomic regions operating in specific tissues or cell types. Here, we have used the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-Seq) to map open chromatin (an index of active regulatory genomic regions) in bulk tissue, NeuN+ and NeuN- nuclei from the prenatal human frontal cortex, and tested enrichment of single-nucleotide polymorphism (SNP) heritability for five neuropsychiatric disorders (autism spectrum disorder, attention deficit hyperactivity disorder [ADHD], bipolar disorder, major depressive disorder, and schizophrenia) within these regions. We observed significant enrichment of SNP heritability for ADHD, major depressive disorder, and schizophrenia within open chromatin regions (OCRs) mapped in bulk fetal frontal cortex, and for all five tested neuropsychiatric conditions when we restricted these sites to those overlapping histone modifications indicative of enhancers (H3K4me1) or promoters (H3K4me3) in fetal brain. SNP heritability for neuropsychiatric disorders was significantly enriched in OCRs identified in fetal frontal cortex NeuN- as well as NeuN+ nuclei overlapping fetal brain H3K4me1 or H3K4me3 sites. We additionally demonstrate the utility of our mapped OCRs for prioritizing potentially functional SNPs at genome-wide significant risk loci for neuropsychiatric disorders. Our data provide evidence for an early neurodevelopmental component to a range of neuropsychiatric conditions and highlight an important role for regulatory genomic regions active within both NeuN+ and NeuN- cells of the prenatal brain.

Three-Dimensional Chromatin Architecture Landscape of Aging and Alzheimer's Disease

Recently, accumulating genetic and epigenetic studies are indicating a close association between dysregulated enhancers and Alzheimer's disease (AD). However, their contributions were poorly defined for lacking direct links to disease genes. To bridge this gap, we presented the maps of three-dimensional chromatin architecture of brain tissue of AD patients, dementia-free aged and young individuals. We found that AD patients were more affected by the dynamic of 3D -chromatin architecture. There were four classes of long-distance interaction changes. Among them, the feedback class was more observed, indicating the prevalence of feedback loops during AD genesis. Integrating with epigenetic data, a new strategy was proposed to improve the resolution of Hi-C data. Hence, we observed more detailed interactions among active genomic regions and found that known AD risk genes were under long-distance promoter-enhancer interactions. As an application, we utilized the Hi-C loops to link the AD risk SNPs to target genes and to identify the dysregulated genes. Functional annotation suggested that long-distance interaction changes contributed to AD genesis by disrupting brain function-related pathways. For future AD and aging studies, our result provides a reference landscape to better interpret signals of association and epigenetic studies.

A discrete subset of epigenetically primed human NK cells mediates antigen-specific immune responses.

Adaptive features of natural killer (NK) cells have been reported in various species with different underlying mechanisms. It is unclear, however, which NK cell populations are capable of mounting antigen-specific recall responses and how such functions are regulated at the molecular level. Here, we identify and characterize a discrete population of CD49a+CD16- NK cells in the human liver that displays increased epigenetic potential to elicit memory responses and has the functional properties to exert antigen-specific immunity in the skin as an effector site. Integrated chromatin-based epigenetic and transcriptomic profiling revealed unique characteristics of hepatic CD49a+CD16- NK cells when compared with conventional CD49a-CD16+ NK cells, thereby defining active genomic regions and molecules underpinning distinct NK cell reactivity. In contrast to conventional NK cells, our results suggest that adaptive CD49a+CD16- NK cells are able to bypass the KIR receptor-ligand system upon antigen-specific stimulation. Furthermore, these cells were highly migratory toward chemokine gradients expressed in epicutaneous patch test lesions as an effector site of adaptive immune responses in the skin. These results define pathways operative in human antigen-specific adaptive NK cells and provide a roadmap for harnessing this NK cell subset for specific therapeutic or prophylactic vaccine strategies.

Multiomics Investigation Revealing the Characteristics of HIV-1-Infected Cells In Vivo

For eradication of HIV-1 infection, it is important to elucidate the detailed features and heterogeneity of HIV-1-infected cells invivo. To reveal multiple characteristics of HIV-1-producing cells invivo, we use a hematopoietic-stem-cell-transplanted humanized mouse model infected with GFP-encoding replication-competent HIV-1. We perform multiomics experiments using recently developed technology to identify the features of HIV-1-infected cells. Genome-wide HIV-1 integration-site analysis reveals that productive HIV-1 infection tends to occur in cells with viral integration into transcriptionally active genomic regions. Bulk transcriptome analysis reveals that a high level of viral mRNA is transcribed in HIV-1-infected cells. Moreover, single-cell transcriptome analysis shows the heterogeneity of HIV-1-infected cells, including CXCL13high cells and a subpopulation with low expression of interferon-stimulated genes, which can contribute to efficient viral spread invivo. Our findings describe multiple characteristics of HIV-1-producing cells invivo, which could provide clues for the development of an HIV-1 cure.

Augmented oxidative stress increases 8-oxoguanine preferentially in the transcriptionally active genomic regions

8-Oxoguanine (8-oxoG) is the most common DNA base modification in the mammalian genome, associated with oxidative stress. Here we analysed the alterations in the distribution of 8-oxoG across the entire murine genome, before and after an elevation of oxidative stress by the use of ferric nitrilotriacetate (Fe-NTA) as an oxidative stress inducer in the renal proximal tubules. We isolated DNA fragments containing 8-oxoGs with immunoprecipitation from the murine genome, and amplified them by PCR for a distribution analysis with microarray-based comparative genomic hybridisation. The distribution profiles revealed that frequencies of 8-oxoG fluctuated with a cycle of 1–10 Mb along the chromosomes and the amplitude of the fluctuation was reduced after Fe-NTA administration. The distributions of 8-oxoG along the entire genome in the control and oxidatively stressed conditions were negatively correlated with that of gene density but positively correlated with that of Lamin B1 interaction, which corresponds to lamina-associated domains. These results on the murine genome were consistent with those on the rat genome we previously reported. We further discovered a negative correlation between the distributions of 8-oxoG and transcriptional activity along the genome. Finally, a comparison of the distributions before and after Fe-NTA administration suggested that 8-oxoGs are generated in response to the augmented oxidative stress preferentially in the transcriptionally active genomic regions, where 8-oxoGs have been less accumulated in the control condition.

Efficient and Reproducible Multigene Expression after Single-Step Transfection Using Improved BAC Transgenesis and Engineering Toolkit.

Achieving stable expression of a single transgene in mammalian cells remains challenging; even more challenging is obtaining simultaneous stable expression of multiple transgenes at reproducible, relative expression levels. Previously, we attained copy-number-dependent, chromosome-position-independent expression of reporter minigenes by embedding them within a BAC "scaffold" containing the mouse Msh3-Dhfr locus (DHFR BAC). Here, we extend this "BAC TG-EMBED" approach. First, we report a toolkit of endogenous promoters capable of driving transgene expression over a 0.01- to 5-fold expression range relative to the CMV promoter, allowing fine-tuning of relative expression levels of multiple reporter genes. Second, we demonstrate little variation in expression level and long-term expression stability of a reporter gene embedded in BACs containing either transcriptionally active or inactive genomic regions, making the choice of BAC scaffolds more flexible. Third, we present a novel BAC assembly scheme, "BAC-MAGIC", for inserting multiple transgenes into BAC scaffolds, which is much more time-efficient than traditional galK-based methods. As a proof-of-principle for our improved BAC TG-EMBED toolkit, we simultaneously fluorescently labeled three nuclear compartments at reproducible, relative intensity levels in 94% of stable clones after a single transfection using a DHFR BAC scaffold containing 4 transgenes assembled with BAC-MAGIC. Our extended BAC TG-EMBED toolkit and BAC-MAGIC method provide an efficient, versatile platform for stable simultaneous expression of multiple transgenes at reproducible, relative levels.

New Pathways of Mutational Change in SARS-CoV-2 Proteomes Involve Regions of Intrinsic Disorder Important for Virus Replication and Release.

The massive worldwide spread of the SARS-CoV-2 virus is fueling the COVID-19 pandemic. Since the first whole-genome sequence was published in January 2020, a growing database of tens of thousands of viral genomes has been constructed. This offers opportunities to study pathways of molecular change in the expanding viral population that can help identify molecular culprits of virulence and virus spread. Here we investigate the genomic accumulation of mutations at various time points of the early pandemic to identify changes in mutationally highly active genomic regions that are occurring worldwide. We used the Wuhan NC_045512.2 sequence as a reference and sampled 15 342 indexed sequences from GISAID, translating them into proteins and grouping them by month of deposition. The per-position amino acid frequencies and Shannon entropies of the coding sequences were calculated for each month, and a map of intrinsic disorder regions and binding sites was generated. The analysis revealed dominant variants, most of which were located in loop regions and on the surface of the proteins. Mutation entropy decreased between March and April of 2020 after steady increases at several sites, including the D614G mutation site of the spike (S) protein that was previously found associated with higher case fatality rates and at sites of the NSP12 polymerase and the NSP13 helicase proteins. Notable expanding mutations include R203K and G204R of the nucleocapsid (N) protein inter-domain linker region and G251V of the viroporin encoded by ORF3a between March and April. The regions spanning these mutations exhibited significant intrinsic disorder, which was enhanced and decreased by the N-protein and viroporin 3a protein mutations, respectively. These results predict an ongoing mutational shift from the spike and replication complex to other regions, especially to encoded molecules known to represent major β-interferon antagonists. The study provides valuable information for therapeutics and vaccine design, as well as insight into mutation tendencies that could facilitate preventive control.

Highly structured homolog pairing reflects functional organization of the Drosophila genome

Trans-homolog interactions have been studied extensively in Drosophila, where homologs are paired in somatic cells and transvection is prevalent. Nevertheless, the detailed structure of pairing and its functional impact have not been thoroughly investigated. Accordingly, we generated a diploid cell line from divergent parents and applied haplotype-resolved Hi-C, showing that homologs pair with varying precision genome-wide, in addition to establishing trans-homolog domains and compartments. We also elucidate the structure of pairing with unprecedented detail, observing significant variation across the genome and revealing at least two forms of pairing: tight pairing, spanning contiguous small domains, and loose pairing, consisting of single larger domains. Strikingly, active genomic regions (A-type compartments, active chromatin, expressed genes) correlated with tight pairing, suggesting that pairing has a functional implication genome-wide. Finally, using RNAi and haplotype-resolved Hi-C, we show that disruption of pairing-promoting factors results in global changes in pairing, including the disruption of some interaction peaks.