Dual-role transcription factors stabilize intermediate expression levels

Precise control of gene expression levels is essential for normal cell functions, yet how they are defined and tightly maintained, particularly at intermediate levels, remains elusive. Here, using a series of newly developed sequencing, imaging, and functional assays, we uncover a class of transcription factors with dual roles as activators and repressors, referred to as condensate-forming level-regulating dual-action transcription factors (TFs). They reduce high expression but increase low expression to achieve stable intermediate levels. Dual-action TFs directly exert activating and repressing functions via condensate-forming domains that compartmentalize core transcriptional unit selectively. Clinically relevant mutations in these domains, which are linked to a range of developmental disorders, impair condensate selectivity and dual-action TF activity. These results collectively address a fundamental question in expression regulation and demonstrate the potential of level-regulating dual-action TFs as powerful effectors for engineering controlled expression levels.

Adeno-to-squamous transition drives resistance to KRAS inhibition in LKB1 mutant lung cancer

KRASG12C inhibitors (adagrasib and sotorasib) have shown clinical promise in targeting KRASG12C-mutated lung cancers; however, most patients eventually develop resistance. In lung patients with adenocarcinoma with KRASG12C and STK11/LKB1 co-mutations, we find an enrichment of the squamous cell carcinoma gene signature in pre-treatment biopsies correlates with a poor response to adagrasib. Studies of Lkb1-deficient KRASG12C and KrasG12D lung cancer mouse models and organoids treated with KRAS inhibitors reveal tumors invoke a lineage plasticity program, adeno-to-squamous transition (AST), that enables resistance to KRAS inhibition. Transcriptomic and epigenomic analyses reveal ΔNp63 drives AST and modulates response to KRAS inhibition. We identify an intermediate high-plastic cell state marked by expression of an AST plasticity signature and Krt6a. Notably, expression of the AST plasticity signature and KRT6A at baseline correlates with poor adagrasib responses. These data indicate the role of AST in KRAS inhibitor resistance and provide predictive biomarkers for KRAS-targeted therapies in lung cancer.

ChIPr: accurate prediction of cohesin-mediated 3D genome organization from 2D chromatin features

The three-dimensional genome organization influences diverse nuclear processes. Here we present Chromatin Interaction Predictor (ChIPr), a suite of regression models based on deep neural networks, random forest, and gradient boosting to predict cohesin-mediated chromatin interaction strength between any two loci in the genome. The predictions of ChIPr correlate well with ChIA-PET data in four cell lines. The standard ChIPr model requires three experimental inputs: ChIP-Seq signals for RAD21, H3K27ac, and H3K27me3 but works well with just RAD21 signal. Integrative analysis reveals novel insights into the role of CTCF motif, its orientation, and CTCF binding on cohesin-mediated chromatin interactions.

Aging-disturbed FUS phase transition impairs hematopoietic stem cells by altering chromatin structure

ABSTRACT

Aged hematopoietic stem cells (HSCs) exhibit compromised reconstitution capacity. The molecular mechanisms behind this phenomenon are not fully understood. Here, we observed that the expression of FUS is increased in aged HSCs, and enforced FUS recapitulates the phenotype of aged HSCs through arginine-glycine-glycine-mediated aberrant FUS phase transition. By using Fus-gfp mice, we observed that FUShigh HSCs exhibit compromised FUS mobility and resemble aged HSCs both functionally and transcriptionally. The percentage of FUShigh HSCs is increased upon physiological aging and replication stress, and FUSlow HSCs of aged mice exhibit youthful function. Mechanistically, FUShigh HSCs exhibit a different global chromatin organization compared with FUSlow HSCs, which is observed in aged HSCs. Many topologically associating domains (TADs) are merged in aged HSCs because of the compromised binding of CCCTC-binding factor with chromatin, which is invoked by aberrant FUS condensates. It is notable that the transcriptional alteration between FUShigh and FUSlow HSCs originates from the merged TADs and is enriched in HSC aging-related genes. Collectively, this study reveals for the first time that aberrant FUS mobility promotes HSC aging by altering chromatin structure.

[Correlation between C-reactive protein to albumin ratio and restenosis after femoral popliteal stenting in patients with lower extremity arteriosclerotic obliterans]

Objective: To investigate the study of the correlation between C-reactive protein to albumin ratio (CAR) and restenosis after stenting in patients with lower extremity atherosclerotic occlusive disease(LEASO). Methods: The clinical data of 95 patients with LEASO admitted to the Department of Vascular Surgery of Nanjing Drum Tower Hospital from June 2020 to December 2022 were retrospectively analyzed. There were 67 males and 28 females,aged (73.1±9.4) years (range:51 to 92 years). The patients were classified into the restenosis group (n=61) and the patency group (n=34) according to the CT angiography results. Independent sample t test,Mann-Whitney U test and χ2 test were used to compare the data between two groups. Risk factors for restenosis after femoropopliteal artery stenting in patients with LEASO were analyzed using multivariate Cox regression. The relationship between preoperative CAR level and restenosis after stent placement was analyzed. Subject operating characteristic(ROC) curves of CAR were plotted to assess the predictive value of CAR for restenosis after stenting,and the results were expressed as area under the curve (AUC). Results: The aortoiliac calcification grade,number of stents,length of stents,C-reactive protein and CAR levels in restenosis group were higher than those in the patency group,and the serum albumin level was lower than that in the patency group(all P<0.05). And the results of multifactorial Cox regression analysis showed that higher pre-procedure CAR level and lower ABI value was an independent risk factor for in-stent restenosis. The AUC of the ROC curve for restenosis was 0.737(95%CI:0.617 to 0.856),the AUC of the ROC curve for 12-month restenosis was 0.709(95%CI:0.602 to 0.815), and the AUC of the ROC curve for 24-month restenosis was 0.702(95%CI:0.594 to 0.811). Conclusion: Higher pre-procedural CAR levels in patients with LEASO is risk factor for in-stent restenosis,and CAR has a predictive value for restenosis after lower extremity arterial stent dilatation and angioplasty.

SPIRAL: integrating and aligning spatially resolved transcriptomics data across different experiments, conditions, and technologies

Properly integrating spatially resolved transcriptomics (SRT) generated from different batches into a unified gene-spatial coordinate system could enable the construction of a comprehensive spatial transcriptome atlas. Here, we propose SPIRAL, consisting of two consecutive modules: SPIRAL-integration, with graph domain adaptation-based data integration, and SPIRAL-alignment, with cluster-aware optimal transport-based coordination alignment. We verify SPIRAL with both synthetic and real SRT datasets. By encoding spatial correlations to gene expressions, SPIRAL-integration surpasses state-of-the-art methods in both batch effect removal and joint spatial domain identification. By aligning spots cluster-wise, SPIRAL-alignment achieves more accurate coordinate alignments than existing methods.

Aging-induced pseudouridine synthase 10 impairs hematopoietic stem cells

Aged hematopoietic stem cells (HSC) exhibit compromised reconstitution capacity and differentiation-bias towards myeloid lineage, however, the molecular mechanism behind it remains not fully understood. In this study, we observed that the expression of pseudouridine (Ψ) synthase 10 is increased in aged hematopoietic stem and progenitor cells (HSPC) and enforced protein of Ψ synthase 10 (PUS10) recapitulates the phenotype of aged HSC, which is not achieved by its Ψ synthase activity. Consistently, we observed no difference of transcribed RNA pseudouridylation profile between young and aged HSPC. No significant alteration of hematopoietic homeostasis and HSC function is observed in young Pus10-/- mice, while aged Pus10-/- mice exhibit mild alteration of hematopoietic homeostasis and HSC function. Moreover, we observed that PUS10 is ubiquitinated by E3 ubiquitin ligase CRL4DCAF1 complex and the increase of PUS10 in aged HSPC is due to aging-declined CRL4DCAF1- mediated ubiquitination degradation signaling. Taken together, this study for the first time evaluated the role of PUS10 in HSC aging and function, and provided a novel insight into HSC rejuvenation and its clinical application.

