Homologous recombination (HR) is a template-driven repair pathway that mends DNA double-stranded breaks (DSBs), and thus helps to maintain genome stability. The RAD51 recombinase facilitates DNA joint formation during HR, but to accomplish this task, RAD51 must be loaded onto the single-stranded DNA. DSS1, a candidate gene for split hand/split foot syndrome, provides the ability to recognize RPA-coated ssDNA to the tumor suppressor BRCA2, which is complexed with RAD51. Together BRCA2-DSS1 displace RPA and load RAD51 onto the ssDNA. In addition, the BRCA2 interacting protein BCCIP normally colocalizes with chromatin bound BRCA2, and upon DSB induction, RAD51 colocalizes with BRCA2-BCCIP foci. Down-regulation of BCCIP reduces DSB repair and disrupts BRCA2 and RAD51 foci formation. While BCCIP is known to interact with BRCA2, the relationship between BCCIP and RAD51 is not known. In this study, we investigated the biochemical role of the β-isoform of BCCIP in relation to the RAD51 recombinase. We demonstrate that BCCIPβ binds DNA and physically and functionally interacts with RAD51 to stimulate its homologous DNA pairing activity. Notably, this stimulatory effect is not the result of RAD51 nucleoprotein filament stabilization; rather, we demonstrate that BCCIPβ induces a conformational change within the RAD51 filament that promotes release of ADP to help maintain an active presynaptic filament. Our findings reveal a functional role for BCCIPβ as a RAD51 accessory factor in HR.

Prdm14 is a sequence-specific transcriptional regulator of embryonic stem cell (ESC) pluripotency and primordial germ cell (PGC) formation. It exerts its function, at least in part, through repressing genes associated with epigenetic modification and cell differentiation. Here, we show that this repressive function is mediated through an ETO-family co-repressor Mtgr1, which tightly binds to the pre-SET/SET domains of Prdm14 and co-occupies its genomic targets in mouse ESCs. We generated two monobodies, synthetic binding proteins, targeting the Prdm14 SET domain and demonstrate their utility, respectively, in facilitating crystallization and structure determination of the Prdm14-Mtgr1 complex, or as genetically encoded inhibitor of the Prdm14-Mtgr1 interaction. Structure-guided point mutants and the monobody abrogated the Prdm14-Mtgr1 association and disrupted Prdm14's function in mESC gene expression and PGC formation in vitro. Altogether, our work uncovers the molecular mechanism underlying Prdm14-mediated repression and provides renewable reagents for studying and controlling Prdm14 functions.

The FaceBase Consortium, funded by the National Institute of Dental and Craniofacial Research, National Institutes of Health, is designed to accelerate understanding of craniofacial developmental biology by generating comprehensive data resources to empower the research community, exploring high-throughput technology, fostering new scientific collaborations among researchers and human/computer interactions, facilitating hypothesis-driven research and translating science into improved health care to benefit patients. The resources generated by the FaceBase projects include a number of dynamic imaging modalities, genome-wide association studies, software tools for analyzing human facial abnormalities, detailed phenotyping, anatomical and molecular atlases, global and specific gene expression patterns, and transcriptional profiling over the course of embryonic and postnatal development in animal models and humans. The integrated data visualization tools, faceted search infrastructure, and curation provided by the FaceBase Hub offer flexible and intuitive ways to interact with these multidisciplinary data. In parallel, the datasets also offer unique opportunities for new collaborations and training for researchers coming into the field of craniofacial studies. Here, we highlight the focus of each spoke project and the integration of datasets contributed by the spokes to facilitate craniofacial research.

Histone H3 lysine 4 (K4) methylation has been linked to the transcriptional activation in a variety of eukaryotic species. Here we show that a common component of MLL1, MLL2, and hSet1 H3 K4 methyltransferase complexes, the WD40-repeat protein WDR5, directly associates with histone H3 di- and trimethylated at K4 and with H3-K4-dimethylated nucleosomes. WDR5 is required for binding of the methyltransferase complex to the K4-dimethylated H3 tail as well as for global H3 K4 trimethylation and HOX gene activation in human cells. WDR5 is essential for vertebrate development, in that WDR5-depleted X. laevis tadpoles exhibit a variety of developmental defects and abnormal spatial Hox gene expression. Our results are the first demonstration that a WD40-repeat protein acts as a module for recognition of a specific histone modification and suggest a mechanism for reading and writing an epigenetic mark for gene activation.

