Identification of pax6-dependent gene regulatory networks in the mouse lens

Single-Cell RNA Sequencing Analysis of the Early Postnatal Mouse Lens Epithelium

Purpose

The lens epithelium maintains the overall health of the organ. We used single-cell RNA sequencing (scRNA-seq) technology to assess transcriptional heterogeneity between cells in the postnatal day 2 (P2) epithelium and identify distinct epithelial cell subtypes. Analysis of these data was used to better understand lens growth, differentiation, and homeostasis on P2.

Methods

scRNA-seq on P2 mouse lenses was performed using the 10x Genomics Chromium Single Cell 3' Kit (v3.1) and short-read Illumina sequencing. Sequence alignment and preprocessing of data were conducted using 10x Genomics Cell Ranger software. Seurat was employed for preprocessing, quality control, dimensionality reduction, and cell clustering, and Monocle was utilized for trajectory analysis to understand the developmental progression of the lens cells. CellChat and GO analyses were used to explore cell-cell communication networks and signaling interactions.

Results

Lens epithelial cells (LECs) were divided into seven subclusters, classified by specific gene markers. The expression of crystallin, cell-cycle, and metabolic genes was not uniform, indicating distinct functional roles of LECs. Trajectory analysis predicted a bifurcation of differentiating and cycling cells from an Igfbp5+ progenitor pool. We also identified heterogeneity in signaling molecules and pathways, suggesting that cycling and progenitor subclusters have prominent roles in coordinating crosstalk.

Conclusions

scRNA-seq corroborated many known markers of epithelial differentiation and proliferation while providing further insight into the pathways and genes directing these processes. Interestingly, we demonstrated that the developing epithelium can be divided into distinct subpopulations. These clusters reflect the transcriptionally diverse roles of the epithelium in proliferation, signaling, and maintenance.

Identification of HSPA8 as an interacting partner of MAB21L2 and an important factor in eye development

BACKGROUND

Pathogenic variants in human MAB21L2 result in microphthalmia, anophthalmia, and coloboma. The exact molecular function of MAB21L2 is currently unknown. We conducted a series of yeast two-hybrid (Y2H) experiments to determine protein interactomes of normal human and zebrafish MAB21L2/mab21l2 as well as human disease-associated variant MAB21L2-p.(Arg51Gly) using human adult retina and zebrafish embryo libraries.

RESULTS

These screens identified klhl31, tnpo1, TNPO2/tnpo2, KLC2/klc2, and SPTBN1/sptbn1 as co-factors of MAB21L2/mab21l2. Several factors, including hspa8 and hspa5, were found to interact with MAB21L2-p.Arg51Gly but not wild-type MAB21L2/mab21l2 in Y2H screens. Further analyses via 1-by-1 Y2H assays, co-immunoprecipitation, and mass spectrometry revealed that both normal and variant MAB21L2 interact with HSPA5 and HSPA8. In situ hybridization detected co-expression of hspa5 and hspa8 with mab21l2 during eye development in zebrafish. Examination of zebrafish mutant hspa8hi138Tg identified reduced hspa8 expression associated with severe ocular developmental defects, including small eye, coloboma, and anterior segment dysgenesis. To investigate the effects of hspa8 deficiency on the mab21l2Arg51_Phe52del allele, corresponding zebrafish double mutants were generated and found to be more severely affected than single mutant lines.

CONCLUSION

This study identifies heat shock proteins as interacting partners of MAB21L2/mab21l2 and suggests a role for this interaction in vertebrate eye development.

Co-expression and Interaction of Pax6 with Genes and Proteins of Immunological Surveillance in the Brain of Mice

The Pax6 binds to promoter sequence elements of genes involved in immunological surveillance and interacts with Iba1, p53, Ras-GAP, and Sparc in the brain of mice. The Pax6 also affects the expression pattern of genes involved in neurogenesis and neurodegeneration. However, the expression and association of Pax6 in the brain under immunologically challenged conditions are still elusive. Therefore, it has been intended to analyze the association of Pax6 in the immunity of the brain using the immune-challenged Dalton's lymphoma (DL) mice model. The expressions of Pax6, Iba1, and Tmem119 decreased, but expressions of Ifn-γ, Tnf-α, Bdnf, and Tgf-β increased in the brain of immune-challenged mice as compared to the control. The level of co-expression of Pax6 decreased in dual positive cells with Iba1, Tmem119, Sparc, p53, Bdnf, and Tgf-β in the brain of immune-challenged mice. Binding of Pax6 to multiple sites of the promoter sequences of Bdnf and Tgf-β indicates their Pax6-associated differential expression and association with immune responsive gene. The levels of binding of Pax6 to Tmem119, Iba1, Ifn-γ, and Tnf-α got altered during the immune-challenged state as compared to control. Results provide the first evidence of the association of Pax6 in brain-specific immunity.

PAX6 disease models for aniridia

Aniridia is a pan-ocular genetic developmental eye disorder characterized by complete or partial iris and foveal hypoplasia, for which there is no treatment currently. Progressive sight loss can arise from cataracts, glaucoma, and aniridia-related keratopathy, which can be managed conservatively or through surgical intervention. The vast majority of patients harbor heterozygous mutations involving the PAX6 gene, which is considered the master transcription factor of early eye development. Over the past decades, several disease models have been investigated to gain a better understanding of the molecular pathophysiology, including several mouse and zebrafish strains and, more recently, human-induced pluripotent stem cells (hiPSCs) derived from aniridia patients. The latter provides a more faithful cellular system to study early human eye development. This review outlines the main aniridia-related animal and cellular models used to study aniridia and highlights the key discoveries that are bringing us closer to a therapy for patients.

Insights into the biochemical and biophysical mechanisms mediating the longevity of the transparent optics of the eye lens

In the human eye, a transparent cornea and lens combine to form the "refracton" to focus images on the retina. This requires the refracton to have a high refractive index "n," mediated largely by extracellular collagen fibrils in the corneal stroma and the highly concentrated crystallin proteins in the cytoplasm of the lens fiber cells. Transparency is a result of short-range order in the spatial arrangement of corneal collagen fibrils and lens crystallins, generated in part by post-translational modifications (PTMs). However, while corneal collagen is remodeled continuously and replaced, lens crystallins are very long-lived and are not replaced and so accumulate PTMs over a lifetime. Eventually, a tipping point is reached when protein aggregation results in increased light scatter, inevitably leading to the iconic protein condensation-based disease, age-related cataract (ARC). Cataracts account for 50% of vision impairment worldwide, affecting far more people than other well-known protein aggregation-based diseases. However, because accumulation of crystallin PTMs begins before birth and long before ARC presents, we postulate that the lens protein PTMs contribute to a "cataractogenic load" that not only increases with age but also has protective effects on optical function by stabilizing lens crystallins until a tipping point is reached. In this review, we highlight decades of experimental findings that support the potential for PTMs to be protective during normal development. We hypothesize that ARC is preventable by protecting the biochemical and biophysical properties of lens proteins needed to maintain transparency, refraction, and optical function.

Cytokine Receptor-Like Factor 3 (CRLF3) Contributes to Early Zebrafish Hematopoiesis

Cytokine receptor-like factor 3 (CRLF3) is an ancient protein conserved across metazoans that contains an archetypal cytokine receptor homology domain (CHD). This domain is found in cytokine receptors present in bilateria, including higher vertebrates, that play key roles in a variety of developmental and homeostatic processes, particularly relating to blood and immune cells. However, understanding of CRLF3 itself remains very limited. This study aimed to investigate this evolutionarily significant protein by studying its embryonic expression and function in early development, particularly of blood and immune cells, using zebrafish as a model. Expression of crlf3 was identified in mesoderm-derived tissues in early zebrafish embryos, including the somitic mesoderm and both anterior and posterior lateral plate mesoderm. Later expression was observed in the thymus, brain, retina and exocrine pancreas. Zebrafish crlf3 mutants generated by genome editing technology exhibited a significant reduction in primitive hematopoiesis and early definitive hematopoiesis, with decreased early progenitors impacting on multiple lineages. No other obvious phenotypes were observed in the crlf3 mutants.

Changes in DNA methylation hallmark alterations in chromatin accessibility and gene expression for eye lens differentiation

Background

Methylation at cytosines (mCG) is a well-known regulator of gene expression, but its requirements for cellular differentiation have yet to be fully elucidated. A well-studied cellular differentiation model system is the eye lens, consisting of a single anterior layer of epithelial cells that migrate laterally and differentiate into a core of fiber cells. Here, we explore the genome-wide relationships between mCG methylation, chromatin accessibility and gene expression during differentiation of eye lens epithelial cells into fiber cells.

