Using genome-wide associations to identify metabolic pathways involved in maize aflatoxin accumulation resistance

Background

Aflatoxin is a potent carcinogen that can contaminate grain infected with the fungus Aspergillus flavus. However, resistance to aflatoxin accumulation in maize is a complex trait with low heritability. Here, two complementary analyses were performed to better understand the mechanisms involved. The first coupled results of a genome-wide association study (GWAS) that accounted for linkage disequilibrium among single nucleotide polymorphisms (SNPs) with gene-set enrichment for a pathway-based approach. The rationale was that the cumulative effects of genes in a pathway would give insight into genetic differences that distinguish resistant from susceptible lines of maize. The second involved finding non-pathway genes close to the most significant SNP-trait associations with the greatest effect on reducing aflatoxin in multiple environments. Unlike conventional GWAS, the latter analysis emphasized multiple aspects of SNP-trait associations rather than just significance and was performed because of the high genotype x environment variability exhibited by this trait.

Results

The most significant metabolic pathway identified was jasmonic acid (JA) biosynthesis. Specifically, there was at least one allelic variant for each step in the JA biosynthesis pathway that conferred an incremental decrease to the level of aflatoxin observed among the inbred lines in the GWAS panel. Several non-pathway genes were also consistently associated with lowered aflatoxin levels. Those with predicted functions related to defense were: leucine-rich repeat protein kinase, expansin B3, reversion-to-ethylene sensitivity1, adaptor protein complex2, and a multidrug and toxic compound extrusion protein.

Conclusions

Our genetic analysis provided strong evidence for several genes that were associated with aflatoxin resistance. Inbred lines that exhibited lower levels of aflatoxin accumulation tended to share similar haplotypes for genes specifically in the pathway of JA biosynthesis, along with several non-pathway genes with putative defense-related functions. Knowledge gained from these two complementary analyses has improved our understanding of population differences in aflatoxin resistance.

Electronic supplementary material

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

Stromelysin-1 and mesothelin are differentially regulated by Wnt-5a and Wnt-1 in C57mg mouse mammary epithelial cells

Background  

The Wnt signal transduction pathway is important in a wide variety of developmental processes as well as in the genesis of human cancer. Vertebrate Wnt pathways can be functionally separated into two classes, the canonical Wnt/beta-catenin pathway and the non-canonical Wnt/Ca²⁺ pathway. Supporting differences in Wnt signaling, gain of function of Wnt-1 in C57mg mouse mammary epithelial cells leads to their morphological transformation while loss of function of Wnt-5a leads to the same transformation. Many downstream target genes of the Wnt/beta-catenin pathway have been identified. In contrast, little is known about the Wnt/Ca²⁺ pathway and whether it regulates gene expression.     

Ror2 modulates the canonical Wnt signaling in lung epithelial cells through cooperation with Fzd2

Background  

Wnt signaling is mediated through 1) the beta-catenin dependent canonical pathway and, 2) the beta-catenin independent pathways. Multiple receptors, including Fzds, Lrps, Ror2 and Ryk, are involved in Wnt signaling. Ror2 is a single-span transmembrane receptor-tyrosine kinase (RTK). The functions of Ror2 in mediating the non-canonical Wnt signaling have been well established. The role of Ror2 in canonical Wnt signaling is not fully understood.     

Integrity of chromatin and replicating DNA in nuclei released from fission yeast by semi-automated grinding in liquid nitrogen

Background

Adenoviruses force quiescent cells to re-enter the cell cycle to replicate their DNA, and for the most part, this is accomplished after they express the E1A protein immediately after infection. In this context, E1A is believed to inactivate cellular proteins (e.g., p130) that are known to be involved in the silencing of E2F-dependent genes that are required for cell cycle entry. However, the potential perturbation of these types of genes by E1A relative to their functions in regulatory networks and canonical pathways remains poorly understood.

Findings

We have used DNA microarrays analyzed with Bayesian ANOVA for microarray (BAM) to assess changes in gene expression after E1A alone was introduced into quiescent cells from a regulated promoter. Approximately 2,401 genes were significantly modulated by E1A, and of these, 385 and 1033 met the criteria for generating networks and functional and canonical pathway analysis respectively, as determined by using Ingenuity Pathway Analysis software. After focusing on the highest-ranking cellular processes and regulatory networks that were responsive to E1A in quiescent cells, we observed that many of the up-regulated genes were associated with DNA replication, the cell cycle and cellular compromise. We also identified a cadre of up regulated genes with no previous connection to E1A; including genes that encode components of global DNA repair systems and DNA damage checkpoints. Among the down-regulated genes, we found that many were involved in cell signalling, cell movement, and cellular proliferation. Remarkably, a subset of these was also associated with p53-independent apoptosis, and the putative suppression of this pathway may be necessary in the viral life cycle until sufficient progeny have been produced.

Conclusions

