GATA-dependent transcriptional and epigenetic control of cardiac lineage specification and differentiation

Extracellular vesicles (EVs): A promising therapeutic tool in the heart tissue regeneration

Mesenchymal stem cells (MSCs) treatment has been widely explored as a therapy for myocardial infarction, peripheral ischemic vascular diseases, dilated cardiomyopathy, and pulmonary hypertension. Latest in vitro studies suggest that MSCs can differentiate into contractile cardiomyocytes. One of the best-characterized MSCs products are MSCs-derived extracellular vesicles (EVs). EVs are crucial paracrine effectors of MSCs. Based on previous works, paracrine effects of MSCs play a primary role in the regenerative ability. Hence, in the current paper, we focused our attention on an alternative approach, exploiting products derived from human dental pulp stem cells (hDPSCs) rather than MSCs themselves, which may denote a cost-effective and safer approach. The focus has been on EVs and the bioactive molecules they contain to evaluate their ability to influence the differentiation process toward cardiomyogenic lineage. The expression of GATA4, ACTC1, CX43, and Nkx2.5 was evaluated using Immunofluorescence, real time-PCR, and Western blotting analyses. Furthermore, the expression profiling analysis of the microRNA hsa-miR-200c-3p, targeting the GATA4 gene, was studied. The hsa-miR-200c-3p was found significantly down-regulated in both c-hDPSCs + EVs-hDPSCs and c-hDPSCs + EVs-HL-1 compared to untreated c-hDPSCs underlying a possible epigenetic mechanism behind the prevalent up-regulation of its targeted GATA4 gene. The aim of the present work was to develop an in vitro model of hDPSCs able to differentiate into cardiomyocytes in order to investigate the role of EVs derived from hDPSCs and derived from HL-1 cardiomyocyte cell line in modulating the differentiation process toward cardiomyogenic lineage.

Dynamic chromatin landscape encodes programs for perinatal transition of cardiomyocytes

The perinatal period occurring immediately before and after birth is critical for cardiomyocytes because they must change rapidly to accommodate the switch from fetal to neonatal circulation after birth. This transition is a well-orchestrated process, and any perturbation leads to unhealthy cardiomyocytes and heart disease. Despite its importance, little is known about how this transition is regulated and controlled. Here, by mapping the genome-wide chromatin accessibility, transcription-centered long-range chromatin interactions and gene expression in cardiomyocytes undergoing perinatal transition, we discovered two key transcription factors, MEF2 and AP1, that are crucial for driving the phenotypic changes within the perinatal window. Thousands of dynamic regulatory elements were found in perinatal cardiomyocytes and we show these elements mediated the transcriptional reprogramming through an elegant chromatin high-order architecture. We recompiled transcriptional program of induced stem cell-derived cardiomyocytes according to our discovered network, and they showed adult cardiomyocyte-like electrophysiological expression. Our work provides a comprehensive regulatory resource of cardiomyocytes perinatal reprogramming, and aids the gap-filling of cardiac translational research.

Towards Understanding the Gene-Specific Roles of GATA Factors in Heart Development: Does GATA4 Lead the Way?

Transcription factors play crucial roles in the regulation of heart induction, formation, growth and morphogenesis. Zinc finger GATA transcription factors are among the critical regulators of these processes. GATA4, 5 and 6 genes are expressed in a partially overlapping manner in developing hearts, and GATA4 and 6 continue their expression in adult cardiac myocytes. Using different experimental models, GATA4, 5 and 6 were shown to work together not only to ensure specification of cardiac cells but also during subsequent heart development. The complex involvement of these related gene family members in those processes is demonstrated through the redundancy among them and crossregulation of each other. Our recent identification at the genome-wide level of genes specifically regulated by each of the three family members and our earlier discovery that gata4 and gata6 function upstream, while gata5 functions downstream of noncanonical Wnt signalling during cardiac differentiation, clearly demonstrate the functional differences among the cardiogenic GATA factors. Such suspected functional differences are worth exploring more widely. It appears that in the past few years, significant advances have indeed been made in providing a deeper understanding of the mechanisms by which each of these molecules function during heart development. In this review, I will therefore discuss current evidence of the role of individual cardiogenic GATA factors in the process of heart development and emphasize the emerging central role of GATA4.

