In this study, the apoptotic mechanism and combinatorial chemotherapeutic effects of the cytotoxic phenylpropanoid compound 1'S-1'-acetoxyeugenol acetate (AEA), extracted from rhizomes of the Malaysian ethnomedicinal plant Alpinia conchigera Griff. (Zingiberaceae), on MCF-7 human breast cancer cells were investigated for the first time. Data from cytotoxic and apoptotic assays such as live and dead and poly-(ADP-ribose) polymerase cleavage assays indicated that AEA was able to induce apoptosis in MCF-7 cells, but not in normal human mammary epithelial cells. A microarray global gene expression analysis of MCF-7 cells, treated with AEA, suggested that the induction of tumor cell death through apoptosis was modulated through dysregulation of the nuclear factor-kappaB (NF-κB) pathway, as shown by the reduced expression of various κB-regulated gene targets. Consequent to this, western blot analysis of proteins corresponding to the NF-κB pathway indicated that AEA inhibited phosphorylation levels of the inhibitor of κB-kinase complex, resulting in the elimination of apoptotic resistance originating from NF-κB activation. This AEA-based apoptotic modulation was elucidated for the first time in this study, and gave rise to the proposal of an NF-κB model termed the 'Switching/Alternating Model.' In addition to this, AEA was also found to synergistically enhance the proapoptotic effects of paclitaxel, when used in combination with MCF-7 cells, presumably by a chemosensitizing role. Therefore, it was concluded that AEA isolated from the Malaysian tropical ginger (A. conchigera) served as a very promising candidate for further in-vivo development in animal models and in subsequent clinical trials involving patients with breast-related malignancies.
FOXP2 shares partially overlapping normal tissue expression and functionality with FOXP1; an established diffuse large B-cell lymphoma (DLBCL) oncogene and marker of poor prognosis. FOXP2 is expressed in the plasma cell malignancy multiple myeloma but has not been studied in DLBCL, where a poor prognosis activated B-cell (ABC)-like subtype display partially blocked plasma cell differentiation. FOXP2 protein expression was detected in ABC-DLBCL cell lines, and in primary DLBCL samples tumoral FOXP2 protein expression was detected in both germinal center B-cell-like (GCB) and non-GCB DLBCL. In biopsies from DLBCL patients treated with immunochemotherapy (R-CHOP), ≥ 20% nuclear tumoral FOXP2-positivity (n = 24/158) correlated with significantly inferior overall survival (OS: P = 0.0017) and progression-free survival (PFS: P = 0.0096). This remained significant in multivariate analysis against either the international prognostic index score or the non-GCB DLBCL phenotype (P < 0.05 for both OS and PFS). Expression of BLIMP1, a marker of plasmacytic differentiation that is commonly inactivated in ABC-DLBCL, did not correlate with patient outcome or FOXP2 expression in this series. Increased frequency of FOXP2 expression significantly correlated with FOXP1-positivity (P = 0.0187), and FOXP1 co-immunoprecipitated FOXP2 from ABC-DLBCL cells indicating that these proteins can co-localize in a multi-protein complex. FOXP2-positive DLBCL had reduced expression of HIP1R (P = 0.0348), which is directly repressed by FOXP1, and exhibited distinct patterns of gene expression. Specifically in ABC-DLBCL these were associated with lower expression of immune response and T-cell receptor signaling pathways. Further studies are warranted to investigate the potential functional cooperativity between FOXP1 and FOXP2 in repressing immune responses during the pathogenesis of high-risk DLBCL.
Nasopharyngeal carcinoma (NPC) is a solid tumor of the head and neck. Multimodal therapy is highly effective when NPC is detected early. However, due to the location of the tumor and the absence of clinical signs, early detection is difficult, making a biomarker for the early detection of NPC a priority. The dysregulation of small non-coding RNAs (miRNAs) during carcinogenesis is the focus of much current biomarker research. Herein, we examine several miRNA discovery methods using two sample matrices to identify circulating miRNAs (c-miRNAs) associated with NPC.
We previously found cytotoxic effects of tomato leaf extract (TLE) on the MCF-7 breast cancer cell line. The aim of this study was to ascertain the molecular mechanisms associated with the usage of TLE as an anticancer agent by microarray analysis using mRNA from MCF-7 breast cancer cells after treatment with TLE for 1 hr and 48 hrs. Approximately 991 genes out of the 30,000 genes in the human genome were significantly (p<0.05) changed after the treatment. Within this gene set, 88 were significantly changed between the TLE treated cells and the untreated MCF-7 cells (control cells) with a cut-off fold change >2.00. In order to focus on genes that were involved in cancer cell growth, only twenty-nine genes were selected, either down-regulated or up-regulated after treatment with TLE. Microarray assay results were confirmed by analyzing 10 of the most up and down regulated genes related to cancer cells progression using real-time PCR. Treatment with TLE induced significant up-regulation in the expression of the CRYAB, PIM1, BTG1, CYR61, HIF1-α and CEBP-β genes after 1 hr and 48 hrs, whereas the TXNIP and THBS1 genes were up-regulated after 1 hr of treatment but down-regulated after 48 hrs. In addition both the HMG1L1 and HIST2H3D genes were down-regulated after 1 hr and 48 hrs of treatment. These results demonstrate the potent activity of TLE as an anticancer agent.
Treatment options for patients with brain metastases (BMs) have limited efficacy and the mortality rate is virtually 100%. Targeted therapy is critically under-utilized, and our understanding of mechanisms underpinning metastatic outgrowth in the brain is limited. To address these deficiencies, we investigated the genomic and transcriptomic landscapes of 36 BMs from breast, lung, melanoma and oesophageal cancers, using DNA copy-number analysis and exome- and RNA-sequencing. The key findings were as follows. (a) Identification of novel candidates with possible roles in BM development, including the significantly mutated genes DSC2, ST7, PIK3R1 and SMC5, and the DNA repair, ERBB-HER signalling, axon guidance and protein kinase-A signalling pathways. (b) Mutational signature analysis was applied to successfully identify the primary cancer type for two BMs with unknown origins. (c) Actionable genomic alterations were identified in 31/36 BMs (86%); in one case we retrospectively identified ERBB2 amplification representing apparent HER2 status conversion, then confirmed progressive enrichment for HER2-positivity across four consecutive metastatic deposits by IHC and SISH, resulting in the deployment of HER2-targeted therapy for the patient. (d) In the ERBB/HER pathway, ERBB2 expression correlated with ERBB3 (r(2) = 0.496; p < 0.0001) and HER3 and HER4 were frequently activated in an independent cohort of 167 archival BM from seven primary cancer types: 57.6% and 52.6% of cases were phospho-HER3(Y1222) or phospho-HER4(Y1162) membrane-positive, respectively. The HER3 ligands NRG1/2 were barely detectable by RNAseq, with NRG1 (8p12) genomic loss in 63.6% breast cancer-BMs, suggesting a microenvironmental source of ligand. In summary, this is the first study to characterize the genomic landscapes of BM. The data revealed novel candidates, potential clinical applications for genomic profiling of resectable BMs, and highlighted the possibility of therapeutically targeting HER3, which is broadly over-expressed and activated in BMs, independent of primary site and systemic therapy.
Salmonella typhimurium is an important biofilm-forming bacteria. It is known to be resistant to a wide range of antimicrobials. The present study was carried out to evaluate the effects of dimethyl sulfoxide (DMSO) against S. typhimurium biofilm and investigate whole-cell protein expression by biofilm cells following treatment with DMSO. Antibiofilm activities were assessed using pellicle assay, crystal violet assay, colony-forming unit counting and extracellular polymeric substance (EPS) matrix assay whilst differential protein expression was investigated using a combination of one dimensional sodium dodecyl sulfate polyacrylamide gel electrophoresis, tandem mass spectrometry and bioinformatics. Treatment with 32% DMSO inhibited pellicle formation, biofilm viability, biofilm biomass and several important components of EPS matrix. Subtractive protein profiling identified two unique protein bands (25.4 and 51.2 kDa) which were present only in control biofilm and not in 32% DMSO-treated biofilm. In turn, 29 and 46 proteins were successfully identified from the protein bands of 25.4 and 51.2 kDa respectively. Protein interaction network analysis identified several biological pathways to be affected, including glycolysis, PhoP-PhoQ phosphorelay signalling and flagellar biosynthesis. The present study suggests that DMSO may inhibit multiple biological pathways to control biofilm formation.
In this day and age, the expectation of cosmetic products to effectively slow down skin photoaging is constantly increasing. However, the detrimental effects of UVB on the skin are not easy to tackle as UVB dysregulates a wide range of molecular changes on the cellular level. In our research, irradiated keratinocyte cells not only experienced a compromise in their redox system, but processes from RNA translation to protein synthesis and folding were also affected. Aside from this, proteins involved in various other processes like DNA repair and maintenance, glycolysis, cell growth, proliferation, and migration were affected while the cells approached imminent cell death. Additionally, the collagen degradation pathway was also activated by UVB irradiation through the upregulation of inflammatory and collagen degrading markers. Nevertheless, with the treatment of Swietenia macrophylla (S. macrophylla) seed extract and fractions, the dysregulation of many genes and proteins by UVB was reversed. The reversal effects were particularly promising with the S. macrophylla hexane fraction (SMHF) and S. macrophylla ethyl acetate fraction (SMEAF). SMHF was able to oppose the detrimental effects of UVB in several different processes such as the redox system, DNA repair and maintenance, RNA transcription to translation, protein maintenance and synthesis, cell growth, migration and proliferation, and cell glycolysis, while SMEAF successfully suppressed markers related to skin inflammation, collagen degradation, and cell apoptosis. Thus, in summary, our research not only provided a deeper insight into the molecular changes within irradiated keratinocytes, but also serves as a model platform for future cosmetic research to build upon. Subsequently, both SMHF and SMEAF also displayed potential photoprotective properties that warrant further fractionation and in vivo clinical trials to investigate and obtain potential novel bioactive compounds against photoaging.
Hepatitis B virus (HBV) barely induces host interferon (IFN)-stimulated genes (ISGs), which allows efficient HBV replication in the immortalized mouse hepatocytes as per human hepatocytes. Here we found that transfection of Isg20 plasmid robustly inhibits the HBV replication in HBV-infected hepatocytes irrespective of IRF3 or IFN promoter activation. Transfection of Isg20 is thus effective to eradicate HBV in the infected hepatocytes. Transfection of HBV genome or ε-stem of HBV pgRNA (active pgRNA moiety) failed to induce Isg20 in the hepatocytes, while control polyI:C (a viral dsRNA analogue mimic) activated MAVS pathway leading to production of type I IFN and then ISGsg20 via the IFN-α/β receptor (IFNAR). Consistently, addition of IFN-α induced Isg20 and partially suppressed HBV replication in hepatocytes. Chasing HBV RNA, DNA and proteins by blotting indicated that ISG20 expression decreased HBV RNA and replicative DNA in HBV-transfected cells, which resulted in low HBs antigen production and virus titer. The exonuclease domains of ISG20 mainly participated in HBV-RNA decay. In vivo hydrodynamic injection, ISG20 was crucial for suppressing HBV replication without degrading host RNA in the liver. Taken together, ISG20 acts as an innate anti-HBV effector that selectively degrades HBV RNA and blocks replication of infectious HBV particles. ISG20 would be a critical effector for ameliorating chronic HBV infection in the IFN therapy.
Plants have a remarkable ability to perceive and respond to various wavelengths of light and initiate regulation of different cascades of light signaling and molecular components. While the perception of red light and the mechanisms of its signaling involving phytochromes are largely known, knowledge of the mechanisms of blue light signaling is still limited. Chemical genetics involves the use of diverse small active or synthetic molecules to evaluate biological processes. By combining chemicals and analyzing the effects they have on plant morphology, we identified a chemical, 3-bromo-7-nitroindazole (3B7N), that promotes hypocotyl elongation of wild-type Arabidopsis only under continuous blue light. Further evaluation with loss-of-function mutants confirmed that 3B7N inhibits photomorphogenesis through cryptochrome-mediated light signaling. Microarray analysis demonstrated that the effect of 3B7N treatment on gene expression in cry1cry2 is considerably smaller than that in the wild type, indicating that 3B7N specifically interrupts cryptochrome function in the control of seedling development in a light-dependent manner. We demonstrated that 3B7N directly binds to CRY1 protein using an in vitro binding assay. These results suggest that 3B7N is a novel chemical that directly inhibits plant cryptochrome function by physical binding. The application of 3B7N can be used on other plants to study further the blue light mechanism and the genetic control of cryptochromes in the growth and development of plant species.
Single-walled carbon nanotubes (SWCNTs) are carbon-based nanomaterials that possess immense industrial potential. Despite accumulating evidence that exposure to SWCNTs might be toxic to humans, our understanding of the mechanisms for cellular toxicity of SWCNTs remain limited. Here, we demonstrated that acute exposure of short (1-3μm) and regular-length (5-30μm) pristine, carboxylated or hydroxylated SWCNTs inhibited cell proliferation in human somatic and human stem cells in a cell type-dependent manner. The toxicity of regular-length pristine SWCNT was most evidenced in NP69>CYT00086>MCF-10A>MRC-5>HaCaT > HEK-293T>HepG2. In contrast, the short pristine SWCNTs were relatively less toxic in most of the cells being tested, except for NP69 which is more sensitive to short pristine SWCNTs as compared to regular-length pristine SWCNTs. Interestingly, carboxylation and hydroxylation of regular-length SWCNTs, but not the short SWCNTs, significantly reduced the cytotoxicity. Exposure of SWCNTs also induced caspase 3 and 9 activities, mitochondrial membrane depolarization, and significant apoptosis and necrosis in MRC-5 embryonic lung fibroblasts. In contrast, SWCNTs inhibited the proliferation of HaCaT human keratinocytes without inducing cell death. Further analyses by gene expression profiling and Connectivity Map analysis showed that SWCNTs induced a gene expression signature characteristic of heat shock protein 90 (HSP90) inhibition in MRC-5 cells, suggesting that SWCNTs may inhibit the HSP90 signaling pathway. Indeed, exposure of MRC-5 cells to SWCNTs results in a dose-dependent decrease in HSP90 client proteins (AKT, CDK4 and BCL2) and a concomitant increase in HSP70 expression. In addition, SWCNTs also significantly inhibited HSP90-dependent protein refolding. Finally, we showed that ectopic expression of HSP90, but not HSP40 or HSP70, completely abrogated the cytotoxic effects of SWCNTs, suggesting that SWCNT-induced cellular toxicity is HSP90 dependent. In summary, our findings suggest that the toxic effects of SWCNTs are mediated through inhibition of HSP90 in human lung fibroblasts and keratinocytes.
Nasopharyngeal carcinoma (NPC) is an epithelial cancer of the nasopharynx which is highly associated with Epstein-Barr virus (EBV). Worldwide, most of the top 20 countries with the highest incidence and mortality rates of NPC are low- and middle-income countries. Many studies had demonstrated that EBV could be detected in the tissue, serum and plasma of NPC patients. In this study, we explored the potential of assays based on non-invasive nasal washings (NW) as a diagnostic and prognostic tool for NPC. A total of 128 patients were evaluated for NW EBV DNA loads and a subset of these samples were also tested for 27 EBV and human miRNAs shortlisted from literature. EBV DNA and seven miRNAs showed area under the receiver operating characteristic curve (AUC) values of more than 0.7, suggestive of their potential utility to detect NPC. Logistic regression analyses suggested that combination of two NW assays that test for EBNA-1 and hsa-miR-21 had the best performance in detecting NPC. The trend of NW EBV DNA load matched with clinical outcome of 71.4% (10 out of 14) NPC patients being followed-up. In summary, the non-invasive NW testing panel may be particularly useful for NPC screening in remote areas where healthcare facilities and otolaryngologists are lacking, and may encourage frequent testing of individuals in the high risk groups who are reluctant to have their blood tested. However, further validation in an independent cohort is required to strengthen the utility of this testing panel as a non-invasive detection tool for NPC.
Cisplatin is the first-line chemotherapeutic agent for the treatment of oral squamous cell carcinoma (OSCC). However, the intrinsic or acquired resistance against cisplatin remains a major obstacle to treatment efficacy in OSCC. Recently, mitochondrial DNA (mtDNA) alterations have been reported in a variety of cancers. However, the role of mtDNA alterations in OSCC has not been comprehensively studied. In this study, we evaluated the correlation between mtDNA alterations (mtDNA content, point mutations, large-scale deletions, and methylation status) and cisplatin sensitivity using two OSCC cell lines, namely SAS and H103, and stem cell-like tumour spheres derived from SAS. By microarray analysis, we found that the tumour spheres profited from aberrant lipid and glucose metabolism and became resistant to cisplatin. By qPCR analysis, we found that the cells with less mtDNA were less responsive to cisplatin (H103 and the tumour spheres). Based on the findings, we theorised that the metabolic changes in the tumour spheres probably resulted in mtDNA depletion, as the cells suppressed mitochondrial respiration and switched to an alternative mode of energy production, i.e. glycolysis. Then, to ascertain the origin of the variation in mtDNA content, we used MinION, a nanopore sequencer, to sequence the mitochondrial genomes of H103, SAS, and the tumour spheres. We found that the lower cisplatin sensitivity of H103 could have been caused by a constellation of genetic and epigenetic changes in its mitochondrial genome. Future work may look into how changes in mtDNA translate into an impact on cell function and therefore cisplatin response.
Long non-coding RNAs (lncRNAs) have been shown to possess a wide range of functions in both cellular and developmental processes including cancers. Although some of the lncRNAs have been implicated in the regulation of the immune response, the exact function of the large majority of lncRNAs still remains unknown. In this study, we characterized the lncRNAs in human primary monocytes, an essential component of the innate immune system. We performed RNA sequencing of monocytes from four individuals and combined our data with eleven other publicly available datasets. Our analysis led to identification of ~8000 lncRNAs of which >1000 have not been previously reported in monocytes. PCR-based validation of a subset of the identified novel long intergenic noncoding RNAs (lincRNAs) revealed distinct expression patterns. Our study provides a landscape of lncRNAs in monocytes, which could facilitate future experimental studies to characterize the functions of these molecules in the innate immune system.
Ornithine decarboxylase (ODC) is an enzyme of one of the two pathways of putrescine biosynthesis in plants. The genes encoding ODC have previously been cloned from Datura stramonium and human. Using differential screening, we isolated ODC cDNA clone from a cDNA library of ripening Capsicum annuum fruit. The cDNA clone designated CUKM10 contains an insert of 1523 bp. The longest open reading frame potentially encodes a peptide of 345 amino acids with an estimated molecular mass of 47 kDa and exhibit striking similarity to other ODCs. Expression analysis showed that the capODC hybridised to a single transcript with a size of 1.7 kb. The capODC transcript was first observed in early ripening and increased steadily until it reached fully ripening stage. From the observation it is suggested that capODC is developmentally regulated especially during later stage of ripening.
The aims of the present study were to undertake gene expression profiling of the blood of glioma patients to determine key genetic components of signaling pathways and to develop a panel of genes that could be used as a potential blood-based biomarker to differentiate between high and low grade gliomas, non-gliomas and control samples. In this study, blood samples were obtained from glioma patients, non-glioma and control subjects. Ten samples each were obtained from patients with high and low grade tumours, respectively, ten samples from non-glioma patients and twenty samples from control subjects. Total RNA was isolated from each sample after which first and second strand synthesis was performed. The resulting cRNA was then hybridized with the Agilent Whole Human Genome (4x44K) microarray chip according to the manufacturer's instructions. Universal Human Reference RNA and samples were labeled with Cy3 CTP and Cy5 CTP, respectively. Microarray data were analyzed by the Agilent Gene Spring 12.1V software using stringent criteria which included at least a 2-fold difference in gene expression between samples. Statistical analysis was performed using the unpaired Student's t-test with a p<0.01. Pathway enrichment was also performed, with key genes selected for validation using droplet digital polymerase chain reaction (ddPCR). The gene expression profiling indicated that were a substantial number of genes that were differentially expressed with more than a 2-fold change (p<0.01) between each of the four different conditions. We selected key genes within significant pathways that were analyzed through pathway enrichment. These key genes included regulators of cell proliferation, transcription factors, cytokines and tumour suppressor genes. In the present study, we showed that key genes involved in significant and well established pathways, could possibly be used as a potential blood-based biomarker to differentiate between high and low grade gliomas, non-gliomas and control samples.
Tissue-specific transcriptional programs control most biological phenotypes, including disease states such as cancer. However, the molecular details underlying transcriptional specificity is largely unknown, hindering the development of therapeutic approaches. Here, we describe novel experimental reporter systems that allow interrogation of the endogenous expression of HIF2A, a critical driver of renal oncogenesis. Using a focused CRISPR-Cas9 library targeting chromatin regulators, we provide evidence that these reporter systems are compatible with high-throughput screening. Our data also suggests redundancy in the control of cancer type-specific transcriptional traits. Reporter systems such as those described here could facilitate large-scale mechanistic dissection of transcriptional programmes underlying cancer phenotypes, thus paving the way for novel therapeutic approaches.
Haloxylon persicum is an endangered western Asiatic desert plant species, which survives under extreme environmental conditions. In this study, we focused on transcriptome analysis of H. persicum to understand the molecular mechanisms associated with drought tolerance. Two different periods of polyethylene glycol (PEG)-induced drought stress (48 h and 72 h) were imposed on H. persicum under in vitro conditions, which resulted in 18 million reads, subsequently assembled by de novo method with more than 8000 transcripts in each treatment. The N50 values were 1437, 1467, and 1524 for the control sample, 48 h samples, and 72 h samples, respectively. The gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis resulted in enrichment of mitogen-activated protein kinase (MAPK) and plant hormone signal transduction pathways under PEG-induced drought conditions. The differential gene expression analysis (DGEs) revealed significant changes in the expression pattern between the control and the treated samples. The KEGG analysis resulted in mapping transcripts with 138 different pathways reported in plants. The differential expression of drought-responsive transcription factors depicts the possible signaling cascades involved in drought tolerance. The present study provides greater insight into the fundamental transcriptome reprogramming of desert plants under drought.
Salt loading (SL) and water deprivation (WD) are experimental challenges that are often used to study the osmotic circuitry of the brain. Central to this circuit is the supraoptic nucleus (SON) of the hypothalamus, which is responsible for the biosynthesis of the hormones, arginine vasopressin (AVP) and oxytocin (OXT), and their transport to terminals that reside in the posterior lobe of the pituitary. On osmotic challenge evoked by a change in blood volume or osmolality, the SON undergoes a function-related plasticity that creates an environment that allows for an appropriate hormone response. Here, we have described the impact of SL and WD compared with euhydrated (EU) controls in terms of drinking and eating behavior, body weight, and recorded physiological data including circulating hormone data and plasma and urine osmolality. We have also used microarrays to profile the transcriptome of the SON following SL and remined data from the SON that describes the transcriptome response to WD. From a list of 2,783 commonly regulated transcripts, we selected 20 genes for validation by qPCR. All of the 9 genes that have already been described as expressed or regulated in the SON by osmotic stimuli were confirmed in our models. Of the 11 novel genes, 5 were successfully validated while 6 were false discoveries.
The Rab GTPase family plays a vital role in several plant physiological processes including fruit ripening. Fruit softening during ripening involves trafficking of cell wall polymers and enzymes between cellular compartments. Mango, an economically important fruit crop, is known for its delicious taste, exotic flavour and nutritional value. So far, there is a paucity of information on the mango Rab GTPase family. In this study, 23 genes encoding Rab proteins were identified in mango by a comprehensive in silico approach. Sequence alignment and similarity tree analysis with the model plant Arabidopsis as a reference enabled the bona fide assignment of the deduced mango proteins to classify into eight subfamilies. Expression analysis by RNA-Sequencing (RNA-Seq) showed that the Rab genes were differentially expressed in ripe and unripe mangoes suggesting the involvement of vesicle trafficking during ripening. Interaction analysis showed that the proteins involved in vesicle trafficking and cell wall softening were interconnected providing further evidence of the involvement of the Rab GTPases in fruit softening. Correlation analyses showed a significant relationship between the expression level of the RabA3 and RabA4 genes and fruit firmness at the unripe stage of the mango varieties suggesting that the differences in gene expression level might be associated with the contrasting firmness of these varieties. This study will not only provide new insights into the complexity of the ripening-regulated molecular mechanism but also facilitate the identification of potential Rab GTPases to address excessive fruit softening.
Due to the limitations in the range of antibodies recognising avian viruses, quantitative real-time PCR (RT-qPCR) is still the most widely used method to evaluate the expression of immunologically related genes in avian viruses. The objective of this study was to identify suitable reference genes for mRNA expression analysis in chicken intraepithelial lymphocyte natural killer (IEL-NK) cells after infection with very-virulent infectious bursal disease virus (vvIBDV). Fifteen potential reference genes were selected based on the references available. The coefficient of variation percentage (CV%) and average count of these 15 genes were determined by NanoString technology for control and infected samples. The M and V values for shortlisted reference genes (ACTB, GAPDH, HMBS, HPRT1, SDHA, TUBB1 and YWHAZ) were calculated using geNorm and NormFinder. GAPDH, YWHAZ and HMBS were the most stably expressed genes. The expression levels of three innate immune response related target genes, CASP8, IL22 and TLR3, agreed in the NanoString and RNA sequencing (RNA-Seq) results using one or two reference genes for normalisation (not HMBS). In conclusion, GAPDH and YWHAZ could be used as reference genes for the normalisation of chicken IEL-NK cell gene responses to infection with vvIBDV.