Joint MiRNA/mRNA expression profiling reveals changes consistent with development of dysfunctional corpus luteum after weight gain

High-Fat Diet and Female Fertility across Lifespan: A Comparative Lesson from Mammal Models

Female reproduction focuses mainly on achieving fully grown follicles and competent oocytes to be successfully fertilized, as well as on nourishing the developing offspring once pregnancy occurs. Current evidence demonstrates that obesity and/or high-fat diet regimes can perturbate these processes, leading to female infertility and transgenerational disorders. Since the mechanisms and reproductive processes involved are not yet fully clarified, the present review is designed as a systematic and comparative survey of the available literature. The available data demonstrate the adverse influences of obesity on diverse reproductive processes, such as folliculogenesis, oogenesis, and embryo development/implant. The negative reproductive impact may be attributed to a direct action on reproductive somatic and germinal compartments and/or to an indirect influence mediated by the endocrine, metabolic, and immune axis control systems. Overall, the present review highlights the fragmentation of the current information limiting the comprehension of the reproductive impact of a high-fat diet. Based on the incidence and prevalence of obesity in the Western countries, this topic becomes a research challenge to increase self-awareness of dietary reproductive risk to propose solid and rigorous preventive dietary regimes, as well as to develop targeted pharmacological interventions.

A Systematic Review of the Effects of High-Fat Diet Exposure on Oocyte and Follicular Quality: A Molecular Point of View

Worldwide, infertility affects between 10 and 15% of reproductive-aged couples. Female infertility represents an increasing health issue, principally in developing countries, as the current inclinations of delaying pregnancy beyond 35 years of age significantly decrease fertility rates. Female infertility, commonly imputable to ovulation disorders, can be influenced by several factors, including congenital malformations, hormonal dysfunction, and individual lifestyle choices, such as smoking cigarettes, stress, drug use and physical activity. Moreover, diet-related elements play an important role in the regulation of ovulation. Modern types of diet that encourage a high fat intake exert a particularly negative effect on ovulation, affecting the safety of gametes and the implantation of a healthy embryo. Identifying and understanding the cellular and molecular mechanisms responsible for diet-associated infertility might help clarify the confounding multifaceted elements of infertility and uncover novel, potentially curative treatments. In this view, this systematic revision of literature will summarize the current body of knowledge of the potential effect of high-fat diet (HFD) exposure on oocyte and follicular quality and consequent female reproductive function, with particular reference to molecular mechanisms and pathways. Inflammation, oxidative stress, gene expression and epigenetics represent the main mechanisms associated with mammal folliculogenesis and oogenesis.

Follicle-stimulating hormone receptors expression in ovine corpora lutea during luteal phase: effect of nutritional plane and follicle-stimulating hormone treatment

