Deletion upstream of MAB21L2 highlights the importance of evolutionarily conserved non-coding sequences for eye development.

Anophthalmia, microphthalmia and coloboma (AMC) comprise a spectrum of developmental eye disorders, accounting for approximately 20% of childhood visual impairment. While non-coding regulatory sequences are increasingly recognised as contributing to disease burden, characterising their impact on gene function and phenotype remains challenging. Furthermore, little is known of the nature and extent of their contribution to AMC phenotypes. We report two families with variants in or near MAB21L2, a gene where genetic variants are known to cause AMC in humans and animal models. The first proband, presenting with microphthalmia and coloboma, has a likely pathogenic missense variant (c.338 G > C; p.[Trp113Ser]), segregating within the family. The second individual, presenting with microphthalmia, carries an ~ 113.5 kb homozygous deletion 19.38 kb upstream of MAB21L2. Modelling of the deletion results in transient small lens and coloboma as well as midbrain anomalies in zebrafish, and microphthalmia and coloboma in Xenopus tropicalis. Using conservation analysis, we identify 15 non-coding conserved elements (CEs) within the deleted region, while ChIP-seq data from mouse embryonic stem cells demonstrates that two of these (CE13 and 14) bind Otx2, a protein with an established role in eye development. Targeted disruption of CE14 in Xenopus tropicalis recapitulates an ocular coloboma phenotype, supporting its role in eye development. Together, our data provides insights into regulatory mechanisms underlying eye development and highlights the importance of non-coding sequences as a source of genetic diagnoses in AMC.

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.

CRISPR-Cas9-mediated functional dissection of the foxc1 genomic region in zebrafish identifies critical conserved cis-regulatory elements.

FOXC1 encodes a forkhead-domain transcription factor associated with several ocular disorders. Correct FOXC1 dosage is critical to normal development, yet the mechanisms controlling its expression remain unknown. Together with FOXQ1 and FOXF2, FOXC1 is part of a cluster of FOX genes conserved in vertebrates. CRISPR-Cas9-mediated dissection of genomic sequences surrounding two zebrafish orthologs of FOXC1 was performed. This included five zebrafish-human conserved regions, three downstream of foxc1a and two remotely upstream of foxf2a/foxc1a or foxf2b/foxc1b clusters, as well as two intergenic regions between foxc1a/b and foxf2a/b lacking sequence conservation but positionally corresponding to the area encompassing a previously reported glaucoma-associated SNP in humans. Removal of downstream sequences altered foxc1a expression; moreover, zebrafish carrying deletions of two or three downstream elements demonstrated abnormal phenotypes including enlargement of the anterior chamber of the eye reminiscent of human congenital glaucoma. Deletions of distant upstream conserved elements influenced the expression of foxf2a/b or foxq1a/b but not foxc1a/b within each cluster. Removal of either intergenic sequence reduced foxc1a or foxc1b expression during late development, suggesting a role in transcriptional regulation despite the lack of conservation at the nucleotide level. Further studies of the identified regions in human patients may explain additional individuals with developmental ocular disorders.

Genetic disruption of zebrafish mab21l1 reveals a conserved role in eye development and affected pathways.

BACKGROUND: The male-abnormal 21 like (MAB21L) genes are important in human ocular development. Homozygous loss of MAB21L1 leads to corneal dystrophy in all affected individuals along with cataracts and buphthalmos in some. The molecular function and downstream pathways of MAB21L factors are largely undefined. RESULTS: We generated the first mab21l1 zebrafish mutant carrying a putative loss-of-function allele, c.107delA p.(Lys36Argfs*7). At the final stages of embryonic development, homozygous mab21l1c.107delA fish displayed enlarged anterior chambers and corneal thinning which progressed with age. Additional studies revealed increased cell death in the mutant corneas, transformation of the cornea into a skin-like epithelium, and progressive lens degeneration with development of fibrous masses in the anterior chamber. RNA-seq of wild-type and mutant ocular transcriptomes revealed significant changes in expression of several genes, including irf1a and b, stat1, elf3, krt17, tlr9, and loxa associated with immunity and/or corneal function. Abnormal expression of lysyl oxidases have been previously linked with corneal thinning, fibrosis, and lens defects in mammals, suggesting a role for loxa misexpression in the progressive mab21l1c.107delA eye phenotype. CONCLUSIONS: Zebrafish mab21l1 is essential for normal corneal development, similar to human MAB21L1. The identified molecular changes in mab21l1c.107delA mutants provide the first clues about possible affected pathways.

De Novo Missense Variants in WDR37 Cause a Severe Multisystemic Syndrome.

