Next Article in Journal
Horticultural Therapy May Reduce Psychological and Physiological Stress in Adolescents with Anorexia Nervosa: A Pilot Study
Next Article in Special Issue
A Pleiotropic Nuclear Hormone Labelled Hundred Years Ago Vitamin D
Previous Article in Journal
Sleep Deprivation and Central Appetite Regulation
Previous Article in Special Issue
Interaction of Vitamin D with Peptide Hormones with Emphasis on Parathyroid Hormone, FGF23, and the Renin-Angiotensin-Aldosterone System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Vitamin D and Systems Biology

by
Shahid Hussain
1,
Clayton Yates
2,3,4 and
Moray J. Campbell
1,*,†
1
Division of Pharmaceutics and Pharmaceutical Chemistry, College of Pharmacy, The Ohio State University, Columbus, OH 43210, USA
2
Department of Biology and Center for Cancer Research, Tuskegee University, Tuskegee, AL 36088, USA
3
Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD 21287, USA
4
Department of Oncology Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD 21287, USA
*
Author to whom correspondence should be addressed.
Current address: Department of Biomedical Sciences, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA.
Nutrients 2022, 14(24), 5197; https://doi.org/10.3390/nu14245197
Submission received: 26 September 2022 / Revised: 28 November 2022 / Accepted: 1 December 2022 / Published: 7 December 2022

Abstract

:
The biological actions of the vitamin D receptor (VDR) have been investigated intensively for over 100 years and has led to the identification of significant insights into the repertoire of its biological actions. These were initially established to be centered on the regulation of calcium transport in the colon and deposition in bone. Beyond these well-known calcemic roles, other roles have emerged in the regulation of cell differentiation processes and have an impact on metabolism. The purpose of the current review is to consider where applying systems biology (SB) approaches may begin to generate a more precise understanding of where the VDR is, and is not, biologically impactful. Two SB approaches have been developed and begun to reveal insight into VDR biological functions. In a top-down SB approach genome-wide scale data are statistically analyzed, and from which a role for the VDR emerges in terms of being a hub in a biological network. Such approaches have confirmed significant roles, for example, in myeloid differentiation and the control of inflammation and innate immunity. In a bottom-up SB approach, current biological understanding is built into a kinetic model which is then applied to existing biological data to explain the function and identify unknown behavior. To date, this has not been applied to the VDR, but has to the related ERα and identified previously unknown mechanisms of control. One arena where applying top-down and bottom-up SB approaches may be informative is in the setting of prostate cancer health disparities.

1. Systems Biology and Biomedicine

1.1. The Opportunities of Applying Systems Biology Approaches in Biomedical Research

Much of the discovery in biomedicine has been centered on the classic paradigms of reductionist biology, in which phenotypes are interpreted as the interaction of either single or small groups of molecules. Across biomedical sciences, there are many remarkable examples of how this approach has led to the discovery of drivers of human disease, and equally remarkable examples of new therapies designed to target these drivers. For example, high-profile examples of this reductionist approach include the discovery of the oncogenic role of the BCR-ABL fusion gene in the etiology of chronic myeloid leukemia [1] and the development of a targeting kinase inhibitor such as Imatinib [2].
This is a striking example of the so-called bench-to-bedside research and represents one of the earliest examples of precision medicine. At times, however, this reductionist approach can appear limited both theoretically and clinically as the etiology of Imatinib-resistant phenotypes only too well demonstrates [3]. Across cancers and other disease phenotypes, a stumbling block can be identifying single strong disease driver mechanisms, which in turn accurately predict drug sensitivities and therapeutic effectiveness. Therefore, delivering a fuller prediction of disease drivers and therapeutic vulnerabilities may require developing different methodologies. Ideally, any methodologies would also take advantage of the ever-increasing stream of high-dimensional biological data to inform diagnosis and prognosis. Perhaps these approaches (SB and reductionist) can actually be highly symbiotic.
Systems biology (SB) aims to apply mathematical approaches to build a predicative and quantitative model of biological systems, with the goal to use model predictions to define specific physiological or pathophysiological states and outcomes. The models derived from such SB approaches applied to experimentally derived biological data aim to identify the dynamic behavior of networks that are the center of cellular behaviors [4,5]. These approaches are readily scalable and not restricted by scope. SB approaches can be applied to discrete cell signaling systems, such as gene regulatory networks and signal transduction cascades, to cell–cell interactions, tissue organization, organismal behavior, and to complex multi-organism interactions as seen within, for example, the function of the gut and even complete ecosystems [6].
Furthermore, the models built by SB approaches aim to define the functioning of living organisms not solely by looking at the constituent molecules such as DNA, RNA, proteins, and metabolites but rather by the process-level biological systems they constitute, such as mitosis and metabolic control [7,8]. A key hallmark of these models is that they capture the states of the system in a predicative and quantitative manner and can be exploited to drive novel understanding. As these in silico models can be interrogated rapidly, multiple components of any given model can be dissected to reveal unintuitive findings [9,10,11].
A consequence of the integrated nature of biological signaling is the emergent complexity, which underpins human health. For example, the dexterous control of transcriptional networks is derived from a high ratio of transcription factors to regulated mRNA or miRNA targets. This ratio, combined with large regulatory regions, results in unparalleled plasticity over the choice, amplitude, and period of transcription [12,13]. More specifically, transcription factor modules recognize, interpret and sustain histone modifications and ultimately establish boundaries between transcriptionally rich euchromatin and transcriptionally restricted heterochromatin. These boundaries are cemented further by the regulation of CpG island methylation; these two processes are dynamically intertwined. Again, from an SB perspective, it is reasoned that modeling these transcriptional processes will help to explain the highly integrated nature and robustness of normal transcriptional control, whereby the processes have redundancies such that they do not radically alter in response to external signals and do not fail when a single component is altered. By contrast, transcriptional networks in cancer cells display a loss of transcriptional plasticity and do not display the full breadth of signaling capacities. The evolution of the malignant transcriptome is seen clearly in the nuclear receptor and MYC superfamilies (reviewed in [14]).
The concept of SB builds on so-called “holistic biology” developed in the 1960s and coupled with informatics theory and modeling approaches developed through the 20th century. However, the application of SB methodologies and the expansion of SB concepts in biomedicine has been boosted since the early 2000s by the technological advances in the development of high dimensional data approaches frequently derived from next-generation sequencing technologies coupled with bioinformatic approaches. This progression from the 1950s to the current state has been profoundly catalyzed by the human genome project, with its draft sequence published in 2001, and a final reference genome published in 2022 [15] alongside other reference genomes across the animal and plant kingdom. In this manner, the combination of high throughput experimental approaches in the wet lab coupled with complex statistical analyses and computational methods in the dry lab have led to a more comprehensive characterization of multiple diverse organisms, and a shift of focus from molecules to their interactions and the networks they form [16]. This sophistication and power of prediction are most likely set to increase when SB models also include a spatio-temporal characterization of cell behavior, including the dynamics of how molecules are exchanged between compartments and exported from the cell. Again, modeling these aspects of cell behavior has been massively impacted by an explosion in single-cell and spatial technologies.
There are several well-justified advantages of applying SB approaches, which have the genuine potential to complement and extend the reductionist paradigm. Expressing the complete interactions within biological systems in mathematical terms reduces the impact of biases introduced by focusing on single proteins, and instead can reveal under-explored control points in the system. Furthermore, expressing biological systems in the context of processes, rather than well-understood individual components, has the potential to assimilate more readily new high dimensional data generated in biological experiments. Consequently, it is reasonably anticipated that significant, novel, and unpredicted strides will be made in understanding the control and responsiveness of biological events, which in turn can generate new insights into disease susceptibility and therapeutic opportunities [17,18].

1.2. Systems Biology Builds upon an Asymptotic Recursion between Wet Lab and Dry Lab

At its core, SB approaches have an asymptotic recursion between research activities in the wet lab and dry lab. Models are built in the dry lab that reflects the biological observations in the wet lab and are used to generate predictions for how the system can behave. Such predictions are then formulated into testable interventions in the wet lab to generate data for model refinement in the dry lab. It is this oscillation between wet and dry labs that is so potentially powerful, but at the same time so daunting. Broadly, two strategies have emerged to develop such models, by either top-down or bottom-up approaches (Figure 1).
In the top-down approach, large datasets are interrogated with statistical methods to find patterns in the data with which to derive predictions on the system organization [19]. The models in top-down SB are phenomenological, meaning they are not directly mechanistically based and do not require knowledge about relationships between different molecular components. Identification of significant associations is used to develop a hypothesis on the nature of the molecule interactions and associations identified. These approaches naturally lend themselves to omics-derived data and develop hypotheses to be tested by wet lab analysis [20]. It is worth noting that the identified associations that are significant may not be causal, and in fact, themselves may not be true. Top-down approaches are often used with subsystems that have not yet been characterized to a high level of mechanistic detail and approaches such as Bayesian modeling are appropriate due to the missingness in the data from biological regulatory networks [21].
By contrast, the bottom-up approach begins with the hypothesis of biological mechanism and formulating equations on how the components of the system interact and then running simulations to generate predictions. This approach relies on experimental studies to determine the kinetic and chemical properties of the components, and starts with formulating system behaviors in rate equations, for example, expressed as differential equations, of the constituting parts of each system. These formulations are then integrated to predict the system behavior, with the goal to combine pathway models into a model for a larger system [21]. The data on the system under study are subjected to perturbations in the cell context and models are refined from the data. An example of a bottom-up approach is the silicon cell program where computational replicas of actual pathways are made to calculate system behavior [22].
The power of these approaches arises from being able to identify and direct experiments to test fundamental questions of a biological system. From the top-down perspective, these questions include identifying in a genome-wide manner the interactions of all components in a system to define the metabolic control that ultimately brings about cell, tissue, or organism behavior. Arising from this approach questions can be asked of a system. For example, what is the interconnectedness of the system and how does that change between health and disease states, or in different development or differentiation states? Which hubs in such networks are central and which are peripheral? Again, how does hub distribution shift in disease and development transitions? How do changes in gene and protein expression combine to control metabolism? How do germline or somatic structural variants change network topology and metabolic flux?
Similarly, from the bottom-up perspective, a model is curated from known biological interactions and expressed in mathematical terms to test behavior and make new predictions. For example, with any signaling system how is activation controlled and silenced? Given the ubiquitous role of enzymes in biology, a common question is how does changing the kinetics of activating/de-activating enzymes lead to altered signal amplitude? To what extent do germline or somatic structural variants in enzymes, or co-factors modulate enzymatic capacity and what effect does that exert on signal strength? More broadly, how are different signaling systems integrated and what is the quantitative effect of convergence? Does crosstalk generate synergistic events in regulatory complexes or are they merely convergent downstream at endpoints? What is the magnitude of control at each step and is there an order to control ranging from find-tuning to signal-independent activation? Finally, questions can be asked of the system in disease or development states and tested for pharmacological relevance including identifying which targets may be impactful but have deleterious side effects.

