Search

Abstract Background Waardenburg syndrome (WS) is a rare genetic disorder mainly characterized by hearing loss and pigmentary abnormalities. Currently, seven causative genes have been identified for WS, but clinical genetic testing results show that 38.9% of WS patients remain molecularly unexplained. In this study, we performed multi-data integration analysis through protein-protein interaction and phenotype-similarity to comprehensively decipher the potential causative factors of undiagnosed WS. In addition, we explored the association between genotypes and phenotypes in WS with the manually collected 443 cases from published literature. Results We predicted two possible WS pathogenic genes (KIT, CHD7) through multi-data integration analysis, which were further supported by gene expression profiles in single cells and phenotypes in gene knockout mouse. We also predicted twenty, seven, and five potential WS pathogenic variations in gene PAX3, MITF, and SOX10, respectively. Genotype-phenotype association analysis showed that white forelock and telecanthus were dominantly present in patients with PAX3 variants; skin freckles and premature graying of hair were more frequently observed in cases with MITF variants; while aganglionic megacolon and constipation occurred more often in those with SOX10 variants. Patients with variations of PAX3 and MITF were more likely to have synophrys and broad nasal root. Iris pigmentary abnormality was more common in patients with variations of PAX3 and SOX10. Moreover, we found that patients with variants of SOX10 had a higher risk of suffering from auditory system diseases and nervous system diseases, which were closely associated with the high expression abundance of SOX10 in ear tissues and brain tissues. Conclusions Our study provides new insights into the potential causative factors of WS and an alternative way to explore clinically undiagnosed cases, which will promote clinical diagnosis and genetic counseling. However, the two potential disease-causing genes (KIT, CHD7) and 32 potential pathogenic variants (PAX3: 20, MITF: 7, SOX10: 5) predicted by multi-data integration in this study are all computational predictions and need to be further verified through experiments in follow-up research.

Abstract Accurate indel calling plays an important role in precision medicine. A benchmarking indel set is essential for thoroughly evaluating the indel calling performance of bioinformatics pipelines. A reference sample with a set of known-positive variants was developed in the FDA-led Sequencing Quality Control Phase 2 (SEQC2) project, but the known indels in the known-positive set were limited. This project sought to provide an enriched set of known indels that would be more translationally relevant by focusing on additional cancer related regions. A thorough manual review process completed by 42 reviewers, two advisors, and a judging panel of three researchers significantly enriched the known indel set by an additional 516 indels. The extended benchmarking indel set has a large range of variant allele frequencies (VAFs), with 87% of them having a VAF below 20% in reference Sample A. The reference Sample A and the indel set can be used for comprehensive benchmarking of indel calling across a wider range of VAF values in the lower range. Indel length was also variable, but the majority were under 10 base pairs (bps). Most of the indels were within coding regions, with the remainder in the gene regulatory regions. Although high confidence can be derived from the robust study design and meticulous human review, this extensive indel set has not undergone orthogonal validation. The extended benchmarking indel set, along with the indels in the previously published known-positive set, was the truth set used to benchmark indel calling pipelines in a community challenge hosted on the precisionFDA platform. This benchmarking indel set and reference samples can be utilized for a comprehensive evaluation of indel calling pipelines. Additionally, the insights and solutions obtained during the manual review process can aid in improving the performance of these pipelines.

Abstract Background Acid sphingomyelinase deficiency (ASMD) disorder, also known as Niemann–Pick disease (NPD) is a rare genetic disease caused by mutations in SMPD1 gene, which encodes sphingomyelin phosphodiesterase (ASM). Except for liver and spleen enlargement and lung disease, two subtypes (Type A and B) of NDP have different onset times, survival times, ASM activities, and neurological abnormalities. To comprehensively explore NPD’s genotype-phenotype association and pathophysiological characteristics, we collected 144 NPD cases with strict quality control through literature mining. Results The difference in ASM activity can differentiate NPD type A from other subtypes, with the ratio of ASM activity to the reference values being lower in type A (threshold 0.045 (4.45%)). Severe variations, such as deletion and insertion, can cause complete loss of ASM function, leading to type A, whereas relatively mild missense mutations generally result in type B. Among reported mutations, the p.Arg3AlafsX76 mutation is highly prevalent in the Chinese population, and the p.R608del mutation is common in Mediterranean countries. The expression profiles of SMPD1 from GTEx and single-cell RNA sequencing data of multiple fetal tissues showed that high expressions of SMPD1 can be observed in the liver, spleen, and brain tissues of adults and hepatoblasts, hematopoietic stem cells, STC2_TLX1-positive cells, mesothelial cells of the spleen, vascular endothelial cells of the cerebellum and the cerebrum of fetuses, indicating that SMPD1 dysfunction is highly likely to have a significant effect on the function of those cell types during development and the clinicians need pay attention to these organs or tissues as well during diagnosis. In addition, we also predicted 21 new pathogenic mutations in the SMPD1 gene that potentially cause the NPD, signifying that more rare cases will be detected with those mutations in SMPD1. Finally, we also analysed the function of the NPD type A cells following the extracellular milieu. Conclusions Our study is the first to elucidate the effects of SMPD1 mutation on cell types and at the tissue level, which provides new insights into the genotype-phenotype association and can help in the precise diagnosis of NPD.

Abstract Background Two main subclasses of macrophages are found in almost all solid tissues: embryo-derived resident tissue macrophages and bone marrow-derived infiltrated macrophages. These macrophage subtypes show transcriptional and functional divergence, and the programs that have shaped the evolution of renal macrophages and related signaling pathways remain poorly understood. To clarify these processes, we performed data analysis based on single-cell transcriptional profiling of renal tissue-resident and infiltrated macrophages in human, mouse and rat. Results In this study, we (i) characterized the transcriptional divergence among species and (ii) illustrated variability in expression among cells of each subtype and (iii) compared the gene regulation network and (iv) ligand-receptor pairs in human and mouse. Using single-cell transcriptomics, we mapped the promoter architecture during homeostasis. Conclusions Transcriptionally divergent genes, such as the differentially TF-encoding genes expressed in resident and infiltrated macrophages across the three species, vary among cells and include distinct promoter structures. The gene regulatory network in infiltrated macrophages shows comparatively better species-wide consistency than resident macrophages. The conserved transcriptional gene regulatory network in infiltrated macrophages among species is uniquely enriched in pathways related to kinases, and TFs associated with largely conserved regulons among species are uniquely enriched in kinase-related pathways.

Cholesterol is an important precursor of many endogenous molecules. Disruption of cholesterol homeostasis can cause many pathological changes, leading to liver and cardiovascular diseases. CYP1A is widely involved in cholesterol metabolic network, but its exact function has not been fully elucidated. Here, we aim to explore how CYP1A regulates cholesterol homeostasis. Our data showed that CYP1A1/2 knockout (KO) rats presented cholesterol deposition in blood and liver. The serum levels of low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and total cholesterol were significantly increased in KO rats. Further studies found that the lipogenesis pathway (LXRα–SREBP1–SCD1) of KO rats was activated, and the key protein of cholesterol ester hydrolysis (CES1) was inhibited. Importantly, lansoprazole can significantly alleviate rat hepatic lipid deposition in hypercholesterolemia models by inducing CYP1A. Our findings reveal the role of CYP1A as a potential regulator of cholesterol homeostasis and provide a new perspective for the treatment of hypercholesterolemia.

