Variants on the UBE2L3/YDJC Autoimmune Disease Risk Haplotype Increase UBE2L3 Expression by Modulating CCCTC-Binding Factor and YY1 Binding.

OBJECTIVE: Genetic variants spanning UBE2L3 are associated with increased expression of the UBE2L3-encoded E2 ubiquitin-conjugating enzyme H7 (UbcH7), which facilitates activation of proinflammatory NF-κB signaling and susceptibility to autoimmune diseases. We undertook this study to delineate how genetic variants carried on the UBE2L3/YDJC autoimmune risk haplotype function to drive hypermorphic UBE2L3 expression. METHODS: We used bioinformatic analyses, electrophoretic mobility shift assays, and luciferase reporter assays to identify and functionally characterize allele-specific effects of risk variants positioned in chromatin accessible regions of immune cells. Chromatin conformation capture with quantitative polymerase chain reaction (3C-qPCR), chromatin immunoprecipitation (ChIP)-qPCR, and small interfering RNA (siRNA) knockdown assays were performed on patient-derived Epstein-Barr virus-transformed B cells homozygous for the UBE2L3/YDJC nonrisk or risk haplotype to determine if the risk haplotype increases UBE2L3 expression by altering the regulatory chromatin architecture in the region. RESULTS: Of the 7 prioritized variants, 5 demonstrated allele-specific increases in nuclear protein binding affinity and regulatory activity. High-throughput sequencing of chromosome conformation capture coupled with ChIP (HiChIP) and 3C-qPCR uncovered a long-range interaction between the UBE2L3 promoter (rs140490, rs140491, rs11089620) and the downstream YDJC promoter (rs3747093) that was strengthened in the presence of the UBE2L3/YDJC risk haplotype, and correlated with the loss of CCCTC-binding factor (CTCF) and gain of YY1 binding at the risk alleles. Depleting YY1 by siRNA disrupted the long-range interaction between the 2 promoters and reduced UBE2L3 expression. CONCLUSION: The UBE2L3/YDJC autoimmune risk haplotype increases UBE2L3 expression through strengthening a YY1-mediated interaction between the UBE2L3 and YDJC promoters.

CTCF chromatin residence time controls three-dimensional genome organization, gene expression and DNA methylation in pluripotent cells.

The 11 zinc finger (ZF) protein CTCF regulates topologically associating domain formation and transcription through selective binding to thousands of genomic sites. Here, we replaced endogenous CTCF in mouse embryonic stem cells with green-fluorescent-protein-tagged wild-type or mutant proteins lacking individual ZFs to identify additional determinants of CTCF positioning and function. While ZF1 and ZF8-ZF11 are not essential for cell survival, ZF8 deletion strikingly increases the DNA binding off-rate of mutant CTCF, resulting in reduced CTCF chromatin residence time. Loss of ZF8 results in widespread weakening of topologically associating domains, aberrant gene expression and increased genome-wide DNA methylation. Thus, important chromatin-templated processes rely on accurate CTCF chromatin residence time, which we propose depends on local sequence and chromatin context as well as global CTCF protein concentration.

Drosophila architectural protein CTCF is not essential for fly survival and is able to function independently of CP190.

CTCF is the most likely ancestor of proteins that contain large clusters of C2H2 zinc finger domains (C2H2) and is conserved among most bilateral organisms. In mammals, CTCF functions as the main architectural protein involved in the organization of topology-associated domains (TADs). In vertebrates and Drosophila, CTCF is involved in the regulation of homeotic genes. Previously, it was found that null mutations in the dCTCF gene died as pharate adults, which failed to eclose from their pupal case, or shortly after hatching of adults. Here, we obtained several new null dCTCF mutations and found that the complete inactivation of dCTCF appears is limited mainly to phenotypic manifestations of the Abd-B gene and fertility of adult flies. Many modifiers that are not associated with an independent phenotypic manifestation can significantly enhance the expressivity of the null dCTCF mutations, indicating that other architectural proteins are able to functionally compensate for dCTCF inactivation in Drosophila. We also mapped the 715-735 aa region of dCTCF as being essential for the interaction with the BTB (Broad-Complex, Tramtrack, and Bric a brac) and microtubule-targeting (M) domains of the CP190 protein, which binds to many architectural proteins. However, the mutational analysis showed that the interaction with CP190 was not important for the functional activity of dCTCF in vivo.

