Understanding human genome regulation through entropic graph-based regularization and submodular optimization

Thesis (Ph.D.)--University of Washington, 2016-06

Oxygen-regulated gene expression in bovine blastocysts

Oxygen concentrations used during in vitro embryo culture can influence embryo development, cell numbers, and gene expression. Here we propose that the preimplantation bovine embryo possesses a molecular mechanism for the detection of, and response to, oxygen, mediated by a family of basic helix-loop-helix transcription factors, the hypoxia-inducible factors (HIFs). Day 5 compacting bovine embryos were cultured under different oxygen tensions (2%, 7%, 20%) and the effect on the expression of oxygen-regulated genes, development, and cell number allocation and HIFalpha protein localization were examined. Bovine in vitro-produced embryos responded to variations in oxygen concentration by altering gene expression. GLUT1 expression was higher following 2% oxygen culture compared with 7% and 20% cultured blastocysts. HIF mRNA expression (HIF1alpha, HIF2alpha) was unaltered by oxygen concentration. HIF2alpha protein was predominantly localized to the nucleus of blastocysts. In contrast, HIF1alpha protein was undetectable at any oxygen concentration or in the presence of the HIF protein stabilizer desferrioxamine (DFO), despite being detectable in cumulus cells following normal maturation conditions, acute anoxic culture, or in the presence of DFO. Oxygen concentration also significantly altered inner cell mass cell proportions at the blastocyst stage. These results suggest that oxygen can influence gene expression in the bovine embryo during postcompaction development and that these effects may be mediated by HIF2alpha.

Inherent size constraints on prokaryote gene networks due to "accelerating" growth

Networks exhibiting accelerating growth have total link numbers growing faster than linearly with network size and either reach a limit or exhibit graduated transitions from nonstationary-to-stationary statistics and from random to scale-free to regular statistics as the network size grows. However, if for any reason the network cannot tolerate such gross structural changes then accelerating networks are constrained to have sizes below some critical value. This is of interest as the regulatory gene networks of single-celled prokaryotes are characterized by an accelerating quadratic growth and are size constrained to be less than about 10,000 genes encoded in DNA sequence of less than about 10 megabases. This paper presents a probabilistic accelerating network model for prokaryotic gene regulation which closely matches observed statistics by employing two classes of network nodes (regulatory and non-regulatory) and directed links whose inbound heads are exponentially distributed over all nodes and whose outbound tails are preferentially attached to regulatory nodes and described by a scale-free distribution. This model explains the observed quadratic growth in regulator number with gene number and predicts an upper prokaryote size limit closely approximating the observed value. (c) 2005 Elsevier GmbH. All rights reserved.

Structure and dynamics of a gene network model incorporating small RNAs

As advances in molecular biology continue to reveal additional layers of complexity in gene regulation, computational models need to incorporate additional features to explore the implications of new theories and hypotheses. It has recently been suggested that eukaryotic organisms owe their phenotypic complexity and diversity to the exploitation of small RNAs as signalling molecules. Previous models of genetic systems are, for several reasons, inadequate to investigate this theory. In this study, we present an artificial genome model of genetic regulatory networks based upon previous work by Torsten Reil, and demonstrate how this model generates networks with biologically plausible structural and dynamic properties. We also extend the model to explore the implications of incorporating regulation by small RNA molecules in a gene network. We demonstrate how, using these signals, highly connected networks can display dynamics that are more stable than expected given their level of connectivity.

Species differential regulation of COX2 can be described by an NFκB-dependent logic AND gate

Cyclooxygenase 2 (COX2), a key regulatory enzyme of the prostaglandin/eicosanoid pathway, is an important target for anti-inflammatory therapy. It is highly induced by pro-inflammatory cytokines in a Nuclear factor kappa B (NFκB)-dependent manner. However, the mechanisms determining the amplitude and dynamics of this important pro-inflammatory event are poorly understood. Furthermore, there is significant difference between human and mouse COX2 expression in response to the inflammatory stimulus tumor necrosis factor alpha (TNFα). Here, we report the presence of a molecular logic AND gate composed of two NFκB response elements (NREs) which controls the expression of human COX2 in a switch-like manner. Combining quantitative kinetic modeling and thermostatistical analysis followed by experimental validation in iterative cycles, we show that the human COX2 expression machinery regulated by NFκB displays features of a logic AND gate. We propose that this provides a digital, noise-filtering mechanism for a tighter control of expression in response to TNFα, such that a threshold level of NFκB activation is required before the promoter becomes active and initiates transcription. This NFκB-regulated AND gate is absent in the mouse COX2 promoter, most likely contributing to its differential graded response in promoter activity and protein expression to TNFα. Our data suggest that the NFκB-regulated AND gate acts as a novel mechanism for controlling the expression of human COX2 to TNFα, and its absence in the mouse COX2 provides the foundation for further studies on understanding species-specific differential gene regulation.

