Heterochromatin protein 1α interacts with parallel RNA and DNA G-quadruplexes

G-quadruplexes in long non-coding RNAs and their interactions with proteins

Long non-coding RNAs (lncRNAs) have emerged as crucial regulators of cellular processes, with their dysregulation linked to various disease states. Among the structural motifs in lncRNAs, RNA G-quadruplexes (rG4s) have gained increasing attention due to their diverse roles in cellular function and disease pathogenesis. This review provides an updated and comprehensive overview of rG4s in lncRNAs, elucidating their formation, interaction with proteins, and distinctive roles in cellular processes. We discuss current methodologies for experimentally probing RNA G4s, including the use of specific small molecules, biomolecular ligands and fluorescent probes. The commonly found RNA G4-interacting protein domains are summarised along with potential strategies for disrupting lncRNA G4-protein interactions from a therapeutic perspective.

Analysis of more than 400,000 women provides case-control evidence for BRCA1 and BRCA2 variant classification

Clinical genetic testing identifies variants causal for hereditary cancer, information that is used for risk assessment and clinical management. Unfortunately, some variants identified are of uncertain clinical significance (VUS), complicating patient management. Case-control data is one evidence type used to classify VUS, and previous findings indicate that case-control likelihood ratios (LRs) outperform odds ratios for variant classification. As an initiative of the Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) Analytical Working Group we analyzed germline sequencing data of BRCA1 and BRCA2 from 96,691 female breast cancer cases and 303,925 unaffected controls from three studies: the BRIDGES study of the Breast Cancer Association Consortium, the Cancer Risk Estimates Related to Susceptibility consortium, and the UK Biobank. We observed 11,227 BRCA1 and BRCA2 variants, with 6,921 being coding, covering 23.4% of BRCA1 and BRCA2 VUS in ClinVar and 19.2% of ClinVar curated (likely) benign or pathogenic variants. Case-control LR evidence was highly consistent with ClinVar assertions for (likely) benign or pathogenic variants; exhibiting 99.1% sensitivity and 95.4% specificity for BRCA1 and 92.2% sensitivity and 86.6% specificity for BRCA2. This approach provides case-control evidence for 785 unclassified variants, that can serve as a valuable element for clinical classification.

Participation of ventricular trabeculae in neonatal cardiac regeneration leads to ectopic recruitment of Purkinje-like cells

Unlike adult mammals, newborn mice can regenerate a functional heart after myocardial infarction; however, the precise origin of the newly formed cardiomyocytes and whether the distal part of the conduction system (the Purkinje fiber (PF) network) is properly formed in regenerated hearts remains unclear. PFs, as well as subendocardial contractile cardiomyocytes, are derived from trabeculae, transient myocardial ridges on the inner ventricular surface. Here, using connexin 40-driven genetic tracing, we uncover a substantial participation of the trabecular lineage in myocardial regeneration through dedifferentiation and proliferation. Concomitantly, regeneration disrupted PF network maturation, resulting in permanent PF hyperplasia and impaired ventricular conduction. Proliferation assays, genetic impairment of PF recruitment, lineage tracing and clonal analysis revealed that PF network hyperplasia results from excessive recruitment of PFs due to increased trabecular fate plasticity. These data indicate that PF network hyperplasia is a consequence of trabeculae participation in myocardial regeneration.

Inverted strand polarity yields thermodynamically stable G-quadruplexes and prevents duplex formation within extended DNA

DNA G-quadruplexes (G4) formed in guanine-rich sequences play a key role in genome function and maintenance, interacting with multiple proteins. However, structural and functional studies of G4s within duplex DNA have been challenging because of the transient nature of G4s and thermodynamic preference of G-rich DNA to form duplexes with their complementary strand rather than G4s. To overcome these challenges, we have incorporated native nucleotides in G-rich sequences using commercially available inverted 3'-O-DMT-5'-O-phosphoramidites of native nucleosides, to give 3'-3' and 5'-5' linkages in the centre of the G-tract. Using circular dichroism and 1H nuclear magnetic resonance spectroscopies and native gel electrophoresis, we demonstrate that these polarity-inverted DNA sequences containing four telomeric repeats form G4s of parallel topology with one lateral or diagonal loop across the face of the quadruplex and two propeller loops across the edges of the quadruplex. These G4s were stable even in the presence of complementary C-rich DNA. As an example, G4 assemblies of inverted polarity were shown to bind to the hinge region of Heterochromatin Protein 1α (HP1α), a known G4-interacting domain. As such, internal polarity inversions in DNA provide a useful tool to control G4 topology while also disrupting the formation of other secondary structures, particularly the canonical duplex.

MOCHA's advanced statistical modeling of scATAC-seq data enables functional genomic inference in large human cohorts

Single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) is being increasingly used to study gene regulation. However, major analytical gaps limit its utility in studying gene regulatory programs in complex diseases. In response, MOCHA (Model-based single cell Open CHromatin Analysis) presents major advances over existing analysis tools, including: 1) improving identification of sample-specific open chromatin, 2) statistical modeling of technical drop-out with zero-inflated methods, 3) mitigation of false positives in single cell analysis, 4) identification of alternative transcription-starting-site regulation, and 5) modules for inferring temporal gene regulatory networks from longitudinal data. These advances, in addition to open chromatin analyses, provide a robust framework after quality control and cell labeling to study gene regulatory programs in human disease. We benchmark MOCHA with four state-of-the-art tools to demonstrate its advances. We also construct cross-sectional and longitudinal gene regulatory networks, identifying potential mechanisms of COVID-19 response. MOCHA provides researchers with a robust analytical tool for functional genomic inference from scATAC-seq data.

Prediction of protein subcellular localization in single cells

The subcellular localization of a protein is important for its function and interaction with other molecules, and its mislocalization is linked to numerous diseases. While atlas-scale efforts have been made to profile protein localization across various cell lines, existing datasets only contain limited pairs of proteins and cell lines which do not cover all human proteins. We present a method that uses both protein sequences and cellular landmark images to perform P redictions of U nseen P roteins' S ubcellular localization ( PUPS ), which can generalize to both proteins and cell lines not used for model training. PUPS combines a protein language model and an image inpainting model to utilize both protein sequence and cellular images for protein localization prediction. The protein sequence input enables generalization to unseen proteins and the cellular image input enables cell type specific prediction that captures single-cell variability. PUPS' ability to generalize to unseen proteins and cell lines enables us to assess the variability in protein localization across cell lines as well as across single cells within a cell line and to identify the biological processes associated with the proteins that have variable localization. Experimental validation shows that PUPS can be used to predict protein localization in newly performed experiments outside of the Human Protein Atlas used for training. Collectively, PUPS utilizes both protein sequences and cellular images to predict protein localization in unseen proteins and cell lines with the ability to capture single-cell variability.

Novel Diagnostic Biomarkers Related to Necroptosis and Immune Infiltration in Coronary Heart Disease

Purpose

Necroptosis, a monitored form of inflammatory cell death, contributes to coronary heart disease (CHD) progression. This study examined the potential of using necroptosis genes as diagnostic markers for CHD and sought to elucidate the underlying roles.

Methods

Through bioinformatic analysis of GSE20680 and GSE20681, we first identified the differentially expressed genes (DEGs) related to necroptosis in CHD. Hub genes were identified using least absolute shrinkage and selection operator (LASSO) regression and random forest analysis after studying immune infiltration and transcription factor-miRNA interaction networks according to the DEGs. Quantitative polymerase chain reaction and immunohistochemistry were used to further investigate hub gene expression in vivo, for which a diagnostic model was constructed and the predictive efficacy was validated. Finally, the CHD group was categorized into high- and low-score groups in accordance with the single-sample gene set enrichment analysis (ssGSEA) score of the necroptosis genes. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, GSEA, and further immune infiltration analyses were performed on the two groups to explore the possible roles of hub genes.

