Dynamic enhancer landscapes in human craniofacial development.

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Novel Association of the NOTCH Pathway Regulator MIB1 Gene With the Development of Bicuspid Aortic Valve.

Importance: Nonsyndromic bicuspid aortic valve (nsBAV) is the most common congenital heart valve malformation. BAV has a heritable component, yet only a few causative genes have been identified; understanding BAV genetics is a key point in developing personalized medicine. Objective: To identify a new gene for nsBAV. Design, Setting, and Participants: This was a comprehensive, multicenter, genetic association study based on candidate gene prioritization in a familial cohort followed by rare and common association studies in replication cohorts. Further validation was done using in vivo mice models. Study data were analyzed from October 2019 to October 2022. Three cohorts of patients with BAV were included in the study: (1) the discovery cohort was a large cohort of inherited cases from 29 pedigrees of French and Israeli origin; (2) the replication cohort 1 for rare variants included unrelated sporadic cases from various European ancestries; and (3) replication cohort 2 was a second validation cohort for common variants in unrelated sporadic cases from Europe and the US. Main Outcomes and Measures: To identify a candidate gene for nsBAV through analysis of familial cases exome sequencing and gene prioritization tools. Replication cohort 1 was searched for rare and predicted deleterious variants and genetic association. Replication cohort 2 was used to investigate the association of common variants with BAV. Results: A total of 938 patients with BAV were included in this study: 69 (7.4%) in the discovery cohort, 417 (44.5%) in replication cohort 1, and 452 (48.2%) in replication cohort 2. A novel human nsBAV gene, MINDBOMB1 homologue MIB1, was identified. MINDBOMB1 homologue (MIB1) is an E3-ubiquitin ligase essential for NOTCH-signal activation during heart development. In approximately 2% of nsBAV index cases from the discovery and replication 1 cohorts, rare MIB1 variants were detected, predicted to be damaging, and were significantly enriched compared with population-based controls (2% cases vs 0.9% controls; P = .03). In replication cohort 2, MIB1 risk haplotypes significantly associated with nsBAV were identified (permutation test, 1000 repeats; P = .02). Two genetically modified mice models carrying Mib1 variants identified in our cohort showed BAV on a NOTCH1-sensitized genetic background. Conclusions and Relevance: This genetic association study identified the MIB1 gene as associated with nsBAV. This underscores the crucial role of the NOTCH pathway in the pathophysiology of BAV and its potential as a target for future diagnostic and therapeutic intervention.

Cell Type- and Tissue-specific Enhancers in Craniofacial Development.

The genetic basis of craniofacial birth defects and general variation in human facial shape remains poorly understood. Distant-acting transcriptional enhancers are a major category of non-coding genome function and have been shown to control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development1-3. However, a lack of accurate maps of the genomic location and cell type-specific in vivo activities of all craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combined histone modification and chromatin accessibility profiling from different stages of human craniofacial development with single-cell analyses of the developing mouse face to create a comprehensive catalogue of the regulatory landscape of facial development at tissue- and single cell-resolution. In total, we identified approximately 14,000 enhancers across seven developmental stages from weeks 4 through 8 of human embryonic face development. We used transgenic mouse reporter assays to determine the in vivo activity patterns of human face enhancers predicted from these data. Across 16 in vivo validated human enhancers, we observed a rich diversity of craniofacial subregions in which these enhancers are active in vivo. To annotate the cell type specificities of human-mouse conserved enhancers, we performed single-cell RNA-seq and single-nucleus ATAC-seq of mouse craniofacial tissues from embryonic days e11.5 to e15.5. By integrating these data across species, we find that the majority (56%) of human craniofacial enhancers are functionally conserved in mice, providing cell type- and embryonic stage-resolved predictions of their in vivo activity profiles. Using retrospective analysis of known craniofacial enhancers in combination with single cell-resolved transgenic reporter assays, we demonstrate the utility of these data for predicting the in vivo cell type specificity of enhancers. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Topologically associating domain boundaries are required for normal genome function.