Targeted Gene Repression Technologies for Regenerative Medicine, Genomics, and Gene Therapy

Gene regulation is a complex and tightly controlled process that defines cell function in physiological and abnormal states. Programmable gene repression technologies enable loss-of-function studies for dissecting gene regulation mechanisms and represent an exciting avenue for gene therapy. Established and recently developed methods now exist to modulate gene sequence, epigenetic marks, transcriptional activity, and post-transcriptional processes, providing unprecedented genetic control over cell phenotype. Our objective was to apply and develop targeted repression technologies for regenerative medicine, genomics, and gene therapy applications. We used RNA interference to control cell cycle regulation in myogenic differentiation and enhance the proliferative capacity of tissue engineered cartilage constructs. These studies demonstrate how modulation of a single gene can be used to guide cell differentiation for regenerative medicine strategies. RNA-guided gene regulation with the CRISPR/Cas9 system has rapidly expanded the targeted repression repertoire from silencing single protein-coding genes to modulation of genes, promoters, and other distal regulatory elements. In order to facilitate its adaptation for basic research and translational applications, we demonstrated the high degree of specificity for gene targeting, gene silencing, and chromatin modification possible with Cas9 repressors. The specificity and effectiveness of RNA-guided transcriptional repressors for silencing endogenous genes are promising characteristics for mechanistic studies of gene regulation and cell phenotype. Furthermore, our results support the use of Cas9-based repressors as a platform for novel gene therapy strategies. We developed an in vivo AAV-based gene repression system for silencing endogenous genes in a mouse model. Together, these studies demonstrate the utility of gene repression tools for guiding cell phenotype and the potential of the RNA-guided CRISPR/Cas9 platform for applications such as causal studies of gene regulatory mechanisms and gene therapy.

Roles of CTCF and YY1 in T Cell Receptor Gene Rearrangement And T Cell Development

Diversity of T cell receptors (TCR) and immunoglobulins (Ig) is generated by V(D)J recombination of antigen receptor (AgR) loci. The Tcra-Tcrd locus is of particular interest because it displays a nested organization of Tcrd and Tcra gene segments and V(D)J recombination follows an intricate developmental program to assemble both TCRδ and TCRα repertoires. However, the mechanisms that dictate the developmental regulation of V(D)J recombination of the Tcra-Tcrd locus remain unclear.

We have previously shown that CCCTC-binding factor (CTCF) regulates Tcra gene transcription and rearrangement through organizing chromatin looping between CTCF- binding elements (CBEs). This study is one of many showing that CTCF functions as a chromatin organizer and transcriptional regulator genome-wide. However, detailed understanding of the impact of specific CBEs is needed to fully comprehend the biological function of CTCF and how CTCF influences the generation of the TCR repertoire during thymocyte development. Thus, we generated several mouse models with genetically modified CBEs to gain insight into the CTCF-dependent regulation of the Tcra-Tcrd locus. We revealed a CTCF-dependent chromatin interaction network at the Tcra-Tcrd locus in double-negative thymocytes. Disruption of a discrete chromatin loop encompassing Dδ, Jδ and Cδ gene segments allowed a single Vδ segment to frequently contact and rearrange to diversity and joining gene segments and dominate the adult TCRδ repertoire. Disruption of this loop also narrowed the TCRα repertoire, which, we believe, followed as a consequence of the restricted TCRδ repertoire. Hence, a single CTCF-mediated chromatin loop directly regulates TCRδ diversity and indirectly regulates TCRα diversity. In addition, we showed that insertion of an ectopic CBE can modify chromatin interactions and disrupt the rearrangement of particular Vδ gene segments. Finally, we investigated the role of YY1 in early T cell development by conditionally deleting YY1 in developing thymocytes. We found that early ablation of YY1 caused severe developmental defects in the DN compartment due to a dramatic increase in DN thymocyte apoptosis. Furthermore, late ablation of YY1 resulted in increased apoptosis of DP thymocytes and a restricted TCRα repertoire. Mechanistically, we showed that p53 was upregulated in both DN and DP YY1-deficient thymocytes. Eliminating p53 in YY1-deficient thymocytes rescued the survival and developmental defects, indicating that these YY1-dependent defects were p53-mediated. We conclude that YY1 is required to maintain cell viability during thymocyte development by thwarting the accumulation of p53.

