Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR

Overexpression of the immediate-early genes Egr1, Egr2, and Egr3 in two strains of rodents susceptible to audiogenic seizures

Genetic animal models of epilepsy are an important tool for further understanding the basic cellular mechanisms underlying epileptogenesis and for developing novel antiepileptic drugs. We conducted a comparative study of gene expression in the inferior colliculus, a nucleus that triggers audiogenic seizures, using two animal models, the Wistar audiogenic rat (WAR) and the genetic audiogenic seizure hamster (GASH:Sal). For this purpose, both models were exposed to high intensity auditory stimulation, and 60min later, the inferior colliculi were collected. As controls, intact Wistar rats and Syrian hamsters were subjected to stimulation and tissue preparation protocols identical to those performed on the experimental animals. Ribonucleic acid was isolated, and microarray analysis comparing the stimulated Wistar and WAR rats showed that the genomic profile of these animals displayed significant (fold change, |FC|≥2.0 and p<0.05) upregulation of 38 genes and downregulation of 47 genes. Comparison of gene expression profiles between stimulated control hamsters and stimulated GASH:Sal revealed the upregulation of 10 genes and the downregulation of 5 genes. Among the common genes that were altered in both models, we identified the zinc finger immediate-early growth response gene Egr3. The Egr3 protein is a transcription factor that is induced by distinct stress-elicited factors. Based on immunohistochemistry, this protein was expressed in the cochlear nucleus complex, the inferior colliculus, and the hippocampus of both animal models as well as in lymphoma tumors of the GASH:Sal. Our results support that the overexpression of the Egr3 gene in both models might contribute to neuronal viability and development of lymphoma in response to stress associated with audiogenic seizures. This article is part of a Special Issue entitled "Genetic and Reflex Epilepsies, Audiogenic Seizures and Strains: From Experimental Models to the Clinic".

Comprehensive processing of high-throughput small RNA sequencing data including quality checking, normalization, and differential expression analysis using the UEA sRNA Workbench

Recently, high-throughput sequencing (HTS) has revealed compelling details about the small RNA (sRNA) population in eukaryotes. These 20 to 25 nt noncoding RNAs can influence gene expression by acting as guides for the sequence-specific regulatory mechanism known as RNA silencing. The increase in sequencing depth and number of samples per project enables a better understanding of the role sRNAs play by facilitating the study of expression patterns. However, the intricacy of the biological hypotheses coupled with a lack of appropriate tools often leads to inadequate mining of the available data and thus, an incomplete description of the biological mechanisms involved. To enable a comprehensive study of differential expression in sRNA data sets, we present a new interactive pipeline that guides researchers through the various stages of data preprocessing and analysis. This includes various tools, some of which we specifically developed for sRNA analysis, for quality checking and normalization of sRNA samples as well as tools for the detection of differentially expressed sRNAs and identification of the resulting expression patterns. The pipeline is available within the UEA sRNA Workbench, a user-friendly software package for the processing of sRNA data sets. We demonstrate the use of the pipeline on a H. sapiens data set; additional examples on a B. terrestris data set and on an A. thaliana data set are described in the Supplemental Information A comparison with existing approaches is also included, which exemplifies some of the issues that need to be addressed for sRNA analysis and how the new pipeline may be used to do this.

Functional Gene Analysis Reveals Cell Cycle Changes and Inflammation in Endothelial Cells Irradiated with a Single X-ray Dose

Background and Purpose: Epidemiological data suggests an excess risk of cardiovascular disease (CVD) at low doses (0.05 and 0.1 Gy) of ionizing radiation (IR). Furthermore, the underlying biological and molecular mechanisms of radiation-induced CVD are still unclear. Because damage to the endothelium could be critical in IR-related CVD, this study aimed to identify the effects of radiation on immortalized endothelial cells in the context of atherosclerosis.

Material and Methods: Microarrays and RT-qPCR were used to compare the response of endothelial cells irradiated with a single X-ray dose (0.05, 0.1, 0.5, 2 Gy) measured after various post-irradiation (repair) times (1 day, 7 days, 14 days). To consolidate and mechanistically support the endothelial cell response to X-ray exposure identified via microarray analysis, DNA repair signaling (γH2AX/TP53BP1-foci quantification), cell cycle progression (BrdU/7AAD flow cytometric analysis), cellular senescence (β-galactosidase assay with CPRG and IGFBP7 quantification) and pro-inflammatory status (IL6 and CCL2) was assessed.

Results: Microarray results indicated persistent changes in cell cycle progression and inflammation. Cells underwent G1 arrest in a dose-dependent manner after high doses (0.5 and 2 Gy), which was compensated by increased proliferation after 1 week and almost normalized after 2 weeks. However, at this point irradiated cells showed an increased β-Gal activity and IGFBP7 secretion, indicative of premature senescence. The production of pro-inflammatory cytokines IL6 and CCL2 was increased at early time points.

Conclusions: IR induces pro-atherosclerotic processes in endothelial cells in a dose-dependent manner. These findings give an incentive for further research on the shape of the dose-response curve, as we show that even low doses of IR can induce premature endothelial senescence at later time points. Furthermore, our findings on the time- and dose-dependent response regarding differentially expressed genes, cell cycle progression, inflammation and senescence bring novel insights into the underlying molecular mechanisms of the endothelial response to X-ray radiation. This may in turn lead to the development of risk-reducing strategies to prevent IR-induced CVD, such as the use of cell cycle modulators and anti-inflammatory drugs as radioprotectors and/or radiation mitigators.

Sodium Octanoate Modulates the Innate Immune Response of Bovine Mammary Epithelial Cells through the TLR2/P38/JNK/ERK1/2 Pathway: Implications during Staphylococcus aureus Internalization

Bovine mammary epithelial cells (bMECs) contribute to mammary gland defense against invading pathogens, such as Staphylococcus aureus (intracellular facultative), which is recognized by TLR2. In a previous report, we showed that sodium octanoate [NaO, a medium chain fatty acid (C8)] induces (0.25 mM) or inhibits (1 mM) S. aureus internalization into bMECs and differentially regulates the innate immune response (IIR). However, the molecular mechanisms have not been described, which was the aim of this study. The results showed that α5β1 integrin membrane abundance (MA) was increased in 0.25 mM NaO-treated cells, but TLR2 or CD36 MA was not modified. When these receptors were blocked individually, 0.25 mM NaO-increased S. aureus internalization was notably reduced. Interestingly, in this condition, the IIR of the bMECs was impaired because MAPK (p38, JNK, and ERK1/2) phosphorylation and the activation of transcription factors related to these pathways were decreased. In addition, the 1 mM NaO treatment induced TLR2 MA, but neither the integrin nor CD36 MA was modified. The reduction in S. aureus internalization induced by 1 mM NaO was increased further when TLR2 was blocked. In addition, the phosphorylation levels of the MAPKs increased, and 13 transcriptional factors related to the IIR were slightly activated (CBF, CDP, c-Myb, AP-1, Ets-1/Pea-3, FAST-1, GAS/ISRE, AP-2, NFAT-1, OCT-1, RAR/DR-5, RXR/DR-1, and Stat-3). Moreover, the 1 mM NaO treatment up-regulated gene expression of IL-8 and RANTES and secretion of IL-1β. Notably, when 1 mM NaO-treated bMECs were challenged with S. aureus, the gene expression of IL-8 and IL-10 increased, while IL-1β secretion was reduced. In conclusion, our results showed that α5β1 integrin, TLR2 and CD36 are involved in 0.25 mM NaO-increased S. aureus internalization in bMECs. In addition, 1 mM NaO activates bMECs via the TLR2 signaling pathways (p38, JNK, and ERK1/2), which improves IIR before S. aureus invasion. Additionally, NaO (1 mM) might exert anti-inflammatory effects after bacterial internalization.

Different genotypes of Trypanosoma cruzi produce distinctive placental environment genetic response in chronic experimental infection

Congenital infection of Trypanosoma cruzi allows transmission of this parasite through generations. Despite the problematic that this entails, little is known about the placenta environment genetic response produced against infection. We performed functional genomics by microarray analysis in C57Bl/6J mice comparing placentas from uninfected animals and from animals infected with two different T. cruzi strains: K98, a clone of the non-lethal myotropic CA-I strain (TcI), and VD (TcVI), isolated from a human case of congenital infection. Analysis of networks by GeneMANIA of differentially expressed genes showed that “Secretory Granule” was a pathway down-regulated in both infected groups, whereas “Innate Immune Response” and “Response to Interferon-gamma” were pathways up-regulated in VD infection but not in K98. Applying another approach, the GSEA algorithm that detects small changes in predetermined gene sets, we found that metabolic processes, transcription and macromolecular transport were down-regulated in infected placentas environment and some pathways related to cascade signaling had opposite regulation: over-represented in VD and down-regulated in K98 group. We also have found a stronger tropism to the placental organ by VD strain, by detection of parasite DNA and RNA, suggesting living parasites. Our study is the first one to describe in a murine model the genetic response of placental environment to T. cruzi infection and suggests the development of a strong immune response, parasite genotype-dependent, to the detriment of cellular metabolism, which may contribute to control infection preventing the risk of congenital transmission.