2. The Opportunities and Challenges of Applying Systems Biology Approaches to Studying the Vitamin D Receptor

A systems-level appreciation for vitamin D signaling has existed for several centuries, given that a description of rickets was first described in the 17th century [23]; rickets arises from impaired bone mineralization due to insufficient signaling via the vitamin D receptor (NR1I1/VDR) (reviewed in [24]).
The VDR is a Type II member of the nuclear receptor superfamily, which binds the active hormone 1α,25(OH)2D3. Several features of how the VDR functions are important from an SB perspective. A feature of Type II receptors such as the VDR is that independent of ligand the receptor is significantly associated with the genome bound to cis-regulatory elements (CRE) and may exert repressive effects, which is reversed by ligand activation leading to genomic redistribution and transactivation [25,26,27,28]. Additionally, of interest from a systems perspective is that VDR interacts in distinct ways with a number of proteins. Firstly, it heterodimerizes with the RXRs (NR2B1/RXRα, NR2B2/RXRβ, NR2B3/RXRγ), which is also a central dimer partner for other Type II receptors including those retinoic acid receptors (NR1B1/RARα, NR1B2/RARβ, NR1B3/RARγ), and the peroxisome proliferator-activated receptors (NR1C1/PPARα, NR1C2/PPARβ, NR1C3/PPARγ) [reviewed in [24]]. Secondly, the VDR interacts with a range of coregulator proteins such as coactivators and corepressors that exert antagonist roles in the control of local chromatin structure, as well as components of the SWI/SNF complexes to remodel nucleosome positioning and other enzymes such as helicases [29,30] and splicing factors [31] to also facilitate transcription. The VDR appears to participate in protein–protein–DNA interactions, for example with other transcription factors, in a trans-regulatory mechanism to regulate transcription. Finally, there is emerging evidence for the VDR, like other nuclear receptors [32,33], to interact functionally with long non-coding RNAs (lncRNA) to control gene expression by directly affecting the DNA environment or other RNA binding proteins [34]. Indeed one such lncRNA, steroid receptor activator (SRA) coimmunoprecipitates with the coactivator NCOA1/SRC1 and may function more broadly as a scaffold for NR complexes (reviewed in [32]).
Finally, of central importance to the complete VDR signaling system is the generation of the ligand, 1α,25(OH)2D3, the levels of which are dynamically controlled in terms of synthesis of the precursor 25(OH)D3 which is held in a relatively tight range in the serum, and signs of deficiency in the VDR system can be seen when this serum level is diminished. The various metabolic steps that lead to the generation of 25(OH)D3, the active ligand 1α,25(OH)2D3, and a range of metabolites that ultimately lead to its catabolism are all tightly controlled by enzymatic reactions both in an endocrine and tissue-specific local intracrine manner.

2.1. Top-Down Approaches Applied to VDR Biology

From a top-down SB perspective, it is interesting to ask in what context the VDR itself, or coregulators, or ligand generation are identified as a central hub in gene networks associated with different cell phenotypes and disease states, and from a bottom-up SB perspective, it is attractive to develop models that accurately capture the cycling interactions of the VDR with different co-factor and predict how this relates to diverse transcriptional outputs, again across cell phenotypes and disease states.
High dimensional data approaches have been applied to ask questions about what the transcriptional effects are of adding 1α,25(OH)2D3, or deleting the VDR. From a top-down SB perspective, the question of VDR function is alternatively phrased to identify biological circumstances where the VDR is identified as a hub in the transcriptional network independent of any knowledge of the system. In this manner, analyses of across myeloid cells [35], granulocytes [36], and megakaryocytes [37] have identified significant control functions for the VDR to regulate specific cell differentiation outcomes. Interestingly, these studies support some of the earliest studies that explored the functions of 1α,25(OH)2D3 and revealed its capacity to initiate differentiation of leukemia cells [38,39]. These findings might suggest a role for the VDR to be mutated or deleted in leukemia and although there were numerous candidate studies of the VDR expression and genomic integrity, large-scale genomic approaches, for example in the TCGA leukemia cohort [40], suggest that the VDR is neither distorted to act as cancer-driver, nor has it been therapeutically exploited to date.
Together, these findings suggest that whilst the VDR is biologically impactful in the normal differentiation, for example, the process to monocytes, it is not itself so frequently disrupted in leukemia such that when it is mutated it acts in an oncogenic manner to disrupt normal progenitor differentiation, and it cannot be targeted to induce leukemia cell differentiation in patients. At first glance, these findings may appear contradictory, but most likely reflect the nature of redundancy in the system, and the role of VDR to undergo diverse protein interactions. Specifically, the VDR is known to interact with CEBPs [41] to induce leukemia cell differentiation [41], and indeed the CEBPs are a master regulator in the physiological transcriptional module containing the VDR to regulate monocyte differentiation [35]. However, CEPPs interact with multiple different nuclear receptors to regulate cell fates [42,43,44,45], and redundancy in the system limits the impact that altered VDR function alone can exert on cell fates. Finally, it is interesting to note that CEBPA mutations are in the top 33 mutational events in the TCGA leukemia cohort [40].
These identified roles for the VDR in monocyte differentiation reflect a broader function for the receptor to signal in differentiation processes [46] and immune phenotypes. Alongside the roles of vitamin D to prevent rickets, in the late 19th century the concept emerged of sunlight therapy, known as heliotherapy, to fight tuberculosis and other infections [47]. Indeed, unbiased transcriptomic analyses of the impact of Mycobacterium tuberculosis demonstrated a role for Toll-like receptors to initiate a signal transduction cascade in macrophages that upregulated the VDR and induced an antimicrobial response. Intriguingly, this response appeared to be dampened in people with African genomic ancestry perhaps associated with low serum levels of 25(OH)D3 [48]. Similarly, Mycobacterium tuberculosis infection of lung cells identified VDR as a key regulated transcription factor [49]. Another infection of interest is the coronaviruses and in 2013 the lung cell response towards severe acute respiratory syndrome (SARS) arising from a coronavirus infection again identified the VDR as a key regulator of response [50]. More recently, COVID-19 infection of T cells from bronchoalveolar lavage demonstrated a role for the VDR to trigger super-enhancer activation and the control of innate immunity [51]. An unbiased identification of the role of the VDR in the regulation of innate immunity and immune phenotypes has also been identified by GWAS [52,53,54,55].
Viewed from the perspective of unbiased top-down statistical analyses has revealed critical roles for the VDR in innate immunity, control of inflammation, regulation of differentiation, and metabolic control [56,57,58], but not in some of the other phenotypes that are investigated at the pre-clinical and candidate level, notably including cancer; this was established in a previous structured literature search (reviewed in [59]). That is, no cancer-associated GWAS identifies structural variants in the VDR associated with cancer risk, and none of the TCGA papers (over 10,000 tumor samples across more than 30 tumor types) identify frequent VDR-associated structural variants or significantly altered expression that associates with cancer clinical phenotypes.
GWAS studies have continued to evolve in terms of the complexity of study design [60] and cohort sizes and different genomic ancestry, and methods have also developed to finding interactions between variants and also considering how non-coding variants withing regulatory regions may determine how efficiently a given transcription factor can function; so-called a cis-expression quantitative trait loci (eQTL). Certainly, there is evidence for significant enrichment of germline structural variants in VDR binding sites associated with genes that regulate immune phenotypes [61]. Likewise, there is evidence for such eQTLs are significantly impacted by genomic ancestry, for example in the case of the ex vivo response towards 1α,25(OH)2D3 in normal colon cells [62].
These methodological developments and the ever-increasing repertoire of high dimensional data are likely to increase the understanding of VDR and vitamin D signaling phenotypes and will most likely become ever more important to develop new hypotheses over how signaling occurs and where it is most biologically relevant. Perhaps this is readily illustrated by the application of Mendelian randomization methodology to test the factors such as levels of serum 25(OH)D3 (either measured or predicted) and their causal relationship with different clinical phenotypes. In this manner 25(OH)D3 levels and components of the VDR axis are significantly casual in multiple sclerosis, hand grip strength, bacterial infection, and type 2 diabetes [63,64,65,66], but not for the prevention of COVID-19 infection, birth weight, or colon cancer [58,67,68].

2.2. Bottom-Up Approaches Applied to VDR Biology

