Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Dominant Distal Myopathy 3 (MPD3) Caused by a Deletion in the HNRNPA1 Gene

Background and Objectives

To determine the genetic cause of the disease in the previously reported family with adult-onset autosomal dominant distal myopathy (myopathy, distal, 3; MPD3).

Methods

Continued clinical evaluation including muscle MRI and muscle pathology. A linkage analysis with single nucleotide polymorphism arrays and genome sequencing were used to identify the genetic defect, which was verified by Sanger sequencing. RNA sequencing was used to investigate the transcriptional effects of the identified genetic defect.

Results

Small hand muscles (intrinsic, thenar, and hypothenar) were first involved with spread to the lower legs and later proximal muscles. Dystrophic changes with rimmed vacuoles and cytoplasmic inclusions were observed in muscle biopsies at advanced stage. A single nucleotide polymorphism array confirmed the previous microsatellite-based linkage to 8p22-q11 and 12q13-q22. Genome sequencing of three affected family members combined with structural variant calling revealed a small heterozygous deletion of 160 base pairs spanning the second last exon 10 of the heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1) gene, which is in the linked region on chromosome 12. Segregation of the mutation with the disease was confirmed by Sanger sequencing. RNA sequencing showed that the mutant allele produces a shorter mutant mRNA transcript compared with the wild-type allele. Immunofluorescence studies on muscle biopsies revealed small p62 and larger TDP-43 inclusions.

Discussion

A small exon 10 deletion in the gene HNRNPA1 was identified as the cause of MPD3 in this family. The new HNRNPA1-related phenotype, upper limb presenting distal myopathy, was thus confirmed, and the family displays the complexities of gene identification.

Missense variant in LOXHD1 is associated with canine nonsyndromic hearing loss

Hearing loss is a common sensory deficit in both humans and dogs. In canines, the genetic basis is largely unknown, as genetic variants have only been identified for a syndromic form of hearing impairment. We observed a congenital or early-onset sensorineural hearing loss in a Rottweiler litter. Assuming an autosomal recessive inheritance, we used a combined approach of homozygosity mapping and genome sequencing to dissect the genetic background of the disorder. We identified a fully segregating missense variant in LOXHD1, a gene that is known to be essential for cochlear hair cell function and associated with nonsyndromic hearing loss in humans and mice. The canine LOXHD1 variant was specific to the Rottweiler breed in our study cohorts of pure-bred dogs. However, it also was present in some mixed-breed dogs, of which the majority showed Rottweiler ancestry. Low allele frequencies in these populations, 2.6% and 0.04%, indicate a rare variant. To summarize, our study describes the first genetic variant for canine nonsyndromic hearing loss, which is clinically and genetically similar to human LOXHD1-related hearing disorder, and therefore, provides a new large animal model for hearing loss. Equally important, the affected breed will benefit from a genetic test to eradicate this LOXHD1-related hearing disorder from the population.

In-frame deletion in canine PITRM1 is associated with a severe early-onset epilepsy, mitochondrial dysfunction and neurodegeneration

We investigated the clinical, genetic, and pathological characteristics of a previously unknown severe juvenile brain disorder in several litters of Parson Russel Terriers. The disease started with epileptic seizures at 6–12 weeks of age and progressed rapidly to status epilepticus and death or euthanasia. Histopathological changes at autopsy were restricted to the brain. There was severe acute neuronal degeneration and necrosis diffusely affecting the grey matter throughout the brain with extensive intraneuronal mitochondrial crowding and accumulation of amyloid-β (Aβ). Combined homozygosity mapping and genome sequencing revealed an in-frame 6-bp deletion in the nuclear-encoded pitrilysin metallopeptidase 1 (PITRM1) encoding for a mitochondrial protease involved in mitochondrial targeting sequence processing and degradation. The 6-bp deletion results in the loss of two amino acid residues in the N-terminal part of PITRM1, potentially affecting protein folding and function. Assessment of the mitochondrial function in the affected brain tissue showed a significant deficiency in respiratory chain function. The functional consequences of the mutation were modeled in yeast and showed impaired growth in permissive conditions and an impaired respiration capacity. Loss-of-function variants in human PITRM1 result in a childhood-onset progressive amyloidotic neurological syndrome characterized by spinocerebellar ataxia with behavioral, psychiatric and cognitive abnormalities. Homozygous Pitrm1-knockout mice are embryonic lethal, while heterozygotes show a progressive, neurodegenerative phenotype characterized by impairment in motor coordination and Aβ deposits. Our study describes a novel early-onset PITRM1-related neurodegenerative canine brain disorder with mitochondrial dysfunction, Aβ accumulation, and lethal epilepsy. The findings highlight the essential role of PITRM1 in neuronal survival and strengthen the connection between mitochondrial dysfunction and neurodegeneration.

