Epigenetic regulation of Wnt7b expression by the cis-acting long noncoding RNA Lnc-Rewind in muscle stem cells

Long non-coding RNAs and their role in muscle regeneration

In mammals, most of the genome is transcribed to generate a large and heterogeneous variety of non-protein coding RNAs, that are broadly grouped according to their size. Long noncoding RNAs include a very large and versatile group of molecules. Despite only a minority of them has been functionally characterized, there is emerging evidence indicating long noncoding RNAs as important regulators of expression at multiple levels. Several of them have been shown to be modulated during myogenic differentiation, playing important roles in the regulation of skeletal muscle development, differentiation and homeostasis, and contributing to neuromuscular diseases. In this chapter, we have summarized the current knowledge about long noncoding RNAs in skeletal muscle and discussed specific examples of long noncoding RNAs (lncRNAs and circRNAs) regulating muscle stem cell biology. We have also discussed selected long noncoding RNAs involved in the most common neuromuscular diseases.

The pancancer overexpressed NFYC Antisense 1 controls cell cycle mitotic progression through in cis and in trans modes of action

Antisense RNAs (asRNAs) represent an underappreciated yet crucial layer of gene expression regulation. Generally thought to modulate their sense genes in cis through sequence complementarity or their act of transcription, asRNAs can also regulate different molecular targets in trans, in the nucleus or in the cytoplasm. Here, we performed an in-depth molecular characterization of NFYC Antisense 1 (NFYC-AS1), the asRNA transcribed head-to-head to NFYC subunit of the proliferation-associated NF-Y transcription factor. Our results show that NFYC-AS1 is a prevalently nuclear asRNA peaking early in the cell cycle. Comparative genomics suggests a narrow phylogenetic distribution, with a probable origin in the common ancestor of mammalian lineages. NFYC-AS1 is overexpressed pancancer, preferentially in association with RB1 mutations. Knockdown of NFYC-AS1 by antisense oligonucleotides impairs cell growth in lung squamous cell carcinoma and small cell lung cancer cells, a phenotype recapitulated by CRISPR/Cas9-deletion of its transcription start site. Surprisingly, expression of the sense gene is affected only when endogenous transcription of NFYC-AS1 is manipulated. This suggests that regulation of cell proliferation is at least in part independent of the in cis transcription-mediated effect on NFYC and is possibly exerted by RNA-dependent in trans effects converging on the regulation of G2/M cell cycle phase genes. Accordingly, NFYC-AS1-depleted cells are stuck in mitosis, indicating defects in mitotic progression. Overall, NFYC-AS1 emerged as a cell cycle-regulating asRNA with dual action, holding therapeutic potential in different cancer types, including the very aggressive RB1-mutated tumors.

Chromatin organization of muscle stem cell

The proper functioning of skeletal muscles is essential throughout life. A crucial crosstalk between the environment and several cellular mechanisms allows striated muscles to perform successfully. Notably, the skeletal muscle tissue reacts to an injury producing a completely functioning tissue. The muscle's robust regenerative capacity relies on the fine coordination between muscle stem cells (MuSCs or "satellite cells") and their specific microenvironment that dictates stem cells' activation, differentiation, and self-renewal. Critical for the muscle stem cell pool is a fine regulation of chromatin organization and gene expression. Acquiring a lineage-specific 3D genome architecture constitutes a crucial modulator of muscle stem cell function during development, in the adult stage, in physiological and pathological conditions. The context-dependent relationship between genome structure, such as accessibility and chromatin compartmentalization, and their functional effects will be analysed considering the improved 3D epigenome knowledge, underlining the intimate liaison between environmental encounters and epigenetics.

LINC02454 promotes thyroid carcinoma progression via upregulating HMGA2 through CREB1

