Qing-Qing Zhang
email: qqzhang@nigpas.ac.cn

André Nel
email: anel@mnhn.fr

Dany Azar
email: azar@mnhn.fr

Bo Wang
email: bowang@nigpas.ac.cn

Nel André

Azar Dany

Bo Wang Zhang

A & Wang

B New

Chinese Psocids

Eocene Fushun

New Chinese psocids from Eocene Fushun amber (Insecta: Psocodea)

published or not. The documents may come

distinctly different from all known Psocodea from Fushun amber in their three-segmented tarsi. S. fushunensis is the second unequivocal fossil of Psyllipsocidae. E. sinica is the first Asiatic and oldest representative of the Neotropical family Ptiloneuridae, and reveals a formerly global distribution of the family. The discovery of these two families in Eocene Fushun amber suggests a rather warm palaeoclimate for the Fushun amber locality.

Key words: Insecta, Psocodea, Fushun amber, China, Paleogene FUSHUN, a city located in northeastern China, Liaoning Province, yields a unique amber biota from lower Eocene deposits (Fig. 1). Fushun amber comes from the layers of the West Opencast Coalmine (41°50' N, 123°54' E), which is the largest opencast coalmine in Asia.

The mining district in the Fushun basin contains the Palaeocene Laohutai and Lizigou formations and the Eocene Guchengzi, Jijuntun, Xilutian, and Gengjiajie formations(see [START_REF] Wang | A diverse paleobiota in Early Eocene Fushun amber from China[END_REF]. The amber-bearing Eocene Guchengzi Formation includes thick coal beds intercalated with carbonaceous shale [START_REF] Wu | Thermodynamic analysis and simulation of coal metamorphism in the Fushun Basin, China[END_REF][START_REF] Meng | Palaeoclimatic evolution during Eocene and its influence on oil shale mineralisation, Fushun basin, China[END_REF][START_REF] Wang | A diverse paleobiota in Early Eocene Fushun amber from China[END_REF].

Fushun amber comes from the middle and upper coal beds of the Guchengzi Formation, which have been assigned an early-middle Ypresian age (50-53 |Ma|) constrained by the paleofloras, paleomagnetic considerations, and isotopic dating [START_REF] Wang | Climatic change during the Palaeocene to Eocene based on fossil plants from Fushun, China. Palaeogeogr, Palaeoclimatol[END_REF][START_REF] Quan | Paleogene evolution of precipitation in northeastern China supporting the middle Eocene intensification of the East Asian monsoon[END_REF]. Fushun amber contains the only Paleogene Asian (Laurasia-originated) arthropod biota that has sufficiently numerous and diverse fossils to permit detailed investigation [START_REF] Wang | A diverse paleobiota in Early Eocene Fushun amber from China[END_REF]. The palaeoclimate in the early Eocene in Fushun is reconstructed as moist subtropical, with variability in seasonal precipitation (very probably a monsoonal system), with dry and wet seasons [START_REF] Quan | Paleogene evolution of precipitation in northeastern China supporting the middle Eocene intensification of the East Asian monsoon[END_REF], 2012[START_REF] Wang | Eocene prevalence of monsoon-like climate over eastern China reflected by hydrological dynamics[END_REF].

The Chinese Fushun amber is known for its entomofauna since the work of [START_REF] Hong | Amber insects of China[END_REF].

However, interest in this fauna has not been great until the recent overview paper of [START_REF] Wang | A diverse paleobiota in Early Eocene Fushun amber from China[END_REF]. The insect fauna proved very diverse with representatives of Diptera, Hymenoptera, Coleoptera, Thysanoptera, Hemiptera and Psocodea. Here we described two new psocids that are of great interest for the biogeography of these groups and palaeoclimate of the Fushun amber locality.

Materials and methods

The specimens are preserved in pieces of relatively clear, yellow amber. The amber pieces containing the inclusions were cut, shaped, and polished. Then they were mounted between two coverslips with a Canada balsam medium as described by [START_REF] Azar | New psychodids from the Cretaceous ambers of Lebanon and France, with a discussion of Eophlebotomus connectens Cockerell, 1920 (Diptera, Psychodidae)[END_REF] before being examined and photographed. The fossils were examined and measured with incident light using an Olympus SZX9 and Leitz Wetzlar binocular microscopes.

The general characteristics of Fushun amber and its geological context were provided by [START_REF] Wang | A diverse paleobiota in Early Eocene Fushun amber from China[END_REF]. We follow the wing venation nomenclature of [START_REF] Nel | Traits and evolution of wing venation pattern in paraneopteran insects[END_REF] and the body part nomenclature of [START_REF] Mockford | North American Psocoptera (Insecta)[END_REF]. Abbreviations for wing venation are as follows: ScP subcosta posterior; RA radius anterior; RP radius posterior; r-rs crossvein between RA and RP; MA median anterior; MP median posterior; CuA cubitus anterior; CuP cubitus posterior, cua-cup brace between CuA and CuP; A1, A2 anal veins.

Systematic palaeontology

Suborder TROGIOMORPHA Roesler, 1940 Family PSYLLIPSOCIDAE Enderlein, 1911 Sinopsyllipsocus gen. nov.

Diagnosis.

Description. Body 2.0 mm long; head 0.4 mm long, 0.65 mm wide; vertex without setae, large, distance between eyes 0.4 mm, ratio of distance between eyes and greatest width of postclypeus 2:1; anterior part of head without setae; vertical epicranial suture distinct but anterior arms less distinct; frons large; postclypeus gibbous but not very prominent; ocelli disposed in triangle far from compound eyes, anterior ocellus as large as others; eye hairless, rather large, 0.15 mm diameter, reaching level of vertex; lacinia present but indistinct; mx2 with a long conical sensillum; mx4 apically broadened, with several short setae but no large sensillum, with smaller sensillum near distal end (Fig. 2C,D); scape 0.08 mm long, pedicel 0.08 mm long, both distinctly broader than flagellomeres; 22 flagellomeres (f), annulated, with two rows of three rather long setae, f1 and f2 0.01 mm long, f3 and f4 0.06 mm long, f5 and more distal flagellomeres nearly all 0.07 mm long (Fig. 3C). Dorsal surface of thorax apparently bare; thorax deformed, 0.5 mm long, 0.45 mm wide; meso-and metanotum not raised above level of head. Abdomen deformed, 1.1 mm long, 0.3 mm wide; genital organs not visible (Fig. 2A,B).

Forewing hyaline, 2.4 mm long, 0.8 mm wide, ratio of forewing length to width 3:1; wing elongate, apex rounded; basal part of ScP not reaching radius and distally evanescent, 0.25 mm long; ScP distinct, 1.4 mm long, directed posteriorly; pterostigmal area not sclerotized; base of RP nearly perpendicular to R and M; r-rs crossvein present; RP and M fusion short, 1.5 mm long; RP divided into two apical branches; M with three branches; CuA divided into two branches CuAa and CuAb, Cu1a 0.8 mm long; vein CuP and A1 distally fused in a nodulus (Fig. 3B); no A2; wing membrane bare, only a few setae disposed in one row on main veins; very few setae on C but no setae on other parts of wing margin (Fig. 3A,D).

Hindwing hyaline, 2.0 mm long, 0.65 mm wide, rather narrow, elongate, but apex rounded; ScP not visible; basiradial cell closed, narrow and elongate, triangular; RP basally fused with M in a common stem; both RP and M with two distal branches; CuA simple (Fig. 3A,E).

Legs with tarsi three-segmented; claws with very small single preapical tooth; short setae on legs, two short apical spurs on posterior tibia.

Discussion. According to the key to the psocopteran families of [START_REF] Smithers | Keys to the families and genera of Psocoptera (Arthropoda, Insecta)[END_REF],

Sinopsyllipsocus fushunensis falls in the couplet coming out as

Psyllipsocidae/Prionoglarididae because of the following combination of characters: macropterous form; tarsi three-segmented; body and wings lacking flattened scales; forewing venation complex; pterostigmal area not opaque; nodulus present; only one anal vein in forewing; lacinia present. The presence of a long, conical sensillum on mx2 favours affiliation with Prionoglarididae, but the wing venation of our fossil better corresponds to that of Psyllipsocidae in that the ScP does not reach the radius, and the areola postica is narrow and elongate (see [START_REF] Smithers | The classification and phylogeny of Psocoptera[END_REF][START_REF] Lienhard | A new genus of Prionoglarididae from a Namibian cave (Insecta: Psocoptera)[END_REF]. According to [START_REF] Mockford | North American Psocoptera (Insecta)[END_REF] and [START_REF] Baz | Archaeatropidae, a new family of Psocoptera from the Cretaceous amber of Álava, northern Spain[END_REF], the presence of more than 20 antennal flagellomeres implies an attribution to Trogiomorpha. The Psocatropetae may possess or lack a sensillum on mx2; they can have annulation or lack it on the antennae. Our fossil differs from Psyllipsocus and all other Psyllipsocidae in the presence of secondary annulation on the antennae plus the presence of a sensillum on mx2 [START_REF] Smithers | The classification and phylogeny of Psocoptera[END_REF].

More precisely, according to the key to genera of Smithers (1990, p25), S. fushunensis would be closely affiliated with Psyllipsocus because of the following characters: macropterous form with numerous veins on wings; forewing margin almost glabrous; and veins with few dorsal setae disposed in one row. In the key to Nearctic genera of Mockford (1993, p46-47), it would also be linked closely to Psyllipsocus because of the forewing being large, with RP and M branched distally, forewing broad with a ratio of forewing length to width 3.44:1, vertex broad, and meso-and metanotum not raised above the level of the head. Dorypteryx Aaron, 1883 has reduced forewings. S. fushunensis cannot be attributed to Psocatropos Ribaga, 1899 and Pseudorypteryx García-Aldrete, 1984 because the forewing CuA is forked into CuAa and CuAb. Furthermore, S. fushunensis does not belong to Pseudorypteryx because it lacks a large sensillum or smaller ones near the distal end of mx4 [START_REF] García-Aldrete | The Trogiomorpha (Psocoptera), of Chamela, Jalisco, Mexico[END_REF], Mockford 1993, p49). Khatangia [START_REF] Vishniakova | Psocoptera in Late Cretaceous insect-bearing resins from the Taimyr[END_REF], which is based on one species, K. inclusa [START_REF] Vishniakova | Psocoptera in Late Cretaceous insect-bearing resins from the Taimyr[END_REF], from the Late Cretaceous Taimyr amber, has a forewing venation similar to that of Psyllipsocus [START_REF] Vishniakova | Psocoptera in Late Cretaceous insect-bearing resins from the Taimyr[END_REF]. Sinopsyllipsocus fushunensis differs from Khatangia in the presence of a nodulus and, in its hindwing, RP being basally fused with M in a common stem and the basiradial cell consequently triangular, as in extant Psyllipsocus.

The lack of precise information on the structure of the lacinia, female genital valvulae and colouration of body and appendages of S. fushunensis prevents any detailed comparison to the modern species of Psyllopsocus. Smithers (1972, p53) Etymology. Named after the Latin for China.

Diagnosis.

Description. Body 1.56 mm long; head elongate, 0.42 mm long, 0.45 mm wide; mandible short; eyes not stalked on a long dorsolateral extension of head capsule; dorsal border of epistomal suture present (Fig. 4C,D); antennae with scape 0.1 mm long, pedicel 0.12 mm long, both distinctly broader than flagellomeres; 11 preserved flagellomeres without any annulation but with long setae disposed into two ventral rows; first flagellomere very long, 0.45 mm long; not placoid sensilla on flagellomeres; pair of setae at base of ocellar group; labial palpus triangular; labral sclerites present but with anterior ends not joined by a sclerotic strip; lacinial tip with outer cusp more developed than inner cusp. Thorax deformed, 0.45 mm long; meso-and metanotum not raised above level of head. Abdomen deformed, 0.75 mm long; female genital organs hardly visible.

Wings hyaline, without flattened scales; forewing 1.8 mm long, 0.59 mm wide; margin and veins glabrous; costa not thickened in margin of pterostigma; pterostigmal area opaque, with microsetae, 0.45 mm long, 0.12 mm wide; pterostigma 3.7 times as long as wide; nodus present; veins Rs and M connected by a short crossvein; areola postica free, 0.25 mm long, 0.14 mm wide, only 1.8 times as long as wide; CuAa and CuAb separated near posterior wing margin; two anal veins (Fig. 5A,C). Hindwing 1.4 mm long, 0.42 mm wide; anal vein long, parallel to wing margin; very small setae on veins Rs and M; row of fine setae on proximal end of costal margin (Fig. 5B,D).

Legs with tarsi three-segmented; pretarsal claw with a preapical denticle; pulvillus of pretarsal claw straight and pointed (Fig. 4E); ctenidiobothria present on t1 and t2 of female mid-and hind leg.

Discussion. Eotriplocania sinica can be attributed to the Psocomorpha due to the presence of 11 flagellomeres without any annulation, the presence of a nodus on the forewing, and a sclerotized pterostigma. It can be attributed to the Epipsocetae on the basis of the following apomorphies, as defined by Casasola-González ( 2006): a pair of setae at the base of the ocellar group; labral sclerites present; lacinial tip with outer cusp more developed than inner cusp; ctenidiobothria present on t1 and t2 of female hind leg; presence of a second anal vein According to the key to the psocopteran families of [START_REF] Smithers | Keys to the families and genera of Psocoptera (Arthropoda, Insecta)[END_REF], Eotriplocania sinica also belongs to Ptiloneuridae because of the following combination of characters: macropterous form, tarsi three-segmented, body and wings lacking flattened scales, forewing venation complex, pterostigmal area thickened, forewing with nodulus and two anal veins; head elongate, labrum with two ridges; forewing veins lacking spurs.

According to the keys to ptiloneurid genera of [START_REF] García-Aldrete | New genera of Psocoptera (Insecta), from Mexico, Belize and Ecuador (Psoquillidae, Ptiloneuridae, Lachesillidae)[END_REF] and [START_REF] Azar | Two new Ptiloneuridae from Colombian copal (Psocodea: Psocomorpha)[END_REF], Eotriplocania sinica would fall near the two Neotropical genera Triplocania [START_REF] Roesler | Neue und wenig bekannte Copeognathengattungen. 1[END_REF]Belicania García-Aldrete, 2006 based on possession of the following characters: hindwing M single-branched; forewing 2A joining wing margin; no crossveins between A1 and wing margin; forewing areola postica high with apex rounded; labral sclerite incomplete, not reaching anterior margin of labrum; and forewing M three-branched. The differences between Eotriplocania and the two Neotropical genera concern the structure of the hypandrium, a structure unknown in our fossil.

Eotriplocania sinica differs from genus Belicania in the Ptiloneuridae in the areola postica being shorter, nearly as long as wide. This fossil differs from all modern species of Triplocania in the Ptiloneuridae by its longer stem of Rs basal of crossvein rs-m. Moreover, most Triplocania species have a longer rs-m than in this fossil, except T. magnifica [START_REF] Roesler | Neue und wenig bekannte Copeognathengattungen. 1[END_REF].

Discussion

Sinocladiopsocus fossilis and Eotriplocania sinica differ from all the Psocodea described by [START_REF] Hong | Amber insects of China[END_REF] from the Fushun amber in their three-segmented rather than two-segmented tarsi. S. fushunensis is the second unequivocal fossil of Psyllipsocidae, and the other is from the Eocene French Oise amber. Modern psyllipsocids live in caves, on rock surfaces, and among dead leaves, and they are restricted to the tropical regions of South America, Africa and Asia, with a few species from North America and one cosmopolitan (García-Aldrete 1997). E. sinica is the oldest record of Ptiloneuridae, of which extant species are all Neotropical. The fossil record of Ptiloneuridae is extremely rare. Until now, only two fossil species have been described from Colombian copal [START_REF] Azar | Two new Ptiloneuridae from Colombian copal (Psocodea: Psocomorpha)[END_REF]. Our fossil is the first Asiatic representative of Ptiloneuridae, and reveals a formerly global distribution of the family (Fig. 6). The poor record of the family, however, precludes any detailed speculation about dispersal patterns. In summary, both families are nowdays distributed in tropical regions, and thus our findings suggest a rather warm palaeoclimate for the Fushun amber locality [START_REF] Wang | Climatic change during the Palaeocene to Eocene based on fossil plants from Fushun, China. Palaeogeogr, Palaeoclimatol[END_REF][START_REF] Quan | Paleogene evolution of precipitation in northeastern China supporting the middle Eocene intensification of the East Asian monsoon[END_REF].

  indicated that the modern Psyllopsocus spp. have antennae with 22 segments or more (see also[START_REF] New | An unusual new species of Psyllipsocidae (Psocoptera) from Malaysia[END_REF], as in S. fushunensis.To date four fossils have been assigned to Psyllipsocidae: Psyllipsocus eocenicusNel et al., 2005 (lowermost Eocene French Oise amber), ?Psyllipsocus banksiCockerell, 1916 (mid- Cretaceous Burmese amber), a "Psyllipsocus sp. (nymph)" from Lower Miocene Mexican amber[START_REF] Mockford | Fossil insects of the Order Psocoptera from Tertiary amber of Chiapas, Mexico[END_REF], and Khatangia inclusa[START_REF] Vishniakova | Psocoptera in Late Cretaceous insect-bearing resins from the Taimyr[END_REF] (Late Cretaceous Taimyr amber) but the last three systematic positions are dubious and require revision. Psyllipsocus eocenicus Nel et al., 2005 differs from S. fushunensis in the presence of only 18 flagellomeres instead of 20. ?Psyllipsocus banksi Cockerell, 1916 from upper Albian Burmese amber, was listed by Spahr (1992) and Ross &York (2000, p12), among others, but has never been critically evaluated. Its original attribution was based on the opinion of Nathan Banks: 'The Psocid you figure must be close to Psyllipsocus de Selys; …' (Cockerell 1916). The original description and figures lack several characters diagnostic of the family and genus. The type specimen needs a redescription before any definite attribution of this taxon can be made. Mockford (1969) described a 'Psyllipsocus sp. (nymph)' from the Miocene Mexican amber. Infraorder Epipsocetae Pearman, 1936 Family PTILONEURIDAE Eertmoed, 1973 Eotriplocania gen. nov. Type species. Eotriplocania sinica gen. et sp. nov. Etymology. Named after the Eocene period and Triplocania Roesler, 1940. Gender feminine. Diagnosis. This fossil genus is characterized by the following combination of characters: 11 flagellomeres; crossvein between Rs and M in forewing very short; areola postica high with apex rounded, nearly as wide as long; forewing 2A joining wing margin; no crossveins between A1 and wing margin; forewing M three-branched; stem of Rs basal of crossvein rs-m long; rs-m very short; hindwing M single-branched; labral sclerite incomplete, not reaching anterior margin of labrum. Eotriplocania sinica sp. nov. (Figs 4, 5) Material. Holotype specimen 12048, deposited in the Nanjing Institute of Geology and Palaeontology, China. Age and occurrence. Early Eocene; Guchengzi Formation, West Opencast Coalmine, Fushun, Liaoning, China.

(

  A2) in forewing; very small setae on veins Rs and M of hindwing; veins Rs and M of forewing with a (short) cross-vein between them. The character state 'row of fine setae on proximal end of costal margin of macropterous hind wing present' is an apomorphy of the clade [Ptiloneuridae, Spurostigma + Cladiopsocus (Cladiopsocidae) and Epipsocidae] whereas the three plesiomorphic character states 'pulvillus of pretarsal claw straight', 'pulvillus apex pointed' and 'pterostigma only 3.7 times as long as wide' would exclude Eotriplocania sinica from the family Dolabellopsocidae. The rather long pterostigma and the anterior ends of labral sclerites not joined by a sclerotic strip are characters of Ptiloneuridae rather than of [Spurostigma+Cladiopsocus (Cladiopsocidae) and Epipsocidae].

Fig. 2 .

 2 Fig. 2. Sinopsyllipsocus fushunensis gen. et sp. nov., holotype 11036. A, photograph of dorsal

Fig. 3 .

 3 Fig. 3. Sinopsyllipsocus fushunensis gen. et sp. nov., holotype 11036. A, photograph of wings;

Fig. 4 .

 4 Fig. 4. Eotriplocania sinica gen. et sp. nov., holotype 12048. A, photograph of dorsal features;

Fig. 5 .

 5 Fig. 5. Eotriplocania sinica gen. et sp. nov., holotype 12048. A, photograph of forewing; B,

  Type species. Sinopsyllipsocus fushunensis gen. et sp. nov.Etymology. Named after the Latin for China, and Psyllipsocus. Gender masculine.

	Diagnosis. Psyllipsocid with secondary annulation on antennae; 22 flagellomeres; long
	sensillum on mx2; no large sensillum or smaller ones near distal end of mx4; macropterous,
	with numerous veins on wings, forewing margin almost glabrous, veins with few dorsal setae
	disposed in one row; ScP not reaching radius; areola postica present, narrow elongate.
	Sinopsyllipsocus fushunensis sp. nov. (Figs 2, 3)
	Material. Holotype specimen 11036, deposited in the Nanjing Institute of Geology and
	Palaeontology, China.
	Etymology. Named after Fushun City.
	Age and occurrence. Early Eocene; Guchengzi Formation, West Opencast Coalmine, Fushun,
	Liaoning, China.

Acknowledgements

12

We are grateful to Prof. E.A. Jarzembowski for improving the English language of the manuscript. This research was supported by the National Basic Research Program of China (2012CB821900), and National Natural Science Foundation of China (41572010, J1210006).

List of figures.


Dr Gerrit Hilgen
email: gerrit.hilgen@ncl.ac.uk

Sahar Pirmoradian

Daniela Pamplona

Pierre Kornprobst

Bruno Cessac

Matthias H Hennig

Evelyne Sernagor

Pan-retinal characterisation of Light Responses from Ganglion Cells in the Developing Mouse Retina Abbreviated Title: Pan-retinal Development of Light Responses

We have investigated the ontogeny of light-driven responses in mouse retinal ganglion cells (RGCs).

Using a large-scale, high-density multielectrode array, we recorded from hundreds to thousands of RGCs simultaneously at pan-retinal level, including dorsal and ventral locations. Responses to different contrasts not only revealed a complex developmental profile for ON, OFF and ON-OFF RGC types, but also unveiled differences between dorsal and ventral RGCs. At eye-opening, dorsal RGCs of all types were more responsive to light, perhaps indicating an environmental priority to nest viewing for pre-weaning pups. The developmental profile of ON and OFF RGCs exhibited antagonistic behavior, with the strongest ON responses shortly after eye-opening, followed by an increase in the strength of OFF responses later on. Further, we found that with maturation receptive field (RF) center sizes decrease, responses to light get stronger, and centers become more circular while seeing differences in all of them between RGC types. These findings show that retinal functionality is not spatially homogeneous, likely reflecting ecological requirements that favour the early development of dorsal retina, and reflecting different roles in vision in the mature animal.

Introduction

The onset of visual experience in mouse occurs around postnatal day (P) 12, at eye opening. Although the retina cannot experience patterned vision beforehand, it is remarkable that RGCs are already capable of encoding information originating from photoreceptors and transmit it to retinal central targets as soon as eyes open. However, these early light responses are far from mature, and they progressively acquire their adult features while the retina develops [START_REF] Masland | Maturation of function in the developing rabbit retina[END_REF][START_REF] Tootle | Early postnatal development of visual function in ganglion cells of the cat retina[END_REF][START_REF] Sernagor | Emergence of complex receptive field properties of ganglion cells in the developing turtle retina[END_REF][START_REF] Sernagor | Retinal Development[END_REF] . In mouse, RGC dendritic stratification in the ON and OFF layers of the inner plexiform layer matures after eye opening [START_REF] Tian | Developmental mechanisms that regulate retinal ganglion cell dendritic morphology[END_REF] and light-driven activity guides the refinement of synaptic connectivity [START_REF] Tian | Visual stimulation is required for refinement on ON and OFF pathways in postnatal retina[END_REF][START_REF] Tian | Visual deprivation alters development of synaptic function in inner retina after eye opening[END_REF] . Consequently, RF sizes [START_REF] Cantrell | Non-Centered Spike-Triggered Covariance Analysis Reveals Neurotrophin-3 as a Developmental Regulator of Receptive Field Properties of ON-OFF Retinal Ganglion Cells[END_REF][START_REF] Koehler | Receptive field center size decreases and firing properties mature in ON and OFF retinal ganglion cells after eye opening in the mouse[END_REF] and complex RF properties such as direction and orientation selectivity [START_REF] Elstrott | Selectivity in the Retina Is Established Independent of Visual Experience and Cholinergic Retinal Waves[END_REF][START_REF] Chen | orientation-selective retinal ganglion cells in mice[END_REF] keep maturing after the onset of visual experience. Yet, despite ongoing maturation after eye opening, longitudinal studies of RF properties have never been fully documented [START_REF] Sernagor | Emergence of complex receptive field properties of ganglion cells in the developing turtle retina[END_REF] .

One important yet often neglected issue is that the retina is not uniformly organised from a functional perspective. Indeed, dorsal, ventral, nasal and temporal domains have evolved to enable optimal not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016; 3 encoding of specific features in the visual scene. For example, mouse cones co-express medium wavelength and short wavelength opsins (M-opsin and S-opsin), with a dorsal-to-ventral increasing gradient in S-opsin (and opposite for M-opsin) [START_REF] Szel | Unique topographic separation of two spectral classes of cones in the mouse retina[END_REF][START_REF] Szél | Distribution of cone photoreceptors in the mammalian retina[END_REF][START_REF]The murine cone photoreceptor: a single cone type expresses both S and M opsins with retinal spatial patterning[END_REF][START_REF] Röhlich | Two different visual pigments in one retinal cone cell[END_REF][START_REF] Lukáts | Photopigment coexpression in mammals: Comparative and developmental aspects[END_REF] . These dorso-ventral gradients affect RGC responses in adult animals with respect to their spectral tuning [START_REF] Chang | Chromatic Coding from Cone-type Unselective Circuits in the Mouse Retina[END_REF][START_REF] Ekesten | Cone and rod inputs to murine retinal ganglion cells: evidence of cone opsin specific channels[END_REF][START_REF] Wang | Spectral and Temporal Sensitivity of Cone-Mediated Responses in Mouse Retinal Ganglion Cells[END_REF] , improving encoding of achromatic contrasts [START_REF] Baden | A tale of two retinal domains: Near-Optimal sampling of achromatic contrasts in natural scenes through asymmetric photoreceptor distribution[END_REF][START_REF] Yin | Physiology and morphology of coloropponent ganglion cells in a retina expressing a dual gradient of S and M opsins[END_REF] and providing evolutionary advantages for visual tasks [START_REF] Peichl | Diversity of mammalian photoreceptor properties: Adaptations to habitat and lifestyle?[END_REF] . The topographical organisation of some RGC subtypes also exhibits dorso-ventral non-uniformity [START_REF] Zhang | The most numerous ganglion cell type of the mouse retina is a selective feature detector[END_REF][START_REF] Bleckert | represented by nonmatching topographies of distinct mouse retinal ganglion cell types[END_REF] . However, nothing is known about the developmental consequences of these inhomogeneities.

Here we present a longitudinal study of RGC RF properties in the developing mouse retina from eye opening up to maturity with emphasis on dorso-ventral topographical differences. We recorded simultaneously from hundreds to thousands of RGCs at near pan-retinal level using the high-density large-scale CMOS-based Active Pixel Sensor multielectrode array (APS-MEA) featuring 4096 electrodes (42 µm pitch) arranged in a 64x64 configuration, covering an active area of 7.12 mm 2 25-27 , allowing us to discriminate topographical differences in light responses. We classified RGCs as ON, OFF and ON-OFF, measuring basic firing properties such as latency, peak amplitude and response duration at different contrast levels. We completed our study by determining the spatio-temporal properties of their RF central areas throughout development using a novel, high resolution stimulus approach.

Results

Simultaneous pan-retinal recording from the dorsal and the ventral retina

The spatial extent (7.12 mm 2 ) of the APS-MEA chip allowed us to record simultaneously from dorsal and ventral RGCs over large retinal areas (Fig. 2A). The small electrode pitch (42 µm) enables sampling from many individual RGCs from these areas, providing us with an unbiased very large analytical sample size (see Table 1) and helps to reduce the amount of experiments/animals needed. Figure 2A shows that it is possible to visualise the outline of the retina, to quantify activity levels on individual channels and to delineate various retinal areas on the APS-MEA chip simply by looking at spiking activity. For each channel the log spike count of a full field stimulus experiment (see Methods) from a P13 (Fig. 2a) and an adult (P38, Fig. 2c) retina were pseudo color-coded and plotted according to their not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016; 4 position on the 64x64 APS-MEA chip. This generates activity maps showing that the outline of the P13 retina is smaller compared to the P38 retina and that the spike rate is overall higher in all channels in that younger retina. Figures 2b and2c illustrate responses to these same full-field stimuli in both retinas following spike sorting and RGC classification (see Methods), yielding spike rasters and histograms for dorsally and ventrally located ON (green), OFF (purple) and ON-OFF (blue) RGCs. The responses of P13 ON cells to the alternating full field stimulus were much stronger and more sustained compared to the responses of OFF and ON-OFF RGCs in the same retina (Fig. 2b) and to all responses in the P38 retina (Fig. 2d). But OFF and ON-OFF cell responses were more clearly defined at P38 than at P13.

Interestingly, the ON cell responses were stronger in the dorsal than in the ventral part of the retina (Fig. 2b, compare left and right). In contrast, at P38, ON responses were somewhat stronger and briefer in the ventral part (Fig. 2d right) compared to the dorsal part, whereas it was the opposite for OFF responses (Fig. 2d left). This simple pan-retinal visualisation demonstrates that the development of basic firing properties varies not only for different RGC types but even for dorsal versus ventral ON and OFF RGCs. Next, we quantified these stimulus-driven responses for different RGC types at different ages (Fig. 3), and we examined how responses change when the full-field stimulus was presented with different Michelson contrasts (Fig. 456).

Ontogeny of dorsal and ventral light response features for different RGC types

Full field stimuli were presented to retinas of different ages, and post-stimulus time histograms (PSTH) were generated for every RGC. The results were classified into 4 age groups: P13 (4 retinas, 4183 RGCs in total), P16/P17 (2 P16 and 2 P17 retinas, 4404 RGCs), P19 (3 retinas, 3210 RGCs) and adult (2 P29 and 2 P38 retinas, 5811 RGCs) and further divided into ON, OFF and ON-OFF (see Methods) dorsal and ventral RGCs (Table 1). For the results in Figures 3456we are referring to these groups. From individual PSTHs, we extracted the peak amplitude, time to peak and response duration for light onset and offset (see Fig. 1a). Fig. 3a summarises results for ON (left), OFF (middle) and ON-OFF (right) RGCs at all ages (x-axis), illustrating the dorsal (green) and ventral (blue) mean peak amplitude A1 (ON cells), A2 (OFF cells) and A1+A2 (ON-OFF) (see Fig. 1a). In line with examples in Fig. 2, dorsal not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016; 5 P13 ON cells were significantly more active (Fig 3a, left) than OFF and ON-OFF cells in both dorsal and ventral areas. There is also a progressive decrease in A1 from P13 to adulthood for dorsal ON cells, whereas A2 increases with age for all OFF cells, dorsal and ventral (Fig. 3a, middle). Except for P13, A2 values were significantly higher in dorsal OFF than in ventral OFF cells. ON-OFF cells, on the other hand, did not show significant changes between the youngest and oldest age groups, and their responses were overall weaker compared to the other cell types (Fig. 3a, right).

How fast these responses to full field stimuli are at different ages was evaluated by measuring the time from stimulus on-or offset to the peak amplitude (T2P1 and T2P2, respectively). We found that T2P1 in ON cells and T2P2 in OFF cells progressively decrease from P13 to adulthood, with slower responses in OFF cells (Fig. 3b, middle). ON-OFF cells had overall slower response onsets than ON and OFF cells. Overall, ventrally located RGCs exhibited slightly more sluggish responses.

How prolonged a RGC response is was defined by measuring the response duration following the peak (RD1 and RD2, respectively). The spike rasters in Fig. 2b clearly show that dorsal ON cells respond in a much more sustained fashion at P13. This observation was confirmed by plotting RD1 and RD2 values (Fig. 3c, left), revealing that dorsal ON responses are significantly more sustained at P13. Subsequently, RD1 values are similar in dorsal and ventral ON cells, gradually dropping down until adulthood (Fig. 3c, left). In contrast, OFF cells become moderately more sustained between P13 and adulthood, eventually producing slightly longer adult responses than ON cells (Fig. 3c, middle). Finally, ON-OFF cells have overall shorter responses to light, and the only developmental change we observed was more prolonged responses in dorsal RGCs at P13 (Fig. 3c,right).

Taken together, dorsal light responses are more prominent after eye-opening and the peak amplitudes and response durations of ON and OFF cells showed an antagonistic behaviour from eye opening to maturity (with the most pronounced changes in ON cells) whereas the ON-OFF cell responses exhibit less developmental variability.

Ventral and dorsal differencepeak amplitude

Next, we investigated how these developmental changes vary with the stimulus contrast. Using the same age groups as above, we recorded full field responses with different Michelson contrasts (MC, see not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016; 6 Methods). Figure 4a illustrates the firing peak amplitudes for all ON cells (A1) from a P13 retina for three different full field contrasts (0.41, 0.53, 0.62) plotted with respect to their electrode position on the array. Overall, the dorsal side clearly shows stronger activity (yellow), especially for MC 0.62. To establish whether this dorsal-ventral trend is present in all retinas of the same age group and/or between the age groups, we calculated A1 and A2 for all full field contrasts for all RGC types (Fig. 4b) but also separately for ON (Fig. 4c), OFF (Fig. 4d) and ON-OFF (Fig. 4e) types. At P13, all dorsal RGC types exhibited significantly higher peak firing rates than ventral cells at all contrast levels (Fig. 4b-e left).

The highest values and most pronounced differences between dorsal and ventral cells were observed at 0.62 Michelson contrast rather than at the highest contrast (MC 0.7). At P16/17 and P19, dorsal OFF RGCs were also significantly more active than their ventral counterparts (Fig. 4d middle), whereas there were no significant differences between ON-OFF cells (Fig. 4e middle). At P19, the trend reversed with ventral ON cells becoming significantly more active than their dorsal counterparts for the three strongest contrast levels (Fig. 4c middle) whereas that trend moved towards lower contrast levels in adult retinas and virtually disappeared at high contrasts (Fig. 4c, d right). Interestingly, and in contrast to P13 (Fig. 4e left), ventral adult ON-OFF cells were significantly more active (for both A1 and A2) than their dorsal counterparts at all but the weakest contrast levels. Taken together, dorsal units are strikingly more active after eye-opening, an effect that disappears in older retinas.

Ventral and dorsal differencetime to peak

Figure 5a shows T2P2 values for all OFF cells from a mature P38 retina for three different full field contrasts (MC 0.19, 0.53, 0.7). Here, dorsal responses appear more sluggish for lower contrasts. Similar trends were observed for all adult RGC types (Fig. 5b-e right), with dorsal RGCs showing significantly longer times to peak than ventral cells, especially prominent at lower contrasts. However, at P13 all ventral RGCs had more sluggish responses at the highest MC levels (Fig. 5b-e left) whereas this was shifted towards lower contrast in OFF P16/P17 and P19 RGCs. There were no differences between dorsal and ventral ON and ON-OFF cells at P19 and P16/P17 (Fig. 5c, e middle). In summary, there is a negative developmental shift of the response time to peak for all RGC types, with the most sluggish responses found in ventral RGCs after eye-opening and slower responses in dorsal units in adults.

not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016;

Ventral and dorsal differenceresponse duration

Figure 6a illustrates response durations for all ON cells in a P13 retina for three different contrast levels (MC 0.19, 0.62, 0.7), with longer responses developing at higher contrasts in the dorsal part of the retina., At P13 dorsal RGCs exhibit more prolonged responses at most contrast levels (Fig. 6b-e left)., and there is a similar trend for ON cells at P16/17 (Fig. 6c, middle). In contrast, the trend was inverted for OFF cells at P16/17, with longer responses in ventral cells (Fig. 6d, middle). As for adult retinas, differences between dorsal and ventral cells became virtually non-existent, except for small differences at higher contrast levels (e.g. longer ventral OFF responses at MC 0.7 at P19; longer dorsal OFF responses MC 0.62, 0.67 and 0.7 in adult retinas).

In summary, there are significant differences between the dorsal and ventral response properties (peak firing, time to peak and response duration) for different RGC types from eye-opening to maturity. Shortly after eye opening, dorsal RGCs exhibit much stronger and more sustained responses to full field flashes than ventral cells. These dorsal-to-ventral differences fade as development progresses, and interestingly, in many cases the gradient even reverses to ventral-to-dorsal.

Maturation of RGC RF centers

We used a novel checkerboard stimulus, shifted white noise (SWN, see Methods), to characterise the maturation of ON and OFF RGC RF centers from eye-opening to maturity. As shown in Pamplona et al. [START_REF] Pamplona | A super-resolution approach for receptive fields estimation of neuronal ensembles[END_REF] , shifted white noise stimulus considerably improves the spatio-temporal resolution of the receptive field computed via STA. This allows therefore a better accuracy in receptive field characterisation. Such an accuracy is necessary to thoroughly describe the variability in receptive fields during development.

The STA is the mean stimulus that precedes a spike and the cell type is inferred from the STA polarity for the phase closest to the spike time 0 (negative for OFF and positive for ON RGCs). Figure 7 shows the temporal STA signal strength (see Methods) and the spatial profile for selected ON (yellow, Fig. 7a) and OFF (blue, Fig. 7b) example cells from all used age groups. The most striking developmental change we observed is that both negative and positive STA peak signal strength values are much weaker not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint ( We did not observe significant differences in STA strength between dorsal and ventral ON or OFF cells at most ages, but overall, we found that responses were significantly stronger in adult RGCs.

We next used the STA to measure RF diameters (see Methods) and quantify their developmental changes from eye-opening to maturity. Figure 7e and show that RF diameters are largest at P13 both in ON (e) and OFF (f) RGCs, dropping down to a minimum at P16, followed by a marginal, nonsignificant increase at P19 and then a more significant increase from P19 to adulthood. Note that P13 ON RF diameters were significantly larger than adult RF diameters but that is not true for P13 OFF RF diameters. Since we observed significant differences in firing properties between dorsal and ventral RGCs, we also looked at RF sizes separately for these same two groups, but we found no major differences (except for adult, where ventral ON RF diameters were significantly larger than dorsal ones).

Finally, we quantified the shape of RFs by measuring their eccentricity (see Methods). As shown in Figures 7g andh, RF eccentricity for both ON and OFF RGCs significantly decreases with development; that is, the shape of RFs becomes closer to circular in later development. We did not observe significant differences between dorsal and ventral eccentricity in ON RGCs except for ventral RGCs being more circular than dorsal RGC at P13. On the opposite, dorsal and ventral RF circularity were significant differet from P16 onwards in OFF RGCs (Fig 7h ,right).

Discussion

In this comprehensive study we have shown that the basic features characterising RGC light responses have a developmental profile that depends both on cell type and retinal location. Furthermore, these properties do not develop synchronously in the dorsal and ventral retina. Finally, using a novel stimulus, we have been able to reliably characterise RGC receptive fields across development, which is difficult in young retinas where light responses to conventional white noise stimuli are weak. We found that ON RGCs have stronger responses to light than other cell types just after eye opening.

As development progresses, these responses become weaker while OFF responses gain strength (antagonistic to ON cells). The relatively high levels of spontaneous activity (residual of spontaneous waves) could not explain these results because peak amplitudes were normalised to baseline activity.

Moreover, if spontaneous activity was indeed affecting our measurements, it would equally affect all RGC types, and not just ON cells. An earlier study reported that RGC light responses reach their peak strength around P28, but this study pooled ON, OFF and ON-OFF responses [START_REF] Tian | Visual deprivation alters development of synaptic function in inner retina after eye opening[END_REF] . Pooling our data would yield similar results because of the antagonistic time course of changes in ON and OFF peak amplitudes throughout development (see Fig. 3a). The different maturation time course of ON and OFF peak firing responses may stem from differences in synaptic connectivity and membrane excitability between these cells. A previous study in mouse reported that light-evoked excitatory postsynaptic potentials are larger in ON than in OFF RGCs at P13 [START_REF] He | Light-evoked synaptic activity of retinal ganglion and amacrine cells is regulated in developing mouse retina[END_REF] , but it also stated that both excitatory and inhibitory light-driven synaptic responses are downregulated with age, which cannot explain why we see differences between ON and OFF cells at later ages. An alternative explanation is that the expression of ionotropic glutamate receptors (AMPA/Kainate and NMDA) differs in ON and OFF RGCs. The NR2A subunit of the NMDA receptor is predominantly found at OFF synapses, while NR2B subunits are preferentially located at ON synapses in rat RGCs [START_REF] Zhang | Subunit-and pathway-specific localization of NMDA receptors and scaffolding proteins at ganglion cell synapses in rat retina[END_REF] . However, Stafford et al 31 found no evidence for a differential localisation of the NR2B subunits at ON and OFF synapses in direction selective RGCs until P28, but no data is yet available for non-direction selective RGCs. Another explanation for our observations is that the time course of inhibitory neurotransmitter receptor expression is different for particular cell types. A good candidate is the GABAc receptor, expressed presynaptically on bipolar cell axon terminals [START_REF] Lukasiewicz | A diversity of GABA receptors in the retina[END_REF] . Retinal GABAC receptor knockout results in stronger and more prolonged spiking activity [START_REF] Lukasiewicz | GABAC receptormediated inhibition in the retina[END_REF] , similar to our observations in young ON RGCs. Further, GABAC-receptor mediated inhibition affects the regulation of glutamatergic synaptic transmission at the ON but not at the OFF bipolar-RGC synapse [START_REF] Sagdullaev | Presynaptic Inhibition Modulates Spillover, Creating Distinct Dynamic Response Ranges of Sensory Output[END_REF] . In contrast, Schubert et al., [START_REF] Schubert | Development of Presynaptic Inhibition Onto Retinal Bipolar Cell Axon not peer-reviewed) is the author/funder[END_REF] stated that ON and OFF bipolar cell spontaneous inhibitory postsynaptic currents (IPSCs) are relatively similar at P12. However, the nature of light-evoked GABAergic inputs onto developing mouse bipolar cells, and how it differs for different RGC types remains to be determined (but see Discussion in [START_REF] He | Light-evoked synaptic activity of retinal ganglion and amacrine cells is regulated in developing mouse retina[END_REF] ).

not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint ( The response durations of ON and OFF RGCs (Fig. 3c) in our study showed a similar antagonistic behaviour as the peak amplitudes (Fig. 3a), but with a decrease much stronger in ON than in OFF RGCs, and a massive drop during the fourth postnatal week (Fig. 3c). In addition to the aforementioned points, GABAC receptors consist of ρ subunits and the ρ1 subunit expression peaks around P9 (start P6) and ρ2, around P15 (start P9) [START_REF] Greka | Expression of GABA C receptor r1 and r2 subunits during development of the mouse retina[END_REF] . The difference in the timing of peak expression of these subunits may explain the differences we observe between ON and OFF RGCs. Response latencies did not exhibit conspicuous differences between different RGC types from eye opening to adulthood and the changes we observed most likely reflect ongoing activity-mediated refinement of the bipolar-RGC synapse [START_REF] Tootle | Early postnatal development of visual function in ganglion cells of the cat retina[END_REF][START_REF] Tian | Developmental mechanisms that regulate retinal ganglion cell dendritic morphology[END_REF][START_REF] Tian | Visual deprivation alters development of synaptic function in inner retina after eye opening[END_REF][START_REF] He | Light-evoked synaptic activity of retinal ganglion and amacrine cells is regulated in developing mouse retina[END_REF][START_REF] Wang | Unique functional properties of on and off pathways in the developing mammalian retina[END_REF][START_REF] Myhr | Mechanisms underlying developmental changes in the firing patterns of ON and OFF retinal ganglion cells during refinement of their central projections[END_REF] .

In mouse, S-opsins are co-expressed with M-opsins in cone photoreceptors and there are more S-opsins expressed in the ventral retina. S-opsin exhibit peak excitation around 360 nm [START_REF] Wang | Spectral and Temporal Sensitivity of Cone-Mediated Responses in Mouse Retinal Ganglion Cells[END_REF][START_REF] Jacobs | Retinal receptors in rodents maximally sensitive to ultraviolet light[END_REF] which is in the UV spectrum. To stimulate retinas in our study, we have used a broad spectrum white light composed of red, green and blue LED lights (~420 -660 nm) but not UV light. Therefore it is highly unlikely that we were able to stimulate S-opsin properly. Under these conditions, dorsal RGCs have an advantage because our white light source is biased towards activation of M-opsin, which is more prevalent in the dorsal retina. Yet, our results do not demonstrate a general increase of response strength in the dorsal side per se, but reveal a differential dorsal/ventral development pattern different RGC types through maturation. Even though a bias in M-opsin activation could potentially explain the stronger responses we recorded at P13, it cannot be responsible for the results we found at later developmental stages, including a switch to stronger ventral responses at low contrast levels in the adult retina. S-and Mopsins are expressed before eye-opening, respectively from ~P1 and ~P8 [START_REF] Szél | Distribution of cone photoreceptors in the mammalian retina[END_REF][START_REF]The murine cone photoreceptor: a single cone type expresses both S and M opsins with retinal spatial patterning[END_REF] . Therefore we can reasonably assume that both S-and M-opsin gradients are already fully established shortly after eye opening, when we sampled the earliest light responses, and differences in opsin expression are unlikely to explain our findings. Further, the light intensity (mean luminance 11 cd/m 2 ) in our experiments was set to co-activate rods and cones. Therefore rod-mediated (by rhodopsin in the visible spectrum) responses are not saturated in our recordings [START_REF] Naarendorp | Dark Light, Rod Saturation, and the Absolute and Incremental Sensitivity of Mouse Cone Vision[END_REF] , and they probably reflect a large proportion of our recordings because rods outnumber cones by 35:1 in the mouse retina [START_REF] Jeon | The Major Cell Populations of the Mouse Retina[END_REF] . Stronger responses in the dorsal retina at the onset of visual experience is not an unreasonable possibility from an ecological and not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016; evolutionary point of view because at that age, pups are still gathered in the nest, near the mother, with no ecological need to look skywards, but they rather concentrate on the nest scenery, using dorsal retinal vision. In the ensuing week, retinal circuits keep developing and refining, leading to maturation of the ventral circuitry as well.

Presenting white noise checkerboard images and performing a post hoc reverse correlation is now a standard approach to estimate RGC RF centers [START_REF] Chichilnisky | A simple white noise analysis of neuronal light responses[END_REF][START_REF] Rieke | Spikes: Exploring the neural code[END_REF] . However, the technique is challenging in very young retinas because of high levels of spontaneous activity and the weak sensitivity to light in young RGCs immediately after eye opening. The problem can be alleviated by increasing the number of trials and/or the pixel size without compromising the minimum resolution for reliable RF estimation [START_REF] Cantrell | Non-Centered Spike-Triggered Covariance Analysis Reveals Neurotrophin-3 as a Developmental Regulator of Receptive Field Properties of ON-OFF Retinal Ganglion Cells[END_REF][START_REF] Koehler | Receptive field center size decreases and firing properties mature in ON and OFF retinal ganglion cells after eye opening in the mouse[END_REF][START_REF] Schwartz | Spike-triggered neural characterization[END_REF] . Here we applied a novel approach, the SWN, which uses large checkerboard pixels (hence evoking strong responses) randomly shifted in space and time by a fraction of the pixel size, yielding better resolution. SWN images are uncorrelated across time, and the bias introduced by the correlation in space due to finite pixel size is negligible. Using this approach, we have shown that STA signal strength defined by the z-score (see Methods) peaks at P19 for ON and OFF RGCs which is in line with Cantrell et al [START_REF] Cantrell | Non-Centered Spike-Triggered Covariance Analysis Reveals Neurotrophin-3 as a Developmental Regulator of Receptive Field Properties of ON-OFF Retinal Ganglion Cells[END_REF] , who showed a signal strength peak at P18. The significantly higher signal strength in adult ventral RGCs is consistent for ON and OFF RGCs but antagonistic to our full field results. Further studies are thus needed to investigate the spatial response inhomogeneities in adult retinas. The main developmental increase in STA signal strength peak correlates well with the end of the programmed cell death period for bipolar cells and outer rods [START_REF] Young | Cell death during differentiation of the retina in the mouse[END_REF] , suggesting that the synaptic connectivity required for mature RF center responses is complete at that time, and only minor synaptic refinements occur later on, resulting in further RF expansion. We found that immature (P13) RF diameters in ON RGCs are significantly larger than later in development, while for OFF cells, RFs are similar in size at P13 and in adults, with a temporary drop during the third postnatal week. Previous studies which used conventional checkerboard images present a different developmental picture. Koehler et al. [START_REF] Koehler | Receptive field center size decreases and firing properties mature in ON and OFF retinal ganglion cells after eye opening in the mouse[END_REF] showed that ON and OFF RFs are smaller in adults than immediately after eye-opening whereas Cantrell et al. [START_REF] Cantrell | Non-Centered Spike-Triggered Covariance Analysis Reveals Neurotrophin-3 as a Developmental Regulator of Receptive Field Properties of ON-OFF Retinal Ganglion Cells[END_REF] stated that OFF, but not ON RFs expand between P15-18 and that there is a significant difference for ON, but not for OFF cells between P18-25. This demonstrates that different checker board approaches can yield different results and will be elaborated in another study [START_REF] Pamplona | A super-resolution approach for receptive fields estimation of neuronal ensembles[END_REF] . Additional factors such as mean luminance, adaptational state not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016; and recording length may account for these differences, and results should therefore be interpreted with caution. From our measurements of RF eccentricity, we found that RFs become more circular between P16-19. These changes are probably caused by the same factors responsible for developmental changes in STA signal strength. Interestingly this goes in line with a previous study in turtle where the authors could show that anisotropic properties of spatiotemporal RFs decrease with maturation [START_REF] Sernagor | Emergence of complex receptive field properties of ganglion cells in the developing turtle retina[END_REF] . The significant difference between dorsal and ventral eccentricity in ON P13 needs further investigation but it seems that higher STA signal strength predicts more isotropic RF centers (data not shown but see dorsal and ventral boxes in Fig 7c, d andFig 7g, h).

In summary, using a large-scale, high-density multielectrode array, we have investigated the ontogeny of light responses in RGCs in the developing mouse retina. We were able to record from hundreds to thousands of RGCs simultaneously at pan-retinal level, and found that the refinement of RF properties strikingly differs between RGC types and between the same RGC types in the dorsal and ventral retina.

These findings suggest that retinal functionality is not spatially uniform and that there might be an ecological advantage to favouring the development of dorsal light responses before the rest of the retina reaches functional maturity.

Materials & Methods

Retina Preparation

Experimental procedures were approved by the UK Home Office, Animals (Scientific procedures) Act 1986. Male and female wild-type mice (C57bl/6), housed under a 12 hour light-dark cycle and aged between postnatal days (P) 13-63 were used for the experiments. Mice were dark-adapted overnight and killed by cervical dislocation. Eyes were enucleated, and following removal of the cornea, lens, and vitreous body, they were placed in artificial cerebrospinal fluid (aCSF) containing the following (in mM): 118 NaCl, 25 NaHCO3, 1 NaH2 PO4, 3 KCl, 1 MgCl2, 2 CaCl2, 10 glucose, and 0.5 l-Glutamine, equilibrated with 95% O2 and 5% CO2. The ventral and dorsal orientation was marked after enucleation and confirmed by using vascular landmarks in the retina [START_REF] Wei | Two-photon targeted recording of GFP-expressing neurons for light responses and live-cell imaging in the mouse retina[END_REF] The retina was isolated from the eye cup and flattened for MEA recordings and the ventral-dorsal, nasal-temporal orientation was noted down. All not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016; 13 procedures were performed in dim red light and the room was maintained in darkness throughout the experiment.

APS-MEA recordings

The isolated retina was placed, RGC layer facing down, onto the APS-MEA and flattened by placing a small piece of translucent polyester membrane filter (Sterlitech Corp., Kent, WA, USA) on the retina followed by a home-made anchor. Similarly as described elsewhere [START_REF] Maccione | Following the ontogeny of retinal waves: pan-retinal recordings of population dynamics in the neonatal mouse[END_REF][START_REF] Portelli | Rank order coding: a retinal information decoding strategy revealed by largescale multielectrode array retinal recordings[END_REF] , throughout recording, retinas were maintained at 33°C using an in-line heater (Warner Instruments LLC, Hamden, CT, USA) and continuously perfused using a peristaltic pump (~1 ml min -1 ). Pan-retinal recordings were performed the BioCam4096 platform with BioChips 4096S+ (3Brain GmbH, Lanquart, Switzerland), integrating 4096 square microelectrodes (21 × 21 μm, pitch 42 μm ) on an active area of 2.67 × 2.67 mm. The platform records at a sampling rate of 7.1 kHz/electrode when using the full 64×64 array and recordings were stored at 12 bits resolution per channel with a 8 kHz low-pass filter/0.8 Khz high-pass filter using 3Brain's BrainWave software application. Data management and analysis for large-scale, high-density APS-MEA recording is particular difficult (as reviewed in [START_REF] Maccione | bioinformatics and neurocomputation for massive neuronal recordings in brain circuits with large scale multielectrode array probes[END_REF] ). To reliably extract spikes from the raw traces we used a quantile-based event detection [START_REF] Maccione | Following the ontogeny of retinal waves: pan-retinal recordings of population dynamics in the neonatal mouse[END_REF][START_REF] Muthmann | Spike Detection for Large Neural Populations Using High Density Multielectrode Arrays[END_REF] and single-unit spikes were sorted using the T-Distribution Expectation-Maximisation algorithm in Offline Sorter (Plexon Inc, Dallas, USA). Sorted units that exhibited at least >0.1 spikes/sec on average over the entire recording session were then verified by visual inspection of the detected clusters in the 2/3D principal component feature space, the calculated cluster interspike intervals with respect to the refractory period and the shape of the spike waveforms in the Offline Sorter GUI. Due to the high density of the electrodes, the same units were sometimes detected on multiple neighbouring channels. These redundant units were removed by comparing coincident spikes between neighbouring units. Briefly, for each unit, spikes occurring within +-2 frames (1 Frame = 1/7.06 ms) were detected in all units on the four closest electrodes and marked. This was done for all units, and units with more than 5% coincident spikes were iteratively removed such that for each coincident group only the one with the largest spike count was retained. On average we finally record from multiple hundreds to thousands individual RGCs in one experiment.

not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016;

14

Light stimuli

Light stimuli were projected onto the retina as described previously [START_REF] Portelli | Rank order coding: a retinal information decoding strategy revealed by largescale multielectrode array retinal recordings[END_REF][START_REF] Maccione | bioinformatics and neurocomputation for massive neuronal recordings in brain circuits with large scale multielectrode array probes[END_REF] by exploiting the temporal and spatial resolution of the experimental platform and were attenuated using neutral density filters to high mesopic light levels (mean luminance 11 cd/m 2 ).

A full field stimulus that switched from light to dark (0.5 Hz, 30 repetitions) was used to define peak response latency, duration and relative amplitude of ON and OFF responses (see Fig. 1A). The luminance contrast for this stimulus measured as the Michelson contrast was defined as (Imax -Imin)/(Imax + Imin) where Imax and Imin are respectively the maximum and minimum luminance and had a maximal value of 0.70. We also used full field stimuli with a series of increasing Michelson contrasts (0.19, 0.41, 0.53, 0.62, 0.67). We estimated each unit's instantaneous firing rate for the different full field intensities by convolving its spike train with a Gaussian kernel smoothing function (standard deviation = 25ms).

We then averaged the trials (Fig. 1A) and extracted several features like the amplitude of ON and OFF responses (A1, A2), the time from stimulus onset or offset to peak of these response (T2P1, T2P2 respectively) and the response duration (RD1, RD2). Statistical significance was evaluated using Oneway ANOVA with a Bonferroni post-hoc test (Prism, GraphPad, CA).

To classify RGCs according to their main response polarity, we measured the relative amplitude of ON and OFF responses and calculated the Bias Index (BI) defined as (A1 -A2)/(A1 + A2) (Carcieri et al., 2003). We used the BI to classify the cells into OFF (BI -1 to -0.33), ON-OFF (BI -0.33 to 0.33) and ON cells (BI 0.33 to 1).

Shifted White Noise

Checkerboard stimuli are routinely used to measure RF areas [START_REF] Chichilnisky | A simple white noise analysis of neuronal light responses[END_REF][START_REF] Rieke | Spikes: Exploring the neural code[END_REF] . Finer resolution can be achieved using smaller unitary checkerboard pixels, but at the same time very small checkerboard pixels may not be able to elicit reliable and repeatable responses, if at all, necessitating to reach a trade-off. Here we use an improved checkerboard stimulus, so-called shifted white noise (SWN), where checkerboard pixels are shifted randomly in space at fixed time steps [START_REF] Pamplona | A super-resolution approach for receptive fields estimation of neuronal ensembles[END_REF] . With this novel approach, the checkerboard pixel size is large enough to reliably evoke significant responses, but at the same time, the RF resolution can be very fine, given by the shift size (Fig. 1B). The SWN at position x, y and time t is defined as: not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint (

Images are composed of B checkerboard pixels with a fixed size P, the area of each pixel is defined by a rectangular function ΠP, b is the block index and its top left corner coordinates are (Xb, Yb). The index of each pixel is given by the random variable ω(b,t) which is taken from a Bernoulli distribution of values -1 (color black) and 1 (color white) with equal probability 0.5 for each block b at each time stamp t. Each image is randomly shifted horizontally (shift x(t)) and vertically (shift y(t)). shifts x(t) and The spike triggered average (STA) [START_REF] Chichilnisky | A simple white noise analysis of neuronal light responses[END_REF] was calculated by computing the average stimulus 500 ms (corresponding to 15 checkerboard frames) before a spike occurred. At the time point of the positive or negative peak maximum of the average temporal STA, a two-dimensional Gaussian was fitted to the corresponding spatial profile frame and an ellipse was drawn around the center with 1 SD of the Gaussian fit. The RF diameter was defined as the diameter of a circle with the same area as the ellipse (2*radius = 2 SD). We also measured the eccentricity of the fitted ellipse to see how 'out of round' an RF is. The eccentricity e is given by:

2 2 1 a b e   , ( 2
)
where a is the length of the semi-major axis (half of the longest diameter of an ellipse) and b is the length of the semi-minor axis (half of the shortest diameter of an ellipse). The eccentricity changes from 0 to 1, where the eccentricity of a circle is zero.

not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016; 16 To measure STA strength at the time of the positive or negative peak maximum of the average temporal STA, we computed the z-score of STA peak amplitudes; that is for a given cell the mean of STA amplitudes was subtracted from the peak amplitude and divided by their standard deviation. This results in a much finer resolution (P/S) for the spatial STA profile. and ON-OFF (e) RGC types and all age groups. All plot conventions are like for Figure 3. The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016;

Tables, Figures and Legends

  at P13-16 than later during development. Quantifying the STA strength (see Methods and legend) for all ON (Fig 7c) and OFF (Fig 7d) cells show that the STA strength significantly increases during development with the main change between P16-19.

  y(t) are random variables taking S possible values with a probability 1/S. Here we take S=4. They are redrawn at each time step. We used 17 x 17 pixels (B = 289) with an edge size of 160 μm (P = 160) and a shift step size of P/4 (40 μm) and we considered 4 steps (S = 0, 40, 80, 120) both for horizontal and vertical shifts. As shown in Pamplona et al, this form of stimulus considerably improves the spatial and temporal resolution of the STA[START_REF] Pamplona | A super-resolution approach for receptive fields estimation of neuronal ensembles[END_REF] . The Michelson contrast was 0.7 with the same mean luminance stated before and SWN images were presented for 33 ms each (30 Hz, ~45 min, ~15 min for adult retinas).

(

Figure 1

 1 Figure 1Response parameter calculations and shifted white noise paradigm. a) Top row is showing a spike raster of a full field stimulus (30 trials, 2 sec on, 2 sec off) where each dot is representing a spike, bottom row shows the averaged trials (see methods, converted to Spikes/Sec). Baseline firing rate was estimated by looking at each unit's activity (Spikes/Sec) before stimuli were presented (>30 sec). A1 (Spikes/Sec) is the first maximum peak (measured from each unit's baseline) after stimulus onset and A2 the first maximum peak after stimulus offset. T2P1 (in ms) is the time from stimulus onset to A1, T2P2 the time from stimulus offset, respectively. RD1 (in ms) is the time, so-called response duration, after A1 where the response is still above A1/e. RD2 is the response duration after A2 to A2/e. b) Top is showing the standard white noise paradigm with static located pixels and the resulting resolution (P) of the spatial STA profile. Bottom is showing the shifted white noise paradigm but with shifted (S) pixel location.

  not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr.26, 2016;

Figure 2

 2 Figure 2 Activity maps and spike raster plots of a P13 and a P38 retina. a, c) The Log spike count (full field stimulus experiment) for each channel from a P13 (a) and a P38 (c) retina is pseudo color-coded and plotted according to electrode coordinates (64 x 64 array). This results in a visualisation of the retina outline and gives an overall estimation of the number of active channels. c, d) Spike raster plots from the same P13 (b) and P38 (d) experiment used for a) and c), respectively, but after spike sorting. Each dot is representing a spike in an alternating full field stimulus experiment and dots are color-coded: green = ON RGCs, purple = OFF RGCs, blue = ON-OFF and the raster/rate plots are divided into dorsal (b,d left) and ventral RGCs (b, d right). The binned (25 ms) average response (Spikes/Sec) of all RGCs is plotted below the raster plots.

  not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr.26, 2016;

