Search

Abstract Single-cell atlases pose daunting computational challenges pertaining to the integration of spatial and temporal information and the visualization of trajectories across large atlases. We introduce StaVia, a computational framework that synergizes multi-faceted single-cell data with higher-order random walks that leverage the memory of cells’ past states, fused with a cartographic Atlas View that offers intuitive graph visualization. This spatially aware cartography captures relationships between cell populations based on their spatial location as well as their gene expression and developmental stage. We demonstrate this using zebrafish gastrulation data, underscoring its potential to dissect complex biological landscapes in both spatial and temporal contexts.

To explore restoration of ovarian function using epigenetically-related, induced pluripotent stem cells (iPSCs), we functionally evaluated the epigenetic memory of novel iPSC lines, derived from mouse and human ovarian granulosa cells (GCs) using c-Myc, Klf4, Sox2 and Oct4 retroviral vectors. The stem cell identity of the mouse and human GC-derived iPSCs (mGriPSCs, hGriPSCs) was verified by demonstrating embryonic stem cell (ESC) antigen expression using immunocytochemistry and RT-PCR analysis, as well as formation of embryoid bodies (EBs) and teratomas that are capable of differentiating into cells from all three germ layers. GriPSCs' gene expression profiles associate more closely with those of ESCs than of the originating GCs as demonstrated by genome-wide analysis of mRNA and microRNA. A comparative analysis of EBs generated from three different mouse cell lines (mGriPSCs; fibroblast-derived iPSC, mFiPSCs; G4 embryonic stem cells, G4 mESCs) revealed that differentiated mGriPSC-EBs synthesize 10-fold more estradiol (E2) than either differentiated FiPSC- or mESC-EBs under identical culture conditions. By contrast, mESC-EBs primarily synthesize progesterone (P4) and FiPSC-EBs produce neither E2 nor P4. Differentiated mGriPSC-EBs also express ovarian markers (AMHR, FSHR, Cyp19a1, ER and Inha) as well as markers of early gametogenesis (Mvh, Dazl, Gdf9, Boule and Zp1) more frequently than EBs of the other cell lines. These results provide evidence of preferential homotypic differentiation of mGriPSCs into ovarian cell types. Collectively, our data support the hypothesis that generating iPSCs from the desired tissue type may prove advantageous due to the iPSCs' epigenetic memory.

Gene regulatory networks (GRNs) play a central role in systems biology, especially in the study of mammalian organ development. One key question remains largely unanswered: Is it possible to infer mammalian causal GRNs using observable gene co-expression patterns alone? We assembled two mouse GRN datasets (embryonic tooth and heart) and matching microarray gene expression profiles to systematically investigate the difficulties of mammalian causal GRN inference. The GRNs were assembled based on > 2,000 pieces of experimental genetic perturbation evidence from manually reading > 150 primary research articles. Each piece of perturbation evidence records the qualitative change of the expression of one gene following knock-down or over-expression of another gene. Our data have thorough annotation of tissue types and embryonic stages, as well as the type of regulation (activation, inhibition and no effect), which uniquely allows us to estimate both sensitivity and specificity of the inference of tissue specific causal GRN edges. Using these unprecedented datasets, we found that gene co-expression does not reliably distinguish true positive from false positive interactions, making inference of GRN in mammalian development very difficult. Nonetheless, if we have expression profiling data from genetic or molecular perturbation experiments, such as gene knock-out or signalling stimulation, it is possible to use the set of differentially expressed genes to recover causal regulatory relationships with good sensitivity and specificity. Our result supports the importance of using perturbation experimental data in causal network reconstruction. Furthermore, we showed that causal gene regulatory relationship can be highly cell type or developmental stage specific, suggesting the importance of employing expression profiles from homogeneous cell populations. This study provides essential datasets and empirical evidence to guide the development of new GRN inference methods for mammalian organ development.

The diagnosis of Systemic Lupus Erythematosus (SLE) is challenging due to its heterogeneous clinical presentation and the lack of robust biomarkers to distinguish it from other autoimmune diseases. Further, currently used laboratory tests do not readily distinguish active and inactive disease. Several groups have attempted to apply emerging high throughput profiling technologies to diagnose and monitor SLE. Despite showing promise, many are expensive and technically challenging for routine clinical use. The goal of this work is to develop a better diagnostic and monitoring tool for SLE. We report a highly customisable antibody microarray that consists of a duplicate arrangement of 82 antibodies directed against surface antigens on peripheral blood mononuclear cells (PMBCs). This high-throughput array was used to profile SLE patients (n = 60) with varying disease activity, compared to healthy controls (n = 24), patients with rheumatoid arthritis (n = 25), and other autoimmune diseases (n = 28). We used a computational algorithm to calculate a score from the entire microarray profile and correlated it with SLE disease activity. Our results demonstrate that leukocyte-capture microarray profiles can readily distinguish active SLE patients from healthy controls (AUROC = 0.84). When combined with the standard laboratory tests (serum anti-dsDNA, complements C3 and C4), the microarrays provide significantly increased discrimination. The antibody microarrays can be enhanced by the addition of other markers for potential application to the diagnosis and stratification of SLE, paving the way for the customised and accurate diagnosis and monitoring of SLE.