Results

Whole genome bisulfite sequencing identified 7621 genomic loci exhibiting significant differences in mCG levels between lens epithelial and fiber cells. Changes in mCG levels were inversely correlated with the differentiation state-specific expression of 1285 genes preferentially expressed in either lens fiber or lens epithelial cells (Pearson correlation r = − 0.37, p < 1 × 10^–42). mCG levels were inversely correlated with chromatin accessibility determined by assay for transposase-accessible sequencing (ATAC-seq) (Pearson correlation r = − 0.86, p < 1 × 10^–300). Many of the genes exhibiting altered regions of DNA methylation, chromatin accessibility and gene expression levels in fiber cells relative to epithelial cells are associated with lens fiber cell structure, homeostasis and transparency. These include lens crystallins (CRYBA4, CRYBB1, CRYGN, CRYBB2), lens beaded filament proteins (BFSP1, BFSP2), transcription factors (HSF4, SOX2, HIF1A), and Notch signaling pathway members (NOTCH1, NOTCH2, HEY1, HES5). Analysis of regions exhibiting cell-type specific alterations in DNA methylation revealed an overrepresentation of consensus sequences of multiple transcription factors known to play key roles in lens cell differentiation including HIF1A, SOX2, and the MAF family of transcription factors.

Conclusions

Collectively, these results link DNA methylation with control of chromatin accessibility and gene expression changes required for eye lens differentiation. The results also point to a role for DNA methylation in the regulation of transcription factors previously identified to be important for lens cell differentiation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-022-00440-z.

Identification of missense MAB21L1 variants in microphthalmia and aniridia

Microphthalmia, coloboma, and aniridia are congenital ocular phenotypes with a strong genetic component but often unknown cause. We present a likely causative novel variant in MAB21L1, c.152G>T p.(Arg51Leu), in two family members with microphthalmia and aniridia, as well as novel or rare compound heterozygous variants of uncertain significance, c.184C>T p.(Arg62Cys)/c.-68T>C, and c.658G>C p.(Gly220Arg)/c.*529A>G, in two additional probands with microphthalmia, coloboma and/or cataracts. All variants were predicted as damaging by in silico programs. In vitro studies of coding variants revealed normal subcellular localization but variable stability for the corresponding mutant proteins. In vivo complementation assays using the zebrafish mab21l2 Q48Sfs*5 loss-of-function line demonstrated that though overexpression of wild-type MAB21L1 messenger RNA (mRNA) compensated for the loss of mab21l2, none of the coding variant mRNAs produced a statistically significant rescue, with p.(Arg51Leu) showing the highest degree of functional deficiency. Dominant variants in a close homolog of MAB21L1, MAB21L2, have been associated with microphthalmia and/or coloboma and repeatedly involved the same Arg51 residue, further supporting its pathogenicity. The possible role of p.(Arg62Cys) and p.(Gly220Arg) in microphthalmia is similarly supported by the observed functional defects, with or without an additional impact from noncoding MAB21L1 variants identified in each patient. This study suggests a broader spectrum of MAB21L1-associated disease.

The Pax6 master control gene initiates spontaneous retinal development via a self-organising Turing network

Understanding how complex organ systems are assembled from simple embryonic tissues is a major challenge. Across the animal kingdom a great diversity of visual organs are initiated by a 'master control gene' called Pax6, which is both necessary and sufficient for eye development. Yet precisely how Pax6 achieves this deeply homologous function is poorly understood. Using the chick as a model organism, we show that vertebrate Pax6 interacts with a pair of morphogen-coding genes, Tgfb2 and Fst, to form a putative Turing network, which we have computationally modelled. Computer simulations suggest that this gene network is sufficient to spontaneously polarise the developing retina, establishing the first organisational axis of the eye and prefiguring its further development. Our findings reveal how retinal self-organisation may be initiated independently of the highly ordered tissue interactions that help to assemble the eye in vivo These results help to explain how stem cell aggregates spontaneously self-organise into functional eye-cups in vitro We anticipate these findings will help to underpin retinal organoid technology, which holds much promise as a platform for disease modelling, drug development and regenerative therapies.

Insights into the regulatory molecules involved in glaucoma pathogenesis

Glaucoma is an important cause of irreversible blindness, characterized by optic nerve anomalies. Increased intraocular pressure (IOP) and aging are major risk factors. Retinal ganglion cells and trabecular meshwork cells are certainly involved in the etiology of glaucoma. Glaucoma is usually a complex disease, and various genes and functions may contribute to its etiology. Among these may be genes that encode regulatory molecules. In this review, regulatory molecules including 18 transcription factors (TFs), 195 microRNAs (miRNAs), 106 long noncoding RNAs (lncRNAs), and two circular RNAs (circRNAs) that are reasonable candidates for having roles in glaucoma pathogenesis are described. The targets of the regulators are reported. Glaucoma-related features including apoptosis, stress responses, immune functions, ECM properties, IOP, and eye development are affected by the targeted genes. The targeted genes that are frequently targeted by multiple regulators most often affect apoptosis and the related features of cell death and cell survival. BCL2, CDKN1A, and TP53 are among the frequent targets of three types of glaucoma-relevant regulators, TFs, miRNAs, and lncRNAs. TP53 was itself identified as a glaucoma-relevant TF. Several of the glaucoma-relevant TFs are themselves among frequent targets of regulatory molecules, which is consistent with existence of a complex network involved in glaucoma pathogenesis.

Pax6 organizes the anterior eye segment by guiding two distinct neural crest waves

Cranial neural crest (NC) contributes to the developing vertebrate eye. By multidimensional, quantitative imaging, we traced the origin of the ocular NC cells to two distinct NC populations that differ in the maintenance of sox10 expression, Wnt signalling, origin, route, mode and destination of migration. The first NC population migrates to the proximal and the second NC cell group populates the distal (anterior) part of the eye. By analysing zebrafish pax6a/b compound mutants presenting anterior segment dysgenesis, we demonstrate that Pax6a/b guide the two NC populations to distinct proximodistal locations. We further provide evidence that the lens whose formation is pax6a/b-dependent and lens-derived TGFβ signals contribute to the building of the anterior segment. Taken together, our results reveal multiple roles of Pax6a/b in the control of NC cells during development of the anterior segment.

Transcriptome analysis of neural progenitor cells derived from Lowe syndrome induced pluripotent stem cells: identification of candidate genes for the neurodevelopmental and eye manifestations

BACKGROUND

Lowe syndrome (LS) is caused by loss-of-function mutations in the X-linked gene OCRL, which codes for an inositol polyphosphate 5-phosphatase that plays a key role in endosome recycling, clathrin-coated pit formation, and actin polymerization. It is characterized by congenital cataracts, intellectual and developmental disability, and renal proximal tubular dysfunction. Patients are also at high risk for developing glaucoma and seizures. We recently developed induced pluripotent stem cell (iPSC) lines from three patients with LS who have hypomorphic variants affecting the 3' end of the gene, and their neurotypical brothers to serve as controls.

METHODS

In this study, we used RNA sequencing (RNA-seq) to obtain transcriptome profiles in LS and control neural progenitor cells (NPCs).

RESULTS

In a comparison of the patient and control NPCs (n = 3), we found 16 differentially expressed genes (DEGs) at the multiple test adjusted p value (padj) < 0.1, with nine at padj < 0.05. Using nominal p value < 0.05, 319 DEGs were detected. The relatively small number of DEGs could be due to the fact that OCRL is not a transcription factor per se, although it could have secondary effects on gene expression through several different mechanisms. Although the number of DEGs passing multiple test correction was small, those that were found are quite consistent with some of the known molecular effects of OCRL protein, and the clinical manifestations of LS. Furthermore, using gene set enrichment analysis (GSEA), we found that genes increased expression in the patient NPCs showed enrichments of several gene ontology (GO) terms (false discovery rate < 0.25): telencephalon development, pallium development, NPC proliferation, and cortex development, which are consistent with a condition characterized by intellectual disabilities and psychiatric manifestations. In addition, a significant enrichment among the nominal DEGs for genes implicated in autism spectrum disorder (ASD) was found (e.g., AFF2, DNER, DPP6, DPP10, RELN, CACNA1C), as well as several that are strong candidate genes for the development of eye problems found in LS, including glaucoma. The most notable example is EFEMP1, a well-known candidate gene for glaucoma and other eye pathologies.

CONCLUSION