Estimates of individual-level genomic ancestry are routinely used in human genetics, and related fields. The analysis of population structure and genomic ancestry can yield insights in terms of modern and ancient populations, allowing us to address questions regarding admixture, and the numbers and identities of the parental source populations. Unrecognized population structure is also an important confounder to correct for in genome-wide association studies. However, it remains challenging to work with heterogeneous datasets from multiple studies collected by different laboratories with diverse genotyping and imputation protocols. This work presents a new approach and an accompanying open-source toolbox that facilitates a robust integrative analysis for population structure and genomic ancestry estimates for heterogeneous datasets. We show robustness against individual outliers and different protocols for the projection of new samples into a reference ancestry space, and the ability to reveal and adjust for population structure in a simulated case-control admixed population. Given that visually evident and easily recognizable patterns of human facial characteristics co-vary with genomic ancestry, and based on the integration of three different sources of genome data, we generate average 3D faces to illustrate genomic ancestry variations within the 1,000 Genome project and for eight ancient-DNA profiles, respectively.

Transposable elements, also known as transposons, are now recognized not only as parasitic DNA, the spread of which in the genome must be controlled by the host, but also as major players in genome evolution and regulation. Long interspersed element-1 (LINE-1, also known as L1), the only currently autonomous mobile transposon in humans, occupies 17% of the genome and generates inter- and intra-individual genetic variation, in some cases resulting in disease. However, how L1 activity is controlled and the function of L1s in host gene regulation are not completely understood. Here we use CRISPR-Cas9 screening strategies in two distinct human cell lines to provide a genome-wide survey of genes involved in the control of L1 retrotransposition. We identify functionally diverse genes that either promote or restrict L1 retrotransposition. These genes, which are often associated with human diseases, control the L1 life cycle at the transcriptional or the post-transcriptional level in a manner that can depend on the endogenous L1 nucleotide sequence, underscoring the complexity of L1 regulation. We further investigate the restriction of L1 by the protein MORC2 and by the human silencing hub (HUSH) complex subunits MPP8 and TASOR. HUSH and MORC2 can selectively bind evolutionarily young, full-length L1s located within transcriptionally permissive euchromatic environments, and promote deposition of histone H3 Lys9 trimethylation (H3K9me3) for transcriptional silencing. Notably, these silencing events often occur within introns of transcriptionally active genes, and lead to the downregulation of host gene expression in a HUSH-, MORC2-, and L1-dependent manner. Together, these results provide a rich resource for studies of L1 retrotransposition, elucidate a novel L1 restriction pathway and illustrate how epigenetic silencing of transposable elements rewires host gene expression programs.

CHARGE syndrome is caused by heterozygous mutations in the chromatin remodeler, CHD7, and is characterized by a set of malformations that, on clinical grounds, were historically postulated to arise from defects in neural crest formation during embryogenesis. To better delineate neural crest defects in CHARGE syndrome, we generated induced pluripotent stem cells (iPSCs) from two patients with typical syndrome manifestations, and characterized neural crest cells differentiated in vitro from these iPSCs (iPSC-NCCs). We found that expression of genes associated with cell migration was altered in CHARGE iPSC-NCCs compared to control iPSC-NCCs. Consistently, CHARGE iPSC-NCCs showed defective delamination, migration and motility in vitro, and their transplantation in ovo revealed overall defective migratory activity in the chick embryo. These results support the historical inference that CHARGE syndrome patients exhibit defects in neural crest migration, and provide the first successful application of patient-derived iPSCs in modeling craniofacial disorders.

Many studies have suggested that developmental instability (DI) could lead to asymmetric development, otherwise known as fluctuating asymmetry (FA). Several attempts to unravel the biological meaning of FA have been made, yet the main step in estimating FA is to remove the effects of directional asymmetry (DA), which is defined as the average bilateral asymmetry at the population level. Here, we demonstrate in a multivariate context that the conventional method of DA correction does not adequately compensate for the effects of DA in other dimensions of asymmetry. This appears to be due to the presence of between-individual variation along the DA dimension. Consequently, we propose to decompose asymmetry into its different orthogonal dimensions, where we introduce a new measure of asymmetry, namely fluctuating directional asymmetry (F-DA). This measure describes individual variation in the dimension of DA, and can be used to adequately correct the asymmetry measurements for the presence of DA. We provide evidence that this measure can be useful in disentangling the different dimensions of asymmetry, and further studies on this measure can provide valuable insight into the underlying biological processes leading to these different asymmetry dimensions.

Dense surface registration, commonly used in computer science, could aid the biological sciences in accurate and comprehensive quantification of biological phenotypes. However, few toolboxes exist that are openly available, non-expert friendly, and validated in a way relevant to biologists. Here, we report a customizable toolbox for reproducible high-throughput dense phenotyping of 3D images, specifically geared towards biological use. Given a target image, a template is first oriented, repositioned, and scaled to the target during a scaled rigid registration step, then transformed further to fit the specific shape of the target using a non-rigid transformation. As validation, we use n = 41 3D facial images to demonstrate that the MeshMonk registration is accurate, with 1.26 mm average error, across 19 landmarks, between placements from manual observers and using the MeshMonk toolbox. We also report no variation in landmark position or centroid size significantly attributable to landmarking method used. Though validated using 19 landmarks, the MeshMonk toolbox produces a dense mesh of vertices across the entire surface, thus facilitating more comprehensive investigations of 3D shape variation. This expansion opens up exciting avenues of study in assessing biological shapes to better understand their phenotypic variation, genetic and developmental underpinnings, and evolutionary history.

Unaffected relatives of individuals with non-syndromic cleft lip with or without cleft palate (NSCL/P) show distinctive facial features. The presence of this facial endophenotype is potentially an expression of underlying genetic susceptibility to NSCL/P in the larger unselected population. To explore this hypothesis, we first partitioned the face into 63 partially overlapping regions representing global-to-local facial morphology and then defined endophenotypic traits by contrasting the 3D facial images from 264 unaffected parents of individuals with NSCL/P versus 3,171 controls. We observed distinct facial features between parents and controls across 59 global-to-local facial segments at nominal significance (p ≤ 0.05) and 52 segments at Bonferroni corrected significance (p < 1.2 × 10-3), respectively. Next, we quantified these distinct facial features as univariate traits in another dataset of 8,246 unaffected European individuals and performed a genome-wide association study. We identified 29 independent genetic loci that were associated (p < 5 × 10-8) with at least one of the tested endophenotypic traits, and nine genetic loci also passed the study-wide threshold (p < 8.47 × 10-10). Of the 29 loci, 22 were in proximity of loci previously associated with normal facial variation, 18 were near genes that show strong evidence in orofacial clefting (OFC), and another 10 showed some evidence in OFC. Additionally, polygenic risk scores for NSCL/P showed associations with the endophenotypic traits. This study thus supports the hypothesis of a shared genetic architecture of normal facial development and OFC.

Evidence from model organisms and clinical genetics suggests coordination between the developing brain and face, but the role of this link in common genetic variation remains unknown. We performed a multivariate genome-wide association study of cortical surface morphology in 19,644 individuals of European ancestry, identifying 472 genomic loci influencing brain shape, of which 76 are also linked to face shape. Shared loci include transcription factors involved in craniofacial development, as well as members of signaling pathways implicated in brain-face cross-talk. Brain shape heritability is equivalently enriched near regulatory regions active in either forebrain organoids or facial progenitors. However, we do not detect significant overlap between shared brain-face genome-wide association study signals and variants affecting behavioral-cognitive traits. These results suggest that early in embryogenesis, the face and brain mutually shape each other through both structural effects and paracrine signaling, but this interplay may not impact later brain development associated with cognitive function.

Pesticide exposure is associated with many long-term health outcomes; the potential underlying mechanisms are not well established for most associations. Epigenetic modifications, such as DNA methylation, may contribute. Individual pesticides may be associated with specific DNA methylation patterns but no epigenome-wide association study (EWAS) has evaluated methylation in relation to individual pesticides.

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how 'Coordinator', a long DNA motif comprised of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, while HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in shared regulation of genes involved in cell-type and positional identities, and ultimately shapes facial morphology and evolution.

Many factors influence human facial morphology, including genetics, age, nutrition, biomechanical forces, and endocrine factors. Moreover, facial features clearly differ between males and females, and these differences are driven primarily by the influence of sex hormones during growth and development. Specific genetic variants are known to influence circulating sex hormone levels in humans, which we hypothesize, in turn, affect facial features. In this study, we investigated the effects of testosterone-related genetic variants on facial morphology. We tested 32 genetic variants across 22 candidate genes related to levels of testosterone, sex hormone-binding globulin (SHGB) and dehydroepiandrosterone sulfate (DHEAS) in three cohorts of healthy individuals for which 3D facial surface images were available (Pittsburgh 3DFN, Penn State and ALSPAC cohorts; total n = 7418). Facial shape was described using a recently developed extension of the dense-surface correspondence approach, in which the 3D facial surface was partitioned into a set of 63 hierarchically organized modules. Each variant was tested against each of the facial surface modules in a multivariate genetic association-testing framework and meta-analyzed. Additionally, the association between these candidate SNPs and five facial ratios was investigated in the Pittsburgh 3DFN cohort. Two significant associations involving intronic variants of SHBG were found: both rs12150660 (p = 1.07E-07) and rs1799941 (p = 6.15E-06) showed an effect on mandible shape. Rs8023580 (an intronic variant of NR2F2-AS1) showed an association with the total and upper facial width to height ratios (p = 9.61E-04 and p = 7.35E-04, respectively). These results indicate that testosterone-related genetic variants affect normal-range facial morphology, and in particular, facial features known to exhibit strong sexual dimorphism in humans.

The theory of mind (ToM) is the ability to attribute mental states to others and is extremely important for social functioning. It has been widely examined in both behavioral and neuroimaging research, usually with the use of the many versions of the false-belief (FB) task. However, there is still not enough evidence from studies on the neurodevelopmental mechanisms of ToM mostly because of methodological limitations: lack of selectivity, mismatch of experimental and control tasks, and focusing on participants older than 6 years old. In the current study, we attempted to develop a computerized tool for ToM assessment suitable for both behavioral and neuriomaging testing in preschoolers. We designed a version of the classic change-of-location task with custom visuals and three fine-tuned conditions: FB, true-belief, and no-belief (NB). The usability of the task for further application in neurodevelopmental research was tested with three methods: first, behaviorally, with the use of a touch screen on a group of 75 children, followed by a functional MRI (fMRI) study on 13 adults, and a functional near-infrared spectroscopy (fNIRS) study on 19 preschool children. In line with our expectations, on the behavioral level, our task elicited the all-or-none performance in preschoolers. There was also a progression of performance with age in the FB condition. On the neural level, we observed the activation of structures involved in the ToM brain network in response to our task in both adults and children. The results therefore suggest that our task can be a useful tool for studying ToM development and its neural underpinnings.

Epithelial to mesenchymal transition (EMT) is a biological process involved in tissue morphogenesis and disease that causes dramatic changes in cell morphology, migration, proliferation, and gene expression. The retinal pigment epithelium (RPE), which supports the neural retina, can undergo EMT, producing fibrous epiretinal membranes (ERMs) associated with vision-impairing clinical conditions, such as macular pucker and proliferative vitreoretinopathy (PVR). We found that co-treatment with TGF-β and TNF-α (TNT) accelerates EMT in adult human RPE stem cell-derived RPE cell cultures. We captured the global epigenomic and transcriptional changes elicited by TNT treatment of RPE and identified putative active enhancers associated with actively transcribed genes, including a set of upregulated transcription factors that are candidate regulators. We found that the vitamin B derivative nicotinamide downregulates these key transcriptional changes, and inhibits and partially reverses RPE EMT, revealing potential therapeutic routes to benefit patients with ERM, macular pucker and PVR.