BackgroundMelanoma is a skin tumor with a high mortality rate, and early diagnosis and effective treatment are the key to reduce its mortality rate. Therefore, more and more attention has been paid for biomarker identification for early diagnosis, prognosis prediction and prognosis evaluation of melanoma. However, there is still a lack of a report that comprehensively and objectively evaluates the research status of melanoma biomarkers. Therefore, this study aims to intuitively analyze the research status and trend of melanoma biomarkers through the methods of bibliometrics and knowledge graph.ObjectiveThis study uses bibliometrics to analyze research in biomarkers in melanoma, summarize the field’s history and current status of research, and predict future research directions.MethodArticles and Reviews related to melanoma biomarkers were retrieved by using Web of Science core collection subject search. Bibliometric analysis was performed in Excel 365, CiteSpace, VOSviewer and Bibliometrix (R-Tool of R-Studio).ResultA total of 5584 documents from 2004 to 2022 were included in the bibliometric analysis. The results show that the number of publications and the frequency of citations in this field are increasing year by year, and the frequency of citations has increased rapidly after 2018. The United States is the most productive and influential country in this field, with the largest number of publications and institutions with high citation frequency. Caroline Robert, F. Stephen Hodi, Suzanne L. Topalian and others are authoritative authors in this field, and The New England Journal of Medicine, Journal of Clinical Oncology and Clinical Cancer Research are the most authoritative journals in this field. Biomarkers related to the diagnosis, treatment and prognosis of melanoma are hot topics and cutting-edge hotspots in this field.ConclusionFor the first time, this study used the bibliometric method to visualize the research in the field of melanoma biomarkers, revealing the trends and frontiers of melanoma biomarkers research, which provides a useful reference for scholars to find key research issues and partners.

An increasing number of studies have reported that microbiome can affect drug response by altering pharmacokinetics and pharmacodynamics of formation of toxic metabolites. With the development of metagenomic sequencing, gut microbial composition as well as the metabolic function are drawing more and more attention for the patient stratification. The established microbiota databases provide useful information about the gut microbe-drug interactions. However, these databases generally lacked the detailed effects on substance and the metabolites, which are helpful in elucidating the mechanisms underlying drug biotransformation and personalized medicine. To address these issues, in this study, we developed Metabolic action of gut Microbiota to Drugs (MagMD), a database and a web-service covering 32, 678 records of interactions between 2,146 gut microbes, 36 enzymes and 219 substrates (mainly drugs). The detailed annotations for each entry, including the taxonomic level of microbes, the molecular form and PubChem ID of drugs from PubChem Compound Database, types of microbial secreted enzymes and the original reference links can also be accessed from the web service. Availability and implementation: MagMD is a publicly available resource, constantly updated. It has an intuitive web interface and can be freely accessed at http://www.unimd.org/magmd.

Abstract Background Reproducible detection of inherited variants with whole genome sequencing (WGS) is vital for the implementation of precision medicine and is a complicated process in which each step affects variant call quality. Systematically assessing reproducibility of inherited variants with WGS and impact of each step in the process is needed for understanding and improving quality of inherited variants from WGS. Results To dissect the impact of factors involved in detection of inherited variants with WGS, we sequence triplicates of eight DNA samples representing two populations on three short-read sequencing platforms using three library kits in six labs and call variants with 56 combinations of aligners and callers. We find that bioinformatics pipelines (callers and aligners) have a larger impact on variant reproducibility than WGS platform or library preparation. Single-nucleotide variants (SNVs), particularly outside difficult-to-map regions, are more reproducible than small insertions and deletions (indels), which are least reproducible when > 5 bp. Increasing sequencing coverage improves indel reproducibility but has limited impact on SNVs above 30×. Conclusions Our findings highlight sources of variability in variant detection and the need for improvement of bioinformatics pipelines in the era of precision medicine with WGS.

Prostate cancer (PCa) is the most common malignancy. New biomarkers are in demand to facilitate the management. The role of the pinin protein (encoded by PNN gene) in PCa has not been thoroughly explored yet. Using The Cancer Genome Atlas (TCGA-PCa) dataset validated with Gene Expression Omnibus (GEO) and protein expression data retrieved from the Human Protein Atlas, the prognostic and diagnostic values of PNN were studied. Highly co-expressed genes with PNN (HCEG) were constructed for pathway enrichment analysis and drug prediction. A prognostic signature based on methylation status using HCEG was constructed. Gene set enrichment analysis (GSEA) and the TISIDB database were utilised to analyse the associations between PNN and tumour-infiltrating immune cells. The upregulated PNN expression in PCa at both transcription and protein levels suggests its potential as an independent prognostic factor of PCa. Analyses of the PNN’s co-expression network indicated that PNN plays a role in RNA splicing and spliceosomes. The prognostic methylation signature demonstrated good performance for progression-free survival. Finally, our results showed that the PNN gene was involved in splicing-related pathways in PCa and identified as a potential biomarker for PCa.

Summary: The MetaSUB Consortium, founded in 2015, is a global consortium with an interdisciplinary team of clinicians, scientists, bioinformaticians, engineers, and designers, with members from more than 100 countries across the globe. This network has continually collected samples from urban and rural sites including subways and transit systems, sewage systems, hospitals, and other environmental sampling. These collections have been ongoing since 2015 and have continued when possible, even throughout the COVID-19 pandemic. The consortium has optimized their workflow for the collection, isolation, and sequencing of DNA and RNA collected from these various sites and processing them for metagenomics analysis, including the identification of SARS-CoV-2 and its variants. Here, the Consortium describes its foundations, and its ongoing work to expand on this network and to focus its scope on the mapping, annotation, and prediction of emerging pathogens, mapping microbial evolution and antibiotic resistance, and the discovery of novel organisms and biosynthetic gene clusters.

Abstract Background Gene copy number variations (CNVs) contribute to genetic diversity and disease prevalence across populations. Substantial efforts have been made to decipher the relationship between CNVs and pathogenesis but with limited success. Results We have developed a novel computational framework X-CNV ( www.unimd.org/XCNV ), to predict the pathogenicity of CNVs by integrating more than 30 informative features such as allele frequency (AF), CNV length, CNV type, and some deleterious scores. Notably, over 14 million CNVs across various ethnic groups, covering nearly 93% of the human genome, were unified to calculate the AF. X-CNV, which yielded area under curve (AUC) values of 0.96 and 0.94 in training and validation sets, was demonstrated to outperform other available tools in terms of CNV pathogenicity prediction. A meta-voting prediction (MVP) score was developed to quantitively measure the pathogenic effect, which is based on the probabilistic value generated from the XGBoost algorithm. The proposed MVP score demonstrated a high discriminative power in determining pathogenetic CNVs for inherited traits/diseases in different ethnic groups. Conclusions The ability of the X-CNV framework to quantitatively prioritize functional, deleterious, and disease-causing CNV on a genome-wide basis outperformed current CNV-annotation tools and will have broad utility in population genetics, disease-association studies, and diagnostic screening.