Overall, the RNA-seq findings present several candidate genes that could help explain the underlying basis for the neurodevelopmental and eye problems seen in boys with LS.

Metabolic and pathologic profiles of human LSS deficiency recapitulated in mice

Skin lesions, cataracts, and congenital anomalies have been frequently associated with inherited deficiencies in enzymes that synthesize cholesterol. Lanosterol synthase (LSS) converts (S)-2,3-epoxysqualene to lanosterol in the cholesterol biosynthesis pathway. Biallelic mutations in LSS have been reported in families with congenital cataracts and, very recently, have been reported in cases of hypotrichosis. However, it remains to be clarified whether these phenotypes are caused by LSS enzymatic deficiencies in each tissue, and disruption of LSS enzymatic activity in vivo has not yet been validated. We identified two patients with novel biallelic LSS mutations who exhibited congenital hypotrichosis and midline anomalies but did not have cataracts. We showed that the blockade of the LSS enzyme reaction occurred in the patients by measuring the (S)-2,3-epoxysqualene/lanosterol ratio in the forehead sebum, which would be a good biomarker for the diagnosis of LSS deficiency. Epidermis-specific Lss knockout mice showed neonatal lethality due to dehydration, indicating that LSS could be involved in skin barrier integrity. Tamoxifen-induced knockout of Lss in the epidermis caused hypotrichosis in adult mice. Lens-specific Lss knockout mice had cataracts. These results confirmed that LSS deficiency causes hypotrichosis and cataracts due to loss-of-function mutations in LSS in each tissue. These mouse models will lead to the elucidation of the pathophysiological mechanisms associated with disrupted LSS and to the development of therapeutic treatments for LSS deficiency.

Identification of genes involved in glaucoma pathogenesis using combined network analysis and empirical studies

Glaucoma is a leading cause of blindness. We aimed in this study to identify genes that may make subtle and cumulative contributions to glaucoma pathogenesis. To this end, we identified molecular interactions and pathways that include transcription factors (TFs) FOXC1, PITX2, PAX6 and NFKB1 and various microRNAs including miR-204 known to have relevance to trabecular meshwork (TM) functions and/or glaucoma. TM tissue is involved in glaucoma pathogenesis. In-house microarray transcriptome results and data sources were used to identify target genes of the regulatory molecules. Bioinformatics analyses were done to filter TM and glaucoma relevant genes. These were submitted to network-creating softwares to define interactions, pathways and a network that would include the genes. The network was stringently scrutinized and minimized, then expanded by addition of microarray data and data on TF and microRNA-binding sites. Selected features of the network were confirmed by empirical studies such as dual luciferase assays, real-time PCR and western blot experiments and apoptosis assays. MYOC, WDR36, LTPBP2, RHOA, CYP1B1, OPA1, SPARC, MEIS2, PLEKHG5, RGS5, BBS5, ALDH1A1, NOMO2, CXCL6, FMNL2, ADAMTS5, CLOCK and DKK1 were among the genes included in the final network. Pathways identified included those that affect ECM properties, IOP, ciliary body functions, retinal ganglion cell viability, apoptosis, focal adhesion and oxidative stress response. The identification of many genes potentially involved in glaucoma pathology is consistent with its being a complex disease. The inclusion of several known glaucoma-related genes validates the approach used.

VCAN Canonical Splice Site Mutation is Associated With Vitreoretinal Degeneration and Disrupts an MMP Proteolytic Site

Purpose

To gain insight into the pathophysiology of vitreoretinal degeneration, the clinical course of three family members with Versican Vitreoretinopathy (VVR) is described, and a canonical splice site mutation in the gene encoding for versican (VCAN) protein was biochemically analyzed.

Methods

A retrospective chart review, human eye histopathology, Sanger DNA sequencing, protein structural modeling, and in vitro proteolysis assays were performed.

Results

The proband (II:1), mother (I:2), and younger sibling (II:2) suffered retinal degeneration with foveal sparing and retinal detachments with proliferative vitreoretinopathy, features that were confirmed on histopathologic analysis. All affected members carried a heterozygous adenine to guanine variant (c.4004-2A>G) predicted to result in exon 8 skipping or the deletion of 13 amino acids at the beginning of the GAGβ chain (VCAN p.1335-1347). This deleted region corresponded to a putative MMP cleavage site, validated using fluorescence resonance energy transfer (FRET)-based proteolysis assays. Proteomic network analysis identified 10 interacting partners in the human vitreous and retina linked to retinal detachment and degeneration.

Conclusions

VVR causes significant ocular disease, including retinal detachment and retinal dystrophy. The intronic VCAN mutation removes an MMP cleavage site, which alters versican structure and results in abnormal vitreous modeling. Disruption of a versican protein network may underlie clinicopathologic disease features and point to targeted therapies.

Loss of Pax6 Causes Regional Changes in Dll1 Expression in Developing Cerebral Cortex

The transcription factor Pax6 controls multiple aspects of forebrain development. Conditional deletion of Pax6 in embryonic mouse cortex causes increased proliferation of cortical progenitor cells and a concomitant decrease in neural differentiation. Notch signaling regulates the balance between proliferation and differentiation of cortical progenitor cells, suggesting a possible connection between Pax6 and Notch signaling. We investigated how expression of the Notch ligand delta-like 1 (Dll1) is altered by loss of Pax6. Acute cortex-specific deletion of Pax6 resulted in a widespread decrease in the density of Dll1+ cells at embryonic days 12.5 and 13.5 (E12.5 and E13.5). In constitutive loss-of-function mutants, decreases in the densities of Dll1+ cells were more limited both spatially and temporally. Controlled over-expression of Pax6 had no detectable effect on Dll1 expression. The proneural transcription factor Neurog2 is a target of Pax6 that can activate Dll1 expression and we found clear co-expression of Neurog2 and Dll1 in radial glial progenitors, suggesting that Pax6’s effect on Dll1 could be mediated through Neurog2. However, we found no change in Dll1+ cells in Neurog2^−/− cortex suggesting either that Neurog2 is not directly involved, or that its loss of function in embryonic cortex can be compensated for.

Transcription factor PAX6 as a novel prognostic factor and putative tumour suppressor in non-small cell lung cancer

Lung cancer is the leading cause of cancer deaths. Novel predictive biomarkers are needed to improve treatment selection and more accurate prognostication. PAX6 is a transcription factor with a proposed tumour suppressor function. Immunohistochemical staining was performed on tissue microarrays from 335 non-small cell lung cancer (NSCLC) patients for PAX6. Multivariate analyses of clinico-pathological variables and disease-specific survival (DSS) was carried out, and phenotypic changes of two NSCLC cell lines with knockdown of PAX6 were characterized. While PAX6 expression was only associated with a trend of better disease-specific survival (DSS) (p = 0.10), the pN+ subgroup (N = 103) showed significant correlation between high PAX6 expression and longer DSS (p = 0.022). Median survival for pN + patients with high PAX6 expression was 127.4 months, versus 22.9 months for patients with low PAX6 expression. In NCI-H661 cells, knockdown of PAX6 strongly activated serum-stimulated migration. In NCI-H460 cells, PAX6 knockdown activated anchorage-independent growth. We did not observe any significant effect of PAX6 on proliferation in either of cell lines. Our findings strongly support the proposition of PAX6 as a valid and positive prognostic marker in NSCLC in node-positive patients. There is a need for further studies, which should provide mechanistical explanation for the role of PAX6 in NSCLC.

Efficacy of Postnatal In Vivo Nonsense Suppression Therapy in a Pax6 Mouse Model of Aniridia

Nonsense mutations leading to premature stop codons are common occurring in approximately 12% of all human genetic diseases. Thus, pharmacological nonsense mutation suppression strategies would be beneficial to a large number of patients if the drugs could be targeted to the affected tissues at the appropriate time. Here, we used nonsense suppression to manipulate Pax6 dosage at different developmental times in the eye of the small eye (Pax6^Sey/+; G194X) mouse model of aniridia. Efficacy was assessed by functional assays for visual capacity, including electroretinography and optokinetic tracking (OKT), in addition to histological and biochemical studies. Malformation defects in the Pax6^Sey/+ postnatal eye responded to topically delivered nonsense suppression in a dose- and time-dependent manner. Elevated levels of Mmp9, a direct downstream target of Pax6 in the cornea, were observed with the different treatment regimens. The lens capsule was particularly sensitive to Pax6 dosage, revealing a potential new role for Pax6 in lens capsule maintenance and development. The remarkable capacity of malformed ocular tissue to respond postnatally to Pax6 dosage in vivo demonstrates that the use of nonsense suppression could be a valuable therapeutic approach for blinding diseases caused by nonsense mutations.

