Detecting novel cell type in single-cell chromatin accessibility data via open-set domain adaptation.

Recent advances in single-cell technologies enable the rapid growth of multi-omics data. Cell type annotation is one common task in analyzing single-cell data. It is a challenge that some cell types in the testing set are not present in the training set (i.e. unknown cell types). Most scATAC-seq cell type annotation methods generally assign each cell in the testing set to one known type in the training set but neglect unknown cell types. Here, we present OVAAnno, an automatic cell types annotation method which utilizes open-set domain adaptation to detect unknown cell types in scATAC-seq data. Comprehensive experiments show that OVAAnno successfully identifies known and unknown cell types. Further experiments demonstrate that OVAAnno also performs well on scRNA-seq data. Our codes are available online at https://github.com/lisaber/OVAAnno/tree/master.

Comprehensive single-cell RNA-seq analysis using deep interpretable generative modeling guided by biological hierarchy knowledge.

Recent advances in microfluidics and sequencing technologies allow researchers to explore cellular heterogeneity at single-cell resolution. In recent years, deep learning frameworks, such as generative models, have brought great changes to the analysis of transcriptomic data. Nevertheless, relying on the potential space of these generative models alone is insufficient to generate biological explanations. In addition, most of the previous work based on generative models is limited to shallow neural networks with one to three layers of latent variables, which may limit the capabilities of the models. Here, we propose a deep interpretable generative model called d-scIGM for single-cell data analysis. d-scIGM combines sawtooth connectivity techniques and residual networks, thereby constructing a deep generative framework. In addition, d-scIGM incorporates hierarchical prior knowledge of biological domains to enhance the interpretability of the model. We show that d-scIGM achieves excellent performance in a variety of fundamental tasks, including clustering, visualization, and pseudo-temporal inference. Through topic pathway studies, we found that d-scIGM-learned topics are better enriched for biologically meaningful pathways compared to the baseline models. Furthermore, the analysis of drug response data shows that d-scIGM can capture drug response patterns in large-scale experiments, which provides a promising way to elucidate the underlying biological mechanisms. Lastly, in the melanoma dataset, d-scIGM accurately identified different cell types and revealed multiple melanin-related driver genes and key pathways, which are critical for understanding disease mechanisms and drug development.

Subgraph extraction and graph representation learning for single cell Hi-C imputation and clustering.

Single-cell Hi-C (scHi-C) technology enables the investigation of 3D chromatin structure variability across individual cells. However, the analysis of scHi-C data is challenged by a large number of missing values. Here, we present a scHi-C data imputation model HiC-SGL, based on Subgraph extraction and graph representation learning. HiC-SGL can also learn informative low-dimensional embeddings of cells. We demonstrate that our method surpasses existing methods in terms of imputation accuracy and clustering performance by various metrics.

Benchmarking deep learning methods for predicting CRISPR/Cas9 sgRNA on- and off-target activities.

In silico design of single guide RNA (sgRNA) plays a critical role in clustered regularly interspaced, short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) system. Continuous efforts are aimed at improving sgRNA design with efficient on-target activity and reduced off-target mutations. In the last 5 years, an increasing number of deep learning-based methods have achieved breakthrough performance in predicting sgRNA on- and off-target activities. Nevertheless, it is worthwhile to systematically evaluate these methods for their predictive abilities. In this review, we conducted a systematic survey on the progress in prediction of on- and off-target editing. We investigated the performances of 10 mainstream deep learning-based on-target predictors using nine public datasets with different sample sizes. We found that in most scenarios, these methods showed superior predictive power on large- and medium-scale datasets than on small-scale datasets. In addition, we performed unbiased experiments to provide in-depth comparison of eight representative approaches for off-target prediction on 12 publicly available datasets with various imbalanced ratios of positive/negative samples. Most methods showed excellent performance on balanced datasets but have much room for improvement on moderate- and severe-imbalanced datasets. This study provides comprehensive perspectives on CRISPR/Cas9 sgRNA on- and off-target activity prediction and improvement for method development.