Transcriptome profiling of hepatocellular carcinoma (HCC) by next-generation sequencing (NGS) technology has been broadly performed by previous studies, which facilitate our understanding of the molecular mechanisms of HCC formation, progression, and metastasis. However, few studies jointly analyze multiple types of noncoding RNAs (ncRNAs), including long noncoding RNAs (lncRNAs), circular RNAs (circRNAs), and micro-RNAs (miRNAs), and further uncover their implications in HCC. In this study, we observed that the circRNA cZRANB1 and lncRNA DUXAP10 were not only significantly upregulated in tumor tissues, but also higher expressed in blood exosomes of HCC as compared with healthy donors. From the analysis of subclass-associated dysregulated ncRNAs, we observed that DLX6-AS1, an antisense RNA of DLX6, and the sense gene DLX6 were highly expressed in S1, a subclass with a more invasive/disseminative phenotype. High correlation between DLX6-AS1 and DLX6 suggested that DLX6-AS1 may function via promoting the transcription of DLX6. Integrative analysis uncovers circRNA–miRNA, lncRNA–miRNA, and competing endogenous RNA networks (ceRNAs). Specifically, cZRANB1, LINC00501, CTD-2008L17.2, and SLC7A11-AS1 may function as ceRNAs that regulate mRNAs by competing the shared miRNAs. Further prognostic analysis demonstrated that the dysregulated ncRNAs had the potential to predict HCC patients’ overall survival. In summary, we identified some novel circRNAs and miRNAs, and dysregulated ncRNAs that could participate in HCC tumorigenesis and progression by inducing transcription of their neighboring genes, increasing their derived miRNAs, or acting as miRNA sponges. Moreover, our systematic analysis provides not only rich data resources for related researchers, but also new insights into the molecular basis of how different ncRNAs coordinately or antagonistically participate in the pathogenesis process of diseases.

Abstract Background Oncopanel genomic testing, which identifies important somatic variants, is increasingly common in medical practice and especially in clinical trials. Currently, there is a paucity of reliable genomic reference samples having a suitably large number of pre-identified variants for properly assessing oncopanel assay analytical quality and performance. The FDA-led Sequencing and Quality Control Phase 2 (SEQC2) consortium analyze ten diverse cancer cell lines individually and their pool, termed Sample A, to develop a reference sample with suitably large numbers of coding positions with known (variant) positives and negatives for properly evaluating oncopanel analytical performance. Results In reference Sample A, we identify more than 40,000 variants down to 1% allele frequency with more than 25,000 variants having less than 20% allele frequency with 1653 variants in COSMIC-related genes. This is 5–100× more than existing commercially available samples. We also identify an unprecedented number of negative positions in coding regions, allowing statistical rigor in assessing limit-of-detection, sensitivity, and precision. Over 300 loci are randomly selected and independently verified via droplet digital PCR with 100% concordance. Agilent normal reference Sample B can be admixed with Sample A to create new samples with a similar number of known variants at much lower allele frequency than what exists in Sample A natively, including known variants having allele frequency of 0.02%, a range suitable for assessing liquid biopsy panels. Conclusion These new reference samples and their admixtures provide superior capability for performing oncopanel quality control, analytical accuracy, and validation for small to large oncopanels and liquid biopsy assays.

Abstract Background Targeted sequencing using oncopanels requires comprehensive assessments of accuracy and detection sensitivity to ensure analytical validity. By employing reference materials characterized by the U.S. Food and Drug Administration-led SEquence Quality Control project phase2 (SEQC2) effort, we perform a cross-platform multi-lab evaluation of eight Pan-Cancer panels to assess best practices for oncopanel sequencing. Results All panels demonstrate high sensitivity across targeted high-confidence coding regions and variant types for the variants previously verified to have variant allele frequency (VAF) in the 5–20% range. Sensitivity is reduced by utilizing VAF thresholds due to inherent variability in VAF measurements. Enforcing a VAF threshold for reporting has a positive impact on reducing false positive calls. Importantly, the false positive rate is found to be significantly higher outside the high-confidence coding regions, resulting in lower reproducibility. Thus, region restriction and VAF thresholds lead to low relative technical variability in estimating promising biomarkers and tumor mutational burden. Conclusion This comprehensive study provides actionable guidelines for oncopanel sequencing and clear evidence that supports a simplified approach to assess the analytical performance of oncopanels. It will facilitate the rapid implementation, validation, and quality control of oncopanels in clinical use.

Summary: Patients with cancer with different molecular characterization and subtypes result in different response to anticancer therapeutics and survival. To identify features that are associated with prognosis is essential to precision medicine by providing clues for target identification, drug discovery. Here, we developed a tumor online prognostic analysis platform (ToPP) which integrated eight multi-omics features and clinical data from 68 cancer projects. It provides multiple approaches for customized prognostic studies, including 1) Prognostic analysis based on multi-omics features and clinical characteristics; 2) Automatic construction of prognostic model; 3) Pancancer prognostic analysis in multi-omics data; 4) Explore the impact of different levels of feature combinations on patient prognosis; 5) More sophisticated prognostic analysis according to regulatory network. ToPP provides a comprehensive source and easy-to-use interface for tumor prognosis research, with one-stop service of multi-omics, subtyping, and online prognostic modeling. The web server is freely available at http://www.biostatistics.online/topp/index.php.

Abstract Background Berardinelli-Seip congenital lipodystrophy (BSCL) is a heterogeneous autosomal recessive disorder characterized by an almost total lack of adipose tissue in the body. Mutations in the AGPAT2, BSCL2, CAV1 and PTRF genes define I-IV subtype of BSLC respectively and clinical data indicate that new causative genes remain to be discovered. Here, we retrieved 341 cases from 60 BSCL-related studies worldwide and aimed to explore genotype-phenotype correlations based on mutations of AGPAT2 and BSCL2 genes from 251 cases. We also inferred new candidate genes for BSCL through protein-protein interaction and phenotype-similarity. Results Analysis results show that BSCL type II with earlier age of onset of diabetes mellitus, higher risk to suffer from premature death and mental retardation, is a more severe disorder than BSCL type I, but BSCL type I patients are more likely to have bone cysts. In BSCL type I, females are at higher risk of developing diabetes mellitus and acanthosis nigricans than males, while in BSCL type II, males suffer from diabetes mellitus earlier than females. In addition, some significant correlations among BSCL-related phenotypes were identified. New candidate genes prediction through protein-protein interaction and phenotype-similarity was conducted and we found that CAV3, EBP, SNAP29, HK1, CHRM3, OBSL1 and DNAJC13 genes could be the pathogenic factors for BSCL. Particularly, CAV3 and EBP could be high-priority candidate genes contributing to pathogenesis of BSCL. Conclusions Our study largely enhances the current knowledge of phenotypic and genotypic heterogeneity of BSCL and promotes the more comprehensive understanding of pathogenic mechanisms for BSCL.

