Discovering monotonic stemness marker genes from time-series stem cell microarray data

Pangenome Data Analysis Reveals Characteristics of Resistance Gene Analogs Associated with Sclerotinia sclerotiorum Resistance in Sunflower

The sunflower, an important oilseed crop and food source across the world, is susceptible to several pathogens, which cause severe losses in sunflower production. The utilization of genetic resistance is the most economical, effective measure to prevent infectious diseases. Based on the sunflower pangenome, in this study, we explored the variability of resistance gene analogs (RGAs) within the species. According to a comparative analysis of RGA candidates in the sunflower pangenome using the RGAugury pipeline, a total of 1344 RGAs were identified, comprising 1107 conserved, 199 varied, and 38 rare RGAs. We also identified RGAs associated with resistance against Sclerotinia sclerotiorum (S. sclerotiorum) in sunflower at the quantitative trait locus (QTL). A total of 61 RGAs were found to be located at four quantitative trait loci (QTLs). Through a detailed expression analysis of RGAs in one susceptible and two tolerant sunflower inbred lines (ILs) across various time points post inoculation, we discovered that 348 RGAs exhibited differential expression in response to Sclerotinia head rot (SHR), with 17 of these differentially expressed RGAs being situated within the QTL regions. In addition, 15 RGA candidates had gene introgression. Our data provide a better understanding of RGAs, which facilitate genomics-based improvements in disease resistance in sunflower.

Generating transcriptional regulatory networks from time-ordered stem cell differentiation RNA sequencing data

We describe steps to 1) identify ascending and descending monotonic key genes from time-ordered stem cell differentiation expression data, 2) construct time-ordered transcriptional regulatory networks, and 3) infer the involvement of transcription factors along the differentiation process.

For complete details on the use and execution of this protocol, please refer to Wong et al. (2020).

•

A protocol to identify ascending and descending monotonic pattern genes

•

Temporal gene regulation in development

•

Pipeline can be adapted to other time-series data sets

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Identification of the Complex Interplay Between Nematode-Related lncRNAs and Their Target Genes in Glycine max L

Soybean (Glycine max) is a major plant protein source and oilseed crop. However, plant-parasitic nematodes (PPNs) affect its annual yield. In the current study, in order to better understand the regulation of defense mechanism against PPNs in soybean, we investigated the role of long non-coding RNAs (lncRNAs) in response to two nematode species, Heterodera glycines (SCN: soybean cyst nematode) and Rotylenchulus reniformis (reniform). To this end, two publicly available RNA-seq data sets (SCN data set and RAD: reniform-associated data set) were employed to discover the lncRNAome profile of soybean under SCN and reniform infection, respectively. Upon identification of unannotated transcripts in these data sets, a seven-step pipeline was utilized to sieve these transcripts, which ended up in 384 and 283 potential lncRNAs in SCN data set and RAD, respectively. These transcripts were then used to predict cis and trans nematode-related targets in soybean genome. Computational prediction of target genes function, some of which were also among differentially expressed genes, revealed the involvement of putative nematode-responsive genes as well as enrichment of multiple stress responses in both data sets. Finally, 15 and six lncRNAs were proposed to be involved in microRNA-mediated regulation of gene expression in soybean in response to SNC and reniform infection, respectively. Collectively, this study provides a novel insight into the signaling and regulatory network of soybean-pathogen interactions and opens a new window for further research.

Unravelling the Complex Interplay of Transcription Factors Orchestrating Seed Oil Content in Brassica napus L

Transcription factors (TFs) and their complex interplay are essential for directing specific genetic programs, such as responses to environmental stresses, tissue development, or cell differentiation by regulating gene expression. Knowledge regarding TF-TF cooperations could be promising in gaining insight into the developmental switches between the cultivars of Brassica napus L., namely Zhongshuang11 (ZS11), a double-low accession with high-oil- content, and Zhongyou821 (ZY821), a double-high accession with low-oil-content. In this regard, we analysed a time series RNA-seq data set of seed tissue from both of the cultivars by mainly focusing on the monotonically expressed genes (MEGs). The consideration of the MEGs enables the capturing of multi-stage progression processes that are orchestrated by the cooperative TFs and, thus, facilitates the understanding of the molecular mechanisms determining seed oil content. Our findings show that TF families, such as NAC, MYB, DOF, GATA, and HD-ZIP are highly involved in the seed developmental process. Particularly, their preferential partner choices as well as changes in their gene expression profiles seem to be strongly associated with the differentiation of the oil content between the two cultivars. These findings are essential in enhancing our understanding of the genetic programs in both cultivars and developing novel hypotheses for further experimental studies.