Systematic comparison of genome information processing and boundary recognition tools used for genomic island detection

Genomic islands are fragments of foreign DNA that are found in bacterial and archaeal genomes, and are typically associated with symbiosis or pathogenesis. While numerous genomic island detection methods have been proposed, there has been limited evaluation of the efficiency of the genome information processing and boundary recognition tools. In this study, we conducted a review of the statistical methods involved in genomic signatures, host signature extraction, informative signature selection, divergence measures, and boundary detection steps in genomic island prediction. We compared the performances of these methods on simulated experiments using alien fragments obtained from both artificial and real genomes. Our results indicate that among the nine genomic signatures evaluated, genomic signature frequency and full probability performed the best. However, their performance declined when normalized to their expectations and variances, such as Z-score and composition vector. Based on our experiments of the E. coli genome, we found that the confidence intervals of the window variances achieved the best performance in the signature extraction of the host, with the best confidence interval being 1.5-2 times the standard error. Ordered kurtosis was most effective in selecting informative signatures from a single genome, without requiring prior knowledge from other datasets. Among the three divergence measures evaluated, the two-sample t-test was the most successful, and a non-overlapping window with a small eye window (size 2) was best suited for identifying compositionally distinct regions. Finally, the maximum of the Markovian Jensen-Shannon divergence score, in terms of GC-content bias, was found to make boundary detection faster while maintaining a similar error rate.

[Clinical analysis of surgical treatment of infection after interventional operation for major iliac artery disease in 6 cases]

Objective: To explore the surgical treatment strategy of stent graft infection after interventional treatment of major iliac artery related diseases. Methods: Retrospective analysis was performed on the clinical data of 6 patients with secondary stent graft infection after interventional treatment for major iliac artery related diseases admitted to the Department of Vascular Surgery,Affiliated Drum Tower Hospital,Medical School of Nanjing University from November 2021 to August 2022.There were 5 males and 1 female,with a mean age of 64 years (range:49 to 79 years).The infection time was 53 to 3 165 days.All the 6 patients received surgical treatment,including 3 patients who underwent anatomic bypass grafting (axillary arterial-femoral artery bypass,femoral arterial-femoral artery bypass) using artificial vessels,and 3 patients who underwent in situ abdominal aorta reconstruction using bovine pericardium.The perioperative situation,postoperative infection and the occurrence of serious adverse events were recorded,and the safety of different treatment methods and materials was evaluated. Results: All patients successfully completed the operation and no death occurred during hospitalization.Intraoperative blood loss was 2 000~5 000 ml,and intraoperative blood transfusion was 1 600 to 5 350 ml.All the patients were followed up for 81 to 395 days after surgery,and the incision healed well,and no reinfection occurred.Postoperative gastrointestinal bleeding occurred in 1 patient,secondary surgery (retroperitoneal hematoma removal) was performed in 1 patient due to postoperative bleeding at the vascular anastomosis,both lower limb amputations were performed in 1 patient due to postoperative lower limb ischemia,and intermittent claudication occurred in 2 patients.All patients were alive at the last follow-up. Conclusion: For patients with aortic stent graft infection,when the infection is not serious and there is enough space to block the proximal and distal aorta,in situ aortic reconstruction is an effective treatment,and different materials can achieve satisfactory results in a short period of time.

Transition cow clusters with distinctive antioxidant ability and their relation to performance and metabolic status in early lactation

Metabolic and oxidative stress have been characterized as risk factors during the transition period from pregnancy to lactation. Although mutual relations between both types of stress have been suggested, they rarely have been studied concomitantly. For this, a total of 99 individual transition dairy cows (117 cases, 18 cows sampled during 2 consecutive lactations) were included in this experiment. Blood samples were taken at -7, 3, 6, 9, and 21 d relative to calving and concentrations of metabolic parameters (glucose, β-hydroxybutyric acid (BHBA), nonesterified fatty acids, insulin, insulin-like growth factor 1, and fructosamine) were determined. In the blood samples of d 21, biochemical profiles related to liver function and parameters related to oxidative status were determined. First, cases were allocated to 2 different BHBA groups (ketotic vs. nonketotic, N:n = 20:33) consisting of animals with an average postpartum BHBA concentration and at least 2 out of 4 postpartum sampling points exceeding 1.2 mmol/L or remaining below 0.8 mmol/L, respectively. Second, oxidative parameters [proportion of oxidized glutathione to total glutathione in red blood cells (%)], activity of glutathione peroxidase, and of superoxide dismutase, concentrations of malondialdehyde and oxygen radical absorbance capacity were used to perform a fuzzy C-means clustering. From this, 2 groups were obtained [i.e., lower antioxidant ability (LAA80%, n = 31) and higher antioxidant ability (HAA80%, n = 19)], with 80% referring to the cutoff value for cluster membership. Increased concentrations of malondialdehyde, decreased superoxide dismutase activity, and impaired oxygen radical absorbance capacity were observed in the ketotic group compared with the nonketotic group, and inversely, the LAA80% group showed increased concentrations of BHBA. In addition, the concentration of aspartate transaminase was higher in the LAA80% group compared with the HAA80% group. Both the ketotic and LAA80% groups showed lower dry matter intake. However, a lower milk yield was observed in the LAA80% group but not in the ketotic group. Only 1 out of 19 (5.3%) and 3 out of 31 (9.7%) cases from the HAA80% and LAA80% clusters belong to the ketotic and nonketotic group, respectively. These findings suggested that dairy cows vary in oxidative status at the beginning of the lactation, and fuzzy C-means clustering allows to classify observations with distinctive oxidative status. Dairy cows with higher antioxidant capacity in early lactation rarely develop ketosis.

Chromatin-associated RNA Dictates the ecDNA Interactome in the Nucleus

Extrachromosomal DNA (ecDNA) promotes cancer by driving copy number heterogeneity and amplifying oncogenes along with functional enhancers. More recent studies suggest two additional mechanisms for further enhancing their oncogenic potential, one via forming ecDNA hubs to augment oncogene expression 1 and the other through acting as portable enhancers to trans-activate target genes 2. However, it has remained entirely elusive about how ecDNA explores the three-dimensional space of the nucleus and whether different ecDNA have distinct interacting mechanisms. Here, by profiling the DNA-DNA and DNA-RNA interactomes in tumor cells harboring different types of ecDNAs in comparison with similarly amplified homogenously staining regions (HSRs) in the chromosome, we show that specific ecDNA interactome is dictated by ecDNA-borne nascent RNA. We demonstrate that the ecDNA co-amplifying PVT1 and MYC utilize nascent noncoding PVT1 transcripts to mediate specific trans-activation of both ecDNA and chromosomal genes. In contrast, the ecDNA amplifying EGFR is weak in this property because of more efficient splicing to remove chromatin-associated nascent RNA. These findings reveal a noncoding RNA-orchestrated program hijacked by cancer cells to enhance the functional impact of amplified oncogenes and associated regulatory elements.

Enterococcus faecalis Strains with Compromised CRISPR-Cas Defense Emerge under Antibiotic Selection for a CRISPR-Targeted Plasmid