The role of CTCF in the organization of the centromeric 11p15 imprinted domain interactome.

DNA methylation, chromatin-binding proteins, and DNA looping are common components regulating genomic imprinting which leads to parent-specific monoallelic gene expression. Loss of methylation (LOM) at the human imprinting center 2 (IC2) on chromosome 11p15 is the most common cause of the imprinting overgrowth disorder Beckwith-Wiedemann Syndrome (BWS). Here, we report a familial transmission of a 7.6 kB deletion that ablates the core promoter of KCNQ1. This structural alteration leads to IC2 LOM and causes recurrent BWS. We find that occupancy of the chromatin organizer CTCF is disrupted proximal to the deletion, which causes chromatin architecture changes both in cis and in trans. We also profile the chromatin architecture of IC2 in patients with sporadic BWS caused by isolated LOM to identify conserved features of IC2 regulatory disruption. A strong interaction between CTCF sites around KCNQ1 and CDKN1C likely drive their expression on the maternal allele, while a weaker interaction involving the imprinting control region element may impede this connection and mediate gene silencing on the paternal allele. We present an imprinting model in which KCNQ1 transcription is necessary for appropriate CTCF binding and a novel chromatin conformation to drive allele-specific gene expression.

Forebrain glutamatergic neuron-specific Ctcf deletion induces reactive microgliosis and astrogliosis with neuronal loss in adult mouse hippocampus.

CCCTC-binding factor (CTCF), a zinc finger protein, is a transcription factor and regulator of chromatin structure. Forebrain excitatory neuron-specific CTCF deficiency contributes to inflammation via enhanced transcription of inflammation-related genes in the cortex and hippocampus. However, little is known about the long-term effect of CTCF deficiency on postnatal neurons, astrocytes, or microglia in the hippocampus of adult mice. To address this, we knocked out the Ctcf gene in forebrain glutamatergic neurons (Ctcf cKO) by crossing Ctcf-floxed mice with Camk2a-Cre mice and examined the hippocampi of 7.5-10-month-old male mice using immunofluorescence microscopy. We found obvious neuronal cell death and reactive gliosis in the hippocampal cornu ammonis (CA)1 in 7.5-10-month-old cKO mice. Prominent rod-shaped microglia that participate in immune surveillance were observed in the stratum pyramidale and radiatum layer, indicating a potential increase in inflammatory mediators released by hippocampal neurons. Although neuronal loss was not observed in CA3, and dentate gyrus (DG) CTCF depletion induced a significant increase in the number of microglia in the stratum oriens of CA3 and reactive microgliosis and astrogliosis in the molecular layer and hilus of the DG in 7.5-10-month-old cKO mice. These results suggest that long-term Ctcf deletion from forebrain excitatory neurons may contribute to reactive gliosis induced by neuronal damage and consequent neuronal loss in the hippocampal CA1, DG, and CA3 in sequence over 7 months of age. [BMB Reports 2021; 54(6): 317-322].

CTCF-binding element regulates ESC differentiation via orchestrating long-range chromatin interaction between enhancers and HoxA.

Proper expression of Homeobox A cluster genes (HoxA) is essential for embryonic stem cell (ESC) differentiation and individual development. However, mechanisms controlling precise spatiotemporal expression of HoxA during early ESC differentiation remain poorly understood. Herein, we identified a functional CTCF-binding element (CBE+47) closest to the 3'-end of HoxA within the same topologically associated domain (TAD) in ESC. CRISPR-Cas9-mediated deletion of CBE+47 significantly upregulated HoxA expression and enhanced early ESC differentiation induced by retinoic acid (RA) relative to wild-type cells. Mechanistic analysis by chromosome conformation capture assay (Capture-C) revealed that CBE+47 deletion decreased interactions between adjacent enhancers, enabling formation of a relatively loose enhancer-enhancer interaction complex (EEIC), which overall increased interactions between that EEIC and central regions of HoxA chromatin. These findings indicate that CBE+47 organizes chromatin interactions between its adjacent enhancers and HoxA. Furthermore, deletion of those adjacent enhancers synergistically inhibited HoxA activation, suggesting that these enhancers serve as an EEIC required for RA-induced HoxA activation. Collectively, these results provide new insight into RA-induced HoxA expression during early ESC differentiation, also highlight precise regulatory roles of the CTCF-binding element in orchestrating high-order chromatin structure.