Pulmonary atresia (PA) is a rare congenital heart defect (CHD) with complex manifestations and a high mortality rate. Since the genetic determinants in the pathogenesis of PA remain elusive, a thorough identification of the genetic factors through whole exome sequencing (WES) will provide novel insights into underlying mechanisms of PA. We performed WES data from PA/VSD (n = 60), PA/IVS (n = 20), TOF/PA (n = 20) and 100 healthy controls. Rare variants and novel genes were identified using variant-based association and gene-based burden analysis. Then we explored the expression pattern of our candidate genes in endothelium cell lines, pulmonary artery tissues, and embryonic hearts. 56 rare damage variants of 7 novel candidate genes (DNAH10, DST, FAT1, HMCN1, HNRNPC, TEP1, and TYK2) were certified to have function in PA pathogenesis for the first time. In our research, the genetic pattern among PA/VSD, PA/IVS and TOF/PA were different to some degree. Taken together, our findings contribute new insights into the molecular basis of this rare congenital birth defect.

Lung adenocarcinoma (LUAD) is the most common subtype of lung cancer with heterogeneous outcomes and diverse therapeutic responses. To classify patients into different groups and facilitate the suitable therapeutic strategy, we first selected eight microRNA (miRNA) signatures in The Cancer Genome Atlas (TCGA)-LUAD cohort based on multi-strategy combination, including differential expression analysis, regulatory relationship, univariate survival analysis, importance clustering, and multivariate combinations analysis. Using the eight miRNA signatures, we further built novel risk scores based on the predefined cutoff and beta coefficients and divided the patients into high-risk and low-risk groups with significantly different overall survival time (p-value < 2 e−16). The risk-score model was confirmed with an independent dataset (p-value = 4.71 e−4). We also observed that the risk scores of early-stage patients were significantly lower than those of late-stage patients. Moreover, our model can also provide new insights into the current clinical staging system and can be regarded as an alternative system for patient stratification. This model unified the variable value as the beta coefficient facilitating the integration of biomarkers obtained from different omics data.

ObjectivesTo investigate the genetic variants that are responsible for peripheral neuroblastic tumors (PNTs) oncogenesis in one family case.Materials and MethodsOne family was recruited, including the healthy parents, sister affected by neuroblastoma (NB), and brother who suffered from ganglioneuroma (GN). Whole-genome sequencing (WGS) of germline DNA from all the family members and RNA-seq of tumor RNA from the siblings were performed. Mutants were validated by Sanger sequencing and co-IP was performed to assess the impact of the mutant on chemosensitivity in the SH-SY5Y cell line.ResultsA novel compound heterozygous mutation of BRCA2 was locked as the cause of carcinogenesis. One allele was BRCA2-S871X (stop-gain) from the siblings’ mother, the other was BRCA2-N372H (missense) from their father. This novel compound heterozygous mutations of the BRCA2 gene associated with PNTs by disordering DNA damage and response (DDR) signal pathway. Moreover, chemosensitivity was reduced in the NB cell line due to the BRCA2-N372H mutant.ConclusionIn summary, these results revealed a novel germline compound heterozygous mutation of the BRCA2 gene associated with familial PNTs.

Summary: The diffuse-type gastric cancer (DGC) constitutes a subgroup of gastric cancer with poor prognosis and no effective molecular therapies. Here, we report a phosphoproteomic landscape of DGC derived from 83 tumors together with their nearby tissues. Based on phosphorylation, DGC could be classified into three molecular subtypes with distinct overall survival (OS) and chemosensitivity. We identified 16 kinases whose activities were associated with poor OS. These activated kinases covered several cancer hallmark pathways, with the MTOR signaling network being the most frequently activated. We proposed a patient-specific strategy based on the hierarchy of clinically actionable kinases for prioritization of kinases for further clinical evaluation. Our global data analysis indicates that in addition to finding activated kinase pathways in DGC, large-scale phosphoproteomics could be used to classify DGCs into subtypes that are associated with distinct clinical outcomes as well as nomination of kinase targets that may be inhibited for cancer treatments. : Biological Sciences; Cancer Systems Biology; Proteomics; Systems Biology Subject Areas: Biological Sciences, Cancer Systems Biology, Proteomics, Systems Biology

Abstract Importance Morbidity and mortality of children are important indicators of the performance of the public health system in any country. In China, the children's disease spectrum has gradually changed in recent years. However, the gender‐ and age‐specific disease spectrum for hospitalized children under 15 years old is still unclear. Objective To explore the gender‐ and age‐based distribution of diseases in hospitalized children under 15 years in China. Methods Medical records home page data for 2016 to 2018 were collected from 18 tertiary children's hospitals in China. The gender‐ and age‐specific disease spectrum was analyzed, using the tenth revision of the International Statistical Classification of Diseases and Related Health Problems (ICD‐10). Results The most common diseases were those of the respiratory system (25.7% of all 2 232 142 hospitalized children). The top three diseases for boys were diseases of the respiratory system (25.6%), diseases of the digestive system (11.4%) and certain conditions originating in the perinatal period (8.6%). The top three diseases for girls were diseases of the respiratory system (25.9%), certain conditions originating in the perinatal period (10.1%), and factors influencing health status and contact with health services (9.4%). The most common diseases for children under 1 year old were certain conditions originating in the perinatal period (38.1%). For all other age groups, the most common conditions were respiratory diseases (33.8% for those aged 1–3 years, 25.2% for those aged 4–6 years, and 12.2% for those aged 7–14 years). Interpretation This study analyzed the medical records home pages of 18 children's hospitals to provide the first overview of the disease spectrum and its gender‐ and age‐specific distribution among children in China.

Abstract Background Glioma is the most commonly diagnosed malignant and aggressive brain cancer in adults. Traditional researches mainly explored the expression profile of glioma at cell-population level, but ignored the heterogeneity and interactions of among glioma cells. Methods Here, we firstly analyzed the single-cell RNA-seq (scRNA-seq) data of 6341 glioma cells using manifold learning and identified neoplastic and healthy cells infiltrating in tumor microenvironment. We systematically revealed cell-to-cell interactions inside gliomas based on corresponding scRNA-seq and TCGA RNA-seq data. Results A total of 16 significantly correlated autocrine ligand-receptor signal pairs inside neoplastic cells were identified based on the scRNA-seq and TCGA data of glioma. Furthermore, we explored the intercellular communications between cancer stem-like cells (CSCs) and macrophages, and identified 66 ligand-receptor pairs, some of which could significantly affect prognostic outcomes. An efficient machine learning model was constructed to accurately predict the prognosis of glioma patients based on the ligand-receptor interactions. Conclusion Collectively, our study not only reveals functionally important cell-to-cell interactions inside glioma, but also detects potentially prognostic markers for predicting the survival of glioma patients.