Results

Based on the results of the LASSO regression and random forest analyses, four genes were used to construct a diagnostic model to establish a nomogram. Additionally, an extensive analysis of all seventeen necroptosis genes revealed notable distinctions in expression between high-risk and low-risk groups. Evaluation of immune infiltration revealed that neutrophils, monocytes, B cells, and activated dendritic cells were highly distributed in the peripheral blood of patients with CHD. Specifically, the high CHD score group exhibited greater neutrophil and monocyte infiltration. Conversely, the high-score group showed lower infiltration of M0 and M2 macrophages, CD8+ T, plasma, and resting mast cells.

Conclusion

TLR3, MLKL, HMGB1, and NDRG2 may be prospective biomarkers for CHD diagnosis. These findings offer plausible explanations for the role of necroptosis in CHD progression through immune infiltration and inflammatory response.

Visfatin (NAMPT) affects global gene expression in porcine anterior pituitary cells during the mid-luteal phase of the oestrous cycle

BACKGROUND

The pituitary belongs to the most important endocrine glands involved in regulating reproductive functions. The proper functioning of this gland ensures the undisturbed course of the oestrous cycle and affects the female's reproductive potential. It is believed that visfatin, a hormone belonging to the adipokine family, may regulate reproductive functions in response to the female's metabolic state. Herein we verified the hypothesis that suggests a modulatory effect of visfatin on the anterior pituitary transcriptome during the mid-luteal phase of the oestrous cycle.

RESULTS

RNA-seq analysis of the porcine anterior pituitary cells revealed changes in the expression of 202 genes (95 up-regulated and 107 down-regulated in the presence of visfatin, when compared to the non-treated controls), assigned to 318 gene ontology terms. We revealed changes in the frequency of alternative splicing events (235 cases), as well as long noncoding RNA expression (79 cases) in the presence of the adipokine. The identified genes were associated, among others, with reproductive system development, epithelial cell proliferation, positive regulation of cell development, gland morphogenesis and cell chemotaxis.

CONCLUSIONS

The obtained results indicate a modulatory influence of visfatin on the regulation of the porcine transcriptome and, in consequence, pituitary physiology during the mid-luteal phase of the oestrous cycle.

Evidence for a role of human blood-borne factors in mediating age-associated changes in molecular circadian rhythms

Aging is associated with a number of physiologic changes including perturbed circadian rhythms; however, mechanisms by which rhythms are altered remain unknown. To test the idea that circulating factors mediate age-dependent changes in peripheral rhythms, we compared the ability of human serum from young and old individuals to synchronize circadian rhythms in culture. We collected blood from apparently healthy young (age 25-30) and old (age 70-76) individuals at 14:00 and used the serum to synchronize cultured fibroblasts. We found that young and old sera are equally competent at initiating robust ~24h oscillations of a luciferase reporter driven by clock gene promoter. However, cyclic gene expression is affected, such that young and old sera promote cycling of different sets of genes. Genes that lose rhythmicity with old serum entrainment are associated with oxidative phosphorylation and Alzheimer's Disease as identified by STRING and IPA analyses. Conversely, the expression of cycling genes associated with cholesterol biosynthesis increased in the cells entrained with old serum. Genes involved in the cell cycle and transcription/translation remain rhythmic in both conditions. We did not observe a global difference in the distribution of phase between groups, but found that peak expression of several clock-controlled genes (PER3, NR1D1, NR1D2, CRY1, CRY2, and TEF) lagged in the cells synchronized ex vivo with old serum. Taken together, these findings demonstrate that age-dependent blood-borne factors affect circadian rhythms in peripheral cells and have the potential to impact health and disease via maintaining or disrupting rhythms respectively.

The LIDPAD Mouse Model Captures the Multisystem Interactions and Extrahepatic Complications in MASLD

Metabolic dysfunction-associated steatotic liver disease (MASLD) represents an impending global health challenge. Current management strategies often face setbacks, emphasizing the need for preclinical models that faithfully mimic the human disease and its comorbidities. The liver disease progression aggravation diet (LIDPAD), a diet-induced murine model, extensively characterized under thermoneutral conditions and refined diets is introduced to ensure reproducibility and minimize species differences. LIDPAD recapitulates key phenotypic, genetic, and metabolic hallmarks of human MASLD, including multiorgan communications, and disease progression within 4 to 16 weeks. These findings reveal gut-liver dysregulation as an early event and compensatory pancreatic islet hyperplasia, underscoring the gut-pancreas axis in MASLD pathogenesis. A robust computational pipeline is also detailed for transcriptomic-guided disease staging, validated against multiple harmonized human hepatic transcriptomic datasets, thereby enabling comparative studies between human and mouse models. This approach underscores the remarkable similarity of the LIDPAD model to human MASLD. The LIDPAD model fidelity to human MASLD is further confirmed by its responsiveness to dietary interventions, with improvements in metabolic profiles, liver histopathology, hepatic transcriptomes, and gut microbial diversity. These results, alongside the closely aligned changing disease-associated molecular signatures between the human MASLD and LIDPAD model, affirm the model's relevance and potential for driving therapeutic development.

Integrative analysis of ultra-deep RNA-seq reveals alternative promoter usage as a mechanism of activating oncogenic programmes during prostate cancer progression

Transcription factor (TF) proteins regulate gene activity by binding to regulatory regions, most importantly at gene promoters. Many genes have alternative promoters (APs) bound by distinct TFs. The role of differential TF activity at APs during tumour development is poorly understood. Here we show, using deep RNA sequencing in 274 biopsies of benign prostate tissue, localized prostate tumours and metastatic castration-resistant prostate cancer, that AP usage increases as tumours progress and APs are responsible for a disproportionate amount of tumour transcriptional activity. Expression of the androgen receptor (AR), the key driver of prostate tumour activity, is correlated with elevated AP usage. We identified AR, FOXA1 and MYC as potential drivers of AP activation. DNA methylation is a likely mechanism for AP activation during tumour progression and lineage plasticity. Our data suggest that prostate tumours activate APs to magnify the transcriptional impact of tumour drivers, including AR and MYC.

Obesity-driven mitochondrial dysfunction in human adipose tissue-derived mesenchymal stem/stromal cells involves epigenetic changes

Obesity exacerbates tissue degeneration and compromises the integrity and reparative potential of mesenchymal stem/stromal cells (MSCs), but the underlying mechanisms have not been sufficiently elucidated. Mitochondria modulate the viability, plasticity, proliferative capacity, and differentiation potential of MSCs. We hypothesized that alterations in the 5-hydroxymethylcytosine (5hmC) profile of mitochondria-related genes may mediate obesity-driven dysfunction of human adipose-derived MSCs. MSCs were harvested from abdominal subcutaneous fat of obese and age/sex-matched non-obese subjects (n = 5 each). The 5hmC profile and expression of nuclear-encoded mitochondrial genes were examined by hydroxymethylated DNA immunoprecipitation sequencing (h MeDIP-seq) and mRNA-seq, respectively. MSC mitochondrial structure (electron microscopy) and function, metabolomics, proliferation, and neurogenic differentiation were evaluated in vitro, before and after epigenetic modulation. hMeDIP-seq identified 99 peaks of hyper-hydroxymethylation and 150 peaks of hypo-hydroxymethylation in nuclear-encoded mitochondrial genes from Obese- versus Non-obese-MSCs. Integrated hMeDIP-seq/mRNA-seq analysis identified a select group of overlapping (altered levels of both 5hmC and mRNA) nuclear-encoded mitochondrial genes involved in ATP production, redox activity, cell proliferation, migration, fatty acid metabolism, and neuronal development. Furthermore, Obese-MSCs exhibited decreased mitochondrial matrix density, membrane potential, and levels of fatty acid metabolites, increased superoxide production, and impaired neuronal differentiation, which improved with epigenetic modulation. Obesity elicits epigenetic changes in mitochondria-related genes in human adipose-derived MSCs, accompanied by structural and functional changes in their mitochondria and impaired fatty acid metabolism and neurogenic differentiation capacity. These observations may assist in developing novel therapies to preserve the potential of MSCs for tissue repair and regeneration in obese individuals.

