Assessment of single-nucleotide variant discovery protocols in RNA-seq data from human cells exposed to mycotoxins.

Food and feed contamination by nonlegislated mycotoxins beauvericin (BEA) and enniatin B (ENB) is a worldwide health concern in the present. The principal objective of this work is to assess some of the existing protocols to discover the single nucleotide variants (SNVs) in transcriptomic data obtained by RNA-seq from Jurkat cells in vitro samples individually exposed to BEA and ENB at three concentration levels (1.5, 3 and 5 µM). Moreover, previous transcriptomic results will be compared with new findings obtained using a different protocol. SNVs rs201003509 in BEA exposed cells and the rs36045790 in ENB were found in the differentially expressed genes in all doses compared to controls by means of the Genome Analysis Toolkit (GATK) Best Practices workflow. SNV-RNA-seq complementary pipeline did not show any SNV. Concerning gene expression, discrepant results were found for 1.5 µM BEA exposed cells compared with previous findings. However, 354 overlapped differentially expressed genes (DEGs) were identified in the three ENB concentrations used, with 147 matches with respect to the 245 DEGs found in the previous results. In conclusion, the two discovery SNVs protocols based on variant calling from RNA-seq used in this work displayed very different results and there were SNVs found manually not identified by any pipeline. Additionally, the new gene expression analysis reported comparable but non identical DEGs to the previous transcriptomic results obtained from these RNA-seq data.

EPDR1 up-regulation in human colorectal cancer is related to staging and favours cell proliferation and invasiveness.

The finding of novel molecular markers for prediction or prognosis of invasiveness in colorectal cancer (CRC) constitutes an appealing challenge. Here we show the up-regulation of EPDR1 in a prospective cohort of 101 CRC patients, in a cDNA array of 43 patients and in in silico analyses. EPDR1 encodes a protein related to ependymins, a family of glycoproteins involved in intercellular contacts. A thorough statistical model allowed us to conclude that the gene is significantly up-regulated in tumour tissues when compared with normal mucosa. These results agree with those obtained by the analysis of three publicly available databases. EPDR1 up-regulation correlates with the TNM staging parameters, especially T and M. Studies with CRC cell lines revealed that the methylation of a CpG island controls EPDR1 expression. siRNA knocking-down and overexpression of the gene following transient plasmid transfection, showed that EPDR1 favours cell proliferation, migration, invasiveness and adhesion to type I collagen fibres, suggesting a role in epithelial to mesenchymal transition. Both statistical and functional analysis correlated EPDR1 overexpression with invasiveness and dissemination of tumour cells, supporting the inclusion of EPDR1 in panels of genes used to improve molecular subtyping of CRC. Eventually, EPDR1 may be an actionable target.