Computational Identification of Master Regulators Influencing Trypanotolerance in Cattle

African Animal Trypanosomiasis (AAT) is transmitted by the tsetse fly which carries pathogenic trypanosomes in its saliva, thus causing debilitating infection to livestock health. As the disease advances, a multistage progression process is observed based on the progressive clinical signs displayed in the host’s body. Investigation of genes expressed with regular monotonic patterns (known as Monotonically Expressed Genes (MEGs)) and of their master regulators can provide important clue for the understanding of the molecular mechanisms underlying the AAT disease. For this purpose, we analysed MEGs for three tissues (liver, spleen and lymph node) of two cattle breeds, namely trypanosusceptible Boran and trypanotolerant N’Dama. Our analysis revealed cattle breed-specific master regulators which are highly related to distinguish the genetic programs in both cattle breeds. Especially the master regulators MYC and DBP found in this study, seem to influence the immune responses strongly, thereby susceptibility and trypanotolerance of Boran and N’Dama respectively. Furthermore, our pathway analysis also bolsters the crucial roles of these master regulators. Taken together, our findings provide novel insights into breed-specific master regulators which orchestrate the regulatory cascades influencing the level of trypanotolerance in cattle breeds and thus could be promising drug targets for future therapeutic interventions.

Identification of Monotonically Differentially Expressed Genes across Pathologic Stages for Cancers

Given the fact that cancer is a multistage progression process resulting from genetic sequence mutations, the genes whose expression values increase or decrease monotonically across pathologic stages are potentially involved in tumor progression. This may provide insightful clues about how human cancers advance, thereby facilitating more personalized treatments. By replacing the expression values of genes with their GeneRanks, we propose a procedure capable of identifying monotonically differentially expressed genes (MEGs) as the disease advances. Using three real-world gene expression data that cover three distinct cancer types-colon, esophageal, and lung cancers-the proposed procedure has demonstrated excellent performance in detecting the potential MEGs. To conclude, the proposed procedure can detect MEGs across pathologic stages of cancers very efficiently and is thus highly recommended.

Decoding the differentiation of mesenchymal stem cells into mesangial cells at the transcriptomic level

Background

Mesangial cells play an important role in the glomerulus to provide mechanical support and maintaine efficient ultrafiltration of renal plasma. Loss of mesangial cells due to pathologic conditions may lead to impaired renal function. Mesenchymal stem cells (MSC) can differentiate into many cell types, including mesangial cells. However transcriptomic profiling during MSC differentiation into mesangial cells had not been studied yet. The aim of this study is to examine the pattern of transcriptomic changes during MSC differentiation into mesangial cells, to understand the involvement of transcription factor (TF) along the differentiation process, and finally to elucidate the relationship among TF-TF and TF-key gene or biomarkers during the differentiation of MSC into mesangial cells.

Results

Several ascending and descending monotonic key genes were identified by Monotonic Feature Selector. The identified descending monotonic key genes are related to stemness or regulation of cell cycle while ascending monotonic key genes are associated with the functions of mesangial cells. The TFs were arranged in a co-expression network in order of time by Time-Ordered Gene Co-expression Network (TO-GCN) analysis. TO-GCN analysis can classify the differentiation process into three stages: differentiation preparation, differentiation initiation and maturation. Furthermore, it can also explore TF-TF-key genes regulatory relationships in the muscle contraction process.

Conclusions

A systematic analysis for transcriptomic profiling of MSC differentiation into mesangial cells has been established. Key genes or biomarkers, TFs and pathways involved in differentiation of MSC-mesangial cells have been identified and the related biological implications have been discussed. Finally, we further elucidated for the first time the three main stages of mesangial cell differentiation, and the regulatory relationships between TF-TF-key genes involved in the muscle contraction process. Through this study, we have increased fundamental understanding of the gene transcripts during the differentiation of MSC into mesangial cells.

Identifying Cattle Breed-Specific Partner Choice of Transcription Factors during the African Trypanosomiasis Disease Progression Using Bioinformatics Analysis

