RNA-Seq SSRs and small RNA-Seq SSRs: New approaches in cancer biomarker discovery

Exploring the Genetic Basis of Calonectria spp. Resistance in Eucalypts

Selecting high-quality varieties with disease resistance by artificial crossbreeding is the most fundamental way to address the damage caused by Calonectria spp. in eucalypt plantations. However, understanding the mechanism of disease-resistant heterosis occurrence in eucalypts is crucial for successful crossbreeding. Two eucalypt hybrids, the susceptible EC333 (H1522 × unknown) and the resistant EC338 (W1767 × P9060), were screened through infection with Calonectria isolates, a pathogen that causes eucalypt leaf blight. RNA-Seq was performed on the susceptible hybrid, the disease-resistant hybrid, and their parents. The gene differential expression analysis showed that there were 3912 differentially expressed genes between EC333 and EC338, with 1631 up-regulated and 2281 down-regulated genes. The expression trends of the differential gene sets in P9060 and EC338 were similar. However, the expression trend of W1767 was opposite that of EC338. The similarity of the expression and the advantage of stress resistance in E. pellita suggested that genes with significant differences in expression likely relate to disease resistance. A GSEA based on GO annotations revealed that the carbohydrate binding pathway genes were differentially expressed between EC338 and EC333. The gene pathways that were differentially expressed between EC338 and EC333 revealed by the GSEA based on KEGG annotations were the sesquiterpenoid and triterpenoid biosynthesis pathways. The alternative splicing analysis demonstrated that an AS event between EC338 and EC333 occurred in LOC104426602. According to our SNP analysis, EC338 had 626 more high-impact mutation loci than the male parent P9060 and 396 more than the female parent W1767; W1767 had 259 more mutation loci in the downstream region than EC338, while P9060 had 3107 fewer mutation loci in the downstream region than EC338. Additionally, EC338 had 9631 more mutation loci in the exon region than EC333. Modules were found via WGCNA that were strongly and oppositely correlated with EC338 and EC333, such as module MEsaddlebrown, likely associated with leaf blight resistance. The present study provides a detailed explanation of the genetic basis of eucalypt leaf blight resistance, providing the foundation for exploring genes related to this phenomenon.

A reference transcriptome for walnut anthracnose pathogen, Ophiognomonia leptostyla, guides the discovery of candidate virulence genes

Despite the economic losses due to the walnut anthracnose, Ophiognomonia leptostyla is an orphan fungus with respect to genomic resources. In the present study, the transcriptome of O. leptostyla was assembled for the first time. RNA sequencing was conducted for the fungal mycelia grown in a liquid media, and the inoculated leaf samples of walnut with the fungal conidia sampled at 48, 96 and 144 h post inoculation (hpi). The completeness, correctness, and contiguity of the de novo transcriptome assemblies generated with Trinity, Oases, SOAPdenovo-Trans and Bridger were compared to identify a single superior reference assembly. In most of the assessment criteria including N50, Transrate score, number of ORFs with known description in gene bank, the percentage of reads mapped back to the transcript (RMBT), BUSCO score, Swiss-Prot coverage bin and RESM-EVAL score, the Bridger assembly was the superior and thus used as a reference for profiling the O. leptostyla transcriptome in liquid media vs. during walnut infection. The k-means clustering of transcripts resulted in four distinct transcription patterns across the three sampling time points. Most of the detected CAZy transcripts had elevated transcription at 96 hpi that is hypothetically concurrent with the start of intracellular growth. The in-silico analysis revealed 103 candidate effectors of which six were members of Necrosis and Ethylene Inducing Like Protein (NLP) gene family belonging to three distinct k-means clusters. This study provided a complex and temporal pattern of the CAZys and candidate effectors transcription during six days post O. leptostyla inoculation on walnut leaves, introducing a list of candidate virulence genes for validation in future studies.

RNA sequencing: new technologies and applications in cancer research