The CCCTC-binding factor CTCF represses hepatitis B virus enhancer I and regulates viral transcription.

Hepatitis B virus (HBV) infection is of global importance with over 2 billion people exposed to the virus during their lifetime and at risk of progressive liver disease, cirrhosis and hepatocellular carcinoma. HBV is a member of the Hepadnaviridae family that replicates via episomal copies of a covalently closed circular DNA (cccDNA) genome. The chromatinization of this small viral genome, with overlapping open reading frames and regulatory elements, suggests an important role for epigenetic pathways to regulate viral transcription. The chromatin-organising transcriptional insulator protein, CCCTC-binding factor (CTCF), has been reported to regulate transcription in a diverse range of viruses. We identified two conserved CTCF binding sites in the HBV genome within enhancer I and chromatin immunoprecipitation (ChIP) analysis demonstrated an enrichment of CTCF binding to integrated or episomal copies of the viral genome. siRNA knock-down of CTCF results in a significant increase in pre-genomic RNA levels in de novo infected HepG2 cells and those supporting episomal HBV DNA replication. Furthermore, mutation of these sites in HBV DNA minicircles abrogated CTCF binding and increased pre-genomic RNA levels, providing evidence of a direct role for CTCF in repressing HBV transcription.

CTCF-mediated genome organization and leukemogenesis.

Recent technological advancements and genome-wide studies provide compelling evidence that dynamic chromatin interaction and three-dimensional genome organization in nuclei play an important role in regulating gene expression. Mammalian genomes consist of many small functional domains termed topologically associated domains (TADs), many of them organized by CCCTC-binding factor (CTCF) and the cohesion complex. Changes in genome TADs might result in inappropriate promoter/enhancer communications leading to activation of oncogenes or suppression of tumor suppressors. During normal hematopoiesis and leukemogenesis, genome structure alters considerably to facilitate normal and malignant hematopoiesis, respectively. Delineating theses normal and abnormal processes will evolve our understanding of disease pathogenesis and development of potential treatment strategies. This review highlights the role of CTCF and its associated protein complexes in three-dimensional genome organization in development and leukemogenesis, as well as the roles of CTCF boundary defined TAD in transcription regulation. We further explore the function of chromatin modulators, such as CTCF, cohesin, and long noncoding RNAs (lncRNAs) in chromosomal interactions and hematopoietic genome organization. Finally, we focus on the implication of 3D genome alteration in the pathogenesis of leukemia and provide a scientific basis for targeted intervention.

CTCF is dispensable for immune cell transdifferentiation but facilitates an acute inflammatory response.

Three-dimensional organization of the genome is important for transcriptional regulation1-7. In mammals, CTCF and the cohesin complex create submegabase structures with elevated internal chromatin contact frequencies, called topologically associating domains (TADs)8-12. Although TADs can contribute to transcriptional regulation, ablation of TAD organization by disrupting CTCF or the cohesin complex causes modest gene expression changes13-16. In contrast, CTCF is required for cell cycle regulation17, embryonic development and formation of various adult cell types18. To uncouple the role of CTCF in cell-state transitions and cell proliferation, we studied the effect of CTCF depletion during the conversion of human leukemic B cells into macrophages with minimal cell division. CTCF depletion disrupts TAD organization but not cell transdifferentiation. In contrast, CTCF depletion in induced macrophages impairs the full-blown upregulation of inflammatory genes after exposure to endotoxin. Our results demonstrate that CTCF-dependent genome topology is not strictly required for a functional cell-fate conversion but facilitates a rapid and efficient response to an external stimulus.

ESCO1 and CTCF enable formation of long chromatin loops by protecting cohesinSTAG1 from WAPL.

Eukaryotic genomes are folded into loops. It is thought that these are formed by cohesin complexes via extrusion, either until loop expansion is arrested by CTCF or until cohesin is removed from DNA by WAPL. Although WAPL limits cohesin's chromatin residence time to minutes, it has been reported that some loops exist for hours. How these loops can persist is unknown. We show that during G1-phase, mammalian cells contain acetylated cohesinSTAG1 which binds chromatin for hours, whereas cohesinSTAG2 binds chromatin for minutes. Our results indicate that CTCF and the acetyltransferase ESCO1 protect a subset of cohesinSTAG1 complexes from WAPL, thereby enable formation of long and presumably long-lived loops, and that ESCO1, like CTCF, contributes to boundary formation in chromatin looping. Our data are consistent with a model of nested loop extrusion, in which acetylated cohesinSTAG1 forms stable loops between CTCF sites, demarcating the boundaries of more transient cohesinSTAG2 extrusion activity.