Control of lens development by Lhx2-regulated neuroretinal FGFs

Fibroblast growth factor (FGF) signaling is an essential regulator of lens epithelial cell proliferation and survival, as well as lens fiber cell differentiation. However, the identities of these FGF factors, their source tissue and the genes that regulate their synthesis are unknown. We have found that Chx10-Cre;Lhx2lox/lox mice, which selectively lack Lhx2 expression in neuroretina from E10.5, showed an early arrest in lens fiber development along with severe microphthalmia. These mutant animals showed reduced expression of multiple neuroretina-expressed FGFs and canonical FGF-regulated genes in neuroretina. When FGF expression was genetically restored in Lhx2-deficient neuroretina of Chx10-Cre;Lhx2lox/lox mice, we observed a partial but nonetheless substantial rescue of the defects in lens cell proliferation, survival and fiber differentiation. These data demonstrate that neuroretinal expression of Lhx2 and neuroretina-derived FGF factors are crucial for lens fiber development in vivo.

Pax6 associates with H3K4-specific histone methyltransferases Mll1, Mll2, and Set1a and regulates H3K4 methylation at promoters and enhancers

BACKGROUND

Pax6 is a key regulator of the entire cascade of ocular lens formation through specific binding to promoters and enhancers of batteries of target genes. The promoters and enhancers communicate with each other through DNA looping mediated by multiple protein-DNA and protein-protein interactions and are marked by specific combinations of histone posttranslational modifications (PTMs). Enhancers are distinguished from bulk chromatin by specific modifications of core histone H3, including H3K4me1 and H3K27ac, while promoters show increased H3K4me3 PTM. Previous studies have shown the presence of Pax6 in as much as 1/8 of lens-specific enhancers but a much smaller fraction of tissue-specific promoters. Although Pax6 is known to interact with EP300/p300 histone acetyltransferase responsible for generation of H3K27ac, a potential link between Pax6 and histone H3K4 methylation remains to be established.

RESULTS

Here we show that Pax6 co-purifies with H3K4 methyltransferase activity in lens cell nuclear extracts. Proteomic studies show that Pax6 immunoprecipitates with Set1a, Mll1, and Mll2 enzymes, and their associated proteins, i.e., Wdr5, Rbbp5, Ash2l, and Dpy30. ChIP-seq studies using chromatin prepared from mouse lens and cultured lens cells demonstrate that Pax6-bound regions are mostly enriched with H3K4me2 and H3K4me1 in enhancers and promoters, though H3K4me3 marks only Pax6-containing promoters. The shRNA-mediated knockdown of Pax6 revealed down-regulation of a set of direct target genes, including Cap2, Farp1, Pax6, Plekha1, Prox1, Tshz2, and Zfp536. Pax6 knockdown was accompanied by reduced H3K4me1 at enhancers and H3K4me3 at promoters, with little or no changes of the H3K4me2 modifications. These changes were prominent in Plekha1, a gene regulated by Pax6 in both lens and retinal pigmented epithelium.

CONCLUSIONS

Our study supports a general model of Pax6-mediated recruitment of histone methyltransferases Mll1 and Mll2 to lens chromatin, especially at distal enhancers. Genome-wide data in lens show that Pax6 binding correlates with H3K4me2, consistent with the idea that H3K4me2 PTMs correlate with the binding of transcription factors. Importantly, partial reduction of Pax6 induces prominent changes in local H3K4me1 and H3K4me3 modification. Together, these data open the field to mechanistic studies of Pax6, Mll1, Mll2, and H3K4me1/2/3 dynamics at distal enhancers and promoters of developmentally controlled genes.

Chromatin remodeling enzyme Snf2h regulates embryonic lens differentiation and denucleation

Ocular lens morphogenesis is a model for investigating mechanisms of cellular differentiation, spatial and temporal gene expression control, and chromatin regulation. Brg1 (Smarca4) and Snf2h (Smarca5) are catalytic subunits of distinct ATP-dependent chromatin remodeling complexes implicated in transcriptional regulation. Previous studies have shown that Brg1 regulates both lens fiber cell differentiation and organized degradation of their nuclei (denucleation). Here, we employed a conditional Snf2h(flox) mouse model to probe the cellular and molecular mechanisms of lens formation. Depletion of Snf2h induces premature and expanded differentiation of lens precursor cells forming the lens vesicle, implicating Snf2h as a key regulator of lens vesicle polarity through spatial control of Prox1, Jag1, p27(Kip1) (Cdkn1b) and p57(Kip2) (Cdkn1c) gene expression. The abnormal Snf2h(-/-) fiber cells also retain their nuclei. RNA profiling of Snf2h(-/) (-) and Brg1(-/-) eyes revealed differences in multiple transcripts, including prominent downregulation of those encoding Hsf4 and DNase IIβ, which are implicated in the denucleation process. In summary, our data suggest that Snf2h is essential for the establishment of lens vesicle polarity, partitioning of prospective lens epithelial and fiber cell compartments, lens fiber cell differentiation, and lens fiber cell nuclear degradation.

PAX6: 25th anniversary and more to learn

The DNA-binding transcription factor PAX6 was cloned 25 years ago by multiple teams pursuing identification of human and mouse eye disease causing genes, cloning vertebrate homologues of pattern-forming regulatory genes identified in Drosophila, or abundant eye-specific transcripts. Since its discovery in 1991, genetic, cellular, molecular and evolutionary studies on Pax6 mushroomed in the mid 1990s leading to the transformative thinking regarding the genetic program orchestrating both early and late stages of eye morphogenesis as well as the origin and evolution of diverse visual systems. Since Pax6 is also expressed outside of the eye, namely in the central nervous system and pancreas, a number of important insights into the development and function of these organs have been amassed. In most recent years, genome-wide technologies utilizing massively parallel DNA sequencing have begun to provide unbiased insights into the regulatory hierarchies of specification, determination and differentiation of ocular cells and neurogenesis in general. This review is focused on major advancements in studies on mammalian eye development driven by studies of Pax6 genes in model organisms and future challenges to harness the technology-driven opportunities to reconstruct, step-by-step, the transition from naïve ectoderm, neuroepithelium and periocular mesenchyme/neural crest cells into the three-dimensional architecture of the eye.

Systems biology of lens development: A paradigm for disease gene discovery in the eye

Over the past several decades, the biology of the developing lens has been investigated using molecular genetics-based approaches in various vertebrate model systems. These efforts, involving target gene knockouts or knockdowns, have led to major advances in our understanding of lens morphogenesis and the pathological basis of cataracts, as well as of other lens related eye defects. In particular, we now have a functional understanding of regulators such as Pax6, Six3, Sox2, Oct1 (Pou2f1), Meis1, Pnox1, Zeb2 (Sip1), Mab21l1, Foxe3, Tfap2a (Ap2-alpha), Pitx3, Sox11, Prox1, Sox1, c-Maf, Mafg, Mafk, Hsf4, Fgfrs, Bmp7, and Tdrd7 in this tissue. However, whether these individual regulators interact or their targets overlap, and the significance of such interactions during lens morphogenesis, is not well defined. The arrival of high-throughput approaches for gene expression profiling (microarrays, RNA-sequencing (RNA-seq), etc.), which can be coupled with chromatin immunoprecipitation (ChIP) or RNA immunoprecipitation (RIP) assays, along with improved computational resources and publically available datasets (e.g. those containing comprehensive protein-protein, protein-DNA information), presents new opportunities to advance our understanding of the lens tissue on a global systems level. Such systems-level knowledge will lead to the derivation of the underlying lens gene regulatory network (GRN), defined as a circuit map of the regulator-target interactions functional in lens development, which can be applied to expedite cataract gene discovery. In this review, we cover the various systems-level approaches such as microarrays, RNA-seq, and ChIP that are already being applied to lens studies and discuss strategies for assembling and interpreting these vast amounts of high-throughput information for effective dispersion to the scientific community. In particular, we discuss strategies for effective interpretation of this new information in the context of the rich knowledge obtained through the application of traditional single-gene focused experiments on the lens. Finally, we discuss our vision for integrating these diverse high-throughput datasets in a single web-based user-friendly tool iSyTE (integrated Systems Tool for Eye gene discovery) - a resource that is already proving effective in the identification and characterization of genes linked to lens development and cataract. We anticipate that application of a similar approach to other ocular tissues such as the retina and the cornea, and even other organ systems, will significantly impact disease gene discovery.

