Comprehensive transcriptome profiling of squamous cell carcinoma of horn in Bos indicus

Assessing the double-edged of extracellular signal-regulated kinase/CCAAT-enhancer-binding protein beta signaling pathway in arsenic-induced skin damage and its potential foodborne interventions

Arsenic exposure is a major environmental public health challenge worldwide. As typical manifestations for arsenic exposure, the pathogenesis of arsenic-induced skin lesions has not been fully elucidated, as well as the lack of effective control measures. In this study, we first determined the short-term and high-dose arsenic exposure can increase the apoptosis rates, while long-term low-dose arsenic exposure decrease the apoptosis rates. Then, the HaCaT cells with knockdown and overexpression of CCAAT-enhancer-binding protein β (CEBPB) and extracellular signal-regulated kinase (ERK) were constructed. The results demonstrate that knockdown of CEBPB and ERK can reduce NaAsO2 -induced cell apoptosis by inhibiting ERK/CEBPB signaling pathway and vice versa. Further cells were treated with Kaji-Ichigoside F1 (KF1). The results clearly show that KF1 can decrease the arsenic-induced cell apoptosis rates and the expression of ERK/CEBPB signaling pathway-related genes. These results provide evidence that ERK/CEBPB signaling pathway acts as a double-edged sword in arsenic-induced skin damage. Another interesting finding was that KF1 can alleviate arsenic-induced skin cell apoptosis by inhibiting the ERK/CEBPB signaling pathway. This study will contribute to a deeper understanding of the mechanisms of arsenic-induced skin cell apoptosis, and our findings will help to identify a potential food-borne intervention in arsenic detoxification.

Bioinformatics analysis of the role of CXC ligands in the microenvironment of head and neck tumor

Chemokines play a significant role in cancer. CXC-motif chemokine ligands (CXCLs) are associated with the tumorigenesis and progression of head and neck squamous cell carcinoma (HNSC); however, their specific functions in the tumor microenvironment remain unclear. Here, we analyzed the molecular networks and transcriptional data of HNSC patients from the Oncomine, GEPIA, String, cBioPortal, Metascape, TISCH, and TIMER databases. To verify immune functions of CXCLs, their expression was analyzed in different immune cell types. To our knowledge, this is the first report on the correlation between CXCL9-12 and 14 expression and advanced tumor stage. CXCL2, 3, 8, 10, 13, and 16 were remarkably related to tumor immunity. Kaplan-Meier and TIMER survival analyses revealed that high expression of CXCL1, 2, 4, and 6-8 is correlated with low survival in HNSC patients, whereas high expression of CXCL9, 10, 13, 14, and 17 predicts high survival. Only CXCL13 and 14 were associated with overall survival in human papilloma virus (HPV)-negative patients. Single-cell datasets confirmed that CXCLs are associated with HNSC-related immune cells. Thus, CXCL1-6, 8-10, 12-14, and 17 could be prognostic targets for HNSC, and CXCL13 and 14 could be novel biomarkers of HPV-negative HNSC.

Molecular portrait of squamous cell carcinoma of the bovine horn evaluated by high-throughput targeted exome sequencing: a preliminary report

Background

Squamous Cell Carcinoma of horn, also known as horn cancer, is a prevailing type of cancer in cattles especially Bos indicus. It is one of the most prevalent disease in Indian bullocks often resulting in death and huge economic losses to farmers. Here, we have reported the use of targeted exome sequencing to identify variants present in horn cancer affected horn mucosa tissue and blood of the same animal to identify some of the prevalent markers of horn cancer.

Results

We have observed higher number of variants present in tissue as compared to blood as well as among cancer samples compared to samples from normal animals. Eighty six and 1437 cancer-specific variants were identified among the predicted variants in blood and tissue samples, respectively. Total 25 missense variants were observed distributed over 18 genes. KRT8 gene coding for Keratin8, one of the key constituents of horn, displayed 5 missense variants. Additionally, three other genes involved in apoptosis pathway and two genes involved in antigen presentation and processing also contained variants.

Conclusions