Corpus luteum (CL), a transient endocrine gland critical for reproductive cyclicity and pregnancy maintenance, is controlled by numerous regulatory factors. Although LH is widely recognized as the major regulator, other factors may also affect luteal functions. It has been demonstrated that FSH receptors (FSHR) are expressed not only in ovarian follicles but also in other tissues within the reproductive tract, including the CL. To evaluate FSHR expression in nontreated (nonsuperovulated; experiment 1) or FSH-treated (superovulated; experiment 2) sheep fed a control (C; maintenance), excess (O; 2 × C), or restricted (U; 0.6 × C) diet, CL were collected at the early, mid and/or late luteal phases (n = 5-7 per group). Protein and messenger RNA (mRNA) expression of FSHR were detected in the CL from all groups using immunohistochemistry followed by image analysis and quantitative RT-PCR, respectively. Follicle-stimulating hormone receptor was immunolocalized to steroidogenic small and large and nonsteroidogenic luteal cells. In both experiments, FSHR protein expression was not affected by stage of luteal development or diet. In experiment 1, expression of mRNA for all FSHR variants was greater (P <0.02 to 0.0003) at the late phase than mid or early luteal phase, and in experiment 2, it was greater (P < 0.001) at the mid than early luteal phase. Plane of nutrition did not affect FSHR mRNA expression. Comparison of FSH-treated with nontreated ewes demonstrated that FSH increased FSHR protein expression by 1.5- to 2-fold (P < 0.0001) in all groups, and mRNA expression by 7- to 30-fold (P < 0.001) for (1) FSHR-1 in all groups except U at the early luteal phase, (2) FSHR-2 in C, O, and U at the mid-phase, but not early luteal phase, and (3) FSHR-3 in U at the mid-luteal phase. Our data demonstrate that (1) FSHRs are expressed in ovine CL at several stages of luteal development, (2) FSHR protein expression does not change during the luteal phase and is not affected by diet, (3) FSHR mRNA expression not only depends on the stage of the estrous cycle but also not affected by diet in nonsuperovulated or superovulated ewes, and (4) in vivo FSH treatment enhanced FSHR protein and/or mRNA expression in the CL depending on diet and phase of the estrous cycle. Presence of FSHR in the CL indicates a regulatory role of FSH in luteal function in sheep. As very little is known about the possible role of FSH and FSHR in luteal functions, further studies should be undertaken to elucidate the endocrine, molecular, and cellular mechanisms of FSH effects on the CL.

Establishment of patient-derived orthotopic xenograft model of 1q+ posterior fossa group A ependymoma

BACKGROUND

Treatment for pediatric posterior fossa group A (PFA) ependymoma with gain of chromosome 1q (1q+) has not improved over the past decade owing partially to lack of clinically relevant models. We described the first 2 1q+ PFA cell lines, which have significantly enhanced our understanding of PFA tumor biology and provided a tool to identify specific 1q+ PFA therapies. However, cell lines do not accurately replicate the tumor microenvironment. Our present goal is to establish patient-derived xenograft (PDX) mouse models.

METHODS

Disaggregated tumors from 2 1q+ PFA patients were injected into the flanks of NSG mice. Flank tumors were then transplanted into the fourth ventricle or lateral ventricle of NSG mice. Characterization of intracranial tumors was performed using imaging, histology, and bioinformatics.

RESULTS

MAF-811_XC and MAF-928_XC established intracranially within the fourth ventricle and retained histological, methylomic, and transcriptomic features of primary patient tumors. We tested the feasibility of treating PDX mice with fractionated radiation or chemotherapy. Mice tolerated radiation despite significant tumor burden, and follow-up imaging confirmed radiation can reduce tumor size. Treatment with fluorouracil reduced tumor size but did not appear to prolong survival.

CONCLUSIONS

MAF-811_XC and MAF-928_XC are novel, authentic, and reliable models for studying 1q+ PFA in vivo. Given the successful response to radiation, these models will be advantageous for testing clinically relevant combination therapies to develop future clinical trials for this high-risk subgroup of pediatric ependymoma.

Transcriptional profiling of murine retinas undergoing semi-synchronous cone photoreceptor differentiation

Uncovering the gene regulatory networks that control cone photoreceptor formation has been hindered because cones only make up a few percent of the retina and form asynchronously during development. To overcome these limitations, we used a γ-secretase inhibitor, DAPT, to disrupt Notch signaling and force proliferating retinal progenitor cells to rapidly adopt neuronal identity. We treated mouse retinal explants at the peak of cone genesis with DAPT and examined tissues at several time-points by histology and bulk RNA-sequencing. We found that this treatment caused supernumerary cone formation in an overwhelmingly synchronized fashion. This analysis revealed several categorical patterns of gene expression changes over time relative to DMSO treated control explants. These were placed in the temporal context of the activation of Otx2, a transcription factor that is expressed at the onset of photoreceptor development and that is required for both rod and cone formation. One group of interest had genes, such as Mybl1, Ascl1, Neurog2, and Olig2, that became upregulated by DAPT treatment before Otx2. Two other groups showed upregulated gene expression shortly after Otx2, either transiently or permanently. This included genes such as Mybl1, Meis2, and Podxl. Our data provide a developmental timeline of the gene expression events that underlie the initial steps of cone genesis and maturation. Applying this strategy to human retinal organoid cultures was also sufficient to induce a massive increase in cone genesis. Taken together, our results provide a temporal framework that can be used to elucidate the gene regulatory logic controlling cone photoreceptor development.