ELISA based quantification of Pax6 expression in the developing Zebrafish embryos

BACKGROUND

Transcription factors are the key regulators of metabolic pathways in the cells, tissues and organ development during embryogenesis. Pax6 is a transcription factor involved in vertebrate eye, brain and central nervous system formation during development.

PURPOSE

A reliable and sensitive assay for the spatiotemporal expression, quantification and detection of Pax6 is not available so far in zebrafish as a developmental model, hence the objective of this work is to develop quantitative assay in zebrafish embryos.

METHODS

The Pax6 transcription factor was purified by heparin agarose affinity chromatography and DEAE cellulose chromatography techniques from the developing zebrafish embryos. The purity was confirmed by SDS-PAGE and western blotting using Pax6 mouse monoclonal antibody. The standard graph was plotted for Pax6 and the expressions in seventeen developmental stages were quantified by indirect ELISA.

RESULTS

The maximum expression of Pax6 was detected at 8 hpf (hours post fertilization) and it was quantified as 179 ng/embryo from the average total protein of 9.5 µg/embryo. The zebrafish Pax6 protein was detected as 48 kDa and confirmed by western blotting.

CONCLUSION

This study paves way to quantify the level of expression of proteins or transcription factors during early embryonic and larval development or embryogenesis using zebrafish as model system.

Regulation of cerebral cortical neurogenesis by the Pax6 transcription factor

Understanding brain development remains a major challenge at the heart of understanding what makes us human. The neocortex, in evolutionary terms the newest part of the cerebral cortex, is the seat of higher cognitive functions. Its normal development requires the production, positioning, and appropriate interconnection of very large numbers of both excitatory and inhibitory neurons. Pax6 is one of a relatively small group of transcription factors that exert high-level control of cortical development, and whose mutation or deletion from developing embryos causes major brain defects and a wide range of neurodevelopmental disorders. Pax6 is very highly conserved between primate and non-primate species, is expressed in a gradient throughout the developing cortex and is essential for normal corticogenesis. Our understanding of Pax6’s functions and the cellular processes that it regulates during mammalian cortical development has significantly advanced in the last decade, owing to the combined application of genetic and biochemical analyses. Here, we review the functional importance of Pax6 in regulating cortical progenitor proliferation, neurogenesis, and formation of cortical layers and highlight important differences between rodents and primates. We also review the pathological effects of PAX6 mutations in human neurodevelopmental disorders. We discuss some aspects of Pax6’s molecular actions including its own complex transcriptional regulation, the distinct molecular functions of its splice variants and some of Pax6’s known direct targets which mediate its actions during cortical development.

Mutations in MAB21L2 result in ocular Coloboma, microcornea and cataracts

Ocular coloboma results from abnormal embryonic development and is often associated with additional ocular and systemic features. Coloboma is a highly heterogeneous disorder with many cases remaining unexplained. Whole exome sequencing from two cousins affected with dominant coloboma with microcornea, cataracts, and skeletal dysplasia identified a novel heterozygous allele in MAB21L2, c.151 C>G, p.(Arg51Gly); the mutation was present in all five family members with the disease and appeared de novo in the first affected generation of the three-generational pedigree. MAB21L2 encodes a protein similar to C. elegans mab-21 cell fate-determining factor; the molecular function of MAB21L2 is largely unknown. To further evaluate the role of MAB21L2, zebrafish mutants carrying a p.(Gln48Serfs*5) frameshift truncation (mab21l2Q48Sfs*5) and a p.(Arg51_Phe52del) in-frame deletion (mab21l2R51_F52del) were developed with TALEN technology. Homozygous zebrafish embryos from both lines developed variable lens and coloboma phenotypes: mab21l2Q48Sfs*5 embryos demonstrated severe lens and retinal defects with complete lethality while mab21l2R51_F52del mutants displayed a milder lens phenotype and severe coloboma with a small number of fish surviving to adulthood. Protein studies showed decreased stability for the human p.(Arg51Gly) and zebrafish p.(Arg51_Phe52del) mutant proteins and predicted a complete loss-of-function for the zebrafish p.(Gln48Serfs*5) frameshift truncation. Additionally, in contrast to wild-type human MAB21L2 transcript, mutant p.(Arg51Gly) mRNA failed to efficiently rescue the ocular phenotype when injected into mab21l2Q48Sfs*5 embryos, suggesting this allele is functionally deficient. Histology, immunohistochemistry, and in situ hybridization experiments identified retinal invagination defects, an increase in cell death, abnormal proliferation patterns, and altered expression of several ocular markers in the mab21l2 mutants. These findings support the identification of MAB21L2 as a novel factor involved in human coloboma and highlight the power of genome editing manipulation in model organisms for analysis of the effects of whole exome variation in humans.

Xenopus pax6 mutants affect eye development and other organ systems, and have phenotypic similarities to human aniridia patients

Mutations in the Pax6 gene cause ocular defects in both vertebrate and invertebrate animal species, and the disease aniridia in humans. Despite extensive experimentation on this gene in multiple species, including humans, we still do not understand the earliest effects on development mediated by this gene. This prompted us to develop pax6 mutant lines in Xenopus tropicalis taking advantage of the utility of the Xenopus system for examining early development and in addition to establish a model for studying the human disease aniridia in an accessible lower vertebrate. We have generated mutants in pax6 by using Transcription Activator-Like Effector Nuclease (TALEN) constructs for gene editing in X. tropicalis. Embryos with putative null mutations show severe eye abnormalities and changes in brain development, as assessed by changes in morphology and gene expression. One gene that we found is downregulated very early in development in these pax6 mutants is myc, a gene involved in pluripotency and progenitor cell maintenance and likely a mediator of some key pax6 functions in the embryo. Changes in gene expression in the developing brain and pancreas reflect other important functions of pax6 during development. In mutations with partial loss of pax6 function eye development is initially relatively normal but froglets show an underdeveloped iris, similar to the classic phenotype (aniridia) seen in human patients with PAX6 mutations. Other eye abnormalities observed in these froglets, including cataracts and corneal defects, are also common in human aniridia. The frog model thus allows us to examine the earliest deficits in eye formation as a result of pax6 lesions, and provides a useful model for understanding the developmental basis for the aniridia phenotype seen in humans.

The cellular and molecular mechanisms of vertebrate lens development

The ocular lens is a model system for understanding important aspects of embryonic development, such as cell specification and the spatiotemporally controlled formation of a three-dimensional structure. The lens, which is characterized by transparency, refraction and elasticity, is composed of a bulk mass of fiber cells attached to a sheet of lens epithelium. Although lens induction has been studied for over 100 years, recent findings have revealed a myriad of extracellular signaling pathways and gene regulatory networks, integrated and executed by the transcription factor Pax6, that are required for lens formation in vertebrates. This Review summarizes recent progress in the field, emphasizing the interplay between the diverse regulatory mechanisms employed to form lens progenitor and precursor cells and highlighting novel opportunities to fill gaps in our understanding of lens tissue morphogenesis.

The homeodomain of Eyeless regulates cell growth and antagonizes the paired domain-dependent retinal differentiation function

Pax6 and its Drosophila homolog Eyeless (Ey) play essential roles during eye development. Ey/Pax6 contains two distinct DNA binding domains, a Paired domain (PD) and a Homeodomain (HD). While Ey/Pax6 PD is required for the expression of key regulators of retinal development, relatively little is known about the HD-dependent Ey function. In this study, we used the UAS/GAL4 system to determine the functions of different Ey domains on cell growth and on retinal development. We showed that Ey can promote cell growth, which requires the HD but not the PD. In contrast, the ability of Ey to activate Ato expression and induce ectopic eye formation requires the PD but not the HD. Interestingly, deletion of the HD enhanced Ey-dependent ectopic eye induction while overexpression of the HD only Ey forms antagonizes ectopic eye induction. These studies revealed a novel function of Ey HD on cell growth and a novel antagonistic effect of Ey HD on Ey PD-dependent eye induction. We further show the third helix of the Ey HD can directly interact with the RED subdomain in Ey PD and that deletion of the HD increased the binding of Ey PD to its target. These results suggest that the direct interaction between the HD and the PD potentially mediates their antagonistic effects. Since different Ey splicing forms are expressed in overlapping regions during normal development, we speculate that the expression ratios of the different Ey splice forms potentially contribute to the regulation of growth and differentiation of these tissues.

Transcriptional regulation of important cellular processes in skeletal myogenesis through interferon-γ

