Search

Background Vestibular schwannoma is the third most common benign intracranial tumor that can occur sporadically or be associated with neurofibromatosis type 2 (neurofibromatosis type 2 vestibular schwannoma [NF2-VS]). The aim of this study is to provide a comprehensive bioinformatic analysis of methylated-differentially expressed genes (MDEGs) in NF2-VS. Methods Transcriptional sequencing datasets (GSE141801 and GSE108524) and gene methylation microarrays (GSE56598) from the Gene Expression Omnibus database were used to identify and analyze MDEGs in NF2-VS. A protein–protein interaction (PPI) network was built, and the hub genes and modules were identified. Finally, potential pharmacotherapy targeting MDEGs were extracted for NF2-VS. Results A total of 57 hypermethylation–low expression genes and 88 hypomethylation–high expression genes were identified. Pathways associated with aberrantly MDEGs included P13K-AKT, MAPK, and Ras, which were also involved in NF2-VS. Six hub genes (EGFR, CCND1, CD53, CSF1R, PLAU, and FGFR1) were identified from the PPI network. Modification of the aforementioned genes altered cell-to-cell communication, response to stimulus, cellular regulation, and membrane and protein bindings. Thirty drugs targeting these pathways were selected based on the hub genes. Conclusions Analysis of MDEGs may enrich the understanding of the molecular mechanisms of NF2-VS pathogenesis and lay the groundwork for potential biomarkers and therapeutic targets for NF2-VS.

Background While most differential coexpression (DC) methods are bound to quantify a single correlation value for a gene pair across multiple samples, a newly devised approach under the name Correlation by Individual Level Product (CILP) revolutionarily projects the summary correlation value to individual product correlation values for separate samples. CILP greatly widened DC analysis opportunities by allowing integration of non-compromised statistical methods. Methods Here, we performed a study to verify our hypothesis that conditional relationships, i.e., gene pairs of remarkable differential coexpression, may be sought as quantitative prognostic markers for human cancers. Alongside the seeking of prognostic gene links in a pan-cancer setting, we also examined whether a trend of global expression correlation loss appeared in a wide panel of cancer types and revisited the controversial subject of mutual relationship between the DE approach and the DC approach. Results By integrating CILP with classical univariate survival analysis, we identified up to 244 conditional gene links as potential prognostic markers in five cancer types. In particular, five prognostic gene links for kidney renal papillary cell carcinoma tended to condense around cancer gene ESPL1, and the transcriptional synchrony between ESPL1 and PTTG1 tended to be elevated in patients of adverse prognosis. In addition, we extended the observation of global trend of correlation loss in more than ten cancer types and empirically proved DC analysis results were independent of gene differential expression in five cancer types. Conclusions Combining the power of CILP and the classical survival analysis, we successfully fetched conditional transcriptional relationships that conferred prognosis power for five cancer types. Despite a general trend of global correlation loss in tumor transcriptomes, most of these prognosis conditional links demonstrated stronger expression correlation in tumors, and their stronger coexpression was associated with poor survival.

Background We previously identified differentially expressed genes on the basis of false discovery rate adjusted P value using empirical Bayes moderated tests. However, that approach yielded a subset of differentially expressed genes without accounting for redundancy between the selected genes. Methods This study is a secondary analysis of a case–control study of the effect of antiretroviral therapy on apoptosis pathway genes comprising of 16 cases (HIV infected with mitochondrial toxicity) and 16 controls (uninfected). We applied the maximum relevance minimum redundancy (mRMR) algorithm on the genes that were differentially expressed between the cases and controls. The mRMR algorithm iteratively selects features (genes) that are maximally relevant for class prediction and minimally redundant. We implemented several machine learning classifiers and tested the prediction accuracy of the two mRMR genes. We next used network analysis to estimate and visualize the association among the differentially expressed genes. We employed Markov Random Field or undirected network models to identify gene networks related to mitochondrial toxicity. The Spinglass model was used to identify clusters of gene communities. Results The mRMR algorithm ranked DFFA and TNFRSF1A, two of the upregulated proapoptotic genes, on the top. The overall prediction accuracy was 86%, the two mRMR genes correctly classified 86% of the participants into their respective groups. The estimated network models showed different patterns of gene networks. In the network of the cases, FASLG was the most central gene. However, instead of FASLG, ABL1 and LTBR had the highest centrality in controls. Conclusion The mRMR algorithm and network analysis revealed a new correlation of genes associated with mitochondrial toxicity.

Background Recurrent spontaneous abortion (RSA) is one of the common complication of pregnancy, bringing heavy burden to the patients and their families. The study aimed to explore the lncRNA-miRNA-mRNA network associated with recurrent spontaneous abortion. Methods By transcriptome sequencing, we detected differences in lncRNA, miRNA and mRNA expression in villus tissue samples collected from 3 patients with RSA and 3 normal abortion patients. Differentially expressed lncRNAs, miRNAs and genes (DELs, DEMs and DEGs, respectively) were identified, and Geno Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used to determine the functions of DELs and DEGs, which were analysed by Fisher’s test. We also observed the regulatory relationships between miRNA-mRNA and lncRNA-miRNA by Cytoscape 3.6.1. Results The results showed that 1008 DELs (523 upregulated and 485 downregulated), 475 DEGs (201 upregulated and 274 downregulated) and 37 DEMs (15 upregulated and 22 downregulated) were identified. And we also constructed a novel lncRNA-related ceRNA network containing 31 lncRNAs, 1 miRNA (hsa-miR-210-5p) and 3 genes (NTNG2, GRIA1 and AQP1). Conclusions lncRNA-related ceRNA network containing 31 lncRNAs, 1 miRNA (hsa-miR-210-5p) and 3 mRNAs (NTNG2, GRIA1 and AQP1) was constructed. The results may provide a basic theory for elucidating the mechanism underlying RSA.

Hui Yu,1 Weikang Guo,2 Yunduo Liu,2 Yaoxian Wang2 1Cardiopulmonary Function Department, Harbin Medical University Cancer Hospital, Harbin, Heilongjiang, 150081, People’s Republic of China; 2Gynecological Department, Harbin Medical University Cancer Hospital, Harbin, Heilongjiang, 150081, People’s Republic of ChinaCorrespondence: Yaoxian WangGynecological Department, Harbin Medical University Cancer Hospital, No. 150 Haping Road, Nangang District, Harbin, Heilongjiang, 150081, People’s Republic of ChinaTel +86-0451-86298578Email wyxxs012@126.comPurpose: The variation in inflammation in chronic obstructive pulmonary disease (COPD) between individuals is genetically determined. This study aimed to identify gene signatures of COPD through bioinformatics analysis based on multiple gene sets and explore their immune characteristics and transcriptional regulation mechanisms.Methods: Data from four microarrays were downloaded from the Gene Expression Omnibus database to screen differentially expressed genes (DEGs) between COPD patients and controls. Weighted gene co-expression network analysis was applied to identify trait-related modules and then select key module-related DEGs. The optimized gene set of signatures was obtained using the least absolute shrinkage and selection operator (LASSO) regression analysis. The CIBERSORT algorithm and Pearson correlation test were used to analyze the relationship between gene signatures and immune cells. Finally, public databases were used to predict the transcription factors (TFs) and upstream miRNAs.Results: A total of 127 DEGs in COPD were identified from the combined dataset. By considering the intersection of DEGs and genes in two trait-related modules, 83 key module-related DEGs were identified, which were mainly enriched in interleukin-related pathways. Seven-gene signatures, including MTHFD2, KANK3, GFPT2, PHLDA1, HS3ST2, FGG, and RPS4Y1, were further selected using the LASSO algorithm. These gene signatures showed the predictive potential for COPD risks and were significantly correlated with 18 types of immune cells. Finally, nine miRNAs and three TFs were predicted to target MTHFD2, GFPT2, PHLDA1, and FGG.Conclusion: We proposed the seven-gene-signature to predict COPD risk and explored its potential immune characteristics and regulatory mechanisms.Keywords: COPD, gene signatures, immune infiltration, transcriptional regulation