Topologically associating domain (TAD) boundaries partition the genome into distinct regulatory territories. Anecdotal evidence suggests that their disruption may interfere with normal gene expression and cause disease phenotypes1-3, but the overall extent to which this occurs remains unknown. Here we demonstrate that targeted deletions of TAD boundaries cause a range of disruptions to normal in vivo genome function and organismal development. We used CRISPR genome editing in mice to individually delete eight TAD boundaries (11-80 kb in size) from the genome. All deletions examined resulted in detectable molecular or organismal phenotypes, which included altered chromatin interactions or gene expression, reduced viability, and anatomical phenotypes. We observed changes in local 3D chromatin architecture in 7 of 8 (88%) cases, including the merging of TADs and altered contact frequencies within TADs adjacent to the deleted boundary. For 5 of 8 (63%) loci examined, boundary deletions were associated with increased embryonic lethality or other developmental phenotypes. For example, a TAD boundary deletion near Smad3/Smad6 caused complete embryonic lethality, while a deletion near Tbx5/Lhx5 resulted in a severe lung malformation. Our findings demonstrate the importance of TAD boundary sequences for in vivo genome function and reinforce the critical need to carefully consider the potential pathogenicity of noncoding deletions affecting TAD boundaries in clinical genetics screening.

Is Bicuspid Aortic Valve Morphology Genetically Determined? A Family-Based Study.

Bicuspid aortic valve (BAV) is a common congenital heart disease, with a 10-fold higher prevalence in first-degree relatives. BAV has different phenotypes based on the morphology of cusp fusion. These phenotypes are associated with different clinical courses and prognoses. Currently, the determinants of the valve phenotype are unknown. In this study we evaluated the role of genetics using familial cohorts. Patients with BAV and their first-degree relatives were evaluated by echocardiography. The concordance in BAV phenotype between pairs of family members was calculated and compared with the concordance expected by chance. We then performed a systematic literature review to identify additional reports and calculated the overall concordance rate. During the study period, 70 cases from 31 families and 327 sporadic cases were identified. BAV was diagnosed in 14% of the screened relatives. The proportions of the morphologies identified was: 12.3% for type 0, 66.2% for type 1-LR, 15.4% for type 1-RN, 4.6% for type 1-NL, and 1.5% for type 2. For the assessment of morphologic concordance, we included 120 pairs of first-degree relatives with BAV from our original cohort and the literature review. Concordance was found only in 62% of the pairs which was not significantly higher than expected by chance. In conclusion, our finding demonstrates intrafamilial variability in BAV morphology, suggesting that morphology is determined by factors other than Mendelian genetics. As prognosis differs by morphology, our findings may suggest that clinical outcomes may vary even between first-degree relatives.

The FABRIC Cancer Portal: A Ranked Catalogue of Gene Selection in Tumors Over the Human Coding Genome.

Contemporary catalogues of cancer driver genes rely primarily on high mutation rates as evidence for gene selection in tumors. Here, we present The Functional Alteration Bias Recovery In Coding-regions Cancer Portal, a comprehensive catalogue of gene selection in cancer based purely on the biochemical functional effects of mutations at the protein level. Gene selection in the portal is quantified by combining genomics data with rich proteomic annotations. Genes are ranked according to the strength of evidence for selection in tumor, based on rigorous and robust statistics. The portal covers the entire human coding genome (∼18,000 protein-coding genes) across 33 cancer types and pan-cancer. It includes a selected set of cross-references to the most relevant resources providing genomics, proteomics, and cancer-related information. We showcase the portal with known and overlooked cancer genes, demonstrating the utility of the portal via its simple visual interface, which allows users to pivot between gene-centric and cancer type views. The portal is available at fabric-cancer.huji.ac.il. SIGNIFICANCE: A new cancer portal quantifies and presents gene selection in tumor over the entire human coding genome across 33 cancer types and pan-cancer.

Comprehensive In Vivo Interrogation Reveals Phenotypic Impact of Human Enhancer Variants.