Overall, this thesis work has shown that CTCF-dependent looping provides a central framework for lineage- and developmental stage-specific regulation of Tcra-Tcrd gene expression and rearrangements. In addition, we identified YY1 as a novel regulator of thymocyte viability.

RNA Localization and Translational Regulation on the Endoplasmic Reticulum

mRNA localization is emerging as a critical cellular mechanism for the spatiotemporal regulation of protein expression and serves important roles in oogenesis, embryogenesis, cell fate specification, and synapse formation. Signal sequence-encoding mRNAs are localized to the endoplasmic reticulum (ER) membrane by either of two mechanisms, a canonical mechanism of translation on ER-bound ribosomes (signal recognition particle pathway), or a poorly understood direct ER anchoring mechanism. In this study, we identify that the ER integral membrane proteins function as RNA-binding proteins and play important roles in the direct mRNA anchoring to the ER. We report that one of the ER integral membrane RNA-binding protein, AEG-1 (astrocyte elevated gene-1), functions in the direct ER anchoring and translational regulation of mRNAs encoding endomembrane transmembrane proteins. HITS-CLIP and PAR-CLIP analyses of the AEG-1 mRNA interactome of human hepatocellular carcinoma cells revealed a high enrichment for mRNAs encoding endomembrane organelle proteins, most notably encoding transmembrane proteins. AEG-1 binding sites were highly enriched in the coding sequence and displayed a signature cluster enrichment downstream of encoded transmembrane domains. In overexpression and knockdown models, AEG-1 expression markedly regulates translational efficiency and protein functions of two of its bound transcripts, MDR1 and NPC1. This study reveals a molecular mechanism for the selective localization of mRNAs to the ER and identifies a novel post-transcriptional gene regulation function for AEG-1 in membrane protein expression.

Phase transitions in gene networks evolved under different selection rules

Mathematical models of gene regulation are a powerful tool for understanding the complex features of genetic control. While various modeling efforts have been successful at explaining gene expression dynamics, much less is known about how evolution shapes the structure of these networks. An important feature of gene regulatory networks is their stability in response to environmental perturbations. Regulatory systems are thought to have evolved to exist near the transition between stability and instability, in order to have the required stability to environmental fluctuations while also being able to achieve a wide variety of functions (corresponding to different dynamical patterns). We study a simplified model of gene network evolution in which links are added via different selection rules. These growth models are inspired by recent work on `explosive' percolation which shows that when network links are added through competitive rather than random processes, the connectivity phase transition can be significantly delayed, and when it is reached, it appears to be first order (discontinuous, e.g., going from no failure at all to large expected failure) instead of second order (continuous, e.g., going from no failure at all to very small expected failure). We find that by modifying the traditional framework for networks grown via competitive link addition to capture how gene networks evolve to avoid damage propagation, we also see significant delays in the transition that depend on the selection rules, but the transitions always appear continuous rather than `explosive'.

The genetic basis of human height : the role of estrogen

Height is a complex physical trait that displays strong heritability. Adult height is related to length of the long bones, which is determined by growth at the epiphyseal growth plate. Longitudinal bone growth occurs via the process of endochondral ossification, where bone forms over the differentiating cartilage template at the growth plate. Estrogen plays a major role in regulating longitudinal bone growth and is responsible for inducing the pubertal growth spurt and fusion of the epiphyseal growth plate. However, the mechanism by which estrogen promotes epiphyseal fusion is poorly understood. It has been hypothesised that estrogen functions to regulate growth plate fusion by stimulating chondrocyte apoptosis, angiogenesis and bone cell invasion in the growth plate. Another theory has suggested that estrogen exposure exhausts the proliferative capacity of growth plate chondrocytes, which accelerates the process of chondrocyte senescence, leading to growth plate fusion. The overall objective of this study was to gain a greater understanding of the molecular mechanisms behind estrogen-mediated growth and height attainment by examining gene regulation in chondrocytes and the role of some of these genes in normal height inheritance. With the heritability of height so well established, the initial hypothesis was that genetic variation in candidate genes associated with longitudinal bone growth would be involved in normal adult height variation. The height-related genes FGFR3, CBFA1, ER and CBFA1 were screened for novel polymorphisms using denaturing HPLC and RFLP analysis. In total, 24 polymorphisms were identified. Two SNPs in ER (rs3757323 C>T and rs1801132 G>C) were strongly associated with adult male height and displayed an 8 cm and 9 cm height difference between homozygous genotypes, respectively. The TC haplotype of these SNPs was associated with a 6 cm decrease in height and remarkably, no homozygous carriers of the TC haplotype were identified in tall subjects. No significant associations with height were found for polymorphisms in the FGFR3, CBFA1 or VDR genes. In the epiphyseal growth plate, chondrocyte proliferation, matrix synthesis and chondrocyte hypertrophy are all major contributors to long bone growth. As estrogen plays such a significant role in both growth and final height attainment, another hypothesis of this study was that estrogen exerted its effects in the growth plate by influencing chondrocyte proliferation and mediating the expression of chondrocyte marker genes. The examination of genes regulated by estrogen in chondrocyte-like cells aimed to identify potential regulators of growth plate fusion, which may further elucidate mechanisms involved in the cessation of linear growth. While estrogen did not dramatically alter the proliferation of the SW1353 cell line, gene expression experiments identified several estrogen regulated genes. Sixteen chondrocyte marker genes were examined in response to estrogen concentrations ranging from 10-12 M to 10-8 M over varying time points. Of the genes analysed, IHH, FGFR3, collagen II and collagen X were not readily detectable and PTHrP, GHR, ER, BMP6, SOX9 and TGF1 mRNAs showed no significant response to estrogen treatments. However, the expression of MMP13, CBFA1, BCL-2 and BAX genes were significantly decreased. Interestingly, the majority of estrogen regulated genes in SW1353 cells are expressed in the hypertrophic zone of the growth plate. Estrogen is also known to regulate systemic GH secretion and local GH action. At the molecular level, estrogen functions to inhibit GH action by negatively regulating GH signalling. GH treated SW1353 cells displayed increases in MMP9 mRNA expression (4.4-fold) and MMP13 mRNA expression (64-fold) in SW1353 cells. Increases were also detected in their respective proteins. Treatment with AG490, an established JAK2 inhibitor, blocked the GH mediated stimulation of both MMP9 and MMP13 mRNA expression. The application of estrogen and GH to SW1353 cells attenuated GH-stimulated MMP13 levels, but did not affect MMP9 levels. Investigation of GH signalling revealed that SW1353 cells have high levels of activated JAK2 and exposure to GH, estrogen, AG490 and other signalling inhibitors did not affect JAK2 phosphorylation. Interestingly, AG490 treatment dramatically decreased ERK2 signalling, although GH did stimulate ERK2 phosphorylation above control levels. AG490 also decreased CBFA1 expression, a transcription factor known to activate MMP9 and MMP13. Finally, GH and estrogen treatment increased expression of SOCS3 mRNA, suggesting that SOCS3 may regulate JAK/STAT signalling in SW1353 cells. The modulation of GH-mediated MMP expression by estrogen in SW1353 cells represents a potentially novel mechanism by which estrogen may regulate longitudinal bone growth. However, further investigation is required in order to elucidate the precise mechanisms behind estrogen and GH regulation of MMP13 expression in SW1353 cells. This study has provided additional evidence that estrogen and the ER gene are major factors in the regulation of growth and the determination of adult height. Newly identified polymorphisms in the ER gene not only contribute to our understanding of the genetic basis of human height, but may also be useful in association studies examining other complex traits. This study also identified several estrogen regulated genes and indicated that estrogen modifies the expression of genes which are primarily expressed in the hypertrophic region of the epiphyseal growth plate. Furthermore, synergistic studies incorporating GH and estrogen have revealed the ability of estrogen to attenuate the effects of GH on MMP13 expression, revealing potential pathways by which estrogen may modulate growth plate fusion, longitudinal bone growth and even arthritis.