DNA methylation status is closely associated with diverse diseases, and is generally more stable than gene expression, thus abnormal DNA methylation could be important biomarkers for tumor diagnosis, treatment and prognosis. However, the signatures regarding DNA methylation changes for pan-cancer diagnosis and prognosis are less explored. Here we systematically analyzed the genome-wide DNA methylation patterns in diverse TCGA cancers with machine learning. We identified seven CpG sites that could effectively discriminate tumor samples from adjacent normal tissue samples for 12 main cancers of TCGA (1216 samples, AUC > 0.99). Those seven potential diagnostic biomarkers were further validated in the other 9 different TCGA cancers and 4 independent datasets (AUC > 0.92). Three out of the seven CpG sites were correlated with cell division, DNA replication and cell cycle. We also identified 12 CpG sites that can effectively distinguish 26 different cancers (7605 samples), and the result was repeatable in independent datasets as well as two disparate tumors with metastases (micro-average AUC > 0.89). Furthermore, a series of potential signatures that could significantly predict the prognosis of tumor patients for 7 different cancer were identified via survival analysis (p-value < 1e-4). Collectively, DNA methylation patterns vary greatly between tumor and adjacent normal tissues, as well as among different types of cancers. Our identified signatures may aid the decision of clinical diagnosis and prognosis for pan-cancer and the potential cancer-specific biomarkers could be used to predict the primary site of metastatic breast and prostate cancers.

The human stomach is divided into seven anatomically distinct regions but their protein composition is largely unknown. Here, the authors present a region-resolved map of the healthy human stomach mucosa as well as mucosa proteomes of tumor and tumor nearby tissue from gastric cancer patients.

Neuroblastoma (NB) is a rare type of cancer but frequently occurred in children. However, it is still unclear whether circular RNAs (circRNAs) play key roles in NB tumorigenesis or progression. In this study, we identified 39,022 circRNAs across the 39 neuroblastoma and 2 normal cell lines. With the gene and circRNA expression data, we classified the NB cell lines, identified and characterized the functional circRNAs in the 3 NB classes. Specifically, 29 circRNAs were found to be dysregulated in the NB classes. Notably, 7 circRNAs located within MYCN-amplified regions were upregulated in cell lines with the high activities of MYC targets and MYCN amplification, and were highly correlated with expression of their parental gene, NBAS . Subsequently, we constructed ceRNA networks for the functional circRNAs. Specifically, hsa_circ_0005379 was identified as a critical regulator in the ceRNA networks because of targeting 13 genes, which formed a complex competing endogenous RNA (ceRNA) network. Moreover, hsa_circ_0002343, which was connected with few genes, might regulate the PI3K/Akt/mTOR signaling via RAC1. Furthermore, 3 genes, including NOTCH2, SERPINH1, and LAMC1, involved in epithelial mesenchymal transition (EMT) were observed to connect with hsa_circ_0001361, suggesting that this circRNA was closely associated with EMT. Consequently, 7 genes, such as DAD1, PPIA, NOTCH2, PGK1, BUB1, EIF2S1, and TCF7L2, were found to be closely associated with both event-free survival (EFS) and overall survival (OS). In conclusion, the present study identified functional circRNAs and predicted their functionality in neuroblastoma cell lines, which not only improved the understanding of circRNAs in neuroblastoma, but also provided the evidences for the related researchers.

Co-occurrence and mutual exclusivity (COME) of DNA methylation refer to two or more genes that tend to be positively or negatively correlated in DNA methylation among different samples. Although COME of gene mutations in pan-cancer have been well explored, little is known about the COME of DNA methylation in pan-cancer. Here, we systematically explored the COME of DNA methylation profile in diverse human cancer. A total of 5,128,332 COME events were identified in 14 main cancers types in The Cancer Genome Atlas (TCGA). We also identified functional epigenetic modules of the zinc finger gene family in six cancer types by integrating the gene expression and DNA methylation data and the frequently occurred COME network. Interestingly, most of the genes in those functional epigenetic modules are epigenetically repressed. Strikingly, those frequently occurred COME events could be used to classify the patients into several subtypes with significant different clinical outcomes in six cancers as well as pan-cancer (p-value ≤ = 0.05). Moreover, we observed significant associations between different COME subtypes and clinical features (e.g., age, gender, histological type, neoplasm histologic grade, and pathologic stage) in distinct cancers. Taken together, we identified millions of COME events of DNA methylation in pan-cancer and detected functional epigenetic COME events that could separate tumor patients into different subtypes, which may benefit the diagnosis and prognosis of pan-cancer.

Diffuse-type gastric cancer (DGC) accounts for 30% of gastric cancers and has few treatment options. Here the authors present a mutation and proteome dataset for 84 patients, identifying three major classes of DGC and indicating potential targets for therapy.

Abstract Environmental chemicals exposure is one of the primary factors for liver toxicity and hepatocarcinoma. Thioacetamide (TAA) is a well-known hepatotoxicant and could be a liver carcinogen in humans. The discovery of early and sensitive microRNA (miRNA) biomarkers in liver injury and tumor progression could improve cancer diagnosis, prognosis, and management. To study this, we performed next generation sequencing of the livers of Sprague-Dawley rats treated with TAA at three doses (4.5, 15 and 45 mg/kg) and four time points (3-, 7-, 14- and 28-days). Overall, 330 unique differentially expressed miRNAs (DEMs) were identified in the entire TAA-treatment course. Of these, 129 DEMs were found significantly enriched for the “liver cancer” annotation. These results were further complemented by pathway analysis (Molecular Mechanisms of Cancer, p53-, TGF-β-, MAPK- and Wnt-signaling). Two miRNAs (rno-miR-34a-5p and rno-miR-455-3p) out of 48 overlapping DEMs were identified to be early and sensitive biomarkers for TAA-induced hepatocarcinogenicity. We have shown significant regulatory associations between DEMs and TAA-induced liver carcinogenesis at an earlier stage than histopathological features. Most importantly, miR-34a-5p is the most suitable early and sensitive biomarker for TAA-induced hepatocarcinogenesis due to its consistent elevation during the entire treatment course.

Urine as a true non-invasive sampling source holds great potential for biomarker discovery. While approximately 2000 proteins can be detected by mass spectrometry in urine from healthy people, the amount of these proteins vary considerably. A systematic evaluation of a large number of samples is needed to determine the range of the variations. Current biomarker studies often measure limited number of urine samples in the discovery phase, which makes it difficult to determine whether proteins differentially expressed between control and disease groups represent actual difference, or are just physiological variations among the individuals, leads to failures in the validation phase with the increased sample numbers. Here, we report a streamlined workflow with capacity of measuring 8 urine proteomes per day at the coverage of >1500 proteins. With this workflow, we evaluated variations in 497 urine proteomes from 167 healthy donors, establishing reference intervals (RIs) that covered urine protein variations. We demonstrated that RIs could be used to monitor physiological changes by detecting transient outlier proteins. Furthermore, we provided a RIs-based algorithm for biomarker discovery and validation to screen for diseases such as cancer. This study provided a proof-of-principle workflow for the use of urine proteome for health monitoring and disease screening.

