Network-based identification of critical regulators as putative drivers of human cleft lip

Single-cell multiomics decodes regulatory programs for mouse secondary palate development

Perturbations in gene regulation during palatogenesis can lead to cleft palate, which is among the most common congenital birth defects. Here, we perform single-cell multiome sequencing and profile chromatin accessibility and gene expression simultaneously within the same cells (n = 36,154) isolated from mouse secondary palate across embryonic days (E) 12.5, E13.5, E14.0, and E14.5. We construct five trajectories representing continuous differentiation of cranial neural crest-derived multipotent cells into distinct lineages. By linking open chromatin signals to gene expression changes, we characterize the underlying lineage-determining transcription factors. In silico perturbation analysis identifies transcription factors SHOX2 and MEOX2 as important regulators of the development of the anterior and posterior palate, respectively. In conclusion, our study charts epigenetic and transcriptional dynamics in palatogenesis, serving as a valuable resource for further cleft palate research.

MicroRNA-Transcription factor regulatory networks in the early strobilar development of Echinococcus granulosus protoscoleces

BACKGROUND

Echinococcus granulosus sensu lato has a complex developmental biology with a variety of factors relating to both intermediate and final hosts. To achieve maximum parasite adaptability, the development of the cestode is dependent on essential changes in transcript regulation. Transcription factors (TFs) and miRNAs are known as master regulators that affect the expression of downstream genes through a wide range of metabolic and signaling pathways. In this study, we aimed to develop a regulatory miRNA-Transcription factor (miRNA-TF) network across early developmental stages of E. granulosus protoscoleces by performing in silico analysis, and to experimentally validate TFs expression in protoscoleces obtained from in vitro culture, and from in vivo experiments.

RESULTS

We obtained list of 394 unique E. granulosus TFs and matched them with 818 differentially expressed genes which identified 41 predicted TFs with differential expression. These TFs were used to predict the potential targets of 31 differentially expressed miRNAs. As a result, eight miRNAs and eight TFs were found, and the predicted network was constructed using Cytoscape. At least four miRNAs (egr-miR-124a, egr-miR-124b-3p, egr-miR-745-3p, and egr-miR-87-3p) and their corresponding differentially expressed TFs (Zinc finger protein 45, Early growth response protein 3, Ecdysone induced protein 78c and ETS transcription factor elf 2) were highlighted in this investigation. The expression of predicted differentially expressed TFs obtained from in vitro and in vivo experiments, were experimentally validated by quantitative polymerase chain reaction. This confirmed findings of RNA-seq data.

CONCLUSION

miRNA-TF networks presented in this study control some of the most important metabolic and signaling pathways in the development and life cycle of E. granulosus, providing a potential approach for disrupting the early hours of dog infection and preventing the development of the helminth in the final host.

Spatiotemporal MicroRNA-Gene Expression Network Related to Orofacial Clefts