Long-Read Sequencing Genome Assembly of Ceratocystis fimbriata Enables Development of Molecular Diagnostics for Sweetpotato Black Rot

Ceratocystis fimbriata is a destructive fungal pathogen of sweetpotato (Ipomoea batatas) that leads to losses at all stages of sweetpotato production. Accurate detection of C. fimbriata would allow for more efficient deployment of management tactics in sweetpotato production. To develop a diagnostic assay, a hybrid genome assembly of C. fimbriata isolate AS236 was generated. The resulting 31.7-MB assembly was near-chromosome level, with 18 contigs, 6,481 predicted genes, and a BUSCO completion score of 98.4% when compared with the fungus-specific lineage database. Additional Illumina DNA reads from C. manginecans, C. platani, and a second C. fimbriata isolate (C1421) were then mapped to the assembled genome using BOWTIE2 and counted using HTSeq, which identified 148 genes present only within C. fimbriata as molecular diagnostic candidates; 6 single-copy and 35 highly multi-copy (>40 BLAST hits), as determined through a self-BLAST-P alignment. Primers for PCR were designed in the 200-bp flanking region of the first exon for each candidate, and the candidates were validated against a diverse DNA panel containing Ceratocystis species, sweetpotato pathogens, and plants. After validation, two diagnostic candidates amplified only C. fimbriata DNA and were considered to be highly specific to the species. These genetic markers will serve as valuable diagnostic tools with multiple applications including the detection of C. fimbriata in seed, soil, and wash water in sweetpotato production.

The molecular evolution of cancer associated genes in mammals

Cancer is a disease that many multicellular organisms have faced for millions of years, and species have evolved various tumour suppression mechanisms to control oncogenesis. Although cancer occurs across the tree of life, cancer related mortality risks vary across mammalian orders, with Carnivorans particularly affected. Evolutionary theory predicts different selection pressures on genes associated with cancer progression and suppression, including oncogenes, tumour suppressor genes and immune genes. Therefore, we investigated the evolutionary history of cancer associated gene sequences across 384 mammalian taxa, to detect signatures of selection across categories of oncogenes (GRB2, FGL2 and CDC42), tumour suppressors (LITAF, Casp8 and BRCA2) and immune genes (IL2, CD274 and B2M). This approach allowed us to conduct a fine scale analysis of gene wide and site-specific signatures of selection across mammalian lineages under the lens of cancer susceptibility. Phylogenetic analyses revealed that for most species the evolution of cancer associated genes follows the species' evolution. The gene wide selection analyses revealed oncogenes being the most conserved, tumour suppressor and immune genes having similar amounts of episodic diversifying selection. Despite BRCA2's status as a key caretaker gene, episodic diversifying selection was detected across mammals. The site-specific selection analyses revealed that the two apoptosis associated domains of the Casp8 gene of bats (Chiroptera) are under opposing forces of selection (positive and negative respectively), highlighting the importance of site-specific selection analyses to understand the evolution of highly complex gene families. Our results highlighted the need to critically assess different types of selection pressure on cancer associated genes when investigating evolutionary adaptations to cancer across the tree of life. This study provides an extensive assessment of cancer associated genes in mammals with highly representative, and substantially large sample size for a comparative genomic analysis in the field and identifies various avenues for future research into the mechanisms of cancer resistance and susceptibility in mammals.

Determinants of gastric cancer immune escape identified from non-coding immune-landscape quantitative trait loci

The landscape of non-coding mutations in cancer progression and immune evasion is largely unexplored. Here, we identify transcrptome-wide somatic and germline 3' untranslated region (3'-UTR) variants from 375 gastric cancer patients from The Cancer Genome Atlas. By performing gene expression quantitative trait loci (eQTL) and immune landscape QTL (ilQTL) analysis, we discover 3'-UTR variants with cis effects on expression and immune landscape phenotypes, such as immune cell infiltration and T cell receptor diversity. Using a massively parallel reporter assay, we distinguish between causal and correlative effects of 3'-UTR eQTLs in immune-related genes. Our approach identifies numerous 3'-UTR eQTLs and ilQTLs, providing a unique resource for the identification of immunotherapeutic targets and biomarkers. A prioritized ilQTL variant signature predicts response to immunotherapy better than standard-of-care PD-L1 expression in independent patient cohorts, showcasing the untapped potential of non-coding mutations in cancer.

Conserved regulatory switches for the transition from natal down to juvenile feather in birds

The transition from natal downs for heat conservation to juvenile feathers for simple flight is a remarkable environmental adaptation process in avian evolution. However, the underlying epigenetic mechanism for this primary feather transition is mostly unknown. Here we conducted time-ordered gene co-expression network construction, epigenetic analysis, and functional perturbations in developing feather follicles to elucidate four downy-juvenile feather transition events. We report that extracellular matrix reorganization leads to peripheral pulp formation, which mediates epithelial-mesenchymal interactions for branching morphogenesis. α-SMA (ACTA2) compartmentalizes dermal papilla stem cells for feather renewal cycling. LEF1 works as a key hub of Wnt signaling to build rachis and converts radial downy to bilateral symmetry. Novel usage of scale keratins strengthens feather sheath with SOX14 as the epigenetic regulator. We show that this primary feather transition is largely conserved in chicken (precocial) and zebra finch (altricial) and discuss the possibility that this evolutionary adaptation process started in feathered dinosaurs.

Cell-type-specific mRNA transcription and degradation kinetics in zebrafish embryogenesis from metabolically labeled single-cell RNA-seq

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the relative contributions of mRNA transcription and degradation to shaping those profiles remains challenging, especially within embryos with diverse cellular identities. Here, we combine single-cell RNA-Seq and metabolic labeling to capture temporal cellular transcriptomes of zebrafish embryos where newly-transcribed (zygotic) and pre-existing (maternal) mRNA can be distinguished. We introduce kinetic models to quantify mRNA transcription and degradation rates within individual cell types during their specification. These models reveal highly varied regulatory rates across thousands of genes, coordinated transcription and destruction rates for many transcripts, and link differences in degradation to specific sequence elements. They also identify cell-type-specific differences in degradation, namely selective retention of maternal transcripts within primordial germ cells and enveloping layer cells, two of the earliest specified cell types. Our study provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

Bento: a toolkit for subcellular analysis of spatial transcriptomics data

The spatial organization of molecules in a cell is essential for their functions. While current methods focus on discerning tissue architecture, cell-cell interactions, and spatial expression patterns, they are limited to the multicellular scale. We present Bento, a Python toolkit that takes advantage of single-molecule information to enable spatial analysis at the subcellular scale. Bento ingests molecular coordinates and segmentation boundaries to perform three analyses: defining subcellular domains, annotating localization patterns, and quantifying gene-gene colocalization. We demonstrate MERFISH, seqFISH + , Molecular Cartography, and Xenium datasets. Bento is part of the open-source Scverse ecosystem, enabling integration with other single-cell analysis tools.