The purpose of the present study was to investigate the effect of interferon (IFN)-γ on the transcriptomic profile of differentiating mouse C2C12 myogenic cells. Global gene expression was evaluated using whole mouse genome oligonucleotide microarrays, and the results were validated through real-time PCR. IFN-γ (1 ng/mL) increased myoblast proliferation but decreased cell respiration and myosin heavy chain content and slightly decreased the fusion index in differentiating C2C12 cell cultures. The genes upregulated through IFN-γ were involved in cell cycle; regulation of cell proliferation; programmed cell death; chemotaxis; and cytokine, growth factor, and peptidase activity, whereas the genes downregulated through IFN-γ primarily contributed to the regulation of transcription, cell-cell signaling, nitrogen compound biosynthesis, ser/thr protein kinase signaling, and regulation of the Wnt pathway. In conclusion, IFN-γ affects the expression of numerous genes associated with the regulation of several processes in myogenesis. The effects of IFN-γ on cellular transcription include (1) alteration of cytokine/growth factor expression, promoting cell proliferation and migration but inhibiting differentiation, (2) impairment of pro-myogenic transcription, (3) disruption of cell adhesion and sarcolemma/cytoskeleton organization, and (4) increased peptidase activity leading to enhanced proteolysis and apoptosis.

In vivo analysis of hyaloid vasculature morphogenesis in zebrafish: A role for the lens in maturation and maintenance of the hyaloid

Two vascular networks nourish the embryonic eye as it develops - the hyaloid vasculature, located at the anterior of the eye between the retina and lens, and the choroidal vasculature, located at the posterior of the eye, surrounding the optic cup. Little is known about hyaloid development and morphogenesis, however. To begin to identify the morphogenetic underpinnings of hyaloid formation, we utilized in vivo time-lapse confocal imaging to characterize morphogenesis of the zebrafish hyaloid through 5 days post fertilization (dpf). Our data segregate hyaloid formation into three distinct morphogenetic stages: Stage I: arrival of hyaloid cells at the lens and formation of the hyaloid loop; Stage II: formation of a branched hyaloid network; Stage III: refinement of the hyaloid network. Utilizing fixed and dissected tissues, distinct Stage II and Stage III aspects of hyaloid formation were quantified over time. Combining in vivo imaging with microangiography, we demonstrate that the hyaloid system becomes fully enclosed by 5dpf. To begin to identify the molecular and cellular mechanisms underlying hyaloid morphogenesis, we identified a recessive mutation in the mab21l2 gene, and in a subset of mab21l2 mutants the lens does not form. Utilizing these "lens-less" mutants, we determined whether the lens was required for hyaloid morphogenesis. Our data demonstrate that the lens is not required for Stage I of hyaloid formation; however, Stages II and III of hyaloid formation are disrupted in the absence of a lens, supporting a role for the lens in hyaloid maturation and maintenance. Taken together, this study provides a foundation on which the cellular, molecular and embryologic mechanisms underlying hyaloid morphogenesis can be elucidated.

PAX6 regulates melanogenesis in the retinal pigmented epithelium through feed-forward regulatory interactions with MITF

During organogenesis, PAX6 is required for establishment of various progenitor subtypes within the central nervous system, eye and pancreas. PAX6 expression is maintained in a variety of cell types within each organ, although its role in each lineage and how it acquires cell-specific activity remain elusive. Herein, we aimed to determine the roles and the hierarchical organization of the PAX6-dependent gene regulatory network during the differentiation of the retinal pigmented epithelium (RPE). Somatic mutagenesis of Pax6 in the differentiating RPE revealed that PAX6 functions in a feed-forward regulatory loop with MITF during onset of melanogenesis. PAX6 both controls the expression of an RPE isoform of Mitf and synergizes with MITF to activate expression of genes involved in pigment biogenesis. This study exemplifies how one kernel gene pivotal in organ formation accomplishes a lineage-specific role during terminal differentiation of a single lineage.

It is currently poorly understood how a single developmental transcription regulator controls early specification as well as a broad range of highly specialized differentiation schemes. PAX6 is one of the most extensively investigated factors in central nervous system development, yet its role in execution of lineage-specific programs remains mostly elusive. Here, we directly investigated the involvement of PAX6 in the differentiation of one lineage, the retinal pigmented epithelium (RPE), a neuroectodermal-derived tissue that is essential for retinal development and function. We revealed that PAX6 accomplishes its role through a unique regulatory interaction with the transcription factor MITF, a master regulator of the pigmentation program. During the differentiation of the RPE, PAX6 regulates the expression of an RPE-specific isoform of Mitf and importantly, at the same time, PAX6 functions together with MITF to directly activate the expression of downstream genes required for pigment biogenesis. These findings provide comprehensive insight into the gene hierarchy that controls RPE development: from a kernel gene (a term referring to the upper-most gene in the gene regulatory network) that is broadly expressed during CNS development through a lineage-specific transcription factor that together with the kernel gene creates cis-regulatory input that contributes to transcriptionally activate a battery of terminal differentiation genes.

Neuronal calcium-binding proteins 1/2 localize to dorsal root ganglia and excitatory spinal neurons and are regulated by nerve injury

Neuronal calcium (Ca(2+))-binding proteins 1 and 2 (NECAB1/2) are members of the phylogenetically conserved EF-hand Ca(2+)-binding protein superfamily. To date, NECABs have been explored only to a limited extent and, so far, not at all at the spinal level. Here, we describe the distribution, phenotype, and nerve injury-induced regulation of NECAB1/NECAB2 in mouse dorsal root ganglia (DRGs) and spinal cord. In DRGs, NECAB1/2 are expressed in around 70% of mainly small- and medium-sized neurons. Many colocalize with calcitonin gene-related peptide and isolectin B4, and thus represent nociceptors. NECAB1/2 neurons are much more abundant in DRGs than the Ca(2+)-binding proteins (parvalbumin, calbindin, calretinin, and secretagogin) studied to date. In the spinal cord, the NECAB1/2 distribution is mainly complementary. NECAB1 labels interneurons and a plexus of processes in superficial layers of the dorsal horn, commissural neurons in the intermediate area, and motor neurons in the ventral horn. Using CLARITY, a novel, bilaterally connected neuronal system with dendrites that embrace the dorsal columns like palisades is observed. NECAB2 is present in cell bodies and presynaptic boutons across the spinal cord. In the dorsal horn, most NECAB1/2 neurons are glutamatergic. Both NECAB1/2 are transported into dorsal roots and peripheral nerves. Peripheral nerve injury reduces NECAB2, but not NECAB1, expression in DRG neurons. Our study identifies NECAB1/2 as abundant Ca(2+)-binding proteins in pain-related DRG neurons and a variety of spinal systems, providing molecular markers for known and unknown neuron populations of mechanosensory and pain circuits in the spinal cord.

Gene regulation by PAX6: structural-functional correlations of missense mutants and transcriptional control of Trpm3/miR-204

PURPOSE

Pax6 is a key regulatory gene for eye, brain, and pancreas development. It acts as a transcriptional activator and repressor. Loss-of-function of Pax6 results in down- and upregulation of a comparable number of genes, although many are secondary targets. Recently, we found a prototype of a Pax6-binding site that acts as a transcriptional repressor. We also identified the Trpm3 gene as a Pax6-direct target containing the miR-204 gene located in intron 6. Thus, there are multiple Pax6-dependent mechanisms of transcriptional repression in the cell. More than 50 Pax6 missense mutations have been identified in humans and mice. Two of these mutations, N50K (Leca4) and R128C (Leca2), were analyzed in depth resulting in different numbers of regulated genes and different ratios of down- and upregulated targets. Thus, additional studies of these mutants are warranted to better understand the molecular mechanisms of the mutants' action.

METHODS

Mutations in PAX6 and PAX6(5a), including G18W, R26G, N50K, G64V, R128C, and R242T, were generated with site-directed mutagenesis. A panel of ten luciferase reporters driven by six copies of Pax6-binding sites representing a spectrum of sites that act as repressors, moderate activators, and strong activators were used. Two additional reporters, including the Pax6-regulated enhancer from mouse Trpm3 and six copies of its individual Pax6-binding site, were also tested in P19 cells.

RESULTS

PAX6 (N50K) acted either as a loss-of-function or neutral mutation. In contrast, PAX6 (R128C) and (R242T) acted as loss-, neutral, and gain-of-function mutations. With three distinct reporters, the PAX6 (N50K) mutation broke the pattern of effects produced by substitutions in the surrounding helices of the N-terminal region of the paired domain. All six mutations tested acted as loss-of-function using the Trpm3 Pax6-binding site.