The Drosophila MSL complex mediates dosage compensation by increasing transcription of the single X chromosome in males approximately two-fold. This is accomplished through recognition of the X chromosome and subsequent acetylation of histone H4K16 on X-linked genes. Initial binding to the X is thought to occur at "entry sites" that contain a consensus sequence motif ("MSL recognition element" or MRE). However, this motif is only ∼2 fold enriched on X, and only a fraction of the motifs on X are initially targeted. Here we ask whether chromatin context could distinguish between utilized and non-utilized copies of the motif, by comparing their relative enrichment for histone modifications and chromosomal proteins mapped in the modENCODE project. Through a comparative analysis of the chromatin features in male S2 cells (which contain MSL complex) and female Kc cells (which lack the complex), we find that the presence of active chromatin modifications, together with an elevated local GC content in the surrounding sequences, has strong predictive value for functional MSL entry sites, independent of MSL binding. We tested these sites for function in Kc cells by RNAi knockdown of Sxl, resulting in induction of MSL complex. We show that ectopic MSL expression in Kc cells leads to H4K16 acetylation around these sites and a relative increase in X chromosome transcription. Collectively, our results support a model in which a pre-existing active chromatin environment, coincident with H3K36me3, contributes to MSL entry site selection. The consequences of MSL targeting of the male X chromosome include increase in nucleosome lability, enrichment for H4K16 acetylation and JIL-1 kinase, and depletion of linker histone H1 on active X-linked genes. Our analysis can serve as a model for identifying chromatin and local sequence features that may contribute to selection of functional protein binding sites in the genome.

To elucidate the transcriptional regulation of Bmp4 expression during organogenesis, we used phylogenetic footprinting and transgenic reporter analyses to identify Bmp4 cis-regulatory modules (CRMs). These analyses identified a regulatory region located ∼46 kb upstream of the mouse Bmp4 transcription start site that had previously been shown to direct expression in lateral plate mesoderm. We refined this regulatory region to a 396-bp minimal enhancer, and show that it recapitulates features of endogenous Bmp4 expression in developing mandibular arch ectoderm and incisor epithelium during the initiation-stage of tooth development. In addition, this enhancer directs expression in the apical ectodermal ridge (AER) of the developing limb and in anterior and posterior limb mesenchyme. Transcript profiling of E11.5 mouse incisor dental lamina, together with protein binding microarray (PBM) analyses, allowed identification of a conserved DNA binding motif in the Bmp4 enhancer for Pitx homeoproteins, which are also expressed in the developing mandibular and incisor epithelium. In vitro electrophoretic mobility shift assays (EMSA) and in vivo transgenic reporter mutational analyses revealed that this site supports Pitx binding and that the site is necessary to recapitulate aspects of endogenous Bmp4 expression in developing craniofacial and limb tissues. Finally, Pitx2 chromatin immunoprecipitation (ChIP) demonstrated direct binding of Pitx2 to this Bmp4 enhancer site in a dental epithelial cell line. These results establish a direct molecular regulatory link between Pitx family members and Bmp4 gene expression in developing incisor epithelium.

Abstract Differential composition analysis — the identification of cell types that have statistically significant changes in abundance between multiple experimental conditions — is one of the most common tasks in single cell omic data analysis. However, it remains challenging to perform differential composition analysis in the presence of flexible experimental designs and uncertainty in cell type assignment. Here, we introduce a statistical model and an open source R package, DCATS, for differential composition analysis based on a beta-binomial regression framework that addresses these challenges. Our empirical evaluation shows that DCATS consistently maintains high sensitivity and specificity compared to state-of-the-art methods.

