Does Background Matter? A Comparative Characterization of Mouse Models of Autosomal Retinitis Pigmentosa rd1 and Pde6b-KO

Many retinal degenerative diseases result in vision impairment or permanent blindness due to photoreceptor loss or dysfunction. It has been observed that Pde6brd1 mice (rd1), which carry a spontaneous nonsense mutation in the pde6b gene, have a strong phenotypic similarity to patients suffering from autosomal recessive retinitis pigmentosa. In this study, we present a novel mouse model of retinitis pigmentosa generated through pde6b gene knockout using CRISPR/Cas9 technology. We compare this Pde6b-KO mouse model to the rd1 mouse model to gain insights into the progression of retinal degeneration. The functional assessment of the mouse retina and the tracking of degeneration dynamics were performed using electrophysiological methods, while retinal morphology was analyzed through histology techniques. Interestingly, the Pde6b-KO mouse model demonstrated a higher amplitude of photoresponse than the rd1 model of the same age. At postnatal day 12, the thickness of the photoreceptor layer in both mouse models did not significantly differ from that of control animals; however, by day 15, a substantial reduction was observed. Notably, the decline in the number of photoreceptors in the rd1 model occurred at a significantly faster rate. These findings suggest that the C3H background may play a significant role in the early stages of retinal degeneration.


