Gene Expression Profiles from Disease Discordant Twins Suggest Shared Antiviral Pathways and Viral Exposures among Multiple Systemic Autoimmune Diseases

Viral agents are of interest as possible autoimmune triggers due to prior reported associations and widely studied molecular mechanisms of antiviral immune responses in autoimmunity. Here we examined new viral candidates for the initiation and/or promotion of systemic autoimmune diseases (SAID), as well as possible related signaling pathways shared in the pathogenesis of those disorders. RNA isolated from peripheral blood samples from 33 twins discordant for SAID and 33 matched, unrelated healthy controls was analyzed using a custom viral-human gene microarray. Paired comparisons were made among three study groups—probands with SAID, their unaffected twins, and matched, unrelated healthy controls—using statistical and molecular pathway analyses. Probands and unaffected twins differed significantly in the expression of 537 human genes, and 107 of those were associated with viral infections. These 537 differentially expressed human genes participate in overlapping networks of several canonical, biologic pathways relating to antiviral responses and inflammation. Moreover, certain viral genes were expressed at higher levels in probands compared to either unaffected twins or unrelated, healthy controls. Interestingly, viral gene expression levels in unaffected twins appeared intermediate between those of probands and the matched, unrelated healthy controls. Of the viruses with overexpressed viral genes, herpes simplex virus-2 (HSV-2) was the only human viral pathogen identified using four distinct oligonucleotide probes corresponding to three HSV-2 genes associated with different stages of viral infection. Although the effects from immunosuppressive therapy on viral gene expression remain unclear, this exploratory study suggests a new approach to evaluate shared viral agents and antiviral immune responses that may be involved in the development of SAID.     

Adipose tissue gene expression analysis reveals changes in inflammatory, mitochondrial respiratory and lipid metabolic pathways in obese insulin-resistant subjects

Background

To elucidate gene expression associated with copy number changes, we performed a genome-wide copy number and expression microarray analysis of 25 pairs of gastric tissues.

Methods

We applied laser capture microdissection (LCM) to obtain samples for microarray experiments and profiled DNA copy number and gene expression using 244K CGH Microarray and Human Exon 1.0 ST Microarray.

Results

Obviously, gain at 8q was detected at the highest frequency (70%) and 20q at the second (63%). We also identified molecular genetic divergences for different TNM-stages or histological subtypes of gastric cancers. Interestingly, the C20orf11 amplification and gain at 20q13.33 almost separated moderately differentiated (MD) gastric cancers from poorly differentiated (PD) type. A set of 163 genes showing the correlations between gene copy number and expression was selected and the identified genes were able to discriminate matched adjacent noncancerous samples from gastric cancer samples in an unsupervised two-way hierarchical clustering. Quantitative RT-PCR analysis for 4 genes (C20orf11, XPO5, PUF60, and PLOD3) of the 163 genes validated the microarray results. Notably, some candidate genes (MCM4 and YWHAZ) and its adjacent genes such as PRKDC, UBE2V2, ANKRD46, ZNF706, and GRHL2, were concordantly deregulated by genomic aberrations.

Conclusions

Taken together, our results reveal diverse chromosomal region alterations for different TNM-stages or histological subtypes of gastric cancers, which is helpful in researching clinicopathological classification, and highlight several interesting genes as potential biomarkers for gastric cancer.     

A metadata approach for clinical data management in translational genomics studies in breast cancer

Background

Paragangliomas of the head and neck are highly vascular and usually clinically benign tumors arising in the paraganglia of the autonomic nervous system. A significant number of cases (10–50%) are proven to be familial. Multiple genes encoding subunits of the mitochondrial succinate-dehydrogenase (SDH) complex are associated with hereditary paraganglioma: SDHB, SDHC and SDHD. Furthermore, a hereditary paraganglioma family has been identified with linkage to the PGL2 locus on 11q13. No SDH genes are known to be located in the 11q13 region, and the exact gene defect has not yet been identified in this family.

Methods