A missense variant in IFT122 associated with a canine model of retinitis pigmentosa

Retinitis pigmentosa (RP) is a blinding eye disease affecting nearly two million people worldwide. Dogs are affected with a similar illness termed progressive retinal atrophy (PRA). Lapponian herders (LHs) are affected with several types of inherited retinal dystrophies, and variants in PRCD and BEST1 genes have been associated with generalized PRA and canine multifocal retinopathy 3 (cmr3), respectively. However, all retinal dystrophy cases in LHs are not explained by these variants, indicating additional genetic causes of disease in the breed. We collected DNA samples from 10 PRA affected LHs, with known PRCD and BEST1 variants excluded, and 34 unaffected LHs. A genome-wide association study identified a locus on CFA20 (p_raw = 2.4 × 10^-7, p_Bonf = 0.035), and subsequent whole-genome sequencing of an affected LH revealed a missense variant, c.3176G>A, in the intraflagellar transport 122 (IFT122) gene. The variant was also found in Finnish Lapphunds, in which its clinical relevancy needs to be studied further. The variant interrupts a highly conserved residue, p.(R1059H), in IFT122 and likely impairs its function. Variants in IFT122 have not been associated with retinal degeneration in mammals, but the loss of ift122 in zebrafish larvae impaired opsin transport and resulted in progressive photoreceptor degeneration. Our study establishes a new spontaneous dog model to study the role of IFT122 in RP biology, while the affected breed will benefit from a genetic test for a recessive condition.

Intronic variant in POU1F1 associated with canine pituitary dwarfism

The anterior pituitary gland secretes several endocrine hormones, essential for growth, reproduction and other basic physiological functions. Abnormal development or function of the pituitary gland leads to isolated or combined pituitary hormone deficiency (CPHD). At least 30 genes have been associated with human CPHD, including many transcription factors, such as POU1F1. CPHD occurs spontaneously also in mice and dogs. Two affected breeds have been reported in dogs: German Shepherds with a splice defect in the LHX3 gene and Karelian Bear Dogs (KBD) with an unknown genetic cause. We obtained samples from five KBDs presenting dwarfism and abnormal coats. A combined analysis of genome-wide association and next-generation sequencing mapped the disease to a region in chromosome 31 and identified a homozygous intronic variant in the fourth exon of the POU1F1 gene in the affected dogs. The identified variant, c.605-3C>A, resided in the splice region and was predicted to affect splicing. The variant's screening in three new prospective cases, related breeds, and ~ 8000 dogs from 207 breeds indicated complete segregation in KBDs with a carrier frequency of 8%, and high breed-specificity as carriers were found at a low frequency only in Lapponian Herders, a related breed. Our study establishes a novel canine model for CPHD with a candidate POU1F1 defect.

Improved Criteria for the Classification of Titin Variants in Inherited Skeletal Myopathies

BACKGROUND

Extensive genetic screening results in the identification of thousands of rare variants that are difficult to interpret. Because of its sheer size, rare variants in the titin gene (TTN) are detected frequently in any individual. Unambiguous interpretation of molecular findings is almost impossible in many patients with myopathies or cardiomyopathies.

OBJECTIVE