While genetic causes are known for many syndromes involving developmental anomalies, a large number of individuals with overlapping phenotypes remain undiagnosed. Using exome-sequencing analysis and review of matchmaker databases, we have discovered four de novo missense variants predicted to affect the N-terminal region of WDR37-p.Ser119Phe, p.Thr125Ile, p.Ser129Cys, and p.Thr130Ile-in unrelated individuals with a previously unrecognized syndrome. Features of WDR37 syndrome include the following: ocular anomalies such as corneal opacity/Peters anomaly, coloboma, and microcornea; dysmorphic facial features; significant neurological impairment with structural brain defects and seizures; poor feeding; poor post-natal growth; variable skeletal, cardiac, and genitourinary defects; and death in infancy in one individual. WDR37 encodes a protein of unknown function with seven predicted WD40 domains and no previously reported human pathogenic variants. Immunocytochemistry and western blot studies showed that wild-type WDR37 is localized predominantly in the cytoplasm and mutant proteins demonstrate similar protein levels and localization. CRISPR-Cas9-mediated genome editing generated zebrafish mutants with novel missense and frameshift alleles: p.Ser129Phe, p.Ser129Cys (which replicates one of the human variants), p.Ser129Tyr, p.Lys127Cysfs, and p.Gln95Argfs. Zebrafish carrying heterozygous missense variants demonstrated poor growth and larval lethality, while heterozygotes with frameshift alleles survived to adulthood, suggesting a potential dominant-negative mechanism for the missense variants. RNA-seq analysis of zebrafish embryos carrying a missense variant detected significant upregulation of cholesterol biosynthesis pathways. This study identifies variants in WDR37 associated with human disease and provides insight into its essential role in vertebrate development and possible molecular functions.

PITX2 deficiency and associated human disease: insights from the zebrafish model.

The PITX2 (paired-like homeodomain 2) gene encodes a bicoid-like homeodomain transcription factor linked with several human disorders. The main associated congenital phenotype is Axenfeld-Rieger syndrome, type 1, an autosomal dominant condition characterized by variable defects in the anterior segment of the eye, an increased risk of glaucoma, craniofacial dysmorphism and dental and umbilical anomalies; in addition to this, one report implicated PITX2 in ring dermoid of the cornea and a few others described cardiac phenotypes. We report three novel PITX2 mutations-c.271C > T, p.(Arg91Trp); c.259T > C, p.(Phe87Leu); and c.356delA, p.(Gln119Argfs*36)-identified in independent families with typical Axenfeld-Rieger syndrome characteristics and some unusual features such as corneal guttata, Wolf-Parkinson-White syndrome, and hyperextensibility. To gain further insight into the diverse roles of PITX2/pitx2 in vertebrate development, we generated various genetic lesions in the pitx2 gene via TALEN-mediated genome editing. Affected homozygous zebrafish demonstrated congenital defects consistent with the range of PITX2-associated human phenotypes: abnormal development of the cornea, iris and iridocorneal angle; corneal dermoids; and craniofacial dysmorphism. In addition, via comparison of pitx2M64* and wild-type embryonic ocular transcriptomes we defined molecular changes associated with pitx2 deficiency, thereby implicating processes potentially underlying disease pathology. This analysis identified numerous affected factors including several members of the Wnt pathway and collagen types I and V gene families. These data further support the link between PITX2 and the WNT pathway and suggest a new role in regulation of collagen gene expression during development.

Functional characterization of zebrafish orthologs of the human Beta 3-Glucosyltransferase B3GLCT gene mutated in Peters Plus Syndrome.

Peters Plus Syndrome (PPS) is a rare autosomal recessive disease characterized by ocular defects, short stature, brachydactyly, characteristic facial features, developmental delay and other highly variable systemic defects. Classic PPS is caused by loss-of-function mutations in the B3GLCT gene encoding for a ß3-glucosyltransferase that catalyzes the attachment of glucose via a ß1-3 glycosidic linkage to O-linked fucose on thrombospondin type 1 repeats (TSRs). B3GLCT was shown to participate in a non-canonical ER quality control mechanism; however, the exact molecular processes affected in PPS are not well understood. Here we report the identification and characterization of two zebrafish orthologs of the human B3GLCT gene, b3glcta and b3glctb. The b3glcta and b3glctb genes encode for 496-aa and 493-aa proteins with 65% and 57% identity to human B3GLCT, respectively. Expression studies demonstrate that both orthologs are widely expressed with strong presence in embryonic tissues affected in PPS. In vitro glucosylation assays demonstrated that extracts from wildtype embryos contain active b3glct enzyme capable of transferring glucose from UDP-glucose to an O-fucosylated TSR, indicating functional conservation with human B3GLCT. To determine the developmental role of the zebrafish genes, single and double b3glct knockouts were generated using TALEN-induced genome editing. Extracts from double homozygous b3glct-/- embryos demonstrated complete loss of in vitro b3glct activity. Surprisingly, b3glct-/- homozygous fish developed normally. Transcriptome analyses of head and trunk tissues of b3glct-/- 24-hpf embryos identified 483 shared differentially regulated transcripts that may be involved in compensation for b3glct function in these embryos. The presented data show that both sequence and function of B3GLCT/b3glct genes is conserved in vertebrates. At the same time, complete b3glct deficiency in zebrafish appears to be inconsequential and possibly compensated for by a yet unknown mechanism.