Over the past few decades, RNA sequencing has significantly progressed, becoming a paramount approach for transcriptome profiling. The revolution from bulk RNA sequencing to single-molecular, single-cell and spatial transcriptome approaches has enabled increasingly accurate, individual cell resolution incorporated with spatial information. Cancer, a major malignant and heterogeneous lethal disease, remains an enormous challenge in medical research and clinical treatment. As a vital tool, RNA sequencing has been utilized in many aspects of cancer research and therapy, including biomarker discovery and characterization of cancer heterogeneity and evolution, drug resistance, cancer immune microenvironment and immunotherapy, cancer neoantigens and so on. In this review, the latest studies on RNA sequencing technology and their applications in cancer are summarized, and future challenges and opportunities for RNA sequencing technology in cancer applications are discussed.

Identification of novel ΔNp63α-regulated miRNAs using an optimized small RNA-Seq analysis pipeline

Advances in high-throughput sequencing have enabled profiling of microRNAs (miRNAs), however, a consensus pipeline for sequencing of small RNAs has not been established. We built and optimized an analysis pipeline using Partek Flow, circumventing the need for analyzing data via scripting languages. Our analysis assessed the effect of alignment reference, normalization method, and statistical model choice on biological data. The pipeline was evaluated using sequencing data from HaCaT cells transfected with either a non-silencing control or siRNA against ΔNp63α, a p53 family member protein which is highly expressed in non-melanoma skin cancer and shown to regulate a number of miRNAs. We posit that 1) alignment and quantification to the miRBase reference provides the most robust quantitation of miRNAs, 2) normalizing sample reads via Trimmed Mean of M-values is the most robust method for accurate downstream analyses, and 3) use of the lognormal with shrinkage statistical model effectively identifies differentially expressed miRNAs. Using our pipeline, we identified previously unrecognized regulation of miRs-149-5p, 18a-5p, 19b-1-5p, 20a-5p, 590-5p, 744-5p and 93-5p by ΔNp63α. Regulation of these miRNAs was validated by RT-qPCR, substantiating our small RNA-Seq pipeline. Further analysis of these miRNAs may provide insight into ΔNp63α's role in cancer progression. By defining the optimal alignment reference, normalization method, and statistical model for analysis of miRNA sequencing data, we have established an analysis pipeline that may be carried out in Partek Flow or at the command line. In this manner, our pipeline circumvents some of the major hurdles encountered during small RNA-Seq analysis.

Construction of synergy networks from gene expression data related to disease

A few methods have been developed to determine whether genes collaborate with each other in relation to a particular disease using an information theoretic measure of synergy. Here, we propose an alternative definition of synergy and justify that our definition improves upon the existing measures of synergy in the context of gene interactions. We use this definition on a prostate cancer data set consisting of gene expression levels in both cancerous and non-cancerous samples and identify pairs of genes which are unable to discriminate between cancerous and non-cancerous samples individually but can do so jointly when we take their synergistic property into account. We also propose a very simple yet effective technique for computation of conditional entropy at a very low cost. The worst case complexity of our method is O(n) while the best case complexity of a state-of-the-art method is O(n(2)). Furthermore, our method can also be extended to find synergistic relation among triplets or even among a larger number of genes. Finally, we validate our results by demonstrating that these findings cannot be due to pure chance and provide the relevance of the synergistic pairs in cancer biology.

Comparative analysis of oncogenes identified by microarray and RNA-sequencing as biomarkers for clinical prognosis

Aims: Although RNA-sequencing has been widely used to identify the differentially expressed genes (DEGs) as biomarkers to guide the therapeutic treatment, it is necessary to investigate the concordance of DEGs identified by microarray and RNA-sequencing for the clinical prognosis. Material & methods: By using The Cancer Genome Atlas data sets, we thoroughly investigated the concordance of DEGs identified from microarray and RNA-sequencing data and their molecular functions. Results: The DEGs identified by both technologies averaged ˜98.6% overlap. The cancer-related gene sets were significantly enriched with the DEGs and consistent between two technologies. Conclusions: The highly consistency of DEGs in their regulation directionality and molecular functions indicated the good reproducibility between microarray and RNA-sequencing in identifying potential oncogenes for clinical prognosis.