Enterococcus faecalis is a Gram-positive bacterium that natively colonizes the human gastrointestinal tract and opportunistically causes life-threatening infections. Multidrug-resistant (MDR) E. faecalis strains have emerged that are replete with mobile genetic elements (MGEs). Non-MDR E. faecalis strains frequently possess CRISPR-Cas systems, which reduce the frequency of MGE acquisition. We demonstrated in previous studies that E. faecalis populations can transiently maintain both a functional CRISPR-Cas system and a CRISPR-Cas target. In this study, we used serial passage and deep sequencing to analyze these populations. In the presence of antibiotic selection for the plasmid, mutants with compromised CRISPR-Cas defense and enhanced ability to acquire a second antibiotic resistance plasmid emerged. Conversely, in the absence of selection, the plasmid was lost from wild-type E. faecalis populations but not E. faecalis populations that lacked the cas9 gene. Our results indicate that E. faecalis CRISPR-Cas can become compromised under antibiotic selection, generating populations with enhanced abilities to undergo horizontal gene transfer. IMPORTANCE Enterococcus faecalis is a leading cause of hospital-acquired infections and disseminator of antibiotic resistance plasmids among Gram-positive bacteria. We have previously shown that E. faecalis strains with an active CRISPR-Cas system can prevent plasmid acquisition and thus limit the transmission of antibiotic resistance determinants. However, CRISPR-Cas is not a perfect barrier. In this study, we observed populations of E. faecalis with transient coexistence of CRISPR-Cas and one of its plasmid targets. Our experimental data demonstrate that antibiotic selection results in compromised E. faecalis CRISPR-Cas function, thereby facilitating the acquisition of additional resistance plasmids by E. faecalis.

[Development in diagnostic criteria for allergic bronchopulmonary aspergillosis]

Allergic bronchopulmonary aspergillosis (ABPA) is an allergic lung disease caused by the sensitization of Aspergillus fumigatus. In recent years, the research into ABPA has progressed, the testing methods have improved and the diagnostic criteria have been continuously updated. There is no gold standard for the diagnosis of the disease. The diagnostic criteria for ABPA include predisposing diseases, fungal-related immunoassay and pathological examination. Understanding the clinical significance of ABPA diagnostic criteria may help to prevent irreversible bronchopulmonary injury, improve respiratory function and improve the prognosis of patients.

[Clinical features of IgG4-related lung disease]

Objective: To explore the clinical features of IgG4-related lung disease. Methods: The clinical data of 60 patients diagnosed with IgG4-related lung disease in Peking University People's Hospital from February 2012 to May 2021 were retrospectively collected. Analysis was made to explore the features of clinical manifestation, laboratory, imaging, prognosis and other characteristics of the disease. Results: A total of 60 patients were included, with 40 males, age of (58.2±12.9) years, an age of onset of (57.1±13.2) years, and 31.7% (19 cases) of the patients had a history of allergic disease. 36.7% (22 cases) of the patients had respiratory symptoms during the disease. 94.6% (53/56) of patients had serum IgG4>1.35 g/L, 24.1% (14/58) of patients had increased eosinophils, 79.2% (38/48) of patients had increased IgE level, and 53.7% (29/54) of patients had decreased C3 or C4. Common imaging findings included nodular changes (38 cases, 63.3%), mediastinal and/or hilar lymphadenopathy (34 cases, 56.7%), and ground glass opacities (31 cases, 51.7%). Fifty-three cases (88.3%) showed two or more imaging changes. The pathological examination of the patient was mainly characterized by lymphoplasmacytic infiltration and fibrosis, with only one case of phlebitis obliterans. Compared with the asymptomatic group (38 cases), patients with respiratory symptoms (22 cases) showed higher level of serum total IgG and eosinophils (43.2 vs 17.8 g/L, 0.30×10⁹/L vs 0.14×10⁹/L, P<0.05), lower proportion of allergic diseases, and higher proportion of consolidation shadows on chest CT (P<0.05). There were no significant differences in serum IgG4, IgE, complement levels, and imaging outcomes after treatment between the two groups (P>0.05). Conclusions: The clinical manifestations of IgG4-related lung disease are atypical, and asymptomatic patients account for a high proportion. The imaging of the disease is highly heterogeneous, and patients are prone to show coexisted multiple imaging changes. The main clinical features and imaging outcomes of patients with and without respiratory symptoms are not significantly different.

SODB facilitates comprehensive exploration of spatial omics data

Spatial omics technologies generate wealthy but highly complex datasets. Here we present Spatial Omics DataBase (SODB), a web-based platform providing both rich data resources and a suite of interactive data analytical modules. SODB currently maintains >2,400 experiments from >25 spatial omics technologies, which are freely accessible as a unified data format compatible with various computational packages. SODB also provides multiple interactive data analytical modules, especially a unique module, Spatial Omics View (SOView). We conduct comprehensive statistical analyses and illustrate the utility of both basic and advanced analytical modules using multiple spatial omics datasets. We demonstrate SOView utility with brain spatial transcriptomics data and recover known anatomical structures. We further delineate functional tissue domains with associated marker genes that were obscured when analyzed using previous methods. We finally show how SODB may efficiently facilitate computational method development. The SODB website is https://gene.ai.tencent.com/SpatialOmics/ . The command-line package is available at https://pysodb.readthedocs.io/en/latest/ .

[Application of meandering mesenteric artery in vascular surgery]

In recent years, the therapeutic concepts and surgical methods of the meandering mesenteric artery involved in atherosclerosis have been reported in the literature, and the importance of this lateral branch circulation in the field of vascular surgery has received more attention. With the improvement of imaging techniques, the discovery rate of this collateral circulation increased. In the presence of major iliac artery occlusion, the meandering mesenteric artery may serve as an important collateral circulation to relieve ischemia in the visceral or lower extremities. The meandering mesenteric artery may also lead to type Ⅱ internal leakage after endovascular repair, while the embolization of inferior mesenteric artery may be performed through the meandering mesenteric artery. The role of meandering mesenteric arteries in vascular surgery needs further study.

Prostate Cancer Transcriptomic Regulation by the Interplay of Germline Risk Alleles, Somatic Mutations, and 3D Genomic Architecture

Prostate cancer is one of the most heritable human cancers. Genome-wide association studies have identified at least 185 prostate cancer germline risk alleles, most noncoding. We used integrative three-dimensional (3D) spatial genomics to identify the chromatin interaction targets of 45 prostate cancer risk alleles, 31 of which were associated with the transcriptional regulation of target genes in 565 localized prostate tumors. To supplement these 31, we verified transcriptional targets for 56 additional risk alleles using linear proximity and linkage disequilibrium analysis in localized prostate tumors. Some individual risk alleles influenced multiple target genes; others specifically influenced only distal genes while leaving proximal ones unaffected. Several risk alleles exhibited widespread germline-somatic interactions in transcriptional regulation, having different effects in tumors with loss of PTEN or RB1 relative to those without. These data clarify functional prostate cancer risk alleles in large linkage blocks and outline a strategy to model multidimensional transcriptional regulation.

SIGNIFICANCE

Many prostate cancer germline risk alleles are enriched in the noncoding regions of the genome and are hypothesized to regulate transcription. We present a 3D genomics framework to unravel risk SNP function and describe the widespread germline-somatic interplay in transcription control. This article is highlighted in the In This Issue feature, p. 2711.

SOTIP is a versatile method for microenvironment modeling with spatial omics data

The rapidly developing spatial omics generated datasets with diverse scales and modalities. However, most existing methods focus on modeling dynamics of single cells while ignore microenvironments (MEs). Here we present SOTIP (Spatial Omics mulTIPle-task analysis), a versatile method incorporating MEs and their interrelationships into a unified graph. Based on this graph, spatial heterogeneity quantification, spatial domain identification, differential microenvironment analysis, and other downstream tasks can be performed. We validate each module's accuracy, robustness, scalability and interpretability on various spatial omics datasets. In two independent mouse cerebral cortex spatial transcriptomics datasets, we reveal a gradient spatial heterogeneity pattern strongly correlated with the cortical depth. In human triple-negative breast cancer spatial proteomics datasets, we identify molecular polarizations and MEs associated with different patient survivals. Overall, by modeling biologically explainable MEs, SOTIP outperforms state-of-art methods and provides some perspectives for spatial omics data exploration and interpretation.