Craniofacial structures change dynamically in morphology during development through the coordinated regulation of various cellular molecules. However, it remains unclear how these complex mechanisms are regulated in a spatiotemporal manner. Here we applied natural cubic splines to model gene and microRNA (miRNA) expression from embryonic day (E) 10.5 to E14.5 in the proximal and distal regions of the maxillary processes to identify spatiotemporal patterns of gene and miRNA expression, followed by constructing corresponding regulatory networks. Three major groups of differentially expressed genes (DEGs) were identified, including 3,927 temporal, 314 spatial, and 494 spatiotemporal DEGs. Unsupervised clustering further resolved these spatiotemporal DEGs into 8 clusters with distinct expression patterns. Interestingly, we found 2 clusters of differentially expressed miRNAs: 1 had 80 miRNAs monotonically decreasing and the other had 97 increasing across developmental stages. To evaluate the phenotypic relevance of these DEGs during craniofacial development, we integrated data from the CleftGeneDB database and constructed the regulatory networks of genes related to orofacial clefts. Our analysis revealed 2 hub miRNAs, mmu-miR-325-3p and mmu-miR-384-5p, that repressed cleft-related genes Adamts3, Runx2, Fgfr2, Acvr1, and Edn2, while their expression increased over time. On the contrary, 2 hub miRNAs, mmu-miR-218-5p and mmu-miR-338-5p, repressed cleft-related genes Pbx2, Ermp1, Snai1, Tbx2, and Bmi1, while their expression decreased over time. Our experiments indicated that these miRNA mimics significantly inhibited cell proliferation in mouse embryonic palatal mesenchymal (MEPM) cells and O9-1 cells through the regulation of genes associated with cleft palate and validated the role of our regulatory networks in orofacial clefts. To facilitate interactive exploration of these data, we developed a user-friendly web tool to visualize the gene and miRNA expression patterns across developmental stages, as well as the regulatory networks (https://fyan.shinyapps.io/facebase_shiny/). Taken together, our results provide a valuable resource that serves as a reference map for future research in craniofacial development.

FaceBase: A Community-Driven Hub for Data-Intensive Research

The FaceBase Consortium, funded by the National Institute of Dental and Craniofacial Research of the National Institutes of Health, was established in 2009 with the recognition that dental and craniofacial research are increasingly data-intensive disciplines. Data sharing is critical for the validation and reproducibility of results as well as to enable reuse of data. In service of these goals, data ought to be FAIR: Findable, Accessible, Interoperable, and Reusable. The FaceBase data repository and educational resources exemplify the FAIR principles and support a broad user community including researchers in craniofacial development, molecular genetics, and genomics. FaceBase demonstrates that a model in which researchers "self-curate" their data can be successful and scalable. We present the results of the first 2.5 y of FaceBase's operations as an open community and summarize the data sets published during this period. We then describe a research highlight from work on the identification of regulatory networks and noncoding RNAs involved in cleft lip with/without cleft palate that both used and in turn contributed new findings to publicly available FaceBase resources. Collectively, FaceBase serves as a dynamic and continuously evolving resource to facilitate data-intensive research, enhance data reproducibility, and perform deep phenotyping across multiple species in dental and craniofacial research.

Identification of Novel Genomic Variations in Susceptibility to Nonsyndromic Cleft Lip and Palate Patients

BACKGROUND

Nonsyndromic cleft lip with or without palate (NSCL/P) is a multifactorial and common birth malformation caused by genetic and environmental factors, as well as by teratogens. Genome-wide association studies found genetic variations with modulatory effects of NSCL/P formation in Chinese and Iranian populations. We aimed to identify the susceptibility of single-nucleotide polymorphisms (SNPs) to nonsyndromic cleft lip with or without palate in the Indian population.

MATERIAL AND METHODS

The present study was conducted on NSCL/P cases and controls. Genomic DNA was extracted from peripheral blood and Axiom- Precision Medicine Research Array (PMRA) was performed. The Axiom-PMRA covers 902,527 markers and several thousand novel risk variants. Quality control-passed samples were included for candidate genetic variation identification, gene functional enrichment, and pathway and network analysis.

RESULTS

The genome-wide association study identified fourteen novel candidate gene SNPs that showed the most significant association with the risk of NSCL/P, and eight were predicted to have regulatory sequences.

CONCLUSION

The GWAS study showed novel candidate genetic variations in NSCL/P formations. These findings contribute to the understanding of genetic predisposition to nonsyndromic cleft lip with or without palate.

Evaluating the Expression of Candidate Homeobox Genes and Their Role in Local-Site Inflammation in Mucosal Tissue Obtained from Children with Non-Syndromic Cleft Lip and Palate

Craniofacial development including palatogenesis is a complex process which requires an orchestrated and spatiotemporal expression of various genes and factors for proper embryogenesis and organogenesis. One such group of genes essential for craniofacial development is the homeobox genes, transcriptional factors that are commonly associated with congenital abnormalities. Amongst these genes, DLX4, HOXB3, and MSX2 have been recently shown to be involved in the etiology of non-syndromic cleft lip and palate. Hence, we investigated the gene and protein expression of these genes in normal and cleft affected mucosal tissue obtained from 22 children, along with analyzing their role in promoting local-site inflammation using NF-κB. Additionally, we investigated the role of PTX3, which plays a critical role in tissue remodeling and wound repair. We found a residual gene and protein expression of DLX4 in cleft mucosa, although no differences in gene expression levels of HOXB3 and MSX2 were noted. However, a significant increase in protein expression for these genes was noted in the cleft mucosa (p < 0.05), indicating increased cellular proliferation. This was coupled with a significant increase in NF-κB protein expression in cleft mucosa (p < 0.05), highlighting the role of these genes in promotion of pro-inflammatory environment. Finally, no differences in gene expression of PTX3 were noted.

miRGTF-net: Integrative miRNA-gene-TF network analysis reveals key drivers of breast cancer recurrence

Analysis of regulatory networks is a powerful framework for identification and quantification of intracellular interactions. We introduce miRGTF-net, a novel tool for construction of miRNA-gene-TF networks. We consider multiple transcriptional and post-transcriptional interaction types, including regulation of gene and miRNA expression by transcription factors, gene silencing by miRNAs, and co-expression of host genes with their intronic miRNAs. The underlying algorithm uses information on experimentally validated interactions as well as integrative miRNA/mRNA expression profiles in a given set of samples. The latter ensures simultaneous tissue-specificity and biological validity of interactions. We applied miRGTF-net to paired miRNA/mRNA-sequencing data of breast cancer samples from The Cancer Genome Atlas (TCGA). Together with topological analysis of the constructed network we showed that considered players can form reliable prognostic gene signatures for ER-positive breast cancer. A number of signatures demonstrated remarkably high accuracy on transcriptomic data obtained by both microarrays and RNA sequencing from several independent patient cohorts. Furthermore, an essential part of prognostic genes were identified as direct targets of transcription factor E2F1. The putative interplay between estrogen receptor alpha and E2F1 was suggested as a potential recurrence factor in patients treated with tamoxifen. Source codes of miRGTF-net are available at GitHub (https://github.com/s-a-nersisyan/miRGTF-net).

A developmental stage-specific network approach for studying dynamic co-regulation of transcription factors and microRNAs during craniofacial development

Craniofacial development is regulated through dynamic and complex mechanisms that involve various signaling cascades and gene regulations. Disruption of such regulations can result in craniofacial birth defects. Here, we propose the first developmental stage-specific network approach by integrating two crucial regulators, transcription factors (TFs) and microRNAs (miRNAs), to study their co-regulation during craniofacial development. Specifically, we used TFs, miRNAs and non-TF genes to form feed-forward loops (FFLs) using genomic data covering mouse embryonic days E10.5 to E14.5. We identified key novel regulators (TFs Foxm1, Hif1a, Zbtb16, Myog, Myod1 and Tcf7, and miRNAs miR-340-5p and miR-129-5p) and target genes (Col1a1, Sgms2 and Slc8a3) expression of which changed in a developmental stage-dependent manner. We found that the Wnt-FoxO-Hippo pathway (from E10.5 to E11.5), tissue remodeling (from E12.5 to E13.5) and miR-129-5p-mediated Col1a1 regulation (from E10.5 to E14.5) might play crucial roles in craniofacial development. Enrichment analyses further suggested their functions. Our experiments validated the regulatory roles of miR-340-5p and Foxm1 in the Wnt-FoxO-Hippo subnetwork, as well as the role of miR-129-5p in the miR-129-5p-Col1a1 subnetwork. Thus, our study helps understand the comprehensive regulatory mechanisms for craniofacial development.

Critical microRNAs and regulatory motifs in cleft palate identified by a conserved miRNA-TF-gene network approach in humans and mice

Cleft palate (CP) is the second most common congenital birth defect. The etiology of CP is complicated, with involvement of various genetic and environmental factors. To investigate the gene regulatory mechanisms, we designed a powerful regulatory analytical approach to identify the conserved regulatory networks in humans and mice, from which we identified critical microRNAs (miRNAs), target genes and regulatory motifs (miRNA-TF-gene) related to CP. Using our manually curated genes and miRNAs with evidence in CP in humans and mice, we constructed miRNA and transcription factor (TF) co-regulation networks for both humans and mice. A consensus regulatory loop (miR17/miR20a-FOXE1-PDGFRA) and eight miRNAs (miR-140, miR-17, miR-18a, miR-19a, miR-19b, miR-20a, miR-451a and miR-92a) were discovered in both humans and mice. The role of miR-140, which had the strongest association with CP, was investigated in both human and mouse palate cells. The overexpression of miR-140-5p, but not miR-140-3p, significantly inhibited cell proliferation. We further examined whether miR-140 overexpression could suppress the expression of its predicted target genes (BMP2, FGF9, PAX9 and PDGFRA). Our results indicated that miR-140-5p overexpression suppressed the expression of BMP2 and FGF9 in cultured human palate cells and Fgf9 and Pdgfra in cultured mouse palate cells. In summary, our conserved miRNA-TF-gene regulatory network approach is effective in detecting consensus miRNAs, motifs, and regulatory mechanisms in human and mouse CP.

MicroRNA and transcription factor co-regulatory networks and subtype classification of seminoma and non-seminoma in testicular germ cell tumors

Recent studies have revealed that feed-forward loops (FFLs) as regulatory motifs have synergistic roles in cellular systems and their disruption may cause diseases including cancer. FFLs may include two regulators such as transcription factors (TFs) and microRNAs (miRNAs). In this study, we extensively investigated TF and miRNA regulation pairs, their FFLs, and TF-miRNA mediated regulatory networks in two major types of testicular germ cell tumors (TGCT): seminoma (SE) and non-seminoma (NSE). Specifically, we identified differentially expressed mRNA genes and miRNAs in 103 tumors using the transcriptomic data from The Cancer Genome Atlas. Next, we determined significantly correlated TF-gene/miRNA and miRNA-gene/TF pairs with regulation direction. Subsequently, we determined 288 and 664 dysregulated TF-miRNA-gene FFLs in SE and NSE, respectively. By constructing dysregulated FFL networks, we found that many hub nodes (12 out of 30 for SE and 8 out of 32 for NSE) in the top ranked FFLs could predict subtype-classification (Random Forest classifier, average accuracy ≥90%). These hub molecules were validated by an independent dataset. Our network analysis pinpointed several SE-specific dysregulated miRNAs (miR-200c-3p, miR-25-3p, and miR-302a-3p) and genes (EPHA2, JUN, KLF4, PLXDC2, RND3, SPI1, and TIMP3) and NSE-specific dysregulated miRNAs (miR-367-3p, miR-519d-3p, and miR-96-5p) and genes (NR2F1 and NR2F2). This study is the first systematic investigation of TF and miRNA regulation and their co-regulation in two major TGCT subtypes.

The International Conference on Intelligent Biology and Medicine (ICIBM) 2018: genomics meets medicine

During June 10–12, 2018, the International Conference on Intelligent Biology and Medicine (ICIBM 2018) was held in Los Angeles, California, USA. The conference included 11 scientific sessions, four tutorials, one poster session, four keynote talks and four eminent scholar talks that covered a wide range of topics ranging from 3D genome structure analysis and visualization, next generation sequencing analysis, computational drug discovery, medical informatics, cancer genomics to systems biology. While medical genomics has always been a main theme in ICIBM, this year we for the first time organized the BMC Medical Genomics Supplement for ICIBM. Here, we describe 15 ICIBM papers selected for publishing in BMC Medical Genomics.