We have performed a RNA expression microarray study in sporadic, SDHD- and PGL2-linked head and neck paragangliomas in order to identify potential differences in gene expression leading to tumorigenesis in these genetically defined paraganglioma subgroups. We have focused our analysis on pathways and functional gene-groups that are known to be associated with SDH function and paraganglioma tumorigenesis, i.e. metabolism, hypoxia, and angiogenesis related pathways. We also evaluated gene clusters of interest on chromosome 11 (i.e. the PGL2 locus on 11q13 and the imprinted region 11p15).

Results

We found remarkable similarity in overall gene expression profiles of SDHD -linked, PGL2-linked and sporadic paraganglioma. The supervised analysis on pathways implicated in PGL tumor formation also did not reveal significant differences in gene expression between these paraganglioma subgroups. Moreover, we were not able to detect differences in gene-expression of chromosome 11 regions of interest (i.e. 11q23, 11q13, 11p15).

Conclusion

The similarity in gene-expression profiles suggests that PGL2, like SDHD, is involved in the functionality of the SDH complex, and that tumor formation in these subgroups involves the same pathways as in SDH linked paragangliomas. We were not able to clarify the exact identity of PGL2 on 11q13. The lack of differential gene-expression of chromosome 11 genes might indicate that chromosome 11 loss, as demonstrated in SDHD-linked paragangliomas, is an important feature in the formation of paragangliomas regardless of their genetic background.     

Comparative salivary gland transcriptomics of sandfly vectors of visceral leishmaniasis

Background

The multi-step process of carcinogenesis can be more fully understood by characterizing gene expression changes induced in cells by carcinogens. In this study, expression microarrays were used to monitor the activity of 18,224 cDNA clones in MCF-7 and HepG2 cells exposed to the carcinogen benzo(a)pyrene (BaP) or its non-carcinogenic isomer benzo(e)pyrene (BeP). Time and concentration gene expression effects of BaP exposure have been assessed and linked to other measures of cellular stress to aid in the identification of novel genes/pathways involved in the cellular response to genotoxic carcinogens.

Results

BaP (0.25–5.0 μM; 6–48 h exposure) modulated 202 clones in MCF-7 cells and 127 in HepG2 cells, including 27 that were altered in both. In contrast, BeP did not induce consistent gene expression changes at the same concentrations. Significant time- and concentration-dependent responses to BaP were seen in both cell lines. Expression changes observed in both cell lines included genes involved in xenobiotic metabolism (e.g., CYP1B1, NQO1, MGST1, AKR1C1, AKR1C3,CPM), cell cycle regulation (e.g., CDKN1A), apoptosis/anti-apoptosis (e.g., BAX, IER3), chromatin assembly (e.g., histone genes), and oxidative stress response (e.g., TXNRD1). RTqPCR was used to validate microarray data. Phenotypic anchoring of the expression data to DNA adduct levels detected by ³²P-postlabelling, cell cycle data and p53 protein expression identified a number of genes that are linked to these biological outcomes, thereby strengthening the identification of target genes. The overall response to BaP consisted of up-regulation of tumour suppressor genes and down-regulation of oncogenes promoting cell cycle arrest and apoptosis. Anti-apoptotic signalling that may increase cell survival and promote tumourigenesis was also evident.

Conclusion

This study has further characterised the gene expression response of human cells after genotoxic insult, induced after exposure to concentrations of BaP that result in minimal cytotoxiCity. We have demonstrated that investigating the time and concentration effect of a carcinogen on gene expression related to other biological end-points gives greater insight into cellular responses to such compounds and strengthens the identification of target genes.     

A genomic survey of the fish parasite Spironucleus salmonicida indicates genomic plastiCity among diplomonads and significant lateral gene transfer in eukaryote genome evolution

Background

Craniosynostosis, the premature fusion of calvarial sutures, is a common craniofacial abnormality. Causative mutations in more than 10 genes have been identified, involving fibroblast growth factor, transforming growth factor beta, and Eph/ephrin signalling pathways. Mutations affect each human calvarial suture (coronal, sagittal, metopic, and lambdoid) differently, suggesting different gene expression patterns exist in each human suture. To better understand the molecular control of human suture morphogenesis we used microarray analysis to identify genes differentially expressed during suture fusion in children with craniosynostosis. Expression differences were also analysed between each unfused suture type, between sutures from syndromic and non-syndromic craniosynostosis patients, and between unfused sutures from individuals with and without craniosynostosis.