EFTUD2 deficiency in vertebrates: Identification of a novel human mutation and generation of a zebrafish model.

BACKGROUND: Congenital microphthalmia and coloboma are severe developmental defects that are frequently associated with additional systemic anomalies and display a high level of genetic heterogeneity. METHODS: To identify the pathogenic variant in a patient with microphthalmia, coloboma, retinal dystrophy, microcephaly, and other features, whole exome sequencing analysis of the patient and parental samples was undertaken. To further explore the identified variant/gene, expression and functional studies in zebrafish were performed. RESULTS: Whole exome sequencing revealed a de novo variant, c.473_474delGA, p.(Arg158Lysfs*4), in EFTUD2 which encodes a component of the spliceosome complex. Dominant mutations in EFTUD2 cause Mandibulofacial Dysostosis, Guion-Almeida type, which does not involve microphthalmia, coloboma, or retinal dystrophy; analysis of genes known to cause these ocular phenotypes identified several variants of unknown significance but no causal alleles in the affected patient. Zebrafish eftud2 demonstrated high sequence conservation with the human gene and broad embryonic expression. TALEN-mediated disruption was employed to generate a c.378_385 del, p.(Ser127Aspfs*23) truncation mutation in eftud2. Homozygous mutants displayed a reduced head size, small eye, curved body, and early embryonic lethality. Apoptosis assays demonstrated a striking increase in terminal deoxynucleotidyl transferase-mediated deoxyuridine triphosphate nick end-labeling (TUNEL)-positive cells in the developing brain, eye, spinal cord, and other tissues starting at 30 hours postfertilization. CONCLUSION: This study reports a novel mutation in EFTUD2 in a Mandibulofacial Dysostosis, Guion-Almeida type patient with unusual ocular features and the generation of a first animal model of eftud2 deficiency. The severe embryonic phenotype observed in eftud2 mutants indicates an important conserved role during development of diverse tissues in vertebrates.

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.

Whole exome sequencing in dominant cataract identifies a new causative factor, CRYBA2, and a variety of novel alleles in known genes.

Pediatric cataracts are observed in 1-15 per 10,000 births with 10-25 % of cases attributed to genetic causes; autosomal dominant inheritance is the most commonly observed pattern. Since the specific cataract phenotype is not sufficient to predict which gene is mutated, whole exome sequencing (WES) was utilized to concurrently screen all known cataract genes and to examine novel candidate factors for a disease-causing mutation in probands from 23 pedigrees affected with familial dominant cataract. Review of WES data for 36 known cataract genes identified causative mutations in nine pedigrees (39 %) in CRYAA, CRYBB1, CRYBB3, CRYGC (2), CRYGD, GJA8 (2), and MIP and an additional likely causative mutation in EYA1; the CRYBB3 mutation represents the first dominant allele in this gene and demonstrates incomplete penetrance. Examination of crystallin genes not yet linked to human disease identified a novel cataract gene, CRYBA2, a member of the ßγ-crystallin superfamily. The p.(Val50Met) mutation in CRYBA2 cosegregated with disease phenotype in a four-generation pedigree with autosomal dominant congenital cataracts with incomplete penetrance. Expression studies detected cryba2 transcripts during early lens development in zebrafish, supporting its role in congenital disease. Our data highlight the extreme genetic heterogeneity of dominant cataract as the eleven causative/likely causative mutations affected nine different genes, and the majority of mutant alleles were novel. Furthermore, these data suggest that less than half of dominant cataract can be explained by mutations in currently known genes.

MIP/Aquaporin 0 represents a direct transcriptional target of PITX3 in the developing lens.

The PITX3 bicoid-type homeodomain transcription factor plays an important role in lens development in vertebrates. PITX3 deficiency results in a spectrum of phenotypes from isolated cataracts to microphthalmia in humans, and lens degeneration in mice and zebrafish. While identification of downstream targets of PITX3 is vital for understanding the mechanisms of normal ocular development and human disease, these targets remain largely unknown. To isolate genes that are directly regulated by PITX3, we performed a search for genomic sequences that contain evolutionarily conserved bicoid/PITX3 binding sites and are located in the proximity of known genes. Two bicoid sites that are conserved from zebrafish to human were identified within the human promoter of the major intrinsic protein of lens fiber, MIP/AQP0. MIP/AQP0 deficiency was previously shown to be associated with lens defects in humans and mice. We demonstrate by both chromatin immunoprecipitation and electrophoretic mobility shift assay that PITX3 binds to MIP/AQP0 promoter region in vivo and is able to interact with both bicoid sites in vitro. In addition, we show that wild-type PITX3 is able to activate the MIP/AQP0 promoter via interaction with the proximal bicoid site in cotransfection experiments and that the introduction of mutations disrupting binding to this site abolishes this activation. Furthermore, mutant forms of PITX3 fail to produce the same levels of transactivation as wild-type when cotransfected with the MIP/AQP0 reporter. Finally, knockdown of pitx3 in zebrafish affects formation of a DNA-protein complex associated with mip1 promoter sequences; and examination of expression in pitx3 morphant and control zebrafish revealed a delay in and reduction of mip1 expression in pitx3-deficient embryos. Therefore, our data suggest that PITX3 is involved in direct regulation of MIP/AQP0 expression and that the alteration of MIP/AQP0 expression is likely to contribute to the lens phenotype in cataract patients with PITX3 mutations.

Three-color cDNA microarrays with prehybridization quality control yield gene expression data comparable to that of commercial platforms.

Despite their lower cost and high content flexibility, a limitation of in-house-prepared arrays has been their susceptibility to quality control (QC) issues and lack of QC standards across laboratories. Therefore, we developed a novel three-color array system that allows prehybridization QC as well as the Matarray software to facilitate acquisition of accurate gene expression data. In this study, we compared performance of our rat cDNA array to the Affymetrix RG-U34A and Agilent G4130A arrays using 2,824 UniGenes represented on all three arrays. Before data filtering, poor interplatform agreement was observed; however, after data filtering, differentially expressed UniGenes exhibited correlation coefficients of 0.91, 0.88, and 0.92 between the Affymetrix vs. Agilent, Affymetrix vs. cDNA, and Agilent vs. cDNA arrays, respectively. The Affymetrix, Agilent, and cDNA arrays agreed well with quantitative RT-PCR conducted on 42 UniGenes, yielding correlation coefficients of 0.90, 0.90, and 0.96, respectively. Each platform underestimated ratios relative to quantitative RT-PCR, possessing respective slopes of 0.86 (R2 = 0.81), 0.65 (R2 = 0.81), and 0.70 (R2 = 0.92). Overall, these data show that the combination of our novel technical and analytic approaches yield an accurate platform for functional genomics that is concordant with commercial discovery arrays in terms of identifying regulated genes and pathways.

Immobilized probe and glass surface chemistry as variables in microarray fabrication.

BACKGROUND: Global gene expression studies with microarrays can offer biological insights never before possible. However, the technology possesses many sources of technical variability that are an obstacle to obtaining high quality data sets. Since spotted microarrays offer design/content flexibility and potential cost savings over commercial systems, we have developed prehybridization quality control strategies for spotted cDNA and oligonucleotide arrays. These approaches utilize a third fluorescent dye (fluorescein) to monitor key fabrication variables, such as print/spot morphology, DNA retention, and background arising from probe redistributed during blocking. Here, our labeled cDNA array platform is used to study, 1) compression of array data using known input ratios of Arabidopsis in vitro transcripts and arrayed serial dilutions of homologous probes; 2) how curing time of in-house poly-L-lysine coated slides impacts probe retention capacity; and 3) the retention characteristics of 13 commercially available surfaces. RESULTS: When array element fluorescein intensity drops below 5,000 RFU/pixel, gene expression measurements become increasingly compressed, thereby validating this value as a prehybridization quality control threshold. We observe that the DNA retention capacity of in-house poly-L-lysine slides decreases rapidly over time (~50% reduction between 3 and 12 weeks post-coating; p < 0.0002) and that there are considerable differences in retention characteristics among commercially available poly-L-lysine and amino silane-coated slides. CONCLUSIONS: High DNA retention rates are necessary for accurate gene expression measurements. Therefore, an understanding of the characteristics and optimization of protocols to an array surface are prerequisites to fabrication of high quality arrays.