Glaucoma Genetic Risk Scores in the Million Veteran Program.

PURPOSE: Primary open-angle glaucoma (POAG) is a degenerative eye disease for which early treatment is critical to mitigate visual impairment and irreversible blindness. POAG-associated loci individually confer incremental risk. Genetic risk score(s) (GRS) could enable POAG risk stratification. Despite significantly higher POAG burden among individuals of African ancestry (AFR), GRS are limited in this population. A recent large-scale, multi-ancestry meta-analysis identified 127 POAG-associated loci and calculated cross-ancestry and ancestry-specific effect estimates, including in European ancestry (EUR) and AFR individuals. We assessed the utility of the 127-variant GRS for POAG risk stratification in EUR and AFR Veterans in the Million Veteran Program (MVP). We also explored the association between GRS and documented invasive glaucoma surgery (IGS). DESIGN: Cross-sectional study. PARTICIPANTS: MVP Veterans with imputed genetic data, including 5830 POAG cases (445 with IGS documented in the electronic health record) and 64 476 controls. METHODS: We tested unweighted and weighted GRS of 127 published risk variants in EUR (3382 cases and 58 811 controls) and AFR (2448 cases and 5665 controls) Veterans in the MVP. Weighted GRS were calculated using effect estimates from the most recently published report of cross-ancestry and ancestry-specific meta-analyses. We also evaluated GRS in POAG cases with documented IGS. MAIN OUTCOME MEASURES: Performance of 127-variant GRS in EUR and AFR Veterans for POAG risk stratification and association with documented IGS. RESULTS: GRS were significantly associated with POAG (P < 5 × 10-5) in both groups; a higher proportion of EUR compared with AFR were consistently categorized in the top GRS decile (21.9%-23.6% and 12.9%-14.5%, respectively). Only GRS weighted by ancestry-specific effect estimates were associated with IGS documentation in AFR cases; all GRS types were associated with IGS in EUR cases. CONCLUSIONS: Varied performance of the GRS for POAG risk stratification and documented IGS association in EUR and AFR Veterans highlights (1) the complex risk architecture of POAG, (2) the importance of diverse representation in genomics studies that inform GRS construction and evaluation, and (3) the necessity of expanding diverse POAG-related genomic data so that GRS can equitably aid in screening individuals at high risk of POAG and who may require more aggressive treatment.

A Discovery with Potential to Revitalize Hammerhead Ribozyme Therapeutics for Treatment of Inherited Retinal Degenerations.

Hammerhead ribozymes (hhRzs), RNA enzymes capable of site-specific cleavage of arbitrary target mRNAs, have faced significant hurdles in development and optimization as gene therapeutics for clinical translation. Chemical and biological barriers must be overcome to realize an effective therapeutic. A new Facilitated ribozyme has been identified with greatly enhanced kinetic properties that lead new insight on the capacity of ribozymes to target mutant genes to treat inherited retinal degenerations.

A cellular high-throughput screening approach for therapeutic trans-cleaving ribozymes and RNAi against arbitrary mRNA disease targets.

Major bottlenecks in development of therapeutic post transcriptional gene silencing (PTGS) agents (e.g. ribozymes, RNA interference, antisense) include the challenge of mapping rare accessible regions of the mRNA target that are open for annealing and cleavage, testing and optimization of agents in human cells to identify lead agents, testing for cellular toxicity, and preclinical evaluation in appropriate animal models of disease. Methods for rapid and reliable cellular testing of PTGS agents are needed to identify potent lead candidates for optimization. Our goal was to develop a means of rapid assessment of many RNA agents to identify a lead candidate for a given mRNA associated with a disease state. We developed a rapid human cell-based screening platform to test efficacy of hammerhead ribozyme (hhRz) or RNA interference (RNAi) constructs, using a model retinal degeneration target, human rod opsin (RHO) mRNA. The focus is on RNA Drug Discovery for diverse retinal degeneration targets. To validate the approach, candidate hhRzs were tested against NUH↓ cleavage sites (N = G,C,A,U; H = C,A,U) within the target mRNA of secreted alkaline phosphatase (SEAP), a model gene expression reporter, based upon in silico predictions of mRNA accessibility. HhRzs were embedded in a larger stable adenoviral VAI RNA scaffold for high cellular expression, cytoplasmic trafficking, and stability. Most hhRz expression plasmids exerted statistically significant knockdown of extracellular SEAP enzyme activity when readily assayed by a fluorescence enzyme assay intended for high throughput screening (HTS). Kinetics of PTGS knockdown of cellular targets is measureable in live cells with the SEAP reporter. The validated SEAP HTS platform was transposed to identify lead PTGS agents against a model hereditary retinal degeneration target, RHO mRNA. Two approaches were used to physically fuse the model retinal gene target mRNA to the SEAP reporter mRNA. The most expedient way to evaluate a large set of potential VAI-hhRz expression plasmids against diverse NUH↓ cleavage sites uses cultured human HEK293S cells stably expressing a dicistronic Target-IRES-SEAP target fusion mRNA. Broad utility of this rational RNA drug discovery approach is feasible for any ophthalmological disease-relevant mRNA targets and any disease mRNA targets in general. The approach will permit rank ordering of PTGS agents based on potency to identify a lead therapeutic compound for further optimization.