DAZL regulates proliferation of human primordial germ cells by direct binding to precursor miRNAs and enhances DICER processing activity

Understanding the molecular and cellular mechanisms of human primordial germ cells (hPGCs) is essential in studying infertility and germ cell tumorigenesis. Many RNA-binding proteins (RBPs) and non-coding RNAs are specifically expressed and functional during hPGC developments. However, the roles and regulatory mechanisms of these RBPs and non-coding RNAs, such as microRNAs (miRNAs), in hPGCs remain elusive. In this study, we reported a new regulatory function of DAZL, a germ cell-specific RBP, in miRNA biogenesis and cell proliferation. First, DAZL co-localized with miRNA let-7a in human PGCs and up-regulated the levels of &gt;100 mature miRNAs, including eight out of nine let-7 family, miR21, miR22, miR125, miR10 and miR199. Purified DAZL directly bound to the loops of precursor miRNAs with sequence specificity of GUU. The binding of DAZL to the precursor miRNA increased the maturation of miRNA by enhancing the cleavage activity of DICER. Furthermore, cell proliferation assay and cell cycle analysis confirmed that DAZL inhibited the proliferation of in vitro PGCs by promoting the maturation of these miRNAs. Evidently, the mature miRNAs up-regulated by DAZL silenced cell proliferation regulators including TRIM71. Moreover, DAZL inhibited germline tumor cell proliferation and teratoma formation. These results demonstrate that DAZL regulates hPGC proliferation by enhancing miRNA processing.

[Periodontitis and systemic diseases: an evidence-based interrelationship and treatment strategies for periodontitis with systemic diseases]

Periodontitis is a common chronic infectious disease, so as to be the leading cause of tooth loss in adults. Numerous studies have confirmed the interrelationship between periodontitis and systemic diseases. However, evidence-based reviews reporting the interrelationship between them and treatment strategies for periodontitis with systemic diseases were still absent currently. Therefore, based on evidence-based medical researches in recent years, this article will summarize the interrelationship between periodontitis and systemic diseases, and briefly state the treatment strategies for periodontitis with systemic diseases.

Avian influenza viruses suppress innate immunity by inducing trans-transcriptional readthrough via SSU72

Innate immunity plays critical antiviral roles. The highly virulent avian influenza viruses (AIVs) H5N1, H7N9, and H5N6 can better escape host innate immune responses than the less virulent seasonal H1N1 virus. Here, we report a mechanism by which transcriptional readthrough (TRT)-mediated suppression of innate immunity occurs post AIV infection. By using cell lines, mouse lungs, and patient PBMCs, we showed that genes on the complementary strand (“trans” genes) influenced by TRT were involved in the disruption of host antiviral responses during AIV infection. The trans-TRT enhanced viral lethality, and TRT abolishment increased cell viability and STAT1/2 expression. The viral NS1 protein directly bound to SSU72, and degradation of SSU72 induced TRT. SSU72 overexpression reduced TRT and alleviated mouse lung injury. Our results suggest that AIVs infection induce TRT by reducing SSU72 expression, thereby impairing host immune responses, a molecular mechanism acting through the NS1-SSU72-trans-TRT-STAT1/2 axis. Thus, restoration of SSU72 expression might be a potential strategy for preventing AIV pandemics.

CTCF functions as an insulator for somatic genes and a chromatin remodeler for pluripotency genes during reprogramming

CTCF mediates chromatin insulation and long-distance enhancer-promoter (EP) interactions; however, little is known about how these regulatory functions are partitioned among target genes in key biological processes. Here, we show that Ctcf expression is progressively increased during induced pluripotency. In this process, CTCF first functions as a chromatin insulator responsible for direct silencing of the somatic gene expression program and, interestingly, elevated Ctcf expression next ensures chromatin accessibility and contributes to increased EP interactions for a fraction of pluripotency-associated genes. Therefore, CTCF functions in a context-specific manner to modulate the 3D genome to enable cellular reprogramming. We further discover that these context-specific CTCF functions also enlist SMARCA5, an imitation switch (ISWI) chromatin remodeler, together rewiring the epigenome to facilitate cell-fate switch. These findings reveal the dual functions of CTCF in conjunction with a key chromatin remodeler to drive reprogramming toward pluripotency.

HPVMD-C: a disease-based mutation database of human papillomavirus in China

Human papillomavirus (HPV) can cause condyloma acuminatum and cervical cancer. Some mutations of these viruses are closely related to the persistent infection of cervical cancer and are ideal cancer vaccine targets. Several databases have been developed to collect HPV sequences, but no HPV mutation database has been published. This paper reports a Chinese HPV mutation database (HPVMD-C), which contains 149 HPV genotypes, 468 HPV mutations, 3409 protein sequences, 4727 domains and 236 epitopes. We analyzed the mutation distribution among HPV genotypes, domains and epitopes. We designed a visualization tool to display these mutations, domains and epitopes and provided more detailed information about the disease, region and related literature. We also proposed an HPV genotype prediction tool, which can predict HPV carcinogenic or non-carcinogenic risk genotypes. We expect that HPVMD-C will complement the existing database and provide valuable resources for HPV vaccine research and cervical cancer treatment. HPVMD-C is freely available at Database URL: http://bioinfo.zstu.edu.cn/hpv.

Integration of single cell data by disentangled representation learning

Recent developments of single cell RNA-sequencing technologies lead to the exponential growth of single cell sequencing datasets across different conditions. Combining these datasets helps to better understand cellular identity and function. However, it is challenging to integrate different datasets from different laboratories or technologies due to batch effect, which are interspersed with biological variances. To overcome this problem, we have proposed Single Cell Integration by Disentangled Representation Learning (SCIDRL), a domain adaption-based method, to learn low-dimensional representations invariant to batch effect. This method can efficiently remove batch effect while retaining cell type purity. We applied it to thirteen diverse simulated and real datasets. Benchmark results show that SCIDRL outperforms other methods in most cases and exhibits excellent performances in two common situations: (i) effective integration of batch-shared rare cell types and preservation of batch-specific rare cell types; (ii) reliable integration of datasets with different cell compositions. This demonstrates SCIDRL will offer a valuable tool for researchers to decode the enigma of cell heterogeneity.

MarkovHC: Markov hierarchical clustering for the topological structure of high-dimensional single-cell omics data with transition pathway and critical point detection

Clustering cells and depicting the lineage relationship among cell subpopulations are fundamental tasks in single-cell omics studies. However, existing analytical methods face challenges in stratifying cells, tracking cellular trajectories, and identifying critical points of cell transitions. To overcome these, we proposed a novel Markov hierarchical clustering algorithm (MarkovHC), a topological clustering method that leverages the metastability of exponentially perturbed Markov chains for systematically reconstructing the cellular landscape. Briefly, MarkovHC starts with local connectivity and density derived from the input and outputs a hierarchical structure for the data. We firstly benchmarked MarkovHC on five simulated datasets and ten public single-cell datasets with known labels. Then, we used MarkovHC to investigate the multi-level architectures and transition processes during human embryo preimplantation development and gastric cancer procession. MarkovHC found heterogeneous cell states and sub-cell types in lineage-specific progenitor cells and revealed the most possible transition paths and critical points in the cellular processes. These results demonstrated MarkovHC's effectiveness in facilitating the stratification of cells, identification of cell populations, and characterization of cellular trajectories and critical points.

MyoD is a 3D genome structure organizer for muscle cell identity