Congenital transmission of Trypanosoma cruzi, the causative agent of Chagas disease, remains a problem of global public health impact in endemic areas where vectorial and transfusional transmission have been controlled and in non-endemic countries due to migration movements. Little is known about how the parasite´s presence and genetic variability affect placental capacity to protect the fetus. This study explores, for the first time, the effects of placental environment infection by analyzing parasite persistence and gene expression using functional genomics and biological network analyses in mice infected by two strains of T. cruzi, with differential capacity of congenital transmission. The infection with the strain with a stronger placental tropism was associated to a higher degree of up-regulation in genes related to innate immunity and response to interferon-gamma. Our findings suggest that the placental environment exerts a strong immune response in detriment of cellular metabolism modulated by the parasite strain. These findings constitute a significant contribution to better understand the mechanisms causing congenital infection of T. cruzi.

Dataset of proinflammatory cytokine and cytokine receptor gene expression in rainbow trout (Oncorhynchus mykiss) measured using a novel GeXP multiplex, RT-PCR assay

A GeXP multiplex, RT-PCR assay was developed and optimized that simultaneously measures expression of a suite of immune-relevant genes in rainbow trout (Oncorhynchus mykiss), concentrating on tumor necrosis factor and interleukin-1 ligand/receptor systems and acute phase response genes. The dataset includes expression values for drpt, il11a, il1b1, il1b2, il1b3, il1r-like-1(e3-5), il1r-like-1(e9-11), il1r1-like-a, il1r1-like-b, il1r2, saa, tnfa1, tnfa2, tnfa3, tnfrsf1a, tnfrsf1a-like-a, tnfrsf1a-like-b, tnfrsf5, and tnfrsf9. Gene expression was measured at four time-points post-challenge in both a resistant line (ARS-Fp-R) and a susceptible line (ARS-Fp-S) of rainbow trout. In addition, fish body weight, spleen index and the Flavobacterium psychrophilum load are reported. These data are an extension of information presented and discussed in “Proinflammatory cytokine and cytokine receptor gene expression kinetics following challenge with Flavobacterium psychrophilum in resistant and susceptible lines of rainbow trout (Oncorhynchus mykiss)” (Kutyrev et al., 2016) [1].

Application of unique sequence index (USI) barcode to gene expression profiling in gastric adenocarcinoma

Accurate expression profiling is imperative for understanding the biological roles of mRNAs. Real-time PCR have been at the forefront of biological innovation in detection and monitoring of gene expression, however, fluorophore-labeled oligonucleotides and double-stranded DNA binding dyes, the two most frequently used dyes in RNA detection, are not very cost effective and have poor specificity, respectively. We have developed a cost effective and specific approach for mRNA expression profiling via added unique sequence index (USI) to cDNAs before amplification. USI is a barcode which enable the detection of each target RNA. Using this method, caudal type homeobox 1 (CDX1) and FAT atypical cadherin 4 (FAT4) expressions were investigated in tumoral and non-tumoral tissues of gastric cancer patients and compared with commercial ABI kit. Both methods indicated that FAT4 and CDX1 expression were significantly reduced in gastric cancer tissues compared with adjacent noncancerous tissues. Moreover, we have shown that this assay is highly sensitive, linear and reproducible. USI barcode not only provides a powerful tool for mRNA detection due to its sensitivity, specificity and cost-effectiveness, but also allows comfortable design for real-time qPCR assays within the least time and empowers the analysis of many transcripts of virtually any organism. Furthermore, USI barcode is highly affordable for large numbers of different samples or small sample sizes without microarray and expensive commercial platforms.

A Systems Approach Implicates a Brain Mitochondrial Oxidative Homeostasis Co-expression Network in Genetic Vulnerability to Alcohol Withdrawal

Genetic factors significantly affect vulnerability to alcohol dependence (alcoholism). We previously identified quantitative trait loci on distal mouse chromosome 1 with large effects on predisposition to alcohol physiological dependence and associated withdrawal following both chronic and acute alcohol exposure in mice (Alcdp1 and Alcw1, respectively). We fine-mapped these loci to a 1.1–1.7 Mb interval syntenic with human 1q23.2-23.3. Alcw1/Alcdp1 interval genes show remarkable genetic variation among mice derived from the C57BL/6J and DBA/2J strains, the two most widely studied genetic animal models for alcohol-related traits. Here, we report the creation of a novel recombinant Alcw1/Alcdp1 congenic model (R2) in which the Alcw1/Alcdp1 interval from a donor C57BL/6J strain is introgressed onto a uniform, inbred DBA/2J genetic background. As expected, R2 mice demonstrate significantly less severe alcohol withdrawal compared to wild-type littermates. Additionally, comparing R2 and background strain animals, as well as reciprocal congenic (R8) and appropriate background strain animals, we assessed Alcw1/Alcdp1 dependent brain gene expression using microarray and quantitative PCR analyses. To our knowledge this includes the first Weighted Gene Co-expression Network Analysis using reciprocal congenic models. Importantly, this allows detection of co-expression patterns limited to one or common to both genetic backgrounds with high or low predisposition to alcohol withdrawal severity. The gene expression patterns (modules) in common contain genes related to oxidative phosphorylation, building upon human and animal model studies that implicate involvement of oxidative phosphorylation in alcohol use disorders (AUDs). Finally, we demonstrate that administration of N-acetylcysteine, an FDA-approved antioxidant, significantly reduces symptoms of alcohol withdrawal (convulsions) in mice, thus validating a phenotypic role for this network. Taken together, these studies support the importance of mitochondrial oxidative homeostasis in alcohol withdrawal and identify this network as a valuable therapeutic target in human AUDs.

Gastrointestinal Transcriptomic Response of Metabolic Vitamin B12 Pathways in Roux-en-Y Gastric Bypass

OBJECTIVES

Vitamin B12 (B12) deficiency after Roux-en-Y gastric bypass (RYGB) is highly prevalent and may contribute to postoperative complications. Decreased production of intrinsic factor owing to gastric fundus removal is thought to have a major role, but other components of B12 metabolism may also be affected. We evaluated changes in the expression levels of multiple B12 pathway-encoding genes in gastrointestinal (GI) tissues to evaluate the potential roles in contributing to post-RYGB B12 deficiency.

METHODS

During double-balloon enteroscopy, serial GI biopsies were collected from 20 obese women (age, 46.9±6.2 years; body mass index, 46.5±5.3 kg/m2) with adult-onset type 2 diabetes (fasting plasma glucose ≥126 mg/dl; hemoglobin A1c≥6.5%) before and, at the same site, 3 months after RYGB. Gene expression levels were assessed by the Affymetrix Human GeneChip 1.0 ST microarray. Findings were validated by real-time quantitative PCR (RT-qPCR).

RESULTS

Gene expression levels with significant changes (P≤0.05) included: transcobalamin I (TCN1) in remnant (-1.914-fold) and excluded (-1.985-fold) gastric regions; gastric intrinsic factor (GIF) in duodenum (-0.725-fold); and cubilin (CUBN) in duodenum (+0.982-fold), jejunum (+1.311-fold), and ileum (+0.685-fold). Validation by RT-qPCR confirmed (P≤0.05) observed changes for TCN1 in the remnant gastric region (-0.132-fold) and CUBN in jejunum (+2.833-fold).

CONCLUSIONS

RYGB affects multiple pathway-encoding genes that may be associated with postoperative B12 deficiency. Decreased TCN1 levels seem to be the main contributing factor. Increased CUBN levels suggest an adaptive genetic reprogramming of intestinal tissue aiming to compensate for impaired intestinal B12 delivery.

Transcriptional Responses to Pre-flowering Leaf Defoliation in Grapevine Berry from Different Growing Sites, Years, and Genotypes