Applying a bottom-up SB approach to VDR-dependent gene regulation requires constructing a kinetic model that captures the function of different VDR-containing complexes and the quantitative regulation patterns of VDR target genes. Such a model would combine spatial-temporal measurements for the generation of 1α,25(OH)2D3 in a target tissue, VDR genomic interactions, such as initial concentrations in the nucleus of the proteins involved, their half-life, and their associations, to explain how these actions predict the downstream transcriptional impact. Model refinement and validation can then be undertaken by empirical measurements of the cycling of VDR genomic interactions and transcriptional outputs.
The working assumption is that the transcription factor complex contains enzymes that remodel chromatin to allow signaling to the RNA polymerase complex and the initiation of transcription. In the case of a transcription factor such as the VDR, which has nuclear residence independent of ligand exposure, the complex exists in two states of genomic interaction, namely with and without ligand, and in each with qualitatively and quantitatively different cofactors, and also in terms of distribution (distal or proximal). The output of the model will then be to predict the kinetics of mRNA accumulation in a manner that reflects the different ligand-dependent and independent states, and spatial association to target genes. Of course, these components are most likely too simple to predict mRNA accumulation accurately. However, even such a simple model for different genes and the VDR could begin to define how much of the difference between mRNA accumulation could be explained by merely considering spatial and temporal factors for transcription factor residence on a gene. Model components could then be added to consider other well-understood biological components such as disease states, and chromatin states, as well as the impact of DNA helicases required for transcription, or even the role of non-coding RNA to impact mRNA accumulation.
To date, such an approach has been applied to several transcription factors in human cells. Within the nuclear receptor superfamily, kinetic modeling of ERα predicted a role of receptor phosphorylation to act as an underexplored feedback mechanism to predict RNA accumulation and was validated by empirical measurements [69]. Outside of nuclear receptors, other human transcription factor-centered models include describing the separate actions of NF-κB [70,71] and p53 [72], and their collective oscillatory behavior [73]. Separately, signal transduction events have also been modeled quite extensively, including the transduction events for RAS [74,75] and its interaction with other regulatory events such as the circadian clock [76] and EGFR signaling [77].
These approaches are powerful, and it is perhaps surprising why they have not been more widely applied to other nuclear receptors beyond the ERα, for example to the clinically relevant VDR. Several impediments no doubt include the time, resources, and interdisciplinary expertise required for model development and refinement. One impediment is establishing collaborations that bring insight from the wet and dry labs. To be successful these collaborations require sufficient data density of dynamic cellular events, such as precise measurements of the concentration of key proteins and RNA molecules in a cell, the frequency of transcription factor genomic interactions, and levels of mRNA accumulation to justify modelling in mathematical terms. Furthermore, these challenges are sometimes impeded by an uneasy relationship between the biological and mathematical communities [5]. The incentive for these intensive modelling efforts, most likely involving interdisciplinary collaborations, is that once such a model is built it has the potential to be translated to other transcription factor genome interactions, or across cell types, or disease states and allow speculation on the biological processes without having to undertake time-intensive and costly experiments in the wet lab.

3. Prostate Cancer and Health Disparities; An Exemplar of the Opportunities Arising from Systems Biology Approaches in Biomedicine

Men of African genomic ancestry in American (African American, (AA)) experience higher risks of developing more aggressive prostate cancer (PCa) than European American (EA) counterparts [78,79,80,81], which reflects underlying genetic [82,83,84,85] and epigenetic [86,87,88,89,90,91] drivers and biopsychosocial processes [89,92,93]. The role of the glucocorticoid receptor (GR) as a primary target for stress response has more recently been examined in the context of cancer (reviewed in [94]) and offers a potential functional explanation for how stress can be an accelerant of PCa in AA men [94,95,96,97]. A role for VDR signaling appears to be more high profile in the etiology of AA rather than EA PCa. Given that UVB radiation degrades folic acid as well as catalyzing vitamin D synthesis, a strong inverse correlation between skin pigmentation and latitude has arisen during ancestral adaptation [98,99], and amongst AA PCa patients there are significant associations between low serum vitamin D3 levels and incidence and progression risks of PCa [100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115]. Furthermore, although vitamin D3 supplementation in the VITAL cohort [116,117] had no overall impact on cancer incidence, the AA participants experienced a suggestive 23% (p = 0.07) reduction in cancer risk, indicating that larger cohorts may be more informative [118,119,120]. Strikingly, vitamin D3 supplementation in AA and EA PCa patients only significantly modulated prostate gene expression in the AA patients [121], associated with the control of inflammation. This was supported by our recent study [122] in AA and EA cell models which have identified qualitatively and quantitatively distinct VDR actions in terms of VDR protein–protein interactions and more frequent VDR genomic interactions associated with significantly distinct target gene expression.
These data support the concept that one cell function that appears to be most highly responsive in a manner that reflects genomic ancestry is the role of VDR to modulate innate immunity [48,123,124,125,126], and reflecting this there are correlative findings that support a relationship between 25(OH)D3 levels and immune-modulatory factors such as interleukin (IL)-6 [127,128]. IL6 release into the serum is associated with activated macrophages and is elevated in Ghanaian and AA men [129,130,131]. Together these studies suggest that the biology of the prostate is the most sensitive VDR signaling in AA men, and more impacted by inadequate VDR signaling either as a result of molecular mechanisms or environments of low 25(OH)D3. One biologically impactful consequence of this is the loss of control of inflammatory signals. These studies also highlight that African ancestry is itself divergent, and is also further modified by admixture, for example in the AA population (refs).

The Potential Application of SB Approaches to Prostate Health and Disease

There are significant knowledge gaps in identifying how genomic ancestry, the environment, and lifestyle choices may combine to impact prostate health and disease. The application of SB modeling could be used to generate biologically plausible hypotheses to test specific relationships in these complex interactions. These knowledge gaps include whether germline structural variants impact the generation of 1α,25(OH)2D3, and the extent of relationships between serum 25(OH)D3 levels and PCa in a manner that reflects genomic ancestry, and if, and to what extent, this is impacted by admixture. In experimental models, it is emerging that the VDR genomic interactions differ between EA and AA prostate cell models (ref), and recently AR genomic interactions also appear different in EA and AA PCa samples. However, it is unclear if either the underlying sequence, for example at enhancer binding sites, or other mechanisms are impacting why these different genomic interactions arise. Finally, given the multi-parameter process of cis-regulatory interactions that drive transcription through to processed RNA and translation to protein, it is also unclear how genomic ancestry impacts these processes potentially in an emergent manner.
Together, these knowledge gaps can be combined in a hypothesis; namely that PCa risks in AA men reflect the interplay of genomic ancestry, including admixture, coupled with altered environmental signals that combine to drive qualitatively and quantitatively distinct functioning of the VDR and potentially other transcription factors and underscore the health disparities in PCa. Potentially, this is also an attractive setting to develop a series of SB models with which to test these interactions in silico and generate predictions to validate in vitro and in vivo (Figure 1).
This provides an opportunity for both top-down and bottom-up SB modeling. From the top-down perspective, high dimensional data are available to test develop models that test how genomic ancestry significantly determines the associations between serum 25(OH)D3 levels [118,119,120] and PCa structural variants [132,133], epigenetic states [86,87,88,89,90,91] and gene expression [87,91,134] and how this significantly relates to serum inflammatory markers. More specifically, Mendelian randomization approaches would be appropriate to test how germline structural variants are associated with serum 25(OH)D3 levels depending on genomic ancestry.
To complement these approaches, other top-down strategies applied to high dimensional data sets could then determine epigenomes such as CpG methylation, histone modification, and chromatin accessibility [135,136] and transcriptome relationships in AA and EA PCa patients. Construction of a spatial matrix from epigenomic data, co-incidence with prostate cancer-specific enhancers [137,138], will be binned by accounting for orientation and distance to gene features. Machine learning approaches can then be applied to identify and test the significance of relationships between this matrix and PCa transcriptomes. For example, bootstrapping approaches can test the associations of the epigenome and patterns of observed gene expression compared to simulated data by random sampling [139]. In parallel, other approaches such as the Pareto optimization algorithm [140] can define the ordered correlation of relationships between the cistrome matrix and genesets, of which the strongest (positive or negative) are of greatest biological significance. Comparing the results of both approaches within and across genomic ancestry will reveal more significant cistrome–transcriptome features that are potentially driving PCa health disparities. Lasso and ridge regression can then be performed on the most significant cistrome–transcriptome relationship genes to identify gene expression patterns that predict clinical outcomes in publicly available data [141,142].
To meet these top-down approaches, bottom-up methods can be applied to build kinetic models that take into consideration the amount of the VDR and its key interacting factors in AA and EA prostate cells to simulate its ligand-dependent and independent genomic states and distribution across epigenetic states. The goal of this modeling would be to provide biological plausibility for why the VDR appears to be significantly transcriptionally divergent between AA and EA PCa. Such predictions would form the justification for validation in the wet lab using AA prostate cell line models and patient-derived xenografts. Finally, such modeling may well justify tailoring recommendations for healthy vitamin D serum levels that reflect genomic ancestry.

4. Summary

The VDR signaling system is physiologically impactful and when its functions are altered it is associated with several disease states. Systems biology approaches to the analyses of VDR functions have been applied in several top-down studies that have identified significant VDR functions. To date, however, no bottom-up approaches to kinetic modeling have been applied. The incentive in biomedicine to continue to apply SB approaches would be to develop insight into what are the most prominent contexts where the VDR system exerts biological impact and in turn how this can be exploited in diagnostic, prognostic, or therapeutic contexts. For example, the application of SB approaches, both from a top-down and bottom-up perspective, has the potential to be impactful in the context of prostate cancer health disparities.

Author Contributions

Conceptualization C.Y., M.J.C.; Writing and Review S.H., C.Y., M.J.C.; Funding acquisition C.Y., M.J.C.; All authors have read and agreed to the published version of the manuscript.

Funding