African Animal Trypanosomiasis (AAT) is a disease caused by pathogenic trypanosomes which affects millions of livestock every year causing huge economic losses in agricultural production especially in sub-Saharan Africa. The disease is spread by the tsetse fly which carries the parasite in its saliva. During the disease progression, the cattle are prominently subjected to anaemia, weight loss, intermittent fever, chills, neuronal degeneration, congestive heart failure, and finally death. According to their different genetic programs governing the level of tolerance to AAT, cattle breeds are classified as either resistant or susceptible. In this study, we focus on the cattle breeds N’Dama and Boran which are known to be resistant and susceptible to trypanosomiasis, respectively. Despite the rich literature on both breeds, the gene regulatory mechanisms of the underlying biological processes for their resistance and susceptibility have not been extensively studied. To address the limited knowledge about the tissue-specific transcription factor (TF) cooperations associated with trypanosomiasis, we investigated gene expression data from these cattle breeds computationally. Consequently, we identified significant cooperative TF pairs (especially DBP−PPARA and DBP−THAP1 in N’Dama and DBP−PAX8 in Boran liver tissue) which could help understand the underlying AAT tolerance/susceptibility mechanism in both cattle breeds.

Identification of monotonically differentially expressed genes for non-small cell lung cancer

Background

Monotonically expressed genes (MEGs) are genes whose expression values increase or decrease monotonically as a disease advances or time proceeds. Non-small cell lung cancer (NSCLC) is a multistage progression process resulting from genetic sequences mutations, the identification of MEGs for NSCLC is important.

Results

With the aid of a feature selection algorithm capable of identifying MEGs – the MFSelector method – two sets of potential MEGs were selected in this study: the MEGs across the different pathologic stages and the MEGs across the risk levels of death for the NSCLC patients at early stages. For the lung adenocarcinoma (AC) subtypes no statistically significant MEGs were identified across pathologic stages, however dozens of MEGs were identified across the risk levels of death. By contrast, for the squamous cell lung carcinoma (SCC) there were no statistically significant MEGs as either stage or risk level advanced.

Conclusions

The pathologic stage of non-small cell lung cancer patients at early stages has no prognostic value, making the identification of prognostic gene signatures for them more meaningful and highly desirable.

Incorporating Pathway Information into Feature Selection towards Better Performed Gene Signatures

To analyze gene expression data with sophisticated grouping structures and to extract hidden patterns from such data, feature selection is of critical importance. It is well known that genes do not function in isolation but rather work together within various metabolic, regulatory, and signaling pathways. If the biological knowledge contained within these pathways is taken into account, the resulting method is a pathway-based algorithm. Studies have demonstrated that a pathway-based method usually outperforms its gene-based counterpart in which no biological knowledge is considered. In this article, a pathway-based feature selection is firstly divided into three major categories, namely, pathway-level selection, bilevel selection, and pathway-guided gene selection. With bilevel selection methods being regarded as a special case of pathway-guided gene selection process, we discuss pathway-guided gene selection methods in detail and the importance of penalization in such methods. Last, we point out the potential utilizations of pathway-guided gene selection in one active research avenue, namely, to analyze longitudinal gene expression data. We believe this article provides valuable insights for computational biologists and biostatisticians so that they can make biology more computable.

Elastic net-based prediction of IFN-β treatment response of patients with multiple sclerosis using time series microarray gene expression profiles

INF-β has been widely used to treat patients with multiple sclerosis (MS) in relapse. Accurate prediction of treatment response is important for effective personalization of treatment. Microarray data have been frequently used to discover new genes and to predict treatment responses. However, conventional analytical methods suffer from three difficulties: high-dimensionality of datasets; high degree of multi-collinearity; and achieving gene identification in time-course data. The use of Elastic net, a sparse modelling method, would decrease the first two issues; however, Elastic net is currently unable to solve these three issues simultaneously. Here, we improved Elastic net to accommodate time-course data analyses. Numerical experiments were conducted using two time-course microarray datasets derived from peripheral blood mononuclear cells collected from patients with MS. The proposed methods successfully identified genes showing a high predictive ability for INF-β treatment response. Bootstrap sampling resulted in an 81% and 78% accuracy for each dataset, which was significantly higher than the 71% and 73% accuracy obtained using conventional methods. Our methods selected genes showing consistent differentiation throughout all time-courses. These genes are expected to provide new predictive biomarkers that can influence INF-β treatment for MS patients.

LncRNAs are altered in lung squamous cell carcinoma and lung adenocarcinoma