Cancer-associated fibroblasts rewire the estrogen receptor response in luminal breast cancer, enabling estrogen independence

Advanced breast cancers represent a major therapeutic challenge due to their refractoriness to treatment. Cancer-associated fibroblasts (CAFs) are the most abundant constituents of the tumor microenvironment and have been linked to most hallmarks of cancer. However, the influence of CAFs on therapeutic outcome remains largely unchartered. Here, we reveal that spatial coincidence of abundant CAF infiltration with malignant cells was associated with reduced estrogen receptor (ER)-α expression and activity in luminal breast tumors. Notably, CAFs mediated estrogen-independent tumor growth by selectively regulating ER-α signaling. Whereas most prototypical estrogen-responsive genes were suppressed, CAFs maintained gene expression related to therapeutic resistance, basal-like differentiation, and invasion. A functional drug screen in co-cultures identified effector pathways involved in the CAF-induced regulation of ER-α signaling. Among these, the Transforming Growth Factor-β and the Janus kinase signaling cascades were validated as actionable targets to counteract the CAF-induced modulation of ER-α activity. Finally, genes that were downregulated in cancer cells by CAFs were predictive of poor response to endocrine treatment. In conclusion, our work reveals that CAFs directly control the luminal breast cancer phenotype by selectively modulating ER-α expression and transcriptional function, and further proposes novel targets to disrupt the crosstalk between CAFs and tumor cells to reinstate treatment response to endocrine therapy in patients.

H3K27ac modification and transcription characteristics of adipose and muscle tissues in Chuxiang Black pig

The establishment of high-quality pork breeds for improving meat quality in the pig industry is needed. The Chuxiang Black (CX) pig is a new breed developed from Chinese local pigs and Western lean pigs that has a high proportion of lean meat and excellent meat quality. However, the characteristics of cis-regulatory elements in CX pigs are still unknown. In this study, cis-regulatory elements of muscle and adipose tissues in CX pigs were investigated using ChIP-seq and RNA sequencing. Compared with the reported cis-regulatory elements of muscle and adipose tissues, 1768 and 1012 highly activated enhancers and 433 and 275 highly activated promoters in CX muscle and adipose tissues were identified, respectively. Motif analysis showed that transcription factors, such as MEF2A and MEF2C, were core regulators of highly activated enhancers and promoters in muscle. Similarly, the transcription factors JUNB and CUX1 were identified as essential for highly activated enhancers and promoters in CX adipose tissue. These results enrich the resources for the analysis of cis-regulatory elements in the pig genome and provide new basic data for further meat quality improvement through breeding in pigs.

Quiescence enables unrestricted cell fate in naive embryonic stem cells

Quiescence in stem cells is traditionally considered as a state of inactive dormancy or with poised potential. Naive mouse embryonic stem cells (ESCs) can enter quiescence spontaneously or upon inhibition of MYC or fatty acid oxidation, mimicking embryonic diapause in vivo. The molecular underpinning and developmental potential of quiescent ESCs (qESCs) are relatively unexplored. Here we show that qESCs possess an expanded or unrestricted cell fate, capable of generating both embryonic and extraembryonic cell types (e.g., trophoblast stem cells). These cells have a divergent metabolic landscape comparing to the cycling ESCs, with a notable decrease of the one-carbon metabolite S-adenosylmethionine. The metabolic changes are accompanied by a global reduction of H3K27me3, an increase of chromatin accessibility, as well as the de-repression of endogenous retrovirus MERVL and trophoblast master regulators. Depletion of methionine adenosyltransferase Mat2a or deletion of Eed in the polycomb repressive complex 2 results in removal of the developmental constraints towards the extraembryonic lineages. Our findings suggest that quiescent ESCs are not dormant but rather undergo an active transition towards an unrestricted cell fate.

Epigenetic regulatory layers in the 3D nucleus

Nearly 7 decades have elapsed since Francis Crick introduced the central dogma of molecular biology, as part of his ideas on protein synthesis, setting the fundamental rules of sequence information transfer from DNA to RNAs and proteins. We have since learned that gene expression is finely tuned in time and space, due to the activities of RNAs and proteins on regulatory DNA elements, and through cell-type-specific three-dimensional conformations of the genome. Here, we review major advances in genome biology and discuss a set of ideas on gene regulation and highlight how various biomolecular assemblies lead to the formation of structural and regulatory features within the nucleus, with roles in transcriptional control. We conclude by suggesting further developments that will help capture the complex, dynamic, and often spatially restricted events that govern gene expression in mammalian cells.

Protein interaction networks in the vasculature prioritize genes and pathways underlying coronary artery disease

Population-based association studies have identified many genetic risk loci for coronary artery disease (CAD), but it is often unclear how genes within these loci are linked to CAD. Here, we perform interaction proteomics for 11 CAD-risk genes to map their protein-protein interactions (PPIs) in human vascular cells and elucidate their roles in CAD. The resulting PPI networks contain interactions that are outside of known biology in the vasculature and are enriched for genes involved in immunity-related and arterial-wall-specific mechanisms. Several PPI networks derived from smooth muscle cells are significantly enriched for genetic variants associated with CAD and related vascular phenotypes. Furthermore, the networks identify 61 genes that are found in genetic loci associated with risk of CAD, prioritizing them as the causal candidates within these loci. These findings indicate that the PPI networks we have generated are a rich resource for guiding future research into the molecular pathogenesis of CAD.

Evolution of lysine-specific demethylase 1 and REST corepressor gene families and their molecular interaction

Lysine-specific demethylase 1A (LSD1) binds to the REST corepressor (RCOR) protein family of corepressors to erase transcriptionally active marks on histones. Functional diversity in these complexes depends on the type of RCOR included, which modulates the catalytic activity of the complex. Here, we studied the duplicative history of the RCOR and LSD gene families and analyzed the evolution of their interaction. We found that RCOR genes are the product of the two rounds of whole-genome duplications that occurred early in vertebrate evolution. In contrast, the origin of the LSD genes traces back before to the divergence of animals and plants. Using bioinformatics tools, we show that the RCOR and LSD1 interaction precedes the RCOR repertoire expansion that occurred in the last common ancestor of jawed vertebrates. Overall, we trace LSD1-RCOR complex evolution and propose that animal non-model species offer advantages in addressing questions about the molecular biology of this epigenetic complex.

Hematopoietic stem cells undergo a lymphoid to myeloid switch in early stages of emergency granulopoiesis

Emergency granulopoiesis is the enhanced and accelerated production of granulocytes that occurs during acute infection. The contribution of hematopoietic stem cells (HSCs) to this process was reported; however, how HSCs participate in emergency granulopoiesis remains elusive. Here, using a mouse model of emergency granulopoiesis we observe transcriptional changes in HSCs as early as 4 h after lipopolysaccharide (LPS) administration. We observe that the HSC identity is changed towards a myeloid-biased HSC and show that CD201 is enriched in lymphoid-biased HSCs. While CD201 expression under steady-state conditions reveals a lymphoid bias, under emergency granulopoiesis loss of CD201 marks the lymphoid-to-myeloid transcriptional switch. Mechanistically, we determine that lymphoid-biased CD201+ HSCs act as a first response during emergency granulopoiesis due to direct sensing of LPS by TLR4 and downstream activation of NF-κΒ signaling. The myeloid-biased CD201- HSC population responds indirectly during an acute infection by sensing G-CSF, increasing STAT3 phosphorylation, and upregulating LAP/LAP* C/EBPβ isoforms. In conclusion, HSC subpopulations support early phases of emergency granulopoiesis due to their transcriptional rewiring from a lymphoid-biased to myeloid-biased population and thus establishing alternative paths to supply elevated numbers of granulocytes.