RTN4IP1-associated non-syndromic optic neuropathy and rod-cone dystrophy.

BACKGROUND: Biallelic variants in RTN4IP1 are a well-established cause of syndromic and nonsyndromic early-onset autosomal recessive optic neuropathy. They have more recently been reported to cause a concomitant but later-onset rod-cone dystrophy with or without syndromic features. METHODS: A comprehensive evaluation was performed that included assessment of visual and retinal function, clinical examination, and retinal imaging. Childhood ophthalmic records as well as the results of genetic testing were evaluated. RESULTS: A 24-year-old female described longstanding reduced visual acuity with more recent subjective impairment of dark adaptation. Visual acuity was subnormal in both eyes. Goldmann kinetic perimetry demonstrated scotomas in a pattern consistent with the presence of both optic neuropathy and rod-cone dystrophy with fundus exam as well as retinal imaging showing corroborating findings. Full-field electroretinography further confirmed the presence of a rod-cone dystrophy. Genetic testing demonstrated biallelic variants in RTN4IP1, one of which was novel, in association with the ocular findings. CONCLUSIONS: RTN4IP1-associated early-onset bilateral optic neuropathy with rod-cone dystrophy is a recently described clinical entity with limited reports available to-date. The present case provides additional support for this dual phenotype and identifies a novel causative variant.

Pharmacokinetics of Monoclonal Antibody and Antibody Fragments in The Mouse Eye Following Intravitreal Administration.

Mice are rarely used in pharmacokinetic (PK) studies of ocular therapeutics due to the small size of their eyes and challenges in drug administration, tissue collection, and analysis of drug concentrations. Therefore, ocular PK of protein therapeutics in mouse eye following intravitreal (IVT) administration is not known. Here, we have presented the first of its kind investigation, to study the PK of 4 different size non-binding protein therapeutics in mouse plasma, cornea/ICB, vitreous humor, retina, and posterior cup (including choroid) following IVT administration. Administered proteins include trastuzumab (150 kDa) and F(ab)2 (100 kDa), Fab, and scFv (27 kDa) fragments of trastuzumab. An imaging and injection apparatus suitable for performing small (50 nL) IVT injections in mice was developed, and techniques for enucleation of the eye and dissection of ocular tissues were developed. Furthermore, a sensitive enzyme-linked immunosorbent assay (ELISA) for detection of proteins in very small amounts of ocular tissues were developed. It was observed that elimination from the vitreous chamber was the primary driver of PK in the cornea/ICB, retina, posterior cup, and plasma. Trastuzumab displays first-order kinetics in the vitreous humor with a half-life of 18.8 h. F(ab)2, Fab, and ScFv show biphasic PK profiles with distribution phases becoming more rapid as molecular weight decreases, and terminal elimination becoming longer as molecular weight decreases, with terminal half-lives of 16.3, 20.6, and 48.9 h, respectively. The mean residence times of trastuzumab, F(ab)2, Fab, and scFv in the vitreous humor were 26.0, 12.2, 10.7, and 8.16 h, respectively. It was found that the mean residence time in vitreous humor doubles with an increase in molecular weight of ∼69 kDa. Interestingly, the PK of proteins measured in the un-injected eye suggest the presence of a pathway for drug transfer between the eyes, which needs to be further validated. Overall, the findings presented here pave the way for drug discovery and development studies of protein therapeutics for ophthalmic indications in mice.