Facial morphology is highly variable, both within and among human populations, and a sizable portion of this variation is attributable to genetics. Previous genome scans have revealed more than 100 genetic loci associated with different aspects of normal-range facial variation. Most of these loci have been detected in Europeans, with few studies focusing on other ancestral groups. Consequently, the degree to which facial traits share a common genetic basis across diverse sets of humans remains largely unknown. We therefore investigated the genetic basis of facial morphology in an East African cohort. We applied an open-ended data-driven phenotyping approach to a sample of 2,595 3D facial images collected on Tanzanian children. This approach segments the face into hierarchically arranged, multivariate features that capture the shape variation after adjusting for age, sex, height, weight, facial size and population stratification. Genome scans of these multivariate shape phenotypes revealed significant (p < 2.5 × 10-8) signals at 20 loci, which were enriched for active chromatin elements in human cranial neural crest cells and embryonic craniofacial tissue, consistent with an early developmental origin of the facial variation. Two of these associations were in highly conserved regions showing craniofacial-specific enhancer activity during embryological development (5q31.1 and 12q21.31). Six of the 20 loci surpassed a stricter threshold accounting for multiple phenotypes with study-wide significance (p < 6.25 × 10-10). Cross-population comparisons indicated 10 association signals were shared with Europeans (seven sharing the same associated SNP), and facilitated fine-mapping of causal variants at previously reported loci. Taken together, these results may point to both shared and population-specific components to the genetic architecture of facial variation.

Transcriptional regulation exhibits extensive robustness, but human genetics indicates sensitivity to transcription factor (TF) dosage. Reconciling such observations requires quantitative studies of TF dosage effects at trait-relevant ranges, largely lacking so far. TFs play central roles in both normal-range and disease-associated variation in craniofacial morphology; we therefore developed an approach to precisely modulate TF levels in human facial progenitor cells and applied it to SOX9, a TF associated with craniofacial variation and disease (Pierre Robin sequence (PRS)). Most SOX9-dependent regulatory elements (REs) are buffered against small decreases in SOX9 dosage, but REs directly and primarily regulated by SOX9 show heightened sensitivity to SOX9 dosage; these RE responses partially predict gene expression responses. Sensitive REs and genes preferentially affect functional chondrogenesis and PRS-like craniofacial shape variation. We propose that such REs and genes underlie the sensitivity of specific phenotypes to TF dosage, while buffering of other genes leads to robust, nonlinear dosage-to-phenotype relationships.

The human face is complex and multipartite, and characterization of its genetic architecture remains challenging. Using a multivariate genome-wide association study meta-analysis of 8,246 European individuals, we identified 203 genome-wide-significant signals (120 also study-wide significant) associated with normal-range facial variation. Follow-up analyses indicate that the regions surrounding these signals are enriched for enhancer activity in cranial neural crest cells and craniofacial tissues, several regions harbor multiple signals with associations to different facial phenotypes, and there is evidence for potential coordinated actions of variants. In summary, our analyses provide insights into the understanding of how complex morphological traits are shaped by both individual and coordinated genetic actions.

Germline specification in mammals occurs through an inductive process whereby competent cells in the post-implantation epiblast differentiate into primordial germ cells (PGC). The intrinsic factors that endow epiblast cells with the competence to respond to germline inductive signals remain unknown. Single-cell RNA sequencing across multiple stages of an in vitro PGC-like cells (PGCLC) differentiation system shows that PGCLC genes initially expressed in the naïve pluripotent stage become homogeneously dismantled in germline competent epiblast like-cells (EpiLC). In contrast, the decommissioning of enhancers associated with these germline genes is incomplete. Namely, a subset of these enhancers partly retain H3K4me1, accumulate less heterochromatic marks and remain accessible and responsive to transcriptional activators. Subsequently, as in vitro germline competence is lost, these enhancers get further decommissioned and lose their responsiveness to transcriptional activators. Importantly, using H3K4me1-deficient cells, we show that the loss of this histone modification reduces the germline competence of EpiLC and decreases PGCLC differentiation efficiency. Our work suggests that, although H3K4me1 might not be essential for enhancer function, it can facilitate the (re)activation of enhancers and the establishment of gene expression programs during specific developmental transitions.