M.J.C. acknowledges support in part from the Department of Defense Congressionally Directed Medical Research Programs (W81XWH-20-1-0373; W81XWH-21-1-0555), C.Y. acknowledges support in part from the Department of Defense Congressionally Directed Medical Research Programs (W81XWH-21-1-0850; W81XWH-18-1-0589); from U54-MD007585-26 (NIH/NIMHD), U54 CA118623 (NIH/NCI). S.H. acknowledges support from a Pelotonia Post-Doctoral Fellowship. M.J.C. and S.H. also acknowledge the National Institute of Health Cancer Center Support Grant (P30CA016058) to the OSUCCC The James.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rowley, J.D. Letter: A new consistent chromosomal abnormality in chronic myelogenous leukaemia identified by quinacrine fluorescence and Giemsa staining. Nature 1973, 243, 290–293. [Google Scholar] [CrossRef] [PubMed]
  2. Druker, B.J.; Lydon, N.B. Lessons learned from the development of an abl tyrosine kinase inhibitor for chronic myelogenous leukemia. J. Clin. Investig. 2000, 105, 3–7. [Google Scholar] [CrossRef] [Green Version]
  3. Bixby, D.; Talpaz, M. Seeking the causes and solutions to imatinib-resistance in chronic myeloid leukemia. Leukemia 2011, 25, 7–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Cassman, M. Barriers to progress in systems biology. Nature 2005, 438, 1079. [Google Scholar] [CrossRef] [PubMed]
  5. Vera, J.; Lischer, C.; Nenov, M.; Nikolov, S.; Lai, X.; Eberhardt, M. Mathematical Modelling in Biomedicine: A Primer for the Curious and the Skeptic. Int. J. Mol. Sci. 2021, 22, 547. [Google Scholar] [CrossRef] [PubMed]
  6. Haran, T.K.; Keren, L. From genes to modules, from cells to ecosystems. Mol. Syst. Biol. 2022, 18, e10726. [Google Scholar]
  7. Angione, C. Human Systems Biology and Metabolic Modelling: A Review-From Disease Metabolism to Precision Medicine. Biomed. Res. Int. 2019, 2019, 8304260. [Google Scholar] [CrossRef]
  8. Kuenzi, B.M.; Ideker, T. A census of pathway maps in cancer systems biology. Nat. Rev. Cancer 2020, 20, 233–246. [Google Scholar] [CrossRef]
  9. Wynn, M.L.; Consul, N.; Merajver, S.D.; Schnell, S. Logic-based models in systems biology: A predictive and parameter-free network analysis method. Integr. Biol. 2012, 4, 1323–1337. [Google Scholar] [CrossRef]
  10. Germain, R.N. Will Systems Biology Deliver Its Promise and Contribute to the Development of New or Improved Vaccines? What Really Constitutes the Study of “Systems Biology” and How Might Such an Approach Facilitate Vaccine Design. Cold Spring Harb. Perspect. Biol. 2018, 10, a033308. [Google Scholar] [CrossRef]
  11. Petrasek, D. Systems biology: The case for a systems science approach to diabetes. J. Diabetes Sci. Technol. 2008, 2, 131–134. [Google Scholar] [CrossRef] [PubMed]
  12. Okamura, H.; Yamaguchi, S.; Yagita, K. Molecular machinery of the circadian clock in mammals. Cell Tissue Res. 2002, 309, 47–56. [Google Scholar] [CrossRef] [PubMed]
  13. Tan, D.; Chen, R.; Mo, Y.; Gu, S.; Ma, J.; Xu, W.; Lu, X.; He, H.; Jiang, F.; Fan, W.; et al. Quantitative control of noise in mammalian gene expression by dynamic histone regulation. Elife 2021, 10, e65654. [Google Scholar] [CrossRef] [PubMed]
  14. Thorne, J.L.; Campbell, M.J.; Turner, B.M. Transcription factors, chromatin and cancer. Int. J. Biochem. Cell Biol. 2009, 41, 164–175. [Google Scholar] [CrossRef] [PubMed]
  15. Nurk, S.; Koren, S.; Rhie, A.; Rautiainen, M.; Bzikadze, A.V.; Mikheenko, A.; Vollger, M.R.; Altemose, N.; Uralsky, L.; Gershman, A.; et al. The complete sequence of a human genome. Science 2022, 376, 44–53. [Google Scholar] [CrossRef] [PubMed]
  16. Barabasi, A.L.; Oltvai, Z.N. Network biology: Understanding the cell’s functional organization. Nat. Rev. Genet. 2004, 5, 101–113. [Google Scholar] [CrossRef] [PubMed]
  17. Carlberg, C.; Dunlop, T.W. An integrated biological approach to nuclear receptor signaling in physiological control and disease. Crit. Rev. Eukaryot. Gene Expr. 2006, 16, 1–22. [Google Scholar] [CrossRef] [Green Version]
  18. Flores, M.; Glusman, G.; Brogaard, K.; Price, N.D.; Hood, L. P4 medicine: How systems medicine will transform the healthcare sector and society. Per. Med. 2013, 10, 565–576. [Google Scholar] [CrossRef] [Green Version]
  19. Sobie, E.A.; Lee, Y.S.; Jenkins, S.L.; Iyengar, R. Systems biology—Biomedical modeling. Sci. Signal 2011, 4, tr2. [Google Scholar] [CrossRef] [Green Version]
  20. Westerhoff, H.V.; Palsson, B.O. The evolution of molecular biology into systems biology. Nat. Biotechnol. 2004, 22, 1249–1252. [Google Scholar] [CrossRef]
  21. Bruggeman, F.J.; Westerhoff, H.V. The nature of systems biology. Trends Microbiol. 2007, 15, 45–50. [Google Scholar] [CrossRef] [PubMed]
  22. Snoep, J.L. The Silicon Cell initiative: Working towards a detailed kinetic description at the cellular level. Curr. Opin. Biotechnol. 2005, 16, 336–343. [Google Scholar] [CrossRef] [PubMed]
  23. Ellis, H. Daniel Whistler: English physician who published the first book on rickets in 1645. Br. J. Hosp. Med. 2019, 80, 51. [Google Scholar] [CrossRef] [PubMed]
  24. Carlberg, C.; Campbell, M.J. Vitamin D receptor signaling mechanisms: Integrated actions of a well-defined transcription factor. Steroids 2013, 78, 127–136. [Google Scholar] [CrossRef] [Green Version]
  25. Seuter, S.; Neme, A.; Carlberg, C. Characterization of genomic vitamin D receptor binding sites through chromatin looping and opening. PLoS ONE 2014, 9, e96184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Satoh, J.; Tabunoki, H. Molecular network of chromatin immunoprecipitation followed by deep sequencing-based vitamin D receptor target genes. Mult. Scler. 2013, 19, 1035–1045. [Google Scholar] [CrossRef] [PubMed]
  27. Ramagopalan, S.V.; Heger, A.; Berlanga, A.J.; Maugeri, N.J.; Lincoln, M.R.; Burrell, A.; Handunnetthi, L.; Handel, A.E.; Disanto, G.; Orton, S.M.; et al. A ChIP-seq defined genome-wide map of vitamin D receptor binding: Associations with disease and evolution. Genome Res. 2010, 20, 1352–1360. [Google Scholar] [CrossRef] [Green Version]
  28. Meyer, M.B.; Goetsch, P.D.; Pike, J.W. VDR/RXR and TCF4/beta-catenin cistromes in colonic cells of colorectal tumor origin: Impact on c-FOS and c-MYC gene expression. Mol. Endocrinol. 2012, 26, 37–51. [Google Scholar] [CrossRef]
  29. Gonzalez-Duarte, R.J.; Cazares-Ordonez, V.; Diaz, L.; Ortiz, V.; Larrea, F.; Avila, E. The expression of RNA helicase DDX5 is transcriptionally upregulated by calcitriol through a vitamin D response element in the proximal promoter in SiHa cervical cells. Mol. Cell Biochem. 2015, 410, 65–73. [Google Scholar] [CrossRef]
  30. Wagner, M.; Rid, R.; Maier, C.J.; Maier, R.H.; Laimer, M.; Hintner, H.; Bauer, J.W.; Onder, K. DDX5 is a multifunctional co-activator of steroid hormone receptors. Mol. Cell Endocrinol. 2012, 361, 80–91. [Google Scholar] [CrossRef]
  31. MacDonald, P.N.; Dowd, D.R.; Zhang, C.; Gu, C. Emerging insights into the coactivator role of NCoA62/SKIP in Vitamin D-mediated transcription. J. Steroid Biochem. Mol. Biol. 2004, 89–90, 179–186. [Google Scholar] [CrossRef] [PubMed]
  32. Foulds, C.E.; Panigrahi, A.K.; Coarfa, C.; Lanz, R.B.; O’Malley, B.W. Long Noncoding RNAs as Targets and Regulators of Nuclear Receptors. Curr. Top Microbiol. Immunol. 2016, 394, 143–176. [Google Scholar] [PubMed]
  33. Lanz, R.B.; Razani, B.; Goldberg, A.D.; O’Malley, B.W. Distinct RNA motifs are important for coactivation of steroid hormone receptors by steroid receptor RNA activator (SRA). Proc. Natl. Acad. Sci. USA 2002, 99, 16081–16086. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Statello, L.; Guo, C.J.; Chen, L.L.; Huarte, M. Gene regulation by long non-coding RNAs and its biological functions. Nat. Rev. Mol. Cell Biol. 2021, 22, 96–118. [Google Scholar] [CrossRef]
  35. Novershtern, N.; Subramanian, A.; Lawton, L.N.; Mak, R.H.; Haining, W.N.; McConkey, M.E.; Habib, N.; Yosef, N.; Chang, C.Y.; Shay, T.; et al. Densely interconnected transcriptional circuits control cell states in human hematopoiesis. Cell 2011, 144, 296–309. [Google Scholar] [CrossRef] [Green Version]
  36. Novikova, S.; Tikhonova, O.; Kurbatov, L.; Farafonova, T.; Vakhrushev, I.; Lupatov, A.; Yarygin, K.; Zgoda, V. Omics Technologies to Decipher Regulatory Networks in Granulocytic Cell Differentiation. Biomolecules 2021, 11, 907. [Google Scholar] [CrossRef]
  37. Fuhrken, P.G.; Chen, C.; Apostolidis, P.A.; Wang, M.; Miller, W.M.; Papoutsakis, E.T. Gene Ontology-driven transcriptional analysis of CD34+ cell-initiated megakaryocytic cultures identifies new transcriptional regulators of megakaryopoiesis. Physiol. Genom. 2008, 33, 159–169. [Google Scholar] [CrossRef] [Green Version]
  38. Yetgin, S.; Ozsoylu, S. Myeloid metaplasia in vitamin D deficiency rickets. Scand. J. Haematol. 1982, 28, 180–185. [Google Scholar] [CrossRef]
  39. Koeffler, H.P. Induction of differentiation of human acute myelogenous leukemia cells: Therapeutic implications. Blood 1983, 62, 709–721. [Google Scholar] [CrossRef] [Green Version]
  40. Tyner, J.W.; Tognon, C.E.; Bottomly, D.; Wilmot, B.; Kurtz, S.E.; Savage, S.L.; Long, N.; Schultz, A.R.; Traer, E.; Abel, M.; et al. Functional genomic landscape of acute myeloid leukaemia. Nature 2018, 562, 526–531. [Google Scholar] [CrossRef]
  41. Ji, Y.; Studzinski, G.P. Retinoblastoma protein and CCAAT/enhancer-binding protein beta are required for 1,25-dihydroxyvitamin D3-induced monocytic differentiation of HL60 cells. Cancer Res. 2004, 64, 370–377. [Google Scholar] [CrossRef] [PubMed]
  42. Shiozaki, Y.; Miyazaki-Anzai, S.; Keenan, A.L.; Miyazaki, M. MEF2D-NR4A1-FAM134B2-mediated reticulophagy contributes to amino acid homeostasis. Autophagy 2022, 18, 1049–1061. [Google Scholar] [CrossRef]
  43. Marchwicka, A.; Marcinkowska, E. Regulation of Expression of CEBP Genes by Variably Expressed Vitamin D Receptor and Retinoic Acid Receptor alpha in Human Acute Myeloid Leukemia Cell Lines. Int. J. Mol. Sci. 2018, 19, 1918. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Mansure, J.J.; Nassim, R.; Chevalier, S.; Szymanski, K.; Rocha, J.; Aldousari, S.; Kassouf, W. A novel mechanism of PPAR gamma induction via EGFR signalling constitutes rational for combination therapy in bladder cancer. PLoS ONE 2013, 8, e55997. [Google Scholar] [CrossRef]
  45. Tang, Q.Q.; Otto, T.C.; Lane, M.D. CCAAT/enhancer-binding protein beta is required for mitotic clonal expansion during adipogenesis. Proc. Natl. Acad. Sci. USA 2003, 100, 850–855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Zolotarenko, A.; Chekalin, E.; Mesentsev, A.; Kiseleva, L.; Gribanova, E.; Mehta, R.; Baranova, A.; Tatarinova, T.V.; Piruzian, E.S.; Bruskin, S. Integrated computational approach to the analysis of RNA-seq data reveals new transcriptional regulators of psoriasis. Exp. Mol. Med. 2016, 48, e268. [Google Scholar] [CrossRef] [Green Version]
  47. Greenhalgh, I.; Butler, A.R. Sanatoria revisited: Sunlight and health. J. R. Coll. Phys. Edinb. 2017, 47, 276–280. [Google Scholar] [CrossRef]
  48. Liu, P.T.; Stenger, S.; Li, H.; Wenzel, L.; Tan, B.H.; Krutzik, S.R.; Ochoa, M.T.; Schauber, J.; Wu, K.; Meinken, C.; et al. Toll-like receptor triggering of a vitamin D-mediated human antimicrobial response. Science 2006, 311, 1770–1773. [Google Scholar] [CrossRef]
  49. Mvubu, N.E.; Pillay, B.; Gamieldien, J.; Bishai, W.; Pillay, M. Canonical pathways, networks and transcriptional factor regulation by clinical strains of Mycobacterium tuberculosis in pulmonary alveolar epithelial cells. Tuberculosis 2016, 97, 73–85. [Google Scholar] [CrossRef]
  50. Sims, A.C.; Tilton, S.C.; Menachery, V.D.; Gralinski, L.E.; Schafer, A.; Matzke, M.M.; Webb-Robertson, B.J.; Chang, J.; Luna, M.L.; Long, C.E.; et al. Release of severe acute respiratory syndrome coronavirus nuclear import block enhances host transcription in human lung cells. J. Virol. 2013, 87, 3885–3902. [Google Scholar] [CrossRef] [Green Version]
  51. Chauss, D.; Freiwald, T.; McGregor, R.; Yan, B.; Wang, L.; Nova-Lamperti, E.; Kumar, D.; Zhang, Z.; Teague, H.; West, E.E.; et al. Autocrine vitamin D signaling switches off pro-inflammatory programs of TH1 cells. Nat. Immunol. 2022, 23, 62–74. [Google Scholar] [CrossRef]
  52. Adolphe, C.; Xue, A.; Fard, A.T.; Genovesi, L.A.; Yang, J.; Wainwright, B.J. Genetic and functional interaction network analysis reveals global enrichment of regulatory T cell genes influencing basal cell carcinoma susceptibility. Genome Med. 2021, 13, 19. [Google Scholar] [CrossRef] [PubMed]
  53. Jostins, L.; Ripke, S.; Weersma, R.K.; Duerr, R.H.; McGovern, D.P.; Hui, K.Y.; Lee, J.C.; Schumm, L.P.; Sharma, Y.; Anderson, C.A.; et al. Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature 2012, 491, 119–124. [Google Scholar] [CrossRef] [Green Version]
  54. Wang, J.; Thingholm, L.B.; Skieceviciene, J.; Rausch, P.; Kummen, M.; Hov, J.R.; Degenhardt, F.; Heinsen, F.A.; Ruhlemann, M.C.; Szymczak, S.; et al. Genome-wide association analysis identifies variation in vitamin D receptor and other host factors influencing the gut microbiota. Nat. Genet. 2016, 48, 1396–1406. [Google Scholar] [CrossRef]
  55. Tavasolian, F.; Hosseini, A.Z.; Soudi, S.; Naderi, M.; Sahebkar, A. A Systems Biology Approach for miRNA-mRNA Expression Patterns Analysis in Rheumatoid Arthritis. Comb. Chem. High Throughput. Screen 2021, 24, 195–212. [Google Scholar] [CrossRef] [PubMed]
  56. Montasser, M.E.; Aslibekyan, S.; Srinivasasainagendra, V.; Tiwari, H.K.; Patki, A.; Bagheri, M.; Kind, T.; Barupal, D.K.; Fan, S.; Perry, J.; et al. An Amish founder population reveals rare-population genetic determinants of the human lipidome. Commun. Biol. 2022, 5, 334. [Google Scholar] [CrossRef] [PubMed]
  57. Skinkyte-Juskiene, R.; Kogelman, L.J.A.; Kadarmideen, H.N. Transcription Factor Co-expression Networks of Adipose RNA-Seq Data Reveal Regulatory Mechanisms of Obesity. Curr. Genom. 2018, 19, 289–299. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Bao, M.H.; Luo, H.Q.; Chen, L.H.; Tang, L.; Ma, K.F.; Xiang, J.; Dong, L.P.; Zeng, J.; Li, G.Y.; Li, J.M. Impact of high fat diet on long non-coding RNAs and messenger RNAs expression in the aortas of ApoE(-/-) mice. Sci. Rep. 2016, 6, 34161. [Google Scholar] [CrossRef] [Green Version]
  59. Campbell, M.J. Bioinformatic approaches to interrogating vitamin D receptor signaling. Mol. Cell Endocrinol. 2017, 453, 3–13. [Google Scholar] [CrossRef]
  60. Vuckovic, D.; Bao, E.L.; Akbari, P.; Lareau, C.A.; Mousas, A.; Jiang, T.; Chen, M.H.; Raffield, L.M.; Tardaguila, M.; Huffman, J.E.; et al. The Polygenic and Monogenic Basis of Blood Traits and Diseases. Cell 2020, 182, 1214–1231.e1211. [Google Scholar] [CrossRef]
  61. Singh, P.K.; van den Berg, P.R.; Long, M.D.; Vreugdenhil, A.; Grieshober, L.; Ochs-Balcom, H.M.; Wang, J.; Delcambre, S.; Heikkinen, S.; Carlberg, C.; et al. Integration of VDR genome wide binding and GWAS genetic variation data reveals co-occurrence of VDR and NF-kappaB binding that is linked to immune phenotypes. BMC Genom. 2017, 18, 132. [Google Scholar] [CrossRef] [PubMed]
  62. Alleyne, D.; Witonsky, D.B.; Mapes, B.; Nakagome, S.; Sommars, M.; Hong, E.; Muckala, K.A.; Di Rienzo, A.; Kupfer, S.S. Colonic transcriptional response to 1alpha,25(OH)2 vitamin D3 in African- and European-Americans. J. Steroid Biochem. Mol. Biol. 2017, 168, 49–59. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Colak, Y.; Nordestgaard, B.G.; Afzal, S. Low vitamin D and risk of bacterial pneumonias: Mendelian randomisation studies in two population-based cohorts. Thorax 2021, 76, 468–478. [Google Scholar] [CrossRef] [PubMed]
  64. Mazidi, M.; Davies, I.G.; Penson, P.; Rikkonen, T.; Isanejad, M. Lifetime serum concentration of 25-hydroxyvitamin D 25(OH) is associated with hand grip strengths: Insight from a Mendelian randomisation. Age Ageing 2022, 51, afac079. [Google Scholar] [CrossRef]
  65. Mokry, L.E.; Ross, S.; Ahmad, O.S.; Forgetta, V.; Smith, G.D.; Goltzman, D.; Leong, A.; Greenwood, C.M.; Thanassoulis, G.; Richards, J.B. Vitamin D and Risk of Multiple Sclerosis: A Mendelian Randomization Study. PLoS Med. 2015, 12, e1001866. [Google Scholar] [CrossRef] [Green Version]
  66. Zheng, J.S.; Luan, J.; Sofianopoulou, E.; Sharp, S.J.; Day, F.R.; Imamura, F.; Gundersen, T.E.; Lotta, L.A.; Sluijs, I.; Stewart, I.D.; et al. The association between circulating 25-hydroxyvitamin D metabolites and type 2 diabetes in European populations: A meta-analysis and Mendelian randomisation analysis. PLoS Med. 2020, 17, e1003394. [Google Scholar] [CrossRef]
  67. He, Y.; Timofeeva, M.; Farrington, S.M.; Vaughan-Shaw, P.; Svinti, V.; Walker, M.; Zgaga, L.; Meng, X.; Li, X.; Spiliopoulou, A.; et al. Exploring causality in the association between circulating 25-hydroxyvitamin D and colorectal cancer risk: A large Mendelian randomisation study. BMC Med. 2018, 16, 142. [Google Scholar] [CrossRef]
  68. Thompson, W.D.; Tyrrell, J.; Borges, M.C.; Beaumont, R.N.; Knight, B.A.; Wood, A.R.; Ring, S.M.; Hattersley, A.T.; Freathy, R.M.; Lawlor, D.A. Association of maternal circulating 25(OH)D and calcium with birth weight: A mendelian randomisation analysis. PLoS Med. 2019, 16, e1002828. [Google Scholar] [CrossRef] [Green Version]
  69. Tian, D.; Solodin, N.M.; Rajbhandari, P.; Bjorklund, K.; Alarid, E.T.; Kreeger, P.K. A kinetic model identifies phosphorylated estrogen receptor-alpha (ERalpha) as a critical regulator of ERalpha dynamics in breast cancer. FASEB J. 2015, 29, 2022–2031. [Google Scholar] [CrossRef] [Green Version]
  70. Cheong, R.; Hoffmann, A.; Levchenko, A. Understanding NF-kappaB signaling via mathematical modeling. Mol. Syst. Biol. 2008, 4, 192. [Google Scholar] [CrossRef]
  71. Bullock, M.E.; Moreno-Martinez, N.; Miller-Jensen, K. A transcriptional cycling model recapitulates chromatin-dependent features of noisy inducible transcription. PLoS Comput. Biol. 2022, 18, e1010152. [Google Scholar] [CrossRef] [PubMed]
  72. Jimenez, A.; Lu, D.; Kalocsay, M.; Berberich, M.J.; Balbi, P.; Jambhekar, A.; Lahav, G. Time-series transcriptomics and proteomics reveal alternative modes to decode p53 oscillations. Mol. Syst. Biol. 2022, 18, e10588. [Google Scholar] [CrossRef] [PubMed]
  73. Ihekwaba, A.E.; Sedwards, S. Communicating oscillatory networks: Frequency domain analysis. BMC Syst. Biol. 2011, 5, 203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Markevich, N.I.; Moehren, G.; Demin, O.V.; Kiyatkin, A.; Hoek, J.B.; Kholodenko, B.N. Signal processing at the Ras circuit: What shapes Ras activation patterns? Syst. Biol. 2004, 1, 104–113. [Google Scholar] [CrossRef]
  75. Stelniec-Klotz, I.; Legewie, S.; Tchernitsa, O.; Witzel, F.; Klinger, B.; Sers, C.; Herzel, H.; Bluthgen, N.; Schafer, R. Reverse engineering a hierarchical regulatory network downstream of oncogenic KRAS. Mol. Syst. Biol. 2012, 8, 601. [Google Scholar] [CrossRef]
  76. Relogio, A.; Thomas, P.; Medina-Perez, P.; Reischl, S.; Bervoets, S.; Gloc, E.; Riemer, P.; Mang-Fatehi, S.; Maier, B.; Schafer, R.; et al. Ras-mediated deregulation of the circadian clock in cancer. PLoS Genet. 2014, 10, e1004338. [Google Scholar] [CrossRef]
  77. Davies, A.E.; Pargett, M.; Siebert, S.; Gillies, T.E.; Choi, Y.; Tobin, S.J.; Ram, A.R.; Murthy, V.; Juliano, C.; Quon, G.; et al. Systems-Level Properties of EGFR-RAS-ERK Signaling Amplify Local Signals to Generate Dynamic Gene Expression Heterogeneity. Cell Syst. 2020, 11, 161–175.e165. [Google Scholar] [CrossRef]
  78. Patel, H.D.; Doshi, C.P.; Koehne, E.L.; Hart, S.; Van Kuiken, M.; Quek, M.L.; Flanigan, R.C.; Gupta, G.N. African American Men have Increased Risk of Prostate Cancer Detection Despite Similar Rates of Anterior Prostatic Lesions and PI-RADS Grade on Multiparametric Magnetic Resonance Imaging. Urology 2022, 163, 132–137. [Google Scholar] [CrossRef]
  79. Powell, I.J. Epidemiology and pathophysiology of prostate cancer in African-American men. J. Urol. 2007, 177, 444–449. [Google Scholar] [CrossRef]
  80. Moul, J.W.; Sesterhenn, I.A.; Connelly, R.R.; Douglas, T.; Srivastava, S.; Mostofi, F.K.; McLeod, D.G. Prostate-specific antigen values at the time of prostate cancer diagnosis in African-American men. JAMA 1995, 274, 1277–1281. [Google Scholar] [CrossRef]
  81. Sanchez-Ortiz, R.F.; Troncoso, P.; Babaian, R.J.; Lloreta, J.; Johnston, D.A.; Pettaway, C.A. African-American men with nonpalpable prostate cancer exhibit greater tumor volume than matched white men. Cancer 2006, 107, 75–82. [Google Scholar] [CrossRef] [PubMed]
  82. Fiorica, P.N.; Schubert, R.; Morris, J.D.; Abdul Sami, M.; Wheeler, H.E. Multi-ethnic transcriptome-wide association study of prostate cancer. PLoS ONE 2020, 15, e0236209. [Google Scholar] [CrossRef] [PubMed]
  83. Hoffmann, T.J.; Van Den Eeden, S.K.; Sakoda, L.C.; Jorgenson, E.; Habel, L.A.; Graff, R.E.; Passarelli, M.N.; Cario, C.L.; Emami, N.C.; Chao, C.R.; et al. A large multiethnic genome-wide association study of prostate cancer identifies novel risk variants and substantial ethnic differences. Cancer Discov. 2015, 5, 878–891. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Helfand, B.T.; Roehl, K.A.; Cooper, P.R.; McGuire, B.B.; Fitzgerald, L.M.; Cancel-Tassin, G.; Cornu, J.N.; Bauer, S.; Van Blarigan, E.L.; Chen, X.; et al. Associations of prostate cancer risk variants with disease aggressiveness: Results of the NCI-SPORE Genetics Working Group analysis of 18,343 cases. Hum. Genet. 2015, 134, 439–450. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Mori, J.O.; White, J.; Elhussin, I.; Duduyemi, B.M.; Karanam, B.; Yates, C.; Wang, H. Molecular and pathological subtypes related to prostate cancer disparities and disease outcomes in African American and European American patients. Front. Oncol. 2022, 12, 928357. [Google Scholar] [CrossRef]
  86. Awasthi, S.; Berglund, A.; Abraham-Miranda, J.; Rounbehler, R.J.; Kensler, K.; Serna, A.; Vidal, A.; You, S.; Freeman, M.R.; Davicioni, E.; et al. Comparative Genomics Reveals Distinct Immune-oncologic Pathways in African American Men with Prostate Cancer. Clin. Cancer Res. 2021, 27, 320–329. [Google Scholar] [CrossRef]
  87. Yuan, J.; Kensler, K.H.; Hu, Z.; Zhang, Y.; Zhang, T.; Jiang, J.; Xu, M.; Pan, Y.; Long, M.; Montone, K.T.; et al. Integrative comparison of the genomic and transcriptomic landscape between prostate cancer patients of predominantly African or European genetic ancestry. PLoS Genet. 2020, 16, e1008641. [Google Scholar] [CrossRef] [Green Version]
  88. Hardiman, G.; Savage, S.J.; Hazard, E.S.; Wilson, R.C.; Courtney, S.M.; Smith, M.T.; Hollis, B.W.; Halbert, C.H.; Gattoni-Celli, S. Systems analysis of the prostate transcriptome in African-American men compared with European-American men. Pharmacogenomics 2016, 17, 1129–1143. [Google Scholar] [CrossRef]
  89. Yamoah, K.; Asamoah, F.A.; Abrahams, A.O.D.; Awasthi, S.; Mensah, J.E.; Dhillon, J.; Mahal, B.A.; Gueye, S.M.; Jalloh, M.; Farahani, S.J.; et al. Prostate tumors of native men from West Africa show biologically distinct pathways-A comparative genomic study. Prostate 2021, 81, 1402–1410. [Google Scholar] [CrossRef]
  90. Koga, Y.; Song, H.; Chalmers, Z.R.; Newberg, J.; Kim, E.; Carrot-Zhang, J.; Piou, D.; Polak, P.; Abdulkadir, S.A.; Ziv, E.; et al. Genomic Profiling of Prostate Cancers from Men with African and European Ancestry. Clin. Cancer Res. 2020, 26, 4651–4660. [Google Scholar] [CrossRef]
  91. Rayford, W.; Beksac, A.T.; Alger, J.; Alshalalfa, M.; Ahmed, M.; Khan, I.; Falagario, U.G.; Liu, Y.; Davicioni, E.; Spratt, D.E.; et al. Comparative analysis of 1152 African-American and European-American men with prostate cancer identifies distinct genomic and immunological differences. Commun. Biol. 2021, 4, 670. [Google Scholar] [CrossRef]
  92. Liu, W.; Zheng, S.L.; Na, R.; Wei, L.; Sun, J.; Gallagher, J.; Wei, J.; Resurreccion, W.K.; Ernst, S.; Sfanos, K.S.; et al. Distinct Genomic Alterations in Prostate Tumors Derived from African American Men. Mol. Cancer Res. 2020, 18, 1815–1824. [Google Scholar] [CrossRef]
  93. Magi-Galluzzi, C.; Tsusuki, T.; Elson, P.; Simmerman, K.; LaFargue, C.; Esgueva, R.; Klein, E.; Rubin, M.A.; Zhou, M. TMPRSS2-ERG gene fusion prevalence and class are significantly different in prostate cancer of Caucasian, African-American and Japanese patients. Prostate 2011, 71, 489–497. [Google Scholar] [CrossRef] [PubMed]
  94. Kerkvliet, C.P.; Truong, T.H.; Ostrander, J.H.; Lange, C.A. Stress sensing within the breast tumor microenvironment: How glucocorticoid receptors live in the moment. Essays Biochem. 2021, 65, 971–983. [Google Scholar] [CrossRef] [PubMed]
  95. Lee, M.J.; Rittschof, C.C.; Greenlee, A.J.; Turi, K.N.; Rodriguez-Zas, S.L.; Robinson, G.E.; Cole, S.W.; Mendenhall, R. Transcriptomic analyses of black women in neighborhoods with high levels of violence. Psychoneuroendocrinology 2021, 127, 105174. [Google Scholar] [CrossRef]
  96. Woods-Burnham, L.; Stiel, L.; Martinez, S.R.; Sanchez-Hernandez, E.S.; Ruckle, H.C.; Almaguel, F.G.; Stern, M.C.; Roberts, L.R.; Williams, D.R.; Montgomery, S.; et al. Psychosocial Stress, Glucocorticoid Signaling, and Prostate Cancer Health Disparities in African American Men. Cancer Health Disparities 2020, 4, e1–e30. [Google Scholar]
  97. Woods-Burnham, L.; Cajigas-Du Ross, C.K.; Love, A.; Basu, A.; Sanchez-Hernandez, E.S.; Martinez, S.R.; Ortiz-Hernandez, G.L.; Stiel, L.; Duran, A.M.; Wilson, C.; et al. Glucocorticoids Induce Stress Oncoproteins Associated with Therapy-Resistance in African American and European American Prostate Cancer Cells. Sci. Rep. 2018, 8, 15063. [Google Scholar] [CrossRef] [Green Version]
  98. Hochberg, Z.; Templeton, A.R. Evolutionary perspective in skin color, vitamin D and its receptor. Hormones 2010, 9, 307–311. [Google Scholar] [CrossRef] [Green Version]
  99. Jablonski, N.G. The evolution of human skin colouration and its relevance to health in the modern world. J. R. Coll. Phys. Edinb. 2012, 42, 58–63. [Google Scholar] [CrossRef] [Green Version]
  100. Yao, S.; Hong, C.C.; Bandera, E.V.; Zhu, Q.; Liu, S.; Cheng, T.D.; Zirpoli, G.; Haddad, S.A.; Lunetta, K.L.; Ruiz-Narvaez, E.A.; et al. Demographic, lifestyle, and genetic determinants of circulating concentrations of 25-hydroxyvitamin D and vitamin D-binding protein in African American and European American women. Am. J. Clin. Nutr. 2017, 105, 1362–1371. [Google Scholar] [CrossRef] [Green Version]
  101. Yao, S.; Haddad, S.A.; Hu, Q.; Liu, S.; Lunetta, K.L.; Ruiz-Narvaez, E.A.; Hong, C.C.; Zhu, Q.; Sucheston-Campbell, L.; Cheng, T.Y.; et al. Genetic variations in vitamin D-related pathways and breast cancer risk in African American women in the AMBER consortium. Int. J. Cancer 2016, 138, 2118–2126. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  102. Mishra, D.K.; Wu, Y.; Sarkissyan, M.; Sarkissyan, S.; Chen, Z.; Shang, X.; Ong, M.; Heber, D.; Koeffler, H.P.; Vadgama, J.V. Vitamin D receptor gene polymorphisms and prognosis of breast cancer among African-American and Hispanic women. PLoS ONE 2013, 8, e57967. [Google Scholar] [CrossRef] [PubMed]
  103. Neuhouser, M.L.; Bernstein, L.; Hollis, B.W.; Xiao, L.; Ambs, A.; Baumgartner, K.; Baumgartner, R.; McTiernan, A.; Ballard-Barbash, R. Serum vitamin D and breast density in breast cancer survivors. Cancer Epidemiol. Biomark. Prev. 2010, 19, 412–417. [Google Scholar] [CrossRef] [Green Version]
  104. John, E.M.; Schwartz, G.G.; Koo, J.; Wang, W.; Ingles, S.A. Sun exposure, vitamin D receptor gene polymorphisms, and breast cancer risk in a multiethnic population. Am. J. Epidemiol. 2007, 166, 1409–1419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Layne, T.M.; Weinstein, S.J.; Graubard, B.I.; Ma, X.; Mayne, S.T.; Albanes, D. Serum 25-hydroxyvitamin D, vitamin D binding protein, and prostate cancer risk in black men. Cancer 2017, 123, 2698–2704. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Nelson, S.M.; Batai, K.; Ahaghotu, C.; Agurs-Collins, T.; Kittles, R.A. Association between Serum 25-Hydroxy-Vitamin D and Aggressive Prostate Cancer in African American Men. Nutrients 2016, 9, 12. [Google Scholar] [CrossRef] [Green Version]
  107. Batai, K.; Murphy, A.B.; Nonn, L.; Kittles, R.A. Vitamin D and Immune Response: Implications for Prostate Cancer in African Americans. Front. Immunol. 2016, 7, 53. [Google Scholar] [CrossRef] [Green Version]
  108. Murphy, A.B.; Nyame, Y.; Martin, I.K.; Catalona, W.J.; Hollowell, C.M.; Nadler, R.B.; Kozlowski, J.M.; Perry, K.T.; Kajdacsy-Balla, A.; Kittles, R. Vitamin D deficiency predicts prostate biopsy outcomes. Clin. Cancer Res. 2014, 20, 2289–2299. [Google Scholar] [CrossRef] [Green Version]
  109. Richards, Z.; Batai, K.; Farhat, R.; Shah, E.; Makowski, A.; Gann, P.H.; Kittles, R.; Nonn, L. Prostatic compensation of the vitamin D axis in African American men. JCI Insight 2017, 2, e91054. [Google Scholar] [CrossRef] [Green Version]
  110. Kidd, L.C.; Paltoo, D.N.; Wang, S.; Chen, W.; Akereyeni, F.; Isaacs, W.; Ahaghotu, C.; Kittles, R. Sequence variation within the 5′ regulatory regions of the vitamin D binding protein and receptor genes and prostate cancer risk. Prostate 2005, 64, 272–282. [Google Scholar] [CrossRef]
  111. Murphy, A.B.; Kelley, B.; Nyame, Y.A.; Martin, I.K.; Smith, D.J.; Castaneda, L.; Zagaja, G.J.; Hollowell, C.M.; Kittles, R.A. Predictors of serum vitamin D levels in African American and European American men in Chicago. Am. J. Mens. Health 2012, 6, 420–426. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  112. Paller, C.J.; Kanaan, Y.M.; Beyene, D.A.; Naab, T.J.; Copeland, R.L.; Tsai, H.L.; Kanarek, N.F.; Hudson, T.S. Risk of prostate cancer in African-American men: Evidence of mixed effects of dietary quercetin by serum vitamin D status. Prostate 2015, 75, 1376–1383. [Google Scholar] [CrossRef] [PubMed]
  113. Batai, K.; Kittles, R.A. Can vitamin D supplementation reduce prostate cancer disparities? Pharmacogenomics 2016, 17, 1117–1120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  114. Taksler, G.B.; Cutler, D.M.; Giovannucci, E.; Smith, M.R.; Keating, N.L. Ultraviolet index and racial differences in prostate cancer incidence and mortality. Cancer 2013, 119, 3195–3203. [Google Scholar] [CrossRef]
  115. Steck, S.E.; Arab, L.; Zhang, H.; Bensen, J.T.; Fontham, E.T.; Johnson, C.S.; Mohler, J.L.; Smith, G.J.; Su, J.L.; Trump, D.L.; et al. Association between Plasma 25-Hydroxyvitamin D, Ancestry and Aggressive Prostate Cancer among African Americans and European Americans in PCaP. PLoS ONE 2015, 10, e0125151. [Google Scholar] [CrossRef] [Green Version]
  116. King, L.; Dear, K.; Harrison, S.L.; van der Mei, I.; Brodie, A.M.; Kimlin, M.G.; Lucas, R.M. Investigating the patterns and determinants of seasonal variation in vitamin D status in Australian adults: The Seasonal D Cohort Study. BMC Public Health 2016, 16, 892. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  117. Diffey, B.L. Modelling the seasonal variation of vitamin D due to sun exposure. Br. J. Dermatol. 2010, 162, 1342–1348. [Google Scholar] [CrossRef]
  118. Manson, J.E.; Cook, N.R.; Lee, I.M.; Christen, W.; Bassuk, S.S.; Mora, S.; Gibson, H.; Albert, C.M.; Gordon, D.; Copeland, T.; et al. Marine n-3 Fatty Acids and Prevention of Cardiovascular Disease and Cancer. N. Engl. J. Med. 2019, 380, 23–32. [Google Scholar] [CrossRef] [PubMed]
  119. Bassuk, S.S.; Chandler, P.D.; Buring, J.E.; Manson, J.E.; Group, V.R. The VITamin D and OmegA-3 TriaL (VITAL): Do Results Differ by Sex or Race/Ethnicity? Am. J. Lifestyle Med. 2021, 15, 372–391. [Google Scholar] [CrossRef]
  120. Chandler, P.D.; Chen, W.Y.; Ajala, O.N.; Hazra, A.; Cook, N.; Bubes, V.; Lee, I.M.; Giovannucci, E.L.; Willett, W.; Buring, J.E.; et al. Effect of Vitamin D3 Supplements on Development of Advanced Cancer: A Secondary Analysis of the VITAL Randomized Clinical Trial. JAMA Netw. Open 2020, 3, e2025850. [Google Scholar] [CrossRef]
  121. Marshall, D.T.; Savage, S.J.; Garrett-Mayer, E.; Keane, T.E.; Hollis, B.W.; Horst, R.L.; Ambrose, L.H.; Kindy, M.S.; Gattoni-Celli, S. Vitamin D3 supplementation at 4000 international units per day for one year results in a decrease of positive cores at repeat biopsy in subjects with low-risk prostate cancer under active surveillance. J. Clin. Endocrinol. Metab. 2012, 97, 2315–2324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  122. Siddappa, M.; Hussain, S.; Wani, S.A.; Tang, H.; Gray, J.S.; Jafari, H.; Wu, H.; Long, M.D.; Elhussin, I.; Karanam, B.; et al. Vitamin D receptor cistrome-transcriptome analyses establishes quantitatively distinct receptor genomic interactions in African American prostate cancer regulated by BAZ1A. bioRxiv 2022. [Google Scholar] [CrossRef]
  123. Singh, A.K.; Prakash, S.; Garg, R.K.; Jain, P.; Kumar, R.; Jain, A. Polymorphisms in vitamin D receptor, toll-like receptor 2 and Toll-Like receptor 4 genes links with Dengue susceptibility. Bioinformation 2021, 17, 506–513. [Google Scholar] [CrossRef] [PubMed]
  124. He, L.; Zhou, M.; Li, Y.C. Vitamin D/Vitamin D Receptor Signaling Is Required for Normal Development and Function of Group 3 Innate Lymphoid Cells in the Gut. iScience 2019, 17, 119–131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  125. Konya, V.; Czarnewski, P.; Forkel, M.; Rao, A.; Kokkinou, E.; Villablanca, E.J.; Almer, S.; Lindforss, U.; Friberg, D.; Hoog, C.; et al. Vitamin D downregulates the IL-23 receptor pathway in human mucosal group 3 innate lymphoid cells. J. Allergy Clin. Immunol. 2018, 141, 279–292. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  126. Gombart, A.F.; Borregaard, N.; Koeffler, H.P. Human cathelicidin antimicrobial peptide (CAMP) gene is a direct target of the vitamin D receptor and is strongly up-regulated in myeloid cells by 1,25-dihydroxyvitamin D3. FASEB J. 2005, 19, 1067–1077. [Google Scholar] [CrossRef] [Green Version]
  127. Liao, S.L.; Lai, S.H.; Tsai, M.H.; Hua, M.C.; Yeh, K.W.; Su, K.W.; Chiang, C.H.; Huang, S.Y.; Kao, C.C.; Yao, T.C.; et al. Maternal Vitamin D Level Is Associated with Viral Toll-Like Receptor Triggered IL-10 Response but Not the Risk of Infectious Diseases in Infancy. Mediat. Inflamm. 2016, 2016, 8175898. [Google Scholar] [CrossRef] [Green Version]
  128. Ojaimi, S.; Skinner, N.A.; Strauss, B.J.; Sundararajan, V.; Woolley, I.; Visvanathan, K. Vitamin D deficiency impacts on expression of toll-like receptor-2 and cytokine profile: A pilot study. J. Transl. Med. 2013, 11, 176. [Google Scholar] [CrossRef] [Green Version]
  129. Darko, S.N.; Yar, D.D.; Owusu-Dabo, E.; Awuah, A.A.; Dapaah, W.; Addofoh, N.; Salifu, S.P.; Awua-Boateng, N.Y.; Adomako-Boateng, F. Variations in levels of IL-6 and TNF-alpha in type 2 diabetes mellitus between rural and urban Ashanti Region of Ghana. BMC Endocr. Disord. 2015, 15, 50. [Google Scholar] [CrossRef] [Green Version]
  130. Mensah, G.I.; Addo, K.K.; Tetteh, J.A.; Sowah, S.; Loescher, T.; Geldmacher, C.; Jackson-Sillah, D. Cytokine response to selected MTB antigens in Ghanaian TB patients, before and at 2 weeks of anti-TB therapy is characterized by high expression of IFN-gamma and Granzyme B and inter- individual variation. BMC Infect. Dis. 2014, 14, 495. [Google Scholar] [CrossRef] [Green Version]
  131. Minas, T.Z.; Candia, J.; Dorsey, T.H.; Baker, F.; Tang, W.; Kiely, M.; Smith, C.J.; Zhang, A.L.; Jordan, S.V.; Obadi, O.M.; et al. Serum proteomics links suppression of tumor immunity to ancestry and lethal prostate cancer. Nat. Commun. 2022, 13, 1759. [Google Scholar] [CrossRef]
  132. Chen, F.; Darst, B.F.; Madduri, R.K.; Rodriguez, A.A.; Sheng, X.; Rentsch, C.T.; Andrews, C.; Tang, W.; Kibel, A.S.; Plym, A.; et al. Validation of a multi-ancestry polygenic risk score and age-specific risks of prostate cancer: A meta-analysis within diverse populations. Elife 2022, 11, e78304. [Google Scholar] [CrossRef] [PubMed]
  133. Conti, D.V.; Darst, B.F.; Moss, L.C.; Saunders, E.J.; Sheng, X.; Chou, A.; Schumacher, F.R.; Olama, A.A.A.; Benlloch, S.; Dadaev, T.; et al. Trans-ancestry genome-wide association meta-analysis of prostate cancer identifies new susceptibility loci and informs genetic risk prediction. Nat. Genet. 2021, 53, 65–75. [Google Scholar] [CrossRef] [PubMed]
  134. Berchuck, J.E.; Adib, E.; Abou Alaiwi, S.; Dash, A.K.; Shin, J.N.; Lowder, D.; McColl, C.; Castro, P.; Carelli, R.; Benedetti, E.; et al. The prostate cancer androgen receptor cistrome in African American men associates with upregulation of lipid metabolism and immune response. Cancer Res. 2022, 82, 2848–2859. [Google Scholar] [CrossRef]
  135. Ernst, J.; Kellis, M. Chromatin-state discovery and genome annotation with ChromHMM. Nat. Protoc. 2017, 12, 2478–2492. [Google Scholar] [CrossRef]
  136. Ernst, J.; Kellis, M. ChromHMM: Automating chromatin-state discovery and characterization. Nat. Methods 2012, 9, 215–216. [Google Scholar] [CrossRef] [Green Version]
  137. Long, M.D.; Jacobi, J.J.; Singh, P.K.; Llimos, G.; Wani, S.A.; Rowsam, A.M.; Rosario, S.R.; Hoogstraat, M.; Linder, S.; Kirk, J.; et al. Reduced NCOR2 expression accelerates androgen deprivation therapy failure in prostate cancer. Cell Rep. 2021, 37, 110109. [Google Scholar] [CrossRef]
  138. Pomerantz, M.M.; Qiu, X.; Zhu, Y.; Takeda, D.Y.; Pan, W.; Baca, S.C.; Gusev, A.; Korthauer, K.D.; Severson, T.M.; Ha, G.; et al. Prostate cancer reactivates developmental epigenomic programs during metastatic progression. Nat. Genet. 2020, 52, 790–799. [Google Scholar] [CrossRef] [PubMed]
  139. Long, M.D.; van den Berg, P.R.; Russell, J.L.; Singh, P.K.; Battaglia, S.; Campbell, M.J. Integrative genomic analysis in K562 chronic myelogenous leukemia cells reveals that proximal NCOR1 binding positively regulates genes that govern erythroid differentiation and Imatinib sensitivity. Nucleic Acids Res. 2015, 43, 7330–7348. [Google Scholar] [CrossRef] [PubMed]
  140. Cao, Y.; Kitanovski, S.; Hoffmann, D. intePareto: An R package for integrative analyses of RNA-Seq and ChIP-Seq data. BMC Genom. 2020, 21, 802. [Google Scholar] [CrossRef]
  141. Liu, J.; Lei, B.; Yu, X.; Li, Y.; Deng, Y.; Yang, G.; Li, Z.; Liu, T.; Ye, L. Combining Immune-Related Genes For Delineating the Extracellular Matrix and Predicting Hormone Therapy and Neoadjuvant Chemotherapy Benefits In Breast Cancer. Front. Immunol. 2022, 13, 888339. [Google Scholar] [CrossRef] [PubMed]
  142. Frost, H.R.; Amos, C.I. Gene set selection via LASSO penalized regression (SLPR). Nucleic Acids Res. 2017, 45, e114. [Google Scholar] [CrossRef] [PubMed]