Diversity is key for cross-ancestry transferability of glaucoma genetic risk scores in Hispanic Veterans in the Million Veteran Program.

A major goal of precision medicine is to stratify patients based on their genetic risk for a disease to inform future screening and intervention strategies. For conditions like primary open-angle glaucoma (POAG), the genetic risk architecture is complicated with multiple variants contributing small effects on risk. Following the tepid success of genome-wide association studies for high-effect disease risk variant discovery, genetic risk scores (GRS), which collate effects from multiple genetic variants into a single measure, have shown promise for disease risk stratification. We assessed the application of GRS for POAG risk stratification in Hispanic-descent (HIS) and European-descent (EUR) Veterans in the Million Veteran Program. Unweighted and cross-ancestry meta-weighted GRS were calculated based on 127 genomic variants identified in the most recent report of cross-ancestry POAG meta-analyses. We found that both GRS types were associated with POAG case-control status and performed similarly in HIS and EUR Veterans. This trend was also seen in our subset analysis of HIS Veterans with less than 50% EUR global genetic ancestry. Our findings highlight the importance of evaluating GRS based on known POAG risk variants in different ancestry groups and emphasize the need for more multi-ancestry POAG genetic studies.

Association Between Fuchs Endothelial Corneal Dystrophy, Diabetes Mellitus, and Multimorbidity.

PURPOSE: The aim of this study was to assess risk for demographic variables and other health conditions that are associated with Fuchs endothelial corneal dystrophy (FECD). METHODS: We developed a FECD case-control algorithm based on structured electronic health record data and confirmed accuracy by individual review of charts at 3 Veterans Affairs (VA) Medical Centers. This algorithm was applied to the Department of VA Million Veteran Program cohort from whom sex, genetic ancestry, comorbidities, diagnostic phecodes, and laboratory values were extracted. Single-variable and multiple variable logistic regression models were used to determine the association of these risk factors with FECD diagnosis. RESULTS: Being a FECD case was associated with female sex, European genetic ancestry, and a greater number of comorbidities. Of 1417 diagnostic phecodes evaluated, 213 had a significant association with FECD, falling in both ocular and nonocular conditions, including diabetes mellitus (DM). Five of 69 laboratory values were associated with FECD, with the direction of change for 4 being consistent with DM. Insulin dependency and type 1 DM raised risk to a greater degree than type 2 DM, like other microvascular diabetic complications. CONCLUSIONS: Female sex, European ancestry, and multimorbidity increased FECD risk. Endocrine/metabolic clinic encounter codes and altered patterns of laboratory values support DM increasing FECD risk. Our results evoke a threshold model in which the FECD phenotype is intensified by DM and potentially other health conditions that alter corneal physiology. Further studies to better understand the relationship between FECD and DM are indicated and may help identify opportunities for slowing FECD progression.

Fine-Scale Adaptations to Environmental Variation and Growth Strategies Drive Phyllosphere Methylobacterium Diversity.