Several genes involved in various apoptotic pathways were found to contain non-synonymous mutations. Keratin8 coding for Keratin, a chief constituent of horn was observed to have the highest number of mutations. In all, we present a preliminary report of mutations observed in horn cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-020-02683-y.

Expression profiling revealed keratins and interleukins as potential biomarkers in squamous cell carcinoma of horn in Indian bullocks (Bos indicus)

Horn cancer is most prevalent in Bos indicus and poorly defined genetic landscape makes disease diagnosis and treatment difficult. In this study, RNA-Seq and data analysis using CLC Genomics Workbench was employed to identify biomarkers associated with horn cancer. As a result, a total of 149 genes were found significant differentially expressed in horn cancer samples compared to horn normal samples. The study revealed 'keratins' and 'interleukins' as apex groups of significant differentially expressed genes (DEGs). Functional analysis showed that the upregulated keratins support metastasis of tumor via cell proliferation, migration, and affecting cell stability, while downregulated interleukins along with other associated chemokine receptors deprive the immune response to tumor posing clear path for metastasis of horn cancer. Combi-action of both the group facilitates the tumor microenvironment to reproduce tumorigenesis. Analysis of pathways enriched in DEGs and exemplified protein-protein interaction network indicated actual role of DEGs in horn cancer at a fine level. Important effect of deregulated expression of keratin and interleukin genes in horn cancer enrolling their candidacy as potential biomarkers for horn cancer prognosis. This study appraises the possibility to mitigate horn cancer at fine resolution to extract attainable identification of prognostic molecular portraits.

Association Between Unstable Angina and CXCL17: a New Potential Biomarker

Atherosclerosis and chemokines are strongly related, but the role of the chemokine CXCL17 in atherogenesis is still poorly understood. We aim to investigate the serum CXCL17 levels in different stages of patients with coronary heart disease and explore whether these differences contribute to atherosclerosis. In the current prospective study, we enrolled 48 patients with unstable angina (UA), 51 patients with stable angina (SA) and 41 patients for the control group (CG). All subjects were diagnosed by coronary angiography and Gensini score was used to evaluate the severity of coronary artery disease. The CXCL17 levels were determined using ELISA, while lipid metabolism indicators and high sensitivity C-reactive protein (hs-CRP) were detected by automatic biochemical analyzer. We observed that the unstable angina group had higher CXCL17 levels compared with the stable angina and the control group. The logistic regression analysis showed that CXCL17 was an independent risk factor for unstable angina. Our results showed that CXCL17 was also statistically correlated with hs-CRP, while it was irrelevant with Gensini score. CXCL17 levels were associated with activity of inflammatory response and the instability of atherosclerotic plaques. These results suggest that CXCL17 elevation may be a potential new biomarker of unstable angina.

Identification and Validation of an Immune-Related RNA Signature to Predict Survival of Patients With Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a heterogeneous disease characterized by different molecular subgroups and clinical features. Therefore, it is important to uncover reliable molecular biomarkers for distinguishing different risk patient subgroup. Here, we conducted a multi-omics analysis to examine the joint predictive power of a multi-type RNA signature in the prognosis of HNSCC patients through integration analysis of mRNA, miRNA, and lncRNA expression profiles and clinical data in a large number of HNSCC patients. A multi-type RNA signature (15SigRS) was constructed which can classify patients into the high-risk group and low-risk group with the significantly different outcome [hazard ratio (HR) = 2.718, 95% confidence interval (CI), 2.258–3.272, p < 0.001] in the discovery data set, and subsequently validated in the Cancer Genome Atlas (TCGA) testing data set (HR = 1.299, 95% CI, 1.170–1.442, p < 0.001) and another independent GSE65858 data set (HR = 1.077, 95% CI, 1.016–1.143, p = 0.013). Further multivariate Cox regression analysis and stratification analysis demonstrated the independence of predictive performance of the 15SigRS relative to conventional clinicopathological factors. Furthermore, the 15SigRS has a prior performance in prognostic prediction than other single RNA type-based signatures. Functional analysis suggested that the 15SigRS are involved in immune- or metabolism-related KEGG pathways. In summary, our study demonstrated the potential application of mixed RNA types as molecular markers for predicting the outcome of cancer patients.