Guided nuclear exploration increases CTCF target search efficiency.

The enormous size of mammalian genomes means that for a DNA-binding protein the number of nonspecific, off-target sites vastly exceeds the number of specific, cognate sites. How mammalian DNA-binding proteins overcome this challenge to efficiently locate their target sites is not known. Here, through live-cell single-molecule tracking, we show that CCCTC-binding factor, CTCF, is repeatedly trapped in small zones that likely correspond to CTCF clusters, in a manner that is largely dependent on an internal RNA-binding region (RBRi). We develop a new theoretical model called anisotropic diffusion through transient trapping in zones to explain CTCF dynamics. Functionally, transient RBRi-mediated trapping increases the efficiency of CTCF target search by ~2.5-fold. Overall, our results suggest a 'guided' mechanism where CTCF clusters concentrate diffusing CTCF proteins near cognate binding sites, thus increasing the local ON-rate. We suggest that local guiding may allow DNA-binding proteins to more efficiently locate their target sites.

CTCF inhibits endoplasmic reticulum stress and apoptosis in cardiomyocytes by upregulating RYR2 via inhibiting S100A1.

AIMS: Pediatric heart failure is a common cardiovascular disease in clinical pediatrics. CCCTC-binding factor (CTCF), a novel transcriptional repressor, was reported to participate in the occurrence of various cardiovascular diseases. The present study focuses on exploring the effects of CTCF on tunicamycin (TM)-induced endoplasmic reticulum (ER) stress, and investigating the underlying mechanisms. MATERIALS AND METHOD: Expression of CTCF in blood samples of heart failure children and TM-induced cardiomyocytes were evaluated by real-time quantitative PCR (RT-qPCR). Apoptotic rate of cardiomyocytes was detected by Annexin v assay. Western blotting and enzyme-linked immunosorbent assay (ELISA) were applied to examine the effect of CTCF on ER stress. Co-immunoprecipitation and western blotting were devoted to explore the mechanism by which CTCF contributes to ER stress. KEY FINDINGS: We proved that CTCF was lowly expressed in blood samples of heart failure children and TM-induced cardiomyocytes, and overexpression of CTCF weaken the TM-induced ER stress. Using co-immunoprecipitation and protein blots, we demonstrated that CTCF upregulates RYR2 by inhibiting S100A1, thus mediating the PERK signaling pathway and regulating ER stress. SIGNIFICANCE: Our data revealed that CTCF protects cardiomyocytes from ER stress through S100A1-RYR2 axis, and can be applied as a therapeutic target for the treatment of pediatric heart failure in future.

CTCF-dependent chromatin boundaries formed by asymmetric nucleosome arrays with decreased linker length.

The CCCTC-binding factor (CTCF) organises the genome in 3D through DNA loops and in 1D by setting boundaries isolating different chromatin states, but these processes are not well understood. Here we investigate chromatin boundaries in mouse embryonic stem cells, defined by the regions with decreased Nucleosome Repeat Length (NRL) for ∼20 nucleosomes near CTCF sites, affecting up to 10% of the genome. We found that the nucleosome-depleted region (NDR) near CTCF is asymmetrically located >40 nucleotides 5'-upstream from the centre of CTCF motif. The strength of CTCF binding to DNA and the presence of cohesin is correlated with the decrease of NRL near CTCF, and anti-correlated with the level of asymmetry of the nucleosome array. Individual chromatin remodellers have different contributions, with Snf2h having the strongest effect on the NRL decrease near CTCF and Chd4 playing a major role in the symmetry breaking. Upon differentiation, a subset of preserved, common CTCF sites maintains asymmetric nucleosome pattern and small NRL. The sites which lost CTCF upon differentiation are characterized by nucleosome rearrangement 3'-downstream, with unchanged NDR 5'-upstream of CTCF motifs. Boundaries of topologically associated chromatin domains frequently contain several inward-oriented CTCF motifs whose effects, described above, add up synergistically.

LINC00346 promotes pancreatic cancer progression through the CTCF-mediated Myc transcription.

Although multiple factors are known to contribute to pancreatic ductal adenocarcinoma (PDAC) progression, the role of long non-coding RNAs (lncRNAs) in PDAC remains largely unknown. In this study, we present data that long intergenic non-coding RNA 346 (LINC00346) functions as a promoting factor for PDAC development. We first show that LINC00346 is highly expressed in pancreatic tumor specimens as compared to normal pancreatic tissue based on interrogation of The Cancer Genome Atlas (TCGA) pancreatic adenocarcinoma dataset. Of significance, this upregulation of LINC00346 is associated with overall survival (OS) and disease-free survival (DFS), respectively. We further show that knockout (KO) of LINC00346 impairs pancreatic cancer cell proliferation, tumorigenesis, migration, and invasion ability. Importantly, these phenotypes can be restored by LINC00346 re-expression in KO cells (i.e., rescue experiment). RNA precipitation assays combined with mass spectrometry analysis indicate that LINC00346 interacts with CCCTC-binding factor (CTCF), a known transcriptional repressor of c-Myc. This interaction between LINC00346 and CTCF prevents the binding of CTCF to c-Myc promoter, relieving the CTCF-mediated repression of c-Myc. Thus, LINC00346 functions as a positive transcriptional regulator of c-Myc. Together, these results suggest that LINC00346 contributes to PDAC pathogenesis by activating c-Myc, and as such, LINC00346 may serve as a potential biomarker and therapeutic target for PDAC.

Acute depletion of CTCF directly affects MYC regulation through loss of enhancer-promoter looping.

Numerous pieces of evidence support the complex, 3D spatial organization of the genome dictates gene expression. CTCF is essential to define topologically associated domain boundaries and to facilitate the formation of insulated chromatin loop structures. To understand CTCF's direct role in global transcriptional regulation, we integrated the miniAID-mClover3 cassette to the endogenous CTCF locus in a human pediatric B-ALL cell line, SEM, and an immortal erythroid precursor cell line, HUDEP-2, to allow for acute depletion of CTCF protein by the auxin-inducible degron system. In SEM cells, CTCF loss notably disrupted intra-TAD loops and TAD integrity in concurrence with a reduction in CTCF-binding affinity, while showing no perturbation to nuclear compartment integrity. Strikingly, the overall effect of CTCF's loss on transcription was minimal. Whole transcriptome analysis showed hundreds of genes differentially expressed in CTCF-depleted cells, among which MYC and a number of MYC target genes were specifically downregulated. Mechanically, acute depletion of CTCF disrupted the direct interaction between the MYC promoter and its distal enhancer cluster residing ∼1.8 Mb downstream. Notably, MYC expression was not profoundly affected upon CTCF loss in HUDEP-2 cells suggesting that CTCF could play a B-ALL cell line specific role in maintaining MYC expression.

An alternative CTCF isoform antagonizes canonical CTCF occupancy and changes chromatin architecture to promote apoptosis.

CTCF plays key roles in gene regulation, chromatin insulation, imprinting, X chromosome inactivation and organizing the higher-order chromatin architecture of mammalian genomes. Previous studies have mainly focused on the roles of the canonical CTCF isoform. Here, we explore the functions of an alternatively spliced human CTCF isoform in which exons 3 and 4 are skipped, producing a shorter isoform (CTCF-s). Functionally, we find that CTCF-s competes with the genome binding of canonical CTCF and binds a similar DNA sequence. CTCF-s binding disrupts CTCF/cohesin binding, alters CTCF-mediated chromatin looping and promotes the activation of IFI6 that leads to apoptosis. This effect is caused by an abnormal long-range interaction at the IFI6 enhancer and promoter. Taken together, this study reveals a non-canonical function for CTCF-s that antagonizes the genomic binding of canonical CTCF and cohesin, and that modulates chromatin looping and causes apoptosis by stimulating IFI6 expression.

CTCF Governs the Identity and Migration of MGE-Derived Cortical Interneurons.

The CCCTC-binding factor (CTCF) is a central regulator of chromatin topology recently linked to neurodevelopmental disorders such as intellectual disability, autism, and schizophrenia. The aim of this study was to identify novel roles of CTCF in the developing mouse brain. We provide evidence that CTCF is required for the expression of the LIM homeodomain factor LHX6 involved in fate determination of cortical interneurons (CINs) that originate in the medial ganglionic eminence (MGE). Conditional Ctcf ablation in the MGE of mice of either sex leads to delayed tangential migration, abnormal distribution of CIN in the neocortex, a marked reduction of CINs expressing parvalbumin and somatostatin (Sst), and an increased number of MGE-derived cells expressing Lhx8 and other markers of basal forebrain projection neurons. Likewise, Ctcf-null MGE cells transplanted into the cortex of wild-type hosts generate fewer Sst-expressing CINs and exhibit lamination defects that are efficiently rescued upon reexpression of LHX6. Collectively, these data indicate that CTCF regulates the dichotomy between Lhx6 and Lhx8 to achieve correct specification and migration of MGE-derived CINs.SIGNIFICANCE STATEMENT This work provides evidence that CCCTC-binding factor (CTCF) controls an early fate decision point in the generation of cortical interneurons mediated at least in part by Lhx6. Importantly, the abnormalities described could reflect early molecular and cellular events that contribute to human neurological disorders previously linked to CTCF, including schizophrenia, autism, and intellectual disability.

Disease-Associated Short Tandem Repeats Co-localize with Chromatin Domain Boundaries.

More than 25 inherited human disorders are caused by the unstable expansion of repetitive DNA sequences termed short tandem repeats (STRs). A fundamental unresolved question is why some STRs are susceptible to pathologic expansion, whereas thousands of repeat tracts across the human genome are relatively stable. Here, we discover that nearly all disease-associated STRs (daSTRs) are located at boundaries demarcating 3D chromatin domains. We identify a subset of boundaries with markedly higher CpG island density compared to the rest of the genome. daSTRs specifically localize to ultra-high-density CpG island boundaries, suggesting they might be hotspots for epigenetic misregulation or topological disruption linked to STR expansion. Fragile X syndrome patients exhibit severe boundary disruption in a manner that correlates with local loss of CTCF occupancy and the degree of FMR1 silencing. Our data uncover higher-order chromatin architecture as a new dimension in understanding repeat expansion disorders.

Dopamine Triggers CTCF-Dependent Morphological and Genomic Remodeling of Astrocytes.

Dopamine is critical for processing of reward and etiology of drug addiction. Astrocytes throughout the brain express dopamine receptors, but consequences of astrocytic dopamine receptor signaling are not well established. We found that extracellular dopamine triggered rapid concentration-dependent stellation of astrocytic processes that was not a result of dopamine oxidation but instead relied on both cAMP-dependent and cAMP-independent dopamine receptor signaling. This was accompanied by reduced duration and increased frequency of astrocytic Ca2+ transients, but little effect on astrocytic voltage-gated potassium channel currents. To isolate possible mechanisms underlying these structural and functional changes, we used whole-genome RNA sequencing and found prominent dopamine-induced enrichment of genes containing the CCCTC-binding factor (CTCF) motif, suggesting involvement of chromatin restructuring in the nucleus. CTCF binding to promoter sites bidirectionally regulates gene transcription and depends on activation of poly-ADP-ribose polymerase 1 (PARP1). Accordingly, antagonism of PARP1 occluded dopamine-induced changes, whereas a PARP1 agonist facilitated dopamine-induced changes on its own. These results indicate that astrocyte response to elevated dopamine involves PARP1-mediated CTCF genomic restructuring and concerted expression of gene networks. Our findings propose epigenetic regulation of chromatin landscape as a critical factor in the rapid astrocyte response to dopamine.SIGNIFICANCE STATEMENT Although dopamine is widely recognized for its role in modulating neuronal responses both in healthy and disease states, little is known about dopamine effects at non-neuronal cells in the brain. To address this gap, we performed whole-genome sequencing of astrocytes exposed to elevated extracellular dopamine and combined it with evaluation of effects on astrocyte morphology and function. We demonstrate a temporally dynamic pattern of genomic plasticity that triggers pronounced changes in astrocyte morphology and function. We further show that this plasticity depends on activation of genes sensitive to DNA-binding protein CTCF. Our results propose that a broad pattern of astrocyte responses to dopamine specifically relies on CTCF-dependent gene networks.