Zebrafish Models in Therapeutic Research of Cardiac Conduction Disease

Malfunction in the cardiac conduction system (CCS) due to congenital anomalies or diseases can cause cardiac conduction disease (CCD), which results in disturbances in cardiac rhythm, leading to syncope and even sudden cardiac death. Insights into development of the CCS components, including pacemaker cardiomyocytes (CMs), atrioventricular node (AVN) and the ventricular conduction system (VCS), can shed light on the pathological and molecular mechanisms underlying CCD, provide approaches for generating human pluripotent stem cell (hPSC)-derived CCS cells, and thus improve therapeutic treatment for such a potentially life-threatening disorder of the heart. However, the cellular and molecular mechanisms controlling CCS development remain elusive. The zebrafish has become a valuable vertebrate model to investigate early development of CCS components because of its unique features such as external fertilization, embryonic optical transparency and the ability to survive even with severe cardiovascular defects during development. In this review, we highlight how the zebrafish has been utilized to dissect the cellular and molecular mechanisms of CCS development, and how the evolutionarily conserved developmental mechanisms discovered in zebrafish could be applied to directing the creation of hPSC-derived CCS cells, therefore providing potential therapeutic strategies that may contribute to better treatment for CCD patients.

ANANSE: an enhancer network-based computational approach for predicting key transcription factors in cell fate determination

Proper cell fate determination is largely orchestrated by complex gene regulatory networks centered around transcription factors. However, experimental elucidation of key transcription factors that drive cellular identity is currently often intractable. Here, we present ANANSE (ANalysis Algorithm for Networks Specified by Enhancers), a network-based method that exploits enhancer-encoded regulatory information to identify the key transcription factors in cell fate determination. As cell type-specific transcription factors predominantly bind to enhancers, we use regulatory networks based on enhancer properties to prioritize transcription factors. First, we predict genome-wide binding profiles of transcription factors in various cell types using enhancer activity and transcription factor binding motifs. Subsequently, applying these inferred binding profiles, we construct cell type-specific gene regulatory networks, and then predict key transcription factors controlling cell fate transitions using differential networks between cell types. This method outperforms existing approaches in correctly predicting major transcription factors previously identified to be sufficient for trans-differentiation. Finally, we apply ANANSE to define an atlas of key transcription factors in 18 normal human tissues. In conclusion, we present a ready-to-implement computational tool for efficient prediction of transcription factors in cell fate determination and to study transcription factor-mediated regulatory mechanisms. ANANSE is freely available at https://github.com/vanheeringen-lab/ANANSE.

GATA4 Controls Epithelial Morphogenesis in the Developing Stomach to Promote Establishment of Glandular Columnar Epithelium

BACKGROUND & AIMS

The transcription factor GATA4 is broadly expressed in nascent foregut endoderm. As development progresses, GATA4 is lost in the domain giving rise to the stratified squamous epithelium of the esophagus and forestomach (FS), while it is maintained in the domain giving rise to the simple columnar epithelium of the hindstomach (HS). Differential GATA4 expression within these domains coincides with the onset of distinct tissue morphogenetic events, suggesting a role for GATA4 in diversifying foregut endoderm into discrete esophageal/FS and HS epithelial tissues. The goal of this study was to determine how GATA4 regulates differential morphogenesis of the mouse gastric epithelium.

METHODS

We used a Gata4 conditional knockout mouse line to eliminate GATA4 in the developing HS and a Gata4 conditional knock-in mouse line to express GATA4 in the developing FS.

RESULTS

We found that GATA4-deficient HS epithelium adopted a FS-like fate, and conversely, that GATA4-expressing FS epithelium adopted a HS-like fate. Underlying structural changes in these epithelia were broad changes in gene expression networks attributable to GATA4 directly activating or repressing expression of HS or FS defining transcripts. Our study implicates GATA4 as having a primary role in suppressing an esophageal/FS transcription factor network during HS development to promote columnar epithelium. Moreover, GATA4-dependent phenotypes in developmental mutants reflected changes in gene expression associated with Barrett's esophagus.

CONCLUSIONS

This study demonstrates that GATA4 is necessary and sufficient to activate the development of simple columnar epithelium, rather than stratified squamous epithelium, in the embryonic stomach. Moreover, similarities between mutants and Barrett's esophagus suggest that developmental biology can provide insight into human disease mechanisms.