CONCLUSIONS

These studies highlight the complexity of Pax6-dependent transcriptional activation and repression mechanisms, and identify the N50K and R128C substitutions as valuable tools for testing interactions between Pax6, Pax6 (N50K), and Pax6 (R128C) with other regulatory proteins, including chromatin remodelers.

Disruption of autoregulatory feedback by a mutation in a remote, ultraconserved PAX6 enhancer causes aniridia

The strictly regulated expression of most pleiotropic developmental control genes is critically dependent on the activity of long-range cis-regulatory elements. This was revealed by the identification of individuals with a genetic condition lacking coding-region mutations in the gene commonly associated with the disease but having a variety of nearby chromosomal abnormalities, collectively described as cis-ruption disease cases. The congenital eye malformation aniridia is caused by haploinsufficiency of the developmental regulator PAX6. We discovered a de novo point mutation in an ultraconserved cis-element located 150 kb downstream from PAX6 in an affected individual with intact coding region and chromosomal locus. The element SIMO acts as a strong enhancer in developing ocular structures. The mutation disrupts an autoregulatory PAX6 binding site, causing loss of enhancer activity, resulting in defective maintenance of PAX6 expression. These findings reveal a distinct regulatory mechanism for genetic disease by disruption of an autoregulatory feedback loop critical for maintenance of gene expression through development.

Histone posttranslational modifications and cell fate determination: lens induction requires the lysine acetyltransferases CBP and p300

Lens induction is a classical embryologic model to study cell fate determination. It has been proposed earlier that specific changes in core histone modifications accompany the process of cell fate specification and determination. The lysine acetyltransferases CBP and p300 function as principal enzymes that modify core histones to facilitate specific gene expression. Herein, we performed conditional inactivation of both CBP and p300 in the ectodermal cells that give rise to the lens placode. Inactivation of both CBP and p300 resulted in the dramatic discontinuation of all aspects of lens specification and organogenesis, resulting in aphakia. The CBP/p300^−/− ectodermal cells are viable and not prone to apoptosis. These cells showed reduced expression of Six3 and Sox2, while expression of Pax6 was not upregulated, indicating discontinuation of lens induction. Consequently, expression of αB- and αA-crystallins was not initiated. Mutant ectoderm exhibited markedly reduced levels of histone H3 K18 and K27 acetylation, subtly increased H3 K27me3 and unaltered overall levels of H3 K9ac and H3 K4me3. Our data demonstrate that CBP and p300 are required to establish lens cell-type identity during lens induction, and suggest that posttranslational histone modifications are integral to normal cell fate determination in the mammalian lens.

Functional dissection of the paired domain of Pax6 reveals molecular mechanisms of coordinating neurogenesis and proliferation

To achieve adequate organ development and size, cell proliferation and differentiation have to be tightly regulated and coordinated. The transcription factor Pax6 regulates patterning, neurogenesis and proliferation in forebrain development. The molecular basis of this regulation is not well understood. As the bipartite DNA-binding paired domain of Pax6 regulates forebrain development, we examined mice with point mutations in its individual DNA-binding subdomains PAI (Pax6(Leca4), N50K) and RED (Pax6(Leca2), R128C). This revealed distinct roles in regulating proliferation in the developing cerebral cortex, with the PAI and RED subdomain mutations reducing and increasing, respectively, the number of mitoses. Conversely, neurogenesis was affected only by the PAI subdomain mutation, phenocopying the neurogenic defects observed in full Pax6 mutants. Genome-wide expression profiling identified molecularly discrete signatures of Pax6(Leca4) and Pax6(Leca2) mutations. Comparison to Pax6 targets identified by chromatin immunoprecipitation led to the identification and functional characterization of distinct DNA motifs in the promoters of target genes dysregulated in the Pax6(Leca2) or Pax6(Leca4) mutants, further supporting the distinct regulatory functions of the DNA-binding subdomains. Thus, Pax6 achieves its key roles in the developing forebrain by utilizing particular subdomains to coordinate patterning, neurogenesis and proliferation simultaneously.

Stimulation of mouse Cyp1b1 during adipogenesis: characterization of promoter activation by the transcription factor Pax6

Cytochrome P4501B1 (Cyp1b1) is expressed specifically in certain neural crest (NC) cells during embryogenesis. Mesenchymal progenitor cells that develop from NC cells are modeled here by mouse C3H10T1/2 and 3T3-L1 cells. Dexamethasone in combination with methylisobutylxanthine (DM) induces Cyp1b1 and a 6.7 kb mouse Cyp1b1 promoter-luciferase reporter in each cell type prior to adipogenesis. An 18 base sequence (at -6.11 kb) (PaxE) which was essential for this reporter stimulation in 3T3-L1 cells bound the transcription factor Pax6. This is shown by gel mobility shifts and sequence mutations. Heterologous vector expression of Pax6 in 3T3-L1 cells enhanced DM stimulated Cyp1b1 promoter activity through cooperation with two Sp1 sites in the proximal promoter region. Chromatin immunoprecipitation showed that DM stimulated binding of Pax6 adjacent to Sp1 in the proximal promoter more than in the PaxE region. The Cyp1b1 induction by DM in C3H10T1/2 cells was more rapid but independent of Pax6. The far upstream enhancer region (FUER) found in rat Cyp1b1 responded to DM but was inactive in the mouse promoter due to key sequence changes. The expression patterns of Pax6 and Cyp1b1 frequently overlap during mouse embryogenesis. The relationship between Pax6 and Cyp1b1 expression warrants further investigation, particularly in the NC.

Pax6 regulates gene expression in the vertebrate lens through miR-204

During development, tissue-specific transcription factors regulate both protein-coding and non-coding genes to control differentiation. Recent studies have established a dual role for the transcription factor Pax6 as both an activator and repressor of gene expression in the eye, central nervous system, and pancreas. However, the molecular mechanism underlying the inhibitory activity of Pax6 is not fully understood. Here, we reveal that Trpm3 and the intronic microRNA gene miR-204 are co-regulated by Pax6 during eye development. miR-204 is probably the best known microRNA to function as a negative modulator of gene expression during eye development in vertebrates. Analysis of genes altered in mouse Pax6 mutants during lens development revealed significant over-representation of miR-204 targets among the genes up-regulated in the Pax6 mutant lens. A number of new targets of miR-204 were revealed, among them Sox11, a member of the SoxC family of pro-neuronal transcription factors, and an important regulator of eye development. Expression of Trpm/miR-204 and a few of its targets are also Pax6-dependent in medaka fish eyes. Collectively, this study identifies a novel evolutionarily conserved mechanism by which Pax6 controls the down-regulation of multiple genes through direct up-regulation of miR-204.

The transcription factor Pax6 is reiteratively employed in space and time for the establishment of progenitor pools and the differentiation of neuronal and non-neuronal lineages of the CNS, pancreas, and eye. Execution of these diverse developmental programs depends on simultaneous activation and repression of gene networks functioning downstream of Pax6. MicroRNAs function as inhibitors of gene expression. Many microRNA genes are transcribed through common promoters of host genes. In this study, using wide-scale analysis of changes in gene expression following Pax6 deletion in the lens, we discover that Pax6 regulates the gene Trpm3 and its hosted microRNA, miR-204. We then show that miR-204 suppresses several target genes in the lens, notably the neuronal gene Sox11. Lastly, by conducting parallel experiments in the medaka fish, we show that Pax6 control of miR-204 and its target genes is evolutionarily conserved between mammals and fish, stressing the biological importance of this pathway. Pax6 regulation of miR-204 explains part of the complex, divergent inhibitory activity of Pax6 in ocular progenitor cells, which is required to establish and maintain the identity and function of ocular tissues.

Pax6 interactions with chromatin and identification of its novel direct target genes in lens and forebrain