Methylobacterium is a prevalent bacterial genus of the phyllosphere. Despite its ubiquity, little is known about the extent to which its diversity reflects neutral processes like migration and drift, versus environmental filtering of life history strategies and adaptations. In two temperate forests, we investigated how phylogenetic diversity within Methylobacterium is structured by biogeography, seasonality, and growth strategies. Using deep, culture-independent barcoded marker gene sequencing coupled with culture-based approaches, we uncovered a considerable diversity of Methylobacterium in the phyllosphere. We cultured different subsets of Methylobacterium lineages depending upon the temperature of isolation and growth (20°C or 30°C), suggesting long-term adaptation to temperature. To a lesser extent than temperature adaptation, Methylobacterium diversity was also structured across large (>100 km; between forests) and small (<1.2 km; within forests) geographical scales, among host tree species, and was dynamic over seasons. By measuring the growth of 79 isolates during different temperature treatments, we observed contrasting growth performances, with strong lineage- and season-dependent variations in growth strategies. Finally, we documented a progressive replacement of lineages with a high-yield growth strategy typical of cooperative, structured communities in favor of those characterized by rapid growth, resulting in convergence and homogenization of community structure at the end of the growing season. Together, our results show how Methylobacterium is phylogenetically structured into lineages with distinct growth strategies, which helps explain their differential abundance across regions, host tree species, and time. This work paves the way for further investigation of adaptive strategies and traits within a ubiquitous phyllosphere genus. IMPORTANCE Methylobacterium is a bacterial group tied to plants. Despite the ubiquity of methylobacteria and the importance to their hosts, little is known about the processes driving Methylobacterium community dynamics. By combining traditional culture-dependent and -independent (metabarcoding) approaches, we monitored Methylobacterium diversity in two temperate forests over a growing season. On the surface of tree leaves, we discovered remarkably diverse and dynamic Methylobacterium communities over short temporal (from June to October) and spatial (within 1.2 km) scales. Because we cultured different subsets of Methylobacterium diversity depending on the temperature of incubation, we suspected that these dynamics partly reflected climatic adaptation. By culturing strains under laboratory conditions mimicking seasonal variations, we found that diversity and environmental variations were indeed good predictors of Methylobacterium growth performances. Our findings suggest that Methylobacterium community dynamics at the surface of tree leaves results from the succession of strains with contrasting growth strategies in response to environmental variations.

Comprehensive Phylogenomics of Methylobacterium Reveals Four Evolutionary Distinct Groups and Underappreciated Phyllosphere Diversity.

Methylobacterium is a group of methylotrophic microbes associated with soil, fresh water, and particularly the phyllosphere, the aerial part of plants that has been well studied in terms of physiology but whose evolutionary history and taxonomy are unclear. Recent work has suggested that Methylobacterium is much more diverse than thought previously, questioning its status as an ecologically and phylogenetically coherent taxonomic genus. However, taxonomic and evolutionary studies of Methylobacterium have mostly been restricted to model species, often isolated from habitats other than the phyllosphere and have yet to utilize comprehensive phylogenomic methods to examine gene trees, gene content, or synteny. By analyzing 189 Methylobacterium genomes from a wide range of habitats, including the phyllosphere, we inferred a robust phylogenetic tree while explicitly accounting for the impact of horizontal gene transfer (HGT). We showed that Methylobacterium contains four evolutionarily distinct groups of bacteria (namely A, B, C, D), characterized by different genome size, GC content, gene content, and genome architecture, revealing the dynamic nature of Methylobacterium genomes. In addition to recovering 59 described species, we identified 45 candidate species, mostly phyllosphere-associated, stressing the significance of plants as a reservoir of Methylobacterium diversity. We inferred an ancient transition from a free-living lifestyle to association with plant roots in Methylobacteriaceae ancestor, followed by phyllosphere association of three of the major groups (A, B, D), whose early branching in Methylobacterium history has been heavily obscured by HGT. Together, our work lays the foundations for a thorough redefinition of Methylobacterium taxonomy, beginning with the abandonment of Methylorubrum.

Development and Evaluation of a Rules-based Algorithm for Primary Open-Angle Glaucoma in the VA Million Veteran Program.

The availability of electronic health record (EHR)-linked biobank data for research presents opportunities to better understand complex ocular diseases. Developing accurate computable phenotypes for ocular diseases for which gold standard diagnosis includes imaging remains inaccessible in most biobank-linked EHRs. The objective of this study was to develop and validate a computable phenotype to identify primary open-angle glaucoma (POAG) through accessing the Department of Veterans Affairs (VA) Computerized Patient Record System (CPRS) and Million Veteran Program (MVP) biobank. Accessing CPRS clinical ophthalmology data from VA Medical Center Eye Clinic (VAMCEC) patients, we developed and iteratively refined POAG case and control algorithms based on clinical, prescription, and structured diagnosis data (ICD-CM codes). Refinement was performed via detailed chart review, initially at a single VAMCEC (n = 200) and validated at two additional VAMCECs (n = 100 each). Positive and negative predictive values (PPV, NPV) were computed as the proportion of CPRS patients correctly classified with POAG or without POAG, respectively, by the algorithms, validated by ophthalmologists and optometrists with access to gold-standard clinical diagnosis data. The final algorithms performed better than previously reported approaches in assuring the accuracy and reproducibility of POAG classification (PPV >83% and NPV >97%) with consistent performance in Black or African American and in White Veterans. Applied to the MVP to identify cases and controls, genetic analysis of a known POAG-associated locus further validated the algorithms. We conclude that ours is a viable approach to use combined EHR-genetic data to study patients with complex diseases that require imaging confirmation.