To refine the current classification framework for TTN-associated skeletal muscle disorders and standardize the interpretation of TTN variants.

METHODS

We used the guidelines issued by the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) to re-analyze TTN genetic findings from our patient cohort.

RESULTS

We identified in the classification guidelines three rules that are not applicable to titin-related skeletal muscle disorders; six rules that require disease-/gene-specific adjustments and four rules requiring quantitative thresholds for a proper use. In three cases, the rule strength need to be modified.

CONCLUSIONS

We suggest adjustments are made to the guidelines. We provide frequency thresholds to facilitate filtering of candidate causative variants and guidance for the use and interpretation of functional data and co-segregation evidence. We expect that the variant classification framework for TTN-related skeletal muscle disorders will be further improved along with a better understanding of these diseases.

webGQT: A Shiny Server for Genotype Query Tools for Model-Based Variant Filtering

Summary

Genotype Query Tools (GQT) were developed to discover disease-causing variations from billions of genotypes and millions of genomes, processes data at substantially higher speed over other existing methods. While GQT has been available to a wide audience as command-line software, the difficulty of constructing queries among non-IT or non-bioinformatics researchers has limited its applicability. To overcome this limitation, we developed webGQT, an easy-to-use tool with a graphical user interface. With pre-built queries across three modules, webGQT allows for pedigree analysis, case-control studies, and population frequency studies. As a package, webGQT allows researchers with less or no applied bioinformatics/IT experience to mine potential disease-causing variants from billions.

Results

webGQT offers a flexible and easy-to-use interface for model-based candidate variant filtering for Mendelian diseases from thousands to millions of genomes at a reduced computation time. Additionally, webGQT provides adjustable parameters to reduce false positives and rescue missing genotypes across all modules. Using a case study, we demonstrate the applicability of webGQT to query non-human genomes. In addition, we demonstrate the scalability of webGQT on large data sets by implementing complex population-specific queries on the 1000 Genomes Project Phase 3 data set, which includes 8.4 billion variants from 2504 individuals across 26 different populations. Furthermore, webGQT supports filtering single-nucleotide variants, short insertions/deletions, copy number or any other variant genotypes supported by the VCF specification. Our results show that webGQT can be used as an online web service, or deployed on personal computers or local servers within research groups.

Availability

webGQT is made available to the users in three forms: 1) as a webserver available at https://vm1138.kaj.pouta.csc.fi/webgqt/, 2) as an R package to install on personal computers, and 3) as part of the same R package to configure on the user's own servers. The application is available for installation at https://github.com/arumds/webgqt.

Recessive missense LAMP3 variant associated with defect in lamellar body biogenesis and fatal neonatal interstitial lung disease in dogs

Neonatal interstitial lung diseases due to abnormal surfactant biogenesis are rare in humans and have never been reported as a spontaneous disorder in animals. We describe here a novel lung disorder in Airedale Terrier (AT) dogs with clinical symptoms and pathology similar to the most severe neonatal forms of human surfactant deficiency. Lethal hypoxic respiratory distress and failure occurred within the first days or weeks of life in the affected puppies. Transmission electron microscopy of the affected lungs revealed maturation arrest in the formation of lamellar bodies (LBs) in the alveolar epithelial type II (AECII) cells. The secretory organelles were small and contained fewer lamellae, often in combination with small vesicles surrounded by an occasionally disrupted common limiting membrane. A combined approach of genome-wide association study and whole exome sequencing identified a recessive variant, c.1159G>A, p.(E387K), in LAMP3, a limiting membrane protein of the cytoplasmic surfactant organelles in AECII cells. The substitution resides in the LAMP domain adjacent to a conserved disulfide bond. In summary, this study describes a novel interstitial lung disease in dogs, identifies a new candidate gene for human surfactant dysfunction and brings important insights into the essential role of LAMP3 in the process of the LB formation.

Novel protective and risk loci in hip dysplasia in German Shepherds