Leaf removal is a grapevine canopy management technique widely used to modify the source-sink balance and/or microclimate around berry clusters to optimize fruit composition. In general, the removal of basal leaves before flowering reduces fruit set, hence achieving looser clusters, and improves grape composition since yield is generally curtailed more than proportionally to leaf area itself. Albeit responses to this practice seem quite consistent, overall vine performance is affected by genotype, environmental conditions, and severity of treatment. The physiological responses of grape varieties to defoliation practices have been widely investigated, and just recently a whole genome transcriptomic approach was exploited showing an extensive transcriptome rearrangement in berries defoliated before flowering. Nevertheless, the extent to which these transcriptomic reactions could be manifested by different genotypes and growing environments is entirely unexplored. To highlight general responses to defoliation vs. different locations, we analyzed the transcriptome of cv. Sangiovese berries sampled at four development stages from pre-flowering defoliated vines in two different geographical areas of Italy. We obtained and validated five markers of the early defoliation treatment in Sangiovese, an ATP-binding cassette transporter, an auxin response factor, a cinnamyl alcohol dehydrogenase, a flavonoid 3-O-glucosyltransferase and an indole-3-acetate beta-glucosyltransferase. Candidate molecular markers were also obtained in another three grapevine genotypes (Nero d'Avola, Ortrugo, and Ciliegiolo), subjected to the same level of selective pre-flowering defoliation (PFD) over two consecutive years in their different areas of cultivation. The flavonol synthase was identified as a marker in the pre-veraison phase, the jasmonate methyltransferase during the transition phase and the abscisic acid receptor PYL4 in the ripening phase. The characterization of transcriptome changes in Sangiovese berry after PFD highlights, on one hand, the stronger effect of environment than treatment on the whole berry transcriptome rearrangement during development and, on the other, expands existing knowledge of the main molecular and biochemical modifications occurring in defoliated vines. Moreover, the identification of candidate genes associated with PFD in different genotypes and environments provides new insights into the applicability and repeatability of this crop practice, as well as its possible agricultural and qualitative outcomes across genetic and environmental variability.

Transcriptome-Wide Identification of Differentially Expressed Genes in Solanum lycopersicon L. in Response to an Alfalfa-Protein Hydrolysate Using Microarrays

An alfalfa-based protein hydrolysate (EM) has been tested in tomato (Solanum lycopersicon L.) plants at two different concentrations (0.1 and 1 mL L^-1) to get insight on its efficacy as biostimulant in this species and to unravel possible metabolic targets and molecular mechanisms that may shed light on its mode of action. EM was efficient in promoting the fresh biomass and content in chlorophyll and soluble sugars of tomato plants, especially when it was applied at the concentration of 1 mL L^-1. This effect on plant productivity was likely related to the EM-dependent up-regulation of genes identified via microarray and involved in primary carbon and nitrogen metabolism, photosynthesis, nutrient uptake and developmental processes. EM also up-regulated a number of genes implied in the secondary metabolism that leads to the synthesis of compounds (phenols and terpenes) functioning in plant development and interaction with the environment. Concomitantly, phenol content was enhanced in EM-treated plants. Several new genes have been identified in tomato as potential targets of EM action, like those involved in detoxification processes from reactive oxygen species and xenobiotic (particularly glutathione/ascorbate cycle-related and ABC transporters), and defense against abiotic and biotic stress. The model hypothesized is that elicitors present in the EM formulation like auxins, phenolics, and amino acids, may trigger a signal transduction pathway via modulation of the intracellular levels of the hormones ethylene, jasmonic acid and abscissic acid, which then further prompt the activation of a cascade events requiring the presence and activity of many kinases and transcription factors to activate stress-related genes. The genes identified suggest these kinases and transcription factors as players involved in a complex crosstalk between biotic and abiotic stress signaling pathways. We conclude that EM acts as a biostimulant in tomato due to its capacity to stimulate plant productivity and up-regulate stress-related responses. Its use in agricultural practices may reduce the need of inorganic fertilizers and pesticides, thereby reducing the environmental impact of productive agriculture.

Finding lost genes in GWAS via integrative-omics analysis reveals novel sub-networks associated with preterm birth

Maternal genome influences associate with up to 40% of spontaneous preterm births (PTB). Multiple genome wide association studies (GWAS) have been completed to identify genetic variants associated with PTB. Disappointingly, no highly significant SNPs have replicated in independent cohorts so far. We developed an approach combining protein-protein interaction (PPI) network data with tissue specific gene expression data to "find" SNPs of modest significance to identify candidate genes of functional importance that would otherwise be overlooked. This approach is based on the assumption that "high-ranking" SNPs falling short of genome wide significance may nevertheless indicate genes that have substantial biological value in understanding PTB. We mapped highly-ranked candidate SNPs from a meta-analysis of PTB-GWAS to coding genes and developed a PPI network enriched with PTB-SNP carrying genes. This network was scored with gene expression data from term and preterm myometrium to identify subnetworks of PTB-SNP associated genes coordinately expressed with labour onset in myometrial tissue. Our analysis consistently identified significant sub-networks associated with the interacting transcription factors MEF2C and TWIST1, genes not previously associated with PTB, both of which regulate processes clearly relevant to birth timing. Other genes in the significant sub-networks were also associated with inflammatory pathways, as well as muscle function and ion channels. Gene expression level dysregulation was confirmed for eight of these networks by qRT-PCR in an independent set of term and pre-term subjects. Our method identifies novel genes dysregulated in PTB and provides a generalized framework to identify GWAS SNPs that would otherwise be overlooked.

miR-152 down-regulation is associated with MET up-regulation in leiomyosarcoma and undifferentiated pleomorphic sarcoma

PURPOSE

Highly aggressive adult soft tissue sarcomas (STS), i.e., leiomyosarcomas (LMS) and undifferentiated pleomorphic sarcomas (UPS), present complex genomic anomalies and overall 5-year survival rates of 20 to 40%. Here, we aimed to identify new biomarkers that may be employed to improve the treatment of non-translocation STS patients. We validated 12 miRNAs implicated in tumor development using primary STS samples and selected miR-152 for further analysis in STS-derived cell lines.

METHODS

59 primary STS samples (27 LMS and 32 UPS) and 10 matched normal control tissues were included in the study, as well as 3 STS-derived cell lines (HT1080, SW872 and SKLMS1) and a normal control mesenchymal cell line (hMSC). miRNA expression analyses were performed using a TaqMan microRNA Array platform and qRT-PCR (miR-152), respectively. The expression levels of the putative miR-152 targets MET and KIT were assessed using qRT-PCR and immunohistochemistry on tissue microarrays, respectively. In addition, various functional analyses were performed before and after miR-152 transfection into SKLMS1 cells.

RESULTS

We found that 12 pre-selected miRNAs were down-regulated in primary STS tumor samples compared to its normal control samples. A statistically significant miR-152 down-regulation was found to be accompanied by high MET and KIT mRNA levels in both the primary samples and the STS-derived cell lines tested. miR-152 transfection in SKLMS1 cells led to a reduction in KIT and MET mRNA and protein levels which, in turn, was associated with a transient down-regulation of the PI3K/AKT pathway, a transient decrease in cell growth, and a transient increase in both apoptotic and S-phase cells.

CONCLUSIONS

Our data indicate that over-expression of MET and KIT in primary STS samples and its derived cell lines is associated with miR-152 down-regulation. This shift may play a role in STS development and, thus, may be used to identify patients at risk. The effect of MET down-regulation on downstream signaling pathways, such as the PI3K/AKT pathway, may provide a basis for the future design of novel STS treatment strategies.

Selection of Suitable Reference Genes for Analysis of Salivary Transcriptome in Non-Syndromic Autistic Male Children

Childhood autism is a severe form of complex genetically heterogeneous and behaviorally defined set of neurodevelopmental diseases, collectively termed as autism spectrum disorders (ASD). Reverse transcriptase quantitative real-time PCR (RT-qPCR) is a highly sensitive technique for transcriptome analysis, and it has been frequently used in ASD gene expression studies. However, normalization to stably expressed reference gene(s) is necessary to validate any alteration reported at the mRNA level for target genes. The main goal of the present study was to find the most stable reference genes in the salivary transcriptome for RT-qPCR analysis in non-syndromic male childhood autism. Saliva samples were obtained from nine drug naïve non-syndromic male children with autism and also sex-, age-, and location-matched healthy controls using the RNA-stabilizer kit from DNA Genotek. A systematic two-phased measurement of whole saliva mRNA levels for eight common housekeeping genes (HKGs) was carried out by RT-qPCR, and the stability of expression for each candidate gene was analyzed using two specialized algorithms, geNorm and NormFinder, in parallel. Our analysis shows that while the frequently used HKG ACTB is not a suitable reference gene, the combination of GAPDH and YWHAZ could be recommended for normalization of RT-qPCR analysis of salivary transcriptome in non-syndromic autistic male children.

MicroRNA-877-5p is involved in the trovafloxacin-induced liver injury

Trovafloxacin develops severe hepatotoxicity; however, the underlying mechanism of the trovafloxacin-induced liver injury has not been cleared. It has been shown that microRNAs (miRNAs) can be involved in the development of drug-induced liver injuries. We performed a miRNA microarray analysis to identify hepatic miRNAs that were induced or reduced by trovafloxacin in mice. It was demonstrated that miR-877-5p was the most increased miRNA in the mouse liver 24h after the trovafloxacin administration. To investigate the role of miR-877-5p in the liver, we established miR-877-5p-overexpressed HepG2 cells. Microarray analysis detected altered expressions in 2077 (>2-fold) and 1547 (<0.5-fold) genes in the miR-877-5p overexpressing cells compared to the mock cells. Especially, SLCO4C1, PEPCK, MT1M, HIST1H2BM, LGI1, and PLA2G2A were markedly increased or decreased in the miR-877-5p overexpressing cells. We conducted a correlation analysis between the expression levels of miR-877-5p and the six genes in eight miR-877-5p stably-expressed clones. It was shown that the PEPCK expression levels were correlated with miR-877-5p expression levels. PEPCK is associated with development of apoptotic cell death; therefore, the increased miR- 877-5p-induced PEPCK can be a trigger that is involved in the development of trovafloxacin-induced liver injury.

Abnormal expression of inflammatory genes in placentas of women with sickle cell anemia and sickle hemoglobin C disease

Sickle cell disease (SCD) is a complex disease that is characterized by the polymerization of deoxyhemoglobin S, altered red blood cell membrane biology, endothelial activation, hemolysis, a procoagulant state, acute and chronic inflammation, and vaso-occlusion. Among the physiological changes that occur during pregnancy, oxygen is consumed by fetal growth, and pregnant women with SCD are more frequently exposed to low oxygen levels. This might lead to red blood cells sickling, and, consequently, to vaso-occlusion. The mechanisms by which SCD affects placental physiology are largely unknown, and chronic inflammation might be involved in this process. This study aimed to evaluate the gene expression profile of inflammatory response mediators in the placentas of pregnant women with sickle cell cell anemia (HbSS) and hemoglobinopathy SC (HbSC). Our results show differences in a number of these genes. For the HbSS group, when compared to the control group, the following genes showed differential expression: IL1RAP (2.76-fold), BCL6 (4.49-fold), CXCL10 (-2.12-fold), CXCR1 (-3.66-fold), and C3 (-2.0-fold). On the other hand, the HbSC group presented differential expressions of the following genes, when compared to the control group: IL1RAP (4.33-fold), CXCL1 (3.05-fold), BCL6 (4.13-fold), CXCL10 (-3.32-fold), C3 (-2.0-fold), and TLR3 (2.38-fold). Taken together, these data strongly suggest a differential expression of several inflammatory genes in both SCD (HbSS and HbSC), indicating that the placenta might become an environment with hypoxia, and increased inflammation, which could lead to improper placental development.

Nutritional Evaluation of an EPA-DHA Oil from Transgenic Camelina sativa in Feeds for Post-Smolt Atlantic Salmon (Salmo salar L.)

Vegetable oils (VO) are possible substitutes for fish oil in aquafeeds but their use is limited by their lack of omega-3 (n-3) long-chain polyunsaturated fatty acids (LC-PUFA). However, oilseed crops can be modified to produce n-3 LC-PUFA such as eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids, representing a potential option to fill the gap between supply and demand of these important nutrients. Camelina sativa was metabolically engineered to produce a seed oil with around 15% total n-3 LC-PUFA to potentially substitute for fish oil in salmon feeds. Post-smolt Atlantic salmon (Salmo salar) were fed for 11-weeks with one of three experimental diets containing either fish oil (FO), wild-type Camelina oil (WCO) or transgenic Camelina oil (DCO) as added lipid source to evaluate fish performance, nutrient digestibility, tissue n-3 LC-PUFA, and metabolic impact determined by liver transcriptome analysis. The DCO diet did not affect any of the performance or health parameters studied and enhanced apparent digestibility of EPA and DHA compared to the WCO diet. The level of total n-3 LC-PUFA was higher in all the tissues of DCO-fed fish than in WCO-fed fish with levels in liver similar to those in fish fed FO. Endogenous LC-PUFA biosynthetic activity was observed in fish fed both the Camelina oil diets as indicated by the liver transcriptome and levels of intermediate metabolites such as docosapentaenoic acid, with data suggesting that the dietary combination of EPA and DHA inhibited desaturation and elongation activities. Expression of genes involved in phospholipid and triacylglycerol metabolism followed a similar pattern in fish fed DCO and WCO despite the difference in n-3 LC-PUFA contents.

Accumulation of Extracellular Matrix in Advanced Lesions of Canine Distemper Demyelinating Encephalitis

In demyelinating diseases, changes in the quality and quantity of the extracellular matrix (ECM) may contribute to demyelination and failure of myelin repair and axonal sprouting, especially in chronic lesions. To characterize changes in the ECM in canine distemper demyelinating leukoencephalitis (DL), histochemical and immunohistochemical investigations of formalin-fixed paraffin-embedded cerebella using azan, picrosirius red and Gomori`s silver stain as well as antibodies directed against aggrecan, type I and IV collagen, fibronectin, laminin and phosphacan showed alterations of the ECM in CDV-infected dogs. A significantly increased amount of aggrecan was detected in early and late white matter lesions. In addition, the positive signal for collagens I and IV as well as fibronectin was significantly increased in late lesions. Conversely, the expression of phosphacan was significantly decreased in early and more pronounced in late lesions compared to controls. Furthermore, a set of genes involved in ECM was extracted from a publically available microarray data set and was analyzed for differential gene expression. Gene expression of ECM molecules, their biosynthesis pathways, and pro-fibrotic factors was mildly up-regulated whereas expression of matrix remodeling enzymes was up-regulated to a relatively higher extent. Summarized, the observed findings indicate that changes in the quality and content of ECM molecules represent important, mainly post-transcriptional features in advanced canine distemper lesions. Considering the insufficiency of morphological regeneration in chronic distemper lesions, the accumulated ECM seems to play a crucial role upon regenerative processes and may explain the relatively small regenerative potential in late stages of this disease.

A 3-methylcrotonyl-CoA carboxylase deficient human skin fibroblast transcriptome reveals underlying mitochondrial dysfunction and oxidative stress

Isolated 3-methylcrotonyl-CoA carboxylase (MCC) deficiency is an autosomal recessive inherited metabolic disease of leucine catabolism with a highly variable phenotype. Apart from extensive mutation analyses of the MCCC1 and MCCC2 genes encoding 3-methylcrotonyl-CoA carboxylase (EC 6.4.1.4), molecular data on MCC deficiency gene expression studies in human tissues is lacking. For IEMs, unbiased '-omics' approaches are starting to reveal the secondary cellular responses to defects in biochemical pathways. Here we present the first whole genome expression profile of immortalized cultured skin fibroblast cells of two clinically affected MCC deficient patients and two healthy individuals generated using Affymetrix(®)HuExST1.0 arrays. There were 16191 significantly differentially expressed transcript IDs of which 3591 were well annotated and present in the predefined knowledge database of Ingenuity Pathway Analysis software used for downstream functional analyses. The most noticeable feature of this MCCA deficient skin fibroblast transcriptome was the typical genetic hallmark of mitochondrial dysfunction, decreased antioxidant response and disruption of energy homeostasis, which was confirmed by mitochondrial functional analyses. The MCC deficient transcriptome seems to predict oxidative stress that could alter the complex secondary cellular response that involve genes of the glycolysis, the TCA cycle, OXPHOS, gluconeogenesis, β-oxidation and the branched-chain fatty acid metabolism. An important emerging insight from this human MCCA transcriptome in combination with previous reports is that chronic exposure to the primary and secondary metabolites of MCC deficiency and the resulting oxidative stress might impact adversely on the quality of life and energy levels, irrespective of whether MCC deficient individuals are clinically affected or asymptomatic.

Transcriptomic Analysis of Drosophila Mushroom Body Neurons Lacking Amyloid-β Precursor-Like Protein Activity

The amyloid-β protein precursor (AβPP) is subjected to sequential intramembrane proteolysis by α-, β-, andγ-secretases, producing secreted amyloid-β (Aβ) peptides and a cytoplasmically released AβPP Intracellular Domain (AICD). AICD complexes with transcription factors in the nucleus, suggesting that this AβPP fragment serves as an active signaling effector that regulates downstream genes, although its nuclear targets are poorly defined. To further understand this potential signaling mechanism mediated by AβPP, we performed a transcriptomic identification of the Drosophila genome that is regulated by the fly AβPP orthologue in fly mushroom body neurons, which control learning- and memory-based behaviors. We find significant changes in expression of 245 genes, representing approximately 1.6% of the Drosophila genome, with the changes ranging from +6 fold to -40 fold. The largest class of responsive targets corresponds to non-protein coding genes and includes microRNAs that have been previously implicated in Alzheimer's disease pathophysiology. Several genes were identified in our Drosophila microarray analyses that have also emerged as putative AβPP targets in similar mammalian transcriptomic studies. Our results also indicate a role for AβPP in cellular pathways involving the regulation of Drosophila Casein Kinase II, mitochondrial oxidative phosphorylation, RNA processing, and innate immunity. Our findings provide insights into the intracellular events that are regulated by AβPP activity in healthy neurons and that might become dysregulated as a result of abnormal AβPP proteolysis in AD.

Maternal Whole Blood Gene Expression at 18 and 28 Weeks of Gestation Associated with Spontaneous Preterm Birth in Asymptomatic Women

The heterogeneity of spontaneous preterm birth (SPTB) requires an interdisciplinary approach to determine potential predictive risk factors of early delivery. The aim of this study was to investigate maternal whole blood gene expression profiles associated with spontaneous preterm birth (SPTB, <37 weeks) in asymptomatic pregnant women. The study population was a matched subgroup of women (51 SPTBs, 114 term delivery controls) who participated in the All Our Babies community based cohort in Calgary (n = 1878). Maternal blood at 17–23 (sampling time point 1, T₁) and 27–33 weeks of gestation (T₂) were collected. Total RNA was extracted and microarray was performed on 326 samples (165 women). Univariate analyses determined significant clinical factors and differential gene expression associated with SPTB. Thirteen genes were validated using qRT-PCR. Three multivariate logistic models were constructed to identify gene expression at T₁ (Model A), T₂ (Model B), and gene expression fold change from T₁ to T₂ (Model C) associated with SPTB. All models were adjusted for clinical factors. Model C can predict SPTB with 65% sensitivity and 88% specificity in asymptomatic women after adjusting for history of abortion and anaemia (occurring before T₂). Clinical data enhanced the sensitivity of the Models to predict SPTB. In conclusion, clinical factors and whole blood gene expression are associated with SPTB in asymptomatic women. An effective screening tool for SPTB during pregnancy would enable targeted preventive approaches and personalised antenatal care.

The defensin from avocado (Persea americana var. drymifolia) PaDef induces apoptosis in the human breast cancer cell line MCF-7

Antimicrobial peptides (AMPs) are cytotoxic to cancer cells; however, mainly the effects of AMPs from animals have been evaluated. In this work, we assessed the cytotoxicity of PaDef defensin from avocado (Persea americana var. drymifolia) on the MCF-7 cancer cell line (a breast cancer cell line) and evaluated its mechanism of action. PaDef inhibited the viability of MCF-7 cells in a concentration-dependent manner, with an IC50=141.62μg/ml. The viability of normal peripheral blood mononuclear cells was unaffected by this AMP. Additionally, PaDef induced apoptosis in MCF-7 cells in a time-dependent manner, but did not affect the membrane potential or calcium flow. In addition, PaDef IC50 induced the expression of cytochrome c, Apaf-1, and the caspase 7 and 9 genes. Likewise, this defensin induced the loss of mitochondrial Δψm and increased the phosphorylation of MAPK p38, which may lead to MCF-7 apoptosis by the intrinsic pathway. This is the first report of an avocado defensin inducing intrinsic apoptosis in cancer cells, which suggests that it could be a potential therapeutic molecule in the treatment of cancer.

Ciliary neurotrophic factor upregulates follistatin and Pak1, causes overexpression of muscle differentiation related genes and downregulation of established atrophy mediators in skeletal muscle

INTRODUCTION

The Ciliary Neurotrophic Factor (CNTF) is a pluripotent cytokine with anorexigenic actions in the hypothalamus that improves insulin sensitivity, increases energy expenditure and induces weight loss. Since CNTF also has an established myotrophic role, we sought to examine whether skeletal muscle contributes to the CNTF-induced metabolic improvement and identify the molecular mechanisms mediating these effects.

METHODS

We used a mouse model of diet-induced obesity, to which high or low CNTF doses were administered for 7days. Whole transcriptome expression levels were analyzed in dissected soleus muscles using microarrays and data were then confirmed using qRT-PCR.

RESULTS

We demonstrate that CNTF administration significantly downregulates leptin, while it upregulates follistatin and Pak1; a molecule associated with insulin sensitization in skeletal muscle. A significant overexpression of muscle differentiation related genes and downregulation of established atrophy mediators was observed.

CONCLUSIONS

The overall gene expression changes suggest an indirect, beneficial effect of CNTF on metabolism, energy expenditure and insulin sensitivity, exerted by the pronounced stimulation of muscle growth, with similarities to the described effect of follistatin and the activation of the Akt pathway in skeletal muscle.

Altered gene expression in the lower respiratory tract of Car6 (-/-) mice

From birth, the respiratory tract mucosa is exposed to various chemical, physical, and microbiological stress factors. Efficient defense mechanisms and strictly regulated renewal systems in the mucosa are thus required. Carbonic anhydrase VI (CA VI) is the only secreted isoenzyme of the α-CA gene family. It is transported in high concentrations in saliva and milk into the alimentary tract where it contributes to optimal pH homeostasis. Earlier study of transcriptomic responses of Car6 (-/-) mice has shown changes in the response to oxidative stress and brown fat cell differentiation in the submandibular gland. It has been suggested that CA VI delivered to the mucosal surface of the bronchiolar epithelium is an essential factor in defense and renewal of the lining epithelium. In this study, the transcriptional effects of CA VI deficiency were investigated in both trachea and lung of Car6 (-/-) mice using a cDNA microarray analysis. Functional clustering of the results indicated significant changes of gene transcription in the lower airways. The altered biological processes included antigen transport by M-cells, potassium transport, muscle contraction, and thyroid hormone synthesis. Immunohistochemical staining confirmed the absence of CA VI in the submandibular gland of Car6 (-/-) mice. Immunostaining of the trachea and lung samples revealed no differences between the knockout and wild type groups nor were any morphological changes observed. The present findings can help us to recognize novel functions for CA VI-one of the major protein constituents of saliva and milk.

Co-expression network analysis reveals transcription factors associated to cell wall biosynthesis in sugarcane

Sugarcane is a hybrid of Saccharum officinarum and Saccharum spontaneum, with minor contributions from other species in Saccharum and other genera. Understanding the molecular basis of cell wall metabolism in sugarcane may allow for rational changes in fiber quality and content when designing new energy crops. This work describes a comparative expression profiling of sugarcane ancestral genotypes: S. officinarum, S. spontaneum and S. robustum and a commercial hybrid: RB867515, linking gene expression to phenotypes to identify genes for sugarcane improvement. Oligoarray experiments of leaves, immature and intermediate internodes, detected 12,621 sense and 995 antisense transcripts. Amino acid metabolism was particularly evident among pathways showing natural antisense transcripts expression. For all tissues sampled, expression analysis revealed 831, 674 and 648 differentially expressed genes in S. officinarum, S. robustum and S. spontaneum, respectively, using RB867515 as reference. Expression of sugar transporters might explain sucrose differences among genotypes, but an unexpected differential expression of histones were also identified between high and low Brix° genotypes. Lignin biosynthetic genes and bioenergetics-related genes were up-regulated in the high lignin genotype, suggesting that these genes are important for S. spontaneum to allocate carbon to lignin, while S. officinarum allocates it to sucrose storage. Co-expression network analysis identified 18 transcription factors possibly related to cell wall biosynthesis while in silico analysis detected cis-elements involved in cell wall biosynthesis in their promoters. Our results provide information to elucidate regulatory networks underlying traits of interest that will allow the improvement of sugarcane for biofuel and chemicals production.

PIWI-interacting RNA (piRNA) signatures in human cardiac progenitor cells

Cardiac progenitors, such as cardiospheres and cardiosphere-derived cells, represent an attractive cell source for cardiac regeneration. The PIWI-interacting RNAs, piRNAs, are an intriguing class of small non-coding RNAs, implicated in the regulation of epigenetic state, maintenance of genomic integrity and stem cell functions. Although non-coding RNAs are an exploiting field in cardiovascular research, the piRNA signatures of cardiac progenitors has not been evaluated yet.We profiled, through microarrays, 15,311 piRNAs expressed in cardiospheres, cardiosphere-derived cells and cardiac fibroblasts. Results showed a set of differentially expressed piRNAs (fold change ≥2, p<0.01): 641 piRNAs were upregulated and 1,301 downregulated in the cardiospheres compared to cardiosphere-derived cells, while 255 and 708 piRNAs resulted up- and down-regulated, respectively, if compared to cardiac fibroblasts. We also identified 181 piRNAs that are overexpressed and 129 are downregulated in cardiosphere-derived cells respect to cardiac fibroblasts.Bioinformatics analysis showed that the deregulated piRNAs were mainly distributed on few chromosomes, suggesting that piRNAs are organized in discrete genomic clusters.Furthermore, the bioinformatics search showed that the most upregulated piRNAs target transposons, especially belonged to LINE-1 class, as validated by qRT-PCR. This reduction is also associated to an activation of AKT signaling, which is beneficial for cardiac regeneration.The present study is the first to show a highly consistent piRNA expression pattern for human cardiac progenitors, likely responsible of their different regenerative power. Moreover, this piRNome analysis may provide new methods for characterize cardiac progenitors and may shed new light on the understanding the complex molecular mechanisms of cardiac regeneration.

The effect of palmitate supplementation on gene expression profile in proliferating myoblasts

High-fat diet, exposure to saturated fatty acids, or the presence of adipocytes in myoblast microenvironment affects skeletal muscle growth and function. The aim of the present study was to investigate the effect of palmitate supplementation on transcriptomic profile of mouse C2C12 myoblasts. Global gene expression was evaluated using whole mouse genome oligonucleotide microarrays, and the results were validated through qPCR. A total of 4047 genes were identified as differentially expressed, including 3492 downregulated and 555 upregulated genes, during a 48-h exposure to palmitate (0.1 mmol/l). Functional classification showed the involvement of these genes in several processes which regulate cell growth. In conclusion, the addition of palmitate modifies the expression of genes associated with (1) myoblast responsiveness to hormones and growth factors, (2) cytokine and growth factor expression, and (3) regulation of cell-cell and cell-matrix communication. Such alterations can affect myoblast growth and differentiation; however, further studies in this field are required.

Oxidative Nanopatterning of Titanium Surface Influences mRNA and MicroRNA Expression in Human Alveolar Bone Osteoblastic Cells

Titanium implants have been extensively used in orthopedic and dental applications. It is well known that micro- and nanoscale surface features of biomaterials affect cellular events that control implant-host tissue interactions. To improve our understanding of how multiscale surface features affect cell behavior, we used microarrays to evaluate the transcriptional profile of osteoblastic cells from human alveolar bone cultured on engineered titanium surfaces, exhibiting the following topographies: nanotexture (N), nano+submicrotexture (NS), and rough microtexture (MR), obtained by modulating experimental parameters (temperature and solution composition) of a simple yet efficient chemical treatment with a H₂SO₄/H₂O₂ solution. Biochemical assays showed that cell culture proliferation augmented after 10 days, and cell viability increased gradually over 14 days. Among the treated surfaces, we observed an increase of alkaline phosphatase activity as a function of the surface texture, with higher activity shown by cells adhering onto nanotextured surfaces. Nevertheless, the rough microtexture group showed higher amounts of calcium than nanotextured group. Microarray data showed differential expression of 716 mRNAs and 32 microRNAs with functions associated with osteogenesis. Results suggest that oxidative nanopatterning of titanium surfaces induces changes in the metabolism of osteoblastic cells and contribute to the explanation of the mechanisms that control cell responses to micro- and nanoengineered surfaces.

Meta-analysis of organ-specific differences in the structure of the immune infiltrate in major malignancies

Anticancer immunosurveillance is one of the major endogenous breaks of tumor progression. Here, we analyzed gene expression pattern indicative of the presence of distinct leukocyte subtypes within four cancer types (breast cancer, colorectal carcinoma, melanoma, and non-small cell lung cancer) and 20 different microarray datasets corresponding to a total of 3471 patients. Multiple metagenes reflecting the presence of such immune cell subtypes were highly reproducible across distinct cohorts. Nonetheless, there were sizable differences in the correlation patterns among such immune-relevant metagenes across distinct malignancies. The reproducibility of the correlations among immune-relevant metagenes was highest in breast cancer (followed by colorectal cancer, non-small cell lung cancer and melanoma), reflecting the fact that mammary carcinoma has an intrinsically better prognosis than the three other malignancies. Among breast cancer patients, we found that the expression of a lysosomal enzyme-related metagene centered around ASAH1 (which codes for N-acylsphingosine amidohydrolase-1, also called acid ceramidase) exhibited a higher correlation with multiple immune-relevant metagenes in patients that responded to neoadjuvant chemotherapy than in non-responders. Altogether, this meta-analysis revealed novel organ-specific features of the immune infiltrate in distinct cancer types, as well as a strategy for defining new prognostic biomarkers.

Profiling Ethylene-Responsive Genes Expressed in the Latex of the Mature Virgin Rubber Trees Using cDNA Microarray

Ethylene is commonly used as a latex stimulant of Hevea brasiliensis by application of ethephon (chloro-2-ethylphosphonic acid); however, the molecular mechanism by which ethylene increases latex production is not clear. To better understand the effects of ethylene stimulation on the laticiferous cells of rubber trees, a latex expressed sequence tag (EST)-based complementary DNA microarray containing 2,973 unique genes (probes) was first developed and used to analyze the gene expression changes in the latex of the mature virgin rubber trees after ethephon treatment at three different time-points: 8, 24 and 48 h. Transcript levels of 163 genes were significantly altered with fold-change values ≥ 2 or ≤ –2 (q-value < 0.05) in ethephon-treated rubber trees compared with control trees. Of the 163 genes, 92 were up-regulated and 71 down-regulated. The microarray results were further confirmed using real-time quantitative reverse transcript-PCR for 20 selected genes. The 163 ethylene-responsive genes were involved in several biological processes including organic substance metabolism, cellular metabolism, primary metabolism, biosynthetic process, cellular response to stimulus and stress. The presented data suggest that the laticifer water circulation, production and scavenging of reactive oxygen species, sugar metabolism, and assembly and depolymerization of the latex actin cytoskeleton might play important roles in ethylene-induced increase of latex production. The results may provide useful insights into understanding the molecular mechanism underlying the effect of ethylene on latex metabolism of H. brasiliensis.

Anti-Inflammatory and Antimicrobial Effects of Estradiol in Bovine Mammary Epithelial Cells during Staphylococcus aureus Internalization

17β-Estradiol (E2), the predominant sexual hormone in females, is associated with the modulation of the innate immune response (IIR), and changes in its levels at parturition are related to intramammary infections, such as mastitis. In bovine mammary epithelial cells (bMECs), E2 regulates differentiation and proliferation, but its immunomodulatory functions have not been explored. Staphylococcus aureus is the predominant pathogen causing mastitis, which can persist intracellularly in bMECs. The aim of this work was to analyze whether E2 modulates the IIR of bMECs during S. aureus internalization. bMECs treated with E2 (50 pg/mL, 24 h) reduced bacteria internalization (~50%). The host receptors α5β1 and TLR2 do not participate in this reduction. However, E2 activates ERα and modulates the IIR reducing the S. aureus induced-mRNA expression of TNF-α (~50%) and IL-1β (90%). E2 also decreased the secretion of these cytokines as well as IL-6 production; however, in infected bMECs, E2 induced the secretion of IL-1β. Furthermore, E2 upregulates the expression of the antimicrobial peptides DEFB1, BNBD5, and psoriasin S100A7 (~5-, 3-, and 6-fold, resp.). In addition, E2 induced the production of antimicrobial compounds in bMEC culture medium, which, together with the modulation of the IIR, could be related to the reduction of S. aureus internalization.

Bilateral gene interaction hierarchy analysis of the cell death gene response emphasizes the significance of cell cycle genes following unilateral traumatic brain injury

BACKGROUND

Delayed or secondary cell death that is caused by a cascade of cellular and molecular processes initiated by traumatic brain injury (TBI) may be reduced or prevented if an effective neuroprotective strategy is employed. Microarray and subsequent bioinformatic analyses were used to determine which genes, pathways and networks were significantly altered 24 h after unilateral TBI in the rat. Ipsilateral hemi-brain, the corresponding contralateral hemi-brain, and naïve (control) brain tissue were used for microarray analysis.

RESULTS

Ingenuity Pathway Analysis showed cell death and survival (CD) to be a top molecular and cellular function associated with TBI on both sides of the brain. One major finding was that the overall gene expression pattern suggested an increase in CD genes in ipsilateral brain tissue and suppression of CD genes contralateral to the injury which may indicate an endogenous protective mechanism. We created networks of genes of interest (GOI) and ranked the genes by the number of direct connections each had in the GOI networks, creating gene interaction hierarchies (GIHs). Cell cycle was determined from the resultant GIHs to be a significant molecular and cellular function in post-TBI CD gene response.

CONCLUSIONS

Cell cycle and apoptosis signalling genes that were highly ranked in the GIHs and exhibited either the inverse ipsilateral/contralateral expression pattern or contralateral suppression were identified and included STAT3, CCND1, CCND2, and BAX. Additional exploration into the remote suppression of CD genes may provide insight into neuroprotective mechanisms that could be used to develop therapies to prevent cell death following TBI.

Characterising granuloma regression and liver recovery in a murine model of schistosomiasis japonica

For hepatic schistosomiasis the egg-induced granulomatous response and the development of extensive fibrosis are the main pathologies. We used a Schistosoma japonicum-infected mouse model to characterise the multi-cellular pathways associated with the recovery from hepatic fibrosis following clearance of the infection with the anti-schistosomal drug, praziquantel. In the recovering liver splenomegaly, granuloma density and liver fibrosis were all reduced. Inflammatory cell infiltration into the liver was evident, and the numbers of neutrophils, eosinophils and macrophages were significantly decreased. Transcriptomic analysis revealed the up-regulation of fatty acid metabolism genes and the identification of Peroxisome proliferator activated receptor alpha as the upstream regulator of liver recovery. The aryl hydrocarbon receptor signalling pathway which regulates xenobiotic metabolism was also differentially up-regulated. These findings provide a better understanding of the mechanisms associated with the regression of hepatic schistosomiasis.

Microarray Analyses Reveal Marked Differences in Growth Factor and Receptor Expression Between 8-Cell Human Embryos and Pluripotent Stem Cells

Previous microarray analyses of RNAs from 8-cell (8C) human embryos revealed a lack of cell cycle checkpoints and overexpression of core circadian oscillators and cell cycle drivers relative to pluripotent human stem cells [human embryonic stem cells/induced pluripotent stem (hES/iPS)] and fibroblasts, suggesting growth factor independence during early cleavage stages. To explore this possibility, we queried our combined microarray database for expression of 487 growth factors and receptors. Fifty-one gene elements were overdetected on the 8C arrays relative to hES/iPS cells, including 14 detected at least 80-fold higher, which annotated to multiple pathways: six cytokine family (CSF1R, IL2RG, IL3RA, IL4, IL17B, IL23R), four transforming growth factor beta (TGFB) family (BMP6, BMP15, GDF9, ENG), one fibroblast growth factor (FGF) family [FGF14(FH4)], one epidermal growth factor member (GAB1), plus CD36, and CLEC10A. 8C-specific gene elements were enriched (73%) for reported circadian-controlled genes in mouse tissues. High-level detection of CSF1R, ENG, IL23R, and IL3RA specifically on the 8C arrays suggests the embryo plays an active role in blocking immune rejection and is poised for trophectoderm development; robust detection of NRG1, GAB1, -2, GRB7, and FGF14(FHF4) indicates novel roles in early development in addition to their known roles in later development. Forty-four gene elements were underdetected on the 8C arrays, including 11 at least 80-fold under the pluripotent cells: two cytokines (IFITM1, TNFRSF8), five TGFBs (BMP7, LEFTY1, LEFTY2, TDGF1, TDGF3), two FGFs (FGF2, FGF receptor 1), plus ING5, and WNT6. The microarray detection patterns suggest that hES/iPS cells exhibit suppressed circadian competence, underexpression of early differentiation markers, and more robust expression of generic pluripotency genes, in keeping with an artificial state of continual uncommitted cell division. In contrast, gene expression patterns of the 8C embryo suggest that it is an independent circadian rhythm-competent equivalence group poised to signal its environment, defend against maternal immune rejection, and begin the rapid commitment events of early embryogenesis.

A Systematic Approach to Discover and Characterize Natural Plant Biostimulants

The use of natural plant biostimulants is proposed as an innovative solution to address the challenges to sustainable agriculture, to ensure optimal nutrient uptake, crop yield, quality, and tolerance to abiotic stress. However, the process of selection and characterization of plant biostimulant matrices is complex and involves a series of rigorous evaluations customized to the needs of the plant. Here, we propose a highly differentiated plant biostimulant development and production platform, which involves a combination of technology, processes, and know-how. Chemistry, biology and omic concepts are combined/integrated to investigate and understand the specific mode(s) of action of bioactive ingredients. The proposed approach allows to predict and characterize the function of natural compounds as biostimulants. By managing and analyzing massive amounts of complex data, it is therefore possible to discover, evaluate and validate new product candidates, thus expanding the uses of existing products to meet the emerging needs of agriculture.

Controlled hydrostatic pressure stress downregulates the expression of ribosomal genes in preimplantation embryos: a possible protection mechanism?

The efficiency of various assisted reproductive techniques can be improved by preconditioning the gametes and embryos with sublethal hydrostatic pressure treatment. However, the underlying molecular mechanism responsible for this protective effect remains unknown and requires further investigation. Here, we studied the effect of optimised hydrostatic pressure treatment on the global gene expression of mouse oocytes after embryonic genome activation. Based on a gene expression microarray analysis, a significant effect of treatment was observed in 4-cell embryos derived from treated oocytes, revealing a transcriptional footprint of hydrostatic pressure-affected genes. Functional analysis identified numerous genes involved in protein synthesis that were downregulated in 4-cell embryos in response to hydrostatic pressure treatment, suggesting that regulation of translation has a major role in optimised hydrostatic pressure-induced stress tolerance. We present a comprehensive microarray analysis and further delineate a potential mechanism responsible for the protective effect of hydrostatic pressure treatment.

Response of Arabidopsis thaliana Roots with Altered Lipid Transfer Protein (LTP) Gene Expression to the Clubroot Disease and Salt Stress

The clubroot disease of Brassicaceae is caused by the obligate biotrophic protist Plasmodiophora brassicae. The disease is characterized by abnormal tumorous swellings of infected roots that result in reduced drought resistance and insufficient distribution of nutrients, leading to reduced crop yield. It is one of the most damaging diseases among cruciferous crops worldwide. The acquisition of nutrients by the protist is not well understood. Gene expression profiles in Arabidopsis thaliana clubroots indicate that lipid transfer proteins (LTPs) could be involved in disease development or at least in adaptation to the disease symptoms. Therefore, the aim of the study was to examine the role of some, of the still enigmatic LTPs during clubroot development. For a functional approach, we have generated transgenic plants that overexpress LTP genes in a root specific manner or show reduced LTP gene expression. Our results showed that overexpression of some of the LTP genes resulted in reduced disease severity whereas the lipid content in clubs of LTP mutants seems to be unaffected. Additional studies indicate a role for some LTPs during salt stress conditions in roots of A. thaliana.

Profiles of microRNA networks in intestinal epithelial cells in a mouse model of colitis

Inflammatory bowel diseases (IBDs) accompany a critical loss of the frontline barrier function that is achieved primarily by intestinal epithelial cells (IECs). Although the gene-regulation pathways underlying these host-defense roles of IECs presumably are deranged during IBD pathogenesis, the quantitative and qualitative alterations of posttranscriptional regulators such as microRNAs (miRNAs) within the cells largely remain to be defined. We aimed to uncover the regulatory miRNA–target gene relationships that arise differentially in inflamed small- compared with large-IECs. Whereas IBD significantly increased the expression of only a few miRNA candidates in small-IECs, numerous miRNAs were upregulated in inflamed large-IECs. These marked alterations might explain why the large, as compared with small, intestine is more sensitive to colitis and shows more severe pathology in this experimental model of IBD. Our in-depth assessment of the miRNA–mRNA expression profiles and the resulting networks prompts us to suggest that miRNAs such as miR-1224, miR-3473a, and miR-5128 represent biomarkers that appear in large-IECs upon IBD development and co-operatively repress the expression of key anti-inflammatory factors. The current study provides insight into gene-regulatory networks in IECs through which dynamic rearrangement of the involved miRNAs modulates the gene expression–regulation machinery between maintaining and disrupting gastrointestinal homeostasis.

Differentiation of oligodendrocyte precursors is impaired in the prefrontal cortex in schizophrenia

The pathophysiology of schizophrenia involves disturbances of information processing across brain regions, possibly reflecting, at least in part, a disruption in the underlying axonal connectivity. This disruption is thought to be a consequence of the pathology of myelin ensheathment, the integrity of which is tightly regulated by oligodendrocytes. In order to gain insight into the possible neurobiological mechanisms of myelin deficit, we determined the messenger RNA (mRNA) expression profile of laser captured cells that were immunoreactive for 2', 3'-cyclic nucleotide 3'-phosphodiesterase (CNPase), a marker for oligodendrocyte progenitor cells (OPCs) in addition to differentiating and myelinating oligodendrocytes, in the white matter of the prefrontal cortex in schizophrenia subjects. Our findings pointed to the hypothesis that OPC differentiation might be impaired in schizophrenia. To address this hypothesis, we quantified cells that were immunoreactive for neural/glial antigen 2 (NG2), a selective marker for OPCs, and those that were immunoreactive for oligodendrocyte transcription factor 2 (OLIG2), an oligodendrocyte lineage marker that is expressed by OPCs and maturing oligodendrocytes. We found that the density of NG2-immunoreactive cells was unaltered, but the density of OLIG2-immunoreactive cells was significantly decreased in subjects with schizophrenia, consistent with the notion that OPC differentiation impairment may contribute to oligodendrocyte disturbances and thereby myelin deficits in schizophrenia.

Development and validation of a mixed-tissue oligonucleotide DNA microarray for Atlantic bluefin tuna, Thunnus thynnus (Linnaeus, 1758)

BACKGROUND

The largest of the tuna species, Atlantic bluefin tuna (Thunnus thynnus), inhabits the North Atlantic Ocean and the Mediterranean Sea and is considered to be an endangered species, largely a consequence of overfishing. T. thynnus aquaculture, referred to as fattening or farming, is a capture based activity dependent on yearly renewal from the wild. Thus, the development of aquaculture practices independent of wild resources can provide an important contribution towards ensuring security and sustainability of this species in the longer-term. The development of such practices is today greatly assisted by large scale transcriptomic studies.

RESULTS

We have used pyrosequencing technology to sequence a mixed-tissue normalised cDNA library, derived from adult T. thynnus. A total of 976,904 raw sequence reads were assembled into 33,105 unique transcripts having a mean length of 893 bases and an N50 of 870. Of these, 33.4% showed similarity to known proteins or gene transcripts and 86.6% of them were matched to the congeneric Pacific bluefin tuna (Thunnus orientalis) genome, compared to 70.3% for the more distantly related Nile tilapia (Oreochromis niloticus) genome. Transcript sequences were used to develop a novel 15 K Agilent oligonucleotide DNA microarray for T. thynnus and comparative tissue gene expression profiles were inferred for gill, heart, liver, ovaries and testes. Functional contrasts were strongest between gills and ovaries. Gills were particularly associated with immune system, signal transduction and cell communication, while ovaries displayed signatures of glycan biosynthesis, nucleotide metabolism, transcription, translation, replication and repair.

CONCLUSIONS

Sequence data generated from a novel mixed-tissue T. thynnus cDNA library provide an important transcriptomic resource that can be further employed for study of various aspects of T. thynnus ecology and genomics, with strong applications in aquaculture. Tissue-specific gene expression profiles inferred through the use of novel oligo-microarray can serve in the design of new and more focused transcriptomic studies for future research of tuna physiology and assessment of the welfare in a production environment.

Co-Transcriptomes of Initial Interactions In Vitro between Streptococcus Pneumoniae and Human Pleural Mesothelial Cells

Streptococcus pneumoniae (Spn) is a major causative organism of empyema, an inflammatory condition occurring in the pleural sac. In this study, we used human and Spn cDNA microarrays to characterize the transcriptional responses occurring during initial contact between Spn and a human pleural mesothelial cell line (PMC) in vitro. Using stringent filtering criteria, 42 and 23 Spn genes were up-and down-regulated respectively. In particular, genes encoding factors potentially involved in metabolic processes and Spn adherence to eukaryotic cells were up-regulated e.g. glnQ, glnA, aliA, psaB, lytB and nox. After Spn initial contact, 870 human genes were differentially regulated and the largest numbers of significant gene expression changes were found in canonical pathways for eukaryotic initiation factor 2 signaling (60 genes out of 171), oxidative phosphorylation (32/103), mitochondrial dysfunction (37/164), eIF4 and p70S6K signaling (28/142), mTOR signaling (27/182), NRF2-mediated oxidative stress response (20/177), epithelial adherens junction remodeling (11/66) and ubiquitination (22/254). The cellular response appeared to be directed towards host cell survival and defense. Spn did not activate NF-kB or phosphorylate p38 MAPK or induce cytokine production from PMC. Moreover, Spn infection of TNF-α pre-stimulated PMC inhibited production of IL-6 and IL-8 secretion by >50% (p<0.01). In summary, this descriptive study provides datasets and a platform for examining further the molecular mechanisms underlying the pathogenesis of empyema.

Genome-wide mapping of Hif-1α binding sites in zebrafish

Background

Hypoxia Inducible Factor (HIF) regulates a cascade of transcriptional events in response to decreased oxygenation, acting from the cellular to the physiological level. This response is evolutionarily conserved, allowing the use of zebrafish (Danio rerio) as a model for studying the hypoxic response. Activation of the hypoxic response can be achieved in zebrafish by homozygous null mutation of the von Hippel-Lindau (vhl) tumour suppressor gene. Previous work from our lab has focused on the phenotypic characterisation of this mutant, establishing the links between vhl mutation, the hypoxic response and cancer. To further develop fish as a model for studying hypoxic signalling, we examine the transcriptional profile of the vhl mutant with respect to Hif-1α. As our approach uses embryos consisting of many cell types, it has the potential to uncover additional HIF regulated genes that have escaped detection in analogous mammalian cell culture studies.

Results

We performed high-density oligonucleotide microarray analysis of the gene expression changes in von Hippel-Lindau mutant zebrafish, which identified up-regulation of well-known hypoxia response genes and down-regulation of genes primarily involved in lipid processing. To identify the dependency of these transcriptional changes on HIF, we undertook Chromatin Immunoprecipitation linked next generation sequencing (ChIP-seq) for the transcription factor Hypoxia Inducible Factor 1α (HIF-1α). We identified HIF-1α binding sites across the genome, with binding sites showing enrichment for an RCGTG motif, showing conservation with the mammalian hypoxia response element.

Conclusions

Transcriptome analysis of vhl mutant embryos detected activation of key hypoxia response genes seen in human cell models of hypoxia, but also suppression of many genes primarily involved in lipid processing. ChIP-seq analysis of Hif-1α binding sites unveiled an unprecedented number of loci, with a high proportion containing a canonical hypoxia response element. Whether these sites are functional remains unknown, nevertheless their frequent location near transcriptional start sites suggests functionality, and will allow for investigation into the potential hypoxic regulation of genes in their vicinity. We expect that our data will be an excellent starting point for analysis of both fish and mammalian gene regulation by HIF.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2169-x) contains supplementary material, which is available to authorized users.

The Effect of Acute and Chronic Social Stress on the Hippocampal Transcriptome in Mice

Psychogenic stress contributes to the formation of brain pathology. Using gene expression microarrays, we analyzed the hippocampal transcriptome of mice subjected to acute and chronic social stress of different duration. The longest period of social stress altered the expression of the highest number of genes and most of the stress-induced changes in transcription were reversible after 5 days of rest. Chronic stress affected genes involved in the functioning of the vascular system (Alas2, Hbb-b1, Hba-a2, Hba-a1), injury response (Vwf, Mgp, Cfh, Fbln5, Col3a1, Ctgf) and inflammation (S100a8, S100a9, Ctla2a, Ctla2b, Lcn2, Lrg1, Rsad2, Isg20). The results suggest that stress may affect brain functions through the stress-induced dysfunction of the vascular system. An important issue raised in our work is also the risk of the contamination of brain tissue samples with choroid plexus. Such contamination would result in a consistent up- or down-regulation of genes, such as Ttr, Igf2, Igfbp2, Prlr, Enpp2, Sostdc1, 1500015O10RIK (Ecrg4), Kl, Clic6, Kcne2, F5, Slc4a5, and Aqp1. Our study suggests that some of the previously reported, supposedly specific changes in hippocampal gene expression, may be a result of the inclusion of choroid plexus in the hippocampal samples.

Metabolomic and Gene Expression Profiles Exhibit Modular Genetic and Dietary Structure Linking Metabolic Syndrome Phenotypes in Drosophila

Genetic and environmental factors influence complex disease in humans, such as metabolic syndrome, and Drosophila melanogaster serves as an excellent model in which to test these factors experimentally. Here we explore the modularity of endophenotypes with an in-depth reanalysis of a previous study by Reed et al. (2014), where we raised 20 wild-type genetic lines of Drosophila larvae on four diets and measured gross phenotypes of body weight, total sugar, and total triglycerides, as well as the endophenotypes of metabolomic and whole-genome expression profiles. We then perform new gene expression experiments to test for conservation of phenotype-expression correlations across different diets and populations. We find that transcript levels correlated with gross phenotypes were enriched for puparial adhesion, metamorphosis, and central energy metabolism functions. The specific metabolites L-DOPA and N-arachidonoyl dopamine make physiological links between the gross phenotypes across diets, whereas leucine and isoleucine thus exhibit genotype-by-diet interactions. Between diets, we find low conservation of the endophenotypes that correlate with the gross phenotypes. Through the follow-up expression study, we found that transcript-trait correlations are well conserved across populations raised on a familiar diet, but on a novel diet, the transcript-trait correlations are no longer conserved. Thus, physiological canalization of metabolic phenotypes breaks down in a novel environment exposing cryptic variation. We cannot predict the physiological basis of disease in a perturbing environment from profiles observed in the ancestral environment. This study demonstrates that variation for disease traits within a population is acquired through a multitude of physiological mechanisms, some of which transcend genetic and environmental influences, and others that are specific to an individual’s genetic and environmental context.

Alteration of Gene Expression Profile in Kidney of Spontaneously Hypertensive Rats Treated with Protein Hydrolysate of Blue Mussel (Mytilus edulis) by DNA Microarray Analysis

Marine organisms are rich sources of bioactive components, which are often reported to have antihypertensive effects. However, the underlying mechanisms have yet to be fully identified. The aim of this study was to investigate the antihypertensive effect of enzymatic hydrolysis of blue mussel protein (HBMP) in rats. Peptides with in vitro ACE inhibitory activity were purified from HBMP by ultrafiltration, gel filtration chromatography and reversed-phase high performance liquid chromatography. And the amino acid sequences of isolated peptides were estimated to be Val-Trp, Leu-Gly-Trp, and Met-Val-Trp-Thr. To study its in vivo action, spontaneously hypertensive rats (SHRs) were orally administration with high- or low-dose of HBMP for 28 days. Major components of the renin-angiotensin (RAS) system in serum of SHRs from different groups were analyzed, and gene expression profiling were performed in the kidney of SHRs, using the Whole Rat Genome Oligonucleotide Microarray. Results indicated although genes involved in RAS system were not significantly altered, those related to blood coagulation system, cytokine and growth factor, and fatty acids metabolism were remarkablely changed. Several genes which were seldom reported to be implicated in pathogenesis of hypertension also showed significant expression alterations after oral administration of HBMP. These data provided valuable information for our understanding of the molecular mechanisms that underlie the potential antihypertensive activities of HBMP, and will contribute towards increased value-added utilization of blue mussel protein.

Development and regulatory application of microRNA biomarkers

MicroRNAs, a class of regulatory small non-coding RNAs, are emerging as promising biomarkers for different health outcomes. Due to their tissue specificity, stability in extracellular space and high conservation between preclinical test species, applications of novel miRNA-based biomarkers for drug safety testing regarding hepatotoxicity and cardiotoxicity are investigated. Furthermore, miRNA expression is altered by environmental exposure such as cigarette smoke or polychlorinated biphenyls. As a consequence, miRNAs potentially influence tumor suppressor genes and oncogenes and may influence carcinogenesis. This has raised the interest in the use of miRNA profiles for health risk assessment. This review summarizes the recent developments in miRNA research with focus on biomarkers for drug safety testing and biomarkers for health outcomes related to environmental exposures.