Rare variants analyses suggest novel cleft genes in the African population.

Non-syndromic orofacial clefts (NSOFCs) are common birth defects with a complex etiology. While over 60 common risk loci have been identified, they explain only a small proportion of the heritability for NSOFCs. Rare variants have been implicated in the missing heritability. Thus, our study aimed to identify genes enriched with nonsynonymous rare coding variants associated with NSOFCs. Our sample included 814 non-syndromic cleft lip with or without palate (NSCL/P), 205 non-syndromic cleft palate only (NSCPO), and 2150 unrelated control children from Nigeria, Ghana, and Ethiopia. We conducted a gene-based analysis separately for each phenotype using three rare-variants collapsing models: (1) protein-altering (PA), (2) missense variants only (MO); and (3) loss of function variants only (LOFO). Subsequently, we utilized relevant transcriptomics data to evaluate associated gene expression and examined their mutation constraint using the gnomeAD database. In total, 13 genes showed suggestive associations (p = E-04). Among them, eight genes (ABCB1, ALKBH8, CENPF, CSAD, EXPH5, PDZD8, SLC16A9, and TTC28) were consistently expressed in relevant mouse and human craniofacial tissues during the formation of the face, and three genes (ABCB1, TTC28, and PDZD8) showed statistically significant mutation constraint. These findings underscore the role of rare variants in identifying candidate genes for NSOFCs.

Effect of scaling on serum high-sensitivity C-reactive protein levels and periodontal parameters in systemically healthy women of reproductive age group with gingivitis.

OBJECTIVE: Varying levels of sex hormones across the menstrual cycle in young systemically healthy females may alter tissue responses to plaque, resulting in increased gingival inflammation. Also, higher severity and prevalence of gingivitis has been demonstrated in adult females than males, attributed to hormonal changes. Further, literature reported that Gingivitis raises the levels of systemic inflammatory markers such as C Reactive Protein. This interventional trial aimed to evaluate the effect of supragingival scaling on serum highsensitivity C-reactive protein (hsCRP) levels along with periodontal parameters in systemically healthy women of reproductive age group with natural gingivitis. MATERIAL AND METHODS: 57 women of reproductive age were enrolled into two groups. Test Group (n=30) comprised of systemically healthy women with gingivitis who received supragingival scaling. Control Group (n=27) included systemically and periodontally healthy females. Periodontal parameters [gingival index (GI), plaque index (PI), pocket probing depth (PPD), bleeding on probing (BOP)], and serum hsCRP levels were recorded at baseline for both the groups. Follow up of Test Group participants was done at 3 and 6 months. RESULTS: Serum hsCRP and periodontal parameters were significantly higher in Test Group than control group at baseline which decreased significantly after treatment in Test Group at 6 months follow up (p≤0.05). GI, BOP and hsCRP in Test Group at 6 months were reduced up to the baseline levels of systemically and periodontally healthy females. CONCLUSION: Treatment of gingival inflammation can help in lowering the systemic and local inflammation up to the levels of systemically and periodontally healthy females.

Shared genetic risk between major orofacial cleft phenotypes in an African population.

Nonsyndromic orofacial clefts (NSOFCs) represent a large proportion (70%-80%) of all OFCs. They can be broadly categorized into nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO). Although NSCL/P and NSCPO are considered etiologically distinct, recent evidence suggests the presence of shared genetic risks. Thus, we investigated the genetic overlap between NSCL/P and NSCPO using African genome-wide association study (GWAS) data on NSOFCs. These data consist of 814 NSCL/P, 205 NSCPO cases, and 2159 unrelated controls. We generated common single-nucleotide variants (SNVs) association summary statistics separately for each phenotype (NSCL/P and NSCPO) under an additive genetic model. Subsequently, we employed the pleiotropic analysis under the composite null (PLACO) method to test for genetic overlap. Our analysis identified two loci with genome-wide significance (rs181737795 [p = 2.58E-08] and rs2221169 [p = 4.5E-08]) and one locus with marginal significance (rs187523265 [p = 5.22E-08]). Using mouse transcriptomics data and information from genetic phenotype databases, we identified MDN1, MAP3k7, KMT2A, ARCN1, and VADC2 as top candidate genes for the associated SNVs. These findings enhance our understanding of genetic variants associated with NSOFCs and identify potential candidate genes for further exploration.