Testing associations between human anxiety and genes previously implicated by mouse anxiety models

Anxiety disorders are common and can be debilitating, with effective treatments remaining hampered by an incomplete understanding of the underlying genetic etiology. Improvements have been made in understanding the genetic influences on mouse behavioral models of anxiety, yet it is unclear the extent to which genes identified in these experimental systems contribute to genetic variation in human anxiety phenotypes. Leveraging new and existing large-scale human genome-wide association studies, we tested whether sets of genes previously identified in mouse anxiety-like behavior studies contribute to a range of human anxiety disorders. When tested as individual genes, 13 mouse-identified genes were associated with human anxiety phenotypes, suggesting an overlap of individual genes contributing to both mouse models of anxiety-like behaviors and human anxiety traits. When genes were tested as sets, we did identify 14 significant associations between mouse gene sets and human anxiety, but the majority of gene sets showed no significant association with human anxiety phenotypes. These few significant associations indicate a need to identify and develop more translatable mouse models by identifying sets of genes that "match" between model systems and specific human phenotypes of interest. We suggest that continuing to develop improved behavioral paradigms and finer-scale experimental data, for instance from individual neuronal subtypes or cell-type-specific expression data, is likely to improve our understanding of the genetic etiology and underlying functional changes in anxiety disorders.

Cellular reprogramming is driven by widespread rewiring of promoter-enhancer interactions

BACKGROUND

Long-range interactions between promoters and cis-regulatory elements, such as enhancers, play critical roles in gene regulation. However, the role of three-dimensional (3D) chromatin structure in orchestrating changes in transcriptional regulation during direct cell reprogramming is not fully understood.

RESULTS

Here, we performed integrated analyses of chromosomal architecture, epigenetics, and gene expression using Hi-C, promoter Capture Hi-C (PCHi-C), ChIP-seq, and RNA-seq during trans-differentiation of Pre-B cells into macrophages with a β-estradiol inducible C/EBPαER transgene. Within 1h of β-estradiol induction, C/EBPα translocated from the cytoplasm to the nucleus, binding to thousands of promoters and putative regulatory elements, resulting in the downregulation of Pre-B cell-specific genes and induction of macrophage-specific genes. Hi-C results were remarkably consistent throughout trans-differentiation, revealing only a small number of TAD boundary location changes, and A/B compartment switches despite significant changes in the expression of thousands of genes. PCHi-C revealed widespread changes in promoter-anchored loops with decreased interactions in parallel with decreased gene expression, and new and increased promoter-anchored interactions in parallel with increased expression of macrophage-specific genes.

CONCLUSIONS

Overall, our data demonstrate that C/EBPα-induced trans-differentiation involves few changes in genome architecture at the level of TADs and A/B compartments, in contrast with widespread reorganization of thousands of promoter-anchored loops in association with changes in gene expression and cell identity.

Polycomb protein binding and looping mediated by Polycomb Response Elements in the ON transcriptional state

Polycomb group proteins (PcG) mediate epigenetic silencing of important developmental genes and other targets. In Drosophila, canonical PcG-target genes contain Polycomb Response Elements (PREs) that recruit PcG protein complexes including PRC2 that trimethylates H3K27 forming large H3K27me3 domains. In the OFF transcriptional state, PREs loop with each other and this looping strengthens silencing. Here we address the question of what PcG proteins bind to PREs when canonical PcG target genes are expressed, and whether PREs loop when these genes are ON. Our data show that the answer to this question is PRE-specific but general conclusions can be made. First, within a PcG-target gene, some regulatory DNA can remain covered with H3K27me3 and PcG proteins remain bound to PREs in these regions. Second, when PREs are within H3K27ac domains, PcG-binding decreases, however, this depends on the protein and PRE. The DNA binding protein GAF, and the PcG protein Ph remain at PREs even when other PcG proteins are greatly depleted. In the ON state, PREs can still loop with each other, but also form loops with presumptive enhancers. These data support the model that, in addition to their role in PcG silencing, PREs can act as "promoter-tethering elements" mediating interactions between promoter proximal PREs and distant enhancers.

Single-cell transcriptome analysis of the germ cells and somatic cells during mitotic quiescence stage in goats

The mitotic quiescence of prospermatogonia is the event known to occur during genesis of the male germline and is tied to the development of the spermatogenic lineage. The regulatory mechanisms and the functional importance of this process have been demonstrated in mice; however, regulation of this process in human and domestic animal is still largely unknown. In this study, we employed single-cell RNA sequencing to identify transcriptional signatures of prospermatogonia and major somatic cell types in testes of goats at E85, E105, and E125. We identified both common and specific Gene Ontology categories, transcription factor regulatory networks, and cell-cell interactions in cell types from goat testis. We also analyzed the transcriptional dynamic changes in prospermatogonia, Sertoli cells, Leydig cells, and interstitial cells. Our datasets provide a useful resource for the study of domestic animal germline development.

scDiffCom: a tool for differential analysis of cell-cell interactions provides a mouse atlas of aging changes in intercellular communication