Antihypertensive treatment during pregnancy induces long-term changes in gut microbiota and the behaviors of the attention deficit hyperactivity disorder offspring.

The pathogenesis of attention-deficit/hyperactivity disorder (ADHD) has not been fully elucidated. Gestational hypertension could double the probability of ADHD in the offspring, while the initial bacterial communication between the mother and offspring has been associated with psychiatric disorders. Thus, we hypothesize that antihypertensive treatment during pregnancy may abate the impairments in neurodevelopment of the offspring. To test this hypothesis, we chose Captopril and Labetalol, to apply to pregnant spontaneously hypertensive rat (SHR) dams and examined the outcomes in the male offspring. Our data demonstrated that maternal treatment with Captopril and Labetalol had long-lasting changes in gut microbiota and behavioral alterations, including decreased hyperactivity and increased curiosity, spatial learning and memory in the male offspring. Increased diversity and composition were identified, and some ADHD related bacteria were found to have the same change in the gut microbiota of both the dam and offspring after the treatments. LC-MS/MS and immunohistochemistry assays suggested elevated expression of brain derived neurotrophic factor (BDNF) and dopamine in the prefrontal cortex and striatum of offspring exposed to Captopril/ Labetalol, which may account for the improvement of the offspring's psychiatric functions. Therefore, our results support the beneficial long-term effects of the intervention of gestational hypertension in the prevention of ADHD.

Faecalibacterium prausnitzii Attenuates CKD via Butyrate-Renal GPR43 Axis.

BACKGROUND: Despite available clinical management strategies, chronic kidney disease (CKD) is associated with severe morbidity and mortality worldwide, which beckons new solutions. Host-microbial interactions with a depletion of Faecalibacterium prausnitzii in CKD are reported. However, the mechanisms about if and how F prausnitzii can be used as a probiotic to treat CKD remains unknown. METHODS: We evaluated the microbial compositions in 2 independent CKD populations for any potential probiotic. Next, we investigated if supplementation of such probiotic in a mouse CKD model can restore gut-renal homeostasis as monitored by its effects on suppression on renal inflammation, improvement in gut permeability and renal function. Last, we investigated the molecular mechanisms underlying the probiotic-induced beneficial outcomes. RESULTS: We observed significant depletion of Faecalibacterium in the patients with CKD in both Western (n=283) and Eastern populations (n=75). Supplementation of F prausnitzii to CKD mice reduced renal dysfunction, renal inflammation, and lowered the serum levels of various uremic toxins. These are coupled with improved gut microbial ecology and intestinal integrity. Moreover, we demonstrated that the beneficial effects in kidney induced by F prausnitzii-derived butyrate were through the GPR (G protein-coupled receptor)-43. CONCLUSIONS: Using a mouse CKD model, we uncovered a novel beneficial role of F prausnitzii in the restoration of renal function in CKD, which is, at least in part, attributed to the butyrate-mediated GPR-43 signaling in the kidney. Our study provides the necessary foundation to harness the therapeutic potential of F prausnitzii for ameliorating CKD.

The landscape of long noncoding RNA-involved and tumor-specific fusions across various cancers.

The majority of the human genome encodes long noncoding RNA (lncRNA) genes, critical regulators of various cellular processes, which largely outnumber protein-coding genes. However, lncRNA-involved fusions have not been surveyed and characterized yet. Here, we present a systematic study of the lncRNA fusion landscape across cancer types and identify >30 000 high-confidence tumor-specific lncRNA fusions (using 8284 tumor and 6946 normal samples). Fusions positively correlated with DNA damage and cancer stemness and were specifically low in microsatellite instable (MSI)-High or virus-infected tumors. Moreover, fusions distribute differently among cancer molecular subtypes, but with shared enrichment in tumors that are microsatellite stable (MSS), with high somatic copy number alterations (SCNA), and with poor survival. Importantly, we find a potentially new mechanism, mediated by enhancer RNAs (eRNA), which generates secondary fusions that form densely connected fusion networks with many fusion hubs targeted by FDA-approved drugs. Finally, we experimentally validate functions of two tumor-promoting chimeric proteins derived from mRNA-lncRNA fusions, KDM4B-G039927 and EPS15L1-lncOR7C2-1. The EPS15L1 fusion protein may regulate (Gasdermin E) GSDME, critical in pyroptosis and anti-tumor immunity. Our study completes the fusion landscape in cancers, sheds light on fusion mechanisms, and enriches lncRNA functions in tumorigenesis and cancer progression.