High-fat diet exposure, regardless of induction of obesity, is associated with altered expression of genes critical to normal ovulatory function

We evaluated the impact of high-fat diet (HFD) on ovarian gene expression. Female 5-week-old C57BL/6J mice were fed a 60% HFD or standard chow for 10 weeks. HFD-fed mice were then separated into obese (HF-Ob) and lean (HF-Ln) based on body weight. HFD exposure led to impairment of the estrous cycle, changes in hormones affecting reproduction, and decreased primordial follicles regardless of the development of obesity. RNA-sequencing of whole ovaries identified multiple genes with altered expression after HFD, with 25 genes displaying decreased expression in both HF-Ln and HF-Ob mice compared to the chow-fed controls (q < 0.05). Several of these 25 genes are involved in normal ovarian functions, including ovulation (Edn2, Tnfaip6, Errfi1, Prkg2, and Nfil3), luteinization (Edn2), and luteolysis (Nr4a1). Taken together, elevated dietary fat intake, regardless of obesity, is associated with impaired estrous cycle, depletion of the ovarian reserve, and altered expression of genes critical to normal ovulatory function.

Elucidating the Role of the Maternal Embryonic Leucine Zipper Kinase in Adrenocortical Carcinoma

Adrenocortical carcinoma (ACC) is an aggressive cancer with a 5-year survival rate <35%. Mortality remains high due to lack of targeted therapies. Using bioinformatic analyses, we identified maternal embryonic leucine zipper kinase (MELK) as 4.1-fold overexpressed in ACC compared with normal adrenal samples. High MELK expression in human tumors correlated with shorter survival and with increased expression of genes involved in cell division and growth. We investigated the functional effects of MELK inhibition using newly developed ACC cell lines with variable MELK expression, CU-ACC1 and CU-ACC2, compared with H295R cells. In vitro treatment with the MELK inhibitor, OTSSP167, resulted in a dose-dependent decrease in rates of cell proliferation, colony formation, and cell survival, with relative sensitivity of each ACC cell line based upon the level of MELK overexpression. To confirm a MELK-specific antitumorigenic effect, MELK was inhibited in H295R cells via multiple short hairpin RNAs. MELK silencing resulted in 1.9-fold decrease in proliferation, and 3- to 10-fold decrease in colony formation in soft agar and clonogenicity assays, respectively. In addition, although MELK silencing had no effect on survival in normoxia, exposure to a hypoxia resulted in a sixfold and eightfold increase in apoptosis as assessed by caspase-3 activation and TUNEL, respectively. Together these data suggest that MELK is a modulator of tumor cell growth and survival in a hypoxic microenvironment in adrenal cancer cells and support future investigation of its role as a therapeutic kinase target in patients with ACC.

Increased FGF8 signaling promotes chondrogenic rather than osteogenic development in the embryonic skull