Rapid, cell-based toxicity screen of potentially therapeutic post-transcriptional gene silencing agents.

Post-transcriptional gene silencing (PTGS) agents such as antisense, ribozymes and RNA interference (RNAi) have great potential as therapeutics for a variety of eye diseases including retinal and macular degenerations, glaucoma, corneal degenerations, inflammatory and viral conditions. Despite their great potential and over thirty years of academic and corporate research only a single PTGS agent is currently approved for human therapy for a single disease. Substantial challenges exist to achieving both efficacious and safe PTGS agents. Efficacy, as measured in specific target mRNA and protein knockdown, depends upon a number of complex factors including the identification of rare regions of target mRNA accessibility, cellular co-localization of the PTGS agent in sufficient concentration with the target mRNA, and stability of the PTGS agent in the target cells in which it is delivered or expressed. Safety is commonly measured by lack of cytotoxicity or other deleterious cellular responses in cells in which the PTGS agent is delivered or expressed. To relieve major bottlenecks in RNA drug discovery novel, efficient, inexpensive, and rapid tools are needed to facilitate lead identification of the most efficacious PTGS agent, rational optimization of efficacy of the lead agent, and lead agent safety determinations. We have developed a technological platform using cell culture expression systems that permits lead identification and efficacy optimization of PTGS agents against arbitrary disease target mRNAs under relatively high throughput conditions. Here, we extend the technology platform to include PTGS safety determinations in cultured human cells that are expected to represent the common cellular housekeeping microenvironment. We developed a high throughput screening (HTS) cytotoxicity assay in 96-well plate format based around the SYTOX Green dye which is excluded from healthy viable cells and becomes substantially fluorescent only after entering cells and binding to nuclear DNA. In this format we can test a number of PTGS agents for cellular toxicity relative to control elements. We also developed an HTS 96-well plate assay that allows us to assess the impact of any given PTGS agent on stimulating a variety of common cellular stress signaling pathways (e.g. CRE, SRE, AP-1, NFκB, Myc, and NFAT) that could indicate possible deleterious effects of PTGS agents either dependent or independent of base pairing complementarity with target mRNAs. To this end we exploited the secreted alkaline phosphatase (SEAP) Pathway Profiling System where the expression of the secreted reporter protein is coupled to transcriptional activation of a variety of promoter elements involved in common cell signaling pathways. We found that a variety of lead hammerhead ribozyme (hhRz) and short hairpin (shRNA) expression constructs did not exert cytotoxicity in human cells when driven by highly active RNA Pol-III promoters. We also found that most of the cell signaling pathways tested (CRE, SRE, Myc, and NFAT) did not significantly couple through upregulation to expression of the set of PTGS agents tested. AP-1 and NFκB upregulation both appear to couple to the expression of some PTGS agents which likely reflect the known properties of these pathways to be stimulated by abundant small structured RNAs.

No evidence for phylosymbiosis in western chipmunk species.

Phylosymbiosis refers to a congruent pattern between the similarity of microbiomes of different species and the branching pattern of the host phylogeny. Phylosymbiosis has been detected in a variety of vertebrate and invertebrate hosts, but has only been assessed in geographically isolated populations. We tested for phylosymbiosis in eight (sub)species of western chipmunks with overlapping ranges and ecological niches; we used a nuclear (Acrosin) and a mitochondrial (CYTB) phylogenetic marker because there are many instances of mitochondrial introgression in chipmunks. We predicted that similarity among microbiomes increases with: (1) increasing host mitochondrial relatedness, (2) increasing host nuclear genome relatedness and (3) decreasing geographic distance among hosts. We did not find statistical evidence supporting phylosymbiosis in western chipmunks. Furthermore, in contrast to studies of other mammalian microbiomes, similarity of chipmunk microbiomes is not predominantly determined by host species. Sampling site explained most variation in microbiome composition, indicating an important role of local environment in shaping microbiomes. Fecal microbiomes of chipmunks were dominated by Bacteroidetes (72.2%), followed by Firmicutes (24.5%), which is one of the highest abundances of Bacteroidetes detected in wild mammals. Future work will need to elucidate the effects of habitat, ecology and host genomics on chipmunk microbiomes.