SANPolyA: a deep learning method for identifying Poly(A) signals.

MOTIVATION: Polyadenylation plays a regulatory role in transcription. The recognition of polyadenylation signal (PAS) motif sequence is an important step in polyadenylation. In the past few years, some statistical machine learning-based and deep learning-based methods have been proposed for PAS identification. Although these methods predict PAS with success, there is room for their improvement on PAS identification. RESULTS: In this study, we proposed a deep neural network-based computational method, called SANPolyA, for identifying PAS in human and mouse genomes. SANPolyA requires no manually crafted sequence features. We compared our method SANPolyA with several previous PAS identification methods on several PAS benchmark datasets. Our results showed that SANPolyA outperforms the state-of-art methods. SANPolyA also showed good performance on leave-one-motif-out evaluation. AVAILABILITY AND IMPLEMENTATION: https://github.com/yuht4/SANPolyA. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Curcumin ameliorates hypertension via gut-brain communication in spontaneously hypertensive rat.

Gut dysbiosis and dysregulation of gut-brain communication have been identified in hypertensive patients and animal models. Previous studies have shown that probiotic or prebiotic treatments exert positive effects on the pathophysiology of hypertension. This study aimed to examine the hypothesis that the microbiota-gut-brain axis is involved in the antihypertensive effects of curcumin, a potential prebiotic obtained from Curcuma longa. Male 8- to 10-week-old spontaneously hypertensive rats (SHRs) and Wistar Kyoto (WKY) rats were divided into four groups: WKY rats and SHRs treated with vehicle and SHRs treated with curcumin in dosage of 100 or 300 mg/kg/day for 12 weeks. Our results show that the elevated blood pressure of SHRs was markedly decreased in both curcumin-treated groups. Curcumin treatment also altered the gut microbial composition and improved intestinal pathology and integrity. These factors were associated with reduced neuroinflammation and oxidative stress in the hypothalamus paraventricular nucleus (PVN). Moreover, curcumin treatment increased butyrate levels in the plasma, which may be the result of increased butyrate-producing gut microorganisms. In addition, curcumin treatment also activated G protein-coupled receptor 43 (GPR 43) in the PVN. These results indicate that curcumin reshapes the composition of the gut microbiota and ameliorates the dysregulation of the gut-brain communication to induce antihypertensive effects.

Pan-cancer analysis identifies telomerase-associated signatures and cancer subtypes.

BACKGROUND: Cancer cells become immortalized through telomere maintenance mechanisms, such as telomerase reverse transcriptase (TERT) activation. In addition to maintaining telomere length, TERT activates manifold cell survival signaling pathways. However, telomerase-associated gene signatures in cancer remain elusive. METHODS: We performed a systematic analysis of TERT high (TERThigh) and low (TERTlow) cancers using multidimensional data from The Cancer Genome Atlas (TCGA). Multidimensional data were analyzed by propensity score matching weight algorithm. Coexpression networks were constructed by weight gene coexpression network analysis (WGCNA). Random forest classifiers were generated to identify cancer subtypes. RESULTS: The TERThigh-specific mRNA expression signature is associated with cell cycle-related coexpression modules across cancer types. Experimental screening of hub genes in the cell cycle module suggested TPX2 and EXO1 as potential regulators of telomerase activity and cell survival. MiRNA analysis revealed that the TERThigh-specific miR-17-92 cluster can target biological processes enriched in TERTlow cancer and that its expression is negatively correlated with the tumor/normal telomere length ratio. Intriguingly, TERThigh cancers tend to have mutations in extracellular matrix organization genes and amplify MAPK signaling. By mining the clinical actionable gene database, we uncovered a number of TERThigh-specific somatic mutations, amplifications and high expression genes containing therapeutic targets. Finally, a random forest classifier integrating telomerase-associated multi-omics signatures identifies two cancer subtypes showed profound differences in telomerase activity and patient survival. CONCLUSIONS: In summary, our results depict a telomerase-associated molecular landscape in cancers and provide therapeutic opportunities for cancer treatment.

Aptamer-Based Fluorescent Sensor Array for Multiplexed Detection of Cyanotoxins on a Smartphone.

Developing easy-to-use and miniaturized detectors is essential for in-field monitoring of environmentally hazardous substances, such as the cyanotoxins. We demonstrated a differential fluorescent sensor array made of aptamers and single-stranded DNA (ssDNA) dyes for multiplexed detection and discrimination of four common cyanotoxins with an ordinary smartphone within 5 min of reaction. The assay reagents were preloaded and dried in a microfluidic chip with a long shelf life over 60 days. Upon the addition of analyte solutions, competitive binding of cyanotoxin to the specific aptamer-dye conjugate occurred. A zone-specific and concentration-dependent reduction in the green fluorescence was observed as a result of the aptamer conformation change. The aptasensors are fully optimized by quantification of their dissociation constants, tuning the stoichiometric ratios of reaction mixtures, and implementation of an internal intensity correction step. The fluorescent sensor array allowed for accurate identification and measurement of four important cyanotoxins, including anatoxin-a (ATX), cylindrospermopsin (CYN), nodularin (NOD), and microcystin-LR (MC-LR), in parallel, with the limit of detection (LOD) down to a few nanomolar (<3 nM), which is close to the World Health Organization's guideline for the maximum concentration allowed in drinking water. The smartphone-based sensor platform also showed remarkable chemical specificity against potential interfering agents in water. The performance of the system was tested and validated with real lake water samples that were contaminated with trace levels of individual cyanotoxins as well as binary, ternary, and quaternary mixtures. Finally, a smartphone app interface has been developed for rapid on-site data processing and result display.

Insights into active intragenic enhancers.

Enhancers can regulate gene transcription from afar. Many enhancers are located in genes. Although the regulatory roles of several individual intragenic enhancers have been elaborated, a genome-wide insight into intragenic enhancers remains to be elucidated. We found that active intragenic enhancers have a preference for being located in expressed genes. Unlike intergenic enhancers, active intragenic enhancers are enriched of H3K79me2 epigenetic signal, and depleted of variant histone H2A.Z. Moreover, eRNAs of active intragenic enhancers show lower degradation rates than those of the other enhancers. Our findings will have implications in understanding functions of intragenic enhancers.

pgRNAFinder: a web-based tool to design distance independent paired-gRNA.

SUMMARY: The CRISPR/Cas System has been shown to be an efficient and accurate genome-editing technique. There exist a number of tools to design the guide RNA sequences and predict potential off-target sites. However, most of the existing computational tools on gRNA design are restricted to small deletions. To address this issue, we present pgRNAFinder, with an easy-to-use web interface, which enables researchers to design single or distance-free paired-gRNA sequences. The web interface of pgRNAFinder contains both gRNA search and scoring system. After users input query sequences, it searches gRNA by 3' protospacer-adjacent motif (PAM), and possible off-targets, and scores the conservation of the deleted sequences rapidly. Filters can be applied to identify high-quality CRISPR sites. PgRNAFinder offers gRNA design functionality for 8 vertebrate genomes. Furthermore, to keep pgRNAFinder open, extensible to any organism, we provide the source package for local use. AVAILABILITY AND IMPLEMENTATION: The pgRNAFinder is freely available at http://songyanglab.sysu.edu.cn/wangwebs/pgRNAFinder/, and the source code and user manual can be obtained from https://github.com/xiexiaowei/pgRNAFinder. CONTACT: songyang@bcm.edu or daizhim@mail.sysu.edu.cn. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Infusion of leukocytes from HLA haplo-identical familial donors as an adjuvant in the HLH-2004 protocol to treat the virus-associated adult hemophagocytic lymphohistiocytosis: a retrospective study of 26 patients.

Adult hemophagocytic lymphohistiocytosis (HLH) is a fatal disease with poor survival and a limited role of drug therapies. To help to recognize virus and enhance survival, we infused leukocytes derived from human leukocyte antigen (HLA) haplo-identical familial donors to patients. We retrospectively investigated 26 adult virus-associated hemophagocytic syndrome (VAHS) patients' medical records from 2006-2017. Eleven of the 26 patients accepted relatives' derived leukocytes infusions in addition to drug therapies recommended in the HLH-2004 protocol. The leukocyte doses ranged from 0.75 to 3.30×108 per kilogram of body weight. The other 15 patients accepted immunosuppressive and supportive therapies referred to in the HLH-2004 protocol. We compared the treatment outcomes of the two groups of patients. Patients in the cell infusion group had a lower viral load (P = 0.023) and better laboratory results and prolonged overall survival (60.44 vs. 20.18 weeks, P = 0.047). A factor that might relate to overall survival is platelet count (P = 0.032), except for the leukocyte infusions (P = 0.012). For patients without acceptable donors, infusions of leukocytes from HLA haplo-identical familial donors could be a feasible treatment to prolong overall survival as an adjuvant to drug therapies.

Horse versus rabbit antithymocyte globulin in immunosuppressive therapy of treatment-naïve aplastic anemia: a systematic review and meta-analysis.

The first-line formulation of antithymocyte globulin (ATG) remains unknown. We aimed to systematically review evidence to compare the efficacy and safety profiles of different ATGs. We did a systematic review and meta-analysis of randomized controlled trials (RCTs) and cohort controlled studies comparing horse and rabbit ATG in immunosuppressive therapy of treatment-naïve aplastic anemia. We searched The Cochrane Library, PubMed, EMBASE, ClinicalTrials.gov , and conference proceedings of American Society of Hematology and European Society for Blood and Marrow Transplantation annual meetings. The outcomes were 3-, 6-, and 12-month response; early mortality; relapse; and evolution. We pooled hazard ratios for relapse and odds ratios (ORs) for other outcomes using fixed-effect or random-effect models based on the heterogeneity. This study was registered with PROSPERO, number CRD42016036945. We included 1636 participants from three RCTs and 11 cohort controlled studies. Allocation to horse ATG increased 6-month response events by 86% compared with rabbit ATG. The benefit of horse ATG was mainly driven by increase in studies with non-Asian (OR 95% CI = 2.39 (1.54-3.69), p < 0.0001) and good partial response criterion (OR 95% CI = 2.73 (1.53-4.89), p = 0.0007). The early mortality and evolution were similar between groups. Compared with rabbit ATG, horse ATG had superior remission by 6 months and equivalent safety profiles in patients with treatment-naïve AA. Evidence for further responses beyond 6 to 12 months was limited.

Lack of association between cytotoxic T-lymphocyte antigen-4 gene polymorphisms and lymphoid malignancy risk: evidence from a meta-analysis.

Cytotoxic T-lymphocyte antigen-4 (CTLA-4) polymorphisms have been associated with susceptibility to lymphoid malignancies. However, results from the published single studies are inconsistent. Therefore, the present meta-analysis was conducted to get a more accurate estimation of the relationship between CTLA-4 gene polymorphisms and the lymphoid malignancy risk. We identified nine independent studies accounting for 3090 subjects up to January 30, 2016. Summary odds ratios (OR) and 95 % confidence intervals (CI) were used to evaluate the risk of lymphoid malignancies. Overall, no significant association was found between +49A/G (rs231775), -318C/T (rs5742909), and +6230A/G (rs3087243) CTLA-4 gene polymorphisms and lymphoid malignancies. Furthermore, ethnicity (Asian and Caucasian) and histopathology subgroup analyses (non-Hodgkin's lymphoma) also failed to detect an association between the studied polymorphisms and lymphoid malignancy risk. Our study shows that common CTLA-4 gene polymorphisms may not contribute to lymphoid malignancy susceptibility based on the current evidence.

Genome-wide analysis of transcription factor binding sites and their characteristic DNA structures.

BACKGROUND: Transcription factors (TF) regulate gene expression by binding DNA regulatory regions. Transcription factor binding sites (TFBSs) are conserved not only in primary DNA sequences but also in DNA structures. However, the global relationship between TFs and their preferred DNA structures remains to be elucidated. RESULTS: In this paper, we have developed a computational method to generate a genome-wide landscape of TFs and their characteristic binding DNA structures in Saccharomyces cerevisiae. We revealed DNA structural features for different TFs. The structural conservation shows positional preference in TFBSs. Structural levels of DNA sequences are correlated with TF-DNA binding affinities. CONCLUSIONS: We provided the genome-wide correspondences of TFs to DNA structures. Our findings will have implications in understanding TF regulatory mechanisms.

Neighboring genes show interchromosomal colocalization after their separation.

The order of genes on eukaryotic chromosomes is nonrandom. Some neighboring genes show order conservation among species, while some neighboring genes separate during evolution. Here, we investigated whether neighboring genes show interactions after their separation. We found that neighboring gene pairs tend to show interchromosomal colocalization (i.e., nuclear colocalization) in the species in which they separate. These nuclear colocalized separated neighboring gene pairs 1) show neighborhood conservation in more species, 2) tend to be regulated by the same transcription factor, and 3) tend to be regulated by the same histone modification. These results suggest a mechanism by which neighboring genes could retain nuclear proximity after their separation.

The pattern and evolution of looped gene bendability.

Gene looping, defined as the physical interaction between the promoter and terminator regions of a RNA polymerase II-transcribed gene, is widespread in yeast and mammalian cells. Gene looping has been shown to play important roles in transcription. Gene-loop formation is dependent on regulatory proteins localized at the 5' and 3' ends of genes, such as TFIIB. However, whether other factors contribute to gene looping remains to be elucidated. Here, we investigated the contribution of intrinsic DNA and chromatin structures to gene looping. We found that Saccharomyces cerevisiae looped genes show high DNA bendability around middle and 3/4 regions in open reading frames (ORFs). This bendability pattern is conserved between yeast species, whereas the position of bendability peak varies substantially among species. Looped genes in human cells also show high DNA bendability. Nucleosome positioning around looped ORF middle regions is unstable. We also present evidence indicating that this unstable nucleosome positioning is involved in gene looping. These results suggest a mechanism by which DNA bendability and unstable nucleosome positioning could assist in the formation of gene loops.

Nuclear colocalization of transcription factor target genes strengthens coregulation in yeast.

Eukaryotic chromosomes are not randomly distributed in the interphase nucleus, but instead occupy distinct territories. Nonetheless, the genome-wide relationships of gene regulation to gene nuclear location remain poorly understood in yeast. In the three-dimensional view of gene regulation, we found that a considerable number of transcription factors (TFs) regulate genes that are colocalized in the nucleus. Colocalized TF target genes are more strongly coregulated compared with the other TF target genes. Target genes of chromatin regulators are also colocalized. These results demonstrate that colocalization of coregulated genes is a common process, and three-dimensional gene positioning is an important part of gene regulation. Our findings will have implications in understanding nuclear architecture and function.