The bones of the cranial vault are formed directly from mesenchymal cells through intramembranous ossification rather than via a cartilage intermediate. Formation and growth of the skull bones involves the interaction of multiple cell-cell signaling pathways, with fibroblast growth factors (FGFs) and their receptors exerting a prominent influence. Mutations within the FGF signaling pathway are the most frequent cause of craniosynostosis, which is a common human craniofacial developmental abnormality characterized by the premature fusion of the cranial sutures. Here, we have developed new mouse models to investigate how different levels of increased FGF signaling can affect the formation of the calvarial bones and associated sutures. Whereas moderate Fgf8 overexpression resulted in delayed ossification followed by craniosynostosis of the coronal suture, higher Fgf8 levels promoted a loss of ossification and favored cartilage over bone formation across the skull. By contrast, endochondral bones were still able to form and ossify in the presence of increased levels of Fgf8, although the growth and mineralization of these bones were affected to varying extents. Expression analysis demonstrated that abnormal skull chondrogenesis was accompanied by changes in the genes required for Wnt signaling. Moreover, further analysis indicated that the pathology was associated with decreased Wnt signaling, as the reduction in ossification could be partially rescued by halving Axin2 gene dosage. Taken together, these findings indicate that mesenchymal cells of the skull are not fated to form bone, but can be forced into a chondrogenic fate through the manipulation of FGF8 signaling. These results have implications for evolution of the different methods of ossification as well as for therapeutic intervention in craniosynostosis.

Development of new preclinical models to advance adrenocortical carcinoma research

Adrenocortical cancer (ACC) is an orphan malignancy that results in heterogeneous clinical phenotypes and molecular genotypes. There are no curative treatments for this deadly cancer with 35% survival at five years. Our understanding of the underlying pathobiology and our ability to test novel therapeutic targets has been limited due to the lack of preclinical models. Here, we report the establishment of two new ACC cell lines and corresponding patient-derived xenograft (PDX) models. CU-ACC1 cell line and PDX were derived from a perinephric metastasis in a patient whose primary tumor secreted aldosterone. CU-ACC2 cell line and PDX were derived from a liver metastasis in a patient with Lynch syndrome. Short tandem repeat profiling confirmed consistent matches between human samples and models. Both exomic and RNA sequencing profiling were performed on the patient samples and the models, and hormonal secretion was evaluated in the new cell lines. RNA sequencing and immunohistochemistry confirmed the expression of adrenal cortex markers in the PDXs and human tumors. The new cell lines replicate two of the known genetic models of ACC. CU-ACC1 cells had a mutation in CTNNB1 and secreted cortisol but not aldosterone. CU-ACC2 cells had a TP53 mutation and loss of MSH2 consistent with the patient's known germline mutation causing Lynch syndrome. Both cell lines can be transfected and transduced with similar growth rates. These new preclinical models of ACC significantly advance the field by allowing investigation of underlying molecular mechanisms of ACC and the ability to test patient-specific therapeutic targets.

AP-2α and AP-2β cooperatively orchestrate homeobox gene expression during branchial arch patterning

The evolution of a hinged moveable jaw with variable morphology is considered a major factor behind the successful expansion of the vertebrates. DLX homeobox transcription factors are crucial for establishing the positional code that patterns the mandible, maxilla and intervening hinge domain, but how the genes encoding these proteins are regulated remains unclear. Herein, we demonstrate that the concerted action of the AP-2α and AP-2β transcription factors within the mouse neural crest is essential for jaw patterning. In the absence of these two proteins, the hinge domain is lost and there are alterations in the size and patterning of the jaws correlating with dysregulation of homeobox gene expression, with reduced levels of Emx, Msx and Dlx paralogs accompanied by an expansion of Six1 expression. Moreover, detailed analysis of morphological features and gene expression changes indicate significant overlap with various compound Dlx gene mutants. Together, these findings reveal that the AP-2 genes have a major function in mammalian neural crest development, influencing patterning of the craniofacial skeleton via the DLX code, an effect that has implications for vertebrate facial evolution, as well as for human craniofacial disorders.

RNA sequencing and pathway analysis identify tumor necrosis factor alpha driven small proline-rich protein dysregulation in chronic rhinosinusitis

BACKGROUND

