|Home | About | Journals | Submit | Contact Us | Français|
Users may view, print, copy, download and text and data- mine the content in such documents, for the purposes of academic research, subject always to the full Conditions of use: http://www.nature.com/authors/editorial_policies/license.html#terms
Individuals with nonsyndromic congenital retinal nonattachment (NCRNA) are totally blind from birth. The disease afflicts ~1% of Kurdish people living in a group of neighboring villages in North Khorasan, Iran. We show NCRNA is caused by a 6523bp deletion that spans a remote cis regulatory element 20 kb upstream from ATOH7 (Math5), a bHLH transcription factor gene required for retinal ganglion cell (RGC) and optic nerve development. In humans, the absence of RGCs stimulates massive neovascular growth of fetal blood vessels within the vitreous, and early retinal detachment. The remote ATOH7 element appears to act as a secondary or ‘shadow’ transcriptional enhancer. It has minimal sequence similarity to the primary enhancer, which is close to the Atoh7 promoter, but drives transgene expression with an identical spatiotemporal pattern in the mouse retina. The human transgene also functions in zebrafish, reflecting deep evolutionary conservation. These dual enhancers may reinforce Atoh7 expression during early critical stages of eye development when retinal neurogenesis is initiated.
The vertebrate neural retina develops during embryogenesis from multipotent progenitors in the optic cup. These cells exit mitosis and differentiate to form photoreceptors, interneurons, Müller glia, and retinal ganglion cells (RGCs) through the actions of intrinsic and extrinsic factors1. Seven major cell types emerge from the progenitor pool in a relatively fixed order, with overlapping birthdates, and differentiate to form a highly laminar structure that detects and processes visual information. Histogenesis starts in the central retina and spreads to the periphery over several days to weeks, depending on the species. Ganglion cells are the first-born retinal neurons in all vertebrates and occupy the innermost nuclear layer. Their axons comprise the optic nerves, which convey all visual information from the eyes to the brain. In mice, RGC genesis begins on embryonic day E11 and peaks around E13. In humans, RGCs are born between the 5th and 11th weeks of gestation. In addition to neurons, the retina contains astrocytes, which originate from a separate progenitor pool in the optic stalk. These enter the eye during late gestation, migrate across the inner retinal surface, and retain mitotic potential2.
Intrinsic factors controlling retinal histogenesis have been defined in large part through analysis of mouse genetic models. An archetypal example is the proneural bHLH (basic helix-loop-helix) transcription factor Atoh7 (Math5), a homolog of Drosophila Atonal3. The Atoh7 gene is transiently expressed in newly post-mitotic precursors, beginning with the onset of neurogenesis4. Atoh7 mutant mice lack optic nerves and have a >95% reduction in ganglion cell abundance5,6. No other cell type is missing. Their eyes detect light and process visual information within the retina but transmit few, if any, signals to the brain7. An orthologous ath5 mutation was identified in the lakritz zebrafish8. In mouse Cre-lox lineage experiments, Atoh7+ precursors give rise to 5% of cells in the normal adult retina, but only a fraction of these precursors develop as RGCs4,9. Atoh7 thus acts as a competence factor for RGC fate. Differentially expressed mRNAs have been identified in mutant retinas10, but the exact mechanism of Atoh7 action is incompletely understood.
The mammalian retina is nourished by two vascular systems. The inner two-thirds is supplied by the central artery, which enters the eye through the optic nerve and branches to form surface and deep capillary beds. The outer retina is oxygenated by diffusion from the surrounding choriocapillaris. Both systems derive from the fetal circulation11,12. In the embryo, inner hyaloid and outer ciliary vessels envelop the optic cup, anastomose across the rim, and fill the vitreous space. Starting at P0 in mice, vessels surrounding the lens regress through a branch-wise process of apoptosis, and the hyaloid trunk is remodeled to form the central retinal artery. In humans, these events occur during late gestation and early postnatal life13.
Many forms of hereditary blindness are caused by mutations in genes required for retinal cell maintenance or physiology, leading to slow or rapidly progressing degeneration14. A smaller number of diseases involve genes such as CRX or NR2E3, which are required for photoreceptor differentiation15. As yet, no human retinal disorder has been specifically associated with defects in fate determination of non-photoreceptor cell types. Congenital retinal vascular disorders are also major causes of blindness worldwide. In retinopathy of prematurity (ROP), familial exudative vitreoretinopathy (FEVR) and persistent hyperplasia of the primary vitreous (PHPV), ischemia and angiogenic signaling defects are the primary stimulus for pathological retinal neovascularization11.
Congenital retinal nonattachment is a severe form of blindness, with no detection of light from birth16. The retina is partially or completely separated from the posterior globe, and there is typically a fibrovascular mass behind the lens. The disease is heterogeneous, with a variety of clinical and histological anomalies involving the retina and blood vessels in one or both eyes. In some patients, the disease occurs with systemic malformations, but there is no consistent association. Sporadic and familial cases have been reported17. A nonsyndromic hereditary form (NCRNA) is common in North Khorasan, Iran, among descendants of a Kurdish founder population, and it is inherited as an autosomal recessive trait with complete penetrance18. The locus was mapped by SSLP linkage and haplotype analysis to a 0.6–1.5 cM region on chromosome 10q21, but the causative mutation and pathogenic mechanism are unknown19,20. In principle, the disease could arise from an intrinsic flaw in neuronal or vascular development.
In this report, we show that blind NCRNA patients lack optic nerves, which prevents light perception and profoundly affects the development of retinal vasculature. Our data suggest that NCRNA is caused by deletion of conserved DNA sequences that appear to control ATOH7 expression in early neuroretinal precursors. Finally, we show that the deleted sequences have properties consistent with a duplicate or ‘shadow’ enhancer21, including:  an identical pattern of transcriptional activation as the primary Atoh7 enhancer,  a remote genomic location,  a comparatively recent evolutionary origin and  conserved biological activity among vertebrates. NCRNA is thus a Mendelian is order that disrupts the intrinsic mechanism of retinal cell fate determination.
This autosomal recessive disorder (MIM221900) was first described in a large consanguineous Kurdish family that traces its ancestry for nine generations (Fig. 1g). The founders were dislocated centuries ago from Kurdistan to North Khorasan province in eastern Iran, during the Safavid Dynasty22, and their descendants remained relatively isolated, linguistically, culturally and genetically, from the main Persian population until recent years. The disease afflicts 1.1% of the population in four neighboring villages, with an estimated carrier frequency of 18% for the NCRNA allele18.
The affected individuals are totally blind, with no perception of light. Their eyes exhibit leukocoria (white pupils), esotropia (misalignment with inward deviation), and a continuous pendular nystagmus (jerky eye movements) (Fig. 1a–f). The pupils have no direct or consensual light reaction. In three-quarters of eyes examined, the anterior chambers were relatively shallow, but there was no consistent elevation of intraocular pressure. During a previous study18, B-scan ultrasonography and serial fundoscopic examinations revealed formation of a dense mass behind the lens in each eye, with complete opacification of the vitreous by four months (Fig. 1e,f). There is also a slowly progressive corneal clouding, consistent with endothelial blood staining23, which was evident to some extent in all adults examined. This develops over several years, is more severe in the inferior part of the cornea (Fig. 1b), and may occur with iridocorneal adhesions. No other sensory or neurological defects have been detected in the blind individuals. Their auditory, olfactory, endocrine and intellectual functions are intact, and they have normal developmental milestones. They have no systemic or other vascular abnormalities. Together, these clinical findings suggest bilateral persistence of the fetal blood vessels, with complete congenital detachment of the neural retina, and intravitreal vascular proliferation during infancy. The obligate heterozygotes are clinically normal. They have good vision, with no structural or functional eye pathology.
Given the complete blindness, congenital nystagmus, and absent pupillary light reflexes, we investigated optic nerve defects as the primary anatomical basis for NCRNA pathology. Nystagmus is common in children with severe optic nerve malformations24, and it is thought to arise from unbalanced sensory input to central autonomic centers controlling gaze, and a failure to calibrate oculomotor systems. The pupillary light reflex is triggered by rare photosensitive RGCs, which relay signals to the cortex of the olivary pretectal nuclei via the optic nerves. Proximity of the NCRNA region to ATOH7 also raised the possibility of optic nerve involvement. Atoh7 mutant mice have essentially no pupillary reflex25, and significant retinovascular abnormalities26, although retinal detachments are rare.
In the original fundoscopic exams of one affected infant prior to vitreal opacification, the status of the optic discs could not be determined18. To characterize the optic nerves anatomically, we performed orbital MRI exams on one adult NCRNA patient (Fig. 2a–c,f). As expected, this study showed heterogeneous consolidation of the vitreous compartments, consistent with complete retinal detachment and a retrolental fibrovascular mass in each eye. The sclerae (white fibrous outer layer of the eye globes), lens, extraocular muscles, retrobulbar fat and adnexa (lacrimal glands and other accessory visual structures) appeared relatively normal, and there were no anatomical defects in the anterior brain. In contrast, the optic nerves and chiasm were absent or greatly attenuated in multiple views. However, because the oculomotor and ciliary nerves and ophthalmic blood vessels were difficult to resolve within the orbital cones, we cannot exclude the presence of small atrophic optic nerves in these patients, similar to residual RGCs reported in some Atoh7 mutant mice27.
The critical region spans 1639kb, from D10S1670 to D10S1418, and overlaps 14 positional candidate genes, including ATOH7 (Fig. 3a). We therefore tested blind (rr) and heterozygous (Rr) family members for mutations in the ATOH7 transcription unit28,29 and known regulatory elements within 5.5kb of the initiation site, by sequencing multiple PCR amplicons (Supplementary Fig. 1, Supplementary Table 1). The proximal 1.8kb upstream segment has been conserved during vertebrate evolution28 and is sufficient to drive reporter transgene expression in embryonic mouse retinas with a pattern similar to the Math5-lacZ knock-in allele30,31.
No sequence changes were found in the ATOH7 gene or proximal upstream region compared to the NCBI reference genome or known variants. Likewise, no mutation was found in any of the 13 other genes in the nonrecombinant interval, including SIRT1, which has a major retinal dysgenesis phenotype in mutant mice (Supplementary Table 2).
To screen for remote structural changes in the ATOH7 region, we amplified single-copy DNA products from the 74 kb intergenic segment delimited by flanking genes PBLD (phenazine biosynthesis-like domain) and MYPN (myopalladin). This analysis revealed a deletion of six adjacent amplicons out of 17 distributed across the region, in samples from blind individuals (rr) compared to non-carrier relatives (RR) (Fig. 3, Supplementary Table 3). To verify this result and define endpoints, we performed multiple PCRs with primers flanking the deletion, and amplified a nested set of junctional products from NCRNA samples (Supplementary Fig. 2). The deletion spans 6523bp, with 5bp microhomology at the endpoints (Fig. 4b). It is located 21.7 to 15.2kb upstream of the ATOH7 transcription start site and 28.8 to 35.3kb downstream from the 3′ end of PBLD (Fig. 3c). We then genotyped all available samples from the pedigree (n =150), including individuals who were not part of the original linkage studies, using a triplex strategy to amplify wild-type (+) and NCRNA (Δ) alleles simultaneously (Fig. 4a). The results were completely concordant with phenotype data. All 31 blind individuals were homozygous for the deletion (ΔΔ), and all 119 sighted relatives were homozygous for the wild-type allele (++, n = 39) or heterozygous for the deletion (Δ+, n = 80, including known obligate carriers). The correlation is highly significant (χ2 = 150.17 for df = 3, giving P = 1.20 × 10−32).
To evaluate the small possibility of a second, coincidental deletion or duplication in the nonrecombinant region, we performed high-resolution SNP (single nucleotide polymorphism) analysis on six DNA samples from the NCRNA pedigree, including two blind patients and two obligate carriers (Supplementary Fig. 3). This analysis detected the deletion 20kb upstream from ATOH7 (spanning rs385145) but no further abnormalities. Unexpectedly, DNA from one blind individual (10877) was heterozygous for SNP markers centromeric to MYPN, revealing a novel recombinant disease haplotype. The underlying crossover narrows the critical region to 713kb and formally excludes four of the original candidate genes. Although an occult point mutation is theoretically possible, the combined genetic, clinical and functional data strongly support a causal relationship between the 6.5 kb deletion and NCRNA disease.
Since heterozygotes are healthy, a loss of gene function is the most logical pathogenic mechanism. The deleted region contains no annotated exons, cDNAs or ESTs (Fig. 3c). In principle, the 6.5 kb deletion could affect expression of ATOH7, PBLD or a more distant gene. Dysregulation of PBLD is an unlikely explanation for NCRNA disease pathology, because it is not expressed in the developing eye (see below).
To rigorously survey the expression of flanking genes, we performed RT-PCR analysis in mice (Supplementary Fig. 4). Although a low level of Mypn mRNA was detected in the fetal eye, only Atoh7 is specifically expressed at high levels. Because EST and 5′RACE data provide no evidence for an alternative upstream exon or mRNA splicing29, these results are most consistent with an enhancer mutation affecting the transcriptional regulation of ATOH7.
The NCRNA deletion removes a cluster of three conserved noncoding elements (CNEs) with a level of nucleotide identity among mammals that is comparable to the ATOH7 coding and proximal regulatory regions (Fig. 3c). To test the function of these CNEs in vivo, we generated transgenic mice in which 3.8 kb human genomic DNA from the deletion was joined to a minimal β-globin promoter (BG) and red fluorescent protein (nuclear Cherry) reporter (Fig. 5). This 3034-BGnCherry transgene was specifically expressed in the developing retina in 5 of 6 founder embryos and 2 of 3 stable lines. No activity was detected from the remaining two integration events. The restricted spatiotemporal pattern matches the Math5-lacZ knock-in allele, beginning at E11.5 with the onset of retinal neurogenesis and continuing to P1 in newly post-mitotic retinal precursors4,5 (Fig. 6a,b,e–). No overlap was observed between nuCherry and the S-phase marker BrdU or the M-phase marker PH3 (Fig. 6c,d,f–h). Together, our data suggest the deleted CNEs may act as a transcriptional enhancer. These CNEs are further upstream from the ATOH7 coding sequence than previously studied regions32. Although not tested directly, the NCRNA deletion may impair the level or timing of ATOH7 expression.
The mammalian ATOH7 region contains two broad regulatory units–a primary enhancer located within 2kb of the promoter, which has been characterized in mouse genomic DNA30,31, and the remote cluster of CNEs that is deleted in the NCRNA pedigree (Fig. 5). To directly compare their activities, we crossed 3034-BGnCherry and Math5-GFP mice, which carry reporters controlled by the human remote and mouse primary enhancers respectively, and examined nuCherry and GFP expression in double transgenic embryos by fluorescence and confocal microscopy, and flow cytometry (Fig. 7, Supplementary Fig. 5). These patterns were essentially overlapping, spatially and temporally, despite the fact that regulatory elements in the transgenes were derived from different species, that the primary and remote regulatory units have little sequence similarity to each other (Supplementary Figs. 6–8), and that Atoh7 is only transcribed for a few hours in each cell, beginning after terminal M phase. Flow analysis of dissociated E14.5 neural retinas revealed a single population of cells, comprising 40% of the total, that coexpress nuCherry and GFP (and presumably Atoh7) (Fig. 7c). There were few single-positive cells, which could arise from differences in folding time and perdurance of the reporter proteins. Thus, despite devastating clinical consequences of the NCRNA deletion (Figs. 1–2), the proximal and remote enhancers appear redundant in their ability to direct transcription in a transgenic context.
Unlike the coding region, which is deeply conserved among vertebrates, there is limited sequence homology between the remote enhancer, in mouse or human DNA, and the corresponding region upstream of Atoh7 (ath5) in teleosts (Supplementary Figs.7–8). Indeed, the syntenic block containing Atoh7 has been significantly disrupted by evolutionary chromosome rearrangements with breakpoints close to the promoter (Supplementary Fig. 6). To test enhancer function in a divergent species, we generated 3034-BGnCherry zebrafish carrying the same nuCherry reporter transgene. In these fish, nuCherry expression was restricted to the developing retina and correlated with the onset of neurogenesis (Fig. 8a,b). Moreover, 3034-BGnCherry fully overlapped expression of the ath5:GFP transgene33 which contains 7kb zebrafish Atoh7 regulatory DNA (Fig. 8c,d and Supplementary Fig. 9). Although distantly related by structure and ancestry, the human remote enhancer is appropriately recognized by transcriptional machinery in the teleost retina and is thus functionally conserved.
In this paper, we identify an ATOH7 regulatory deletion as the probable cause of nonsyndromic congenital retinal nonattachment, a Mendelian disorder of cell fate determination. Although circumstantial, this conclusion is based on the most parsimonious interpretation of concordant phenotypic, genetic, transgene expression and evolutionary data. Since Atoh7 controls development of the first-born retinal neurons (RGCs), the deleted enhancer, and the proximal regulatory region, appear to be the principal sites in the mammalian genome responsible for initiating retinal neurogenesis.
Although relatively rare worldwide, NCRNA is common in certain geographic areas, and other families have been reported with this disease16,17. Clinically related malformations, such as retinal dysplasia, congenital falciform folds and PHPV, have been described histologically, and may coexist with retinal nonattachment in the same pedigree or contralateral eye. To our knowledge, optic nerve aplasia has not been previously associated with retinal nonattachment. Our findings provide a basis for molecular genetic analysis of these disorders and carrier detection in the Khorasani pedigree.
The ocular pathology in NCRNA patients resembles the phenotype of Atoh7 mutant mice beyond the absence of optic nerves. In normal mice, retinal vascular development occurs between birth and P14, with spreading of surface vessels from the central artery to the ora serrata, elaboration of two deep vascular plexi, and coordinated regression of the fetal hyaloid network12. In Atoh7 mutants, there is complete agenesis of the intrinsic retinal vasculature26. Instead, the fetal vessels persist after birth and collapse upon the retinal surface, where they penetrate and ramify in an irregular pattern. These abnormal vessels are perfused by cilioretinal arteries and the remnant hyaloid trunk, and are prone to hemorrhage. A comparable but less severe neovascularization has been noted in Bst/+ mice and some human sporadic optic nerve aplasia cases34, and may reflect a general deficiency of RGC paracrine factors that control retinal angiogenesis.
In NCRNA patients, the pathogenic sequence of events is likely to be similar to that in Atoh7 mutant mice, except that the persistent fetal blood vessels appear to proliferate more dramatically during infancy, causing complete perinatal detachment of the retina. This may occur during or shortly after neurogenesis, since the pathology is apparent at birth, and it may have a tractional basis, involving fibrovascular strands between the retina, posterior lens and hyaloid vessels. Indeed, vascular traction is thought to underlie the prevalence of juvenile detachments in patients with ROP, PHPV and FEVR. RGC axons comprising the nerve fiber layer may also stabilize the retina mechanically. The difference between mouse Atoh7 and human NCRNA phenotypes may reflect the large disparity in scale between species for several parameters:  the size of the eye globe compared to the fixed distance for O2 diffusion and thickness of the neuroretinal shell;  the relative volume of the lens and vitreous chamber;  the exposure of the retina to mechanical forces; and  the surface area requiring perfusion. Consequently, the hypoxia experienced by cells in the newborn retina, and the extent of the proliferative neovascular response, may be significantly greater in NCRNA patients compared to Atoh7 mutant mice, and this may promote retinal detachment. There are also important differences in retinal structure and timing of vascular development between humans and mice, which may affect NCRNA pathogenesis13.
Regardless of these differences, the vascular phenotype in both species is a secondary consequence of the absence of RGCs, since Atoh7 expression is limited to neurons. Although the molecular pathways are not fully defined, ganglion cells are known to stimulate migration of retinal astrocytes from the optic stalk, which promotes the normal program of vascularization12. These events are mediated by signaling factors such as Shh, Pdgfa, Vegf and Ndp2. In this way, the pathogenesis of NCRNA further highlights the interdependence between neural and vascular development.
The profusion of vitreal blood vessels in NCRNA resembles the vascular events in retinopathy of prematurity, a major form of childhood blindness35. Although multiple risk factors have been identified, ROP is primarily stimulated by retinal ischemia during postnatal development, with an excessive angiogenic response. The new vessels create a fibrovascular ridge, near the equator of the globe, with a risk of folding and tractional detachment. The time course for retrolental fibroplasia in NCRNA and ROP progression are similar18,36.
In children and adults with NCRNA, recurrent intravitreal hemorrhages may communicate with the anterior chamber, causing chronic corneal blood staining, a local hemosiderin reaction. Although this effect has not been demonstrated histologically in Atoh7 mutants, hyphemas do occur in these mice26.
The role of atonal bHLH proteins in generating the first retinal neurons has been conserved during vertebrate and invertebrate evolution5,37. In the Drosophila eye imaginal disc, Atonal specifies R8 cells, which induce formation of photoreceptors R1–R7 through a succession of cell-cell interactions. Consequently, in ato mutant flies, all of these dependent cell types are missing. In Atoh7 mutant mice and lakritz zebrafish, the absence of ganglion cells does not trigger a similar neurogenic catastrophe. However, the loss of RGC inductive signals does cause a major retinovascular failure in humans with NCRNA disease.
The signals that initiate retinal neurogenesis converge upon Atoh7 transcriptional enhancers. Several extrinsic factors have been implicated, which may stimulate (FGF), suppress (GDF11) or modulate (SHH) Atoh7 expression (Supplementary Notes). Likewise, a variety of homeodomain and bHLH proteins appear to act directly upstream, including Pax6 and Neurog23,32,38. Functional binding sites have been identified within the bipartite 1.8 kb primary enhancer for mouse, Xenopus and chicken Atoh730–32,39 and are predicted in the most conserved portions of the remote enhancer (Supplementary Fig. 7d,e). These two DNA segments drive expression in the developing retinas of multiple species in an appropriate spatiotemporal pattern for Atoh7, despite differences in timing with respect to cell cycle exit. Mouse Atoh7 is expressed after terminal M phase4, whereas zebrafish ath5 is activated during G2 phase of the last cell division, such that every ath5+ cell gives rise to two daughters33 (Supplementary Fig. 9d,e). In each case, regulatory protein complexes are likely to assemble on the enhancer chromatin in mitotic progenitors prior to Atoh7 transcription, in a latent configuration, defining the RGC competence state1. The factors that mark these cells for terminal division are unknown.
On the basis of our findings, we speculate that the remote CNEs may function as a ‘shadow’ or secondary enhancer for ATOH7. This class of regulatory elements are thought to reinforce, refine or stabilize the effects of primary enhancers and may serve as vehicles for evolutionary change21,40. First, the NCRNA element is located ~20kb from the primary ATOH7 enhancer, which is close to the mRNA initiation site, and therefore presumably interacts with the promoter during retinogenesis via chromatin looping (Supplementary Fig. 10). Second, the remote and primary enhancers drive nearly identical patterns of expression, and are thus at least partly redundant (Fig. 7). The two enhancers have no obvious sequence similarity, suggesting an independent evolutionary origin, and may have different selective constraints. They appear to have evolved as separate entities in the tetrapod lineage, by progressive displacement of 5′ regulatory modules from a single ancestral cluster (Supplementary Fig. 7c) with refinement by accumulation of point mutations. Third, the remote enhancer may be important for evolutionary fitness.
Why does ATOH7 have duplicate enhancers? Each is sufficient for appropriate reporter expression in transgenic mice and, although not formally tested, both are likely to be necessary for full development of RGCs in vivo. The remote enhancer may greatly boost the level of gene expression, acting via a cooperative mechanism, or stabilize the onset of transcription at a critical early stage of retinal histogenesis. In principle, delayed development of RGCs, due to a disruption in Atoh7 expression or cellular metabolism, could prevent ganglion cell axons from exiting the optic cup before closure of the choroid fissure, leading to optic nerve hypoplasia and/or coloboma. These effects may be compounded because retinal neurogenesis in some vertebrates appears to propagate as a wave, suggesting a feed-forward mechanism. This may involve direct binding of Atoh7 protein to E-box sites in the proximal enhancer32, although there is no evidence for Atoh7 autoregulation in mammals5. Alternatively, an internal neurogenic timer may be cued by a gradient of signaling molecules41. In either model, dual enhancers would assure robust firing of Atoh7 transcripts.
Congenital optic nerve malformations are relatively common causes of severe visual impairment in humans and have environmental and genetic risk factors. Optic nerve hypoplasia (ONH) affects 12% of blind children and often occurs with hypothalamic, pituitary or other CNS abnormalities (Supplementary Notes). Mutations in several genes are known to cause ONH in humans or mice, including Pax6, Pou4f2, Rx, Hesx1, Otx2, Six6, Chx10, Sox2, Dcc, Ntn1 and Rpl24. Bilateral optic nerve aplasia is the most extreme ocular manifestation of this disease spectrum. There is also considerable variation in optic disc size among healthy people42 and RGC number among inbred strains of mice43. In humans, two recent genome-wide association studies (GWAS) implicate the ATOH7 region as the most important determinant of variation in optic disc area, which is directly related to RGC number44,45. Remarkably, the highest statistical significance for this papillometric trait is centered over the NCRNA deletion, including one genotyped SNP (rs3858145) in the shadow enhancer (Fig. 3c, Supplementary Fig. 8b).
The earliest stages of retinal neurogenesis may be relatively precarious, particularly the brief but essential burst of Atoh7 transcription, which is required for vision and thus evolutionary survival. In this way, the remote enhancer may buffer this critical process against environmental and genetic stresses, as recently demonstrated for the svb and snail loci in Drosophila46,47. The dual enhancers appear to function equivalently but have no obvious sequence similarity, apart from short motifs (Supplementary Fig. 7). In future studies, it will be important to identify the molecular factors that activate these enhancers in vivo, triggering the onset of retinal histogenesis, and to determine the transcriptional and phenotypic effects of deleting the orthologous remote enhancer in mice.
The NCRNA deletion may silence ATOH7 expression entirely, as a null allele, or prevent high-level transcription, similar to deletions of the HBB locus control region or 3′ PAX6 enhancer. A gain-of-function mechanism is unlikely, since NCRNA heterozygotes are normal, and transgenic mice with broad ectopic retinal expression of Atoh7 have no obvious histopathology48. Although remote regulatory mutations are known in human genetics, the NCRNA deletion provides a clear example of a Mendelian disease involving a duplicate or ‘shadow’ enhancer.
Detailed eye exams and photography were performed on 12 patients for this study at Emam Ali Hospital in Bojnourd, North Khorasan. In addition, the findings reflect observations from over 40 NCRNA patients examined at regional clinics during the past two decades18. Magnetic resonance images (MRI) were obtained for an 18 year-old male subject in Tehran, using a Siemens Magnetom Expert-plus instrument. The orbital MRI protocol included T1W SE-90, T2 TIR-180, T1 TSE-180 and T1 SE sequences, with fat suppression, before and after administration of gadolinium contrast. Images were collected in horizontal, coronal and oblique sagittal planes. Approval was obtained from Shahid Beheshti University Medical School and University of Michigan institutional review boards in accordance with accepted guidelines.
Genomic DNA samples from 150 pedigree members were collected as part of a previous linkage study19 and include 28 blind subjects. We amplified 184 exonic and 25 intergenic fragments from the 14 positional candidates by PCR using custom oligonucleotide primers (Supplementary Tables 1–3, Supplementary Fig. 1) and a 2:1 mixture of Platinum Taq (Invitrogen) and EXPAND Hi-Fidelity (Roche) polymerases. Segments spanning the GC-rich ATOH7 gene and 5′ flanking DNA were amplified in the presence of betaine29 as detailed (Supplementary Table 1). Screening was performed on a core panel of samples from blind (rr), obligate carrier (Rr) and wild-type (RR) individuals. To minimize potential effects of chemical inhibitors and DNA damage, a multiple displacement amplification (MDA) step (Repli-g, Qiagen) was included before PCR. Products were sequenced directly after agarose gel purification (Wizard SV gel system, Promega). Endpoints of the NCRNA deletion were determined by junctional PCR, using five different primer pairs. RR homozygotes and Rr family members were identified based on previous 10q SSLP haplotype analysis20. As expected, blind individuals were homozygous for all SNPs surveyed within the nonrecombinant interval.
To determine NCRNA genotypes in family members throughout the pedigree, and identify mutation carriers, we used a triplex PCR strategy. Diagnostic reactions (25μl) were performed in 1X PCR buffer (50 mM KCl, 20 mM Tris pH 8.4, 1.5 mM MgCl2) containing 0.3 M betaine (1X MasterAmp™, Epicentre)29, 200 μM dNTPs, 0.25 U Taq polymerase (Invitrogen), 75 ng genomic DNA and three primers (5 pmol each). The reverse primers were 30R (5′-CCCAGTGTTTCCTAGCTTTGTTAGG) and 37R (5′-TTCAGAACAGGAGCAGCGTG), and the common forward primer was 30Fref (5′-TGCAAGGACTTTGTGGCTCTTG). PCRs were performed in a 96-well MJ Thermocycler with the following protocol: 95°C × 5 min; followed by 40 cycles of 95°C denaturation × 45 sec, 55°C annealing × 45 sec, and 72°C extension × 60 sec; followed by a final 72°C extension step for 7 min, with termination at 4°C. Allelic wild-type (727 bp) and NCRNA deletion (462 bp) products were resolved on a 1% agarose gel.
Four archival samples from the NCRNA pedigree, and fresh blood and lymphoblastoid cell line DNA from one blind (rr) individual, were genotyped for 1.1 million informative SNP (single nucleotide polymorphism) markers, distributed with 2.4 kb average spacing genome-wide, including 390 SNPs in the NCRNA critical interval. Genotypes were determined using the Illumina Omni1-quad Infinium BeadChip platform and analyzed for dosage (logR ratio) and heterozygosity (B allele frequency) using GenomeStudio software.
Total RNA was isolated from adult retina, liver, kidney and heart tissues and E14.5 fetal eyes of C57BL/6 mice by Trizol extraction (Invitrogen) and treated with RNase-free DNase (Roche). cDNA was generated from 5 μg RNA using d(N)6 random hexamers, recombinant RT (Superscript II, Invitrogen) and RNaseH in a 40 μl reaction, and 1 μl of the resulting cDNA was used as template per 25 μl PCR. Products were amplified for 34 cycles with Platinum Taq (Invitrogen) using primers specific for Atoh7, Mypn, Pbld and Actb (Supplementary Fig. 3). The PCR cycles were 95°C denaturation × 45 sec, 58°C annealing × 45 sec, and 72°C extension × 60 sec.
Plasmid p3034-BGnCherry was created by inserting a 3834 bp human genomic fragment into the XbaI and BamHI sites of pBGnCherry (kindly provided by Jane Johnson), in the normal transcriptional orientation. This vector contains a minimal human β-globin promoter (BG) and a monomeric Cherry red fluorescent protein reporter cassette (nuCherry) with a synthetic amino-terminal nuclear localization signal (MAPKKKRKVEDV) downstream of the BamHI site, and it has no intrinsic activity in transgenic mice. The fragment was amplified from human BAC clone RP24-1114I5 with primers 30Fref-SpeI (CATTTAACTaGTCAGCTTGCAAGGACTTTGTGGC) and 34R (5′-GTGAATCAAAGTACAGCATGATGG) and EXPAND Hi-Fidelity polymerase, creating a terminal SpeI restriction site. The fragment was purified after digestion with SpeI and BamHI, subcloned, and verified by sequencing. It extends between −21.9 kb and −18.1 kb, relative to the ATOH7 transcription start site (TSS), and spans three CNEs within the NCRNA deletion.
Founder embryos and one stable transgenic line were generated by injecting the 5.0 kb SalI-PspOMI fragment of plasmid p3034-BGnCherry into fertilized (C57 x SJL) F2 or (B6SJL × CD-1 Atoh7tm1Gla) F1 oocytes. To directly compare nuCherry and Math5-lacZ expression, we examined embryos carrying the Atoh7tm1Gla lacZ knock-in allele5. Math5-GFP mice carry an EGFP reporter transgene regulated by a 2.1 kb mouse DNA segment that contains the Atoh7 promoter, 5′UTR and 1.8 kb upstream genomic DNA49,31. Embryos from timed pregnancies were evaluated for nuCherry and GFP expression by stereofluorescence microscopy, or were fixed in 4% paraformaldehyde (PFA), cryosectioned and immunostained. To assess cell division, gravid dams were given an intraperitoneal injection of BrdU (100 mg/kg) one hour before harvest. Wholemount GFP and nuCherry images were obtained using 480/40 or 560/40 nm excitation, and 510 or 610 nm LP emission filters, respectively, with a Leica MZ FLIII stereomicroscope.
Fixed sections or dissociated cells were blocked with 0.5% Triton X-100, 0.1 M NaPO4 pH 7.3 (PBTx) with 1% bovine serum albumin (BSA) and 10% normal donkey serum (NDS). Primary antibodies were applied overnight at 4°C in PBTx with 3% NDS 1% BSA. After washing in PBS 0.5% BSA, minimally cross-reactive donkey anti-rabbit IgG sera conjugated with biotin (1:500), or donkey anti-chicken IgY or anti-rat IgG sera conjugated with Dylight-488 (1:1000, Jackson Immunoresearch, West Grove, PA), was applied for 2 hr at room temperature. After secondary antibody staining and washing in PBS 0.5% BSA, tertiary streptavidin-conjugated Texas Red (1:400) or Dylight-594 (1:1000) reagents were applied for 1 hr in PBTx with 1 μg/ml 4′,6-diamidino-2- phenylindole (DAPI). Slides were then washed in PBS and visualized by epifluorescence (Olympus BX-51) or laser confocal microscopy (Zeiss LSM-510). Primary antibodies were anti-DsRed (Clontech rabbit polyclonal, 1:500) for Discoma sp. nuCherry; anti-GFP (Abcam chicken polyclonal, 1:2000); anti-β-galactosidase (rat polyclonal, 1:1000); and anti-BrdU (Harlan BU/175 rat monoclonal, 1:100). Genetic and serologic controls showed no cross-detection of nuCherry and GFP antigens. BrdU immunostaining was performed as a terminal procedure, following application of other 1° and 2° antibodies and an antigen retrieval step (2.4 N HCl in PBS 0.4% Triton X-100 for 1 hr). PH3 (phosphohistone H3) was detected using polyclonal rabbit antibody directly conjugated to Alexa-448 (Cell Signaling 9708, 1:50), after nuCherry immunostaining and 10% normal rabbit serum blocking steps, or a mouse monoclonal (6G3, Cell Signaling, 1:100) with biotinylated anti-mouse IgG and Alexa-647 streptavidin.
Double-positive E14.5 littermates carrying nuCherry and GFP reporters were generated, along with single-positive and double-negative controls, by crossing Math5-GFP and 3034-BGnCherry transgenic mice. Genotypes were provisionally assigned by retinal fluorescence and verified by PCR of amniotic membrane DNA. Neural retinas from each embryo were dissected in cold PBS and dissociated by digestion with 1 μg/ml trypsin (Gibco) for 15 min at 37°C, with intermittant trituration, in a 400 μl volume. The reaction was terminated by addition of 1 μg/ml SBTI (soybean trypsin inhibitor, igma), 2.5 mM MgCl2 and 100 μg/ml DNase I (Sigma). Freshly dissociated cells were kept on ice, protected from light, and sieved immediately prior to flow cytometric analysis by centrifugation at 100×g for 5 min into 70 μm mesh strainer collection tubes (BD Biosciences).
Two-dimensional flow cytometry was performed using a BD FACS Diva instrument equipped with 488 nm argon and 561 nm solid-state excitation lasers (Spectraphysics) and 510/21 nm and 610/20 nm bandpass emission filters (Omega optics) to detect GFP and nuCherry fluorescence, respectively. No spectral compensation was required. Two embryos of each genotype (n = 8) were analyzed. For each cell suspension, the primary population was gated using FSC-A (forward scatter) and SSC-A (side scatter), and multiparticle events were excluded using the FSC-W parameter. Living cells (≥95%) were further gated by DAPI (1 μg/ml) exclusion, which was detected using a 351 nm UV laser (Coherent 300) and 424/44 nm bandpass filter. Flow data from ≥10,000 events were analyzed for each embryo using BD FACS Diva software.
To visualize dissociated retinal cells, 20 μl aliquots from each suspension were plated on 8-well chamber glass slides coated with poly-D-lysine, and the adherent cells were fixed with 4% PFA. Immunocytochemistry was performed as noted above, with rabbit anti-DsRed (1:2000) and chicken anti-GFP (1:5000). Cells were viewed using an Olympus BX-51 microscope with widefield fluorescence optics and a 100X oil immersion objective.
Zebrafish in our breeding colony originated from a local pet store and are maintained at 28°C with a 10/14 hour light/dark cycle. Transgenic fish were produced by Tol2 transposase-mediated integration (Supplementary Notes). Expression cassettes (3034-BGnCherry and BGnCherry control) were subcloned between the XhoI and BsrGI sites of pT2AL200R150G, a plasmid vector that harbors Tol2 transposable elements flanking the cloning site. Capped Tol2 transposase mRNA was synthesized from NotI-linearized pCS-TP template using the Mmessage Mmachine SP6 kit (Ambion), purified by phenol-chloroform extraction, precipitated with isopropanol, and resuspended in nuclease-free water at 1 μg/ml. Approximately 0.3 μl of a mixture containing 25 ng/μl supercoiled expression plasmid DNA and 25 ng/μl Tol2 transposase mRNA was injected into fertilized zebrafish eggs at the 1–2 cell stage. Injected embryos were surveyed for nuCherry expression 2–5 days post-fertilization (dpf) using a Leica MZ FLIII stereofluorescence microscope. No expression was detected in embryos injected with the control BGnCherry construct, whereas most embryos injected with the 3034-BGnCherry exhibited bright red fluorescence. Approximately 100 embryos expressing the reporter were grown to adulthood.
To identify germline integration of the transgene, injected fish were bred in groups of 6–8, and their offspring were examined for nuCherry expression. Groups containing fluorescent progeny were outcrossed to identify founders (F0) with transgenes integrated in the germline. The resulting N1 progeny were then bred to wild- type fish to perpetuate heterozygous lines harboring solitary transgenes flanked by LTRs. Three 3034-BGnCherry lines were generated in this way. To visualize retinal nuCherry expression in transparent embryos, 0.003% 1-phenyl-2-thiourea (PTU) was added to the water at 12 hpf to inhibit pigment formation. 3034-BGnCherry progeny were also crossed to the ath5:GFP line (kindly provided by Ichiro Masai), which carries a GFP transgene regulated by 7 kb genomic DNA from zebrafish ath5 (Supplementary Notes). Mitotic cells were labeled by immersing zebrafish embryos for 20 min in 1 mM EdU (ethynyl deoxyuridine) 15% DMSO prior to fixation in 4% PFA and cryopreservation. EdU-labeled nuclei were detected using CuSO4 and Alexafluor 647 azide (Click-iT reagents, Invitrogen).
Sequence comparisons were performed between remote and primary ATOH7 enhancer elements in human and mouse genomes using pairwise BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Upstream ATOH7 regulatory regions from diverse taxa were aligned using pairwise (BLASTZ), MultiZ and PhastCons features of the UCSC browser (http://genome.ucsc.edu), for panels containing 12 eutherian mammals and 37 vertebrate species (Supplementary notes) and mVISTA (LAGAN). To search for homologous amphibian and fish sequences, we interrogated Xenopus, teleost (Danio, Oryzias, Gasterosteus, Tetraodon) and lamprey (Petromyzon) genomes and upstream regions with individual CNE segments from the NCRNA deletion, using high- and low-stringency search parameters. PCR sequence analysis and optimal primer design for mutation screening were performed using MacVector 7.2 software (Accelrys). For mutation screening, PCR sequences from the NCRNA pedigree were compared to the NCBI reference genome (NCBI36/hg18, Mar 2006 assembly) and known variants (dbSNP, including the March 2010 interim release from the 1000 genomes project). Potential regulatory motifs for DNA-binding factors were identified using multiple analysis tools, including the JASPAR database (http://jaspar.genereg.net/).
The authors are grateful to Anand Swaroop for facilitating this collaborative project; Avisa Adelavood for logistical support and encouragement; Hemant Parmar, Jonathan Trobe and Edward Oliver for help interpreting MRI studies; Lev Prasov for assistance with mouse tissue cDNAs and retinal dissections; Jane Johnson for pBGn-Cherry plasmid DNA; Nadean Brown for Math5-GFP transgenic mice; Ichiro Masai for ath5:GFP transgenic fish; Mike Pihalja and Mark Chiang for flow cytometry advice; Steve Philips and Tehmina Masud for screening candidate genes; Susan Dagenais and Robert Lyons for SNP genotyping analysis; Scott Barolo, Nadean Brown, Lev Prasov, Chris Chou, Miriam Meisler and Brian Link for helpful suggestions; and to Thom Saunders, Maggie van Keuren and the UM transgenic animal, flow cytometry and DNA sequencing cores for technical support. We are profoundly grateful to NCRNA family members for their involvement and dedication to this study over many years. The research was funded by grants from The Glaucoma Foundation (TGF) and the UM Center for Genetics in Health and Medicine (CGHM) to TG, and from the NIH to TG (EY14259), DG (EY18132) and JAB (T32EY13934).
AUTHOR CONTRIBUTIONSNMG and MM collected clinical data; DDR, JAB and TG performed genomic and functional experiments; DG developed and bred transgenic fish; NMG, DDR, JAB, DG and TG analyzed data and wrote the manuscript.