Canine hip dysplasia is a common, non-congenital, complex and hereditary disorder. It can inflict severe pain via secondary osteoarthritis and lead to euthanasia. An analogous disorder exists in humans. The genetic background of hip dysplasia in both species has remained ambiguous despite rigorous studies. We aimed to investigate the genetic causes of this disorder in one of the high-risk breeds, the German Shepherd. We performed genetic analyses with carefully phenotyped case-control cohorts comprising 525 German Shepherds. In our genome-wide association studies we identified four suggestive loci on chromosomes 1 and 9. Targeted resequencing of the two loci on chromosome 9 from 24 affected and 24 control German Shepherds revealed deletions of variable sizes in a putative enhancer element of the NOG gene. NOG encodes for noggin, a well-described bone morphogenetic protein inhibitor affecting multiple developmental processes, including joint development. The deletion was associated with the healthy controls and mildly dysplastic dogs suggesting a protective role against canine hip dysplasia. Two enhancer variants displayed a decreased activity in a dual luciferase reporter assay. Our study identifies novel loci and candidate genes for canine hip dysplasia, with potential regulatory variants in the NOG gene. Further research is warranted to elucidate how the identified variants affect the expression of noggin in canine hips, and what the potential effects of the other identified loci are.

A novel KRT71 variant in curly-coated dogs

Curly fur is a common phenotype in many dog breeds, known to result from a missense variant (c.451C>T) in exon 2 of the keratin 71 (KRT71) gene. During screening for this variant across various breeds, we found that Curly Coated Retrievers (CCRs) fixed with the trait did not carry the known variant. By analysis of whole-genome sequencing data of one CCR we identified a novel genetic cause for curly fur. We found a novel structural variant in exon 7 of the KRT71 gene (c.1266_1273delinsACA) that was predicted to result in a frameshift and stop loss, therefore significantly affecting the structure of the protein, if translated. The variant was also found at lower frequencies in five other breeds, including Lagotto Romagnolo, Bichon Frise, Spanish Water Dog, Chesapeake Bay Retriever and Irish Terrier. One curly-coated Lagotto carried neither of the two KRT71 variants. These results identify a second variant for curly coat in KRT71 and suggest the existence of additional alleles. This study enables the development of an additional KRT71 gene test for breeders to understand and manage coat types.

Author Correction: Assembly and Analysis of Unmapped Genome Sequence Reads Reveal Novel Sequence and Variation in Dogs

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Copy number variation analysis increases the diagnostic yield in muscle diseases

Objective

Copy number variants (CNVs) were analyzed from next-generation sequencing data, with the aim of improving diagnostic yield in skeletal muscle disorder cases.

Methods

Four publicly available bioinformatic analytic tools were used to analyze CNVs from sequencing data from patients with muscle diseases. The patients were previously analyzed with a targeted gene panel for single nucleotide variants and small insertions and deletions, without achieving final diagnosis. Variants detected by multiple CNV analysis tools were verified with either array comparative genomic hybridization or PCR. The clinical significance of the verified CNVs was interpreted, considering previously identified variants, segregation studies, and clinical information of the patient cases.

Results

Combining analysis of all different mutation types enabled integration of results and identified the final cause of the disease in 9 myopathy cases. Complex effects like compound heterozygosity of different mutation types and compound disease arising from variants of different genes were unraveled. We identified the first large intragenic deletion of the titin (TTN) gene implicated in the pathogenesis of a severe form of myopathy. Our work also revealed a "double-trouble" effect in a patient carrying a single heterozygous insertion/deletion mutation in the TTN gene and a Becker muscular dystrophy causing deletion in the dystrophin gene.

Conclusions

Causative CNVs were identified proving that analysis of CNVs is essential for increasing the diagnostic yield in muscle diseases. Complex severe muscular dystrophy phenotypes can be the result of different mutation types but also of the compound effect of 2 different genetic diseases.

Nonsense variant in COL7A1 causes recessive dystrophic epidermolysis bullosa in Central Asian Shepherd dogs