Secreted proteins (SPs) play important roles in diverse important biological processes; however, a comprehensive and high-quality list of human SPs is still lacking. Here we identified 6,943 high-confidence human SPs (3,522 of them are novel) based on 330,427 human proteins derived from databases of UniProt, Ensembl, AceView, and RefSeq. Notably, 6,267 of 6,943 (90.3%) SPs have the supporting evidences from a large amount of mass spectrometry (MS) and RNA-seq data. We found that the SPs were broadly expressed in diverse tissues as well as human body fluid, and a significant portion of them exhibited tissue-specific expression. Moreover, 14 cancer-specific SPs that their expression levels were significantly associated with the patients’ survival of eight different tumors were identified, which could be potential prognostic biomarkers. Strikingly, 89.21% of 6,943 SPs (2,927 novel SPs) contain known protein domains. Those novel SPs we mainly enriched with the known domains regarding immunity, such as Immunoglobulin V-set and C1-set domain. Specifically, we constructed a user-friendly and freely accessible database, SPRomeDB (www.unimd.org/SPRomeDB), to catalog those SPs. Our comprehensive SP identification and characterization gain insights into human secretome and provide valuable resource for future researches.

Osteogenesis imperfecta (OI), mainly caused by structural abnormalities of type I collagen, is a hereditary rare disease characterized by increased bone fragility and reduced bone mass. Clinical manifestations of OI mostly include multiple repeated bone fractures, thin skin, blue sclera, hearing loss, cardiovascular and pulmonary system abnormalities, triangular face, dentinogenesis imperfecta (DI), and walking with assistance. Currently, 20 causative genes with 18 subtypes have been identified for OI, of them, variations in COL1A1 and COL1A2 have been demonstrated to be major causative factors to OI. However, the complexity of the bone formation process indicates that there are potential new pathogenic genes associated with OI. To comprehensively explore the underlying mechanism of OI, we conducted association analysis between genotypes and phenotypes of OI diseases and found that mutations in COL1A1 and COL1A2 contributed to a large proportion of the disease phenotypes. We categorized the clinical phenotypes and the genotypes based on the variation types for those 155 OI patients collected from literature, and association study revealed that three phenotypes (bone deformity, DI, walking with assistance) were enriched in two variation types (the Gly-substitution missense and groups of frameshift, nonsense, and splicing variations). We also identified four novel variations (c.G3290A (p.G1097D), c.G3289C (p.G1097R), c.G3289A (p.G1097S), c.G3281A (p.G1094D)) in gene COL1A1 and two novel variations (c.G2332T (p.G778C), c.G2341T (p.G781C)) in gene COL1A2, which could potentially contribute to the disease. In addition, we identified several new potential pathogenic genes (ADAMTS2, COL5A2, COL8A1) based on the integration of protein–protein interaction and pathway enrichment analysis. Our study provides new insights into the association between genotypes and phenotypes of OI and novel information for dissecting the underlying mechanism of the disease.

The rare autosomal dominant disorder acute intermittent porphyria (AIP) is caused by the deficient activity of hydroxymethylbilane synthase (HMBS). The symptoms of AIP are acute neurovisceral attacks which are induced by the dysfunction of heme biosynthesis. To better interpret the underlying mechanism of clinical phenotypes, we collected 117 HMBS gene mutations from reported individuals with AIP and evaluated the mutations’ impacts on the corresponding protein structure and function. We found that several mutations with most severe clinical symptoms are located at dipyromethane cofactor (DPM) binding domain of HMBS. Mutations on these residues likely significantly influence the catalytic reaction. To infer new pathogenic mutations, we evaluated the pathogenicity for all the possible missense mutations of HMBS gene with different bioinformatic prediction algorithms, and identified 34 mutations with serious pathogenicity and low allele frequency. In addition, we found that gene PPARA may also play an important role in the mechanisms of AIP attacks. Our analysis about the distribution frequencies of the 23 variations revealed different distribution patterns among eight ethnic populations, which could help to explain the genetic basis that may contribute to population disparities in AIP prevalence. Our systematic analysis provides a better understanding for this disease and helps for the diagnosis and treatment of AIP.

Combinatorial drug therapy can improve the therapeutic effect and reduce the corresponding adverse events. In silico strategies to classify synergistic vs. antagonistic drug pairs is more efficient than experimental strategies. However, most of the developed methods have been applied only to cancer therapies. In this study, we introduce a novel method, XGBoost, based on five features of drugs and biomolecular networks of their targets, to classify synergistic vs. antagonistic drug combinations from different drug categories. We found that XGBoost outperformed other classifiers in both stratified fivefold cross-validation (CV) and independent validation. For example, XGBoost achieved higher predictive accuracy than other models (0.86, 0.78, 0.78, and 0.83 for XGBoost, logistic regression, naïve Bayesian, and random forest, respectively) for an independent validation set. We also found that the five-feature XGBoost model is much more effective at predicting combinatorial therapies that have synergistic effects than those with antagonistic effects. The five-feature XGBoost model was also validated on TCGA data with accuracy of 0.79 among the 61 tested drug pairs, which is comparable to that of DeepSynergy. Among the 14 main anatomical/pharmacological groups classified according to WHO Anatomic Therapeutic Class, for drugs belonging to five groups, their prediction accuracy was significantly increased (odds ratio < 1) or reduced (odds ratio > 1) (Fisher’s exact test, p < 0.05). This study concludes that our five-feature XGBoost model has significant benefits for classifying synergistic vs. antagonistic drug combinations.

Single-cell RNA sequencing (scRNA-seq) technologies allow the dissection of gene expression at single-cell resolution, which greatly revolutionizes transcriptomic studies. A number of scRNA-seq protocols have been developed, and these methods possess their unique features with distinct advantages and disadvantages. Due to technical limitations and biological factors, scRNA-seq data are noisier and more complex than bulk RNA-seq data. The high variability of scRNA-seq data raises computational challenges in data analysis. Although an increasing number of bioinformatics methods are proposed for analyzing and interpreting scRNA-seq data, novel algorithms are required to ensure the accuracy and reproducibility of results. In this review, we provide an overview of currently available single-cell isolation protocols and scRNA-seq technologies, and discuss the methods for diverse scRNA-seq data analyses including quality control, read mapping, gene expression quantification, batch effect correction, normalization, imputation, dimensionality reduction, feature selection, cell clustering, trajectory inference, differential expression calling, alternative splicing, allelic expression, and gene regulatory network reconstruction. Further, we outline the prospective development and applications of scRNA-seq technologies.