Development of lead hammerhead ribozyme candidates against human rod opsin mRNA for retinal degeneration therapy.

To identify lead candidate allele-independent hammerhead ribozymes (hhRz) for the treatment of autosomal dominant mutations in the human rod opsin (RHO) gene, we tested a series of hhRzs for potential to significantly knockdown human RHO gene expression in a human cell expression system. Multiple computational criteria were used to select target mRNA regions likely to be single stranded and accessible to hhRz annealing and cleavage. Target regions are tested for accessibility in a human cell culture expression system where the hhRz RNA and target mRNA and protein are coexpressed. The hhRz RNA is embedded in an adenoviral VAI RNA chimeric RNA of established structure and properties which are critical to the experimental paradigm. The chimeric hhRz-VAI RNA is abundantly transcribed so that the hhRzs are expected to be in great excess over substrate mRNA. HhRz-VAI traffics predominantly to the cytoplasm to colocalize with the RHO mRNA target. Colocalization is essential for second-order annealing reactions. The VAI chimera protects the hhRz RNA from degradation and provides for a long half-life. With cell lines chosen for high transfection efficiency and a molar excess of hhRz plasmid over target plasmid, the conditions of this experimental paradigm are specifically designed to evaluate for regions of accessibility of the target mRNA in cellulo. Western analysis was used to measure the impact of hhRz expression on RHO protein expression. Three lead candidate hhRz designs were identified that significantly knockdown target protein expression relative to control (p<0.05). Successful lead candidates (hhRz CUC [see in text downward arrow] 266, hhRz CUC [see in text downward arrow] 1411, hhRz AUA [see in text downward arrow] 1414) targeted regions of human RHO mRNA that were predicted to be accessible by a bioinformatics approach, whereas regions predicted to be inaccessible supported no knockdown. The maximum opsin protein level knockdown is approximately 30% over a 48h paradigm of testing. These results validate a rigorous computational bioinformatics approach to detect accessible regions of target mRNAs in cellulo. The opsin knockdown effect could prove to be clinically significant when integrated over longer periods in photoreceptors. Further optimization and animal testing are the next step in this stratified RNA drug discovery program. A recently developed novel and efficient screening assay based upon expression of a dicistronic mRNA (RHO-IRES-SEAP) containing both RHO and reporter (SEAP) cDNAs was used to compare the hhRz 266 lead candidate to another agent (Rz525/hhRz485) already known to partially rescue retinal degeneration in a rodent model. Lead hhRz 266 CUC [see in text downward arrow] proved more efficacious than Rz525/hhRz485 which infers viability for rescue of retinal degeneration in appropriate preclinical models of disease.

Systematic Screening, Rational Development, and Initial Optimization of Efficacious RNA Silencing Agents for Human Rod Opsin Therapeutics.

PURPOSE: To systematically evaluate human rod opsin (hRHO) mRNA for potential target sites sensitive to posttranscriptional gene silencing (PTGS) by hammerhead ribozyme (hhRz) or RNA interference (RNAi) in human cells. To develop a comprehensive strategy to identify and optimize lead candidate agents for PTGS gene therapeutics. METHODS: In multidisciplinary RNA drug discovery, computational mRNA accessibility and in vitro experimental methods using reverse transcription-polymerase chain reaction (RT-PCR) were used to map accessibility in full-length hRHO transcripts. HhRzs targeted predicted accessible and inaccessible sites and were screened for cellular knockdown using a bicistronic reporter construct. Lead hhRz and RNAi PTGS agents were rationally optimized for target knockdown in human cells. RESULTS: Systematic screening of hRHO mRNA targeting agents resulted in lead candidate identification of a novel hhRz embedded in an RNA scaffold. Rational optimization strategies identified a minimal 725 hhRz as the most active agent. Recently identified tertiary accessory elements did not enhance activity. A 725-short-hairpin RNA (shRNA) agent exerts log-order knockdown. Silent modulation of the 725-hhRz target site in hRHO mRNA resulted in resistance to knockdown. CONCLUSIONS: Combining rational RNA drug design with cell-based screening allowed rapid identification of lead agents targeting hRHO. Optimization strategies identified the agent with highest intracellular activity. These agents have therapeutic potential in a mutation-independent strategy for adRP, or other degenerations where hRHO is a target. This approach can be broadly applied to any validated target mRNA, regardless of the disease. TRANSLATIONAL RELEVANCE: This work establishes a platform approach to develop RNA biologicals for the treatment of human disease.

Phylogeography, genetic diversity, and connectivity of brown bear populations in Central Asia.

Knowledge of genetic diversity and population structure is critical for conservation and management planning at the population level within a species' range. Many brown bear populations in Central Asia are small and geographically isolated, yet their phylogeographic relationships, genetic diversity, and contemporary connectivity are poorly understood. To address this knowledge gap, we collected brown bear samples from the Gobi Desert (n = 2360), Altai, Sayan, Khentii, and Ikh Khyangan mountains of Mongolia (n = 79), and Deosai National Park in the Himalayan Mountain Range of Pakistan (n = 5) and generated 927 base pairs of mitochondrial DNA (mtDNA) sequence data and genotypes at 13 nuclear DNA microsatellite loci. We documented high levels of mtDNA and nDNA diversity in the brown bear populations of northern Mongolia (Altai, Sayan, Buteeliin nuruu and Khentii), but substantially lower diversity in brown bear populations in the Gobi Desert and Himalayas of Pakistan. We detected 3 brown bear mtDNA phylogeographic groups among bears of the region, with clade 3a1 in Sayan, Khentii, and Buteeliin nuruu mountains, clade 3b in Altai, Sayan, Buteeliin nuruu, Khentii, and Ikh Khyangan, and clade 6 in Gobi and Pakistan. Our results also clarified the phylogenetic relationships and divergence times with other brown bear mtDNA clades around the world. The nDNA genetic structure analyses revealed distinctiveness of Gobi bears and different population subdivisions compared to mtDNA results. For example, genetic distance for nDNA microsatellite loci between the bears in Gobi and Altai (FST = 0.147) was less than that of the Gobi and Pakistan (FST = 0.308) suggesting more recent male-mediated nuclear gene flow between Gobi and Altai than between Gobi and the Pakistan bears. Our results provide valuable information for conservation and management of bears in this understudied region of Central Asia and highlight the need for special protection and additional research on Gobi brown bears.

Genetically-guided algorithm development and sample size optimization for age-related macular degeneration cases and controls in electronic health records from the VA Million Veteran Program.

Electronic health records (EHRs) linked to extensive biorepositories and supplemented with lifestyle, behavioral, and environmental exposure data, have enormous potential to contribute to genomic discovery, a necessary step in the pathway towards translational or precision medicine. A major bottleneck in incorporating EHRs into genomic studies is the extraction of research-grade variables for analysis, particularly when gold-standard measurements are not available or accessible. Here we develop algorithms for age-related macular degeneration (AMD), a common cause of blindness among the elderly, and controls free of AMD. These computable phenotypes were developed using billing codes (ICD-9-CM and ICD-10-CM) and Current Procedural Terminology (CPT) codes and evaluated in two study sites of the Veterans Affairs Million Veteran Program: Louis Stokes Cleveland VA Medical Center and the Providence VA Medical Center. After establishing a high overall positive and negative predictive values (93% and 95%, respectively) through manual chart review, the candidate algorithm was deployed in the full VA MVP dataset of >500,000 participants. The algorithm was then optimized in a data cube using a variety of approaches including adjusting inclusion age thresholds by examining previously-reported genetic associations for CFH (rs10801555, a proxy for rs1061170) and ARMS2 (rs10490924). The algorithm with the smallest p-values for the known genetic associations was selected for downstream and on-going AMD genomic discovery efforts. This two-phase approach to developing research-grade case/control variables for AMD genomic studies capitalizes on established genetic associations resulting in high precision and optimized sample sizes, an approach that can be applied to other large-scale biobanks linked to EHRs for precision medicine research.