A rare hereditary mechanobullous disorder called epidermolysis bullosa (EB) causes blistering in the skin and the mucosal membranes. To date, nineteen EB-related genes have been discovered in human and other species. We describe here a novel EB variant in dogs. Two newborn littermates of Central Asian Shepherd dogs with severe signs of skin blistering were brought to a veterinary clinic and euthanized due to poor prognosis. In post-mortem examination, the puppies were shown to have findings in the skin and the mucosal membranes characteristic of EB. A whole-genome sequencing of one of the affected puppies was performed to identify the genetic cause. The resequencing data were filtered under a recessive model against variants from 31 other dog genomes, revealing a homozygous case-specific nonsense variant in one of the known EB-causing genes, COL7A1 (c.4579C>T, p.R1527*). The variant results in a premature stop codon and likely absence of the functional protein in the basement membrane of the skin in the affected dogs. This was confirmed by immunohistochemistry using a COL7A1 antibody. Additional screening of the variant indicated full penetrance and breed specificity at ~28% carrier frequency. In summary, this study reveals a novel COL7A1 variant causing recessive dystrophic EB and provides a genetic test for the eradication of the disease from the breed.

Association study reveals novel risk loci for sporadic inclusion body myositis

BACKGROUND AND PURPOSE

The aim was to identify potential genetic risk factors associated with sporadic inclusion body myositis (sIBM).

METHODS

An association based case-control approach was utilized on whole exome sequencing data of 30 Finnish sIBM patients and a control cohort (n = 193). A separate Italian cohort of sIBM patients (n = 12) was used for evaluation of the results.

RESULTS

Seven single nucleotide polymorphisms were identified in five genes that have a considerably higher observed frequency in Finnish sIBM patients compared to the control population, and the previous association of the genetic human leukocyte antigen region was confirmed.

CONCLUSIONS

All seven identified variants could individually or in combination increase the susceptibility for sIBM.

Molecular Characterization of Three Canine Models of Human Rare Bone Diseases: Caffey, van den Ende-Gupta, and Raine Syndromes

One to two percent of all children are born with a developmental disorder requiring pediatric hospital admissions. For many such syndromes, the molecular pathogenesis remains poorly characterized. Parallel developmental disorders in other species could provide complementary models for human rare diseases by uncovering new candidate genes, improving the understanding of the molecular mechanisms and opening possibilities for therapeutic trials. We performed various experiments, e.g. combined genome-wide association and next generation sequencing, to investigate the clinico-pathological features and genetic causes of three developmental syndromes in dogs, including craniomandibular osteopathy (CMO), a previously undescribed skeletal syndrome, and dental hypomineralization, for which we identified pathogenic variants in the canine SLC37A2 (truncating splicing enhancer variant), SCARF2 (truncating 2-bp deletion) and FAM20C (missense variant) genes, respectively. CMO is a clinical equivalent to an infantile cortical hyperostosis (Caffey disease), for which SLC37A2 is a new candidate gene. SLC37A2 is a poorly characterized member of a glucose-phosphate transporter family without previous disease associations. It is expressed in many tissues, including cells of the macrophage lineage, e.g. osteoclasts, and suggests a disease mechanism, in which an impaired glucose homeostasis in osteoclasts compromises their function in the developing bone, leading to hyperostosis. Mutations in SCARF2 and FAM20C have been associated with the human van den Ende-Gupta and Raine syndromes that include numerous features similar to the affected dogs. Given the growing interest in the molecular characterization and treatment of human rare diseases, our study presents three novel physiologically relevant models for further research and therapy approaches, while providing the molecular identity for the canine conditions.