Long non-coding RNAs (lncRNAs) have been implicated in pathogenesis of various cancers, including lung squamous cell carcinoma (LUSC) and lung adenocarcinoma (LUAD). We used cBioPortal to analyze lncRNA alteration frequencies and their ability to predict overall survival (OS) using 504 LUSC and 522 LUAD samples from The Cancer Genome Atlas (TCGA) database. In LUSC, 624 lncRNAs had alteration rates > 1% and 64 > 10%. In LUAD 625 lncRNAs had alteration rates > 1% and 36 > 10%. Among those, 620 lncRNAs had alteration frequencies > 1% in both LUSC and LUAD, while 22 were LUSC-specific and 23 were LUAD-specific. Twenty lncRNAs had alteration frequencies > 10% in both LUSC and LUAD, while 44 were LUSC-specific and 16 were LUAD specific. Genome ontology and pathway analyses produced similar results for LUSC and LUAD. Two lncRNAs (IGF2BP2-AS1 and DGCR5) correlated with better OS in LUSC, and three (MIR31HG, CDKN2A-AS1 and LINC01600) predicted poor OS in LUAD. Chip-seq and luciferase reporter assays identified potential IGF2BP2-AS1, DGCR5 and LINC01600 promoters and enhancers. This study presented lncRNA landscapes and revealed differentially expressed, highly altered lncRNAs in LUSC and LUAD. LncRNAs that act as oncogenes and lncRNA-regulating transcription factors provide novel targets for anti-lung cancer therapeutics.

Classification of early-stage non-small cell lung cancer by weighing gene expression profiles with connectivity information

Pathway-based feature selection algorithms, which utilize biological information contained in pathways to guide which features/genes should be selected, have evolved quickly and become widespread in the field of bioinformatics. Based on how the pathway information is incorporated, we classify pathway-based feature selection algorithms into three major categories-penalty, stepwise forward, and weighting. Compared to the first two categories, the weighting methods have been underutilized even though they are usually the simplest ones. In this article, we constructed three different genes' connectivity information-based weights for each gene and then conducted feature selection upon the resulting weighted gene expression profiles. Using both simulations and a real-world application, we have demonstrated that when the data-driven connectivity information constructed from the data of specific disease under study is considered, the resulting weighted gene expression profiles slightly outperform the original expression profiles. In summary, a big challenge faced by the weighting method is how to estimate pathway knowledge-based weights more accurately and precisely. Only until the issue is conquered successfully will wide utilization of the weighting methods be impossible.

Systematic identification of long non-coding RNAs with cancer-testis expression patterns in 14 cancer types

Cancer-testis (CT) genes are a group of genes that are potential targets of immunotherapy and candidate epi-drivers participating in the development of cancers. Previous studies mainly focused on protein-coding genes, neglecting long non-coding RNAs with the same expression patterns. In this study, we performed a systematic investigation of cancer-testis long non-coding RNAs (CT-lncRNAs) with multiple independent open-access databases.We identified 1,325 extremely highly expressed CT-lncRNAs (EECT-lncRNAs) in 14 cancer types. Functional annotation revealed that CT-lncRNAs reactivated in cancers could promote genome instability and the malignant potential of cancers. We observed a mutually exclusive pattern of EECT-lncRNA activation and mutation in known oncogenes, suggesting their potential role as drivers of cancer that complement known mut-driver genes. Additionally, we provided evidence that testis-specific regulatory elements and promoter hypo-methylation may be EECT-lncRNA activation mechanisms, and EECT-lncRNAs may regulate CT gene reactivation. Taken together, our study puts forth a new hypothesis in the research field of CT genes, whereby CT-lncRNAs/EECT-lncRNAs play important roles in the progression and maintenance of tumorigenesis, expanding candidate CT epi-driver genes from coding genes to non-coding RNAs.

Identification of potential prognostic ceRNA module biomarkers in patients with pancreatic adenocarcinoma

Accumulating evidence suggested that long non-coding RNAs (lncRNAs) can function as competing endogenous RNAs (ceRNAs) to interact with other RNA transcripts and ceRNAs perturbation play important roles in cancer initiation and progression including pancreatic adenocarcinoma (PAAD). In this study, we constructed a PAAD-specific hallmark gene-related ceRNA network (HceNet) using paired genome-wide expression profiles of mRNA, lncRNA and miRNA and regulatory relationships between them. Based on “ceRNA hypothesis”, we analyzed the characteristics of HceNet and identified a ceRNA module comprising of 29 genes (12 lncRNAs, two miRNAs and 15 mRNAs) as potential prognostic biomarkers related to overall survival of patients with PAAD. The prognostic value of ceRNA module biomarkers was further validated in the train (Hazard Ratio (HR) =1.661, 95% CI: 1.275–2.165, p<1.00e-4), test (HR=1.546, 95% CI: 1.238-1.930, p<1.00e-4), and entire (HR=1.559, 95% CI: 1.321-1.839, p<1.00e-4) datasets. Our study provides candidate prognostic biomarkers for PAAD and increases our understanding of ceRNA-related regulatory mechanism in PAAD pathogenesis.