BRD4 (Bromodomain-Containing Protein 4) Interacts with GATA4 (GATA Binding Protein 4) to Govern Mitochondrial Homeostasis in Adult Cardiomyocytes

BACKGROUND

Gene regulatory networks control tissue homeostasis and disease progression in a cell type-specific manner. Ubiquitously expressed chromatin regulators modulate these networks, yet the mechanisms governing how tissue specificity of their function is achieved are poorly understood. BRD4 (bromodomain-containing protein 4), a member of the BET (bromo- and extraterminal domain) family of ubiquitously expressed acetyl-lysine reader proteins, plays a pivotal role as a coactivator of enhancer signaling across diverse tissue types in both health and disease and has been implicated as a pharmacological target in heart failure. However, the cell-specific role of BRD4 in adult cardiomyocytes remains unknown.

METHODS

We combined conditional mouse genetics, unbiased transcriptomic and epigenomic analyses, and classic molecular biology and biochemical approaches to understand the mechanism by which BRD4 in adult cardiomyocyte homeostasis.

RESULTS

Here, we show that cardiomyocyte-specific deletion of Brd4 in adult mice leads to acute deterioration of cardiac contractile function with mutant animals demonstrating a transcriptomic signature characterized by decreased expression of genes critical for mitochondrial energy production. Genome-wide occupancy data show that BRD4 enriches at many downregulated genes (including the master coactivators Ppargc1a, Ppargc1b, and their downstream targets) and preferentially colocalizes with GATA4 (GATA binding protein 4), a lineage-determining cardiac transcription factor not previously implicated in regulation of adult cardiac metabolism. BRD4 and GATA4 form an endogenous complex in cardiomyocytes and interact in a bromodomain-independent manner, revealing a new functional interaction partner for BRD4 that can direct its locus and tissue specificity.

CONCLUSIONS

These results highlight a novel role for a BRD4-GATA4 module in cooperative regulation of a cardiomyocyte-specific gene program governing bioenergetic homeostasis in the adult heart.

Combination of 5-azaytidine and hanging drop culture convert fat cell into cardiac cell

One of the promising approaches for the treatment of cardiac disease is stem cell therapy. In this study, we compared the cardiomyogenic differentiation rate, from human adipose-derived stem cells (hADSCs) in a three-dimensional (3D) hanging drop (HD) spheroid culture system, versus a two-dimensional (2D) culture condition at different concentrations of 5-azacytidine (5-Aza). 5-Azaytidine (5-Aza) is a pyrimidine nucleoside analogue of cytidine that initiates cell differentiation programs through DNA demethylation. The hADSCs were isolated and cultured both in 2D and 3D HD conditions, with either 10 or 50 μM concentrations of 5-Aza. Then DNA content, gene expression, and protein content were analyzed. 3D HD culture resulted in a higher percentage of cells in G0/G1 and S phase in the cell division cycle, whereas 2D culture led to a greater percentage of cells in the G2/M phase. A significantly higher gene expression rate of HAND1, HAND2, cTnI, Cx43, βMHC, GATA4, NKX2.5, and MLC2V was observed in HD treated with 50 μM 5-Aza. This was confirmed by immunocytochemistry. These findings suggest that 50 μM concentration of 5-Aza can induce hADSCs to differentiate into cardiomyocytes. The differentiation rate was significantly higher when accompanied by the 3D HD culture system. This work provides a new culture system for cell differentiation for cardiovascular tissue engineering.

Variation in a Left Ventricle-Specific Hand1 Enhancer Impairs GATA Transcription Factor Binding and Disrupts Conduction System Development and Function

RATIONALE

The ventricular conduction system (VCS) rapidly propagates electrical impulses through the working myocardium of the ventricles to coordinate chamber contraction. GWAS (Genome-wide association studies) have associated nucleotide polymorphisms, most are located within regulatory intergenic or intronic sequences, with variation in VCS function. Two highly correlated polymorphisms (r²>0.99) associated with VCS functional variation (rs13165478 and rs13185595) occur 5' to the gene encoding the basic helix-loop-helix transcription factor HAND1 (heart- and neural crest derivatives-expressed protein 1).

OBJECTIVE

Here, we test the hypothesis that these polymorphisms influence HAND1 transcription thereby influencing VCS development and function.