Rare Variants Analyses Suggest Novel Cleft Genes in the African Population.

Non-syndromic orofacial clefts (NSOFCs) are common birth defects with a complex etiology. While over 60 common risk loci have been identified, they explain only a small proportion of the heritability for NSOFC. Rare variants have been implicated in the missing heritability. Thus, our study aimed to identify genes enriched with nonsynonymous rare coding variants associated with NSOFCs. Our sample included 814 non-syndromic cleft lip with or without palate (NSCL/P), 205 non-syndromic cleft palate only (NSCPO), and 2150 unrelated control children from Nigeria, Ghana, and Ethiopia. We conducted a gene-based analysis separately for each phenotype using three rare-variants collapsing models: (1) protein-altering (PA), (2) missense variants only (MO); and (3) loss of function variants only (LOFO). Subsequently, we utilized relevant transcriptomics data to evaluate associated gene expression and examined their mutation constraint using the gnomeAD database. In total, 13 genes showed suggestive associations (p = E-04). Among them, eight genes (ABCB1, ALKBH8, CENPF, CSAD, EXPH5, PDZD8, SLC16A9, and TTC28) were consistently expressed in relevant mouse and human craniofacial tissues during the formation of the face, and three genes (ABCB1, TTC28, and PDZD8) showed statistically significant mutation constraint. These findings underscore the role of rare variants in identifying candidate genes for NSOFCs.

Misconduct in Biomedical Research: A Meta-Analysis and Systematic Review.

Aims and Objectives: This study aimed to assess the nature and prevalence of misconduct in self and nonself-reported biomedical research. Materials and Methods: A detailed review of previously conducted studies was conducted through PubMed Central, PubMed, and Google Scholar using MeSH terms: "scientific misconduct," "Publications," "plagiarism," and "authorship," and keywords: scientific misconduct, gift authorship, ghost authorship, and duplicate publication. MeSH terms and keywords were searched in combinations using Boolean operators "AND" and "OR." Of 7771 articles that appeared in the search, 107 were selected for inspection. The articles were screened for their quality and inclusion criteria. Finally, 16 articles were selected for meta-analysis. Data analysis was conducted using an Open-Source, Open Meta Analyst, statistical software using the package "metaphor." Results: Plagiarism, data fabrication, and falsification were prevalent in most articles reviewed. The prevalence of research misconduct for plagiarism was 4.2% for self-reported and 27.9% for nonself-reported studies. Data fabrication was 4.5% in self-reported and 21.7% in nonself-reported studies. Data falsification was 9.7% in self-reported and 33.4% in nonself-reported studies, with significant heterogeneity. Conclusion: This meta-analysis gives a pooled estimate of the misconduct in research done in biomedical fields such as medicine, dental, pharmacy, and others across the world. We found that there is an alarming rate of misconduct in recent nonself-reported studies, and they were higher than that in the self-reported studies.

Proteomic profiling of retina and retinal pigment epithelium combined embryonic tissue to facilitate ocular disease gene discovery.

