Understanding and treating human diseases require valid animal models. Leveraging the genetic diversity in rhesus macaque populations across eight primate centers in the United States, we conduct targeted-sequencing on 1845 individuals for 374 genes linked to inherited human retinal and neurodevelopmental diseases. We identify over 47,000 single nucleotide variants, a substantial proportion of which are shared with human populations. By combining rhesus and human allele frequencies with established variant prediction methods, we develop a machine learning-based score that outperforms established methods in predicting missense variant pathogenicity. Remarkably, we find a marked number of loss-of-function variants and putative deleterious variants, which may lead to the development of rhesus disease models. Through phenotyping of macaques carrying a pathogenic OPA1:p.A8S variant, we identify a genetic model of autosomal dominant optic atrophy. Finally, we present a public website housing variant and genotype data from over two thousand rhesus macaques.

The vast majority of protein-coding genes in the human genome produce multiple mRNA isoforms through alternative splicing, significantly enhancing the complexity of the transcriptome and proteome. To establish an efficient method for characterizing transcript isoforms within tissue samples, we conducted a systematic comparison between single-cell long-read and conventional short-read RNA sequencing techniques. The transcriptome of approximately 30,000 mouse retina cells was profiled using 1.54 billion Illumina short reads and 1.40 billion Oxford Nanopore Technologies long reads. Consequently, we identify 44,325 transcript isoforms, with a notable 38% previously uncharacterized and 17% expressed exclusively in distinct cellular subclasses. We observe that long-read sequencing not only matches the gene expression and cell-type annotation performance of short-read sequencing but also excel in the precise identification of transcript isoforms. While transcript isoforms are often shared across various cell types, their relative abundance shows considerable cell type–specific variation. The data generated from our study significantly enhance the existing repertoire of transcript isoforms, thereby establishing a resource for future research into the mechanisms and implications of alternative splicing within retinal biology and its links to related diseases.

Inherited retinal diseases (IRDs) are a diverse set of visual disorders that collectively represent a major cause of early-onset blindness. With the reduction in sequencing costs in recent years, whole-genome sequencing (WGS) is being used more frequently, particularly when targeted gene panels and whole-exome sequencing (WES) fail to detect pathogenic mutations in patients. In this study, we performed mutation screens using WGS for a cohort of 311 IRD patients whose mutations were undetermined. A total of nine putative pathogenic mutations in six IRD patients were identified, including six novel mutations. Among them, four were deep intronic mutations that affected mRNA splicing, while the other five affected protein-coding sequences. Our results suggested that the rate of resolution of unsolved cases via targeted gene panels and WES can be further enhanced with WGS; however, the overall improvement may be limited.

Despite increasing success in determining genetic diagnosis for patients with inherited retinal diseases (IRDs), mutations in about 30% of the IRD cases remain unclear or unsettled after targeted gene panel or whole exome sequencing. In this study, we aimed to investigate the contributions of structural variants (SVs) to settling the molecular diagnosis of IRD with whole-genome sequencing (WGS). A cohort of 755 IRD patients whose pathogenic mutations remain undefined were subjected to WGS. Four SV calling algorithms including include MANTA, DELLY, LUMPY and CNVnator were used to detect SVs throughout the genome. All SVs identified by any one of these four algorithms were included for further analysis. AnnotSV was used to annotate these SVs. SVs that overlap with known IRD-associated genes were examined with sequencing coverage, junction reads and discordant read pairs. Polymerase Chain Reaction (PCR) followed by Sanger sequencing was used to further confirm the SVs and identify the breakpoints. Segregation of the candidate pathogenic alleles with the disease was performed when possible. A total of 16 candidate pathogenic SVs were identified in 16 families, including deletions and inversions, representing 2.1% of patients with previously unsolved IRDs. Autosomal dominant, autosomal recessive and X-linked inheritance of disease-causing SVs were observed in 12 different genes. Among these, SVs in CLN3, EYS and PRPF31 were found in multiple families. Our study suggests that the contribution of SVs detected by short-read WGS is about 0.25% of our IRD patient cohort and is significantly lower than that of single nucleotide changes and small insertions and deletions.

The vast majority of protein-coding genes in the human genome produce multiple mRNA isoforms through alternative splicing, significantly enhancing the complexity of the transcriptome and proteome. To establish an efficient method for characterizing transcript isoforms within tissue samples, we conducted a systematic comparison between single-cell long-read and conventional short-read RNA sequencing techniques. The transcriptome of approximately 30,000 mouse retina cells was profiled using 1.54 billion Illumina short reads and 1.40 billion Oxford Nanopore long reads. Consequently, we identified 44,325 transcript isoforms, with a notable 38% previously uncharacterized and 17% expressed exclusively in distinct cellular subclasses. We observed that long-read sequencing not only matched the gene expression and cell-type annotation performance of short-read sequencing but also excelled in the precise identification of transcript isoforms. While transcript isoforms are often shared across various cell types, their relative abundance shows considerable cell-type-specific variation. The data generated from our study significantly enhance the existing repertoire of transcript isoforms, thereby establishing a foundational resource for future research into the mechanisms and implications of alternative splicing within retinal biology and its links to related diseases.

Background: Leber congenital amaurosis (LCA) is a congenital onset severe form of inherited retinal dystrophy (IRD) and a common cause of pediatric blindness. Disease-causing variants in at least 14 genes are reported to predispose LCA phenotype. LCA is inherited as an autosomal recessive disease. It can be an isolated eye disorder or as part of a syndrome, such as Senior Loken or Joubert syndrome. Sequencing studies from consanguineous populations have proven useful for novel variants identification; thus, the present study aimed to explore the genetic heterogeneity of 15 consanguineous Pakistani families, each segregating a severe IRD phenotype using targeted next generation sequencing. Methods: This study enrolled 15 consanguineous families, each with multiple affected cases of retinal dystrophy phenotype. DNA was extracted from blood samples. Targeted panel sequencing of 344 known genes for IRDs was performed, followed by Sanger sequencing for segregation analysis. Results: Data analysis revealed a total of eight reported (c.316C>T and c.506G>A in RDH12; c.864dup and c.1012C>T in SPATA7, as well as c.1459T>C, c.1062_1068del, c.1495+1G>A, c.998G>A in the CRB1, LCA5, TULP1, and IFT140 genes, respectively) and four novel homozygous (c.720+1G>T in LCA5, c.196G>C in LRAT, c.620_625del in PRPH2, and c.3411_3414del in CRB1) variants segregating with disease phenotype in each respective family. Furthermore, a novel heterozygous variant of CRB1 gene, i.e., c.1935delC in compound heterozygous condition was found segregating with disease phenotype in one large family with multiple consanguinity loops. Conclusion: Comprehensive molecular diagnosis of 15 consanguineous Pakistani families led to the identification of a total of 5 novel variants contributing to genetic heterogeneity of LCA-associated genes and helped to provide genetic counseling to the affected families.