Chronic rhinosinusitis (CRS) is a heterogeneous inflammatory disorder in which many pathways contribute to end-organ disease. Small proline-rich proteins (SPRR) are polypeptides that have recently been shown to contribute to epithelial biomechanical properties relevant in T-helper type 2 inflammation. There is evidence that genetic polymorphism in SPRR genes may predict the development of asthma in children with atopy and, correlatively, that expression of SPRRs is increased under allergic conditions, which leads to epithelial barrier dysfunction in atopic disease.

METHODS

RNAs from uncinate tissue specimens from patients with CRS and control subjects were compared by RNA sequencing by using Ingenuity Pathway Analysis (n = 4 each), and quantitative polymerase chain reaction (PCR) (n = 15). A separate cohort of archived sinus tissue was examined by immunohistochemistry (n = 19).

RESULTS

A statistically significant increase of SPRR expression in CRS sinus tissue was identified that was not a result of atopic presence. SPRR1 and SPRR2A expressions were markedly increased in patients with CRS (p < 0.01) on RNA sequencing, with confirmation by using real-time PCR. Immunohistochemistry of archived surgical samples demonstrated staining of SPRR proteins within squamous epithelium of both groups. Pathway analysis indicated tumor necrosis factor (TNF) alpha as a master regulator of the SPRR gene products.

CONCLUSION

Expression of SPRR1 and of SPRR2A is increased in mucosal samples from patients with CRS and appeared as a downstream result of TNF alpha modulation, which possibly resulted in epithelial barrier dysfunction.

Chi-miR-4110 promotes granulosa cell apoptosis by targeting Sma- and Mad-related protein 2 (Smad2) in the caprine ovary

Follicular atresia mainly results from the apoptosis of granulosa cells (GCs). Whilst our previous investigations examined the role of chi-miR-4110 in regulating ovarian function, the present study detected the role of chi-miR-4110 in GC development. We transfected caprine GCs cultured in vitro with chi-miR-4110 mimics. Results revealed that chi-miR-4110 decreased mRNA and protein levels of Smad2 by targeting its 3'-untranslated region (3'UTR). FoxC1 and Sp1 mRNA and protein levels markedly increased, whereas those of bHLHe22 significantly decreased (P<0.01 or 0.05) in GCs transfected with the chi-miR-4110 mimics. Further studies revealed a significantly higher number of apoptotic cells in GCs transfected with the chi-miR-4110 mimics (P< 0.05) than in GCs transfected with mimics negative control. GCs transfected with the chi-miR-4110 mimics exhibited significantly increased mRNA and protein levels of the pro-apoptotic gene Bax (P<0.01) and significantly decreased expression levels of the anti-apoptotic gene BCL-2 (P<0.01). Smad2 interference (Si-1282) results were consistent with those of the chi-miR-4110 mimics. Previous reports and our results showed that chi-miR-4110 increases Sp1 expression by repressing Smad2. The increase in Sp1 induces p53-upregulated modulator of apoptosis, which increases the relative abundance of Bax and causes caprine GC apoptosis. Our findings may provide relevant data for the investigation of miRNA-mediated regulation of ovarian functions.

The RNA-binding protein caper is required for sensory neuron development in Drosophila melanogaster

BACKGROUND

Alternative splicing mediated by RNA-binding proteins (RBPs) is emerging as a fundamental mechanism for the regulation of gene expression. Alternative splicing has been shown to be a widespread phenomenon that facilitates the diversification of gene products in a tissue-specific manner. Although defects in alternative splicing are rooted in many neurological disorders, only a small fraction of splicing factors have been investigated in detail.

RESULTS

We find that the splicing factor Caper is required for the development of multiple different mechanosensory neuron subtypes at multiple life stages in Drosophila melanogaster. Disruption of Caper function causes defects in dendrite morphogenesis of larval dendrite arborization neurons and neuronal positioning of embryonic proprioceptors, as well as the development and maintenance of adult mechanosensory bristles. Additionally, we find that Caper dysfunction results in aberrant locomotor behavior in adult flies. Transcriptome-wide analyses further support a role for Caper in alternative isoform regulation of genes that function in neurogenesis.