REGγ is a component of ubiquitin-independent 20S proteasome that targets many regulatory proteins for degradation. Here the authors show that REGγ is induced in DSS colitis and promotes degradation of IκBɛ, and that REGγ-deficient mice have less NFκB activation and are more resistant to the disease.

Hepatoblastoma (HB), a leading primary hepatic malignancy in children, originates from primitive hepatic stem cells. This study aimed to uncover the genetic variants that are responsible for HB oncogenesis. One family, which includes the healthy parents, and two brothers affected by HB, was recruited. Whole-genome sequencing (WGS) of germline DNA from all the family members identified two maternal variants, located within APC gene and X-linked WAS gene, which were harbored by the two brothers. The mutation of APC (rs137854573, c.C1606T, p.R536X) could result in HB carcinogenesis by activating Wnt signaling. The WAS variant (c.G3T, p.M1-P5del) could promote HB cell proliferation and inhibit T-cell-based immunity by activating PLK1 signaling and inactivating TCR signaling. Further analysis reflected that WAS deficiency might affect the antitumor activity of natural killer and dendritic cells. In summary, the obtained results imply that an APC mutant together with an X-linked WAS mutant, could lead to HB tumorigenesis by activating Wnt and PLK1 signaling, inhibiting TCR signaling, and reducing the antitumor activity of natural killer and dendritic cells.

DNA sequencing has allowed for the discovery of the genetic cause for a considerable number of diseases, paving the way for new disease diagnostics. However, due to the lack of clinical samples and records, the molecular cause for rare diseases is always hard to identify, significantly limiting the number of rare Mendelian diseases diagnosed through sequencing technologies. Clinical phenotype information therefore becomes a major resource to diagnose rare diseases. In this article, we adopted both a phenotypic similarity method and a machine learning method to build four diagnostic models to support rare disease diagnosis. All the diagnostic models were validated using the real medical records from RAMEDIS. Each model provides a list of the top 10 candidate diseases as the prediction outcome and the results showed that all models had a high diagnostic precision (≥98%) with the highest recall reaching up to 95% while the models with machine learning methods showed the best performance. To promote effective diagnosis for rare disease in clinical application, we developed the phenotype-based Rare Disease Auxiliary Diagnosis system (RDAD) to assist clinicians in diagnosing rare diseases with the above four diagnostic models. The system is freely accessible through http://www.unimd.org/RDAD/.

High-risk neuroblastoma is a very aggressive disease, with excessive tumor growth and poor outcomes. A proper stratification of the high-risk patients by prognostic outcome is important for treatment. However, there is still a lack of survival stratification for the high-risk neuroblastoma. To fill the gap, we adopt a deep learning algorithm, Autoencoder, to integrate multi-omics data, and combine it with K-means clustering to identify two subtypes with significant survival differences. By comparing the Autoencoder with PCA, iCluster, and DGscore about the classification based on multi-omics data integration, Autoencoder-based classification outperforms the alternative approaches. Furthermore, we also validated the classification in two independent datasets by training machine-learning classification models, and confirmed its robustness. Functional analysis revealed that MYCN amplification was more frequently occurred in the ultra-high-risk subtype, in accordance with the overexpression of MYC/MYCN targets in this subtype. In summary, prognostic subtypes identified by deep learning-based multi-omics integration could not only improve our understanding of molecular mechanism, but also help the clinicians make decisions.

Leaf growth and senescence are controlled by tight genetic factors involved regulation at multiple levels. Circular RNAs (circRNAs) have recently been reported as the microRNA sponge to accomplish corresponding regulatory roles. This study aims to explore the expression profile and functional role of circRNAs in Arabidopsis leaf growth and senescence. We used publically available RNA-seq data of Arabidopsis leaves to identify the circular RNA expression profile and used quantitative real-time PCR to validate our identified circRNAs. The functions of circRNAs were explored using distinct bioinformatics methods including analysis of network, gene ontology and KEGG pathway. We identified 168 circRNAs, including 40 novel circRNAs, in Arabidopsis thaliana leaves, with 158 (94.1%) circRNAs arising from the exons of genes. Real-time PCRs were used to verify 4 highly expressed circRNAs and they all showed consistent expression patterns with the RNA-seq results. Interestingly, 6 and 35 circRNAs were differentially expressed at G- to -M stage and M- to -S stage, respectively. The circRNAs display an upregulation trend during the lifespan of Arabidopsis leaves. Moreover, the expression of circRNAs during senescence is independent of host gene expression to a certain degree. The gene ontology (GO) and KEGG pathway analysis of the targeted mRNA of circRNA–miRNA–mRNA network showed that the circRNAs may be involved in plant hormone signal transduction, Porphyrin and chlorophyll metabolism during leaves senescence. Our comprehensive analysis of the expression profile of circRNAs and their potential functions during leaf growth and senescence suggest that circRNAs may function as new post-transcriptional regulators in the senescence of Arabidopsis leaves.

The original version of this Article contained an error in the email address of the corresponding author Jun Qin. The correct email is jqin1965@126.com. The error has been corrected in the HTML and PDF versions of the Article.

Si-Wu-Tang (SWT) is a Traditional Chinese Medicine (TCM) formula widely used for the treatments of gynecological diseases. To explore the pharmacological mechanism of SWT, we incorporated microarray data of SWT with our herbal target database TCMID to analyze the potential activity mechanism of SWT's herbal ingredients and targets. We detected 2,405 differentially expressed genes in the microarray data, 20 of 102 proteins targeted by SWT were encoded by these DEGs and can be targeted by 2 FDA-approved drugs and 39 experimental drugs. The results of pathway enrichment analysis of the 20 predicted targets were consistent with that of 2,405 differentially expressed genes, elaborating the potential pharmacological mechanisms of SWT. Further study from a perspective of protein-protein interaction (PPI) network showed that the predicted targets of SWT function cooperatively to perform their multi-target effects. We also constructed a network to combine herbs, ingredients, targets and drugs together which bridges the gap between SWT and conventional medicine, and used it to infer the potential mechanisms of herbal ingredients. Moreover, based on the hypothesis that the same or similar effects between different TCM formulae may result from targeting the same proteins, we analyzed 27 other TCM formulae which can also treat the gynecological diseases, the subsequent result provides additional insight to understand the potential mechanisms of SWT in treating amenorrhea. Our bioinformatics approach to detect the pharmacology of SWT may shed light on drug discovery for gynecological diseases and could be utilized to investigate other TCM formulae as well.

Human protein complexes play crucial roles in various biological processes as the functional module. However, the expression features of human protein complexes at the transcriptome cascade are poorly understood. Here, we used the RNA-Seq data from 16 disparate tissues and four types of human cancers to explore the characteristics and dynamics of human protein complexes. We observed that many individual components of human protein complexes can be generated by multiple distinct transcripts. Similar with yeast, the human protein complex constituents are inclined to co-express in diverse tissues. The dominant isoform of the genes involved in protein complexes tend to encode the complex constituents in each tissue. Our results indicate that the protein complex dynamics not only correlate with the presence or absence of complexes, but may also be related to the major isoform switching for complex subunits. Between any two cancers of breast, colon, lung and prostate, we found that only a few of the differentially expressed transcripts associated with complexes were identical, but 5-10 times more protein complexes involved in differentially expressed transcripts were common. Collectively, our study reveals novel properties and dynamics of human protein complexes at the transcriptome cascade in diverse normal tissues and different cancers.

Recent progress in stem cell biology, notably cell fate conversion, calls for novel theoretical understanding for cell differentiation. The existing qualitative concept of Waddington's "epigenetic landscape" has attracted particular attention because it captures subsequent fate decision points, thus manifesting the hierarchical ("tree-like") nature of cell fate diversification. Here, we generalized a recent work and explored such a developmental landscape for a two-gene fate decision circuit by integrating the underlying probability landscapes with different parameters (corresponding to distinct developmental stages). The change of entropy production rate along the parameter changes indicates which parameter changes can represent a normal developmental process while other parameters' change can not. The transdifferentiation paths over the landscape under certain conditions reveal the possibility of a direct and reversible phenotypic conversion. As the intensity of noise increases, we found that the landscape becomes flatter and the dominant paths more straight, implying the importance of biological noise processing mechanism in development and reprogramming. We further extended the landscape of the one-step fate decision to that for two-step decisions in central nervous system (CNS) differentiation. A minimal network and dynamic model for CNS differentiation was firstly constructed where two three-gene motifs are coupled. We then implemented the SDEs (Stochastic Differentiation Equations) simulation for the validity of the network and model. By integrating the two landscapes for the two switch gene pairs, we constructed the two-step development landscape for CNS differentiation. Our work provides new insights into cellular differentiation and important clues for better reprogramming strategies.

Subsyndromal symptomatic depression (SSD) is a subtype of subthreshold depressive and also lead to significant psychosocial functional impairment as same as major depressive disorder (MDD). Several studies have suggested that SSD is a transitory phenomena in the depression spectrum and is thus considered a subtype of depression. However, the pathophysioloy of depression remain largely obscure and studies on SSD are limited. The present study compared the expression profile and made the classification with the leukocytes by using whole-genome cRNA microarrays among drug-free first-episode subjects with SSD, MDD, and matched controls (8 subjects in each group). Support vector machines (SVMs) were utilized for training and testing on candidate signature expression profiles from signature selection step. Firstly, we identified 63 differentially expressed SSD signatures in contrast to control (P< = 5.0E-4) and 30 differentially expressed MDD signatures in contrast to control, respectively. Then, 123 gene signatures were identified with significantly differential expression level between SSD and MDD. Secondly, in order to conduct priority selection for biomarkers for SSD and MDD together, we selected top gene signatures from each group of pair-wise comparison results, and merged the signatures together to generate better profiles used for clearly classify SSD and MDD sets in the same time. In details, we tried different combination of signatures from the three pair-wise compartmental results and finally determined 48 gene expression signatures with 100% accuracy. Our finding suggested that SSD and MDD did not exhibit the same expressed genome signature with peripheral blood leukocyte, and blood cell-derived RNA of these 48 gene models may have significant value for performing diagnostic functions and classifying SSD, MDD, and healthy controls.

In their expression process, different genes can generate diverse functional products, including various protein-coding or noncoding RNAs. Here, we investigated the protein-coding capacities and the expression levels of their isoforms for human known genes, the conservation and disease association of long noncoding RNAs (ncRNAs) with two transcriptome sequencing datasets from human brain tissues and 10 mixed cell lines. Comparative analysis revealed that about two-thirds of the genes expressed between brain and cell lines are the same, but less than one-third of their isoforms are identical. Besides those genes specially expressed in brain and cell lines, about 66% of genes expressed in common encoded different isoforms. Moreover, most genes dominantly expressed one isoform and some genes only generated protein-coding (or noncoding) RNAs in one sample but not in another. We found 282 human genes could encode both protein-coding and noncoding RNAs through alternative splicing in the two samples. We also identified more than 1,000 long ncRNAs, and most of those long ncRNAs contain conserved elements across either 46 vertebrates or 33 placental mammals or 10 primates. Further analysis showed that some long ncRNAs differentially expressed in human breast cancer or lung cancer, several of those differentially expressed long ncRNAs were validated by RT-PCR. In addition, those validated differentially expressed long ncRNAs were found significantly correlated with certain breast cancer or lung cancer related genes, indicating the important biological relevance between long ncRNAs and human cancers. Our findings reveal that the differences of gene expression profile between samples mainly result from the expressed gene isoforms, and highlight the importance of studying genes at the isoform level for completely illustrating the intricate transcriptome.

Hepatocellular carcinoma (HCC) is a primary malignancy of the liver, which is closely related to hepatitis C and cirrhosis. The molecular mechanisms underlying the hepatocarcinogenesis induced by HCV infection remain clarified from a standpoint of systems biology. By integrating data from protein-protein interactions, transcriptional regulation, and disease related microarray analysis, we carried out a dynamic biological network analysis on the progression of HCV induced hepatocarcinogenesis, and systematically explored the potentially disease-related mechanisms through a network view. The dysfunctional interactions among proteins and deregulatory relationships between transcription factors and their target genes could be causes for the occurrence and progression of this disease. The six pathologically defined disease stages in the development and progression of HCC after HCV infection were included in this study. We constructed disease-related biological networks for each disease stage, and identified progression-related sub-networks that potentially play roles in the developmental stage of the corresponding disease and participate in the later stage of cancer progression. In addition, we identified novel risk factors related to HCC based on the analysis of the progression-related sub-networks. The dynamic characteristics of the network reflect important features of the disease development and progression, which provide important information for us to further explore underlying mechanisms of the disease.

Mitochondria are major players on the production of energy, and host several key reactions involved in basic metabolism and biosynthesis of essential molecules. Currently, the majority of nucleus-encoded mitochondrial proteins are unknown even for model plant Arabidopsis. We reported a computational framework for predicting Arabidopsis mitochondrial proteins based on a probabilistic model, called Naive Bayesian Network, which integrates disparate genomic data generated from eight bioinformatics tools, multiple orthologous mappings, protein domain properties and co-expression patterns using 1,027 microarray profiles. Through this approach, we predicted 2,311 candidate mitochondrial proteins with 84.67% accuracy and 2.53% FPR performances. Together with those experimental confirmed proteins, 2,585 mitochondria proteins (named CoreMitoP) were identified, we explored those proteins with unknown functions based on protein-protein interaction network (PIN) and annotated novel functions for 26.65% CoreMitoP proteins. Moreover, we found newly predicted mitochondrial proteins embedded in particular subnetworks of the PIN, mainly functioning in response to diverse environmental stresses, like salt, draught, cold, and wound etc. Candidate mitochondrial proteins involved in those physiological acitivites provide useful targets for further investigation. Assigned functions also provide comprehensive information for Arabidopsis mitochondrial proteome.