METHODS AND RESULTS

We employed transgenic mouse models to identify an enhancer that is sufficient for left ventricle (LV) cis-regulatory activity. Two evolutionarily conserved GATA transcription factor cis-binding elements within this enhancer are bound by GATA4 and are necessary for cis-regulatory activity, as shown by in vitro DNA binding assays. CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9-mediated deletion of this enhancer dramatically reduces Hand1 expression solely within the LV but does not phenocopy previously published mouse models of cardiac Hand1 loss-of-function. Electrophysiological and morphological analyses reveals that mice homozygous for this deleted enhancer display a morphologically abnormal VCS and a conduction system phenotype consistent with right bundle branch block. Using 1000 Genomes Project data, we identify 3 additional single nucleotide polymorphisms (SNPs), located within the Hand1 LV enhancer, that compose a haplotype with rs13165478 and rs13185595. One of these SNPs, rs10054375, overlaps with a critical GATA cis-regulatory element within the Hand1 LV enhancer. This SNP, when tested in electrophoretic mobility shift assays, disrupts GATA4 DNA-binding. Modeling 2 of these SNPs in mice causes diminished Hand1 expression and mice present with abnormal VCS function.

CONCLUSIONS

Together, these findings reveal that SNP rs10054375, which is located within a necessary and sufficient LV-specific Hand1 enhancer, exhibits reduces GATA DNA-binding in electrophoretic mobility shift assay, and this enhancer in total, is required for VCS development and function in mice and perhaps humans.

Drosophila Mediator Subunit Med1 Is Required for GATA-Dependent Developmental Processes: Divergent Binding Interfaces for Conserved Coactivator Functions

DNA-bound transcription factors (TFs) governing developmental gene regulation have been proposed to recruit polymerase II machinery at gene promoters through specific interactions with dedicated subunits of the evolutionarily conserved Mediator (MED) complex. However, whether such MED subunit-specific functions and partnerships have been conserved during evolution has been poorly investigated. To address this issue, we generated the first Drosophila melanogaster loss-of-function mutants for Med1, known as a specific cofactor for GATA TFs and hormone nuclear receptors in mammals. We show that Med1 is required for cell proliferation and hematopoietic differentiation depending on the GATA TF Serpent (Srp). Med1 physically binds Srp in cultured cells and in vitro through its conserved GATA zinc finger DNA-binding domain and the divergent Med1 C terminus. Interestingly, GATA-Srp interaction occurs through the longest Med1 isoform, suggesting a functional diversity of MED complex populations. Furthermore, we show that Med1 acts as a coactivator for the GATA factor Pannier during thoracic development. In conclusion, the Med1 requirement for GATA-dependent regulatory processes is a common feature in insects and mammals, although binding interfaces have diverged. Further work in Drosophila should bring valuable insights to fully understand GATA-MED functional partnerships, which probably involve other MED subunits depending on the cellular context.

pouC Regulates Expression of bmp4 During Atrioventricular Canal Formation in Zebrafish

BACKGROUND

Many human gene mutations have been linked to congenital heart disease (CHD), yet CHD remains a major health issue worldwide due in part to an incomplete understanding of the molecular basis for cardiac malformation.

RESULTS

Here we identify the orthologous mouse Pou6f1 and zebrafish pouC as POU homeodomain transcription factors enriched in the developing heart. We find that pouC is a multi-functional transcriptional regulator containing separable activation, repression, protein-protein interaction, and DNA binding domains. Using zebrafish heart development as a model system, we demonstrate that pouC knockdown impairs cardiac morphogenesis and affects cardiovascular function. We also find that levels of pouC expression must be fine-tuned to enable proper heart formation. At the cellular level, we demonstrate that pouC knockdown disrupts atrioventricular canal (AVC) cardiomyocyte maintenance, although chamber myocyte specification remains intact. Mechanistically, we show that pouC binds a bmp4 intronic regulatory element to mediate transcriptional activation.

CONCLUSIONS

Taken together, our study establishes pouC as a novel transcriptional input into the regulatory hierarchy that drives AVC morphogenesis in zebrafish. We anticipate that these findings will inform future efforts to explore functional conservation in mammals and potential association with atrioventricular septal defects in humans. Developmental Dynamics 248:173-188, 2019. © 2018 Wiley Periodicals, Inc.