To expedite gene discovery in eye development and its associated defects, we previously developed a bioinformatics resource-tool iSyTE (integrated Systems Tool for Eye gene discovery). However, iSyTE is presently limited to lens tissue and is predominantly based on transcriptomics datasets. Therefore, to extend iSyTE to other eye tissues on the proteome level, we performed high-throughput tandem mass spectrometry (MS/MS) on mouse embryonic day (E)14.5 retina and retinal pigment epithelium combined tissue and identified an average of 3300 proteins per sample (n = 5). High-throughput expression profiling-based gene discovery approaches-involving either transcriptomics or proteomics-pose a key challenge of prioritizing candidates from thousands of RNA/proteins expressed. To address this, we used MS/MS proteome data from mouse whole embryonic body (WB) as a reference dataset and performed comparative analysis-termed "in silico WB-subtraction"-with the retina proteome dataset. In silico WB-subtraction identified 90 high-priority proteins with retina-enriched expression at stringency criteria of ≥ 2.5 average spectral counts, ≥ 2.0 fold-enrichment, false discovery rate < 0.01. These top candidates represent a pool of retina-enriched proteins, several of which are associated with retinal biology and/or defects (e.g., Aldh1a1, Ank2, Ank3, Dcn, Dync2h1, Egfr, Ephb2, Fbln5, Fbn2, Hras, Igf2bp1, Msi1, Rbp1, Rlbp1, Tenm3, Yap1, etc.), indicating the effectiveness of this approach. Importantly, in silico WB-subtraction also identified several new high-priority candidates with potential regulatory function in retina development. Finally, proteins exhibiting expression or enriched-expression in the retina are made accessible in a user-friendly manner at iSyTE ( https://research.bioinformatics.udel.edu/iSyTE/ ), to allow effective visualization of this information and facilitate eye gene discovery.

High-Throughput Transcriptomics of Celf1 Conditional Knockout Lens Identifies Downstream Networks Linked to Cataract Pathology.

Defects in the development of the ocular lens can cause congenital cataracts. To understand the various etiologies of congenital cataracts, it is important to characterize the genes linked to this developmental defect and to define their downstream pathways that are relevant to lens biology and pathology. Deficiency or alteration of several RNA-binding proteins, including the conserved RBP Celf1 (CUGBP Elav-like family member 1), has been described to cause lens defects and early onset cataracts in animal models and/or humans. Celf1 is involved in various aspects of post-transcriptional gene expression control, including regulation of mRNA stability/decay, alternative splicing and translation. Celf1 germline knockout mice and lens conditional knockout (Celf1cKO) mice develop fully penetrant cataracts in early postnatal stages. To define the genome-level changes in RNA transcripts that result from Celf1 deficiency, we performed high-throughput RNA-sequencing of Celf1cKO mouse lenses at postnatal day (P) 0. Celf1cKO lenses exhibit 987 differentially expressed genes (DEGs) at cut-offs of >1.0 log2 counts per million (CPM), ≥±0.58 log2 fold-change and <0.05 false discovery rate (FDR). Of these, 327 RNAs were reduced while 660 were elevated in Celf1cKO lenses. The DEGs were subjected to various downstream analyses including iSyTE lens enriched-expression, presence in Cat-map, and gene ontology (GO) and representation of regulatory pathways. Further, a comparative analysis was done with previously generated microarray datasets on Celf1cKO lenses P0 and P6. Together, these analyses validated and prioritized several key genes mis-expressed in Celf1cKO lenses that are relevant to lens biology, including known cataract-linked genes (e.g., Cryab, Cryba2, Cryba4, Crybb1, Crybb2, Cryga, Crygb, Crygc, Crygd, Cryge, Crygf, Dnase2b, Bfsp1, Gja3, Pxdn, Sparc, Tdrd7, etc.) as well as novel candidates (e.g., Ell2 and Prdm16). Together, these data have defined the alterations in lens transcriptome caused by Celf1 deficiency, in turn uncovering downstream genes and pathways (e.g., structural constituents of eye lenses, lens fiber cell differentiation, etc.) associated with lens development and early-onset cataracts.

Proteomic profiling of retina and retinal pigment epithelium combined embryonic tissue to facilitate ocular disease gene discovery.

To expedite gene discovery in eye development and its associated defects, we previously developed a bioinformatics resource-tool iSyTE (integrated Systems Tool for Eye gene discovery). However, iSyTE is presently limited to lens tissue and is predominantly based on transcriptomics datasets. Therefore, to extend iSyTE to other eye tissues on the proteome level, we performed high-throughput tandem mass spectrometry (MS/MS) on mouse embryonic day (E)14.5 retina and retinal pigment epithelium combined tissue and identified an average of 3,300 proteins per sample (n=5). High-throughput expression profiling-based gene discovery approaches-involving either transcriptomics or proteomics-pose a key challenge of prioritizing candidates from thousands of RNA/proteins expressed. To address this, we used MS/MS proteome data from mouse whole embryonic body (WB) as a reference dataset and performed comparative analysis-termed "in silico WB-subtraction"-with the retina proteome dataset. In silico WB-subtraction identified 90 high-priority proteins with retina-enriched expression at stringency criteria of ³2.5 average spectral counts, ³2.0 fold-enrichment, False Discovery Rate <0.01. These top candidates represent a pool of retina-enriched proteins, several of which are associated with retinal biology and/or defects (e.g., Aldh1a1, Ank2, Ank3, Dcn, Dync2h1, Egfr, Ephb2, Fbln5, Fbn2, Hras, Igf2bp1, Msi1, Rbp1, Rlbp1, Tenm3, Yap1, etc.), indicating the effectiveness of this approach. Importantly, in silico WB-subtraction also identified several new high-priority candidates with potential regulatory function in retina development. Finally, proteins exhibiting expression or enriched-expression in the retina are made accessible in a user-friendly manner at iSyTE (https://research.bioinformatics.udel.edu/iSyTE/), to allow effective visualization of this information and facilitate eye gene discovery.

Lens Epithelial Explants Treated with Vitreous Humor Undergo Alterations in Chromatin Landscape with Concurrent Activation of Genes Associated with Fiber Cell Differentiation and Innate Immune Response.

Lens epithelial explants are comprised of lens epithelial cells cultured in vitro on their native basement membrane, the lens capsule. Biologists have used lens epithelial explants to study many different cellular processes including lens fiber cell differentiation. In these studies, fiber differentiation is typically measured by cellular elongation and the expression of a few proteins characteristically expressed by lens fiber cells in situ. Chromatin and RNA was collected from lens epithelial explants cultured in either un-supplemented media or media containing 50% bovine vitreous humor for one or five days. Chromatin for ATAC-sequencing and RNA for RNA-sequencing was prepared from explants to assess regions of accessible chromatin and to quantitatively measure gene expression, respectively. Vitreous humor increased chromatin accessibility in promoter regions of genes associated with fiber differentiation and, surprisingly, an immune response, and this was associated with increased transcript levels for these genes. In contrast, vitreous had little effect on the accessibility of the genes highly expressed in the lens epithelium despite dramatic reductions in their mRNA transcripts. An unbiased analysis of differentially accessible regions revealed an enrichment of cis-regulatory motifs for RUNX, SOX and TEAD transcription factors that may drive differential gene expression in response to vitreous.

Deficiency of the bZIP transcription factors Mafg and Mafk causes misexpression of genes in distinct pathways and results in lens embryonic developmental defects.

Deficiency of the small Maf proteins Mafg and Mafk cause multiple defects, namely, progressive neuronal degeneration, cataract, thrombocytopenia and mid-gestational/perinatal lethality. Previous data shows Mafg -/-:Mafk +/- compound knockout (KO) mice exhibit cataracts age 4-months onward. Strikingly, Mafg -/-:Mafk -/- double KO mice develop lens defects significantly early in life, during embryogenesis, but the pathobiology of these defects is unknown, and is addressed here. At embryonic day (E)16.5, the epithelium of lens in Mafg -/-:Mafk -/- animals appears abnormally multilayered as demonstrated by E-cadherin and nuclear staining. Additionally, Mafg -/-:Mafk -/- lenses exhibit abnormal distribution of F-actin near the "fulcrum" region where epithelial cells undergo apical constriction prior to elongation and reorientation as early differentiating fiber cells. To identify the underlying molecular changes, we performed high-throughput RNA-sequencing of E16.5 Mafg -/-:Mafk -/- lenses and identified a cohort of differentially expressed genes that were further prioritized using stringent filtering criteria and validated by RT-qPCR. Several key factors associated with the cytoskeleton, cell cycle or extracellular matrix (e.g., Cdk1, Cdkn1c, Camsap1, Col3a1, Map3k12, Sipa1l1) were mis-expressed in Mafg -/-:Mafk -/- lenses. Further, the congenital cataract-linked extracellular matrix peroxidase Pxdn was significantly overexpressed in Mafg -/-:Mafk -/- lenses, which may cause abnormal cell morphology. These data also identified the ephrin signaling receptor Epha5 to be reduced in Mafg -/-:Mafk -/- lenses. This likely contributes to the Mafg -/-:Mafk -/- multilayered lens epithelium pathology, as loss of an ephrin ligand, Efna5 (ephrin-A5), causes similar lens defects. Together, these findings uncover a novel early function of Mafg and Mafk in lens development and identify their new downstream regulatory relationships with key cellular factors.

Whole-genome sequencing reveals de-novo mutations associated with nonsyndromic cleft lip/palate.

The majority (85%) of nonsyndromic cleft lip with or without cleft palate (nsCL/P) cases occur sporadically, suggesting a role for de novo mutations (DNMs) in the etiology of nsCL/P. To identify high impact protein-altering DNMs that contribute to the risk of nsCL/P, we conducted whole-genome sequencing (WGS) analyses in 130 African case-parent trios (affected probands and unaffected parents). We identified 162 high confidence protein-altering DNMs some of which are based on available evidence, contribute to the risk of nsCL/P. These include novel protein-truncating DNMs in the ACTL6A, ARHGAP10, MINK1, TMEM5 and TTN genes; as well as missense variants in ACAN, DHRS3, DLX6, EPHB2, FKBP10, KMT2D, RECQL4, SEMA3C, SEMA4D, SHH, TP63, and TULP4. Many of these protein-altering DNMs were predicted to be pathogenic. Analysis using mouse transcriptomics data showed that some of these genes are expressed during the development of primary and secondary palate. Gene-set enrichment analysis of the protein-altering DNMs identified palatal development and neural crest migration among the few processes that were significantly enriched. These processes are directly involved in the etiopathogenesis of clefting. The analysis of the coding sequence in the WGS data provides more evidence of the opportunity for novel findings in the African genome.

Oral mesotherapy technique for the treatment of physiologic gingival melanin hyperpigmentation using locally injectable vitamin C: a clinical and histologic cases series.

OBJECTIVES: Vitamin C (ascorbic acid) is widely used in dermatology for skin depigmentation. However, there are very few clinical studies on the efficacy of vitamin C in gingival depigmentation. This preliminary case series aims to present the clinical effectiveness, histologic changes, and patient-reported outcomes of intra-epidermal vitamin C injections (oral mesotherapy) for managing patients with gingival melanin hyperpigmentation. METHOD AND MATERIALS: Five patients were administered locally injectable vitamin C (once per week for 4 to 5 visits) in maxillary or mandibular anterior pigmented gingiva. The depigmentation effect was evaluated using the Dummett Oral Pigmentation Index (DOPI) and Gingival Pigmentation Index (GPI). Digital photographs were used to assess gingival luminescence (L*) and pigmented surface area (PSA). Parameters were recorded at baseline and at 1 and 3 months. Melanocyte histopathologic count was determined at baseline and at 3 months. Pain, gingival color change, and patient satisfaction scores were also assessed. RESULTS: Median GPI, DOPI, and PSA were significantly reduced (P ≤ .05) from baseline to 1 month. There was no statistically significant change from 1 month to 3 months. L* value significantly increased from baseline to 3 months. A median pain score of 3 (scale of 0 to 10) was observed on the day of the procedure. A score of 3 (scale of 0 to 4) was reported for the gingival color and scores 3 and 4 (scale of 0 to 4) for the overall patient satisfaction. CONCLUSION: Oral mesotherapy using locally injectable vitamin C is a nonsurgical, minimally invasive, and efficient technique for gingival depigmentation. Indian patients were satisfied with the gingival color obtained and the overall treatment experience. CLINICAL IMPORTANCE: As all the branches of medicine, specifically dentistry, direct to minimally invasive approaches, mesotherapy shows great importance to reduce the surgical interventions, especially when esthetic outcomes are needed. Oral mesotherapy using locally injectable vitamin C can be a useful nonsurgical technique for achieving gingival depigmentation in the esthetic zone.

Variant analyses of candidate genes in orofacial clefts in multi-ethnic populations.

OBJECTIVES: Cleft lip with/without cleft palate and cleft palate only is congenital birth defects where the upper lip and/or palate fail to fuse properly during embryonic facial development. Affecting ~1.2/1000 live births worldwide, these orofacial clefts impose significant social and financial burdens on affected individuals and their families. Orofacial clefts have a complex etiology resulting from genetic variants combined with environmental covariates. Recent genome-wide association studies and whole-exome sequencing for orofacial clefts identified significant genetic associations and variants in several genes. Of these, we investigated the role of common/rare variants in SHH, RORA, MRPL53, ACVR1, and GDF11. MATERIALS AND METHODS: We sequenced these five genes in 1255 multi-ethnic cleft lip with/without palate and cleft palate only samples in order to find variants that may provide potential explanations for the missing heritability of orofacial clefts. Rare and novel variants were further analyzed using in silico predictive tools. RESULTS: Ninteen total variants of interest were found, with variant types including stop-gain, missense, synonymous, intronic, and splice-site variants. Of these, 3 novel missense variants were found, one in SHH, one in RORA, and one in GDF11. CONCLUSION: This study provides evidence that variants in SHH, RORA, MRPL53, ACVR1, and GDF11 may contribute to risk of orofacial clefts in various populations.

A large multiethnic GWAS meta-analysis of cataract identifies new risk loci and sex-specific effects.

Cataract is the leading cause of blindness among the elderly worldwide and cataract surgery is one of the most common operations performed in the United States. As the genetic etiology of cataract formation remains unclear, we conducted a multiethnic genome-wide association meta-analysis, combining results from the GERA and UK Biobank cohorts, and tested for replication in the 23andMe research cohort. We report 54 genome-wide significant loci, 37 of which were novel. Sex-stratified analyses identified CASP7 as an additional novel locus specific to women. We show that genes within or near 80% of the cataract-associated loci are significantly expressed and/or enriched-expressed in the mouse lens across various spatiotemporal stages as per iSyTE analysis. Furthermore, iSyTE shows 32 candidate genes in the associated loci have altered gene expression in 9 different gene perturbation mouse models of lens defects/cataract, suggesting their relevance to lens biology. Our work provides further insight into the complex genetic architecture of cataract susceptibility.

Genome-Wide Analysis of Differentially Expressed miRNAs and Their Associated Regulatory Networks in Lenses Deficient for the Congenital Cataract-Linked Tudor Domain Containing Protein TDRD7.

Mutations/deficiency of TDRD7, encoding a tudor domain protein involved in post-transcriptional gene expression control, causes early onset cataract in humans. While Tdrd7 is implicated in the control of key lens mRNAs, the impact of Tdrd7 deficiency on microRNAs (miRNAs) and how this contributes to transcriptome misexpression and to cataracts, is undefined. We address this critical knowledge-gap by investigating Tdrd7-targeted knockout (Tdrd7-/-) mice that exhibit fully penetrant juvenile cataracts. We performed Affymetrix miRNA 3.0 microarray analysis on Tdrd7-/- mouse lenses at postnatal day (P) 4, a stage preceding cataract formation. This analysis identifies 22 miRNAs [14 over-expressed (miR-15a, miR-19a, miR-138, miR-328, miR-339, miR-345, miR-378b, miR-384, miR-467a, miR-1224, miR-1935, miR-1946a, miR-3102, miR-3107), 8 reduced (let-7b, miR-34c, miR-298, miR-382, miR-409, miR-1198, miR-1947, miR-3092)] to be significantly misexpressed (fold-change ≥ ± 1.2, p-value < 0.05) in Tdrd7-/- lenses. To understand how these misexpressed miRNAs impact Tdrd7-/- cataract, we predicted their mRNA targets and examined their misexpression upon Tdrd7-deficiency by performing comparative transcriptomics analysis on P4 and P30 Tdrd7-/- lens. To prioritize these target mRNAs, we used various stringency filters (e.g., fold-change in Tdrd7-/- lens, iSyTE-based lens-enriched expression) and identified 98 reduced and 89 elevated mRNA targets for overexpressed and reduced miRNAs, respectively, which were classified as "top-priority" "high-priority," and "promising" candidates. For Tdrd7-/- lens overexpressed miRNAs, this approach identified 18 top-priority reduced target mRNAs: Alad, Ankrd46, Ceacam10, Dgat2, Ednrb, H2-Eb1, Klhl22, Lin7a, Loxl1, Lpin1, Npc1, Olfm1, Ppm1e, Ppp1r1a, Rgs8, Shisa4, Snx22 and Wnk2. Majority of these targets were also altered in other gene-specific perturbation mouse models (e.g., Brg1, E2f1/E2f2/E2f3, Foxe3, Hsf4, Klf4, Mafg/Mafk, Notch) of lens defects/cataract, suggesting their importance to lens biology. Gene ontology (GO) provided further insight into their relevance to lens pathology. For example, the Tdrd7-deficient lens capsule defect may be explained by reduced mRNA targets (e.g., Col4a3, Loxl1, Timp2, Timp3) associated with "basement membrane". GO analysis also identified new genes (e.g., Casz1, Rasgrp1) recently linked to lens biology/pathology. Together, these analyses define a new Tdrd7-downstream miRNA-mRNA network, in turn, uncovering several new mRNA targets and their associated pathways relevant to lens biology and offering molecular insights into the pathology of congenital cataract.

The Tudor-domain protein TDRD7, mutated in congenital cataract, controls the heat shock protein HSPB1 (HSP27) and lens fiber cell morphology.

Mutations of the RNA granule component TDRD7 (OMIM: 611258) cause pediatric cataract. We applied an integrated approach to uncover the molecular pathology of cataract in Tdrd7-/- mice. Early postnatal Tdrd7-/- animals precipitously develop cataract suggesting a global-level breakdown/misregulation of key cellular processes. High-throughput RNA sequencing integrated with iSyTE-bioinformatics analysis identified the molecular chaperone and cytoskeletal modulator, HSPB1, among high-priority downregulated candidates in Tdrd7-/- lens. A protein fluorescence two-dimensional difference in-gel electrophoresis (2D-DIGE)-coupled mass spectrometry screen also identified HSPB1 downregulation, offering independent support for its importance to Tdrd7-/- cataractogenesis. Lens fiber cells normally undergo nuclear degradation for transparency, posing a challenge: how is their cell morphology, also critical for transparency, controlled post-nuclear degradation? HSPB1 functions in cytoskeletal maintenance, and its reduction in Tdrd7-/- lens precedes cataract, suggesting cytoskeletal defects may contribute to Tdrd7-/- cataract. In agreement, scanning electron microscopy (SEM) revealed abnormal fiber cell morphology in Tdrd7-/- lenses. Further, abnormal phalloidin and wheat germ agglutinin (WGA) staining of Tdrd7-/- fiber cells, particularly those exhibiting nuclear degradation, reveals distinct regulatory mechanisms control F-actin cytoskeletal and/or membrane maintenance in post-organelle degradation maturation stage fiber cells. Indeed, RNA immunoprecipitation identified Hspb1 mRNA in wild-type lens lysate TDRD7-pulldowns, and single-molecule RNA imaging showed co-localization of TDRD7 protein with cytoplasmic Hspb1 mRNA in differentiating fiber cells, suggesting that TDRD7-ribonucleoprotein complexes may be involved in optimal buildup of key factors. Finally, Hspb1 knockdown in Xenopus causes eye/lens defects. Together, these data uncover TDRD7's novel upstream role in elevation of stress-responsive chaperones for cytoskeletal maintenance in post-nuclear degradation lens fiber cells, perturbation of which causes early-onset cataracts.

Six2 regulates Pax9 expression, palatogenesis and craniofacial bone formation.

In this study, we investigated the role of the transcription factor Six2 in palate development. Six2 was selected using the SysFACE tool to predict genes from the 2p21 locus, a region associated with clefting in humans by GWAS, that are likely to be involved in palatogenesis. We functionally validated the predicted role of Six2 in palatogenesis by showing that 22% of Six2 null embryos develop cleft palate. Six2 contributes to palatogenesis by promoting mesenchymal cell proliferation and regulating bone formation. The clefting phenotype in Six2-/- embryos is similar to Pax9 null embryos, so we examined the functional relationship of these two genes. Mechanistically, SIX2 binds to a PAX9 5' upstream regulatory element and activates PAX9 expression. In addition, we identified a human SIX2 coding variant (p.Gly264Glu) in a proband with cleft palate. We show this missense mutation affects the stability of the SIX2 protein and leads to decreased PAX9 expression. The low penetrance of clefting in the Six2 null mouse combined with the mutation in one patient with cleft palate underscores the potential combinatorial interactions of other genes in clefting. Our study demonstrates that Six2 interacts with the developmental gene regulatory network in the developing palate.

MS/MS in silico subtraction-based proteomic profiling as an approach to facilitate disease gene discovery: application to lens development and cataract.

While the bioinformatics resource-tool iSyTE (integrated Systems Tool for Eye gene discovery) effectively identifies human cataract-associated genes, it is currently based on just transcriptome data, and thus, it is necessary to include protein-level information to gain greater confidence in gene prioritization. Here, we expand iSyTE through development of a novel proteome-based resource on the lens and demonstrate its utility in cataract gene discovery. We applied high-throughput tandem mass spectrometry (MS/MS) to generate a global protein expression profile of mouse lens at embryonic day (E)14.5, which identified 2371 lens-expressed proteins. A major challenge of high-throughput expression profiling is identification of high-priority candidates among the thousands of expressed proteins. To address this problem, we generated new MS/MS proteome data on mouse whole embryonic body (WB). WB proteome was then used as a reference dataset for performing "in silico WB-subtraction" comparative analysis with the lens proteome, which effectively identified 422 proteins with lens-enriched expression at ≥ 2.5 average spectral counts, ≥ 2.0 fold enrichment (FDR < 0.01) cut-off. These top 20% candidates represent a rich pool of high-priority proteins in the lens including known human cataract-linked genes and many new potential regulators of lens development and homeostasis. This rich information is made publicly accessible through iSyTE (https://research.bioinformatics.udel.edu/iSyTE/), which enables user-friendly visualization of promising candidates, thus making iSyTE a comprehensive tool for cataract gene discovery.

High-throughput transcriptome analysis reveals that the loss of Pten activates a novel NKX6-1/RASGRP1 regulatory module to rescue microphthalmia caused by Fgfr2-deficient lenses.

FGFR signaling is critical to development and disease pathogenesis, initiating phosphorylation-driven signaling cascades, notably the RAS-RAF-MEK-ERK and PI3 K-AKT cascades. PTEN antagonizes FGFR signaling by reducing AKT and ERK activation. Mouse lenses lacking FGFR2 exhibit microphakia and reduced ERK and AKT phosphorylation, widespread apoptosis, and defective lens fiber cell differentiation. In contrast, simultaneous deletion of both Fgfr2 and Pten restores ERK and AKT activation levels as well as lens size, cell survival and aspects of fiber cell differentiation; however, the molecular basis of this "rescue" remains undefined. We performed transcriptomic analysis by RNA sequencing of mouse lenses with conditional deletion of Fgfr2, Pten or both Fgfr2 and Pten, which reveal new molecular mechanisms that uncover how FGFR2 and PTEN signaling interact during development. The FGFR2-deficient lens transcriptome demonstrates overall loss of fiber cell identity with deregulated expression of 1448 genes. We find that ~ 60% of deregulated genes return to normal expression levels in lenses lacking both Fgfr2 and Pten. Further, application of customized filtering parameters to these RNA-seq data sets identified 68 high-priority candidate genes. Bioinformatics analyses showed that the cis-binding motif of a high-priority homeodomain transcription factor, NKX6-1, was present in the putative promoters of ~ 78% of these candidates. Finally, biochemical reporter assays demonstrate that NKX6-1 activated the expression of the high-priority candidate Rasgrp1, a RAS-activating protein. Together, these data define a novel regulatory module in which NKX6-1 directly activates Rasgrp1 expression to restore the balance of ERK and AKT activation, thus providing new insights into alternate regulation of FGFR downstream events.

Missense Pathogenic variants in KIF4A Affect Dental Morphogenesis Resulting in X-linked Taurodontism, Microdontia and Dens-Invaginatus.

The etiology of dental anomalies is multifactorial; and genetic and environmental factors that affect the dental lamina have been implicated. We investigated two families of European ancestry in which males were affected by taurodontism, microdontia and dens invaginatus. In both families, males were related to each other via unaffected females. A linkage analysis was conducted in a New Zealand family, followed by exome sequencing and focused analysis of the X-chromosome. In a US family, exome sequencing of the X-chromosome was followed by Sanger sequencing to conduct segregation analyses. We identified two independent missense variants in KIF4A that segregate in affected males and female carriers. The variant in a New Zealand family (p.Asp371His) predicts the substitution of a residue in the motor domain of the protein while the one in a US family (p.Arg771Lys) predicts the substitution of a residue in the domain that interacts with Protein Regulator of Cytokinesis 1 (PRC1). We demonstrated that the gene is expressed in the developing tooth bud during development, and that the p.Arg771Lys variant influences cell migration in an in vitro assay. These data implicate missense variations in KIF4A in a pathogenic mechanism that causes taurodontism, microdontia and dens invaginatus phenotypes.