Thyroid carcinoma (THCA) is the most common malignancy in the endocrine system. Long intergenic non-coding RNA 2454 (LINC02454) exhibits an HMGA2-like expression pattern, but their relationship and roles in THCA are largely unknown. The present purpose was to delineate the roles of LINC02454 in THCA progression and its molecular mechanisms. We collected THCA tissues from patients and monitored patient survival. THCA cell colony formation, migration, and invasion were evaluated. Metastasis was evaluated by examining EMT markers through Western blotting. Gene interaction was determined with ChIP, RIP, RNA pull-down, and luciferase activity assays. A mouse model of a subcutaneous tumor was used to determine the activity of LINC02454 knockdown in vivo. We found that LINC02454 was highly expressed in THCA, and its upregulation was associated with poor survival. The knockdown of LINC02454 repressed colony formation, migration, and invasion. Moreover, loss of LINC02454 inhibited tumor growth and metastasis in mice. HMGA2 promoted LINC02454 transcription via binding to the LINC02454 promoter, and silencing of HMGA2 suppressed malignant behaviors through downregulation of LINC02454. HMGA2 was a novel functional target of LINC02454 in THCA cells, and knockdown of LINC02454-mediated anti-tumor effects was reversed by HMGA2 overexpression. Mechanically, LINC02454 promoted CREB1 phosphorylation and nuclear translocation, and CREB1 was subsequently bound to the HMGA2 promoter to facilitate its expression. LINC02454 cis-regulates HMGA2 transcription via facilitating CREB1 phosphorylation and nuclear translocation, and, in turn, HMGA2 promotes LINC02454 expression, thus accelerating thyroid carcinoma progression. Our results support therapeutic targets of LINC02454 and HMGA2 for THCA.

Preparation of Cardiac Extracts from Embryonal Hearts to Capture RNA-protein Interactions by CLIP

The interaction of RNA with specific RNA-binding proteins (RBP) leads to the establishment of complex regulatory networks through which gene expression is controlled. Careful consideration should be given to the exact environment where a given RNA/RBP interplay occurs, as the functional responses might depend on the type of organism as well as the specific cellular or subcellular contexts. This requisite becomes particularly crucial for the study of long non-coding RNAs (lncRNA), as a consequence of their peculiar tissue-specificity and timely regulated expression. The functional characterization of lncRNAs has traditionally relied on the use of established cell lines that, although useful, are unable to fully recapitulate the complexity of a tissue or organ. Here, we detail an optimized protocol, with comments and tips, to identify the RNA interactome of given RBPs by performing cross-linking immunoprecipitation (CLIP) from mouse embryonal hearts. We tested the efficiency of this protocol on the murine pCharme, a muscle-specific lncRNA interacting with Matrin3 (MATR3) and forming RNA-enriched condensates of biological significance in the nucleus. Key features • The protocol refines previous methods of cardiac extracts preparation to use for CLIP assays. • The protocol allows the quantitative RNA-seq analysis of transcripts interacting with selected proteins. • Depending on the embryonal stage, a high number of hearts can be required as starting material. • The steps are adaptable to other tissues and biochemical assays.

Down-regulation of OIP5-AS1 inhibits obesity-induced myocardial pyroptosis and miR-22/NLRP3 inflammasome axis

BACKGROUND

Obesity can induce myocardial pyroptosis, but the exact mechanism is still unknown. A recent study reported the association of opa-interacting protein 5-antisense transcript 1 (OIP5-AS1), an evolutionarily conserved long noncoding RNA, with pyroptosis. Therefore, this study aimed to investigate the role of OIP5-AS1 in obesity-induced myocardial pyroptosis.

METHODS

OIP5-AS1 was downregulated in H9c2 cells, followed by treatment with 400 μM palmitic acid (PA). Propidium iodide (PI) staining, lactic dehydrogenase (LDH) release assay, caspase-1 activity assay, IL-1β, and IL-18 activity assay were performed to detect pyroptotic phenotype. The interaction between OIP5-AS1 and microRNAs (miRNAs) was analyzed using RNA pull-down and luciferase assay. The effect of OIP5-AS1 knockdown in high-fat diet (HFD)-induced obesity rat on cardiac function, myocardial hypertrophy, fibrosis, and remodeling was evaluated.

RESULTS

Fat deposition was observed in cardiomyocytes 24 h after PA treatment; moreover, PA-treated cardiomyocytes showed significant increase in the rate of pyroptotic cells, release of LDH, protein expressions of NLRP3 and cleaved caspase-1, and the activity of caspase-1, IL-1β, and IL-18 as well as OIP5-AS1 expression. These findings suggested that PA activated pyroptosis and induced OIP5-AS1 expression in cardiomyocytes. Moreover, OIP5-AS1 knockdown inhibited PA-induced pyroptosis. Mechanistically, OIP5-AS1 was found to specifically bind to miR-22 and to regulate NLRP3 inflammasome-mediated pyroptosis via miR-22. Furthermore, OIP5-AS1 knockdown ameliorated HFD-induced cardiac dysfunction, myocardial hypertrophy, fibrosis, remodeling, and pyroptosis.

CONCLUSION

Our results revealed that downregulation of OIP5-AS1 can inhibit obesity-induced myocardial pyroptosis via miR-22/NLRP3 inflammasome axis. This finding lays a foundation of gene therapy for heart disease targeting OIP5-AS1.

Exploring the landscape of tools and resources for the analysis of long non-coding RNAs

In recent years, research on long non-coding RNAs (lncRNAs) has gained considerable attention due to the increasing number of newly identified transcripts. Several characteristics make their functional evaluation challenging, which called for the urgent need to combine molecular biology with other disciplines, including bioinformatics. Indeed, the recent development of computational pipelines and resources has greatly facilitated both the discovery and the mechanisms of action of lncRNAs. In this review, we present a curated collection of the most recent computational resources, which have been categorized into distinct groups: databases and annotation, identification and classification, interaction prediction, and structure prediction. As the repertoire of lncRNAs and their analysis tools continues to expand over the years, standardizing the computational pipelines and improving the existing annotation of lncRNAs will be crucial to facilitate functional genomics studies.

Transcriptome Analysis Reveals the Profile of Long Non-Coding RNAs during Myogenic Differentiation in Goats

The long non-coding RNAs (lncRNAs) are emerging as essential regulators of the growth and development of skeletal muscles. However, little is known about the expression profiles of lncRNAs during the proliferation and differentiation of skeletal muscle satellite cells (MuSCs) in goats. In this study, we investigate potential regulatory lncRNAs that govern muscle development by performing lncRNA expression profiling analysis during the proliferation (cultured in the growth medium, GM) and differentiation (cultured in the differentiation medium, DM1/DM5) of MuSCs. In total, 1001 lncRNAs were identified in MuSC samples, and 314 differentially expressed (DE) (FDR < 0.05, |log2FC| > 1) lncRNAs were screened by pairwise comparisons from three comparison groups (GM-vs-DM1, GM-vs-DM5, DM1-vs-DM5). Moreover, we identified the cis-, trans-, and antisense-regulatory target genes of DE lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that these target genes were significantly enriched in muscle development-related GO terms and KEGG pathways. In addition, the network of interactions between DE lncRNAs and their target genes was identified, which included well-known myogenesis regulators such as Myogenic differentiation 1 (MyoD), Myogenin (MyoG), and Myosin heavy chain (MyHC). Meanwhile, competing endogenous RNA (ceRNA) network analysis showed that 237 DE lncRNAs could bind to 329 microRNAs (miRNAs), while miRNAs could target 564 mRNAs. Together, our results provide a genome-wide resource of lncRNAs that may contribute to myogenic differentiation in goats and lay the groundwork for future investigation into their functions during skeletal muscle development.

Clinical Prognostic Implications of Wnt Hub Genes Expression in Medulloblastoma

Medulloblastoma is the most common type of pediatric malignant primary brain tumor, and about one-third of patients die due to disease recurrence and most survivors suffer from long-term side effects. MB is clinically, genetically, and epigenetically heterogeneous and subdivided into at least four molecular subgroups: WNT, SHH, Group 3, and Group 4. We evaluated common differentially expressed genes between a Brazilian RNA-seq GSE181293 dataset and microarray GSE85217 dataset cohort of pediatric MB samples using bioinformatics methodology in order to identify hub genes of the molecular subgroups based on PPI network construction, survival and functional analysis. The main finding was the identification of five hub genes from the WNT subgroup that are tumor suppressors, and whose lower expression is related to a worse prognosis for MB patients. Furthermore, the common genes correlated with the five tumor suppressors participate in important pathways and processes for tumor initiation and progression, as well as development and differentiation, and some of them control cell stemness and pluripotency. These genes have not yet been studied within the context of MB, representing new important elements for investigation in the search for therapeutic targets, prognostic markers or for understanding of MB biology.

The chromatin-associated RNAs in gene regulation and cancer

Eukaryotic genomes are prevalently transcribed into many types of RNAs that translate into proteins or execute gene regulatory functions. Many RNAs associate with chromatin directly or indirectly and are called chromatin-associated RNAs (caRNAs). To date, caRNAs have been found to be involved in gene and transcriptional regulation through multiple mechanisms and have important roles in different types of cancers. In this review, we first present different categories of caRNAs and the modes of interaction between caRNAs and chromatin. We then detail the mechanisms of chromatin-associated nascent RNAs, chromatin-associated noncoding RNAs and emerging m⁶A on caRNAs in transcription and gene regulation. Finally, we discuss the roles of caRNAs in cancer as well as epigenetic and epitranscriptomic mechanisms contributing to cancer, which could provide insights into the relationship between different caRNAs and cancer, as well as tumor treatment and intervention.

Regulation of muscle stem cell fate

Skeletal muscle plays a critical role in human health. Muscle stem cells (MuSCs) serve as the major cell type contributing to muscle regeneration by directly differentiating to mature muscle cells. MuSCs usually remain quiescent with occasionally self-renewal and are activated to enter cell cycle for proliferation followed by differentiation upon muscle injury or under pathological conditions. The quiescence maintenance, activation, proliferation, and differentiation of MuSCs are tightly regulated. The MuSC cell-intrinsic regulatory network and the microenvironments work coordinately to orchestrate the fate transition of MuSCs. The heterogeneity of MuSCs further complicates the regulation of MuSCs. This review briefly summarizes the current progress on the heterogeneity of MuSCs and the microenvironments, epigenetic, and transcription regulations of MuSCs.

The RNA helicase DDX5 cooperates with EHMT2 to sustain alveolar rhabdomyosarcoma growth

Rhabdomyosarcoma (RMS) is the most common soft-tissue sarcoma of childhood characterized by the inability to exit the proliferative myoblast-like stage. The alveolar fusion positive subtype (FP-RMS) is the most aggressive and is mainly caused by the expression of PAX3/7-FOXO1 oncoproteins, which are challenging pharmacological targets. Here, we show that the DEAD box RNA helicase 5 (DDX5) is overexpressed in alveolar RMS cells and that its depletion and pharmacological inhibition decrease FP-RMS viability and slow tumor growth in xenograft models. Mechanistically, we provide evidence that DDX5 functions upstream of the EHMT2/AKT survival signaling pathway, by directly interacting with EHMT2 mRNA, modulating its stability and consequent protein expression. We show that EHMT2 in turns regulates PAX3-FOXO1 activity in a methylation-dependent manner, thus sustaining FP-RMS myoblastic state. Together, our findings identify another survival-promoting loop in FP-RMS and highlight DDX5 as a potential therapeutic target to arrest RMS growth.

Long non-coding RNAs in <i>Sus scrofa</i> ileum under starvation stress

Objective

In this study, we aimed to identify long non-coding RNAs (lncRNAs) that play important roles in starvation stress, analyze their functions, and discover potential molecular targets to alleviate starvation stress to provide a theoretical reference for subsequent in-depth research.

Methods

We generated a piglet starvation stress animal model. Nine Yorkshire weaned piglets were randomly divided into a long-term starvation stress group (starved for 72 h), short-term starvation stress group (starved for 48 h), and the control group. LncRNA libraries were constructed using high-throughput sequencing of piglet ileums.

Results

We obtained 11,792 lncRNAs, among which, 2,500 lncRNAs were novel. In total, 509 differentially expressed (DE)lncRNAs were identified in this study. Target genes of DElncRNAs were predicted via cis and trans interactions, and functional and pathway analyses were performed. Gene ontology functions and Kyoto encyclopedia of genes and genomes analysis revealed that lncRNA-targeted genes mainly participated in metabolic pathways, cellular processes, immune system processes, digestive systems, and transport activities. To reveal the mechanism underlying starvation stress, the interaction network between lncRNAs and their targets was constructed based on 26 DElncRNAs and 72 DEmRNAs. We performed an interaction network analysis of 121 DElncRNA–DEmRNA pairs with a Pearson correlation coefficient greater than 0.99.

Conclusion

We found that MSTRG.19894.13, MSTRG.16726.3, and MSTRG.12176.1 might play important roles in starvation stress. This study not only generated a library of enriched lncRNAs in piglets, but its outcomes also provide a strong foundation to screen key lncRNAs involved in starvation stress and a reference for subsequent in-depth research.

LncRNA UCA1 mediates Cetuximab resistance in Colorectal Cancer via the MiR-495 and HGF/c-MET Pathways

Background: Cetuximab is one of the most widely used monoclonal antibodies to treat patients with RAS/BRAF wild-type metastatic colorectal cancer (mCRC). Unfortunately, cetuximab resistance often occurs during targeted therapy. However, the underlying epigenetic mechanisms remain unclear. Our previous study demonstrated that the exosomal transfer of urothelial carcinoma-associated 1 (UCA1) confers cetuximab resistance to CRC cells. The goal of this study was to elucidate the detailed role of UCA1 in cetuximab resistance in CRC and the underlying molecular mechanism.

Methods:In vitro and in vivo functional studies were performed to assess the role of UCA1 in cetuximab resistance in CRC cell lines and xenograft models. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was used to examine UCA1 localization and expression. Bioinformatics analysis was performed to predict the potential mechanism of UCA1, which was further validated by the dual-luciferase reporter assay and the RNA immunoprecipitation (RIP) assay. Cells treated with indicators were subjected to Cell Counting Kit-8 (CCK-8) and western blotting to investigate the role of hepatocyte growth factor (HGF)/c-mesenchymal-epithelial transition (c-MET) signalling in UCA1-mediated cetuximab resistance.

Results: We showed that UCA1 decreased CRC cell sensitivity to cetuximab by suppressing apoptosis. Mechanistic studies revealed that UCA1 promoted cetuximab resistance by competitively binding miR-495 to facilitate HGF and c-MET expression in CRC cells. Moreover, HGF was shown to attenuate the cetuximab-induced inhibition of cell proliferation by activating the HGF/c-MET pathway in CRC cells.

Conclusion: We provide the first evidence of a UCA1-miR-495-HGF/c-MET regulatory network involved in cetuximab resistance in CRC. Therefore, UCA1 has potential as a predictor and therapeutic target for cetuximab resistance.

Targeting the Expression of Long Noncoding RNAs in Murine Satellite Cells from Single Myofibers

LncRNAs have been recently implicated in the epigenetic control of muscle differentiation and their functional characterization has traditionally relied upon in vitro models of myogenic differentiation. However, the use of experimental paradigms to specifically target lncRNAs expression in muscle stem cells (MuSCs), also known as satellite cells, represents an important requisite to interrogate their function in more physiological contexts. Since isolation and culture of single myofibers preserves satellite cells within their physiological niche underneath the surrounding basal lamina, this procedure represents the optimal approach to follow satellite cell dynamics ex-vivo, such as activation from quiescence, expansion of committed progenitors, differentiation, and self-renewal. Here, we detail an optimized protocol to isolate viable single myofibers from the extensor digitorum longus (EDL) skeletal muscle of adult mice and to manipulate the expression of lncRNAs by antisense LNA GapmeRs-mediated knock-down (KD). Furthermore, we describe a method of EdU incorporation that, coupled to lncRNA KD and subsequent immunofluorescence analysis of proliferating, differentiating, and satellite cell-specific markers, permits the inference of lncRNAs function on muscle stem cells dynamics. Graphic abstract: Graphical representation of the single myofiber isolation method. Experimental workflow showing the main steps of the protocol procedure: EDL muscle harvesting from the mouse hindlimb; EDL digestion into single myofibers; transfection with antisense oligos and culture for 96h; immunofluorescence protocol and image outcome.

Muscle Regeneration and RNA: New Perspectives for Ancient Molecules

The ability of the ribonucleic acid (RNA) to self-replicate, combined with a unique cocktail of chemical properties, suggested the existence of an RNA world at the origin of life. Nowadays, this hypothesis is supported by innovative high-throughput and biochemical approaches, which definitively revealed the essential contribution of RNA-mediated mechanisms to the regulation of fundamental processes of life. With the recent development of SARS-CoV-2 mRNA-based vaccines, the potential of RNA as a therapeutic tool has received public attention. Due to its intrinsic single-stranded nature and the ease with which it is synthesized in vitro, RNA indeed represents the most suitable tool for the development of drugs encompassing every type of human pathology. The maximum effectiveness and biochemical versatility is achieved in the guise of non-coding RNAs (ncRNAs), which are emerging as multifaceted regulators of tissue specification and homeostasis. Here, we report examples of coding and ncRNAs involved in muscle regeneration and discuss their potential as therapeutic tools. Small ncRNAs, such as miRNA and siRNA, have been successfully applied in the treatment of several diseases. The use of longer molecules, such as lncRNA and circRNA, is less advanced. However, based on the peculiar properties discussed below, they represent an innovative pool of RNA biomarkers and possible targets of clinical value.

Deciphering the chromatin organization and dynamics for muscle stem cell function

The chromatin landscape represents a critical regulatory layer for precise transcriptional control. Chromosome architecture restrains the physical access to the DNA elements and is one of the determinants that specifies cell identity. Adult stem cells possess the unique ability to differentiate into a specific lineage. One of the underexplored areas in skeletal muscle biology is the molecular mechanism guiding the chromatin organization changes in muscle stem cell specification, myogenic determination, and differentiation. In this review, we focus on the regulatory network guiding the progression of muscle stem cells to differentiated progeny. We summarize recent findings regarding the mechanisms directing myogenic cell fate decision and differentiation, with a particular focus on three-dimensional chromosome architecture and long noncoding RNA-associated chromatin accessibility changes.

LINC02381 contributes to cell proliferation and hinders cell apoptosis in glioma by transcriptionally enhancing CBX5

Glioma, featured with high incidence and low survival rate, is the most common type of primary brain tumor, severely affecting human life worldwide. LINC02381 is an interesting lncRNA functioning as oncogenic lncRNA in some cancers but as tumor-suppressor in others, but no report demonstrates its association with and function in glioma. Intriguingly, we found in a bioinformatics website LncRNADisease that LINC02381 was closely related to malignant glioma, so this study aimed to figure out the expression and function of LINC02381 in glioma. By RT-qPCR, we confirmed LINC02381 upregulation in glioma cells. Functional experiments demonstrated that LINC02381 knockdown repressed glioma cell proliferation and induced apoptosis. Boinformatics tools and RT-qPCR revealed the positive correlation between LINC02381 and CBX5 in glioma cells. More importantly, we confirmed that LINC02381 could interact and work synergistically with CEBPβ to bind to CBX5 promoter and activate CBX5 transcriptionally. Additionally, rescue experiments indicated that CBX5 up-regulation reversed the decline in cell proliferation and the augment in cell apoptosis caused by LINC02381 knockdown. To conclude, LINC02381 could facilitate CBX5 transcription via interaction with CEBPβ, thus exerting its oncogenic role in glioma cells, which could contribute to better understanding of glioma.