LncSubpathway: a novel approach for identifying dysfunctional subpathways associated with risk lncRNAs by integrating lncRNA and mRNA expression profiles and pathway topologies

Long non-coding RNAs (lncRNAs) play important roles in various biological processes, including the development of many diseases. Pathway analysis is a valuable aid for understanding the cellular functions of these transcripts. We have developed and characterized LncSubpathway, a novel method that integrates lncRNA and protein coding gene (PCG) expression with interactome data to identify disease risk subpathways that functionally associated with risk lncRNAs. LncSubpathway identifies the most relevance regions which are related with risk lncRNA set and implicated with study conditions through simultaneously considering the dysregulation extent of lncRNAs, PCGs and their correlations. Simulation studies demonstrated that the sensitivity and false positive rates of LncSubpathway were within acceptable ranges, and that LncSubpathway could accurately identify dysregulated regions that related with disease risk lncRNAs within pathways. When LncSubpathway was applied to colorectal carcinoma and breast cancer subtype datasets, it identified cancer type- and breast cancer subtype-related meaningful subpathways. Further, analysis of its robustness and reproducibility indicated that LncSubpathway was a reliable means of identifying subpathways that functionally associated with lncRNAs. LncSubpathway is freely available at http://www.bio-bigdata.com/lncSubpathway/.

Systemically identifying and prioritizing risk lncRNAs through integration of pan-cancer phenotype associations

LncRNAs have emerged as a major class of regulatory molecules involved in normal cellular physiology and disease, our knowledge of lncRNAs is very limited and it has become a major research challenge in discovering novel disease-related lncRNAs in cancers. Based on the assumption that diverse diseases with similar phenotype associations show similar molecular mechanisms, we presented a pan-cancer network-based prioritization approach to systematically identify disease-specific risk lncRNAs by integrating disease phenotype associations. We applied this strategy to approximately 2800 tumor samples from 14 cancer types for prioritizing disease risk lncRNAs. Our approach yielded an average area under the ROC curve (AUC) of 80.66%, with the highest AUC (98.14%) for medulloblastoma. When evaluated using leave-one-out cross-validation (LOOCV) for prioritization of disease candidate genes, the average AUC score of 97.16% was achieved. Moreover, we demonstrated the robustness as well as the integrative importance of this approach, including disease phenotype associations, known disease genes and the numbers of cancer types. Taking glioblastoma multiforme as a case study, we identified a candidate lncRNA gene SNHG1 as a novel disease risk factor for disease diagnosis and prognosis. In summary, we provided a novel lncRNA prioritization approach by integrating pan-cancer phenotype associations that could help researchers better understand the important roles of lncRNAs in human cancers.

Subpathway-LNCE: Identify dysfunctional subpathways competitively regulated by lncRNAs through integrating lncRNA-mRNA expression profile and pathway topologies

Recently, studies have reported that long noncoding RNAs (lncRNAs) can act as modulators of mRNAs through competitively binding to microRNAs (miRNAs) and have relevance to tumorigenesis as well as other diseases. Identify lncRNA competitively regulated subpathway not only can gain insight into the initiation and progression of disease, but also help for understanding the functional roles of lncRNAs in the disease context. Here, we present an effective method, Subpathway-LNCE, which was specifically designed to identify lncRNAs competitively regulated functions and the functional roles of these competitive regulation lncRNAs have not be well characterized in diseases. Moreover, the method integrated lncRNA-mRNA expression profile and pathway topologies. Using prostate cancer datasets and LUAD data sets, we confirmed the effectiveness of our method in identifying disease associated dysfunctional subpathway that regulated by lncRNAs. By analyzing kidney renal clear cell carcinoma related lncRNA competitively regulated subpathway network, we show that Subpathway-LNCE can help uncover disease key lncRNAs. Furthermore, we demonstrated that our method is reproducible and robust. Subpathway-LNCE provide a flexible tool to identify lncRNA competitively regulated signal subpathways underlying certain condition, and help to expound the functional roles of lncRNAs in various status. Subpathway-LNCE has been developed as an R package freely available at https://cran.rstudio.com/web/packages/SubpathwayLNCE/.