A novel corepressor, BCoR-L1, represses transcription through an interaction with CtBP

Corepressors play a crucial role in negative gene regulation and are defective in several diseases. BCoR is a corepressor for the BCL6 repressor protein. Here we describe and functionally characterize BCoR-L1, a homolog of BCoR. When tethered to a heterologous promoter, BCoR-L1 is capable of strong repression. Like other corepressors, BCoR-L1 associates with histone deacetylase (HDAC) activity. Specifically, BCoR-L1 coprecipitates with the Class II HDACs, HDAC4, HDAC5, and HDAC7, suggesting that they are involved in its role as a transcriptional repressor. BCoR-L1 also interacts with the CtBP corepressor through a CtBP-interacting motif in its amino terminus. Abrogation of the CtBP binding site within BCoR-L1 partially relieves BCoR-L1-mediated transcriptional repression. Furthermore, BCoR-L1 is located on the E-cadherin promoter, a known CtBP-regulated promoter, and represses the E-cadherin promoter activity in a reporter assay. The inhibition of BCoR-L1 expression by RNA-mediated interference results in derepression of E-cadherin in cells that do not normally express E-cadherin, indicating that BCoR-L1 contributes to the repression of an authentic endogenous CtBP target.

Myocyte enhancer factor 2C : an osteoblast transcription factor identified by DMSO enhanced mineralization

Rapid mineralization of cultured osteoblasts could be a useful characteristic in stem-cell mediated therapies for fracture and other orthopaedic problems. Dimethyl sulfoxide (DMSO) is a small amphipathic solvent molecule capable of simulating cell differentiation. We report that, in primary human osteoblasts, DMSO dose-dependently enhanced the expression of osteoblast differentiation markers alkaline phosphatase (ALP) activity and extracellular matrix mineralization. Furthermore, similar DMSO mediated mineralization enhancement was observed in primary osteoblast-like cells differentiated from mouse mesenchymal cells derived from fat, a promising source of starter cells for cell-based therapy. Using a convenient mouse pre-osteoblast model cell line MC3T3-E1 we further investigated this phenomenon showing that numerous osteoblast-expressed genes were elevated in response to DMSO treatment and correlated with enhanced mineralization. Myocyte enhancer factor 2c (Mef2c) was identified as the transcription factor most induced by DMSO, among numerous DMSO-induced genes, suggesting a role for Mef2c in osteoblast gene regulation. Immunohistochemistry confirmed expression of Mef2c in osteoblast-like cells in mouse mandible, cortical and trabecular bone. shRNAi-mediated Mef2c gene silencing resulted in defective osteoblast differentiation, decreased ALP activity and matrix mineralization and knockdown of osteoblast specific gene expression, including osteocalcin and bone sialoprotein. Flow on knockdown of bone specific transcription factors, Runx2 and osterix by shRNAi knockdown of Mef2c suggests that Mef2c lies upstream of these two important factors in the cascade of gene expression in osteoblasts.

Stochastic simulation for spatial modelling of dynamic process in a living cell

One of the fundamental motivations underlying computational cell biology is to gain insight into the complicated dynamical processes taking place, for example, on the plasma membrane or in the cytosol of a cell. These processes are often so complicated that purely temporal mathematical models cannot adequately capture the complex chemical kinetics and transport processes of, for example, proteins or vesicles. On the other hand, spatial models such as Monte Carlo approaches can have very large computational overheads. This chapter gives an overview of the state of the art in the development of stochastic simulation techniques for the spatial modelling of dynamic processes in a living cell.

Characterisation of the TaNF-Y family of transcription factors in wheat