The genome exists as an organized, three-dimensional (3D) dynamic architecture, and each cell type has a unique 3D genome organization that determines its cell identity. An unresolved question is how cell type-specific 3D genome structures are established during development. Here, we analyzed 3D genome structures in muscle cells from mice lacking the muscle lineage transcription factor (TF), MyoD, versus wild-type mice. We show that MyoD functions as a “genome organizer” that specifies 3D genome architecture unique to muscle cell development, and that H3K27ac is insufficient for the establishment of MyoD-induced chromatin loops in muscle cells. Moreover, we present evidence that other cell lineage-specific TFs might also exert functional roles in orchestrating lineage-specific 3D genome organization during development.

Pioneer transcription factors (TFs) have been proposed to act as protein anchors to orchestrate cell type-specific 3D genome architecture. MyoD is a pioneer TF for myogenic lineage specification. Here the authors provide further support for the role of MyoD in 3D genome architecture in muscle stem cells by comparing MyoD knockout and wild-type mice.

A data-driven method to learn a jump diffusion process from aggregate biological gene expression data

Dynamic models of gene expression are urgently required. In this paper, we describe the time evolution of gene expression by learning a jump diffusion process to model the biological process directly. Our algorithm needs aggregate gene expression data as input and outputs the parameters of the jump diffusion process. The learned jump diffusion process can predict population distributions of gene expression at any developmental stage, obtain long-time trajectories for individual cells, and offer a novel approach to computing RNA velocity. Moreover, it studies biological systems from a stochastic dynamic perspective. Gene expression data at a time point, which is a snapshot of a cellular process, is treated as an empirical marginal distribution of a stochastic process. The Wasserstein distance between the empirical distribution and predicted distribution by the jump diffusion process is minimized to learn the dynamics. For the learned jump diffusion process, its trajectories correspond to the development process of cells, the stochasticity determines the heterogeneity of cells, its instantaneous rate of state change can be taken as "RNA velocity", and the changes in scales and orientations of clusters can be noticed too. We demonstrate that our method can recover the underlying nonlinear dynamics better compared to previous parametric models and the diffusion processes driven by Brownian motion for both synthetic and real world datasets. Our method is also robust to perturbations of data because the computation involves only population expectations.

ScaffComb: A Phenotype-Based Framework for Drug Combination Virtual Screening in Large-Scale Chemical Datasets

Combinational therapy is used for a long time in cancer treatment to overcome drug resistance related to monotherapy. Increased pharmacological data and the rapid development of deep learning methods have enabled the construction of models to predict and screen drug pairs. However, the size of drug libraries is restricted to hundreds to thousands of compounds. The ScaffComb framework, which aims to bridge the gaps in the virtual screening of drug combinations in large-scale databases, is proposed here. Inspired by phenotype-based drug design, ScaffComb integrates phenotypic information into molecular scaffolds, which can be used to screen the drug library and identify potent drug combinations. First, ScaffComb is validated using the US food and drug administration dataset and known drug combinations are successfully reidentified. Then, ScaffComb is applied to screen the ZINC and ChEMBL databases, which yield novel drug combinations and reveal an ability to discover new synergistic mechanisms. To our knowledge, ScaffComb is the first method to use phenotype-based virtual screening of drug combinations in large-scale chemical datasets.

[Transient spontaneous regression of pulmonary lymphoma: two case report and review of the literature]

Objective: To investigate the clinical features and pathogenesis of spontaneous regression of lymphoma involving the respiratory system. Methods: Two cases of pulmonary lymphoma which regressed spontaneously were reported. Literatures with"spontaneous regression, lymphoma""spontaneous remission, lymphoma"published before February 2020 were searched in Wanfang, CNKI and Pubmed database. And cases with respiratory system involvement of lymphoma which regressed spontaneously were analyzed. A total of 22 cases were finally retrieved. There were 6 males and 16 females, with an average age of (62.5±16.8) years. Results: Common symptoms included dyspnea, cough, expectoration, fever and weakness. Imaging examination showed that any parts in respiratory system could be involved. The proportion of invasive and indolent lymphomas was approximately similar. The time interval between diagnosis of lymphoma and first sign of spontaneous regression of the disease was from 2 weeks to 1 year. Spontaneous regression of the disease could sustain from 50 days to 60 months. Moreover, 78.6% of patients achieved complete remission in their courses of treatment, and some patients showed wax and wane phenomena. Conclusions: Spontaneous regression of lymphoma may occur in any part of the respiratory system. When a patient exhibits spontaneous regression of pulmonary lesions, lymphoma should be considered in the differential diagnoses. Pathology study is needed to achieve a definitive diagnosis. Misdiagnosis and delayed diagnosis related to empirical treatment of antibiotics and corticosteroid should be avoided.

SEAM is a spatial single nuclear metabolomics method for dissecting tissue microenvironment

Spatial metabolomics can reveal intercellular heterogeneity and tissue organization. Here we report on the spatial single nuclear metabolomics (SEAM) method, a flexible platform combining high-spatial-resolution imaging mass spectrometry and a set of computational algorithms that can display multiscale and multicolor tissue tomography together with identification and clustering of single nuclei by their in situ metabolic fingerprints. We first applied SEAM to a range of wild-type mouse tissues, then delineated a consistent pattern of metabolic zonation in mouse liver. We further studied the spatial metabolic profile in the human fibrotic liver. We discovered subpopulations of hepatocytes with special metabolic features associated with their proximity to the fibrotic niche, and validated this finding by spatial transcriptomics with Geo-seq. These demonstrations highlighted SEAM's ability to explore the spatial metabolic profile and tissue histology at the single-cell level, leading to a deeper understanding of tissue metabolic organization.

Characterizing microRNA-mediated modulation of gene expression noise and its effect on synthetic gene circuits

MicroRNAs (miRNAs) have been shown to modulate gene expression noise, but less is known about how miRNAs with different properties may regulate noise differently. Here, we investigate the role of competing RNAs and the composition of miRNA response elements (MREs) in modulating noise. We find that weak competing RNAs could introduce lower noise than strong competing RNAs. In comparison with a single MRE, both repetitive and composite MREs can reduce the noise at low expression, but repetitive MREs can elevate the noise remarkably at high expression. We further observed the behavior of a synthetic cell-type classifier with miRNAs as inputs and find that miRNAs and MREs that could introduce higher noise tend to enhance cell state transition. These results provide a systematic and quantitative understanding of the function of miRNAs in controlling gene expression noise and the utilization of miRNAs to modulate the behavior of synthetic gene circuits.

Quantifying the phase separation property of chromatin-associated proteins under physiological conditions using an anti-1,6-hexanediol index

BACKGROUND

Liquid-liquid phase separation (LLPS) is an important organizing principle for biomolecular condensation and chromosome compartmentalization. However, while many proteins have been reported to undergo LLPS, quantitative and global analysis of chromatin LLPS property remains absent.

RESULTS

Here, by combining chromatin-associated protein pull-down, quantitative proteomics and 1,6-hexanediol (1,6-HD) treatment, we develop Hi-MS and define an anti-1,6-HD index of chromatin-associated proteins (AICAP) to quantify 1,6-HD sensitivity of chromatin-associated proteins under physiological conditions. Compared with known physicochemical properties involved in phase separation, we find that proteins with lower AICAP are associated with higher content of disordered regions, higher hydrophobic residue preference, higher mobility and higher predicted LLPS potential. We also construct BL-Hi-C libraries following 1,6-HD treatment to study the sensitivity of chromatin conformation to 1,6-HD treatment. We find that the active chromatin and high-order structures, as well as the proteins enriched in corresponding regions, are more sensitive to 1,6-HD treatment.

CONCLUSIONS

Our work provides a global quantitative measurement of LLPS properties of chromatin-associated proteins and higher-order chromatin structure. Hi-MS and AICAP data provide an experimental tool and quantitative resources valuable for future studies of biomolecular condensates.

Single-cell alternative polyadenylation analysis delineates GABAergic neuron types

Background