Bottlenecks in development of retinal therapeutic post-transcriptional gene silencing agents.

Development of post-transcriptional gene silencing (PTGS) agents for therapeutic purposes is an immense challenge in modern biology. Established technologies used to knockdown a specific target RNA and its cognate protein: antisense, ribozyme, RNAi, all conditionally depend upon an initial, critical annealing event of the PTGS ligand to a target RNA. In this review we address the nature of the bottlenecks, emphasizing the biocomplexity of target RNA structure, that currently limit PTGS therapeutic development. We briefly review existing and emerging technologies designed to release these constraints to realize the potential of PTGS agents in gene based therapies.

Ultrahigh Resolution Mouse Optical Coherence Tomography to Aid Intraocular Injection in Retinal Gene Therapy Research.

HR-SD-OCT is utilized to monitor the progression of photoreceptor degeneration in live mouse models, assess the delivery of therapeutic agents into the subretinal space, and to evaluate toxicity and efficacy in vivo. HR-SD-OCT uses near infrared light (800-880 nm) and has optics specifically designed for the unique optics of the mouse eye with sub-2-micron axial resolution. Transgenic mouse models of outer retinal (photoreceptor) degeneration and controls were imaged to assess the disease progression. Pulled glass microneedles were used to deliver sub retinal injections of adeno-associated virus (AAV) or nanoparticles (NP) via a trans-scleral and trans-choroidal approach. Careful positioning of the needle into the subretinal space was required prior to a calibrated pressure injection, which delivers fluid into the sub retinal space. Real time subretinal surgery was conducted on our retinal imaging system (RIS). HR-SD-OCT demonstrated progressive uniform retinal degeneration due to expression of a toxic mutant human mutant rhodopsin (P347S) (RHOP347S) transgene in mice. HR-SD-OCT allows rigorous quantification of all the retinal layers. Outer nuclear layer (ONL) thickness and photoreceptor outer segment length (OSL) measurements correlate with photoreceptor vitality, degeneration, or rescue. The RIS delivery system allows real-time visualization of subretinal injections in neonatal (~P10-14) or adult mice, and HR-SD-OCT immediately determines success of delivery and maps areal extent. HR-SD-OCT is a powerful tool that can evaluate the success of subretinal surgery in mice, in addition to measuring vitality of photoreceptors in vivo. HR-SD-OCT can also be used to identify uniform animal cohorts to evaluate the extent of retinal degeneration, toxicity, and therapeutic rescue in preclinical gene therapy research studies.

Effects of Pathogenic Variations in the Human Rhodopsin Gene (hRHO) on the Predicted Accessibility for a Lead Candidate Ribozyme.

Purpose: The mutation-independent strategy for hammerhead ribozyme (hhRz) or RNA interference (RNAi)-based gene therapeutics to treat autosomal dominant diseases is predicated on the hypothesis that a single therapeutic would equivalently suppress all/most of the diverse mutant mRNAs in patients with the disease phenotype. However, the hypothesis has not been formally tested. We address this through a comprehensive bioinformatics study of how mutations affect target mRNA structure accessibility for a single lead hhRz therapeutic (725GUC↓), designed against human rod rhodopsin mRNA (hRHO), for patients with hRHO mutations that cause autosomal dominant retinitis pigmentosa. Methods: A total of 199 in silico coding region mutations (missense, nonsense, insert, deletion, indel) were made in hRHO mRNA based on Human Gene Mutation Database and Database of Single Nucleotide Polymorphisms. Each mRNA was folded with MFold, SFold, and OligoWalk algorithms and subjected to a bioinformatics model called multiparameter prediction of RNA accessibility. Predicted accessibility of each mutant over both a broad local region and the explicit lead ribozyme annealing site were compared quantitatively to wild-type hRHO mRNA. Results: Accessibility of the 725GUC↓ site is sensitive to some mutations. For single nucleotide missense mutations, proximity of the mutation to the hhRz annealing site increases the impact on predicted accessibility, but some distant mutations also influence accessibility. Conclusions: A mutation-independent strategy appears viable in this specific context but certain mutations could significantly influence ribozyme or RNAi efficacy through impact on accessibility at the target annealing site/region. This possibility must be considered in applications of this gene therapy strategy.