Dysregulation of intercellular communication is a hallmark of aging. To better quantify and explore changes in intercellular communication, we present scDiffCom and scAgeCom. scDiffCom is an R package, relying on approximately 5,000 curated ligand-receptor interactions, that performs differential intercellular communication analysis between two conditions from single-cell transcriptomics data. Built upon scDiffCom, scAgeCom is an atlas of age-related cell-cell communication changes covering 23 mouse tissues from 58 single-cell RNA sequencing datasets from Tabula Muris Senis and the Calico murine aging cell atlas. It offers a comprehensive resource of tissue-specific and sex-specific aging dysregulations and highlights age-related intercellular communication changes widespread across the whole body, such as the upregulation of immune system processes and inflammation, the downregulation of developmental processes, angiogenesis and extracellular matrix organization and the deregulation of lipid metabolism. Our analysis emphasizes the relevance of the specific ligands, receptors and cell types regulating these processes. The atlas is available online ( https://scagecom.org ).

Unique repetitive nucleic acid structures mirror switch regions in the human IgH locus

Immunoglobulin (Ig) genes carry the unique ability to be reshaped in peripheral B lymphocytes after these cells encounter a specific antigen. B cells can then further improve their affinity, acquire new functions as memory cells and eventually end up as antibody-secreting cells. Ig class switching is an important change that occurs in this context, thanks to local DNA lesions initiated by the enzyme activation-induced deaminase (AID). Several cis-acting elements of the Ig heavy (IgH) chain locus make it accessible to the AID-mediated lesions that promote class switch recombination (CSR). DNA repeats, with a non-template strand rich in G-quadruplexes (G4)-DNA, are prominent cis-targets of AID and define the so-called "switch" (S) regions specifically targeted for CSR. By analyzing the structure of the human IgH locus, we uncover that abundant DNA repeats, some with a putative G4-rich template strand, are additionally present in downstream portions of the IgH coding genes. These like-S (LS) regions stand as 3' mirror-images of S regions and also show analogies to some previously reported repeats associated with the IgH locus 3' super-enhancer. A regulatory role of LS repeats is strongly suggested by their specific localization close to exons encoding the membrane form of Ig molecules, and by their conservation during mammalian evolution.

Immune checkpoint inhibitor-induced colitis is mediated by polyfunctional lymphocytes and is dependent on an IL23/IFNγ axis

Immune checkpoint inhibitors (CPIs) are a relatively newly licenced cancer treatment, which make a once previously untreatable disease now amenable to a potential cure. Combination regimens of anti-CTLA4 and anti-PD-1 show enhanced efficacy but are prone to off-target immune-mediated tissue injury, particularly at the barrier surfaces. To probe the impact of immune checkpoints on intestinal homoeostasis, mice are challenged with anti-CTLA4 and anti-PD-1 immunotherapy and manipulation of the intestinal microbiota. The immune profile of the colon of these mice with CPI-colitis is analysed using bulk RNA sequencing, single-cell RNA sequencing and flow cytometry. CPI-colitis in mice is dependent on the composition of the intestinal microbiota and by the induction of lymphocytes expressing interferon-γ (IFNγ), cytotoxicity molecules and other pro-inflammatory cytokines/chemokines. This pre-clinical model of CPI-colitis could be attenuated following blockade of the IL23/IFNγ axis. Therapeutic targeting of IFNγ-producing lymphocytes or regulatory networks, may hold the key to reversing CPI-colitis.

TCMBank: bridges between the largest herbal medicines, chemical ingredients, target proteins, and associated diseases with intelligence text mining

Traditional Chinese Medicine (TCM) has long been viewed as a precious source of modern drug discovery. AI-assisted drug discovery (AIDD) has been investigated extensively. However, there are still two challenges in applying AIDD to guide TCM drug discovery: the lack of a large amount of standardized TCM-related information and AIDD is prone to pathological failures in out-of-domain data. We have released TCM Database@Taiwan in 2011, and it has been widely disseminated and used. Now, we developed TCMBank, the largest systematic free TCM database, which is an extension of TCM Database@Taiwan. TCMBank contains 9192 herbs, 61 966 ingredients (unduplicated), 15 179 targets, 32 529 diseases, and their pairwise relationships. By integrating multiple data sources, TCMBank provides 3D structure information of ingredients and provides a standard list and detailed information on herbs, ingredients, targets and diseases. TCMBank has an intelligent document identification module that continuously adds TCM-related information retrieved from the literature in PubChem. In addition, driven by TCMBank big data, we developed an ensemble learning-based drug discovery protocol for identifying potential leads and drug repurposing. We take colorectal cancer and Alzheimer's disease as examples to demonstrate how to accelerate drug discovery by artificial intelligence. Using TCMBank, researchers can view literature-driven relationship mapping between herbs/ingredients and genes/diseases, allowing the understanding of molecular action mechanisms for ingredients and identification of new potentially effective treatments. TCMBank is available at https://TCMBank.CN/.

Insights into the genetics of menopausal vasomotor symptoms: genome-wide analyses of routinely-collected primary care health records

BACKGROUND

Vasomotor symptoms (VMS) can often significantly impact women's quality of life at menopause. In vivo studies have shown that increased neurokinin B (NKB) / neurokinin 3 receptor (NK3R) signalling contributes to VMS, with previous genetic studies implicating the TACR3 gene locus that encodes NK3R. Large-scale genomic analyses offer the possibility of biological insights but few such studies have collected data on VMS, while proxy phenotypes such as hormone replacement therapy (HRT) use are likely to be affected by changes in clinical practice. We investigated the genetic basis of VMS by analysing routinely-collected health records.

METHODS

We performed a GWAS of VMS derived from linked primary-care records and cross-sectional self-reported HRT use in up to 153,152 women from UK Biobank, a population-based cohort. In a subset of this cohort (n = 39,356), we analysed exome-sequencing data to test the association with VMS of rare deleterious genetic variants. Finally, we used Mendelian randomisation analysis to investigate the reasons for HRT use over time.

RESULTS

Our GWAS of health-records derived VMS identified a genetic signal near TACR3 associated with a lower risk of VMS (OR=0.76 (95% CI 0.72,0.80) per A allele, P=3.7x10-27), which was consistent with previous studies, validating this approach. Conditional analyses demonstrated independence of genetic signals for puberty timing and VMS at the TACR3 locus, including a rare variant predicted to reduce functional NK3R levels that was associated with later menarche (P = 5 × 10-9) but showed no association with VMS (P = 0.6). Younger menopause age was causally-associated with greater HRT use before 2002 but not after.

CONCLUSIONS

We provide support for TACR3 in the genetic basis of VMS but unexpectedly find that rare genomic variants predicted to lower NK3R levels did not modify VMS, despite the proven efficacy of NK3R antagonists. Using genomics we demonstrate changes in genetic associations with HRT use over time, arising from a change in clinical practice since the early 2000s, which is likely to reflect a switch from preventing post-menopausal complications in women with earlier menopause to primarily treating VMS. Our study demonstrates that integrating routinely-collected primary care health records and genomic data offers great potential for exploring the genetic basis of symptoms.

A lamprey neural cell type atlas illuminates the origins of the vertebrate brain

The vertebrate brain emerged more than ~500 million years ago in common evolutionary ancestors. To systematically trace its cellular and molecular origins, we established a spatially resolved cell type atlas of the entire brain of the sea lamprey-a jawless species whose phylogenetic position affords the reconstruction of ancestral vertebrate traits-based on extensive single-cell RNA-seq and in situ sequencing data. Comparisons of this atlas to neural data from the mouse and other jawed vertebrates unveiled various shared features that enabled the reconstruction of cell types, tissue structures and gene expression programs of the ancestral vertebrate brain. However, our analyses also revealed key tissues and cell types that arose later in evolution. For example, the ancestral brain was probably devoid of cerebellar cell types and oligodendrocytes (myelinating cells); our data suggest that the latter emerged from astrocyte-like evolutionary precursors in the jawed vertebrate lineage. Altogether, our work illuminates the cellular and molecular architecture of the ancestral vertebrate brain and provides a foundation for exploring its diversification during evolution.

PTEN deficiency exposes a requirement for an ARF GTPase module for integrin-dependent invasion in ovarian cancer

Dysregulation of the PI3K/AKT pathway is a common occurrence in high-grade serous ovarian carcinoma (HGSOC), with the loss of the tumour suppressor PTEN in HGSOC being associated with poor prognosis. The cellular mechanisms of how PTEN loss contributes to HGSOC are largely unknown. We here utilise time-lapse imaging of HGSOC spheroids coupled to a machine learning approach to classify the phenotype of PTEN loss. PTEN deficiency induces PI(3,4,5)P3 -rich and -dependent membrane protrusions into the extracellular matrix (ECM), resulting in a collective invasion phenotype. We identify the small GTPase ARF6 as a crucial vulnerability of HGSOC cells upon PTEN loss. Through a functional proteomic CRISPR screen of ARF6 interactors, we identify the ARF GTPase-activating protein (GAP) AGAP1 and the ECM receptor β1-integrin (ITGB1) as key ARF6 interactors in HGSOC regulating PTEN loss-associated invasion. ARF6 functions to promote invasion by controlling the recycling of internalised, active β1-integrin to maintain invasive activity into the ECM. The expression of the CYTH2-ARF6-AGAP1 complex in HGSOC patients is inversely associated with outcome, allowing the identification of patient groups with improved versus poor outcome. ARF6 may represent a therapeutic vulnerability in PTEN-depleted HGSOC.

An extended transcription factor regulatory network controls hepatocyte identity

Cell identity is specified by a core transcriptional regulatory circuitry (CoRC), typically limited to a small set of interconnected cell-specific transcription factors (TFs). By mining global hepatic TF regulons, we reveal a more complex organization of the transcriptional regulatory network controlling hepatocyte identity. We show that tight functional interconnections controlling hepatocyte identity extend to non-cell-specific TFs beyond the CoRC, which we call hepatocyte identity (Hep-ID)CONNECT TFs. Besides controlling identity effector genes, Hep-IDCONNECT TFs also engage in reciprocal transcriptional regulation with TFs of the CoRC. In homeostatic basal conditions, this translates into Hep-IDCONNECT TFs being involved in fine tuning CoRC TF expression including their rhythmic expression patterns. Moreover, a role for Hep-IDCONNECT TFs in the control of hepatocyte identity is revealed in dedifferentiated hepatocytes where Hep-IDCONNECT TFs are able to reset CoRC TF expression. This is observed upon activation of NR1H3 or THRB in hepatocarcinoma or in hepatocytes subjected to inflammation-induced loss of identity. Our study establishes that hepatocyte identity is controlled by an extended array of TFs beyond the CoRC.

Heuristic-enabled active machine learning: A case study of predicting essential developmental stage and immune response genes in Drosophila melanogaster

Computational prediction of absolute essential genes using machine learning has gained wide attention in recent years. However, essential genes are mostly conditional and not absolute. Experimental techniques provide a reliable approach of identifying conditionally essential genes; however, experimental methods are laborious, time and resource consuming, hence computational techniques have been used to complement the experimental methods. Computational techniques such as supervised machine learning, or flux balance analysis are grossly limited due to the unavailability of required data for training the model or simulating the conditions for gene essentiality. This study developed a heuristic-enabled active machine learning method based on a light gradient boosting model to predict essential immune response and embryonic developmental genes in Drosophila melanogaster. We proposed a new sampling selection technique and introduced a heuristic function which replaces the human component in traditional active learning models. The heuristic function dynamically selects the unlabelled samples to improve the performance of the classifier in the next iteration. Testing the proposed model with four benchmark datasets, the proposed model showed superior performance when compared to traditional active learning models (random sampling and uncertainty sampling). Applying the model to identify conditionally essential genes, four novel essential immune response genes and a list of 48 novel genes that are essential in embryonic developmental condition were identified. We performed functional enrichment analysis of the predicted genes to elucidate their biological processes and the result evidence our predictions. Immune response and embryonic development related processes were significantly enriched in the essential immune response and embryonic developmental genes, respectively. Finally, we propose the predicted essential genes for future experimental studies and use of the developed tool accessible at http://heal.covenantuniversity.edu.ng for conditional essentiality predictions.

A role for SETD2 loss in tumorigenesis through DNA methylation dysregulation

SETD2-dependent H3 Lysine-36 trimethylation (H3K36me3) has been recently linked to the deposition of de-novo DNA methylation. SETD2 is frequently mutated in cancer, however, the functional impact of SETD2 loss and depletion on DNA methylation across cancer types and tumorigenesis is currently unknown. Here, we perform a pan-cancer analysis and show that both SETD2 mutation and reduced expression are associated with DNA methylation dysregulation across 21 out of the 24 cancer types tested. In renal cancer, these DNA methylation changes are associated with altered gene expression of oncogenes, tumour suppressors, and genes involved in neoplasm invasiveness, including TP53, FOXO1, and CDK4. This suggests a new role for SETD2 loss in tumorigenesis and cancer aggressiveness through DNA methylation dysregulation. Moreover, using a robust machine learning methodology, we develop and validate a 3-CpG methylation signature which is sufficient to predict SETD2 mutation status with high accuracy and correlates with patient prognosis.

Comprehensive mRNA-sequencing-based characterization of three HEK-293 cell lines during an rAAV production process for gene therapy applications

Human embryonal kidney cells (HEK-293) are the most common host cells used for transient recombinant adeno-associated virus (rAAV) production in pharmaceutical industry. To better cover the expected gene therapy product demands in the future, different traditional strategies such as cell line sub-cloning and/or addition of chemical substances to the fermentation media have been used to maximize titers and improve product quality. A more effective and advanced approach to boost yield can be envisaged by characterizing the transcriptome of different HEK-293 cell line pedigrees with distinct rAAV productivity patterns to subsequently identify potential gene targets for cell engineering. In this work, the mRNA expression profile of three HEK-293 cell lines, resulting in various yields during a fermentation batch process for rAAV production, was investigated to gain basic insight into cell variability and eventually to identify genes that correlate with productivity. Mock runs using only transfection reagents were performed in parallel as a control. It finds significant differences in gene regulatory behaviors between the three cell lines at different growth and production stages. The evaluation of these transcriptomics profiles combined with collected in-process control parameters and titers shed some light on potential cell engineering targets to maximize transient production of rAAV in HEK-293 cells.

ORC1 binds to cis-transcribed RNAs for efficient activation of replication origins

Cells must coordinate the activation of thousands of replication origins dispersed throughout their genome. Active transcription is known to favor the formation of mammalian origins, although the role that RNA plays in this process remains unclear. We show that the ORC1 subunit of the human Origin Recognition Complex interacts with RNAs transcribed from genes with origins in their transcription start sites (TSSs), displaying a positive correlation between RNA binding and origin activity. RNA depletion, or the use of ORC1 RNA-binding mutant, result in inefficient activation of proximal origins, linked to impaired ORC1 chromatin release. ORC1 RNA binding activity resides in its intrinsically disordered region, involved in intra- and inter-molecular interactions, regulation by phosphorylation, and phase-separation. We show that RNA binding favors ORC1 chromatin release, by regulating its phosphorylation and subsequent degradation. Our results unveil a non-coding function of RNA as a dynamic component of the chromatin, orchestrating the activation of replication origins.

Regulation of thrombin activity by ligand-induced topological alteration in a thrombin-binding aptamer

Thrombin-binding aptamer (TBA), which forms a G-quadruplex (G4) structure with anti-parallel topology, interacts with thrombin to inhibit its enzymatic activity. Here we show that the G4-topology-altering ligand L2H2-2M2EA-6LCO (6LCO) changes the anti-parallel topology of TBA G4 to the parallel topology, thereby abrogating the thrombin-inhibitory activity of TBA. This finding suggests that G4 ligands that alter topology may be promising drug candidates for diseases involving G4-binding proteins.

Unveiling the role of G-quadruplex structure in promoter region: Regulation of ABCA1 expression in macrophages possibly via NONO protein recruitment

ABCA1 has been found to be critical for cholesterol efflux in macrophages. Understanding the mechanism regulating ABCA1 expression is important for the prevention and treatment of atherosclerosis. In the present study, a G-quadruplex (G4) structure was identified in the ABCA1 promoter region. This G4 was shown to be essential for ABCA1 transcription. Stabilizing the G4 by ligands surprisingly upregulated ABCA1 expression in macrophages. Knocking out the G4 remarkably reduced ABCA1 expression, and abolished the increase of ABCA1 expression induced by the G4 ligand. By pull-down assays, the protein NONO was identified as an ABCA1 G4 binder. Overexpression or repression of NONO significantly induced upregulation and downregulation of ABCA1 expression, respectively. ChIP and EMSA experiments showed that the G4 ligand promoted the binding between the ABCA1 G4 and NONO, which led to more recruitment of NONO to the promoter region and enhanced ABCA1 transcription. Finally, the G4 ligand was shown to significantly reduce the accumulation of cholesterol in macrophages. This study showed a new insight into the regulation of gene expression by G4, and provided a new molecular mechanism regulating ABCA1 expression in macrophages. Furthermore, the study showed a possible novel application of the G4 ligand: preventing and treating atherosclerosis.

Exosomal miRNA-223-3p derived from tumor associated macrophages promotes pulmonary metastasis of breast cancer 4T1 cells

Research about the effect of exosomes derived from tumor associated macrophages (TAM-exos) in the distant organ metastasis of breast cancer is limited. In this study, we found that TAM-exos could promote the migration of 4T1 cells. Through comparing the expression of microRNAs in 4T1 cells, TAM-exos, and exosomes from bone marrow derived macrophages (BMDM-exos) by sequencing, miR-223-3p and miR-379-5p were screened out as two noteworthy differentially expressed microRNAs. Furthermore, miR-223-3p was confirmed to be the reason for the improved migration and metastasis of 4T1 cells. The expression of miR-223-3p was also increased in 4T1 cells isolated from the lung of tumor-bearing mice. Cbx5, which has been reported to be closely related with metastasis of breast cancer, was identified to be the target of miR-223-3p. Based on the information of breast cancer patients from online databases, miR-223-3p had a negative correlation with the overall survival rate of breast cancer patients within a three-year follow-up, while Cbx5 showed an opposite relationship. Taken together, miR-223-3p in TAM-exos can be delivered into 4T1 cells and exosomal miR-223-3p promotes pulmonary metastasis of 4T1 cells by targeting Cbx5.

South Asian medical cohorts reveal strong founder effects and high rates of homozygosity

The benefits of large-scale genetic studies for healthcare of the populations studied are well documented, but these genetic studies have traditionally ignored people from some parts of the world, such as South Asia. Here we describe whole genome sequence (WGS) data from 4806 individuals recruited from the healthcare delivery systems of Pakistan, India and Bangladesh, combined with WGS from 927 individuals from isolated South Asian populations. We characterize population structure in South Asia and describe a genotyping array (SARGAM) and imputation reference panel that are optimized for South Asian genomes. We find evidence for high rates of reproductive isolation, endogamy and consanguinity that vary across the subcontinent and that lead to levels of rare homozygotes that reach 100 times that seen in outbred populations. Founder effects increase the power to associate functional variants with disease processes and make South Asia a uniquely powerful place for population-scale genetic studies.

DNA based neoepitope vaccination induces tumor control in syngeneic mouse models

Recent findings have positioned tumor mutation-derived neoepitopes as attractive targets for cancer immunotherapy. Cancer vaccines that deliver neoepitopes via various vaccine formulations have demonstrated promising preliminary results in patients and animal models. In the presented work, we assessed the ability of plasmid DNA to confer neoepitope immunogenicity and anti-tumor effect in two murine syngeneic cancer models. We demonstrated that neoepitope DNA vaccination led to anti-tumor immunity in the CT26 and B16F10 tumor models, with the long-lasting presence of neoepitope-specific T-cell responses in blood, spleen, and tumors after immunization. We further observed that engagement of both the CD4+ and CD8+ T cell compartments was essential to hamper tumor growth. Additionally, combination therapy with immune checkpoint inhibition provided an additive effect, superior to either monotherapy. DNA vaccination offers a versatile platform that allows the encoding of multiple neoepitopes in a single formulation and is thus a feasible strategy for personalized immunotherapy via neoepitope vaccination.

G-Quadruplexes in Nuclear Biomolecular Condensates

G-quadruplexes (G4s) have long been implicated in the regulation of chromatin packaging and gene expression. These processes require or are accelerated by the separation of related proteins into liquid condensates on DNA/RNA matrices. While cytoplasmic G4s are acknowledged scaffolds of potentially pathogenic condensates, the possible contribution of G4s to phase transitions in the nucleus has only recently come to light. In this review, we summarize the growing evidence for the G4-dependent assembly of biomolecular condensates at telomeres and transcription initiation sites, as well as nucleoli, speckles, and paraspeckles. The limitations of the underlying assays and the remaining open questions are outlined. We also discuss the molecular basis for the apparent permissive role of G4s in the in vitro condensate assembly based on the interactome data. To highlight the prospects and risks of G4-targeting therapies with respect to the phase transitions, we also touch upon the reported effects of G4-stabilizing small molecules on nuclear biomolecular condensates.

Aberrant splicing prediction across human tissues

Aberrant splicing is a major cause of genetic disorders but its direct detection in transcriptomes is limited to clinically accessible tissues such as skin or body fluids. While DNA-based machine learning models can prioritize rare variants for affecting splicing, their performance in predicting tissue-specific aberrant splicing remains unassessed. Here we generated an aberrant splicing benchmark dataset, spanning over 8.8 million rare variants in 49 human tissues from the Genotype-Tissue Expression (GTEx) dataset. At 20% recall, state-of-the-art DNA-based models achieve maximum 12% precision. By mapping and quantifying tissue-specific splice site usage transcriptome-wide and modeling isoform competition, we increased precision by threefold at the same recall. Integrating RNA-sequencing data of clinically accessible tissues into our model, AbSplice, brought precision to 60%. These results, replicated in two independent cohorts, substantially contribute to noncoding loss-of-function variant identification and to genetic diagnostics design and analytics.

Alternative CDC20 translational isoforms tune mitotic arrest duration

Mitotic defects activate the spindle-assembly checkpoint, which inhibits the anaphase-promoting complex co-activator CDC20 to induce a prolonged cell cycle arrest1,2. Once errors are corrected, the spindle-assembly checkpoint is silenced, allowing anaphase onset to occur. However, in the presence of persistent unresolvable errors, cells can undergo 'mitotic slippage', exiting mitosis into a tetraploid G1 state and escaping the cell death that results from a prolonged arrest. The molecular logic that enables cells to balance these duelling mitotic arrest and slippage behaviours remains unclear. Here we demonstrate that human cells modulate the duration of their mitotic arrest through the presence of conserved, alternative CDC20 translational isoforms. Downstream translation initiation results in a truncated CDC20 isoform that is resistant to spindle-assembly-checkpoint-mediated inhibition and promotes mitotic exit even in the presence of mitotic perturbations. Our study supports a model in which the relative levels of CDC20 translational isoforms control the duration of mitotic arrest. During a prolonged mitotic arrest, new protein synthesis and differential CDC20 isoform turnover create a timer, with mitotic exit occurring once the truncated Met43 isoform achieves sufficient levels. Targeted molecular changes or naturally occurring cancer mutations that alter CDC20 isoform ratios or its translational control modulate mitotic arrest duration and anti-mitotic drug sensitivity, with potential implications for the diagnosis and treatment of human cancers.

Single-cell temporal dynamics reveals the relative contributions of transcription and degradation to cell-type specific gene expression in zebrafish embryos

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the underlying regulation of mRNA transcription and degradation remains a challenge, especially within whole embryos with diverse cellular identities. Here, we collect temporal cellular transcriptomes of zebrafish embryos, and decompose them into their newly-transcribed (zygotic) and pre-existing (maternal) mRNA components by combining single-cell RNA-Seq and metabolic labeling. We introduce kinetic models capable of quantifying regulatory rates of mRNA transcription and degradation within individual cell types during their specification. These reveal different regulatory rates between thousands of genes, and sometimes between cell types, that shape spatio-temporal expression patterns. Transcription drives most cell-type restricted gene expression. However, selective retention of maternal transcripts helps to define the gene expression profiles of germ cells and enveloping layer cells, two of the earliest specified cell-types. Coordination between transcription and degradation restricts expression of maternal-zygotic genes to specific cell types or times, and allows the emergence of spatio-temporal patterns when overall mRNA levels are held relatively constant. Sequence-based analysis links differences in degradation to specific sequence motifs. Our study reveals mRNA transcription and degradation events that control embryonic gene expression, and provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.