Background: identification of molecular mechanisms that determine tumour progression in cancer patients is a prerequisite for developing new disease treatment guidelines. Even though the predictive performance of current machine learning models is promising, extracting significant and meaningful knowledge from the data simultaneously during the learning process is a difficult task considering the high-dimensional and highly correlated nature of genomic datasets. Thus, there is a need for models that not only predict tumour volume from gene expression data of patients but also use prior information coming from pathway/gene sets during the learning process, to distinguish molecular mechanisms which play crucial role in tumour progression and therefore, disease prognosis. Results: in this study, instead of initially choosing several pathways/gene sets from an available set and training a model on this previously chosen subset of genomic features, we built a novel machine learning algorithm, PrognosiT, that accomplishes both tasks together. We tested our algorithm on thyroid carcinoma patients using gene expression profiles and cancer-specific pathways/gene sets. Predictive performance of our novel multiple kernel learning algorithm (PrognosiT) was comparable or even better than random forest (RF) and support vector regression (SVR). It is also notable that, to predict tumour volume, PrognosiT used gene expression features less than one-tenth of what RF and SVR algorithms used. Conclusions: PrognosiT was able to obtain comparable or even better predictive performance than SVR and RF. Moreover, we demonstrated that during the learning process, our algorithm managed to extract relevant and meaningful pathway/gene sets information related to the studied cancer type, which provides insights about its progression and aggressiveness. We also compared gene expressions of the selected genes by our algorithm in tumour and normal tissues, and we then discussed up- and down-regulated genes selected by our algorithm while learning, which could be beneficial for determining new biomarkers., Scientific and Technological Research Council of Turkey (TÜBİTAK); Turkish Academy of Science (TÜBA-GEBİP) The Young Scientist Award Program; Science Academy of Turkey (BAGEP) The Young Scientist Award Program

Background Amyloids are insoluble fibrillar aggregates that are highly associated with complex human diseases, such as Alzheimer’s disease, Parkinson’s disease, and type II diabetes. Recently, many studies reported that some specific regions of amino acid sequences may be responsible for the amyloidosis of proteins. It has become very important for elucidating the mechanism of amyloids that identifying the amyloidogenic regions. Accordingly, several computational methods have been put forward to discover amyloidogenic regions. The majority of these methods predicted amyloidogenic regions based on the physicochemical properties of amino acids. In fact, position, order, and correlation of amino acids may also influence the amyloidosis of proteins, which should be also considered in detecting amyloidogenic regions. Results To address this problem, we proposed a novel machine-learning approach for predicting amyloidogenic regions, called ReRF-Pred. Firstly, the pseudo amino acid composition (PseAAC) was exploited to characterize physicochemical properties and correlation of amino acids. Secondly, tripeptides composition (TPC) was employed to represent the order and position of amino acids. To improve the distinguishability of TPC, all possible tripeptides were analyzed by the binomial distribution method, and only those which have significantly different distribution between positive and negative samples remained. Finally, all samples were characterized by PseAAC and TPC of their amino acid sequence, and a random forest-based amyloidogenic regions predictor was trained on these samples. It was proved by validation experiments that the feature set consisted of PseAAC and TPC is the most distinguishable one for detecting amyloidosis. Meanwhile, random forest is superior to other concerned classifiers on almost all metrics. To validate the effectiveness of our model, ReRF-Pred is compared with a series of gold-standard methods on two datasets: Pep-251 and Reg33. The results suggested our method has the best overall performance and makes significant improvements in discovering amyloidogenic regions. Conclusions The advantages of our method are mainly attributed to that PseAAC and TPC can describe the differences between amyloids and other proteins successfully. The ReRF-Pred server can be accessed at http://106.12.83.135:8080/ReRF-Pred/.

The discovery of various noncoding RNAs (ncRNAs) and their biological implications is a growing area in cell biology. Increasing evidence has revealed canonical and noncanonical functions of long and small ncRNAs, including microRNAs, long ncRNAs (lncRNAs), circular RNAs, PIWI-interacting RNAs, and tRNA-derived fragments. These ncRNAs have the ability to regulate gene expression and modify metabolic pathways. Thus, they may have important roles as diagnostic biomarkers or therapeutic targets in various diseases, including neurodegenerative disorders, especially Parkinson’s disease. Recently, through diverse sequencing technologies and a wide variety of bioinformatic analytical tools, such as reverse transcriptase quantitative PCR, microarrays, next-generation sequencing and long-read sequencing, numerous ncRNAs have been shown to be associated with neurodegenerative disorders, including Parkinson’s disease. In this review article, we will first introduce the biogenesis of different ncRNAs, including microRNAs, PIWI-interacting RNAs, circular RNAs, long noncoding RNAs, and tRNA-derived fragments. The pros and cons of the detection platforms of ncRNAs and the reproducibility of bioinformatic analytical tools will be discussed in the second part. Finally, the recent discovery of numerous PD-associated ncRNAs and their association with the diagnosis and pathophysiology of PD are reviewed, and microRNAs and long ncRNAs that are transported by exosomes in biofluids are particularly emphasized. Supplementary Information The online version contains supplementary material available at 10.1186/s12929-021-00775-x.

Background Single-cell RNA sequencing (scRNA-seq) provides new insights to address biological and medical questions, and it will benefit more from the ultralow input RNA or subcellular sequencing. Results Here, we present a highly sensitive library construction protocol for ultralow input RNA sequencing (ulRNA-seq). We systematically evaluate experimental conditions of this protocol, such as reverse transcriptase, template-switching oligos (TSO), and template RNA structure. It was found that Maxima H Minus reverse transcriptase and rN modified TSO, as well as all RNA templates capped with m7G improved the sequencing sensitivity and low abundance gene detection ability. RNA-seq libraries were successfully prepared from total RNA samples as low as 0.5 pg, and more than 2000 genes have been identified. Conclusions The ability of low abundance gene detection and sensitivity were largely enhanced with this optimized protocol. It was also confirmed in single-cell sequencing, that more genes and cell markers were identified compared to conventional sequencing method. We expect that ulRNA-seq will sequence and transcriptome characterization for the subcellular of disease tissue, to find the corresponding treatment plan.

Glioblastoma is the most aggressive type of brain tumors resistant to a number of antitumor drugs. The problem of therapy and drug treatment course is complicated by extremely high heterogeneity in the benign cell populations, the random arrangement of tumor cells, and polymorphism of their nuclei. The pathogenesis of gliomas needs to be studied using modern cellular technologies, genome- and transcriptome-wide technologies of high-throughput sequencing, analysis of gene expression on microarrays, and methods of modern bioinformatics to find new therapy targets. Functional annotation of genes related to the disease could be retrieved based on genetic databases and cross-validated by integrating complementary experimental data. Gene network reconstruction for a set of genes (proteins) proved to be effective approach to study mechanisms underlying disease progression. We used online bioinformatics tools for annotation of gene list for glioma, reconstruction of gene network and comparative analysis of gene ontology categories. The available tools and the databases for glioblastoma gene analysis are discussed together with the recent progress in this field.

Background CRISPR-Cas9 genome-wide screens are being increasingly performed, allowing systematic explorations of cancer dependencies at unprecedented accuracy and scale. One of the major computational challenges when analysing data derived from such screens is to identify genes that are essential for cell survival invariantly across tissues, conditions, and genomic-contexts (core-fitness genes), and to distinguish them from context-specific essential genes. This is of paramount importance to assess the safety profile of candidate therapeutic targets and for elucidating mechanisms involved in tissue-specific genetic diseases. Results We have developed CoRe: an R package implementing existing and novel methods for the identification of core-fitness genes (at two different level of stringency) from joint analyses of multiple CRISPR-Cas9 screens. We demonstrate, through a fully reproducible benchmarking pipeline, that CoRe outperforms state-of-the-art tools, yielding more reliable and biologically relevant sets of core-fitness genes. Conclusions CoRe offers a flexible pipeline, compatible with many pre-processing methods for the analysis of CRISPR data, which can be tailored onto different use-cases. The CoRe package can be used for the identification of high-confidence novel core-fitness genes, as well as a means to filter out potentially cytotoxic hits while analysing cancer dependency datasets for identifying and prioritising novel selective therapeutic targets.