Introduction
Retinitis pigmentosa (RP) is a group of hereditary diseases that lead to visual impairment or permanent blindness due to the loss or dysfunction of photoreceptors and the retinal pigment epithelium.RP is characterized by pronounced genetic and clinical heterogeneity; therefore, the clinical age of the onset of symptoms varies greatly and ranges from childhood to adulthood [1].To date, autosomal recessive RP has been associated with mutations in at least 74 genes (RetNet: https://sph.uth.edu;accessed 5 July 2023).Recessive biallelic mutations in the pde6b gene are a common cause of autosomal recessive RP in both humans and animals [2,3].According to the Ensembl genome browser database, more than a thousand SNPs in the pde6b gene have been found in humans, leading to the development of RP.Since the mouse pde6b gene has an orthologous intron-exon relationship comparable to that of the human pde6b gene [4], the common phenotypic characteristics of Pde6b inactivation appear, which include the acute degeneration of photoreceptors and their nuclei in the outer nuclear layer, damage to retinal vessels, pigment spots on the fundus, and a decrease in electroretinogram (ERG) amplitudes [5].The pde6b gene encodes the catalytic β-subunit of cyclic guanidine monophosphate (cGMP) phosphodiesterase (Pde6b) expressed specifically in rods' outer segments.Since Pde6b plays a key role in the conversion of light into electrical response of rod photoreceptor cells, mutations in the pde6b gene result in the permanent opening of cGMP-gated cation channels in the photoreceptor membrane, allowing excess extracellular calcium ions to enter these cells, ultimately leading to cell death through apoptosis.
Currently, there are two in vivo mouse models with spontaneous mutations in the pde6b gene that are most studied: rd1 (Pde6b rd1 ) and rd10 (Pde6b rd10 ) [3].Both models exhibit different phenotypes (allelic heterogeneity).The rd1 phenotype is associated with two changes in the pde6b locus of chromosome 5: murine leukemia virus insertion (Xmv-28) in reverse orientation (3 to 5 ) in intron 1 and nonsense mutation (p.Y347X, c.C1041A) in exon 7 of the Pde6β subunit.The rd10 phenotype carries a missense mutation (p.R560C, c.C1678T) in exon 13 of the pde6b gene and is another model of recessive retinal degeneration.Unlike rd1, in which Pde6b protein production and activity are not detected, rd10 mice are characterized by reduced Pde6b activity and relatively slower onset of retinal degeneration, starting at 4 weeks [6,7].Residual Pde6b activity in the rd10 retina may prevent the toxic increase in cGMP early in mouse life, thus delaying the onset of photoreceptor cell death and widening the window for therapeutic interventions.
In this work, we established a new mouse model of RP disease with the knockout for the pde6b gene.In contrast to the rd1 mouse line, in the pde6b gene of our model, a frameshift occurred in exon 1, and there was no insertion of the mouse leukemia virus.

The Development and Biochemical Characterization of the Pde6b Knockout Mouse Model in Comparison with the rd1 Mouse Model
As we sought to establish a new knockout mouse line lacking the Pde6b protein, we designed sgRNA targeting the first exon of the corresponding gene.A combination of sgRNA and the Cas9 protein was microinjected into mouse zygotes, and 16 F0 mice were genotyped.We selected the F0 mouse model harboring two-nucleotide deletion in the first exon of the pde6b gene, leading to a nonsense codon in position 104 (Figure 1).This mouse model was crossed with a wild-type C57Bl/6 mouse model in order to obtain heterozygous (F1) mice.F1 mice were then crossed to obtain homozygotes and subsequently genotyped via PCR and DNA sequencing analysis.
Furthermore, to eliminate potential off-targets, six backcrosses of homozygotes with wild-type C57Bl/6 mice were performed.Each backcross involved successively crossing homozygous mice with wild-type C57Bl/6 mice and then crossing the resulting heterozygotes with each other to produce homozygotes.
Next, we compared the level of pde6b mRNA in the retinas of wild-type, rd1, and Pde6b-KO mice.We found that on the ninth postnatal day, the amount of pde6b mRNA in both rd1 and Pde6b-KO mice was drastically decreased (Figure 1).We did not find any detectable amounts of the Pde6b protein both in rd1 and Pde6b-KO retinas on day 12.This corresponds with previously published results on rd1 mice [8].detectable amounts of the Pde6b protein both in rd1 and Pde6b-KO retinas on day 12.This corresponds with previously published results on rd1 mice [8].

Electroretinography of Pde6b Knockout Mice in Comparison with rd1 Mice
To assess the functional state of the retina, an electroretinogram (ERG) recording is usually performed in order to measure the electrical response of retinal cells to a light stimulus [9].The ERG method allows for the calculation of the amplitude and time to the peak of the two primary response components known as the a-wave and the b-wave.The a-wave represents the initial response of photoreceptors to a short flash of light [10], whereas the b-wave reflects the response of underlying retinal neurons, including ON bipolar cells to photoreceptor stimulation.The loss of the a-or b-wave amplitude can be attributed to a number of retinal dystrophies, while "supernormal" waves with increased amplitude have been attributed to cone dystrophy [11].An analysis of photoresponses demonstrated that the shapes of the ERG on the 18th day of life in Pde6b-KO and rd1 mice were significantly different from that in wild-type animals (Figure 2a,b).Since the rod signaling cascade in rd1 and Pde6b-KO mice lacks a key effector enzyme, they never respond to low-intensity stimuli.Nevertheless, dark-adapted Pde6b-KO mice still exhibited a b-wave response to the bright flash on days 18-21, unlike rd1 mice.However, by the 24th day of life, neither the b-waves nor a-waves were detectable in both strains of mice, indicating a complete loss of visual function (Figure 2c).The responses of both strains on the 27th day were hardly detectable over noise level, and no differences were revealed between them (Figure S1).

Electroretinography of Pde6b Knockout Mice in Comparison with rd1 Mice
To assess the functional state of the retina, an electroretinogram (ERG) recording is usually performed in order to measure the electrical response of retinal cells to a light stimulus [9].The ERG method allows for the calculation of the amplitude and time to the peak of the two primary response components known as the a-wave and the b-wave.The a-wave represents the initial response of photoreceptors to a short flash of light [10], whereas the b-wave reflects the response of underlying retinal neurons, including ON bipolar cells to photoreceptor stimulation.The loss of the a-or b-wave amplitude can be attributed to a number of retinal dystrophies, while "supernormal" waves with increased amplitude have been attributed to cone dystrophy [11].An analysis of photoresponses demonstrated that the shapes of the ERG on the 18th day of life in Pde6b-KO and rd1 mice were significantly different from that in wild-type animals (Figure 2a,b).Since the rod signaling cascade in rd1 and Pde6b-KO mice lacks a key effector enzyme, they never respond to low-intensity stimuli.Nevertheless, dark-adapted Pde6b-KO mice still exhibited a b-wave response to the bright flash on days 18-21, unlike rd1 mice.However, by the 24th day of life, neither the b-waves nor a-waves were detectable in both strains of mice, indicating a complete loss of visual function (Figure 2c).The responses of both strains on the 27th day were hardly detectable over noise level, and no differences were revealed between them (Figure S1).Pde6b-KO and rd1 mice on the 18th day of the postnatal period demonstrated reduced photoresponse in both lines compared to the wild type.The response was measured using a stimulating 10 ms flash with an intensity of 31.6 cd*s/m 2 ; (c) comparison of b-wave amplitudes in Pde6b-KO and rd1 mice (n = 6) on days 18, 21, and 24 of the postnatal period revealed that the Pde6b-KO photoresponse was higher than that of the rd1 model.The Mann-Whitney test revealed statistically significant differences (* p < 0.05).

Retinal Morphology of Pde6b Knockout Mice in Comparison with rd1 Mice
In wild-type mice, normal retinal development was observed: the outer nuclear layer (ONL), the inner nuclear layer (INL), and retinal lamination were maintained at all experimental time points (12, 15, 18, and 21 days of postnatal age), and retinal thickness increased, especially in the ONL.
There was a dramatic decrease in retinal ONL thickness in both models compared to the wild type on the 15th postnatal day.On the 18th postnatal day in both models, there was a decrease in the thickness of the ONL, but in rd1 mice, the ONL consisted of one row of cell nuclei (3 ± 1 μm).In Pde6b-KO mice, the layer thickness was 2-3 cells (10 ± 3 μm).On the 21st day of postnatal age, the ONL was reduced to one row of cell nuclei with the relative preservation of the INL in both rd1 and Pde6b-KO mice (Figure 3a).
The thickness of the ONL gradually decreased in rd1 and Pde6b-KO mice but did not significantly change over time in wild-type mice.On the 15th postnatal day, the ONL was significantly thicker in Pde6b-KO mice at all measured points, compared with those of rd1 mice (p < 0.05 for all).By the 21st postnatal day, the ONL thickness in rd1 mice and Pde6b-KO mice was not significantly different (Figure 3b).(n = 6) on days 18, 21, and 24 of the postnatal period revealed that the Pde6b-KO photoresponse was higher than that of the rd1 model.The Mann-Whitney test revealed statistically significant differences (* p < 0.05).

Retinal Morphology of Pde6b Knockout Mice in Comparison with rd1 Mice
In wild-type mice, normal retinal development was observed: the outer nuclear layer (ONL), the inner nuclear layer (INL), and retinal lamination were maintained at all experimental time points (12, 15, 18, and 21 days of postnatal age), and retinal thickness increased, especially in the ONL.
There was a dramatic decrease in retinal ONL thickness in both models compared to the wild type on the 15th postnatal day.On the 18th postnatal day in both models, there was a decrease in the thickness of the ONL, but in rd1 mice, the ONL consisted of one row of cell nuclei (3 ± 1 µm).In Pde6b-KO mice, the layer thickness was 2-3 cells (10 ± 3 µm).On the 21st day of postnatal age, the ONL was reduced to one row of cell nuclei with the relative preservation of the INL in both rd1 and Pde6b-KO mice (Figure 3a).

Functional Analysis of Genes
In our previous experiments, differences in retina activity and degradation between Pde6b-KO mice and rd1 mice were observed.The Pde6b-KO mouse line had a C57Bl/6 background, and rd1 mice originated from the C3H mouse line, so we hypothesized that this slight difference in retinopathy development could be due to alterations in the genetic background of these mice.
We analyzed 502 genes related to eye development and function that have polymorphisms in C3H mice to assess their role in vision perception with GO enrichment analysis.The genes coding proteins specific to photoreceptors are the most over-represented in this gene set (Figure 4b).Also, the selected gene subcluster is most enriched by the genes described with GO molecular function terms "GO:0015267: Channel activity" and "GO:0008140: Passive transmembrane transporters activity" (Fig- ure 4c).The thickness of the ONL gradually decreased in rd1 and Pde6b-KO mice but did not significantly change over time in wild-type mice.On the 15th postnatal day, the ONL was significantly thicker in Pde6b-KO mice at all measured points, compared with those of rd1 mice (p < 0.05 for all).By the 21st postnatal day, the ONL thickness in rd1 mice and Pde6b-KO mice was not significantly different (Figure 3b).

Functional Analysis of Genes
In our previous experiments, differences in retina activity and degradation between Pde6b-KO mice and rd1 mice were observed.The Pde6b-KO mouse line had a C57Bl/6 background, and rd1 mice originated from the C3H mouse line, so we hypothesized that this slight difference in retinopathy development could be due to alterations in the genetic background of these mice.
We analyzed 502 genes related to eye development and function that have polymorphisms in C3H mice to assess their role in vision perception with GO enrichment analysis.
The genes coding proteins specific to photoreceptors are the most over-represented in this gene set (Figure 4b).Also, the selected gene subcluster is most enriched by the genes described with GO molecular function terms "GO:0015267: Channel activity" and "GO:0008140: Passive transmembrane transporters activity" (Figure 4c).tions of 2 amino acid residues in positions 971 and 972.The significance of amino acid substitutions was estimated using web tools SIFT [12], PolyPhen2 [13], and PhD-SNP [14].All substitutions were identified as "neutral" with PhD-SNP and "tolerated" with SIFT.Only with PolyPhen2 were we able to characterize several polymorphisms in the gene Rpgrip1 as possibly damaging or probably damaging mutations (refer to Table S2).Taking into account that other tools did not reveal any damaging effect of these mutations, the results obtained with PolyPhen2 should be considered with caution, as this tool is characterized by a low specificity rate, with 15-20% of false-positive results [15,16].The deletion of two glutamic acid residues was localized on the interdomain spacer but could change the protein charge.We assessed the level of Rpgrip1 mRNA in rd1 and Pde6b-KO mice and found no statistical difference.The results are presented in Figure S2.One amino acid substitution with a neutral effect was identified in the Cacna2d4 protein.The analysis demonstrates that the genes that are associated with morphogenesis, especially with neurogenesis, and cell-cell interactions (refer to Figure 4a) are significantly over-represented in the identified cluster of polymorphic genes harboring SNPs in the C3H mouse line compared to the C57Bl/6 mouse model.The significant overrepresentation of the genes annotated with GO terms associated with retina functioning, including "GO:0060041 retina development in camera-type eye", "GO:0060042 retina morphogenesis in camera-type eye", "GO:0031290 retinal ganglion cell axon guidance", "GO:0003407 neural retina development", and "GO:0060221 retinal rod cell differentiation", was also observed.
Rpgrip1 has 13 amino acid substitutions (listed in Table S2) and involves the deletions of 2 amino acid residues in positions 971 and 972.The significance of amino acid substitutions was estimated using web tools SIFT [12], PolyPhen2 [13], and PhD-SNP [14].All substitutions were identified as "neutral" with PhD-SNP and "tolerated" with SIFT.Only with PolyPhen2 were we able to characterize several polymorphisms in the gene Rpgrip1 as possibly damaging or probably damaging mutations (refer to Table S2).Taking into account that other tools did not reveal any damaging effect of these mutations, the results obtained with PolyPhen2 should be considered with caution, as this tool is characterized by a low specificity rate, with 15-20% of false-positive results [15,16].The deletion of two glutamic acid residues was localized on the interdomain spacer but could change the protein charge.We assessed the level of Rpgrip1 mRNA in rd1 and Pde6b-KO mice and found no statistical difference.The results are presented in Figure S2.One amino acid substitution with a neutral effect was identified in the Cacna2d4 protein.

Discussion
A number of hereditary retinal diseases such as RP lead to vision loss due to the degeneration of photoreceptor cells.The most extensively studied model of RP is the rd1 mouse model with a nonsense mutation in exon 7 of the pde6b gene.The Pde6b-KO mouse model presented in this paper differs from the rd1 model in the localization of a frameshift mutation in the pde6b gene.The homozygous Pde6b-KO mouse model demonstrated clinical signs of autosomal recessive retinitis pigmentosa due to a two-nucleotide deletion in the pde6b gene, leading to a frameshift and nonsense codon in the 104 position.It is assumed that the main cause of photoreceptors' death in rd1 and Pde6b-KO mice is the increased intracellular calcium influx through continuously opened cGMP-gated channels [17].Since the stop codon in the rd1 mutant mouse model is located in exon 7, unlike Pde6b-KO, one would expect that the shortened protein would accumulate.And yet, the results of Western blot analysis showed the absence of protein in both models.The functional state of the retina was monitored via in vivo electroretinogram (ERG) recordings.Histological analysis of eyecup cryosections at different time points of postnatal age was performed in order to track morphological changes.The analysis of the ERG revealed that the amplitude of the photoresponse in the Pde6b-KO mouse model at postnatal day 18 was higher than in the rd1 model.The amplitude of the responses continued to decrease until it became indistinguishable from the background noise at 24-27th day.The morphological analysis revealed that in both mouse models, the thickness of the photoreceptor layer at postnatal day 12 did not significantly differ from that of control animals but significantly reduced by day 15.It appears that the decrease in the number of photoreceptors in the rd1 model occurs considerably faster.This difference might arise from different genetic backgrounds because rd1 mice had a mixed origin (C3H/C57Bl/6), whereas Pde6b-KO was obtained on a pure C57Bl/6 background.Using bioinformatic analysis, we found more than 800,000 SNPs that distinguish the C3H and C57Bl/6 mouse lines, and some of them are associated with the development of RP.These data may indicate that the C3H background significantly contributes to the process of retinal degeneration at the early stages of RP development.Further research is needed to establish the exact cause of this phenomenon.

Production of Pde6b KO Mouse Model
All animal studies were approved by the Ethical Committee in the animal research field at Saint Petersburg State University (No 131-03-5, 11 October 2022).All procedures were conducted strictly in accordance with the Guide for the Care and Use of Laboratory Animals.All surgery was performed under general anesthesia using an intraperitoneal injection of TZ (30 mg/kg; Zoletil, Virbac Lab, Carros, France) and xylazine (6 mg/kg; Interchemie Werken de Adelaar, Venray, The Netherlands).
Standard methods of molecular biology and bioinformatics were used to generate CRISPR guide RNA constructs for the pde6b gene.The selection of the spacer sequence for the guide RNA was performed using the CHOPCHOP program "http://chopchop.cbu.uib.no/ (accessed on 3 November 2020)".In vitro transcription was performed using T7 RNA polymerase (Invitrogen, Carlsbad, CA, USA), and the gRNA sequence is listed in Table S1.RNA isolation and purification were performed using a MEGAclear™ Transcription Clean-Up Kit (Invitrogen, Carlsbad, CA, USA).
The injected zygotes were cultured in a KSOM medium (Sigma, St. Louis, MO, USA) at 37 • C under 5% CO 2 for around one hour.Afterward, approximately 8-12 successfully microinjected embryos were transferred into the oviducts of pseudopregnant CD1 mouse females with copulation plugs.In order to obtain CD1 pseudopregnant females (recipients), CD1 females were mated with vasectomized CD1 male mice the day before injection.On the 21st day after birth, F0 pups were genotyped via PCR followed by DNA sequencing analysis.Genomic DNA was extracted from mice's ears and analyzed with PCR using the primers listed in Table S1.Primers for amplification and sequencing were selected using the PerlPrimer program "http://perlprimer.sourceforge.net/(accessed on 2 December 2020)".

Electroretinography
To assess the functional state of the retina, the in vivo electroretinography method was used [19].This technique enables the measurement of changes in the potential on the cornea that occur as a result of retinal responses to light stimuli.The recordings started from the moment when the pup's eyes were fully opened (18th postnatal day, see [19]).To prepare the mice for ERG recordings, they were kept in darkness for 12 h prior to the procedure.All animal procedures were conducted under dim red lights to maintain darkness.During the recording, each mouse was anesthetized using a Minor Vet anesthesia apparatus (Zoomed, Moscow, Russia) with a gas mixture of isoflurane (2.5%, Laboratorios Karizoo, Spain) and oxygen.Additionally, the local anesthetic oxybuprocaine (inocaine 0.4% solution, Sentiss Rus, Moscow, Russia) and the pupil dilator tropicamide (midrimax 0.8% solution, Sentiss Russ, Moscow, Russia) were applied to both eyes of the mouse.For ERG recording, silver-embedded thread electrodes (Ocuscience, Henderson, NV, USA) were placed on the cornea of the left eye (stimulated by light) and the cornea of the right eye (kept in darkness under the black cap).The ground electrode was securely attached underneath the animal's skin on the back.The light stimulation system was equipped with a light-emitting diode that emitted light at a peak wavelength of λmax = 525 nm.The intensity of the light was controlled with a photodiode current and a set of gray neutral-density filters.The intensity of the LED was calibrated with an OPT-301 optosensor (Burr-Brown, Tucson, AZ, USA).Stimuli were presented as brief (10-50 ms) flashes in darkness.Since both rd1 and Pde6b-KO mice lacked rod-derived ERG components, only bright flashes (0.6 cd*s/m 2 and higher) were able to elicit a response.Signals from electrodes were passed through an amplifier (Model 3000, A-M Systems, Sequim, WA, USA) with a 3 kHz low-frequency analog filter and were recorded at a sampling frequency of 1 ms/point.For data acquisition and controlling the timing, duration, and intensity of the light stimulus during the experiments, the customwritten LabView 7.1 software (National Instruments, Austin, TX, USA) was employed.After baseline correction, the amplitudes of ERG b-waves were calculated and analyzed.All data are presented as the mean ± SEM.A Mann-Whitney U test with a p-value < 0.05 was used for statistical analysis.

Histology
Retinal morphology was compared between rd1 (3 males, 3 eyecups), Pde6b-KO (3 males, 3 eyecups), and age-matched wild-type mice (3 males, 3 eyecups) at each time point.At 12, 15, 18, and 21 days of postnatal age, the eyes were enucleated after mice were sacrificed via cervical dislocation.The enucleated whole eyes were prefixed for 5 min in 4% buffered paraformaldehyde, after which the corneas were cut out, and the lens and vitreous body were removed.The remaining eyecups were further fixed for at least 48 h in the same fixative at 4 • C.
Samples were then embedded in paraffin wax for sections (Histomix, Bio-Optica, Saint-Petersburg, Russia).Vertical sections of 6 µm thickness were cut, passing through or close to the optic disc.Specimens were stained using the standard hematoxylin and eosin (H&E) method and were captured using a light microscope fitted with a digital camera (Leica DM 6000, Leica Biosystems Nussloch GmbH, Nussloch, Germany).
For measuring the thickness of the outer nuclear layer (ONL), each hemisphere of the eyecup was divided into 3 equal segments from the optic nerve.At each segment, three independent measurements were performed and averaged.Measurements excluded the first 100 µm from the optic nerve due to the thinness of the retina at this location.At least three sections from each eye were counted, and the averages were used to represent the eye.
The morphometry of retinal layers was determined on the captured digital images of sections using a FIJI v1.53 (Fiji, RRID:SciRes_000137, SCR_002285) image processing package.All data are presented as the mean ± SEM.A Mann-Whitney U test with a p-value <0.05 was used for statistical analysis.

Total RNA Extraction and RT-PCR
At the 9th day of postnatal age, the eyes were enucleated under general anesthesia using an intraperitoneal injection of TZ (30 mg/kg; Zoletil, Virbac Lab, Carros, France) and xylazine (6 mg/kg; Interchemie Werken de Adelaar, Venray, The Netherlands).The corneas were cut out, and the lens and vitreous body were removed.The remaining eyecups were treated with RNAprotect reagent (Promega, Madison, WI, USA), and the total RNA was extracted with an RNeasy Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions.The RNA concentration and quality were measured after gel electrophoresis using the ChemiDoc XRS+ Gel Imaging System (Bio-Rad, Hercules, CA, USA) and Image Lab 6.0.1.software (Bio-Rad, Hercules, CA, USA).cDNA was synthesized using an iScript cDNA Synthesis Kit (Bio-Rad, Hercules, CA, USA), and the primers for RT-PCR are listed in Table S1.The 2 −∆∆Ct method was used to calculate the difference in the relative expression levels of the target genes normalized on the housekeeping gene GAPDH.

Protein Extraction and Western Blot Hybridization
On day 12 of postnatal age, the eyes were enucleated under general anesthesia using an intraperitoneal injection of TZ (30 mg/kg; Zoletil, Virbac Lab, Carros, France) and xylazine (6 mg/kg; Interchemie Werken de Adelaar, Venray, The Netherlands).The corneas were cut out, and the lens and vitreous body were removed.The remaining eyecups were homogenized using 0.5 µm glass beads (Sigma, St. Louis, MO, USA) in lysis buffer (25 mM Tris-HCl, pH 7.4, 100 mM NaCl, 1 mM dithiothreitol, 10 mM EDTA, and 2 mM PMSF) and a Complete™ protease inhibitor mixture (Roche Applied Science, Penzberg, Germany).The total protein concentration was measured using a Qubit 2.0 fluorometer and a protein assay kit (Invitrogen, Carlsbad, CA, USA), 35 µg of the total protein extracts was used per lane.The Pde6b protein was detected with the primary Anti-PDE6B Polyclonal antibody PA1-722 (Invitrogen, Carlsbad, CA, USA) and secondary HRP-conjugated Goat Anti-Rabbit IgG H&L ab205718 (Abcam, Cambridge, UK).The GAPDH protein was detected with the primary Anti-GAPDH polyclonal antibody AF7021 (Affinity Biosciences, Zhenjiang, China) and secondary HRP-conjugated Goat Anti-Rabbit IgG H&L ab205718 (Abcam, Cambridge, UK).Chemiluminescent detection was performed using an Amersham ECL Prime Western Blotting Detection Reagent (GE Healthcare, Chicago, IL, USA).

The Identification and Description of Genetic Polymorphisms
The C3H/HeJ mouse genome provided in the NCBI database (NCBI GenBankID GCA_921997125.2) was compared with the reference C57BL/6 genome (NCBI RefSeqID GCF_000001635.27) using the minimap2 tool (version 2.26) [20] with default settings.Postprocessing was performed using a custom Python script.

Function and Enrichment Analysis
Gene Ontology (GO) term enrichment analysis (i.e., the identification of GO terms [22] that are significantly enriched by the genes of the selected set) was performed in the gene cluster, which included the genes that harbor nucleotide substitutions in C3H mice compared to the C57Bl/6 mice genome.The enriched terms were identified and visualized using the clusterProfiler Bioconductor package [23].Biologic process (BP) ontology was used [22].
Then, the GO enrichment analysis was performed in the same way in the cluster, which includes the selected genes associated with retina development and functioning.Cell component (CC) and molecular function (MF) ontologies were used [22].

Functional Analysis of Genes Harboring Nucleotide Polymorphisms in C3H
The alignment of C3H and C57Bl/6 mouse genomes led to the identification of 221,946 nucleotide substitutions in 15,597 genes in C3H compared to the reference C57Bl/6 genome sequence.To explain the specific biological function of the genes in which polymorphisms were identified, we performed the biological process GO terms' enrichment analysis.

Patents
The patent for the proposed method for the generation and application of a mouse model with pde6b gene knockout was obtained with Patent Number RU 2 791 687 C1.

Figure 1 .
Figure 1.The mutations in the pde6b gene in both rd1 and Pde6b-KO mice lead to a decrease in corresponding mRNA levels and a complete absence of the mutant proteins: (a) the positions of mutations in rd1 and Pde6b-KO mice; the sequence of sgRNA spacer is underlined, and amino acid substitution and nucleotide substitutions are shown in green; (b) the level of pde6b mRNA in retinas of wt, rd1, and Pde6b-KO mice on the 9th day; (c) the amount of Pde6b protein in mouse retina; the size of Pde6b wt protein is 99 kDa, and the size of Pde6b rd1 and Pde6bKO hypothetical proteins is 40 and 11 kDa, respectively.

Figure 1 .
Figure 1.The mutations in the pde6b gene in both rd1 and Pde6b-KO mice lead to a decrease in corresponding mRNA levels and a complete absence of the mutant proteins: (a) the positions of mutations in rd1 and Pde6b-KO mice; the sequence of sgRNA spacer is underlined, and amino acid substitution and nucleotide substitutions are shown in green; (b) the level of pde6b mRNA in retinas of wt, rd1, and Pde6b-KO mice on the 9th day; (c) the amount of Pde6b protein in mouse retina; the size of Pde6b wt protein is 99 kDa, and the size of Pde6b rd1 and Pde6bKO hypothetical proteins is 40 and 11 kDa, respectively.

Figure 2 .
Figure 2. Pde6b-KO and rd1 mice demonstrate different levels of photoresponse amplitude in the electrophysiological analysis: (a) comparison between the amplitudes of b-wave ERGs of C57Bl/6 (wild type, wt), Pde6b-KO and rd1 mice on the 18th day of life (n = 4-6); the Mann-Whitney test revealed statistically significant differences (* p < 0.05); (b) comparison of the response shapes inPde6b-KO and rd1 mice on the 18th day of the postnatal period demonstrated reduced photoresponse in both lines compared to the wild type.The response was measured using a stimulating 10 ms flash with an intensity of 31.6 cd*s/m 2 ; (c) comparison of b-wave amplitudes in Pde6b-KO and rd1 mice (n = 6) on days 18, 21, and 24 of the postnatal period revealed that the Pde6b-KO photoresponse was higher than that of the rd1 model.The Mann-Whitney test revealed statistically significant differences (* p < 0.05).

Figure 2 .
Figure 2. Pde6b-KO and rd1 mice demonstrate different levels of photoresponse amplitude in the electrophysiological analysis: (a) comparison between the amplitudes of b-wave ERGs of C57Bl/6 (wild type, wt), Pde6b-KO and rd1 mice on the 18th day of life (n = 4-6); the Mann-Whitney test revealed statistically significant differences (* p < 0.05); (b) comparison of the response shapes in Pde6b-KO and rd1 mice on the 18th day of the postnatal period demonstrated reduced photoresponse in both lines compared to the wild type.The response was measured using a stimulating 10 ms flash with an intensity of 31.6 cd*s/m 2 ; (c) comparison of b-wave amplitudes in Pde6b-KO and rd1 mice (n = 6) on days 18, 21, and 24 of the postnatal period revealed that the Pde6b-KO photoresponse was higher than that of the rd1 model.The Mann-Whitney test revealed statistically significant differences (* p < 0.05).

Figure 3 .
Figure 3. Retinal morphology and morphometry of retinal thickness: (a) H&E staining of sections from the central part of the retina of wild-type mice, as well as rd1 and Pde6b-KO mice at 12, 15, 18, and 21 days of postnatal age.ONL, outer nuclear layer (photoreceptor cell bodies); INL, inner nuclear layer; GCL, ganglion cell layer; H&E, hematoxylin and eosin; pd, postnatal day.Scale bars: 50 μm.(b) The thickness of the retinal outer nuclear layer of wild-type mice, and rd1 and Pde6b-KO mice at 12, 15, 18, and 21 days of postnatal age.OD, optic disk; L1-L3, R1-R3, measured segments of left and right hemispheres of the eyecup, respectively.Normally, the center is thicker than the periphery, but in models on the 15th day, the thickness of the central parts of the retina sharply decreases.

Figure 3 .
Figure 3. Retinal morphology and morphometry of retinal thickness: (a) H&E staining of sections from the central part of the retina of wild-type mice, as well as rd1 and Pde6b-KO mice at 12, 15, 18, and 21 days of postnatal age.ONL, outer nuclear layer (photoreceptor cell bodies); INL, inner nuclear layer; GCL, ganglion cell layer; H&E, hematoxylin and eosin; pd, postnatal day.Scale bars: 50 µm.(b) The thickness of the retinal outer nuclear layer of wild-type mice, and rd1 and Pde6b-KO mice at 12, 15, 18, and 21 days of postnatal age.OD, optic disk; L1-L3, R1-R3, measured segments of left and right hemispheres of the eyecup, respectively.Normally, the center is thicker than the periphery, but in models on the 15th day, the thickness of the central parts of the retina sharply decreases.

Figure 4 .
Figure 4. Functional analysis of genes that harbor polymorphisms in C3H mice: (a) the dot plot representing the Gene Ontology (GO) biological process terms' enrichment in the cluster of genes for which nucleotide substitutions in C3H mice were identified (top-20 terms are represented); (b) the dot plot representing the Gene Ontology (GO) cell component terms' enrichment in the cluster of eye development and function-related genes for which nucleotide substitutions in C3H mice were identified (top-10 terms are presented); (c) the dot plot representing the Gene Ontology (GO)

Figure 4 .
Figure 4. Functional analysis of genes that harbor polymorphisms in C3H mice: (a) the dot plot representing the Gene Ontology (GO) biological process terms' enrichment in the cluster of genes for which nucleotide substitutions in C3H mice were identified (top-20 terms are represented); (b) the dot plot representing the Gene Ontology (GO) cell component terms' enrichment in the cluster of eye development and function-related genes for which nucleotide substitutions in C3H mice were identified (top-10 terms are presented); (c) the dot plot representing the Gene Ontology (GO) molecular function terms' enrichment in the cluster of eye development and function-related genes for which nucleotide substitutions in C3H mice were identified (top-10 terms are presented).
Institutional Review Board Statement: The animal study protocol was approved by the Ethics Committee in the field of animal research at Saint Petersburg State University (Protocol Code 131-03-5, 11 October 2022).Not applicable.All data are available upon reasonable request.