Results

We identified genes with increased expression in unfused sutures compared to fusing/fused sutures that may be pivotal to the maintenance of suture patency or in controlling early osteoblast differentiation (i.e. RBP4, GPC3, C1QTNF3, IL11RA, PTN, POSTN). In addition, we have identified genes with increased expression in fusing/fused suture tissue that we suggest could have a role in premature suture fusion (i.e. WIF1, ANXA3, CYFIP2). Proteins of two of these genes, glypican 3 and retinol binding protein 4, were investigated by immunohistochemistry and localised to the suture mesenchyme and osteogenic fronts of developing human calvaria, respectively, suggesting novel roles for these proteins in the maintenance of suture patency or in controlling early osteoblast differentiation. We show that there is limited difference in whole genome expression between sutures isolated from patients with syndromic and non-syndromic craniosynostosis and confirmed this by quantitative RT-PCR. Furthermore, distinct expression profiles for each unfused suture type were noted, with the metopic suture being most disparate. Finally, although calvarial bones are generally thought to grow without a cartilage precursor, we show histologically and by identification of cartilage-specific gene expression that cartilage may be involved in the morphogenesis of lambdoid and posterior sagittal sutures.

Conclusion

This study has provided further insight into the complex signalling network which controls human calvarial suture morphogenesis and craniosynostosis. Identified genes are candidates for targeted therapeutic development and to screen for craniosynostosis-causing mutations.     

Expression profiling of formalin-fixed paraffin-embedded primary breast tumors using cancer-specific and whole genome gene panels on the DASL® platform

Background

The c D NA-mediated A nnealing, extension, S election and L igation (DASL) assay has become a suitable gene expression profiling system for degraded RNA from paraffin-embedded tissue. We examined assay characteristics and the performance of the DASL 502-gene Cancer Panel^v1 (1.5K) and 24,526-gene panel (24K) platforms at differentiating nine human epidermal growth factor receptor 2- positive (HER2+) and 11 HER2-negative (HER2-) paraffin-embedded breast tumors.

Methods

Bland-Altman plots and Spearman correlations evaluated intra/inter-panel agreement of normalized expression values. Unequal-variance t -statistics tested for differences in expression levels between HER2 + and HER2 - tumors. Regulatory network analysis was performed using Metacore (GeneGo Inc., St. Joseph, MI).

Results

Technical replicate correlations ranged between 0.815-0.956 and 0.986-0.997 for the 1.5K and 24K panels, respectively. Inter-panel correlations of expression values for the common 498 genes across the two panels ranged between 0.485-0.573. Inter-panel correlations of expression values of 17 probes with base-pair sequence matches between the 1.5K and 24K panels ranged between 0.652-0.899. In both panels, erythroblastic leukemia viral oncogene homolog 2 (ERBB2) was the most differentially expressed gene between the HER2 + and HER2 - tumors and seven additional genes had p-values < 0.05 and log2 -fold changes > |0.5| in expression between HER2 + and HER2 - tumors: topoisomerase II alpha (TOP2A), cyclin a2 (CCNA2), v-fos fbj murine osteosarcoma viral oncogene homolog (FOS), wingless-type mmtv integration site family, member 5a (WNT5A), growth factor receptor-bound protein 7 (GRB7), cell division cycle 2 (CDC2), and baculoviral iap repeat-containing protein 5 (BIRC5). The top 52 discriminating probes from the 24K panel are enriched with genes belonging to the regulatory networks centered around v-myc avian myelocytomatosis viral oncogene homolog (MYC), tumor protein p53 (TP53), and estrogen receptor α (ESR1). Network analysis with a two-step extension also showed that the eight discriminating genes common to the 1.5K and 24K panels are functionally linked together through MYC, TP53, and ESR1.

Conclusions

The relative RNA abundance obtained from two highly differing density gene panels are correlated with eight common genes differentiating HER2 + and HER2 - breast tumors. Network analyses demonstrated biological consistency between the 1.5K and 24K gene panels.     