Variants in cardiac GATA genes associated with bicuspid aortic valve

BACKGROUND

Bicuspid aortic valve (BAV) is a heterogeneous and still not fully understood condition, with diverse genetic aetiology and associated phenotypes ranging from aortic stenosis or regurgitation to aneurysm and dissection. Several genes have been associated with the presence of BAV, notably some members of the GATA family of transcription factors that play important roles in cardiac embryogenesis.

METHODS

A case-control study with 122 unrelated and ethnically matched patients with bicuspid and 154 with tricuspid aortic valve was performed. All exons of GATA4, GATA5, and GATA6 genes were sequenced searching for variants. Frequencies were compared and functional effects of rare variants were analysed by informatic prediction tools.

RESULTS

Four rare and potentially pathogenic variants were identified in only five sporadic cases: a missense p.Arg202Gln (rs782614097) in GATA5 and three synonymous variants, p.Cys274= (rs55980825) and p.His302= (rs201516339) in GATA4, and p.Asn458= (rs143026087) in GATA6. Minor alleles of p.His302=, p.Arg202Gln and rs3764962 are enriched in BAV patients compared to ExAC database (P = 0.01, 0.03, and 0.01 respectively). In addition, a common polymorphism in GATA5 (p.Asp203=, rs41305803) is associated with BAV showing a protective effect in carriers of the minor allele (OR [95%CI] = 0.45[0.25-0.81]; P = 0.004).

CONCLUSION

This study associates additional genetic variants in GATA4 and GATA5 with BAV, supporting the implication of these genes in the development of this valvulopathy. The discovery of all the genetic factors involved will contribute to a better understanding of the process and, therefore, to detect a genetic predisposition and even to the identification of therapeutic targets.

Genome and epigenome analysis of monozygotic twins discordant for congenital heart disease

BACKGROUND

Congenital heart disease (CHD) is the leading non-infectious cause of death in infants. Monozygotic (MZ) twins share nearly all of their genetic variants before and after birth. Nevertheless, MZ twins are sometimes discordant for common complex diseases. The goal of this study is to identify genomic and epigenomic differences between a pair of twins discordant for a form of congenital heart disease, double outlet right ventricle (DORV).

RESULTS

A monoamniotic monozygotic (MZ) twin pair discordant for DORV were subjected to genome-wide sequencing and methylation analysis. We identified few genomic differences but 1566 differentially methylated regions (DMRs) between the MZ twins. Twenty percent (312/1566) of the DMRs are located within 2 kb upstream of transcription start sites (TSS), containing 121 binding sites of transcription factors. Particularly, ZIC3 and NR2F2 are found to have hypermethylated promoters in both the diseased twin and additional patients suffering from DORV.

CONCLUSIONS

The results showed a high correlation between hypermethylated promoters at ZIC3 and NR2F2 and down-regulated gene expression levels of these two genes in patients with DORV compared to normal controls, providing new insight into the potential mechanism of this rare form of CHD.

HAND2 Target Gene Regulatory Networks Control Atrioventricular Canal and Cardiac Valve Development

The HAND2 transcriptional regulator controls cardiac development, and we uncover additional essential functions in the endothelial to mesenchymal transition (EMT) underlying cardiac cushion development in the atrioventricular canal (AVC). In Hand2-deficient mouse embryos, the EMT underlying AVC cardiac cushion formation is disrupted, and we combined ChIP-seq of embryonic hearts with transcriptome analysis of wild-type and mutants AVCs to identify the functionally relevant HAND2 target genes. The HAND2 target gene regulatory network (GRN) includes most genes with known functions in EMT processes and AVC cardiac cushion formation. One of these is Snai1, an EMT master regulator whose expression is lost from Hand2-deficient AVCs. Re-expression of Snai1 in mutant AVC explants partially restores this EMT and mesenchymal cell migration. Furthermore, the HAND2-interacting enhancers in the Snai1 genomic landscape are active in embryonic hearts and other Snai1-expressing tissues. These results show that HAND2 directly regulates the molecular cascades initiating AVC cardiac valve development.

Novel mutation of GATA4 gene in Kurdish population of Iran with nonsyndromic congenital heart septals defects

BACKGROUND

The mutations in GATA4 gene induce inherited atrial and ventricular septation defects, which is the most frequent forms of congenital heart defects (CHDs) constituting about half of all cases.