These studies have identified for the first time a large number of genes that are relevant to E1A's activities in promoting quiescent cells to re-enter the cell cycle in order to create an optimum environment for adenoviral replication.     

Ancient intron insertion sites and palindromic genomic duplication evolutionally shapes an elementally functioning membrane protein family

Background  

In spite of the recent accumulation of genomic data, the evolutionary pathway in the individual genes of present-day living taxa is still elusive for most genes. Among ion channels, inward K⁺ rectifier (IRK) channels are the fundamental and well-defined protein group. We analyzed the genomic structures of this group and compared them among a phylogenetically wide range with our sequenced Halocynthia roretzi, a tunicate, IRK genomic genes.     

Extending pathways and processes using molecular interaction networks to analyse cancer genome data

Background  

Cellular processes and pathways, whose deregulation may contribute to the development of cancers, are often represented as cascades of proteins transmitting a signal from the cell surface to the nucleus. However, recent functional genomic experiments have identified thousands of interactions for the signalling canonical proteins, challenging the traditional view of pathways as independent functional entities. Combining information from pathway databases and interaction networks obtained from functional genomic experiments is therefore a promising strategy to obtain more robust pathway and process representations, facilitating the study of cancer-related pathways.     

Estimation and testing for the effect of a genetic pathway on a disease outcome using logistic kernel machine regression via logistic mixed models

Background  

Growing interest on biological pathways has called for new statistical methods for modeling and testing a genetic pathway effect on a health outcome. The fact that genes within a pathway tend to interact with each other and relate to the outcome in a complicated way makes nonparametric methods more desirable. The kernel machine method provides a convenient, powerful and unified method for multi-dimensional parametric and nonparametric modeling of the pathway effect.     

The androgen receptor can signal through Wnt/β-Catenin in prostate cancer cells as an adaptation mechanism to castration levels of androgens

Background  

A crucial event in Prostate Cancer progression is the conversion from a hormone-sensitive to a hormone-refractory disease state. Correlating with this transition, androgen receptor (AR) amplification and mutations are often observed in patients failing hormonal ablation therapies. β-Catenin, an essential component of the canonical Wnt signaling pathway, was shown to be a coactivator of the AR signaling in the presence of androgens. However, it is not yet clear what effect the increased levels of the AR could have on the Wnt signaling pathway in these hormone-refractory prostate cells.     

Male-like sexual behavior of female mouse lacking fucose mutarotase

Background

Quantitative traits often underlie risk for complex diseases. For example, weight and body mass index (BMI) underlie the human abdominal obesity-metabolic syndrome. Many attempts have been made to identify quantitative trait loci (QTL) over the past decade, including association studies. However, a single QTL is often capable of affecting multiple traits, a quality known as gene pleiotropy. Gene pleiotropy may therefore cause a loss of power in association studies focused only on a single trait, whether based on single or multiple markers.

Results

We propose using principal-component-based multivariate regression (PCBMR) to test for gene pleiotropy with comprehensive evaluation. This method generates one or more independent canonical variables based on the principal components of original traits and conducts a multivariate regression to test for association with these new variables. Systematic simulation studies have shown that PCBMR has great power. PCBMR-based pleiotropic association studies of abdominal obesity-metabolic syndrome and its possible linkage to chromosomal band 3q27 identified 11 susceptibility genes with significant associations. Whereas some of these genes had been previously reported to be associated with metabolic traits, others had never been identified as metabolism-associated genes.

Conclusions

PCBMR is a computationally efficient and powerful test for gene pleiotropy. Application of PCBMR to abdominal obesity-metabolic syndrome indicated the existence of gene pleiotropy affecting this syndrome.     

Comprehensive Expression of Wnt Signaling Pathway Genes during Development and Maturation of the Mouse Cochlea

Background

In the inner ear Wnt signaling is necessary for proliferation, cell fate determination, growth of the cochlear duct, polarized orientation of stereociliary bundles, differentiation of the periotic mesenchyme, and homeostasis of the stria vascularis. In neonatal tissue Wnt signaling can drive proliferation of cells in the sensory region, suggesting that Wnt signaling could be used to regenerate the sensory epithelium in the damaged adult inner ear. Manipulation of Wnt signaling for regeneration will require an understanding of the dynamics of Wnt pathway gene expression in the ear. We present a comprehensive screen for 84 Wnt signaling related genes across four developmental and postnatal time points.

Results

We identified 72 Wnt related genes expressed in the inner ear on embryonic day (E) 12.5, postnatal day (P) 0, P6 and P30. These genes included secreted Wnts, Wnt antagonists, intracellular components of canonical signaling and components of non-canonical signaling/planar cell polarity.

Conclusion

A large number of Wnt signaling molecules were dynamically expressed during cochlear development and in the early postnatal period, suggesting complex regulation of Wnt transduction. The data revealed several potential key regulators for further study.     