Figure 1. A workflow for top-down and bottoms-up system biology (SB) approaches to address prostate cancer health disparities. For top-down SB there is a plethora of high dimensional data including germline structural variants, population characteristics as serum-borne, and other clinical data that can be integrated through a range of approaches to deliver insights into how the prostate could either be either sensitive or resistant to the growth restraint properties of the VDR. From the bottoms-up perspective, more focused kinetic models can be developed for how the VDR functions and is disrupted between normal and tumor cells. Combining these approaches has the potential to develop a holistic understanding of how the VDR functions in the prostate gland and how this is impacted by genomic ancestry.
Figure 1. A workflow for top-down and bottoms-up system biology (SB) approaches to address prostate cancer health disparities. For top-down SB there is a plethora of high dimensional data including germline structural variants, population characteristics as serum-borne, and other clinical data that can be integrated through a range of approaches to deliver insights into how the prostate could either be either sensitive or resistant to the growth restraint properties of the VDR. From the bottoms-up perspective, more focused kinetic models can be developed for how the VDR functions and is disrupted between normal and tumor cells. Combining these approaches has the potential to develop a holistic understanding of how the VDR functions in the prostate gland and how this is impacted by genomic ancestry.
Nutrients 14 05197 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hussain, S.; Yates, C.; Campbell, M.J. Vitamin D and Systems Biology. Nutrients 2022, 14, 5197. https://doi.org/10.3390/nu14245197

AMA Style

Hussain S, Yates C, Campbell MJ. Vitamin D and Systems Biology. Nutrients. 2022; 14(24):5197. https://doi.org/10.3390/nu14245197

Chicago/Turabian Style

Hussain, Shahid, Clayton Yates, and Moray J. Campbell. 2022. "Vitamin D and Systems Biology" Nutrients 14, no. 24: 5197. https://doi.org/10.3390/nu14245197

APA Style

Hussain, S., Yates, C., & Campbell, M. J. (2022). Vitamin D and Systems Biology. Nutrients, 14(24), 5197. https://doi.org/10.3390/nu14245197

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop