Correction: Genetic variation affects morphological retinal phenotypes extracted from UK Biobank optical coherence tomography images.

[This corrects the article DOI: 10.1371/journal.pgen.1009497.].

Genetic variation affects morphological retinal phenotypes extracted from UK Biobank optical coherence tomography images.

Optical Coherence Tomography (OCT) enables non-invasive imaging of the retina and is used to diagnose and manage ophthalmic diseases including glaucoma. We present the first large-scale genome-wide association study of inner retinal morphology using phenotypes derived from OCT images of 31,434 UK Biobank participants. We identify 46 loci associated with thickness of the retinal nerve fibre layer or ganglion cell inner plexiform layer. Only one of these loci has been associated with glaucoma, and despite its clear role as a biomarker for the disease, Mendelian randomisation does not support inner retinal thickness being on the same genetic causal pathway as glaucoma. We extracted overall retinal thickness at the fovea, representative of foveal hypoplasia, with which three of the 46 SNPs were associated. We additionally associate these three loci with visual acuity. In contrast to the Mendelian causes of severe foveal hypoplasia, our results suggest a spectrum of foveal hypoplasia, in part genetically determined, with consequences on visual function.

Transient MicroRNA Expression Enhances Myogenic Potential of Mouse Embryonic Stem Cells.

MicroRNAs (miRNAs) are known regulators of various cellular processes, including pluripotency and differentiation of embryonic stem cells (ESCs). We analyzed differentiation of two ESC lines-D3 and B8, and observed significant differences in the expression of miRNAs and genes involved in pluripotency and differentiation. We also examined if transient miRNA overexpression could serve as a sufficient impulse modulating differentiation of mouse ESCs. ESCs were transfected with miRNA Mimics and differentiated in embryoid bodies and embryoid body outgrowths. miRNAs involved in differentiation of mesodermal lineages, such as miR145 and miR181, as well as miRNAs regulating myogenesis (MyomiRs)-miR1, miR133a, miR133b, and miR206 were tested. Using such approach, we proved that transient overexpression of molecules selected by us modulated differentiation of mouse ESCs. Increase in miR145 levels upregulated Pax3, Pax7, Myod1, Myog, and MyHC2, while miR181 triggered the expression of such crucial myogenic factors as Myf5 and MyHC2. As a result, the ability of ESCs to initiate myogenic differentiation and form myotubes was enhanced. Premature expression of MyomiRs had, however, an adverse effect on myogenic differentiation of ESCs. Stem Cells 2018;36:655-670.

Role of lysophosphatidic acid in the retinal pigment epithelium and photoreceptors.

The human retina is a complex structure of organised layers of specialised cells that support the transmission of light signals to the visual cortex. The outermost layer of the retina, the retinal pigment epithelium (RPE), forms part of the blood retina barrier and is implicated in many retinal diseases. Lysophosphatidic acid (LPA) is a bioactive lipid exerting pleiotropic effects in various cell types, during development, normal physiology and disease. Its producing enzyme, AUTOTAXIN (ATX), is highly expressed by the pigmented epithelia of the human eye, including the RPE. Using human pluripotent stem cell (hPSC)-derived retinal cells, we interrogated the role of LPA in the human RPE and photoreceptors. hPSC-derived RPE cells express and synthesize functional ATX, which is predominantly secreted apically of the RPE, suggesting it acts in a paracrine manner to regulate photoreceptor function. In RPE cells, LPA regulates tight junctions, in a receptor-dependent mechanism, with an increase in OCCLUDIN and ZONULA OCCLUDENS (ZO)-1 expression at the cell membrane, accompanied by an increase in the transepithelial resistance of the epithelium. High concentration of LPA decreases phagocytosis of photoreceptor outer segments by the RPE. In hPSC-derived photoreceptors, LPA induces morphological rearrangements by modulating the actin myosin cytoskeleton, as evidenced by Myosin Light Chain l membrane relocation. Collectively, our data suggests an important role of LPA in the integrity and functionality of the healthy retina and blood retina barrier.

Organelle mapping in dendrites of human iPSC-derived neurons reveals dynamic functional dendritic Golgi structures.

Secretory pathways within dendrites of neurons have been proposed for local transport of newly synthesized proteins. However, little is known about the dynamics of the local secretory system and whether the organelles are transient or stable structures. Here, we quantify the spatial and dynamic behavior of dendritic Golgi and endosomes during differentiation of human neurons generated from induced pluripotent stem cells (iPSCs). In early neuronal development, before and during migration, the entire Golgi apparatus transiently translocates from the soma into dendrites. In mature neurons, dynamic Golgi elements, containing cis and trans cisternae, are transported from the soma along dendrites, in an actin-dependent process. Dendritic Golgi outposts are dynamic and display bidirectional movement. Similar structures were observed in cerebral organoids. Using the retention using selective hooks (RUSH) system, Golgi resident proteins are transported efficiently into Golgi outposts from the endoplasmic reticulum. This study reveals dynamic, functional Golgi structures in dendrites and a spatial map for investigating dendrite trafficking in human neurons.

A semi-automated pipeline for quantifying drusen-like deposits in human induced pluripotent stem cell-derived retinal pigment epithelium cells.

Age-Related Macular Degeneration (AMD) is a highly prevalent form of retinal disease amongst Western communities over 50 years of age. A hallmark of AMD pathogenesis is the accumulation of drusen underneath the retinal pigment epithelium (RPE), a biological process also observable in vitro. The accumulation of drusen has been shown to predict the progression to advanced AMD, making accurate characterisation of drusen in vitro models valuable in disease modelling and drug development. More recently, deposits above the RPE in the subretinal space, called reticular pseudodrusen (RPD) have been recognized as a sub-phenotype of AMD. While in vitro imaging techniques allow for the immunostaining of drusen-like deposits, quantification of these deposits often requires slow, low throughput manual counting of images. This further lends itself to issues including sampling biases, while ignoring critical data parameters including volume and precise localization. To overcome these issues, we developed a semi-automated pipeline for quantifying the presence of drusen-like deposits in vitro, using RPE cultures derived from patient-specific induced pluripotent stem cells (iPSCs). Using high-throughput confocal microscopy, together with three-dimensional reconstruction, we developed an imaging and analysis pipeline that quantifies the number of drusen-like deposits, and accurately and reproducibly provides the location and composition of these deposits. Extending its utility, this pipeline can determine whether the drusen-like deposits locate to the apical or basal surface of RPE cells. Here, we validate the utility of this pipeline in the quantification of drusen-like deposits in six iPSCs lines derived from patients with AMD, following their differentiation into RPE cells. This pipeline provides a valuable tool for the in vitro modelling of AMD and other retinal disease, and is amenable to mid and high throughput screenings.

AAV capsid bioengineering in primary human retina models.

Adeno-associated viral (AAV) vector-mediated retinal gene therapy is an active field of both pre-clinical as well as clinical research. As with other gene therapy clinical targets, novel bioengineered AAV variants developed by directed evolution or rational design to possess unique desirable properties, are entering retinal gene therapy translational programs. However, it is becoming increasingly evident that predictive preclinical models are required to develop and functionally validate these novel AAVs prior to clinical studies. To investigate if, and to what extent, primary retinal explant culture could be used for AAV capsid development, this study performed a large high-throughput screen of 51 existing AAV capsids in primary human retina explants and other models of the human retina. Furthermore, we applied transgene expression-based directed evolution to develop novel capsids for more efficient transduction of primary human retina cells and compared the top variants to the strongest existing benchmarks identified in the screening described above. A direct side-by-side comparison of the newly developed capsids in four different in vitro and ex vivo model systems of the human retina allowed us to identify novel AAV variants capable of high transgene expression in primary human retina cells.

A village in a dish model system for population-scale hiPSC studies.

The mechanisms by which DNA alleles contribute to disease risk, drug response, and other human phenotypes are highly context-specific, varying across cell types and different conditions. Human induced pluripotent stem cells are uniquely suited to study these context-dependent effects but cell lines from hundreds or thousands of individuals are required. Village cultures, where multiple induced pluripotent stem lines are cultured and differentiated in a single dish, provide an elegant solution for scaling induced pluripotent stem experiments to the necessary sample sizes required for population-scale studies. Here, we show the utility of village models, demonstrating how cells can be assigned to an induced pluripotent stem line using single-cell sequencing and illustrating that the genetic, epigenetic or induced pluripotent stem line-specific effects explain a large percentage of gene expression variation for many genes. We demonstrate that village methods can effectively detect induced pluripotent stem line-specific effects, including sensitive dynamics of cell states.

Culture Variabilities of Human iPSC-Derived Cerebral Organoids Are a Major Issue for the Modelling of Phenotypes Observed in Alzheimer's Disease.

Apolipoprotein E (APOE) is the most important susceptibility gene for late onset of Alzheimer's disease (AD), with the presence of APOE-ε4 associated with increased risk of developing AD. Here, we reprogrammed human fibroblasts from individuals with different APOE-ε genotypes into induced pluripotent stem cells (iPSCs), and generated isogenic lines with different APOE profiles. Following characterisation of the newly established iPSC lines, we used an unguided/unpatterning differentiation method to generate six-month-old cerebral organoids from all iPSC lines to assess the suitability of this in vitro system to measure APOE, ß amyloid, and Tau phosphorylation levels. We identified variabilities in the organoids' cell composition between cell lines, and between batches of differentiation for each cell line. We observed more homogenous cerebral organoids, and similar levels of APOE, ß amyloid, and Tau when using the CRISPR-edited APOE isogenic lines, with the exception of one site of Tau phosphorylation which was higher in the APOE-ε4/ε4 organoids. These data describe that pathological hallmarks of AD are observed in cerebral organoids, and that their variation is mainly independent of the APOE-ε status of the cells, but associated with the high variability of cerebral organoid differentiation. It demonstrates that the cell-line-to-cell-line and batch-to-batch variabilities need to be considered when using cerebral organoids.

AAV2-mediated gene therapy for Bietti crystalline dystrophy provides functional CYP4V2 in multiple relevant cell models.

Bietti crystalline dystrophy (BCD) is an inherited retinal disease (IRD) caused by mutations in the CYP4V2 gene. It is a relatively common cause of IRD in east Asia. A number of features of this disease make it highly amenable to gene supplementation therapy. This study aims to validate a series of essential precursor in vitro experiments prior to developing a clinical gene therapy for BCD. We demonstrated that HEK293, ARPE19, and patient induced pluripotent stem cell (iPSC)-derived RPE cells transduced with AAV2 vectors encoding codon optimization of CYP4V2 (AAV2.coCYP4V2) resulted in elevated protein expression levels of CYP4V2 compared to those transduced with AAV2 vectors encoding wild type CYP4V2 (AAV2.wtCYP4V2), as assessed by immunocytochemistry and western blot. Similarly, we observed significantly increased CYP4V2 enzyme activity in cells transduced with AAV2.coCYP4V2 compared to those transduced with AAV2.wtCYP4V2. We also showed CYP4V2 expression in human RPE/choroid explants transduced with AAV2.coCYP4V2 compared to those transduced with AAV2.wtCYP4V2. These preclinical data support the further development of a gene supplementation therapy for a currently untreatable blinding condition-BCD. Codon-optimized CYP4V2 transgene was superior to wild type in terms of protein expression and enzyme activity. Ex vivo culture of human RPE cells provided an effective approach to test AAV-mediated transgene delivery.

Transcriptomic and proteomic retinal pigment epithelium signatures of age-related macular degeneration.

There are currently no treatments for geographic atrophy, the advanced form of age-related macular degeneration. Hence, innovative studies are needed to model this condition and prevent or delay its progression. Induced pluripotent stem cells generated from patients with geographic atrophy and healthy individuals were differentiated to retinal pigment epithelium. Integrating transcriptional profiles of 127,659 retinal pigment epithelium cells generated from 43 individuals with geographic atrophy and 36 controls with genotype data, we identify 445 expression quantitative trait loci in cis that are asssociated with disease status and specific to retinal pigment epithelium subpopulations. Transcriptomics and proteomics approaches identify molecular pathways significantly upregulated in geographic atrophy, including in mitochondrial functions, metabolic pathways and extracellular cellular matrix reorganization. Five significant protein quantitative trait loci that regulate protein expression in the retinal pigment epithelium and in geographic atrophy are identified - two of which share variants with cis- expression quantitative trait loci, including proteins involved in mitochondrial biology and neurodegeneration. Investigation of mitochondrial metabolism confirms mitochondrial dysfunction as a core constitutive difference of the retinal pigment epithelium from patients with geographic atrophy. This study uncovers important differences in retinal pigment epithelium homeostasis associated with geographic atrophy.

Retinal ganglion cell-specific genetic regulation in primary open-angle glaucoma.

To assess the transcriptomic profile of disease-specific cell populations, fibroblasts from patients with primary open-angle glaucoma (POAG) were reprogrammed into induced pluripotent stem cells (iPSCs) before being differentiated into retinal organoids and compared with those from healthy individuals. We performed single-cell RNA sequencing of a total of 247,520 cells and identified cluster-specific molecular signatures. Comparing the gene expression profile between cases and controls, we identified novel genetic associations for this blinding disease. Expression quantitative trait mapping identified a total of 4,443 significant loci across all cell types, 312 of which are specific to the retinal ganglion cell subpopulations, which ultimately degenerate in POAG. Transcriptome-wide association analysis identified genes at loci previously associated with POAG, and analysis, conditional on disease status, implicated 97 statistically significant retinal ganglion cell-specific expression quantitative trait loci. This work highlights the power of large-scale iPSC studies to uncover context-specific profiles for a genetically complex disease.

Automation of Organoid Cultures: Current Protocols and Applications.

Organoids are three-dimensional, functional structures that mimic in vivo organs. They offer new opportunities for the modeling of cancer and infectious and rare hereditary diseases. Furthermore, the advent of organoid biobanks opens new avenues for drug screening in a personalized fashion and holds much promise for personalized regenerative medicine. Thus, there is a need for reproducible, large-scale organoid generation with minimal variability, making manual approaches impracticable. Here, we review the current use of automation in organoid culture and analysis, using cerebral and retinal organoids as illustrations of current applications. An increased demand for automated organoid platforms is anticipated.

Generation of a gene-corrected human isogenic iPSC line from an Alzheimer's disease iPSC line carrying the London mutation in APP (V717I).

We report the genome-editing of an existing iPSC line carrying the London mutation in APP (V717I) into an iPSC line in which the pathogenic mutation was corrected. The resulting isogenic iPSC line maintained pluripotent stem cell morphology, a normal karyotype, expression of pluripotency markers and the ability to differentiate into the three germ-layers in vitro.

Single cell eQTL analysis identifies cell type-specific genetic control of gene expression in fibroblasts and reprogrammed induced pluripotent stem cells.

BACKGROUND: The discovery that somatic cells can be reprogrammed to induced pluripotent stem cells (iPSCs) has provided a foundation for in vitro human disease modelling, drug development and population genetics studies. Gene expression plays a critical role in complex disease risk and therapeutic response. However, while the genetic background of reprogrammed cell lines has been shown to strongly influence gene expression, the effect has not been evaluated at the level of individual cells which would provide significant resolution. By integrating single cell RNA-sequencing (scRNA-seq) and population genetics, we apply a framework in which to evaluate cell type-specific effects of genetic variation on gene expression. RESULTS: Here, we perform scRNA-seq on 64,018 fibroblasts from 79 donors and map expression quantitative trait loci (eQTLs) at the level of individual cell types. We demonstrate that the majority of eQTLs detected in fibroblasts are specific to an individual cell subtype. To address if the allelic effects on gene expression are maintained following cell reprogramming, we generate scRNA-seq data in 19,967 iPSCs from 31 reprogramed donor lines. We again identify highly cell type-specific eQTLs in iPSCs and show that the eQTLs in fibroblasts almost entirely disappear during reprogramming. CONCLUSIONS: This work provides an atlas of how genetic variation influences gene expression across cell subtypes and provides evidence for patterns of genetic architecture that lead to cell type-specific eQTL effects.

Comparative performance of the BGI and Illumina sequencing technology for single-cell RNA-sequencing.

The libraries generated by high-throughput single cell RNA-sequencing (scRNA-seq) platforms such as the Chromium from 10× Genomics require considerable amounts of sequencing, typically due to the large number of cells. The ability to use these data to address biological questions is directly impacted by the quality of the sequence data. Here we have compared the performance of the Illumina NextSeq 500 and NovaSeq 6000 against the BGI MGISEQ-2000 platform using identical Single Cell 3' libraries consisting of over 70 000 cells generated on the 10× Genomics Chromium platform. Our results demonstrate a highly comparable performance between the NovaSeq 6000 and MGISEQ-2000 in sequencing quality, and the detection of genes, cell barcodes, Unique Molecular Identifiers. The performance of the NextSeq 500 was also similarly comparable to the MGISEQ-2000 based on the same metrics. Data generated by both sequencing platforms yielded similar analytical outcomes for general single-cell analysis. The performance of the NextSeq 500 and MGISEQ-2000 were also comparable for the deconvolution of multiplexed cell pools via variant calling, and detection of guide RNA (gRNA) from a pooled CRISPR single-cell screen. Our study provides a benchmark for high-capacity sequencing platforms applied to high-throughput scRNA-seq libraries.

Automated Cell Culture Systems and Their Applications to Human Pluripotent Stem Cell Studies.

Pluripotent stem cells are an extremely powerful tool in modeling human diseases and hold much promise for personalized regenerative or cell replacement therapies. There is an increasing need for reproducible large-scale stem cell and differentiated progeny production, with minimal variation, rendering manual approaches impracticable. Here, we provide an overview of systems currently available for automated stem cell culture, and undertake a review of their capacities, capabilities, and relative limitations. With the merging of modern technology and stem cell biology, an increased demand and implementation of automated platforms for stem cell studies is anticipated.

Single cell RNA sequencing of stem cell-derived retinal ganglion cells.

We used single cell sequencing technology to characterize the transcriptomes of 1,174 human embryonic stem cell-derived retinal ganglion cells (RGCs) at the single cell level. The human embryonic stem cell line BRN3B-mCherry (A81-H7), was differentiated to RGCs using a guided differentiation approach. Cells were harvested at day 36 and prepared for single cell RNA sequencing. Our data indicates the presence of three distinct subpopulations of cells, with various degrees of maturity. One cluster of 288 cells showed increased expression of genes involved in axon guidance together with semaphorin interactions, cell-extracellular matrix interactions and ECM proteoglycans, suggestive of a more mature RGC phenotype.

Single-Cell Profiling Identifies Key Pathways Expressed by iPSCs Cultured in Different Commercial Media.

We assessed the pluripotency of human induced pluripotent stem cells (iPSCs) maintained on an automated platform using StemFlex and TeSR-E8 media. Analysis of transcriptome of single cells revealed similar expression of core pluripotency genes, as well as genes associated with naive and primed states of pluripotency. Analysis of individual cells from four samples consisting of two different iPSC lines each grown in the two culture media revealed a shared subpopulation structure with three main subpopulations different in pluripotency states. By implementing a machine learning approach, we estimated that most cells within each subpopulation are very similar between all four samples. The single-cell RNA sequencing analysis of iPSC lines grown in both media reports the molecular signature in StemFlex medium and how it compares to that observed in the TeSR-E8 medium.

Generation of a human induced pluripotent stem cell line CERAi001-A-6 using episomal vectors.

We report the generation of the hiPSC line CERAi001-A-6 from primary human dermal fibroblasts. Reprogramming was performed using episomal vector delivery of OCT4, SOX2, KLF4, L-MYC, LIN28 and shRNA for p53.