Figure 3

 3 Figure 3 Response properties of different RGC types from P13 to adult. a, b, c: respectively illustrate peak amplitude (A1, A2), time to peak (T2P1, T2P2) and response duration (RD1, RD2) for dorsal (green) and ventral (blue) ON (left), OFF (middle) and ON-OFF (right) RGCs for each age group (mean valueswith 95% confidence interval, n: see Table1). For ON-OFF cells the light on-and offset values were calculated as ON-OFF 1 (solid line) and ON-OFF 2 (dotted line), respectively. Significance: * = p < 0.05; ** = p < 0.01; *** = p < 0.001. Dark asterisks are for onset values, grey asterisks for offset values, no asterisks mean not significant.

24 Figure 4

 244 Figure 4 Dorsal-ventral gradient of peak amplitudes to different contrasts after eye-opening. a) ON peak amplitudes (A1 in Spikes/Sec) for three different full field contrasts (0.41, 0.53, 0.62) from a single P13 retina are plotted in pseudo colors according to firing strength and electrode position. For visualisation, the individual x,y electrode position for the peak values was slightly randomly shifted (+/-0.25) because after spike sorting multiple RGC units are assigned to the same electrode position. b-e) Mean peak response amplitudes to different full field Michelson contrasts (0.19, 0.41, 0.53, 0.62, 0.67) for all (b), ON (c), OFF (d) and ON-OFF (e) RGC types for all age groups (ascending from left to right) with respect to their dorsal (green) or ventral (blue) location. All plot conventions are like for Figure 3.