Alternative polyadenylation (APA) is emerging as an important mechanism in the post-transcriptional regulation of gene expression across eukaryotic species. Recent studies have shown that APA plays key roles in biological processes, such as cell proliferation and differentiation. Single-cell RNA-seq technologies are widely used in gene expression heterogeneity studies; however, systematic studies of APA at the single-cell level are still lacking.

Results

Here, we described a novel computational framework, SAPAS, that utilizes 3′-tag-based scRNA-seq data to identify novel poly(A) sites and quantify APA at the single-cell level. Applying SAPAS to the scRNA-seq data of phenotype characterized GABAergic interneurons, we identified cell type-specific APA events for different GABAergic neuron types. Genes with cell type-specific APA events are enriched for synaptic architecture and communications. In further, we observed a strong enrichment of heritability for several psychiatric disorders and brain traits in altered 3′ UTRs and coding sequences of cell type-specific APA events. Finally, by exploring the modalities of APA, we discovered that the bimodal APA pattern of Pak3 could classify chandelier cells into different subpopulations that are from different laminar positions.

Conclusions

We established a method to characterize APA at the single-cell level. When applied to a scRNA-seq dataset of GABAergic interneurons, the single-cell APA analysis not only identified cell type-specific APA events but also revealed that the modality of APA could classify cell subpopulations. Thus, SAPAS will expand our understanding of cellular heterogeneity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01076-3.

DAGM: A novel modelling framework to assess the risk of HER2-negative breast cancer based on germline rare coding mutations

BACKGROUND

Breast cancers can be divided into HER2-negative and HER2-positive subtypes according to different status of HER2 gene. Despite extensive studies connecting germline mutations with possible risk of HER2-negative breast cancer, the main category of breast cancer, it remains challenging to obtain accurate risk assessment and to understand the potential underlying mechanisms.

METHODS

We developed a novel framework named Damage Assessment of Genomic Mutations (DAGM), which projects rare coding mutations and gene expressions into Activity Profiles of Signalling Pathways (APSPs).

FINDINGS

We characterized and validated DAGM framework at multiple levels. Based on an input of germline rare coding mutations, we obtained the corresponding APSP spectrum to calculate the APSP risk score, which was capable of distinguish HER2-negative from HER2-positive cases. These findings were validated using breast cancer data from TCGA (AUC = 0.7). DAGM revealed that HER2 signalling pathway was up-regulated in germline of HER2-negative patients, and those with high APSP risk scores had exhibited immune suppression. These findings were validated using RNA sequencing, phosphoproteome analysis, and CyTOF. Moreover, using germline mutations, DAGM could evaluate the risk for HER2-negative breast cancer, not only in women carrying BRCA1/2 mutations, but also in those without known disease-associated mutations.

INTERPRETATION

The DAGM can facilitate the screening of subjects at high risk of HER2-negative breast cancer for primary prevention. This study also provides new insights into the potential mechanisms of developing HER2-negative breast cancer. The DAGM has the potential to be applied in the prevention, diagnosis, and treatment of HER2-negative breast cancer.

FUNDING

This work was supported by the National Key Research and Development Program of China (grant no. 2018YFC0910406 and 2018AAA0103302 to CZ); the National Natural Science Foundation of China (grant no. 81202076 and 82072939 to MY, 81871513 to KW); the Guangzhou Science and Technology Program key projects (grant no. 2014J2200007 to MY, 202002030236 to KW); the National Key R&D Program of China (grant no. 2017YFC1309100 to CL); Shenzhen Science and Technology Planning Project (grant no. JCYJ20170817095211560 574 to YN); and the Natural Science Foundation of Guangdong Province (grant no. 2017A030313882 to KW and S2013010012048 to MY); Hefei National Laboratory for Physical Sciences at the Microscale (grant no. KF2020009 to GN); and RGC General Research Fund (grant no. 17114519 to YQS).

The loss of heterochromatin is associated with multiscale three-dimensional genome reorganization and aberrant transcription during cellular senescence

Heterochromatin remodeling is critical for various cell processes. In particular, the "loss of heterochromatin" phenotype in cellular senescence is associated with the process of aging and age-related disorders. Although biological processes of senescent cells, including senescence-associated heterochromatin foci (SAHF) formation, chromosome compaction, and redistribution of key proteins, have been closely associated with high-order chromatin structure, the relationship between the high-order chromatin reorganization and the loss of heterochromatin phenotype during senescence has not been fully understood. By using senescent and deep senescent fibroblasts induced by DNA damage harboring the "loss of heterochromatin" phenotype, we observed progressive 3D reorganization of heterochromatin during senescence. Facultative and constitutive heterochromatin marked by H3K27me3 and H3K9me3, respectively, show different alterations. Facultative heterochromatin tends to switch from the repressive B-compartment to the active A-compartment, whereas constitutive heterochromatin shows no significant changes at the compartment level but enhanced interactions between themselves. Both types of heterochromatin show increased chromatin accessibility and gene expression leakage during senescence. Furthermore, increased chromatin accessibility in potential CTCF binding sites accompanies the establishment of novel loops in constitutive heterochromatin. Finally, we also observed aberrant expression of repetitive elements, including LTR (long terminal repeat) and satellite classes. Overall, facultative and constitutive heterochromatin show both similar and distinct multiscale alterations in the 3D map, chromatin accessibility, and gene expression leakage. This study provides an epigenomic map of heterochromatin reorganization during senescence.

Deciphering hierarchical organization of topologically associated domains through change-point testing

Background

The nucleus of eukaryotic cells spatially packages chromosomes into a hierarchical and distinct segregation that plays critical roles in maintaining transcription regulation. High-throughput methods of chromosome conformation capture, such as Hi-C, have revealed topologically associating domains (TADs) that are defined by biased chromatin interactions within them.

Results

We introduce a novel method, HiCKey, to decipher hierarchical TAD structures in Hi-C data and compare them across samples. We first derive a generalized likelihood-ratio (GLR) test for detecting change-points in an interaction matrix that follows a negative binomial distribution or general mixture distribution. We then employ several optimal search strategies to decipher hierarchical TADs with p values calculated by the GLR test. Large-scale validations of simulation data show that HiCKey has good precision in recalling known TADs and is robust against random collisions of chromatin interactions. By applying HiCKey to Hi-C data of seven human cell lines, we identified multiple layers of TAD organization among them, but the vast majority had no more than four layers. In particular, we found that TAD boundaries are significantly enriched in active chromosomal regions compared to repressed regions.

Conclusions

HiCKey is optimized for processing large matrices constructed from high-resolution Hi-C experiments. The method and theoretical result of the GLR test provide a general framework for significance testing of similar experimental chromatin interaction data that may not fully follow negative binomial distributions but rather more general mixture distributions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04113-8.

Integrative Genomic Analysis of Gemcitabine Resistance in Pancreatic Cancer by Patient-derived Xenograft Models

PURPOSE

Gemcitabine is most commonly used for pancreatic cancer. However, the molecular features and mechanisms of the frequently occurring resistance remain unclear. This work aims at exploring the molecular features of gemcitabine resistance and identifying candidate biomarkers and combinatorial targets for the treatment.

EXPERIMENTAL DESIGN

In this study, we established 66 patient-derived xenografts (PDXs) on the basis of clinical pancreatic cancer specimens and treated them with gemcitabine. We generated multiomics data (including whole-exome sequencing, RNA sequencing, miRNA sequencing, and DNA methylation array) of 15 drug-sensitive and 13 -resistant PDXs before and after the gemcitabine treatment. We performed integrative computational analysis to identify the molecular networks related to gemcitabine intrinsic and acquired resistance. Then, short hairpin RNA-based high-content screening was implemented to validate the function of the deregulated genes.

RESULTS

The comprehensive multiomics analysis and functional experiment revealed that MRPS5 and GSPT1 had strong effects on cell proliferation, and CD55 and DHTKD1 contributed to gemcitabine resistance in pancreatic cancer cells. Moreover, we found miR-135a-5p was significantly associated with the prognosis of patients with pancreatic cancer and could be a candidate biomarker to predict gemcitabine response. Comparing the molecular features before and after the treatment, we found that PI3K-Akt, p53, and hypoxia-inducible factor-1 pathways were significantly altered in multiple patients, providing candidate target pathways for reducing the acquired resistance.

CONCLUSIONS

This integrative genomic study systematically investigated the predictive markers and molecular mechanisms of chemoresistance in pancreatic cancer and provides potential therapy targets for overcoming gemcitabine resistance.

Abstract PS8-29: Rare variants in the germline genome holistically determine receptor-independent Her2 signaling pathway activation and immune suppression, shaping pathological type and risk of HER2-negative breast cancer

Background: Pathogenic factors embedded in the germline genome are widely recognized as being crucial to breast cancer development. However, current knowledge is either concentrated on the pathogenic variants of a few individual genes or SNPs distributed sparsely across the genome in non-coding regions. Methods: We developed a multi-layered framework, DAGG, which converts somatic mutations or germline rare coding variants (gRCVs) into a functional spectrum of dozens of cellular functions and signaling pathways to identify potential pathogenic factors.Findings: We analyzed whole-exome sequencing (WES) data of 726 germline DNA samples and 169 breast tumor DNA samples from breast cancer patients with various pathological types and cancer-free female subjects, we found that germline pathogens of breast cancers were (1) mainly distributed in HER2-negative subtypes, and (2) involved Her2 signaling pathway activation and immune suppression. These computational discoveries were experimentally validated and can provide digital features to explain the germline differences between diseased and healthy genome (AUC = 0.76). Furthermore, an individual’s risk for breast cancer can be estimated by calculating the combined effects of these identified germline pathogens. Carriers of BRCA1/2 pathogenic variants were found to have a significantly higher average risk (p = 0.02). Interpretation: The results demonstrated that the identified pathogenic mechanisms by DAGG were compatible with our current understanding of the causes of breast cancer. Moreover, DAGG provides improved performance over currently used polygenic risk score method of measuring complex disease risks. Our framework promises possible future applications for the prevention, diagnosis, and treatment of breast cancer.

Citation Format: Mei Yang, Yanhui Fan, Zhi-Yong Wu, Zhendong Feng, Qiangzu Zhang, Shunhua Han, Zhonghai Zhang, Xu Li, Yiqing Xue, Xiaoling Li, Meixia Hu, Jieqing Li, Weiping Li, Hongfei Gao, Ciqiu Yang, Chunming Zhang, Liulu Zhang, Teng Zhu, Minyi Cheng, Fei Ji, Juntao Xu, Hening Cui, Guangming Tan, Michael Q Zhang, Changhong Liang, Zaiyi Liu, You-Qiang Song, Gang Niu, Kun Wang. Rare variants in the germline genome holistically determine receptor-independent Her2 signaling pathway activation and immune suppression, shaping pathological type and risk of HER2-negative breast cancer [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS8-29.

The landscape of RNA polymerase II-associated chromatin interactions in prostate cancer

Transcriptional dysregulation is a hallmark of prostate cancer (PCa). We mapped the RNA polymerase II-associated (RNA Pol II-associated) chromatin interactions in normal prostate cells and PCa cells. We discovered thousands of enhancer-promoter, enhancer-enhancer, as well as promoter-promoter chromatin interactions. These transcriptional hubs operate within the framework set by structural proteins - CTCF and cohesins - and are regulated by the cooperative action of master transcription factors, such as the androgen receptor (AR) and FOXA1. By combining analyses from metastatic castration-resistant PCa (mCRPC) specimens, we show that AR locus amplification contributes to the transcriptional upregulation of the AR gene by increasing the total number of chromatin interaction modules comprising the AR gene and its distal enhancer. We deconvoluted the transcription control modules of several PCa genes, notably the biomarker KLK3, lineage-restricted genes (KRT8, KRT18, HOXB13, FOXA1, ZBTB16), the drug target EZH2, and the oncogene MYC. By integrating clinical PCa data, we defined a germline-somatic interplay between the PCa risk allele rs684232 and the somatically acquired TMPRSS2-ERG gene fusion in the transcriptional regulation of multiple target genes - VPS53, FAM57A, and GEMIN4. Our studies implicate changes in genome organization as a critical determinant of aberrant transcriptional regulation in PCa.

A novel cytogenetic method to image chromatin interactions at subkilobase resolution: Tn5 transposase-based fluorescence in situ hybridization

There is an increasing interest in understanding how three-dimensional organization of the genome is regulated. Different strategies have been used to identify genome-wide chromatin interactions. However, owing to current limitations in resolving genomic contacts, visualization and validation of these genomic loci at subkilobase resolution remain unsolved to date. Here, we describe Tn5 transposase-based fluorescence in situ hybridization (Tn5-FISH), a polymerase chain reaction-based, cost-effective imaging method, which can colocalize the genomic loci at subkilobase resolution, dissect genome architecture, and verify chromatin interactions detected by chromatin configuration capture-derived methods. To validate this method, short-range interactions in the keratin-encoding gene (KRT) locus in the topologically associated domain were imaged by triple-color Tn5-FISH, indicating that Tn5-FISH is very useful to verify short-range chromatin interactions inside the contact domain and TAD. Therefore, Tn5-FISH can be a powerful molecular tool for clinical detection of cytogenetic changes in numerous genetic diseases such as cancers.

Molecular basis for histone H3 "K4me3-K9me3/2" methylation pattern readout by Spindlin1

Histone recognition by "reader" modules serves as a fundamental mechanism in epigenetic regulation. Previous studies have shown that Spindlin1 is a reader of histone H3K4me3 as well as "K4me3-R8me2a" and promotes transcription of rDNA or Wnt/TCF4 target genes. Here we show that Spindlin1 also acts as a potent reader of histone H3 "K4me3-K9me3/2" bivalent methylation pattern. Calorimetric titration revealed a binding affinity of 16 nm between Spindlin1 and H3 "K4me3-K9me3" peptide, which is one to three orders of magnitude stronger than most other histone readout events at peptide level. Structural studies revealed concurrent recognition of H3K4me3 and H3K9me3/2 by aromatic pockets 2 and 1 of Spindlin1, respectively. Epigenomic profiling studies showed that Spindlin1 colocalizes with both H3K4me3 and H3K9me3 peaks in a subset of genes enriched in biological processes of transcription and its regulation. Moreover, the distribution of Spindlin1 peaks is primarily associated with H3K4me3 but not H3K9me3, which suggests that Spindlin1 is a downstream effector of H3K4me3 generated in heterochromatic regions. Collectively, our work calls attention to an intriguing function of Spindlin1 as a potent H3 "K4me3-K9me3/2" bivalent mark reader, thereby balancing gene expression and silencing in H3K9me3/2-enriched regions.

Highly diversified core promoters in the human genome and their effects on gene expression and disease predisposition

Background

Core promoter controls transcription initiation. However, little is known for core promoter diversity in the human genome and its relationship with diseases. We hypothesized that as a functional important component in the genome, the core promoter in the human genome could be under evolutionary selection, as reflected by its highly diversification in order to adjust gene expression for better adaptation to the different environment.

Results

Applying the “Exome-based Variant Detection in Core-promoters” method, we analyzed human core-promoter diversity by using the 2682 exome data sets of 25 worldwide human populations sequenced by the 1000 Genome Project. Collectively, we identified 31,996 variants in the core promoter region (− 100 to + 100) of 12,509 human genes (https://dbhcpd.fhs.um.edu.mo). Analyzing the rich variation data identified highly ethnic-specific patterns of core promoter variation between different ethnic populations, the genes with highly variable core promoters, the motifs affected by the variants, and their involved functional pathways. eQTL test revealed that 12% of core promoter variants can significantly alter gene expression level. Comparison with GWAS data we located 163 variants as the GWAS identified traits associated with multiple diseases, half of these variants can alter gene expression.

Conclusion

Data from our study reals the highly diversified nature of core promoter in the human genome, and highlights that core promoter variation could play important roles not only in gene expression regulation but also in disease predisposition.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07222-5.

LncRNA Platr22 promotes super-enhancer activity and stem cell pluripotency

Super-enhancers (SEs) comprise large clusters of enhancers, which are co-occupied by multiple lineage-specific and master transcription factors, and play pivotal roles in regulating gene expression and cell fate determination. However, it is still largely unknown whether and how SEs are regulated by the noncoding portion of the genome. Here, through genome-wide analysis, we found that long noncoding RNA (lncRNA) genes preferentially lie next to SEs. In mouse embryonic stem cells (mESCs), depletion of SE-associated lncRNA transcripts dysregulated the activity of their nearby SEs. Specifically, we revealed a critical regulatory role of the lncRNA gene Platr22 in modulating the activity of a nearby SE and the expression of the nearby pluripotency regulator ZFP281. Through these regulatory events, Platr22 contributes to pluripotency maintenance and proper differentiation of mESCs. Mechanistically, Platr22 transcripts coat chromatin near the SE region and interact with DDX5 and hnRNP-L. DDX5 further recruits p300 and other factors related to active transcription. We propose that these factors assemble into a transcription hub, thus promoting an open and active epigenetic chromatin state. Our study highlights an unanticipated role for a class of lncRNAs in epigenetically controlling the activity and vulnerability to perturbation of nearby SEs for cell fate determination.

Model-based analysis of chromatin interactions from dCas9-Based CAPTURE-3C-seq

Deciphering long-range chromatin interactions is critical for understanding temporal and tissue-specific gene expression regulated by cis- and trans-acting factors. By combining the chromosome conformation capture (3C) and biotinylated dCas9 system, we previously established a method CAPTURE-3C-seq to unbiasedly identify high-resolution and locus-specific long-range DNA interactions. Here we present the statistical model and a flexible pipeline, C3S, for analysing CAPTURE-3C-seq or similar experimental data from raw sequencing reads to significantly interacting chromatin loci. C3S provides all steps for data processing, quality control and result illustration. It can automatically define the bin size based on the binding peak of the dCas9-targeted regions. Furthermore, it supports the analysis of intra- and inter-chromosomal interactions for different mammalian cell types. We successfully applied C3S across multiple datasets in human K562 cells and mouse embryonic stem cells (mESC) for detecting known and new chromatin interactions at multiple scales. Integrative and topological analysis of the interacted loci at the human β-globin gene cluster provides new insights into mechanisms in developmental gene regulation and network structure in local chromosomal architecture. Furthermore, computational results in mESCs reveal a role for chromatin interacting loops between enhancers and promoters in regulating alternative transcripts of the pluripotency gene OCT4.

A large-scale CRISPR screen and identification of essential genes in cellular senescence bypass

Cellular senescence is an important mechanism of autonomous tumor suppression, while its consequence such as the senescence-associated secretory phenotype (SASP) may drive tumorigenesis and age-related diseases. Therefore, controlling the cell fate optimally when encountering senescence stress is helpful for anti-cancer or anti-aging treatments. To identify genes essential for senescence establishment or maintenance, we carried out a CRISPR-based screen with a deliberately designed single-guide RNA (sgRNA) library. The library comprised of about 12,000 kinds of sgRNAs targeting 1378 senescence-associated genes selected by integrating the information of literature mining, protein-protein interaction network, and differential gene expression. We successfully detected a dozen gene deficiencies potentially causing senescence bypass, and their phenotypes were further validated with a high true positive rate. RNA-seq analysis showed distinct transcriptome patterns of these bypass cells. Interestingly, in the bypass cells, the expression of SASP genes was maintained or elevated with CHEK2, HAS1, or MDK deficiency; but neutralized with MTOR, CRISPLD2, or MORF4L1 deficiency. Pathways of some age-related neurodegenerative disorders were also downregulated with MTOR, CRISPLD2, or MORF4L1 deficiency. The results demonstrated that disturbing these genes could lead to distinct cell fates as a consequence of senescence bypass, suggesting that they may play essential roles in cellular senescence.

[Clinical features of 2019 novel coronavirus pneumonia in the early stage from a fever clinic in Beijing]

Objective: To summarize and analyze the clinical and imaging characteristics of patients with 2019 novel coronavirus pneumonia in the early stage in Beijing. Methods: A retrospective analysis of clinical and imaging data of 9 patients with 2019 novel coronavirus infection diagnosed in one fever clinicic in Beijing from January 18, 2020 to February 3, 2020. Results: 5 male and 4 female was included in those 9 patients, whose median age was 36 years, and the age range from 15 to 49 years. 8 of these patients had no underlying disease and one suffered from diabetes. 7 patients had a history of travel to Wuhan City or Hubei Province, and one patient was a medical staff. Two family clustered was found. The incubation period was 1 to 6 days. The clinical manifestations were fever in 8 cases (8/9) , dry cough in 5 cases (5/9) , pharyngalgia in 4 cases (4/9) , fatigue in 4 cases (4/9) , body soreness in 4 cases (4/9) , and blocked or watery nose in 1 case (1/9) . Six patients (6/9) had abnormal cell peripheral blood, of which 3 (3/9) had an increased monocyte count, 2 (2/9) had a reduced lymphocyte, and 1 (1/9) had an increased leukocyte count, while the 3 patients had normal cell blood routines. The median of CRP was 16.3 mg/L, including 5 patients with slightly elevated (5/9) , 4 patients with normal values (4/9) . the results of procalcitonin test were negative in5 patients. Three patients were examined by chest X-ray examination, one of which was normal, one case showed infiltrates of right upper lung, and another showed in right lower lung. All patients underwent chest HRCT. And 7 cases (7/9) showed multiple ground glass exudation, including 5 cases (5/7) involved bilateral lungs, 2 cases (2/7) involved unilateral lung, 3 cases (3/7) with patchy consolidation, and 2 cases (2/9) showed no abnormality. Conclusions: The patents with 2019 novel coronavirus pneumonia in this study generally have an epidemiological history. The clinical manifestations are fever and cough. Peripheral white blood cell counts were most normal And PCT were all negative. Chest HRCT manifested as multiple ground-glass opacities with partly consolidation. Some patients had normal chest radiographs but HRCT showed pneumonia. Some patients had no pneumonia on chest HRCT.

Single‑cell RNA sequencing of t(8;21) acute myeloid leukemia for risk prediction

Single-cell RNA sequencing (scRNA-seq) of bone marrow or peripheral blood samples from patients with acute myeloid leukemia (AML) enables the characterization of heterogeneous malignant cells. A total of 87 cells from two patients with t(8;21) AML were analyzed using scRNA-seq. Clustering methods were used to separate leukemia cells into different sub-populations, and the expression patterns of specific marker genes were used to annotate these populations. Among the 31 differentially expressed genes in the cells of a patient who relapsed after hematopoietic stem cell transplantation, 13 genes were identified to be associated with leukemia. Furthermore, three genes, namely AT-rich interaction domain 2, lysine methyltransferase 2A and synaptotagmin binding cytoplasmic RNA interacting protein were validated as possible prognostic biomarkers using two bulk expression datasets. Taking advantage of scRNA-seq, the results of the present study may provide clinicians with several possible biomarkers to predict the prognostic outcomes of t(8;21) AML.