Canonical Wnt signaling differently modulates osteogenic differentiation of mesenchymal stem cells derived from bone marrow and from periodontal ligament under inflammatory conditions

Background

Cellular plasticity and complex functional requirements of the periodontal ligament (PDL) assume a local stem cell (SC) niche to maintain tissue homeostasis and repair. Here, pathological alterations caused by inflammatory insults might impact the regenerative capacities of these cells. As bone homeostasis is fundamentally controlled by Wnt-mediated signals, it was the aim of this study to characterize the SC-like capacities of cells derived from PDL and to investigate their involvement in bone pathophysiology especially regarding the canonical Wnt pathway.

Methods

PDLSCs were investigated for their SC characteristics via analysis of cell surface marker expression, colony forming unit efficiency, proliferation, osteogenic differentiation and adipogenic differentiation, and compared to bone marrow derived mesenchymal SCs (BMMSCs). To determine the impact of both inflammation and the canonical Wnt pathway on osteogenic differentiation, cells were challenged with TNF-α, maintained with or without Wnt3a or DKK-1 under osteogenic induction conditions and investigated for p-IκBα, p-NF-κB, p-Akt, β-catenin, p-GSK-3β, ALP and Runx2.

Results

PDLSCs exhibit weaker adipogenic and osteogenic differentiation capacities compared to BMMSCs. TNF-α inhibited osteogenic differentiation of PDLSCs more than BMMSCs mainly through regulating canonical Wnt pathway. Blocking the canonical Wnt pathway by DKK-1 reconstituted osteogenic differentiation of PDLSCs under inflammatory conditions, whereas activation by Wnt3a increased osteogenic differentiation of BMMSCs.

Conclusions

Our results suggest a diverse regulation of the inhibitory effect of TNF-α in BMMSCs and PDLSCs via canonical Wnt pathway modulation.

General significance

These findings provide novel insights on PDLSC SC-like capacities and their involvement in bone pathophysiology under the impact of the canonical Wnt pathway.     

Identifying differential correlation in gene/pathway combinations

Background  

An important emerging trend in the analysis of microarray data is to incorporate known pathway information a priori. Expression level "summaries" for pathways, obtained from the expression data for the genes constituting the pathway, permit the inclusion of pathway information, reduce the high dimensionality of microarray data, and have the power to elucidate gene-interaction dependencies which are not already accounted for through known pathway identification.     

Components of the Canonical and Non-Canonical Wnt Pathways Are Not Mis-Expressed in Pituitary Tumors

Introduction

Canonical and non-canonical Wnt pathways are involved in the genesis of multiple tumors; however, their role in pituitary tumorigenesis is mostly unknown.

Objective

This study evaluated gene and protein expression of Wnt pathways in pituitary tumors and whether these expression correlate to clinical outcome.

Materials and Methods

Genes of the Wnt canonical pathway: activating ligands (WNT11, WNT4, WNT5A), binding inhibitors (DKK3, sFRP1), β-catenin (CTNNB1), β-catenin degradation complex (APC, AXIN1, GSK3β), inhibitor of β-catenin degradation complex (AKT1), sequester of β-catenin (CDH1), pathway effectors (TCF7, MAPK8, NFAT5), pathway mediators (DVL-1, DVL-2, DVL-3, PRICKLE, VANGL1), target genes (MYB, MYC, WISP2, SPRY1, TP53, CCND1); calcium dependent pathway (PLCB1, CAMK2A, PRKCA, CHP); and planar cell polarity pathway (PTK7, DAAM1, RHOA) were evaluated by QPCR, in 19 GH-, 18 ACTH-secreting, 21 non-secreting (NS) pituitary tumors, and 5 normal pituitaries. Also, the main effectors of canonical (β-catenin), planar cell polarity (JNK), and calcium dependent (NFAT5) Wnt pathways were evaluated by immunohistochemistry.

Results

There are no differences in gene expression of canonical and non-canonical Wnt pathways between all studied subtypes of pituitary tumors and normal pituitaries, except for WISP2, which was over-expressed in ACTH-secreting tumors compared to normal pituitaries (4.8x; p = 0.02), NS pituitary tumors (7.7x; p = 0.004) and GH-secreting tumors (5.0x; p = 0.05). β-catenin, NFAT5 and JNK proteins showed no expression in normal pituitaries and in any of the pituitary tumor subtypes. Furthermore, no association of the studied gene or protein expression was observed with tumor size, recurrence, and progressive disease. The hierarchical clustering showed a regular pattern of genes of the canonical and non-canonical Wnt pathways randomly distributed throughout the dendrogram.

Conclusions

Our data reinforce previous reports suggesting no activation of canonical Wnt pathway in pituitary tumorigenesis. Moreover, we describe, for the first time, evidence that non-canonical Wnt pathways are also not mis-expressed in the pituitary tumors.