Prdm16-mediated H3K9 methylation controls fibro-adipogenic progenitors identity during skeletal muscle repair

H3K9 methylation maintains cell identity orchestrating stable silencing and anchoring of alternate fate genes within the heterochromatic compartment underneath the nuclear lamina (NL). However, how cell type-specific genomic regions are specifically targeted to the NL is still elusive. Using fibro-adipogenic progenitors (FAPs) as a model, we identified Prdm16 as a nuclear envelope protein that anchors H3K9-methylated chromatin in a cell-specific manner. We show that Prdm16 mediates FAP developmental capacities by orchestrating lamina-associated domain organization and heterochromatin sequestration at the nuclear periphery. We found that Prdm16 localizes at the NL where it cooperates with the H3K9 methyltransferases G9a/GLP to mediate tethering and silencing of myogenic genes, thus repressing an alternative myogenic fate in FAPs. Genetic and pharmacological disruption of this repressive pathway confers to FAP myogenic competence, preventing fibro-adipogenic degeneration of dystrophic muscles. In summary, we reveal a druggable mechanism of heterochromatin perinuclear sequestration exploitable to reprogram FAPs in vivo.

Cis-Repression of Foxq1 Expression Affects Foxf2-Mediated Gene Expression in Palate Development

Disruption of FOXF2, encoding a member of the Forkhead family transcription factors, has been associated with cleft palate in humans and mice. FOXF2 is located in a conserved gene cluster containing FOXQ1, FOXF2, and FOXC1. We found that expression of Foxq1 is dramatically upregulated in the embryonic palatal mesenchyme in Foxf2 ^-/- mouse embryos. We show here that the Foxf2 promoter-deletion mutation caused dramatically increased expression of the cis-linked Foxq1 allele but had little effect on the Foxq1 allele in trans. We analyzed effects of the Foxf2 mutation on the expression of other neighboring genes and compared those effects with the chromatin domain structure and recently identified enhancer-promoter associations as well as H3K27ac ChIP-seq data. We show that the Foxf2 mutation resulted in significantly increased expression of the Foxq1 and Exoc2 genes located in the same topologically associated domain with Foxf2 but not the expression of the Foxc1 and Gmds genes located in the adjacent chromatin domain. We inactivated the Foxq1 gene in mice homozygous for a Foxf2 conditional allele using CRISPR genome editing and generated (Foxf2/Foxq1)^+/- mice with loss-of-function mutations in Foxf2 and Foxq1 in cis. Whereas the (Foxf2/Foxq1)^-/- mice exhibited cleft palate at birth similar as in the Foxf2 ^-/- mice, systematic expression analyses of a large number of Foxf2-dependent genes revealed that the (Foxf2/Foxq1)^-/- embryos exhibited distinct effects on the domain-specific expression of several important genes, including Foxf1, Shox2, and Spon1, in the developing palatal shelves compared with Foxf2 ^-/- embryos. These results identify a novel cis-regulatory effect of the Foxf2 mutation and demonstrate that cis-regulation of Foxq1 contributed to alterations in palatal gene expression in Foxf2 ^-/- embryos. These results have important implications for interpretation of results and mechanisms from studies of promoter- or gene-deletion alleles. In addition, the unique mouse lines generated in this study provide a valuable resource for understanding the cross-regulation and combinatorial functions of the Foxf2 and Foxq1 genes in development and disease.