25 Figure 5

 255 Figure 5 Dorsal-ventral differences for the time to peak values at lower contrast in the adult retina. a) OFF RGC color-coded response times to peak (T2P2 in ms) from a single P38 at three different full field contrasts (0.19, 0.53, 0.7) plotted according to their electrode position. b-e) Mean times to peak plotted for the different full field Michelson contrasts (x axis, 0.19, 0.41, 0.53, 0.62, 0.67) for all (b), ON (c), OFF (d)

  not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.The copyright holder for this preprint (which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr.26, 2016;

Figure 6

 6 Figure 6 Responses to light are more sustained for all dorsal RGC types than for their ventral counterparts. a) The ON RGC response duration times (RD1 in ms) from a single P13 retina of three different full field contrasts (0.19, 0.62, 0.7) were pseudo colour-coded plotted according to their electrode position. b-e) The mean response duration times to different full field Michelson contrasts (x axis, 0.19, 0.41, 0.53, 0.62, 0.67) for all (b), ON (c), OFF (d) and ON-OFF (e) RGC types and all age groups. All plot conventions are like for Figure 3.

27 Figure 7

 277 Figure 7 Shifted White Noise used to study the development of RGC RF central areas. a, b: STA signal strength (top) and spatial profiles (bottom) for selected example RF central areas from ON (a) and OFF (c) RGCs from all 4 age groups (see below). c-h: box plots (whiskers: 10-90 percentile, mean indicated by + symbol) of ON (yellow boxes; P13-N=2 retinas, n = 449 cells; P16-N=2, n=425; P19-N=2, n=926; adult-N=2; n=698) and OFF (blue boxes; P13-N=2, n=309; P16-N=2, n=239; P19-N=2, n=794; adult-N=2, n=514) RGCs from all retinal areas for STA signal strength (c, d), RF diameters (e, f), and RF (g, h). All values are also shown separately for dorsal and ventral located RGCs (small black boxes). Significance: * = p < 0.05; ** = p < 0.01; *** = p < 0.001; ns = not significant.

  not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

  not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

	bioRxiv preprint first posted online Apr. 26, 2016;	doi:	. http://dx.doi.org/10.1101/050393	The copyright holder for this preprint (which was
			9

  which was . http://dx.doi.org/10.1101/050393 doi: bioRxiv preprint first posted online Apr. 26, 2016;

	10

  𝐼(𝑥, 𝑦, 𝑡) = ∑𝛱 𝑃 (𝑥 -𝑋 𝑏 -𝜀 𝑆 (𝑥, 𝑡))𝛱 𝑃 (𝑦 -𝑌 𝑏 -𝜀 𝑆 (𝑦, 𝑡))𝜔(𝑏, 𝑡)

	bioRxiv preprint first posted online Apr. 26, 2016;	doi:	. http://dx.doi.org/10.1101/050393	which was
	𝐵 𝑏=1

Table 1 Numbers

 1 peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

	(4 retinas)

(n) of dorsal and ventral ON, OFF and ON-OFF RGCs for the different age groups used for Fig.

3-6

. not

The research received financial support from the 7th Framework Program for Research of the European Commission (Grant agreement no 600847: RENVISION, project of the Future and Emerging Technologies (FET) program Neuro-bio-inspired systems (NBIS) FET-Proactive Initiative)). We thank the IIT-NetS3 laboratory (member of the RENVISION consortium) for providing the light-stimulation system integrated with the CMOS-MEA platform used in this work.

not peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.

Conflict of Interest

The authors declare no competing financial interests.