Pax6 encodes a specific DNA-binding transcription factor that regulates the development of multiple organs, including the eye, brain and pancreas. Previous studies have shown that Pax6 regulates the entire process of ocular lens development. In the developing forebrain, Pax6 is expressed in ventricular zone precursor cells and in specific populations of neurons; absence of Pax6 results in disrupted cell proliferation and cell fate specification in telencephalon. In the pancreas, Pax6 is essential for the differentiation of α-, β- and δ-islet cells. To elucidate molecular roles of Pax6, chromatin immunoprecipitation experiments combined with high-density oligonucleotide array hybridizations (ChIP-chip) were performed using three distinct sources of chromatin (lens, forebrain and β-cells). ChIP-chip studies, performed as biological triplicates, identified a total of 5,260 promoters occupied by Pax6. 1,001 (133) of these promoter regions were shared between at least two (three) distinct chromatin sources, respectively. In lens chromatin, 2,335 promoters were bound by Pax6. RNA expression profiling from Pax6^+/− lenses combined with in vivo Pax6-binding data yielded 76 putative Pax6-direct targets, including the Gaa, Isl1, Kif1b, Mtmr2, Pcsk1n, and Snca genes. RNA and ChIP data were validated for all these genes. In lens cells, reporter assays established Kib1b and Snca as Pax6 activated and repressed genes, respectively. In situ hybridization revealed reduced expression of these genes in E14 cerebral cortex. Moreover, we examined differentially expressed transcripts between E9.5 wild type and Pax6^−/− lens placodes that suggested Efnb2, Fat4, Has2, Nav1, and Trpm3 as novel Pax6-direct targets. Collectively, the present studies, through the identification of Pax6-direct target genes, provide novel insights into the molecular mechanisms of Pax6 gene control during mouse embryonic development. In addition, the present data demonstrate that Pax6 interacts preferentially with promoter regions in a tissue-specific fashion. Nevertheless, nearly 20% of the regions identified are accessible to Pax6 in multiple tissues.

Origins of shared genetic variation in African cichlids

Cichlid fishes have evolved tremendous morphological and behavioral diversity in the waters of East Africa. Within each of the Great Lakes Tanganyika, Malawi, and Victoria, the phenomena of hybridization and retention of ancestral polymorphism explain allele sharing across species. Here, we explore the sharing of single nucleotide polymorphisms (SNPs) between the major East African cichlid assemblages. A set of approximately 200 genic and nongenic SNPs was ascertained in five Lake Malawi species and genotyped in a diverse collection of ∼160 species from across Africa. We observed segregating polymorphism outside of the Malawi lineage for more than 50% of these loci; this holds similarly for genic versus nongenic SNPs, as well as for SNPs at putative CpG versus non-CpG sites. Bayesian and principal component analyses of genetic structure in the data demonstrate that the Lake Malawi endemic flock is not monophyletic and that river species have likely contributed significantly to Malawi genomes. Coalescent simulations support the hypothesis that river cichlids have transported polymorphism between lake assemblages. We observed strong genetic differentiation between Malawi lineages for approximately 8% of loci, with contributions from both genic and nongenic SNPs. Notably, more than half of these outlier loci between Malawi groups are polymorphic outside of the lake. Cichlid fishes have evolved diversity in Lake Malawi as new mutations combined with standing genetic variation shared across East Africa.

Development of embryonic stem cells in recombinant kidneys

Embryonic stem cells (ESC) are self-renewing and can generate all cell types during normal development. Previous studies have begun to explore fates of ESCs and their mesodermal derivatives after injection into explanted intact metanephric kidneys and neonatal kidneys maturing in vivo. Here, we exploited a recently described recombinant organ culture model, mixing fluorescent quantum dot labeled mouse exogenous cells with host metanephric cells. We compared abilities of undifferentiated ESCs with ESC-derived mesodermal or non-mesodermal cells to contribute to tissue compartments within recombinant, chimeric metanephroi. ESC-derived mesodermal cells downregulated Oct4, a marker of undifferentiated cells, and, as assessed by locations of quantum dots, contributed to Wilms’ tumor 1-expressing forming nephrons, synaptopodin-expressing glomeruli, and organic ion-transporting tubular epithelia. Similar results were observed when labeled native metanephric cells were recombined with host cells. In striking contrast, non-mesodermal ESC-derived cells strongly inhibited growth of embryonic kidneys, while undifferentiated ESCs predominantly formed Oct4 expressing colonies between forming nephrons and glomeruli. These findings clarify the conclusion that ESC-derived mesodermal cells have functional nephrogenic potential, supporting the idea that they could potentially replace damaged epithelia in diseased kidneys. On the other hand, undifferentiated ESCs and non-mesodermal precursors derived from ESCs would appear to be less suitable materials for use in kidney cell therapies.

Patterns of gene expression in microarrays and expressed sequence tags from normal and cataractous lenses

In this contribution, we have examined the patterns of gene expression in normal and cataractous lenses as presented in five different papers using microarrays and expressed sequence tags. The purpose was to evaluate unique and common patterns of gene expression during development, aging and cataracts.

Analysis of the transcriptomes downstream of Eyeless and the Hedgehog, Decapentaplegic and Notch signaling pathways in Drosophila melanogaster

Tissue-specific transcription factors are thought to cooperate with signaling pathways to promote patterned tissue specification, in part by co-regulating transcription. The Drosophila melanogaster Pax6 homolog Eyeless forms a complex, incompletely understood regulatory network with the Hedgehog, Decapentaplegic and Notch signaling pathways to control eye-specific gene expression. We report a combinatorial approach, including mRNAseq and microarray analyses, to identify targets co-regulated by Eyeless and Hedgehog, Decapentaplegic or Notch. Multiple analyses suggest that the transcriptomes resulting from co-misexpression of Eyeless+signaling factors provide a more complete picture of eye development compared to previous efforts involving Eyeless alone: (1) Principal components analysis and two-way hierarchical clustering revealed that the Eyeless+signaling factor transcriptomes are closer to the eye control transcriptome than when Eyeless is misexpressed alone; (2) more genes are upregulated at least three-fold in response to Eyeless+signaling factors compared to Eyeless alone; (3) based on gene ontology analysis, the genes upregulated in response to Eyeless+signaling factors had a greater diversity of functions compared to Eyeless alone. Through a secondary screen that utilized RNA interference, we show that the predicted gene CG4721 has a role in eye development. CG4721 encodes a neprilysin family metalloprotease that is highly up-regulated in response to Eyeless+Notch, confirming the validity of our approach. Given the similarity between D. melanogaster and vertebrate eye development, the large number of novel genes identified as potential targets of Ey+signaling factors will provide novel insights to our understanding of eye development in D. melanogaster and humans.

Pax6: a multi-level regulator of ocular development

Eye development has been a paradigm for the study of organogenesis, from the demonstration of lens induction through epithelial tissue morphogenesis, to neuronal specification and differentiation. The transcription factor Pax6 has been shown to play a key role in each of these processes. Pax6 is required for initiation of developmental pathways, patterning of epithelial tissues, activation of tissue-specific genes and interaction with other regulatory pathways. Herein we examine the data accumulated over the last few decades from extensive analyses of biochemical modules and genetic manipulation of the Pax6 gene. Specifically, we describe the regulation of Pax6's expression pattern, the protein's DNA-binding properties, and its specific roles and mechanisms of action at all stages of lens and retinal development. Pax6 functions at multiple levels to integrate extracellular information and execute cell-intrinsic differentiation programs that culminate in the specification and differentiation of a distinct ocular lineage.

Mab21l2 is essential for embryonic heart and liver development

During mouse embryogenesis, proper formation of the heart and liver is especially important and is crucial for embryonic viability. In this study, we showed that Mab21l2 was expressed in the trabecular and compact myocardium, and that deletion of Mab21l2 resulted in a reduction of the trabecular myocardium and thinning of the compact myocardium. Mab21l2-deficient embryonic hearts had decreased expression of genes that regulate cell proliferation and apoptosis of cardiomyocytes. These results show that Mab21l2 functions during heart development by regulating the expression of such genes. Mab21l2 was also expressed in the septum transversum mesenchyme (STM). Epicardial progenitor cells are localized to the anterior surface of the STM (proepicardium), and proepicardial cells migrate onto the surface of the heart and form the epicardium, which plays an important role in heart development. The rest of the STM is essential for the growth and survival of hepatoblasts, which are bipotential progenitors for hepatocytes and cholangiocytes. Proepicardial cells in Mab21l2-deficient embryos had defects in cell proliferation, which led to a small proepicardium, in which α4 integrin expression, which is essential for the migration of proepicardial cells, was down-regulated, suggesting that defects occurred in its migration. In Mab21l2-deficient embryos, epicardial formation was defective, suggesting that Mab21l2 plays important roles in epicardial formation through the regulation of the cell proliferation of proepicardial cells and the migratory process of proepicardial cells. Mab21l2-deficient embryos also exhibited hypoplasia of the STM surrounding hepatoblasts and decreased hepatoblast proliferation with a resultant loss of defective morphogenesis of the liver. These findings demonstrate that Mab21l2 plays a crucial role in both heart and liver development through STM formation.