CONCLUSIONS

Our results provide the first evidence for a fundamental and broad requirement for the highly conserved splicing factor Caper in the development and maintenance of the nervous system and provide a framework for future studies on the detailed mechanism of Caper-mediated RNA regulation. Developmental Dynamics 246:610-624, 2017. © 2017 Wiley Periodicals, Inc.

Gsg1, Trnp1, and Tmem215 Mark Subpopulations of Bipolar Interneurons in the Mouse Retina

Purpose

How retinal bipolar cell interneurons are specified and assigned to specialized subtypes is only partially understood. In part, this is due to a lack of early pan- and subtype-specific bipolar cell markers. To discover these factors, we identified genes that were upregulated in Blimp1 (Prdm1) mutant retinas, which exhibit precocious bipolar cell development.

Methods

Postnatal day (P)2 retinas from Blimp1 conditional knock-out (CKO) mice and controls were processed for RNA sequencing. Genes that increased at least 45% and were statistically different between conditions were considered candidate bipolar-specific factors. Candidates were further evaluated by RT-PCR, in situ hybridization, and immunohistochemistry. Knock-in Tmem215-LacZ mice were used to better trace retinal expression.

Results

A comparison between Blimp1 CKO and control RNA-seq datasets revealed approximately 40 significantly upregulated genes. We characterized the expression of three genes that have no known function in the retina, Gsg1 (germ cell associated gene), Trnp1 (TMF-regulated nuclear protein), and Tmem215 (a predicted transmembrane protein). Germ cell associated gene appeared restricted to a small subset of cone bipolars while Trnp1 was seen in all ON type bipolar cells. Using Tmem215-LacZ heterozygous knock-in mice, we observed that β-galactosidase expression started early in bipolar cell development. In adults, Tmem215 was expressed by a subset of ON and OFF cone bipolar cells.

Conclusions

We have identified Gsg1, Tmem215, and Trnp1 as novel bipolar subtype-specific genes. The spatial and temporal pattern of their expression is consistent with a role in controlling bipolar subtype fate choice, differentiation, or physiology.

An autonomous metabolic role for Spen

Preventing obesity requires a precise balance between deposition into and mobilization from fat stores, but regulatory mechanisms are incompletely understood. Drosophila Split ends (Spen) is the founding member of a conserved family of RNA-binding proteins involved in transcriptional regulation and frequently mutated in human cancers. We find that manipulating Spen expression alters larval fat levels in a cell-autonomous manner. Spen-depleted larvae had defects in energy liberation from stores, including starvation sensitivity and major changes in the levels of metabolic enzymes and metabolites, particularly those involved in β-oxidation. Spenito, a small Spen family member, counteracted Spen function in fat regulation. Finally, mouse Spen and Spenito transcript levels scaled directly with body fat in vivo, suggesting a conserved role in fat liberation and catabolism. This study demonstrates that Spen is a key regulator of energy balance and provides a molecular context to understand the metabolic defects that arise from Spen dysfunction.

All animals need energy to fuel development and survive as adults. Excess energy stored as fat provides a means to endure periods when external energy is unavailable, but there is a delicate balance between accumulating sufficient fat stores and becoming obese. While the enzymes that mediate energy deposition into and mobilization from fat stores are well studied, the complex upstream regulatory pathways have not been fully worked out. We report here that two members of a conserved family of RNA-binding proteins, Spen and Nito, operate in fat storage cells in fruit fly larvae to control the expression of genes that mediate energy liberation from fat stores. Manipulating Spen or Spenito function grossly perturbs larval energy metabolism, including imbalances in the amounts of stored fats, key metabolites, and metabolic enzymes, and resulting in defects in survival under starvation conditions. Interestingly, Nito opposes Spen functions, indicative of a regulatory mechanism that helps keep energy balance in check. We find that the mouse homologs of Spen and Nito, which were known to regulate gene expression in other pathways, respond similarly to changes in body fat induced by a high-fat diet, suggesting that the balancing effect of these two proteins also prevents mammalian obesity.