Light plays a unique role for plants as it is both a source of energy for growth and a signal for development. Light captured by the pigments in the light harvesting complexes is used to drive the synthesis of the chemical energy required for carbon assimilation. The light perceived by photoreceptors activates effectors, such as transcription factors (TFs), which modulate the expression of light-responsive genes. Recently, it has been speculated that increasing the photosynthetic rate could further improve the yield potential of three carbon (C3) crops such as wheat. However, little is currently known about the transcriptional regulation of photosynthesis genes, particularly in crop species. Nuclear factor Y (NF-Y) TF is a functionally diverse regulator of growth and development in the model plant species, with demonstrated roles in embryo development, stress response, flowering time and chloroplast biogenesis. Furthermore, a light-responsive NF-Y binding site (CCAAT-box) is present in the promoter of a spinach photosynthesis gene. As photosynthesis genes are co-regulated by light and co-regulated genes typically have similar regulatory elements in their promoters, it seems likely that other photosynthesis genes would also have light-responsive CCAAT-boxes. This provided the impetus to investigate the NF-Y TF in bread wheat. This thesis is focussed on wheat NF-Y members that have roles in light-mediated gene regulation with an emphasis on their involvement in the regulation of photosynthesis genes. NF-Y is a heterotrimeric complex, comprised of the three subunits NF-YA, NF-YB and NF-YC. Unlike the mammalian and yeast counterparts, each of the three subunits is encoded by multiple genes in Arabidopsis. The initial step taken in this study was the identification of the wheat NF-Y family (Chapter 3). A search of the current wheat nucleotide sequence databases identified 37 NF-Y genes (10 NF-YA, 11 NF-YB, 14 NF-YC & 2 Dr1). Phylogenetic analysis revealed that each of the three wheat NF-Y (TaNF-Y) subunit families could be divided into 4-5 clades based on their conserved core regions. Outside of the core regions, eleven motifs were identified to be conserved between Arabidopsis, rice and wheat NF-Y subunit members. The expression profiles of TaNF-Y genes were constructed using quantitative real-time polymerase chain reaction (RT-PCR). Some TaNF-Y subunit members had little variation in their transcript levels among the organs, while others displayed organ-predominant expression profiles, including those expressed mainly in the photosynthetic organs. To investigate their potential role in light-mediated gene regulation, the light responsiveness of the TaNF-Y genes were examined (Chapters 4 and 5). Two TaNF-YB and five TaNF-YC members were markedly upregulated by light in both the wheat leaves and seedling shoots. To identify the potential target genes of the light-upregulated NF-Y subunit members, a gene expression correlation analysis was conducted using publically available Affymetrix Wheat Genome Array datasets. This analysis revealed that the transcript expression levels of TaNF-YB3 and TaNF-YC11 were significantly correlated with those of photosynthesis genes. These correlated express profiles were also observed in the quantitative RT-PCR dataset from wheat plants grown under light and dark conditions. Sequence analysis of the promoters of these wheat photosynthesis genes revealed that they were enriched with potential NF-Y binding sites (CCAAT-box). The potential role of TaNF-YB3 in the regulation of photosynthetic genes was further investigated using a transgenic approach (Chapter 5). Transgenic wheat lines constitutively expressing TaNF-YB3 were found to have significantly increased expression levels of photosynthesis genes, including those encoding light harvesting chlorophyll a/b-binding proteins, photosystem I reaction centre subunits, a chloroplast ATP synthase subunit and glutamyl-tRNA reductase (GluTR). GluTR is a rate-limiting enzyme in the chlorophyll biosynthesis pathway. In association with the increased expression of the photosynthesis genes, the transgenic lines had a higher leaf chlorophyll content, increased photosynthetic rate and had a more rapid early growth rate compared to the wild-type wheat. In addition to its role in the regulation of photosynthesis genes, TaNF-YB3 overexpression lines flower on average 2-days earlier than the wild-type (Chapter 6). Quantitative RT-PCR analysis showed that there was a 13-fold increase in the expression level of the floral integrator, TaFT. The transcript levels of other downstream genes (TaFT2 and TaVRN1) were also increased in the transgenic lines. Furthermore, the transcript levels of TaNF-YB3 were significantly correlated with those of constans (CO), constans-like (COL) and timing of chlorophyll a/b-binding (CAB) expression 1 [TOC1; (CCT)] domain-containing proteins known to be involved in the regulation of flowering time. To summarise the key findings of this study, 37 NF-Y genes were identified in the crop species wheat. An in depth analysis of TaNF-Y gene expression profiles revealed that the potential role of some light-upregulated members was in the regulation of photosynthetic genes. The involvement of TaNF-YB3 in the regulation of photosynthesis genes was supported by data obtained from transgenic wheat lines with increased constitutive expression of TaNF-YB3. The overexpression of TaNF-YB3 in the transgenic lines revealed this NF-YB member is also involved in the fine-tuning of flowering time. These data suggest that the NF-Y TF plays an important role in light-mediated gene regulation in wheat.