Abstract Transient gut microbiota alterations have been reported after antibiotic therapy for Helicobacter pylori. However, alteration in the gut virome after H. pylori eradication remains uncertain. Here, we apply metagenomic sequencing to fecal samples of 44 H. pylori-infected patients at baseline, 6-week (N = 44), and 6-month (N = 33) after treatment. Following H. pylori eradication, we discover contraction of the gut virome diversity, separation of virome community with increased community difference, and shifting towards a higher proportion of core virus. While the gut microbiota is altered at 6-week and restored at 6-month, the virome community shows contraction till 6-month after the treatment with enhanced phage-bacteria interactions at 6-week. Multiple courses of antibiotic treatments further lead to lower virus community diversity when compared with treatment naive patients. Our results demonstrate that H. pylori eradication therapies not only result in transient alteration in gut microbiota but also significantly alter the previously less known gut virome community.

Read alignment is an important step in RNA-seq analysis as the result of alignment forms the basis for downstream analyses. However, recent studies have shown that published alignment tools have variable mapping sensitivity and do not necessarily align all the reads which should have been aligned, a problem we termed as the false-negative non-alignment problem. Here we present Scavenger, a python-based bioinformatics pipeline for recovering unaligned reads using a novel mechanism in which a putative alignment location is discovered based on sequence similarity between aligned and unaligned reads. We showed that Scavenger could recover unaligned reads in a range of simulated and real RNA-seq datasets, including single-cell RNA-seq data. We found that recovered reads tend to contain more genetic variants with respect to the reference genome compared to previously aligned reads, indicating that divergence between personal and reference genomes plays a role in the false-negative non-alignment problem. Even when the number of recovered reads is relatively small compared to the total number of reads, the addition of these recovered reads can impact downstream analyses, especially in terms of estimating the expression and differential expression of lowly expressed genes, such as pseudogenes.

Read alignment is an important step in RNA-seq analysis as the result of alignment forms the basis for downstream analyses. However, recent studies have shown that published alignment tools have variable mapping sensitivity and do not necessarily align all the reads which should have been aligned, a problem we termed as the false-negative non-alignment problem. Here we present Scavenger, a python-based bioinformatics pipeline for recovering unaligned reads using a novel mechanism in which a putative alignment location is discovered based on sequence similarity between aligned and unaligned reads. We showed that Scavenger could recover unaligned reads in a range of simulated and real RNA-seq datasets, including single-cell RNA-seq data. We found that recovered reads tend to contain more genetic variants with respect to the reference genome compared to previously aligned reads, indicating that divergence between personal and reference genomes plays a role in the false-negative non-alignment problem. Even when the number of recovered reads is relatively small compared to the total number of reads, the addition of these recovered reads can impact downstream analyses, especially in terms of estimating the expression and differential expression of lowly expressed genes, such as pseudogenes.

Abstract Most existing dimensionality reduction and clustering packages for single-cell RNA-seq (scRNA-seq) data deal with dropouts by heavy modeling and computational machinery. Here, we introduce CIDR (Clustering through Imputation and Dimensionality Reduction), an ultrafast algorithm that uses a novel yet very simple implicit imputation approach to alleviate the impact of dropouts in scRNA-seq data in a principled manner. Using a range of simulated and real data, we show that CIDR improves the standard principal component analysis and outperforms the state-of-the-art methods, namely t-SNE, ZIFA, and RaceID, in terms of clustering accuracy. CIDR typically completes within seconds when processing a data set of hundreds of cells and minutes for a data set of thousands of cells. CIDR can be downloaded at https://github.com/VCCRI/CIDR .

To elucidate the transcriptional regulation of Bmp4 expression during organogenesis, we used phylogenetic footprinting and transgenic reporter analyses to identify Bmp4 cis-regulatory modules (CRMs). These analyses identified a regulatory region located ,46 kb upstream of the mouse Bmp4 transcription start site that had previously been shown to direct expression in lateral plate mesoderm. We refined this regulatory region to a 396-bp minimal enhancer, and show that it recapitulates features of endogenous Bmp4 expression in developing mandibular arch ectoderm and incisor epithelium during the initiation-stage of tooth development. In addition, this enhancer directs expression in the apical ectodermal ridge (AER) of the developing limb and in anterior and posterior limb mesenchyme. Transcript profiling of E11.5 mouse incisor dental lamina, together with protein binding microarray (PBM) analyses, allowed identification of a conserved DNA binding motif in the Bmp4 enhancer for Pitx homeoproteins, which are also expressed in the developing mandibular and incisor epithelium. In vitro electrophoretic mobility shift assays (EMSA) and in vivo transgenic reporter mutational analyses revealed that this site supports Pitx binding and that the site is necessary to recapitulate aspects of endogenous Bmp4 expression in developing craniofacial and limb tissues. Finally, Pitx2 chromatin immunoprecipitation (ChIP) demonstrated direct binding of Pitx2 to this Bmp4 enhancer site in a dental epithelial cell line. These results establish a direct molecular regulatory link between Pitx family members and Bmp4 gene expression in developing incisor epithelium. [ABSTRACT FROM AUTHOR]

Screening combinatorial mutants is too massive for wet-lab experiment alone. Here the authors present a machine learning-coupled combinatorial mutagenesis approach to vastly reduce experimental burden for engineering Cas9 genome editing enzymes.

Mitochondrial variants are informative endogenous barcodes for clonal substructure. Here, the authors developed a computational method MQuad to effectively detect these clonal informed mtDNA variants from single-cell RNA, DNA or ATAC sequencing data.

Scalable trajectory inference for multi-omic single cell datasets is challenging in terms of capturing non-tree complex topologies. Here the authors present a method, VIA, that scales to millions of cells across multiple omic modalities using lazy-teleporting random walks.

Abstract Triaging and prioritising patients for RT-PCR test had been essential in the management of COVID-19 in resource-scarce countries. In this study, we applied machine learning (ML) to the task of detection of SARS-CoV-2 infection using basic laboratory markers. We performed the statistical analysis and trained an ML model on a retrospective cohort of 5148 patients from 24 hospitals in Hong Kong to classify COVID-19 and other aetiology of pneumonia. We validated the model on three temporal validation sets from different waves of infection in Hong Kong. For predicting SARS-CoV-2 infection, the ML model achieved high AUCs and specificity but low sensitivity in all three validation sets (AUC: 89.9–95.8%; Sensitivity: 55.5–77.8%; Specificity: 91.5–98.3%). When used in adjunction with radiologist interpretations of chest radiographs, the sensitivity was over 90% while keeping moderate specificity. Our study showed that machine learning model based on readily available laboratory markers could achieve high accuracy in predicting SARS-CoV-2 infection.

Abstract Background The human gut microbiome plays a critical role in the carcinogenesis of colorectal cancer (CRC). However, a comprehensive analysis of the interaction between the host and microbiome is still lacking. Results We found correlations between the change in abundance of microbial taxa, butyrate-related colonic metabolites, and methylation-associated host gene expression in colonic tumour mucosa tissues compared with the adjacent normal mucosa tissues. The increase of genus Fusobacterium abundance was correlated with a decrease in the level of 4-hydroxybutyric acid (4-HB) and expression of immune-related peptidase inhibitor 16 (PI16), Fc Receptor Like A (FCRLA) and Lymphocyte Specific Protein 1 (LSP1). The decrease in the abundance of another potentially 4-HB-associated genus, Prevotella 2, was also found to be correlated with the down-regulated expression of metallothionein 1 M (MT1M). Additionally, the increase of glutamic acid-related family Halomonadaceae was correlated with the decreased expression of reelin (RELN). The decreased abundance of genus Paeniclostridium and genus Enterococcus were correlated with increased lactic acid level, and were also linked to the expression change of Phospholipase C Beta 1 (PLCB1) and Immunoglobulin Superfamily Member 9 (IGSF9) respectively. Interestingly, 4-HB, glutamic acid and lactic acid are all butyrate precursors, which may modify gene expression by epigenetic regulation such as DNA methylation. Conclusions Our study identified associations between previously reported CRC-related microbial taxa, butyrate-related metabolites and DNA methylation-associated gene expression in tumour and normal colonic mucosa tissues from CRC patients, which uncovered a possible mechanism of the role of microbiome in the carcinogenesis of CRC. In addition, these findings offer insight into potential new biomarkers, therapeutic and/or prevention strategies for CRC.

Abstract Background Read alignment and transcript assembly are the core of RNA-seq analysis for transcript isoform discovery. Nonetheless, current tools are not designed to be scalable for analysis of full-length bulk or single cell RNA-seq (scRNA-seq) data. The previous version of our cloud-based tool Falco only focuses on RNA-seq read counting, but does not allow for more flexible steps such as alignment and read assembly. Results The Falco framework can harness the parallel and distributed computing environment in modern cloud platforms to accelerate read alignment and transcript assembly of full-length bulk RNA-seq and scRNA-seq data. There are two new modes in Falco: alignment-only and transcript assembly. In the alignment-only mode, Falco can speed up the alignment process by 2.5–16.4x based on two public scRNA-seq datasets when compared to alignment on a highly optimised standalone computer. Furthermore, it also provides a 10x average speed-up compared to alignment using published cloud-enabled tool for read alignment, Rail-RNA. In the transcript assembly mode, Falco can speed up the transcript assembly process by 1.7–16.5x compared to performing transcript assembly on a highly optimised computer. Conclusion Falco is a significantly updated open source big data processing framework that enables scalable and accelerated alignment and assembly of full-length scRNA-seq data on the cloud. The source code can be found at https://github.com/VCCRI/Falco.

Abstract Background The application of next-generation sequencing in cancer has revealed the genomic landscape of many tumour types and is nowadays routinely used in research and clinical settings. Multiple algorithms have been developed to detect somatic variation from sequencing data using either paired tumour-blood or tumour-only samples. Most of these methods have been developed and evaluated for the identification of somatic variation using Illumina sequencing datasets of moderate coverage. However, a comprehensive evaluation of somatic variant detection algorithms on Ion Torrent targeted deep sequencing data has not been performed. Methods We have applied three somatic detection algorithms, Torrent Variant Caller, MuTect2 and VarScan2, on a large cohort of ovarian cancer patients comprising of 208 paired tumour-blood samples and 253 tumour-only samples sequenced deeply on Ion Torrent Proton platform across 330 amplicons. Subsequently, the concordance and performance of the three somatic variant callers were assessed. Results We have observed low concordance across the algorithms with only 0.5% of SNV and 0.02% of INDEL calls in common across all three methods. The intersection of all methods showed better performance when assessed using correlation with known mutational signatures, overlap with COSMIC variation and by examining the variant characteristics. The Torrent Variant Caller also performed well with the advantage of not eliminating a high number of variants that could lead to high type II error. Conclusions Our results suggest that caution should be taken when applying state-of-the-art somatic variant algorithms to Ion Torrent targeted deep sequencing data. Better quality control procedures and strategies that combine results from multiple methods should ensure that higher accuracy is achieved. This is essential to ensure that results from bioinformatics pipelines using Ion Torrent deep sequencing can be robustly applied in cancer research and in the clinic.

Abstract DamID, in which a protein of interest is fused to Dam methylase, enables mapping of protein-DNA binding through readout of adenine methylation in genomic DNA. DamID offers a compelling alternative to chromatin immunoprecipitation sequencing (ChIP-Seq), particularly in cases where cell number or antibody availability is limiting. This comes at a cost, however, of high non-specific signal and a lowered spatial resolution of several kb, limiting its application to transcription factor-DNA binding. Here we show that mutations in Dam, when fused to the transcription factor Tcf7l2, greatly reduce non-specific methylation. Combined with a simplified DamID sequencing protocol, we find that these Dam mutants allow for accurate detection of transcription factor binding at a sensitivity and spatial resolution closely matching that seen in ChIP-seq.

The cell elasticity has attracted extensive research interests since it not only provides new insights into cell biology but also is an emerging mechanical marker for the diagnosis of some diseases. This paper reports the method for the precise measurement of mechanical properties of single cells deformed to a large extent using a novel microfluidic system integrated with a pressure feedback system and small particle separation unit. The particle separation system was employed to avoid the blockage of the cell deformation channel to enhance the measurement throughput. This system is of remarkable application potential in the precise evaluation of cell mechanical properties. In brief, this paper reports: • The manufacturing of the chip using standard soft lithography; • The methods to deform single cells in a microchannel and measure the relevant pressure drop using a pressure sensor connecting to the microfluidic chip; • Calculation of the mechanical properties including stiffness and fluidity of each cell based on a power-law rheology model describing the viscoelastic behaviors of cells; • Automatic and real-time measurement of the mechanical properties using video processing software.

Abstract Background It has been observed that many transcription factors (TFs) can bind to different genomic loci depending on the cell type in which a TF is expressed in, even though the individual TF usually binds to the same core motif in different cell types. How a TF can bind to the genome in such a highly cell-type specific manner, is a critical research question. One hypothesis is that a TF requires co-binding of different TFs in different cell types. If this is the case, it may be possible to observe different combinations of TF motifs – a motif grammar – located at the TF binding sites in different cell types. In this study, we develop a bioinformatics method to systematically identify DNA motifs in TF binding sites across multiple cell types based on published ChIP-seq data, and address two questions: (1) can we build a machine learning classifier to predict cell-type specificity based on motif combinations alone, and (2) can we extract meaningful cell-type specific motif grammars from this classifier model. Results We present a Random Forest (RF) based approach to build a multi-class classifier to predict the cell-type specificity of a TF binding site given its motif content. We applied this RF classifier to two published ChIP-seq datasets of TF (TCF7L2 and MAX) across multiple cell types. Using cross-validation, we show that motif combinations alone are indeed predictive of cell types. Furthermore, we present a rule mining approach to extract the most discriminatory rules in the RF classifier, thus allowing us to discover the underlying cell-type specific motif grammar. Conclusions Our bioinformatics analysis supports the hypothesis that combinatorial TF motif patterns are cell-type specific.