Transcriptome analysis identified aberrant gene expression in pollen developmental pathways leading to CGMS in cotton (Gossypium hirsutum L.)

Male sterility (induced or natural) is a potential tool for commercial hybrid seed production in different crops. Despite numerous endeavors to understand the physiological, hereditary, and molecular cascade of events governing CMS in cotton, the exact biological process controlling sterility and fertility reconstruction remains obscure. During current study, RNA-Seq using Ion Torrent S5 platform is carried out to identify 'molecular portraits' in floral buds among the Cytoplasmic Genic Male Sterility (CGMS) line, its near-isogenic maintainer, and restorer lines. A total of 300, 438 and 455 genes were differentially expressed in CGMS, Maintainer, and Restorer lines respectively. The functional analysis using AgriGo revealed suppression in the pathways involved in biogenesis and metabolism of secondary metabolites which play an important role in pollen and anther maturation. Enrichment analysis showed dearth related to pollen and anther's development in sterile line, including anomalous expression of genes and transcription factors that have a role in the development of the reproductive organ, abnormal cytoskeleton formation, defects in cell wall formation. The current study found aberrant expression of DYT1, AMS and cytochrome P450 genes involved in tapetum formation, pollen development, pollen exine and anther cuticle formation associated to male sterility as well as fertility restoration of CGMS. In the current study, more numbers of DEGs were found on Chromosome D05 and A05 as compared to other chromosomes. Expression pattern analysis of fourteen randomly selected genes using qRT-PCR showed high concurrence with gene expression profile of RNA-Seq analysis accompanied by a strong correlation of 0.82. The present study provides an important support for future studies in identifying interaction between cyto-nuclear molecular portraits, to accelerate functional genomics and molecular breeding related to cytoplasmic male sterility studies in cotton.

Identification of Differentially Expressed Genes between Original Breast Cancer and Xenograft Using Machine Learning Algorithms

Breast cancer is one of the most common malignancies in women. Patient-derived tumor xenograft (PDX) model is a cutting-edge approach for drug research on breast cancer. However, PDX still exhibits differences from original human tumors, thereby challenging the molecular understanding of tumorigenesis. In particular, gene expression changes after tissues are transplanted from human to mouse model. In this study, we propose a novel computational method by incorporating several machine learning algorithms, including Monte Carlo feature selection (MCFS), random forest (RF), and rough set-based rule learning, to identify genes with significant expression differences between PDX and original human tumors. First, 831 breast tumors, including 657 PDX and 174 human tumors, were collected. Based on MCFS and RF, 32 genes were then identified to be informative for the prediction of PDX and human tumors and can be used to construct a prediction model. The prediction model exhibits a Matthews coefficient correlation value of 0.777. Seven interpretable interactions within the informative gene were detected based on the rough set-based rule learning. Furthermore, the seven interpretable interactions can be well supported by previous experimental studies. Our study not only presents a method for identifying informative genes with differential expression but also provides insights into the mechanism through which gene expression changes after being transplanted from human tumor into mouse model. This work would be helpful for research and drug development for breast cancer.

Transcriptomic comparison of primary bovine horn core carcinoma culture and parental tissue at early stage

Aim:

Squamous cell carcinoma or SCC of horn in bovines (bovine horn core carcinoma) frequently observed in Bos indicus affecting almost 1% of cattle population. Freshly isolated primary epithelial cells may be closely related to the malignant epithelial cells of the tumor. Comparison of gene expression in between horn’s SCC tissue and its early passage primary culture using next generation sequencing was the aim of this study.

Materials and Methods:

Whole transcriptome sequencing of horn’s SCC tissue and its early passage cells using Ion Torrent PGM were done. Comparative expression and analysis of different genes and pathways related to cancer and biological processes associated with malignancy, proliferating capacity, differentiation, apoptosis, senescence, adhesion, cohesion, migration, invasion, angiogenesis, and metabolic pathways were identified.

Results:

Up-regulated genes in SCC of horn’s early passage cells were involved in transporter activity, catalytic activity, nucleic acid binding transcription factor activity, biogenesis, cellular processes, biological regulation and localization and the down-regulated genes mainly were involved in focal adhesion, extracellular matrix receptor interaction and spliceosome activity.

Conclusion:

The experiment revealed similar transcriptomic nature of horn’s SCC tissue and its early passage cells.

Identification and quantification of novel RNA isoforms in horn cancer of Bos indicus by comprehensive RNA-Seq

Horn cancer (HC) is a squamous cell carcinoma of horn, commonly observed in Bos indicus of the Asian countries. To elucidate the complexity of alternative splicing present in the HC, high-throughput sequencing and analysis of HC and matching horn normal (HN) tissue were carried out. A total of 535,067 and 849,077 reads were analysed after stringent quality filtering for HN and HC, respectively. Cufflinks pipeline for transcriptome analysis revealed 4786 novel splice isoforms comprising 2432 exclusively in HC, 2055 exclusively in HN and 298 in both the conditions. Based on pathway clustering and in silico verification, 102 novel splice isoforms were selected and further analysed with respect to change in protein sequence using Blastp. Finally, fourteen novel splicing events supported both by Cufflinks and UCSC genome browser were selected and confirmed expression by RT-qPCR. Future studies targeted at in-depth characterization of these potential candidate splice isoforms might be helpful in the development of relevant biomarkers for early diagnosis of HC. The results reported in this study refine the available information on transcriptome repertoire of bovine species and boost the research in the line of development of relevant biomarkers for early diagnosis of HC.

Identification of novel SNPs in differentially expressed genes and its association with horn cancer of Bos indicus bullocks by next-generation sequencing

The use of polymorphic markers like SNPs promises to provide comprehensive tool for analysing genome and identifying genomic regions that contribute to cancer phenotype. Horn cancer is the most common cancer among Bos indicus animals. Increased expression of some genes due to polymorphisms increases risk of HC incidence. We successfully amplified 91 SNPs located in 69 genes in 52 samples, each of HC and HN. Equimolar concentration of amplicons from 69 PCR products of each sample was pooled and subjected to sequencing using Ion Torrent PGM. Data obtained were analysed using DNASTAR software package and case control analysis using SAS software. We found SNP present in BPIFA1 gene of B. indicus shows association with event of HC which reflects its potential to be a genetic marker. Bioinformatic analysis to detect structural and functional impact nsSNP of BPIFA1 added another layer of confirmation to our result. We successfully identified SNP associated with HC as well as demonstrated efficient approach for limited number of SNP discovery and validation in targeted genomics regions in large number of samples combining PCR amplification and Ion Torrent PGM sequencing which suits small-scale laboratories with limited budget.

Upregulation of Bpifb1 expression in the parotid glands of non-obese diabetic mice

OBJECTIVE

To define the increased mRNA expression of Bpifb1, a member of the bactericidal/permeability-increasing protein family, in parotid acinar cells from non-obese diabetic (NOD) mice, an animal model for Sjögren's syndrome.

MATERIALS AND METHODS

Parotid acinar cells were prepared from female NOD (NOD/ShiJcl) mice with or without diabetes, as well as from control (C57BL/6JJcl) mice. Total RNA and homogenate were prepared from the parotid acinar cells. Embryonic cDNA from a Mouse MTC(™) Panel I kit was used. The expression of Bpifb1 was determined by cDNA microarray analysis, RT-PCR, real-time PCR, northern blotting and in situ hybridization.

RESULTS

The expression of Bpifb1 mRNA was high in parotid acinar cells from diabetic and non-diabetic NOD mice at 5-50 weeks of age. Acinar cells in the C57BL/6 mice had a low expression of Bpifb1 mRNA at an age >8 weeks, but had a relatively high expression in the foetus and infantile stages.

CONCLUSIONS

Bpifb1 mRNA is upregulated in parotid acinar cells in NOD mice, but its expression is not related to the onset of diabetes. These findings suggest that high expression levels of Bpifb1 might predict disease traits before the onset of autoimmunity.