Pathway Analysis Based on a Genome-Wide Association Study of Polycystic Ovary Syndrome

Background

Polycystic ovary syndrome (PCOS) is one of the most common endocrine disorders in women of reproductive age, and it is affected by both environmental and genetic factors. Although the genetic component of PCOS is evident, studies aiming to identify susceptibility genes have shown controversial results. This study conducted a pathway-based analysis using a dataset obtained through a genome-wide association study (GWAS) to elucidate the biological pathways that contribute to PCOS susceptibility and the associated genes.

Methods

We used GWAS data on 636,797 autosomal single nucleotide polymorphisms (SNPs) from 1,221 individuals (432 PCOS patients and 789 controls) for analysis. A pathway analysis was conducted using meta-analysis gene-set enrichment of variant associations (MAGENTA). Top-ranking pathways or gene sets associated with PCOS were identified, and significant genes within the pathways were analyzed.

Results

The pathway analysis of the GWAS dataset identified significant pathways related to oocyte meiosis and the regulation of insulin secretion by acetylcholine and free fatty acids (all nominal gene-set enrichment analysis (GSEA) P-values < 0.05). In addition, INS, GNAQ, STXBP1, PLCB3, PLCB2, SMC3 and PLCZ1 were significant genes observed within the biological pathways (all gene P-values < 0.05).

Conclusions

By applying MAGENTA pathway analysis to PCOS GWAS data, we identified significant pathways and candidate genes involved in PCOS. Our findings may provide new leads for understanding the mechanisms underlying the development of PCOS.     

Expression analyses of nuclear receptor genes in breast cancer cell lines exposed to soy phytoestrogens after BRCA2 knockdown by TaqMan Low-Density Array (TLDA)

Background

Most of breast cancers are considered sporadic and modulation of the two major genes BRCA1 and BRCA2 expressions caused by tissue-specific somatic mutations lead to this pathology. The nutritional intake of phytoestrogens seems to reduce the risk of breast cancer and investigation of their potential as anticancer agents has increased. However, the possible mechanisms and signalling pathways of phytoestrogen action in breast cancer prevention remains unknown.

Results

Using Taqman Low Density Array technology, we investigated the BRCA2 loss of function role in sporadic breast cancers and the links existing with soy isoflavones on a panel of nuclear receptor expression. Human breast cell lines (MCF-7, MDA-MB-231, and MCF-10a) were transfected by BRCA2-siRNA and treated with genistein (18.5 μM) or daidzein (78.5 μM) for 72 h. Generating the transitory knockdown of BRCA2 oncosuppressor, we observed different modulations in several nuclear receptor genes such as ER, RAR and RXR, as well as PPARs and VDR according to the studied breast cell line. Additional isoflavone treatments showed different nuclear receptor gene modulation profiles.

Conclusion

Our results seemed to implicate the oncosuppressor BRCA2 and the phytoestrogen pathways in different nuclear gene expressions via an ER-independent manner.     

Ayurvedic herbal medicine and lead poisoning

Background

Myelodysplastic syndrome with isolated chromosome 5q deletion (5q- syndrome) is a clonal stem cell disorder characterized by ineffective hematopoiesis. MicroRNAs (miRNAs) are important regulators of hematopoiesis and their aberrant expression was detected in some clonal hematopoietic disorders. We thus analyzed miRNA expressions in bone marrow CD34+ cells of 5q- syndrome patients. Further, we studied gene expressions of miR-143, miR-145, miR-378 and miR-146a mapped within the 5q deletion.

Results

Using microarrays we identified 21 differently expressed miRNAs in 5q- patients compared to controls. Especially, miR-34a was markedly overexpressed in 5q- patients, suggesting its role in an increased apoptosis of bone marrow progenitors. Out of four miRNAs at del(5q), only miR-378 and miR-146a showed reduced gene expression in the patients. An integrative analysis of mRNA profiles and predicted putative targets defined potential downstream targets of the deregulated miRNAs. The list of targets included several genes that play an important role in the regulation of hematopoiesis (e.g. KLF4, LEF1, SPI1).

Conclusions

The study demonstrates global overexpression of miRNAs is associated with 5q- phenotype. Identification of hematopoiesis-relevant target genes indicates that the deregulated miRNAs may be involved in the pathogenesis of 5q- syndrome by a modulation of these targets. The expression data on miRNAs at del(5q) suggest the presence of mechanisms for compensation of a gene dosage.     

Expression patterns of WRKY genes in di-haploid Populus simonii × P. nigra in response to salinity stress revealed by quantitative real-time PCR and RNA sequencing

Key message

Spatio-temporal expression patterns of 13 out of 119 poplar WRKY genes indicated dynamic and tissue-specific roles of WRKY family proteins in salinity stress tolerance.

Abstract

To understand the expression patterns of poplar WRKY genes under salinity stress, 51 of the 119 WRKY genes were selected from di-haploid Populus simonii × P. nigra by quantitative real-time PCR (qRT-PCR). We used qRT-PCR to profile the expression of the top 13 genes under salinity stress across seven time points, and employed RNA-Seq platforms to cross-validate it. Results demonstrated that all the 13 WRKY genes were expressed in root, stem, and leaf tissues, but their expression levels and overall patterns varied notably in these tissues. Regarding overall gene expression in roots, the 13 genes were significantly highly expressed at all six time points after the treatment, reaching the plateau of expression at hour 9. In leaves, the 13 genes were similarly up-regulated from 3 to 12 h in response to NaCl treatment. In stems, however, expression levels of the 13 genes did not show significant changes after the NaCl treatment. Regarding individual gene expression across the time points and the three tissues, the 13 genes can be classified into three clusters: the lowly expressed Cluster 1 containing PthWRKY28, 45 and 105; intermediately expressed Clusters 2 including PthWRKY56, 88 and 116; and highly expressed Cluster 3 consisting of PthWRKY41, 44, 51, 61, 62, 75 and 106. In general, genes in Cluster 2 and 3 displayed a dynamic pattern of “induced amplification—recovering”, suggesting that these WRKY genes and corresponding pathways may play a critical role in mediating salt response and tolerance in a dynamic and tissue-specific manner.     

Spontaneous inflammatory pain model from a mouse line with N-ethyl-N-nitrosourea mutagenesis

Background

Most filarial nematodes contain Wolbachia symbionts. The purpose of this study was to examine the effects of doxycycline on gene expression in Wolbachia and adult female Brugia malayi.

Methods

Brugia malayi infected gerbils were treated with doxycycline for 6-weeks. This treatment largely cleared Wolbachia and arrested worm reproduction. RNA recovered from treated and control female worms was labeled by random priming and hybridized to the Version 2- filarial microarray to obtain expression profiles.

Results and discussion

Results showed significant changes in expression for 200 Wolbachia (29% of Wolbachia genes with expression signals in untreated worms) and 546 B. malayi array elements after treatment. These elements correspond to known genes and also to novel genes with unknown biological functions. Most differentially expressed Wolbachia genes were down-regulated after treatment (98.5%). In contrast, doxycycline had a mixed effect on B. malayi gene expression with many more genes being significantly up-regulated after treatment (85% of differentially expressed genes). Genes and processes involved in reproduction (gender-regulated genes, collagen, amino acid metabolism, ribosomal processes, and cytoskeleton) were down-regulated after doxycycline while up-regulated genes and pathways suggest adaptations for survival in response to stress (energy metabolism, electron transport, anti-oxidants, nutrient transport, bacterial signaling pathways, and immune evasion).

Conclusions

Doxycycline reduced Wolbachia and significantly decreased bacterial gene expression. Wolbachia ribosomes are believed to be the primary biological target for doxycycline in filarial worms. B. malayi genes essential for reproduction, growth and development were also down-regulated; these changes are consistent with doxycycline effects on embryo development and reproduction. On the other hand, many B. malayi genes involved in energy production, electron-transport, metabolism, anti-oxidants, and others with unknown functions had increased expression signals after doxycycline treatment. These results suggest that female worms are able to compensate in part for the loss of Wolbachia so that they can survive, albeit without reproductive capacity. This study of doxycycline induced changes in gene expression has provided new clues regarding the symbiotic relationship between Wolbachia and B. malayi.