Rare developmental disorders make a major contribution to pediatric hospital admissions and mortality. There is a growing interest in the development of therapeutics for these conditions, but that requires understanding of the genetic cause and pathology. Research can be facilitated by physiologically relevant models, such as dogs with corresponding disorders. We have characterized the clinical features and genetic causes of three developmental syndromes in dogs, including craniomandibular osteopathy (CMO), a previously undescribed skeletal syndrome, and dental hypomineralization, for which we identified mutations in the canine SLC37A2, SCARF2 and FAM20C genes, respectively. CMO is a clinical equivalent to an infantile cortical hyperostosis (Caffey disease) for which SLC37A2 is a new candidate gene. SLC37A2 is a glucose-phosphate transporter in osteoclasts, and its defect suggests an impaired glucose homeostasis in developing bone, leading to hyperostosis. Mutations in the SCARF2 and FAM20C genes have been associated with the human van den Ende-Gupta and Raine syndromes. Our study provides molecular identity for the canine conditions and presents three novel physiologically relevant models of human rare diseases.

Identification of a common risk haplotype for canine idiopathic epilepsy in the ADAM23 gene

Background

Idiopathic epilepsy is a common neurological disease in human and domestic dogs but relatively few risk genes have been identified to date. The seizure characteristics, including focal and generalised seizures, are similar between the two species, with gene discovery facilitated by the reduced genetic heterogeneity of purebred dogs. We have recently identified a risk locus for idiopathic epilepsy in the Belgian Shepherd breed on a 4.4 megabase region on CFA37.

Results

We have expanded a previous study replicating the association with a combined analysis of 157 cases and 179 controls in three additional breeds: Schipperke, Finnish Spitz and Beagle (p_c = 2.9e–07, p_GWAS = 1.74E-02). A targeted resequencing of the 4.4 megabase region in twelve Belgian Shepherd cases and twelve controls with opposite haplotypes identified 37 case-specific variants within the ADAM23 gene. Twenty-seven variants were validated in 285 cases and 355 controls from four breeds, resulting in a strong replication of the ADAM23 locus (p_raw = 2.76e–15) and the identification of a common 28 kb-risk haplotype in all four breeds. Risk haplotype was present in frequencies of 0.49–0.7 in the breeds, suggesting that ADAM23 is a low penetrance risk gene for canine epilepsy.

Conclusions

These results implicate ADAM23 in common canine idiopathic epilepsy, although the causative variant remains yet to be identified. ADAM23 plays a role in synaptic transmission and interacts with known epilepsy genes, LGI1 and LGI2, and should be considered as a candidate gene for human epilepsies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1651-9) contains supplementary material, which is available to authorized users.

Increased expression of MERTK is associated with a unique form of canine retinopathy

Progressive retinal degenerations are among the most common causes of blindness both in human and in dogs. Canine progressive retinal atrophy (PRA) resembles human retinitis pigmentosa (RP) and is typically characterized by a progressive loss of rod photoreceptors followed by a loss of cone function. The disease gradually progress from the loss of night and day vision to a complete blindness. We have recently described a unique form of retinopathy characterized by the multifocal gray/brown discoloration and thinning of the retina in the Swedish Vallhund (SV) breed. We aimed to identify the genetic cause by performing a genome wide association analysis in a cohort of 18 affected and 10 healthy control dogs using Illumina's canine 22k SNP array. We mapped the disease to canine chromosome 17 (p = 7.7×10⁻⁵) and found a 6.1 Mb shared homozygous region in the affected dogs. A combined analysis of the GWAS and replication data with additional 60 dogs confirmed the association (p = 4.3×10⁻⁸, OR = 11.2 for homozygosity). A targeted resequencing of the entire associated region in four cases and four controls with opposite risk haplotypes identified several variants in the coding region of functional candidate genes, such as a known retinopathy gene, MERTK. However, none of the identified coding variants followed a compelling case- or breed-specific segregation pattern. The expression analyses of four candidate genes in the region, MERTK, NPHP1, ANAPC1 and KRCC1, revealed specific upregulation of MERTK in the retina of the affected dogs. Collectively, these results indicate that the retinopathy is associated with overexpression of MERTK, however further investigation is needed to discover the regulatory mutation for the better understanding of the disease pathogenesis. Our study establishes a novel gain-of-function model for the MERTK biology and provides a therapy model for retinopathy MERTK inhibitors. Meanwhile, a marker-based genetic counseling can be developed to revise breeding programs.