METHOD

We have performed High resolution melting (HRM) mutation scanning of GATA4 coding exons of nonsyndrome 100 patients as a case group including 39 atrial septal defects (ASD), 57 ventricular septal defects (VSD) and four patients with both above defects and 50 healthy individuals as a control group. Our samples are categorized according to their HRM graph. The genome sequencing has been done for 15 control samples and 25 samples of patients whose HRM analysis were similar to healthy subjects for each exon. The PolyPhen-2 and MUpro have been used to determine the causative possibility and structural stability prediction of GATA4 sequence variation.

RESULTS

The HRM curve analysis exhibit that 21 patients and 3 normal samples have deviated curves for GATA4 coding exons. Sequencing analysis has revealed 12 nonsynonymous mutations while all of them resulted in stability structure of protein 10 of them are pathogenic and 2 of them are benign. Also we found two nucleotide deletions which one of them was novel and one new indel mutation resulting in frame shift mutation, and 4 synonymous variations or polymorphism in 6 of patients and 3 of normal individuals. Six or about 50% of these nonsynonymous mutations have not been previously reported.

CONCLUSION

Our results show that there is a spectrum of GATA4 mutations resulting in septal defects.

Cardiac reprogramming factor Gata4 reduces postinfarct cardiac fibrosis through direct repression of the profibrotic mediator snail

OBJECTIVE

The administration of a variety of reprogramming factor cocktails has now been shown to reprogram cardiac fibroblasts into induced cardiomyocyte-like cells. However, reductions in ventricular fibrosis observed after reprogramming factor administration seem to far exceed the extent of induced cardiomyocyte-like cell generation in vivo. We investigated whether reprogramming factor administration might primarily play a role in activating antifibrotic molecular pathways.

METHODS

Adult rat cardiac fibroblasts were infected with lentivirus encoding the transcription factors Gata4, Mef2c, or Tbx5, all 3 vectors, or a green fluorescent protein control vector. Gene and protein expression assays were performed to identify relevant antifibrotic targets of these factors. The antifibrotic effects of these factors were then investigated in a rat coronary ligation model.

RESULTS

Gata4, Mef2c, or Tbx5 administration to rat cardiac fibroblasts in vitro significantly downregulated expression of Snail and the profibrotic factors connective tissue growth factor, collagen1a1, and fibronectin. Of these factors, Gata4 was shown to be the one responsible for the downregulation of the profibrotic factors and Snail (mRNA expression fold change relative to green fluorescent protein for Snail, Gata4: 0.5 ± 0.3, Mef2c: 1.3 ± 1.0, Tbx5: 0.9 ± 0.5, Gata4, Mef2c, or Tbx5: 0.6 ± 0.2, P < .05). Chromatin immunoprecipitation quantitative polymerase chain reaction identified Gata4 binding sites in the Snail promoter. In a rat coronary ligation model, only Gata4 administration alone improved postinfarct ventricular function and reduced the extent of postinfarct fibrosis.

CONCLUSIONS

Gata4 administration reduces postinfarct ventricular fibrosis and improves ventricular function in a rat coronary ligation model, potentially as a result of Gata4-mediated downregulation of the profibrotic mediator Snail.

Mechanisms of retinoic acid signaling during cardiogenesis

Substantial experimental and epidemiological data have highlighted the interplay between nutritional and genetic factors in the development of congenital heart defects. Retinoic acid (RA), a derivative of vitamin A, plays a key role during vertebrate development including the formation of the heart. Retinoids bind to RA and retinoid X receptors (RARs and RXRs) which then regulate tissue-specific genes. Here, we will focus on the roles of RA signaling and receptors in gene regulation during cardiogenesis, and the consequence of deregulated retinoid signaling on heart formation and congenital heart defects.

The formation and function of the cardiac conduction system

The cardiac conduction system (CCS) consists of distinctive components that initiate and conduct the electrical impulse required for the coordinated contraction of the cardiac chambers. CCS development involves complex regulatory networks that act in stage-, tissue- and dose-dependent manners, and recent findings indicate that the activity of these networks is sensitive to common genetic variants associated with cardiac arrhythmias. Here, we review how these findings have provided novel insights into the regulatory mechanisms and transcriptional networks underlying CCS formation and function.