Background Birth defects pose a major challenge to infant health. Thus far, however, the causes of most birth defects remain cryptic. Over the past few decades, considerable effort has been expended on disclosing the underlying mechanisms related to birth defects, yielding myriad treatises and data. To meet the increasing requirements for data resources, we developed a freely accessible birth defect multi-omics database (BDdb, http://t21omics.cngb.org) consisting of multi-omics data and potential disease biomarkers. Results In total, omics datasets from 136 Gene Expression Omnibus (GEO) Series records, including 5245 samples, as well as 869 biomarkers of 22 birth defects in six different species, were integrated into the BDdb. The database provides a user-friendly interface for searching, browsing, and downloading data of interest. The BDdb also enables users to explore the correlations among different sequencing methods, such as chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq) from different studies, to obtain the information on gene expression patterns from diverse aspects. Conclusion To the best of our knowledge, the BDdb is the first comprehensive database associated with birth defects, which should benefit the diagnosis and prevention of birth defects.

Refereed/Peer-reviewed Cancer treatment is constantly evolving from a one-size-fits-all towards bespoke approaches for each patient. In certain solid cancers, including breast and lung, tumor genome profiling has been incorporated into therapeutic decision-making. For chronic phase chronic myeloid leukemia (CML), while tyrosine kinase inhibitor therapy is the standard treatment, current clinical scoring systems cannot accurately predict the heterogeneous treatment outcomes observed in patients. Biomarkers capable of segregating patients according to outcome at diagnosis are needed to improve management, and facilitate enrollment in clinical trials seeking to prevent blast crisis transformation and improve the depth of molecular responses. To this end, gene expression (GE) profiling studies have evaluated whether GE signatures at diagnosis are clinically informative. Patient material from a variety of sources has been profiled using microarrays, RNA sequencing and, more recently, single-cell RNA sequencing. However, differences in the cell types profiled, the technologies used, and the inherent complexities associated with the interpretation of genomic data pose challenges in distilling GE datasets into biomarkers with clinical utility. The goal of this paper is to review previous studies evaluating GE profiling in CML, and explore their potential as risk assessment tools for individualized CML treatment. We also review the contribution that acquired mutations, including those seen in clonal hematopoiesis, make to GE profiles, and how a model integrating contributions of genetic and epigenetic factors in resistance to tyrosine kinase inhibitors and blast crisis transformation can define a route to GE-based biomarkers. Finally, we outline a four-stage approach for the development of GE-based biomarkers in CML.

Background Several computational tools for predicting protein Ubiquitylation and SUMOylation sites have been proposed to study their regulatory roles in gene location, gene expression, and genome replication. However, existing methods generally rely on feature engineering, and ignore the natural similarity between the two types of protein translational modification. This study is the first all-in-one deep network to predict protein Ubiquitylation and SUMOylation sites from protein sequences as well as their crosstalk sites simultaneously. Our deep learning architecture integrates several meta classifiers that apply deep neural networks to protein sequence information and physico-chemical properties, which were trained on multi-label classification mode for simultaneously identifying protein Ubiquitylation and SUMOylation as well as their crosstalk sites. Results The promising AUCs of our method on Ubiquitylation, SUMOylation and crosstalk sites achieved 0.838, 0.888, and 0.862 respectively on tenfold cross-validation. The corresponding APs reached 0.683, 0.804 and 0.552, which also validated our effectiveness. Conclusions The proposed architecture managed to classify ubiquitylated and SUMOylated lysine residues along with their crosstalk sites, and outperformed other well-known Ubiquitylation and SUMOylation site prediction tools.

Background Phylogenetic profiling is widely used to predict novel members of large protein complexes and biological pathways. Although methods combined with phylogenetic trees have significantly improved prediction accuracy, computational efficiency is still an issue that limits its genome-wise application. Results Here we introduce a new tree-based phylogenetic profiling algorithm named GFICLEE, which infers common single and continuous loss (SCL) events in the evolutionary patterns. We validated our algorithm with human pathways from three databases and compared the computational efficiency with current tree-based with 10 different scales genome dataset. Our algorithm has a better predictive performance with high computational efficiency. Conclusions The GFICLEE is a new method to infers genome-wide gene function. The accuracy and computational efficiency of GFICLEE make it possible to explore gene functions at the genome-wide level on a personal computer.

Background Accurate copy number variant (CNV) detection is especially challenging for both targeted sequencing (TS) and whole‐exome sequencing (WES) data. To maximize the performance, the parameters of the CNV calling algorithms should be optimized for each specific dataset. This requires obtaining validated CNV information using either multiplex ligation-dependent probe amplification (MLPA) or array comparative genomic hybridization (aCGH). They are gold standard but time-consuming and costly approaches. Results We present isoCNV which optimizes the parameters of DECoN algorithm using only NGS data. The parameter optimization process is performed using an in silico CNV validated dataset obtained from the overlapping calls of three algorithms: CNVkit, panelcn.MOPS and DECoN. We evaluated the performance of our tool and showed that increases the sensitivity in both TS and WES real datasets. Conclusions isoCNV provides an easy-to-use pipeline to optimize DECoN that allows the detection of analysis-ready CNV from a set of DNA alignments obtained under the same conditions. It increases the sensitivity of DECoN without the need for orthogonal methods. isoCNV is available at https://gitlab.com/sequentiateampublic/isocnv.

Background Paralogs formed through gene duplication and isoforms formed through alternative splicing have been important processes for increasing protein diversity and maintaining cellular homeostasis. Despite their recognized importance and the advent of large-scale genomic and transcriptomic analyses, paradoxically, accurate annotations of all gene loci to allow the identification of paralogs and isoforms remain surprisingly incomplete. In particular, the global analysis of the transcriptome of a non-model organism for which there is no reference genome is especially challenging. Results To reliably discriminate between the paralogs and isoforms in RNA-seq data, we redefined the pre-existing sequence features (sequence similarity, inverse count of consecutive identical or non-identical blocks, and match-mismatch fraction) previously derived from full-length cDNAs and EST sequences and described newly discovered genomic and transcriptomic features (twilight zone of protein sequence alignment and expression level difference). In addition, the effectiveness and relevance of the proposed features were verified with two widely used support vector machine (SVM) and random forest (RF) models. From nine RNA-seq datasets, all AUC (area under the curve) scores of ROC (receiver operating characteristic) curves were over 0.9 in the RF model and significantly higher than those in the SVM model. Conclusions In this study, using an RF model with five proposed RNA-seq features, we implemented our method called Paralogs and Isoforms Classifier based on Machine-learning approaches (PIC-Me) and showed that it outperformed an existing method. Finally, we envision that our tool will be a valuable computational resource for the genomics community to help with gene annotation and will aid in comparative transcriptomics and evolutionary genomics studies, especially those on non-model organisms.

Background Atherosclerosis is the leading cause of cardiovascular disease with a high mortality worldwide. Understanding the atherosclerosis pathogenesis and identification of efficient diagnostic signatures remain major problems of modern medicine. This study aims to screen the potential diagnostic genes for atherosclerosis. Methods We downloaded the gene chip data of 135 peripheral blood samples, including 57 samples with atherosclerosis and 78 healthy subjects from GEO database (Accession Number: GSE20129). The weighted gene co-expression network analysis was applied to identify atherosclerosis-related genes. Functional enrichment analysis was conducted by using the clusterProfiler R package. The interaction pairs of proteins encoded by atherosclerosis-related genes were screened using STRING database, and the interaction network was further optimized with the cytoHubba plug-in of Cytoscape software. Results The logistic regression diagnostic model was constructed to predict normal and atherosclerosis samples. A gene module which included 532 genes related to the occurrence of atherosclerosis were screened. Functional enrichment analysis basing on the 532 genes identified 235 significantly enriched GO terms and 44 significantly enriched KEGG pathways. The top 50 hub genes of the protein–protein interaction network were identified. The final logistic regression diagnostic model was established by the optimal 10 key genes, which could distinguish atherosclerosis samples from normal samples. Conclusions A predictive model based on 10 potential atherosclerosis-related genes was obtained, which should shed light on the diagnostic research of atherosclerosis.

Background Generating high-quality de novo genome assemblies is foundational to the genomics study of model and non-model organisms. In recent years, long-read sequencing has greatly benefited genome assembly and scaffolding, a process by which assembled sequences are ordered and oriented through the use of long-range information. Long reads are better able to span repetitive genomic regions compared to short reads, and thus have tremendous utility for resolving problematic regions and helping generate more complete draft assemblies. Here, we present LongStitch, a scalable pipeline that corrects and scaffolds draft genome assemblies exclusively using long reads. Results LongStitch incorporates multiple tools developed by our group and runs in up to three stages, which includes initial assembly correction (Tigmint-long), followed by two incremental scaffolding stages (ntLink and ARKS-long). Tigmint-long and ARKS-long are misassembly correction and scaffolding utilities, respectively, previously developed for linked reads, that we adapted for long reads. Here, we describe the LongStitch pipeline and introduce our new long-read scaffolder, ntLink, which utilizes lightweight minimizer mappings to join contigs. LongStitch was tested on short and long-read assemblies of Caenorhabditis elegans, Oryza sativa, and three different human individuals using corresponding nanopore long-read data, and improves the contiguity of each assembly from 1.2-fold up to 304.6-fold (as measured by NGA50 length). Furthermore, LongStitch generates more contiguous and correct assemblies compared to state-of-the-art long-read scaffolder LRScaf in most tests, and consistently improves upon human assemblies in under five hours using less than 23 GB of RAM. Conclusions Due to its effectiveness and efficiency in improving draft assemblies using long reads, we expect LongStitch to benefit a wide variety of de novo genome assembly projects. The LongStitch pipeline is freely available at https://github.com/bcgsc/longstitch.

Background Liver cancer is a common malignant tumor in China, with high mortality. Its occurrence and development were thoroughly studied by high-throughput expression microarray, which produced abundant data on gene expression, mRNA quantification and the clinical data of liver cancer. However, the hub genes, which can be served as biomarkers for diagnosis and treatment of early liver cancer, are not well screened. Results Here we present a new method for getting 6 key genes, aiming to diagnose and treat the early liver cancer. We firstly analyzed the different expression microarrays based on TCGA database, and a total of 1564 differentially expressed genes were obtained, of which 1400 were up-regulated and 164 were down-regulated. Furthermore, these differentially expressed genes were studied by using GO and KEGG enrichment analysis, a PPI network was constructed based on the STRING database, and 15 hub genes were obtained. Finally, 15 hub genes were verified by applying the survival analysis method on Oncomine database, and 6 key genes were ultimately identified, including PLK1, CDC20, CCNB2, BUB1, MAD2L1 and CCNA2. The robustness analysis of four independent data sets verifies the accuracy of the key gene’s classification of the data set. Conclusions Although there are complicated differences between cancer and normal cells in gene functions, cancer cells could be differentiated in case that a group of special genes expresses abnormally. Here we presented a new method to identify the 6 key genes for diagnosis and treatment of early liver cancer, and these key genes can help us understand the pathogenesis of liver cancer more deeply.

Background Comparisons of the molecular framework among organisms can be done on both structural and functional levels. One of the most common top-down approaches for functional comparisons is RNA sequencing. This estimation of organismal transcriptional responses is of interest for understanding evolution of molecular activity, which is used for answering a diversity of questions ranging from basic biology to pre-clinical species selection and translation. However, direct comparison between species is often hindered by evolutionary divergence in structure of molecular framework, as well as large difference in the depth of our understanding of the genetic background between humans and other species. Here, we focus on the latter. We attempt to understand how differences in transcriptome annotation affect direct gene abundance comparisons between species. Results We examine and suggest some straightforward approaches for direct comparison given the current available tools and using a sample dataset from human, cynomolgus monkey, dog, rat and mouse with a common quantitation and normalization approach. In addition, we examine how variation in genome annotation depth and quality across species may affect these direct comparisons. Conclusions Our findings suggest that further efforts for better genome annotation or computational normalization tools may be of strong interest.

Background Reference sequences play a vital role in next-generation sequencing (NGS), impacting mapping quality during genome analyses. However, reference genomes usually do not represent the full range of genetic diversity of a species as a result of geographical divergence and independent demographic events of different populations. For the mitochondrial genome (mitogenome), which occurs in high copy numbers in cells and is strictly maternally inherited, an optimal reference sequence has the potential to make mitogenome alignment both more accurate and more efficient. In this study, we used three different types of reference sequences for mitogenome mapping, i.e., the commonly used reference sequence (CU-ref), the breed-specific reference sequence (BS-ref) and the sample-specific reference sequence (SS-ref), respectively, and compared the accuracy of mitogenome alignment and SNP calling among them, for the purpose of proposing the optimal reference sequence for mitochondrial DNA (mtDNA) analyses of specific populations Results Four pigs, representing three different breeds, were high-throughput sequenced, subsequently mapping reads to the reference sequences mentioned above, resulting in a largest mapping ratio and a deepest coverage without increased running time when aligning reads to a BS-ref. Next, single nucleotide polymorphism (SNP) calling was carried out by 18 detection strategies with the three tools SAMtools, VarScan and GATK with different parameters, using the bam results mapping to BS-ref. The results showed that all eighteen strategies achieved the same high specificity and sensitivity, which suggested a high accuracy of mitogenome alignment by the BS-ref because of a low requirement for SNP calling tools and parameter choices. Conclusions This study showed that different reference sequences representing different genetic relationships to sample reads influenced mitogenome alignment, with the breed-specific reference sequences being optimal for mitogenome analyses, which provides a refined processing perspective for NGS data.

Background Various epigenetic factors are responsible for the non-genetic regulation on gene expression. The epigenetically dysregulated oncogenes or tumor suppressors by miRNA and/or DNA methylation are often observed in cancer cells. Each of these epigenetic regulators has been studied well in cancer progressions; however, their mutual regulatory relationship in cancer still remains unclear. In this study, we propose an integrative framework to systematically investigate epigenetic interactions between miRNA and methylation at the alternatively spliced mRNA level in bladder cancer. Each of these epigenetic regulators has been studied well in cancer progressions; however, their mutual regulatory relationship in cancer still remains unclear. Results The integrative analyses yielded 136 significant combinations (methylation, miRNA and isoform). Further, overall survival analysis on the 136 combinations based on methylation and miRNA, high and low expression groups resulted in 13 combinations associated with survival. Additionally, different interaction patterns were examined. Conclusions Our study provides a higher resolution of molecular insight into the crosstalk between two epigenetic factors, DNA methylation and miRNA. Given the importance of epigenetic interactions and alternative splicing in cancer, it is timely to identify and understand the underlying mechanisms based on epigenetic markers and their interactions in cancer, leading to alternative splicing with primary functional impact.

Background Differential expression (DE) analysis of RNA-seq data typically depends on gene annotations. Different sets of gene annotations are available for the human genome and are continually updated–a process complicated with the development and application of high-throughput sequencing technologies. However, the impact of the complexity of gene annotations on DE analysis remains unclear. Results Using “mappability”, a metric of the complexity of gene annotation, we compared three distinct human gene annotations, GENCODE, RefSeq, and NONCODE, and evaluated how mappability affected DE analysis. We found that mappability was significantly different among the human gene annotations. We also found that increasing mappability improved the performance of DE analysis, and the impact of mappability mainly evident in the quantification step and propagated downstream of DE analysis systematically. Conclusions We assessed how the complexity of gene annotations affects DE analysis using mappability. Our findings indicate that the growth and complexity of gene annotations negatively impact the performance of DE analysis, suggesting that an approach that excludes unnecessary gene models from gene annotations improves the performance of DE analysis.

Background Lately, high-throughput RNA sequencing has been extensively used to elucidate the transcriptome landscape and dynamics of cell types of different species. In particular, for most non-model organisms lacking complete reference genomes with high-quality annotation of genetic information, reference-free (RF) de novo transcriptome analyses, rather than reference-based (RB) approaches, are widely used, and RF analyses have substantially contributed toward understanding the mechanisms regulating key biological processes and functions. To date, numerous bioinformatics studies have been conducted for assessing the workflow, production rate, and completeness of transcriptome assemblies within and between RF and RB datasets. However, the degree of consistency and variability of results obtained by analyzing gene expression levels through these two different approaches have not been adequately documented. Results In the present study, we evaluated the differences in expression profiles obtained with RF and RB approaches and revealed that the former tends to be satisfactorily replaced by the latter with respect to transcriptome repertoires, as well as from a gene expression quantification perspective. In addition, we urge cautious interpretation of these findings. Several genes that are lowly expressed, have long coding sequences, or belong to large gene families must be validated carefully, whenever gene expression levels are calculated using the RF method. Conclusions Our empirical results indicate important contributions toward addressing transcriptome-related biological questions in non-model organisms.

Background Circular RNA (circRNA) is an emerging class of RNA molecules attracting researchers due to its potential for serving as markers for diagnosis, prognosis, or therapeutic targets of cancer, cardiovascular, and autoimmune diseases. Current methods for detection of circRNA from RNA sequencing (RNA-seq) focus mostly on improving mapping quality of reads supporting the back-splicing junction (BSJ) of a circRNA to eliminate false positives (FPs). We show that mapping information alone often cannot predict if a BSJ-supporting read is derived from a true circRNA or not, thus increasing the rate of FP circRNAs. Results We have developed Circall, a novel circRNA detection method from RNA-seq. Circall controls the FPs using a robust multidimensional local false discovery rate method based on the length and expression of circRNAs. It is computationally highly efficient by using a quasi-mapping algorithm for fast and accurate RNA read alignments. We applied Circall on two simulated datasets and three experimental datasets of human cell-lines. The results show that Circall achieves high sensitivity and precision in the simulated data. In the experimental datasets it performs well against current leading methods. Circall is also substantially faster than the other methods, particularly for large datasets. Conclusions With those better performances in the detection of circRNAs and in computational time, Circall facilitates the analyses of circRNAs in large numbers of samples. Circall is implemented in C++ and R, and available for use at https://www.meb.ki.se/sites/biostatwiki/circall and https://github.com/datngu/Circall.

Circular RNAs (circRNAs) are novel biomarkers of various cancers. CircRNAs can sponge miRNAs and regulate target mRNAs, which was called competing endogenous RNAs (ceRNA). This study was designed to identify circRNAs related to patients with clear cell renal cell carcinoma (ccRCC) and the first to select three independent Gene Expression Omnibus microarrays covering circRNAs, miRNAs, and mRNAs for multiple analyses. The data of clinical cases applied in our study were obtained from The Cancer Genome Atlas. We successfully conducted a circRNA/miRNA/mRNA ceRNA network related to ccRCC patients via R software and Cytoscape including 8 circRNAs, 6 miRNAs, and 49 mRNAs. The prognosis‐associated subnet covered 8 circRNAs, 6 miRNAs, and 22 mRNAs. Quantitative real‐time PCR was applied to measure our prediction in three renal cell lines and 23 pairs of tissues. Small interfering RNA targeting the back‐splice region of hsa_circ_0001167 was further implied to confirm the regulation. Ultimately, hsa_circ_0001167/hsa‐miR‐595/CCDC8 regulatory axis was identified in this study, which may serve as prognostic indicators. Lower levels of hsa_circ_0001167 and CCDC8 were potentially correlated with worse patient survival., Based on the analyzation of three Gene Expression Omnibus microarrays and The Cancer Genome Atlas cases, we were the first to construct a ccRCC centred ceRNA network and its prognosis‐associated subnet. We further validated a hsa_circ_0001167/ has‐miR‐595/CCDC8 axis in renal cell lines, which may serve as prognostic indictors.

Recent advances in next-generation sequencing, deep networks and other bioinformatic tools have enabled us to mine huge amount of genomic information about living organisms in the post-microarray era. However, these tools do not explicitly factor in the role of the underlying DNA biochemistry (particularly, DNA hybridization) essential to life processes. Here, we focus more precisely on the role that DNA hybridization plays in determining properties of biological organisms at the macro-level. We illustrate its role with solutions to challenging problems in human disease. These solutions are made possible by novel structural properties of DNA hybridization landscapes revealed by a metric model of oligonucleotides of a common length that makes them reminiscent of some planets in our solar system, particularly Earth and Saturn. They allow a judicious selection of so-called noncrosshybridizing (nxh) bases that offer substantial reduction of DNA sequences of arbitrary length into a few informative features. The quality assessment of the information extracted by them is high because of their very low Shannon Entropy, i.e. they minimize the degree of uncertainty in hybridization that makes results on standard microarrays irreproducible. For example, SNP classification (pathogenic/non-pathogenic) and pathogen identification can be solved with high sensitivity (~77%/100%) and specificity (~92%/100%, respectively) for combined taxa on a sample of over 264 fully coding sequences in whole bacterial genomes and fungal mitochondrial genomes using machine learning (ML) models. These methods can be applied to several other interesting research questions that could be addressed with similar genomic analyses.

Background With the broad application of high-throughput sequencing and its reduced cost, simple sequence repeat (SSR) genotyping by sequencing (SSR-GBS) has been widely used for interpreting genetic data across different fields, including population genetic diversity and structure analysis, the construction of genetic maps, and the investigation of intraspecies relationships. The development of accurate and efficient typing strategies for SSR-GBS is urgently needed and several tools have been published. However, to date, no suitable accurate genotyping method can tolerate single nucleotide variations (SNVs) in SSRs and flanking regions. These SNVs may be caused by PCR and sequencing errors or SNPs among varieties, and they directly affect sequence alignment and genotyping accuracy. Results Here, we report a new integrated strategy named the accurate microsatellite genotyping tool based on targeted sequencing (AMGT-TS) and provide a user-friendly web-based platform and command-line version of AMGT-TS. To handle SNVs in the SSRs or flanking regions, we developed a broad matching algorithm (BMA) that can quickly and accurately achieve SSR typing for ultradeep coverage and high-throughput analysis of loci with SNVs compatibility and grouping of typed reads for further in-depth information mining. To evaluate this tool, we tested 21 randomly sampled loci in eight maize varieties, accompanied by experimental validation on actual and simulated sequencing data. Our evaluation showed that, compared to other tools, AMGT-TS presented extremely accurate typing results with single base resolution for both homozygous and heterozygous samples. Conclusion This integrated strategy can achieve accurate SSR genotyping based on targeted sequencing, and it can tolerate single nucleotide variations in the SSRs and flanking regions. This method can be readily applied to divergent sequencing platforms and species and has excellent application prospects in genetic and population biology research. The web-based platform and command-line version of AMGT-TS are available at https://amgt-ts.plantdna.site:8445 and https://github.com/plantdna/amgt-ts, respectively.

Background Zebrafish is a popular animal model used for high-throughput screening of chemical hazards, however, investigations of transcriptomic mechanisms of toxicity are still needed. Here, our goal was to identify genes and biological pathways that Aryl Hydrocarbon Receptor 2 (AHR2) Activators and flame retardant chemicals (FRCs) alter in developing zebrafish. Taking advantage of a compendium of phenotypically-anchored RNA sequencing data collected from 48-h post fertilization (hpf) zebrafish, we inferred a co-expression network that grouped genes based on their transcriptional response. Results Genes responding to the FRCs and AHR2 Activators localized to distinct regions of the network, with FRCs inducing a broader response related to neurobehavior. AHR2 Activators centered in one region related to chemical stress responses. We also discovered several highly co-expressed genes in this module, including cyp1a, and we subsequently show that these genes are definitively within the AHR2 signaling pathway. Systematic removal of the two chemical types from the data, and analysis of network changes identified neurogenesis associated with FRCs, and regulation of vascular development associated with both chemical classes. We also identified highly connected genes responding specifically to each class that are potential biomarkers of exposure. Conclusions Overall, we created the first zebrafish chemical-specific gene co-expression network illuminating how chemicals alter the transcriptome relative to each other. In addition to our conclusions regarding FRCs and AHR2 Activators, our network can be leveraged by other studies investigating chemical mechanisms of toxicity.

Background Variation in mitochondrial DNA (mtDNA) identified by genotyping microarrays or by sequencing only the hypervariable regions of the genome may be insufficient to reliably assign mitochondrial genomes to phylogenetic lineages or haplogroups. This lack of resolution can limit functional and clinical interpretation of a substantial body of existing mtDNA data. To address this limitation, we developed and evaluated a large, curated reference alignment of complete mtDNA sequences as part of a pipeline for imputing missing mtDNA single nucleotide variants (mtSNVs). We call our reference alignment and pipeline MitoImpute. Results We aligned the sequences of 36,960 complete human mitochondrial genomes downloaded from GenBank, filtered and controlled for quality. These sequences were reformatted for use in imputation software, IMPUTE2. We assessed the imputation accuracy of MitoImpute by measuring haplogroup and genotype concordance in data from the 1000 Genomes Project and the Alzheimer’s Disease Neuroimaging Initiative (ADNI). The mean improvement of haplogroup assignment in the 1000 Genomes samples was 42.7% (Matthew’s correlation coefficient = 0.64). In the ADNI cohort, we imputed missing single nucleotide variants. Conclusion These results show that our reference alignment and panel can be used to impute missing mtSNVs in existing data obtained from using microarrays, thereby broadening the scope of functional and clinical investigation of mtDNA. This improvement may be particularly useful in studies where participants have been recruited over time and mtDNA data obtained using different methods, enabling better integration of early data collected using less accurate methods with more recent sequence data.

MotivationComputationally predicting the specificity of T cell receptors can be a powerful tool to shed light on the immune response against infectious diseases and cancers, autoimmunity, cancer immunotherapy, and immunopathology. With more T cell receptor sequence data becoming available, the need for bioinformatics approaches to tackle this problem is even more pressing. Here we present SwarmTCR, a method that uses labeled sequence data to predict the specificity of T cell receptors using a nearest-neighbor approach. SwarmTCR works by optimizing the weights of the individual CDR regions to maximize classification performance.ResultsWe compared the performance of SwarmTCR against a state-of-the-art method (TCRdist) and showed that SwarmTCR performed significantly better on epitopes EBV-BRLF1300, EBV-BRLF1109, NS4B214–222 with single cell data and epitopes EBV-BRLF1300, EBV-BRLF1109, IAV-M158 with bulk sequencing data (α and β chains). In addition, we show that the weights returned by SwarmTCR are biologically interpretable.AvailabilitySwarmTCR is distributed freely under the terms of the GPL-3 license. The source code and all sequencing data are available at GitHub (https://github.com/thecodingdoc/SwarmTCR)Contactdghersi@unomaha.edu

BackgroundOne of the major challenges in precision medicine is accurate prediction of individual patient’s response to drugs. A great number of computational methods have been developed to predict compounds activity using genomic profiles or chemical structures, but more exploration is yet to be done to combine genetic mutation, gene expression, and cheminformatics in one machine learning model.ResultsWe presented here a novel deep-learning model that integrates gene expression, genetic mutation, and chemical structure of compounds in a multi-task convolutional architecture. We applied our model to the Genomics of Drug Sensitivity in Cancer (GDSC) and Cancer Cell Line Encyclopedia (CCLE) datasets. We selected relevant cancer-related genes based on oncology genetics database and L1000 landmark genes, and used their expression and mutations as genomic features in model training. We obtain the cheminformatics features for compounds from PubChem or ChEMBL. Our finding is that combining gene expression, genetic mutation, and cheminformatics features greatly enhances the predictive performance.ConclusionWe implemented an extended Graph Neural Network for molecular graphs and Convolutional Neural Network for gene features. With the employment of multi-tasking and self-attention functions to monitor the similarity between compounds, our model outperforms recently published methods using the same training and testing datasets.

Background Huntington’s disease is a kind of chronic progressive neurodegenerative disease with complex pathogenic mechanisms. To data, the pathogenesis of Huntington’s disease is still not fully understood, and there has been no effective treatment. The rapid development of high-throughput sequencing technologies makes it possible to explore the molecular mechanisms at the transcriptome level. Our previous studies on Huntington’s disease have shown that it is difficult to distinguish disease-associated genes from non-disease genes. Meanwhile, recent progress in bio-medicine shows that the molecular origin of chronic complex diseases may not exist in the diseased tissue, and differentially expressed genes between different tissues may be helpful to reveal the molecular origin of chronic diseases. Therefore, developing integrative analysis computational methods for the multi-tissues gene expression data, exploring the relationship between differentially expressed genes in different tissues and the disease, can greatly accelerate the molecular discovery process. Methods For analysis of the intra- and inter- tissues’ differentially expressed genes, we designed an integrative enrichment analysis method based on an artificial neuron (IEAAN). Firstly, we calculated the differential expression scores of genes which are seen as features of the corresponding gene, using fold-change approach with intra- and inter- tissues’ gene expression data. Then, we weighted sum all the differential expression scores through a sigmoid function to get differential expression enrichment score. Finally, we ranked the genes according to the enrichment score. Top ranking genes are supposed to be the potential disease-associated genes. Results In this study, we conducted large amounts of experiments to analyze the differentially expressed genes of intra- and inter- tissues. Experimental results showed that genes differentially expressed between different tissues are more likely to be Huntington’s disease-associated genes. Five disease-associated genes were selected out in this study, two of which have been reported to be implicated in Huntington’s disease. Conclusions We proposed a novel integrative enrichment analysis method based on artificial neuron (IEAAN), which displays better prediction precision of disease-associated genes in comparison with the state-of-the-art statistical-based methods. Our comprehensive evaluation suggests that genes differentially expressed between striatum and liver tissues of health individuals are more likely to be Huntington’s disease-associated genes.

BackgroundClinically, behavior, cognitive, and mental functions are affected during the neurodegenerative disease progression. To date, the molecular pathogenesis of these complex disease is still unclear. With the rapid development of sequencing technologies, it is possible to delicately decode the molecular mechanisms corresponding to different clinical phenotypes at the genome-wide transcriptomic level using computational methods. Our previous studies have shown that it is difficult to distinguish disease genes from non-disease genes. Therefore, to precisely explore the molecular pathogenesis under complex clinical phenotypes, it is better to identify biomarkers corresponding to different disease stages or clinical phenotypes. So, in this study, we designed a label propagation-based semi-supervised feature selection approach (LPFS) to prioritize disease-associated genes corresponding to different disease stages or clinical phenotypes.MethodsIn this study, we pioneering put label propagation clustering and feature selection into one framework and proposed label propagation-based semi-supervised feature selection approach. LPFS prioritizes disease genes related to different disease stages or phenotypes through the alternative iteration of label propagation clustering based on sample network and feature selection with gene expression profiles. Then the GO and KEGG pathway enrichment analysis were carried as well as the gene functional analysis to explore molecular mechanisms of specific disease phenotypes, thus to decode the changes in individual behavioral and mental characteristics during neurodegenerative disease progression.ResultsLarge amounts of experiments were conducted to verify the performance of LPFS with Huntington’s gene expression data. Experimental results shown that LPFS performs better in comparison with the-state-of-art methods. GO and KEGG enrichment analysis of key gene sets shown that TGF-beta signaling pathway, cytokine-cytokine receptor interaction, immune response, and inflammatory response were gradually affected during the Huntington’s disease progression. In addition, we found that the expression of SLC4A11, ZFP474, AMBP, TOP2A, PBK, CCDC33, APSL, DLGAP5, and Al662270 changed seriously by the development of the disease.ConclusionsIn this study, we designed a label propagation-based semi-supervised feature selection model to precisely selected key genes of different disease phenotypes. We conducted experiments using the model with Huntington’s disease mice gene expression data to decode the mechanisms of it. We found many cell types, including astrocyte, microglia, and GABAergic neuron, could be involved in the pathological process.

Background Sorafenib, the first approved targeted therapy for advanced hepatocellular carcinoma (HCC), is often reported to comprised survival-benefit due to resistance. An underlying mechanism of resistance was proposed using bioinformatics analysis based on differentially expressed genes (DEGs) from microarrays. However, most DEGs were invalidated at both the expression level, and the role in causing resistance. Therefore, we conducted a bioinformatics analysis based on experimentally determined sorafenib-resistance-related genes (SRRGs) to elucidate the mechanism of sorafenib resistance. Methods The SRRGs, which have been experimentally determined to promote or inhibit resistance, were collected from published studies. The Database for Annotation, Visualization and Integrated Discovery (DAVID) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to perform Gene Ontology (GO) and pathway enrichment analysis, respectively. A corresponding protein-protein interaction network (PPI) was created using the Cytoscape software program, and network hub genes were proposed. Results A total of 145 SRRGs, with 117 promoting and 28 inhibiting resistance, were identified. Cell proliferation, migration, development, response to oxygen levels, epithelial-to-mesenchymal transition (EMT), cell skeleton, protein function, and autophagy were all proposed as crucial gene functions related to resistance. The pathways related to cell proliferation or apoptosis, immune function, endocrine metabolism, stem cell function, and differentiation were identified as key resistance-related pathways. A total of 81 hub genes were proposed, including the following top 10 genes: TP53, AKT1, EGFR, STAT3, VEGFA, JUN, MAPK1, IL6, PTEN, and CTNNB1. Conclusions In conclusion, this study gathered experimentally validated genes that determine sorafenib resistance in HCC, provided an overview of the underlying mechanisms of resistance, and further validated sorafenib resistance in HCC.

Background Genetic heterogeneity of a cancer tumor that develops during clonal evolution is one of the reasons for cancer treatment failure, by increasing the chance of drug resistance. Clones are cell populations with different genotypes, resulting from differences in somatic mutations that occur and accumulate during cancer development. An appropriate approach for identifying clones is determining the variant allele frequency of mutations that occurred in the tumor. Although bulk sequencing data can be used to provide that information, the frequencies are not informative enough for identifying different clones with the same prevalence and their evolutionary relationships. On the other hand, single-cell sequencing data provides valuable information about branching events in the evolution of a cancerous tumor. However, the temporal order of mutations may be determined with ambiguities using only single-cell data, while variant allele frequencies from bulk sequencing data can provide beneficial information for inferring the temporal order of mutations with fewer ambiguities. Result In this study, a new method called Conifer (ClONal tree Inference For hEterogeneity of tumoR) is proposed which combines aggregated variant allele frequency from bulk sequencing data with branching event information from single-cell sequencing data to more accurately identify clones and their evolutionary relationships. It is proven that the accuracy of clone identification and clonal tree inference is increased by using Conifer compared to other existing methods on various sets of simulated data. In addition, it is discussed that the evolutionary tree provided by Conifer on real cancer data sets is highly consistent with information in both bulk and single-cell data. Conclusions In this study, we have provided an accurate and robust method to identify clones of tumor heterogeneity and their evolutionary history by combining single-cell and bulk sequencing data.

Background Ovarian cancer is one of the most common gynecological tumors, and among gynecological tumors, its incidence and mortality rates are fairly high. However, the pathogenesis of ovarian cancer is not clear. The present study aimed to investigate the differentially expressed genes and signaling pathways associated with ovarian cancer by bioinformatics analysis. Methods The data from three mRNA expression profiling microarrays (GSE14407, GSE29450, and GSE54388) were obtained from the Gene Expression Omnibus (GEO) database. Differentially expressed genes between ovarian cancer tissues and normal tissues were identified using R software. The overlapping genes from the three GEO datasets were identified, and profound analysis was performed. The overlapping genes were used for pathway and Gene Ontology (GO) functional enrichment analysis using the Metascape online tool. Protein–protein interactions were analyzed with the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING). Subnetwork models were selected using the plugin molecular complex detection (MCODE) application in Cytoscape. Kaplan–Meier curves were used to analyze the univariate survival outcomes of the hub genes. The Human Protein Atlas (HPA) database and Gene Expression Profiling Interactive Analysis (GEPIA) were used to validate hub genes. Results In total, 708 overlapping genes were identified through analyses of the three microarray datasets (GSE14407, GSE29450, and GSE54388). These genes mainly participated in mitotic sister chromatid segregation, regulation of chromosome segregation and regulation of the cell cycle process. High CCNA2 expression was associated with poor overall survival (OS) and tumor stage. The expression of CDK1, CDC20, CCNB1, BUB1B, CCNA2, KIF11, CDCA8, KIF2C, NDC80 and TOP2A was increased in ovarian cancer tissues compared with normal tissues according to the Oncomine database. Higher expression levels of these seven candidate genes in ovarian cancer tissues compared with normal tissues were observed by GEPIA. The protein expression levels of CCNA2, CCNB1, CDC20, CDCA8, CDK1, KIF11 and TOP2A were high in ovarian cancer tissues, which was further confirmed via the HPA database. Conclusion Taken together, our study provided evidence concerning the altered expression of genes in ovarian cancer tissues compared with normal tissues. In vivo and in vitro experiments are required to verify the results of the present study.

To screen the key genes in the development of gastric cancer and their influence on prognosis. The GEO database was used to screen gastric cancer-related gene chips as a training set, and the R packages limma tool was used to analyze the differential genes expressed in gastric cancer tissues compared to normal tissues, and then the selected genes were verified in the validation set. The String database was used to calculate their Protein–protein interaction (PPI) network, using Cytoscape software's Centiscape and other plug-ins to analyze key genes in the PPI network. The DAVID database was used to enrich and annotate gene functions of differential genes and PPI key module genes, and further explore correlation between expression level and clinical stage and prognosis. Based on clinical data and patient samples, differential expression of key node genes was verified by immunohistochemistry. The 63 characteristic differential genes screened had good discrimination between gastric cancer and normal tissues, and are mainly involved in regulating extracellular matrix receptor interactions and the PI3k-AKT signaling pathway. Key nodes in the PPI network regulate tumor proliferation and metastasis. Analysis of the expression levels of key node genes found that relative to normal tissues, the expression of ITGB1 and COL1A2 was significantly increased in gastric cancer tissues, and patients with late clinical stages of tumors had higher expression of ITGB1 and COL1A2 in tumor tissues, and their survival time was longer (P

Background Gmelina arborea Roxb is a fast-growing tree species of commercial importance for tropical countries due to multiple industrial uses of its wood. Wood is primarily composed of thick secondary cell walls of xylem cells which imparts the strength to the wood. Identification of the genes involved in the secondary cell wall biosynthesis as well as their cognate regulators is crucial to understand how the production of wood occurs and serves as a starting point for developing breeding strategies to produce varieties with improved wood quality, better paper pulping or new potential uses such as biofuel production. In order to gain knowledge on the molecular mechanisms and gene regulation related with wood development in white teak, a de novo sequencing and transcriptome assembly approach was used employing secondary cell wall synthesizing cells from young white teak trees. Results For generation of transcriptome, RNA-seq reads were assembled into 110,992 transcripts and 49,364 genes were functionally annotated using plant databases; 5071 GO terms and 25,460 SSR markers were identified within xylem transcripts and 10,256 unigenes were assigned to KEGG database in 130 pathways. Among transcription factor families, C2H2, C3H, bLHLH and MYB were the most represented in xylem. Differential gene expression analysis using leaves as a reference was carried out and a total of 20,954 differentially expressed genes were identified including monolignol biosynthetic pathway genes. The differential expression of selected genes (4CL, COMT, CCoAOMT, CCR and NST1) was validated using qPCR. Conclusions We report the very first de novo transcriptome of xylem-related genes in this tropical timber species of commercial importance and constitutes a valuable extension of the publicly available transcriptomic resource aimed at fostering both basic and breeding studies.

Background Antimicrobial peptides (AMPs) are oligopeptides that act as crucial components of innate immunity, naturally occur in all multicellular organisms, and are involved in the first line of defense function. Recent studies showed that AMPs perpetuate great potential that is not limited to antimicrobial activity. They are also crucial regulators of host immune responses that can modulate a wide range of activities, such as immune regulation, wound healing, and apoptosis. However, a microorganism's ability to adapt and to resist existing antibiotics triggered the scientific community to develop alternatives to conventional antibiotics. Therefore, to address this issue, we proposed Co-AMPpred, an in silico-aided AMP prediction method based on compositional features of amino acid residues to classify AMPs and non-AMPs. Results In our study, we developed a prediction method that incorporates composition-based sequence and physicochemical features into various machine-learning algorithms. Then, the boruta feature-selection algorithm was used to identify discriminative biological features. Furthermore, we only used discriminative biological features to develop our model. Additionally, we performed a stratified tenfold cross-validation technique to validate the predictive performance of our AMP prediction model and evaluated on the independent holdout test dataset. A benchmark dataset was collected from previous studies to evaluate the predictive performance of our model. Conclusions Experimental results show that combining composition-based and physicochemical features outperformed existing methods on both the benchmark training dataset and a reduced training dataset. Finally, our proposed method achieved 80.8% accuracies and 0.871 area under the receiver operating characteristic curve by evaluating on independent test set. Our code and datasets are available at https://github.com/onkarS23/CoAMPpred.

Background Next-Generation-Sequencing (NGS) enables detection of microorganisms present in biological and other matrices of various origin and nature, allowing not only the identification of known phyla and strains but also the discovery of novel ones. The large amount of metagenomic shotgun data produced by NGS require comprehensive and user-friendly pipelines for data analysis, that speed up the bioinformatics steps, relieving the users from the need to manually perform complex and time-consuming tasks. Results We describe here HOME-BIO (sHOtgun MEtagenomic analysis of BIOlogical entities), an exhaustive pipeline for metagenomics data analysis, comprising three independent analytical modules designed for an inclusive analysis of large NGS datasets. Conclusions HOME-BIO is a powerful and easy-to-use tool that can be run also by users with limited computational expertise. It allows in-depth analyses by removing low-complexity/ problematic reads, integrating the analytical steps that lead to a comprehensive taxonomy profile of each sample by querying different source databases, and it is customizable according to specific users’ needs.

The application of array-based oligonucleotides in biological studies is described. These oligonucleotides are mainly used to design large libraries of various nucleotide sequences, which are applied to study protein-nucleic acid interactions, splicing, transcription, translation, and other regulatory processes in mammalian, yeast, and bacterial systems. The application of gene libraries generated by array-based nucleotides along with advanced methods of the combination of DNA duplexes will make it possible to obtain complex genetic designs for synthetic biology.

Background New breeding technologies (NBT) using CRISPR/Cas9-induced homology directed repair (HDR) has the potential to expedite genetic improvement in aquaculture. The long generation time in Atlantic salmon makes breeding an unattractive solution to obtain homozygous mutants and improving the rates of perfect HDR in founder (F0) fish is thus required. Genome editing can represent small DNA changes down to single nucleotide replacements (SNR). This enables edits such as premature stop codons or single amino acid changes and may be used to obtain fish with traits favorable to aquaculture, e.g. disease resistance. A method for SNR has not yet been demonstrated in salmon. Results Using CRISPR/Cas9 and asymmetrical ODNs, we were able to perform precise SNR and introduce a premature stop codon in dnd in F0 salmon. Deep sequencing demonstrated up to 59.2% efficiency in single embryos. In addition, using the same asymmetrical ODN design, we inserted a FLAG element into slc45a2 and dnd, showing high individual perfect HDR efficiencies (up to 36.7 and 32.7%, respectively). Conclusions In this work, we demonstrate that precise SNR and knock-in (KI) can be performed in F0 salmon embryos using asymmetrical oligonucleotide (ODN) donors. We suggest that HDR-induced SNR can be applied as a powerful NBT, allowing efficient introgression of favorable alleles and bypassing challenges associated with traditional selective breeding.

Background Tumors are composed by a number of cancer cell subpopulations (subclones), characterized by a distinguishable set of mutations. This phenomenon, known as intra-tumor heterogeneity (ITH), may be studied using Copy Number Aberrations (CNAs). Nowadays ITH can be assessed at the highest possible resolution using single-cell DNA (scDNA) sequencing technology. Additionally, single-cell CNA (scCNA) profiles from multiple samples of the same tumor can in principle be exploited to study the spatial distribution of subclones within a tumor mass. However, since the technology required to generate large scDNA sequencing datasets is relatively recent, dedicated analytical approaches are still lacking. Results We present PhyliCS, the first tool which exploits scCNA data from multiple samples from the same tumor to estimate whether the different clones of a tumor are well mixed or spatially separated. Starting from the CNA data produced with third party instruments, it computes a score, the Spatial Heterogeneity score, aimed at distinguishing spatially intermixed cell populations from spatially segregated ones. Additionally, it provides functionalities to facilitate scDNA analysis, such as feature selection and dimensionality reduction methods, visualization tools and a flexible clustering module. Conclusions PhyliCS represents a valuable instrument to explore the extent of spatial heterogeneity in multi-regional tumour sampling, exploiting the potential of scCNA data.

Assays for the molecular detection of nucleic acids are typically constrained by the level of multiplexing (this is the case for the quantitative polymerase chain reaction (qPCR) and for isothermal amplification), turnaround times (as with microarrays and next-generation sequencing), quantification accuracy (isothermal amplification, microarrays and nanopore sequencing) or specificity for single-nucleotide differences (microarrays and nanopore sequencing). Here we show that a portable and battery-powered PCR assay performed in a toroidal convection chamber housing a microarray of fluorescently quenched oligonucleotide probes allows for the rapid and sensitive quantification of multiple DNA targets with single-nucleotide discrimination. The assay offers a limit of detection of 10 DNA copies within 30 min of turnaround time and a dynamic range spanning 4 orders of magnitude of DNA concentration, and we show its performance by detecting 20 genomic loci and 30 single-nucleotide polymorphisms in human genomic DNA samples, and 15 bacterial species in clinical isolates. Portable devices for the fast and highly multiplexed detection of nucleic acids may offer advantages in point-of-care diagnostics. A portable PCR assay performed in a toroidal convection chamber housing an array of fluorescently quenched oligonucleotide probes allows for the rapid quantification of multiple DNA targets with single-nucleotide discrimination.

Background Identifying differentially expressed genes between the same or different species is an urgent demand for biological and medical research. For RNA-seq data, systematic technical effects and different sequencing depths are usually encountered when conducting experiments. Normalization is regarded as an essential step in the discovery of biologically important changes in expression. The present methods usually involve normalization of the data with a scaling factor, followed by detection of significant genes. However, more than one scaling factor may exist because of the complexity of real data. Consequently, methods that normalize data by a single scaling factor may deliver suboptimal performance or may not even work.The development of modern machine learning techniques has provided a new perspective regarding discrimination between differentially expressed (DE) and non-DE genes. However, in reality, the non-DE genes comprise only a small set and may contain housekeeping genes (in same species) or conserved orthologous genes (in different species). Therefore, the process of detecting DE genes can be formulated as a one-class classification problem, where only non-DE genes are observed, while DE genes are completely absent from the training data. Results In this study, we transform the problem to an outlier detection problem by treating DE genes as outliers, and we propose a scaling-free minimum enclosing ball (SFMEB) method to construct a smallest possible ball to contain the known non-DE genes in a feature space. The genes outside the minimum enclosing ball can then be naturally considered to be DE genes. Compared with the existing methods, the proposed SFMEB method does not require data normalization, which is particularly attractive when the RNA-seq data include more than one scaling factor. Furthermore, the SFMEB method could be easily extended to different species without normalization. Conclusions Simulation studies demonstrate that the SFMEB method works well in a wide range of settings, especially when the data are heterogeneous or biological replicates. Analysis of the real data also supports the conclusion that the SFMEB method outperforms other existing competitors. The R package of the proposed method is available at https://bioconductor.org/packages/MEB.

Background Recent scientific research has enabled the identification of macrophages related-genes (MaRG), which play a key role in the control of the immune microenvironment in many human cancers. However, the functional role of MaRGs in human tumors is ill-defined. Herein, we aimed at bioinformatically exploring the molecular signatures of MaRGs in colorectal cancer. Methods A list of MaRGs was generated and their differential expression was analyzed across multiple datasets downloaded from the publicly available functional genomics database Gene Expression Omnibus. The weighted gene co-expression network analysis (WGCNA) was also applied to identify the partner genes of these MaRGs in colorectal cancer. Results After integration of the results from analyses of different datasets, we found that 29 differentially expressed MaRGs (DE-MaRGs) could be considered as CRC-related genes as obtained from the WGCNA analysis. These genes were functionally involved in positive regulation of DNA biosynthetic process and glutathione metabolism. Protein–protein interaction network analysis indicated that PDIA6, PSMA1, PRC1, RRM2, HSP90AB1, CDK4, MCM7, RFC4, and CCT5 were the hub MaRGs. The LASSO approach was used for validating the 29 MaRGs in TCGA-COAD and TCGA-READ data and the results showed that ten among the 29 genes could be considered as MaRGs significantly involved in CRC. The maftools analysis showed that MaRGs were mutated at varying degrees. The nomogram analysis indicated the correlation of these MaRGs with diverse clinical features of CRC patients. Conclusions Conclusively, the present disclosed a signature of MaRGs as potential key regulators involved in CRC pathogenesis and progression. These findings contribute not only to the understanding of the molecular mechanism of CRC pathogenesis but also to the development of adequate immunotherapies for CRC patients.

Background Network models are well-established as very useful computational-statistical tools in cell biology. However, a genomic network model based only on gene expression data can, by definition, only infer gene co-expression networks. Hence, in order to infer gene regulatory patterns, it is necessary to also include data related to binding of regulatory factors to DNA. Results We propose a new dynamic genomic network model, for inferring patterns of genomic regulatory influence in dynamic processes such as development. Our model fuses experiment-specific gene expression data with publicly available DNA-binding data. The method we propose is computationally efficient, and can be applied to genome-wide data with tens of thousands of transcripts. Thus, our method is well suited for use as an exploratory tool for genome-wide data. We apply our method to data from human fetal cortical development, and our findings confirm genomic regulatory patterns which are recognised as being fundamental to neuronal development. Conclusions Our method provides a mathematical/computational toolbox which, when coupled with targeted experiments, will reveal and confirm important new functional genomic regulatory processes in mammalian development.