Luteal function during the estrous cycle in arginine-treated ewes fed different planes of nutrition

Functions of corpus luteum (CL) are influenced by numerous factors including hormones, growth and angiogenic factors, nutritional plane and dietary supplements such as arginine (Arg), a semi-essential amino acid and precursor for proteins, polyamines and nitric oxide (NO). The aim of this study was to determine if Arg supplementation to ewes fed different planes of nutrition influences: (1) progesterone (P4) concentrations in serum and luteal tissue, (2) luteal vascularity, cell proliferation, endothelial NO synthase (eNOS) and receptor (R) soluble guanylate cyclase β protein and mRNA expression and (3) luteal mRNA expression for selected angiogenic factors during the estrous cycle. Ewes (n = 111) were categorized by weight and randomly assigned to one of three nutritional planes: maintenance control (C), overfed (2× C) and underfed (0.6× C) beginning 60 days prior to onset of estrus. After estrus synchronization, ewes from each nutritional plane were assigned randomly to one of two treatments: Arg or saline. Serum and CL were collected at the early, mid and late luteal phases. The results demonstrated that: (1) nutritional plane affected ovulation rates, luteal vascularity, cell proliferation andNOS3,GUCY1B3, vascular endothelial growth factor (VEGF) andVEGFR2mRNA expression, (2) Arg affected luteal vascularity, cell proliferation andNOS3,GUCY1B3,VEGFandVEGFR2mRNA expression and (3) luteal vascularity, cell proliferation and the VEGF and NO systems depend on the stage of the estrous cycle. These data indicate that plane of nutrition and/or Arg supplementation can alter vascularization and expression of selected angiogenic factors in luteal tissue during the estrous cycle in sheep.

Genomewide Association Study of African Children Identifies Association of SCHIP1 and PDE8A with Facial Size and Shape

The human face is a complex assemblage of highly variable yet clearly heritable anatomic structures that together make each of us unique, distinguishable, and recognizable. Relatively little is known about the genetic underpinnings of normal human facial variation. To address this, we carried out a large genomewide association study and two independent replication studies of Bantu African children and adolescents from Mwanza, Tanzania, a region that is both genetically and environmentally relatively homogeneous. We tested for genetic association of facial shape and size phenotypes derived from 3D imaging and automated landmarking of standard facial morphometric points. SNPs within genes SCHIP1 and PDE8A were associated with measures of facial size in both the GWAS and replication cohorts and passed a stringent genomewide significance threshold adjusted for multiple testing of 34 correlated traits. For both SCHIP1 and PDE8A, we demonstrated clear expression in the developing mouse face by both whole-mount in situ hybridization and RNA-seq, supporting their involvement in facial morphogenesis. Ten additional loci demonstrated suggestive association with various measures of facial shape. Our findings, which differ from those in previous studies of European-derived whites, augment understanding of the genetic basis of normal facial development, and provide insights relevant to both human disease and forensics.

The human face is made up of distinct yet related anatomic structures that together make both individuals and families recognizable. It is clear there is a strong genetic component to the human face, and though the genetics of the face have been studied for several years, there are relatively few genes known to impact normal human facial development and facial shape. We report here a large-scale human genetic study in which we successfully identify and replicate genetic markers associated with normal facial variation using advanced 3D facial imaging in African children. We identified two significant replicated genes associated with measures of human facial size, SCHIP1 and PDE8A, demonstrated their clear expression in the developing face in the mouse, and identified 10 additional candidate genetic loci for human facial shape. Gene discovery for human facial development is an important first step for both diagnosing and treating craniofacial syndromes and for developing forensic modeling of the human face.