Establishing causal links between non-coding variants and human phenotypes is an increasing challenge. Here, we introduce a high-throughput mouse reporter assay for assessing the pathogenic potential of human enhancer variants in vivo and examine nearly a thousand variants in an enhancer repeatedly linked to polydactyly. We show that 71% of all rare non-coding variants previously proposed as causal lead to reporter gene expression in a pattern consistent with their pathogenic role. Variants observed to alter enhancer activity were further confirmed to cause polydactyly in knockin mice. We also used combinatorial and single-nucleotide mutagenesis to evaluate the in vivo impact of mutations affecting all positions of the enhancer and identified additional functional substitutions, including potentially pathogenic variants hitherto not observed in humans. Our results uncover the functional consequences of hundreds of mutations in a phenotype-associated enhancer and establish a widely applicable strategy for systematic in vivo evaluation of human enhancer variants.

Emergence of online communities: Empirical evidence and theory.

Online communities, which have become an integral part of the day-to-day life of people and organizations, exhibit much diversity in both size and activity level; some communities grow to a massive scale and thrive, whereas others remain small, and even wither. In spite of the important role of these proliferating communities, there is limited empirical evidence that identifies the dominant factors underlying their dynamics. Using data collected from seven large online platforms, we observe a relationship between online community size and its activity which generally repeats itself across platforms: First, in most platforms, three distinct activity regimes exist-one of low-activity and two of high-activity. Further, we find a sharp activity phase transition at a critical community size that marks the shift between the first and the second regime in six out of the seven online platforms. Essentially, we argue that it is around this critical size that sustainable interactive communities emerge. The third activity regime occurs above a higher characteristic size in which community activity reaches and remains at a constant and higher level. We find that there is variance in the steepness of the slope of the second regime, that leads to the third regime of saturation, but that the third regime is exhibited in six of the seven online platforms. We propose that the sharp activity phase transition and the regime structure stem from the branching property of online interactions.

Centralized nucleation in online networks leads to high social inequality.

Even though the heterogeneity of centrality in social networks is well documented, its role and effect on network stability in real life remains unclear. The literature roughly suggests that network structure is such that networks have an "inner" highly-connected nucleus and, in contrast, sparse outer shells. But to what extent is the existence of this nucleus crucial for the survival of a network? To what extent is the outer shells' much larger population essential to the longevity of the network? Furthermore, as a network grows and forms, theoretically speaking, network structure should be dependent on the patterns of change of degree centrality, i.e., social mobility between centrality shells. What is the role of social mobility in the formation of the nucleus-to-periphery profile, and is it related to network lifetime? Here, we explore these questions using data collected covering over a decade of activity from more than 10, 000 networked communities, with more than 134,000 users. We find that: (i) social mobility is, on average, negative but that, (ii) the higher the social mobility of the members of the network, the more stable and long-living the network is. Further, (iii) the network is, indeed, composed of two phases - a large but ephemeral sparsely connected "cloud" of actors, that nucleates around a highly stable nucleus of users. Lastly, (iv) networked communities which maintain a specific nucleus-to-periphery ratio η, i.e., a ratio of the size of the nucleus to periphery of around η = 1 4 , have a greater chance of survival. We find that deviations from this nucleus-to-periphery ratio predict a collapse of network activity, especially in the case of younger communities.

Gene ORGANizer: linking genes to the organs they affect.

One of the biggest challenges in studying how genes work is understanding their effect on the physiology and anatomy of the body. Existing tools try to address this using indirect features, such as expression levels and biochemical pathways. Here, we present Gene ORGANizer (geneorganizer.huji.ac.il), a phenotype-based tool that directly links human genes to the body parts they affect. It is built upon an exhaustive curated database that links >7000 genes to ∼150 anatomical parts using >150 000 gene-organ associations. The tool offers user-friendly platforms to analyze the anatomical effects of individual genes, and identify trends within groups of genes. We demonstrate how Gene ORGANizer can be used to make new discoveries, showing that chromosome X is enriched with genes affecting facial features, that positive selection targets genes with more constrained phenotypic effects, and more. We expect Gene ORGANizer to be useful in a variety of evolutionary, medical and molecular studies aimed at understanding the phenotypic effects of genes.