Repeated evolution of eye loss in Mexican cavefish: Evidence of similar developmental mechanisms in independently evolved populations.

Evolution in similar environments often leads to convergence of behavioral and anatomical traits. A classic example of convergent trait evolution is the reduced traits that characterize many cave animals: reduction or loss of pigmentation and eyes. While these traits have evolved many times, relatively little is known about whether these traits repeatedly evolve through the same or different molecular and developmental mechanisms. The small freshwater fish, Astyanax mexicanus, provides an opportunity to investigate the repeated evolution of cave traits. A. mexicanus exists as two forms, a sighted, surface-dwelling form and at least 29 populations of a blind, cave-dwelling form that initially develops eyes that subsequently degenerate. We compared eye morphology and the expression of eye regulatory genes in developing surface fish and two independently evolved cavefish populations, Pachón and Molino. We found that many of the previously described molecular and morphological alterations that occur during eye development in Pachón cavefish are also found in Molino cavefish. However, for many of these traits, the Molino cavefish have a less severe phenotype than Pachón cavefish. Further, cave-cave hybrid fish have larger eyes and lenses during early development compared with fish from either parental population, suggesting that some different changes underlie eye loss in these two populations. Together, these data support the hypothesis that these two cavefish populations evolved eye loss independently, yet through some of the same developmental and molecular mechanisms.

Transplantation of iPSC-derived TM cells rescues glaucoma phenotypes in vivo.

Glaucoma is a common cause of vision loss or blindness and reduction of intraocular pressure (IOP) has been proven beneficial in a large fraction of glaucoma patients. The IOP is maintained by the trabecular meshwork (TM) and the elevation of IOP in open-angle glaucoma is associated with dysfunction and loss of the postmitotic cells residing within this tissue. To determine if IOP control can be maintained by replacing lost TM cells, we transplanted TM-like cells derived from induced pluripotent stem cells into the anterior chamber of a transgenic mouse model of glaucoma. Transplantation led to significantly reduced IOP and improved aqueous humor outflow facility, which was sustained for at least 9 wk. The ability to maintain normal IOP engendered survival of retinal ganglion cells, whose loss is ultimately the cause for reduced vision in glaucoma. In vivo and in vitro analyses demonstrated higher TM cellularity in treated mice compared with littermate controls and indicated that this increase is primarily because of a proliferative response of endogenous TM cells. Thus, our study provides in vivo demonstration that regeneration of the glaucomatous TM is possible and points toward novel approaches in the treatment of this disease.

Complementary feature selection from alternative splicing events and gene expression for phenotype prediction.

MOTIVATION: A central task of bioinformatics is to develop sensitive and specific means of providing medical prognoses from biomarker patterns. Common methods to predict phenotypes in RNA-Seq datasets utilize machine learning algorithms trained via gene expression. Isoforms, however, generated from alternative splicing, may provide a novel and complementary set of transcripts for phenotype prediction. In contrast to gene expression, the number of isoforms increases significantly due to numerous alternative splicing patterns, resulting in a prioritization problem for many machine learning algorithms. This study identifies the empirically optimal methods of transcript quantification, feature engineering and filtering steps using phenotype prediction accuracy as a metric. At the same time, the complementary nature of gene and isoform data is analyzed and the feasibility of identifying isoforms as biomarker candidates is examined. RESULTS: Isoform features are complementary to gene features, providing non-redundant information and enhanced predictive power when prioritized and filtered. A univariate filtering algorithm, which selects up to the N highest ranking features for phenotype prediction is described and evaluated in this study. An empirical comparison of pipelines for isoform quantification is reported by performing cross-validation prediction tests with datasets from human non-small cell lung cancer (NSCLC) patients, human patients with chronic obstructive pulmonary disease (COPD) and amyotrophic lateral sclerosis (ALS) transgenic mice, each including samples of diseased and non-diseased phenotypes. AVAILABILITY AND IMPLEMENTATION: https://github.com/clabuzze/Phenotype-Prediction-Pipeline.git CONTACT: clabuzze@iastate.edu, antoniom@bc.edu, watsondk@musc.edu, andersonpe2@cofc.edu.

Transcription factor Olig2 defines subpopulations of retinal progenitor cells biased toward specific cell fates.

Previous lineage analyses have shown that retinal progenitor cells (RPCs) are multipotent throughout development, and expression-profiling studies have shown a great deal of molecular heterogeneity among RPCs. To determine if the molecular heterogeneity predicts that an RPC will produce particular types of progeny, clonal lineage analysis was used to investigate the progeny of a subset of RPCs, those that express the basic helix-loop-helix transcription factor, Olig2. The embryonic Olig2(+) RPCs underwent terminal divisions, producing small clones with primarily two of the five cell types being made by the pool of RPCs at that time. The later, postnatal Olig2(+) RPCs also made terminal divisions, which were biased toward production of rod photoreceptors and amacrine cell interneurons. These data indicate that the multipotent progenitor pool is made up of distinctive types of RPCs, which have biases toward producing subsets of retinal neurons in a terminal division, with the types of neurons produced varying over time. This strategy is similar to that of the developing Drosophila melanogaster ventral nerve cord, with the Olig2(+) cells behaving as ganglion mother cells.

Analysis of gene expression in wild-type and Notch1 mutant retinal cells by single cell profiling.

BACKGROUND: The vertebrate retina comprises sensory neurons, the photoreceptors, as well as many other types of neurons and one type of glial cell. These cells are generated by multipotent and restricted retinal progenitor cells (RPCs), which express Notch1. Loss of Notch1 in RPCs late during retinal development results in the overproduction of rod photoreceptors at the expense of interneurons and glia. RESULTS: To examine the molecular underpinnings of this observation, microarray analysis of single retinal cells from wild-type or Notch1 conditional knockout retinas was performed. In situ hybridization was carried out to validate some of the findings. CONCLUSIONS: The majority of Notch1-mutant cells lost expression of known Notch target genes. These cells also had low levels of RPC and cell cycle genes, and robustly up-regulated rod precursor genes. In addition, single wild-type cells, in which cell cycle marker genes were down-regulated, expressed markers of both rod photoreceptors and interneurons.

Development and diversification of retinal amacrine interneurons at single cell resolution.

The vertebrate retina uses diverse neuronal cell types arrayed into complex neural circuits to extract, process, and relay information from the visual scene to the higher order processing centers of the brain. Amacrine cells, a class of interneurons, are thought to mediate much of the processing of the visual signal that occurs within the retina. Although amacrine cells display extensive morphological diversity, the molecular nature of this diversity is largely unknown. Furthermore, it is not known how this diversity arises during development. Here, we have combined in vivo genetic labeling, single cell genome-wide expression profiling, and classical birthdating to (i) identify specific molecular types of amacrine cells, (ii) demonstrate the molecular diversity of the amacrine cell class, and (iii) show that amacrine cell diversity arises at least in part through temporal patterning.

Identification of genes expressed preferentially in the developing peripheral margin of the optic cup.

Specification of the peripheral optic cup by Wnt signaling is critical for formation of the ciliary body/iris. Identification of marker genes for this region during development provides a starting point for functional analyses. During transcriptional profiling of single cells from the developing eye, two cells were identified that expressed genes not found in most other single cell profiles. In situ hybridizations demonstrated that many of these genes were expressed in the peripheral optic cup in both early mouse and chicken development, and in the ciliary body/iris at subsequent developmental stages. These analyses indicate that the two cells probably originated from the developing ciliary body/iris. Changes in expression of these genes were assayed in embryonic chicken retinas when canonical Wnt signaling was ectopically activated by CA-beta-catenin. Twelve ciliary body/iris genes were identified as upregulated following induction, suggesting they are excellent candidates for downstream effectors of Wnt signaling in the optic cup.

AAV vectors engineered to target insulin receptor greatly enhance intramuscular gene delivery.

Adeno-associated virus (AAV) is one of the most commonly used vectors for gene therapy, and the applications for AAV-delivered therapies are numerous. However, the current state of technology is limited by the low efficiency with which most AAV vectors transduce skeletal muscle tissue. We demonstrate that vector efficiency can be enhanced by modifying the AAV capsid with a peptide that binds a receptor highly expressed in muscle tissue. When an insulin-mimetic peptide, S519, previously characterized for its high affinity to insulin receptor (IR), was inserted into the capsid, the AAV9 transduction efficiency of IR-expressing cell lines as well as differentiated primary human muscle cells was dramatically enhanced. This vector also exhibited efficient transduction of mouse muscle in vivo, resulting in up to 18-fold enhancement over AAV9. Owing to its superior transduction efficiency in skeletal muscle, we named this vector "enhanced AAV9" (eAAV9). We also found that the modification enhanced the transduction efficiency of several other AAV serotypes. Together, these data show that AAV transduction of skeletal muscle can be improved by targeting IR. They also show the broad utility of this modular strategy and suggest that it could also be applied to next-generation vectors that have yet to be engineered.

Molecular signatures of retinal ganglion cells revealed through single cell profiling.

Retinal ganglion cells can be classified into more than 40 distinct subtypes, whether by functional classification or transcriptomics. The examination of these subtypes in relation to their physiology, projection patterns, and circuitry would be greatly facilitated through the identification of specific molecular identifiers for the generation of transgenic mice. Advances in single cell transcriptomic profiling have enabled the identification of molecular signatures for cellular subtypes that are only rarely found. Therefore, we used single cell profiling combined with hierarchical clustering and correlate analyses to identify genes expressed in distinct populations of Parvalbumin-expressing cells and functionally classified RGCs. RGCs were manually isolated based either upon fluorescence or physiological distinction through cell-attached recordings. Microarray hybridization and RNA-Sequencing were employed for the characterization of transcriptomes and in situ hybridization was utilized to further characterize gene candidate expression. Gene candidates were identified based upon cluster correlation, as well as expression specificity within physiologically distinct classes of RGCs. Further, we identified Prph, Ctxn3, and Prkcq as potential candidates for ipRGC classification in the murine retina. The use of these genes, or one of the other newly identified subset markers, for the generation of a transgenic mouse would enable future studies of RGC-subtype specific function, wiring, and projection.

Thyroid hormone components are expressed in three sequential waves during development of the chick retina.

BACKGROUND: Thyroid hormone (TH) is an important developmental regulator in many tissues, including the retina. TH is activated locally via deiodinase 2 (Dio2), and it is destroyed by deiodinase 3 (Dio3). The TH receptors, TRa and TRb, mediate TH activity through hormone and DNA binding, and interactions with transcription regulators. RESULTS: In the current work, the expression of these TH components was examined in the chick retina over time. Three waves of expression were characterized and found to be correlated with critical developmental events. The first wave occurred as progenitor cells began to make photoreceptors, the second as some cell types adopted a more mature location and differentiation state, and the third as Müller glia were generated. The cell types expressing the components, as well as the kinetics of expression within the cell cycle, were defined. TRb expression initiated during G2 in progenitor cells, concomitant with NeuroD and Otx2, which are expressed in early photoreceptor cells. TRb was expressed in photoreceptor cells for several days and then was reduced in expression level, as the expression of Crx, a later photoreceptor gene, became more evident. Dio3 was expressed throughout the cell cycle in progenitor cells. TRa was in most, if not all, retinal cells. Dio2 appeared transiently in a ventral (high) to dorsal gradient, likely in a subset of photoreceptor cells. CONCLUSION: Multiple TH components were expressed in dynamic patterns in cycling progenitor cells and photoreceptors cells across the developing chick retina. These dynamic patterns suggest that TH is playing several roles in retinal development, both within the cycling progenitor cells and possibly with respect to the timing of differentiation of photoreceptor cells.

The Trim family of genes and the retina: Expression and functional characterization.

To better understand the mechanisms that govern the development of retinal neurons, it is critical to gain additional insight into the specific intrinsic factors that control cell fate decisions and neuronal maturation. In the developing mouse retina, Atoh7, a highly conserved transcription factor, is essential for retinal ganglion cell development. Moreover, Atoh7 expression in the developing retina occurs during a critical time period when progenitor cells are in the process of making cell fate decisions. We performed transcriptome profiling of Atoh7+ individual cells isolated from mouse retina. One of the genes that we found significantly correlated with Atoh7 in our transcriptomic data was the E3 ubiquitin ligase, Trim9. The correlation between Trim9 and Atoh7 coupled with the expression of Trim9 in the early mouse retina led us to hypothesize that this gene may play a role in the process of cell fate determination. To address the role of Trim9 in retinal development, we performed a functional analysis of Trim9 in the mouse and did not detect any morphological changes in the retina in the absence of Trim9. Thus, Trim9 alone does not appear to be involved in cell fate determination or early ganglion cell development in the mouse retina. We further hypothesize that the reason for this lack of phenotype may be compensation by one of the many additional TRIM family members we find expressed in the developing retina.

ALS-associated genes display CNS expression in the developing zebrafish.

Amyotrophic lateral sclerosis (ALS) is characterized by progressive muscle atrophy resulting from the deterioration of motor neurons in the central nervous system (CNS). Recent genome-wide association studies have revealed several genes linked to ALS, further demonstrating the complexity of the disease. The zebrafish (Danio rerio) is an attractive model organism to study the function of the rapidly expanding number of ALS-associated genes, in part, due to the development of genome editing techniques that have facilitated specific gene targeting. Before investing in the manipulation and phenotypic examination of these genes, however, it is important to ascertain the localization of expression in this organism. We performed an expression analysis of 29 total ALS-linked genes in the developing zebrafish, specifically focusing on those genes that displayed robust and reproducible expression at multiple different timepoints. First, we classified a subset of the most robustly expressed genes into three distinct groups: head-only expression, head and weak trunk expression, and head and robust trunk expression. Then, we defined the characteristic pattern of each gene at 2, 3, and 4 days post fertilization. This analysis will facilitate improved mutant phenotype assessment in zebrafish by focusing researchers on the areas of expression.

Molecular heterogeneity of developing retinal ganglion and amacrine cells revealed through single cell gene expression profiling.

During development of the central nervous system (CNS), cycling uncommitted progenitor cells give rise to a variety of distinct neuronal and glial cell types. As these different cell types are born they progress from newly specified cells to fully differentiated neurons and glia. In order to define the developmental processes of individual cell types, single cell expression profiling was carried out on developing ganglion and amacrine cells of the murine retina. Individual cells from multiple developmental stages were isolated and profiled on Affymetrix oligonucleotide arrays. Two-color fluorescent in situ hybridization on dissociated retinas was used to verify and extend the microarray results by allowing quantitative measurements of a large number of cells coexpressing two genes. Together, these experiments have yielded an expanded view of the processes underway in developing retinal ganglion and amacrine cells, as well as several hundred new marker genes for these cell types. In addition, this study has allowed for the definition of some of the molecular heterogeneity both between developing ganglion and amacrine cells and among subclasses of each cell type.

Single-cell RNA-Seq of Defined Subsets of Retinal Ganglion Cells.

The discovery of cell type-specific markers can provide insight into cellular function and the origins of cellular heterogeneity. With a recent push for the improved understanding of neuronal diversity, it is important to identify genes whose expression defines various subpopulations of cells. The retina serves as an excellent model for the study of central nervous system diversity, as it is composed of multiple major cell types. The study of each major class of cells has yielded genetic markers that facilitate the identification of these populations. However, multiple subtypes of cells exist within each of these major retinal cell classes, and few of these subtypes have known genetic markers, although many have been characterized by morphology or function. A knowledge of genetic markers for individual retinal subtypes would allow for the study and mapping of brain targets related to specific visual functions and may also lend insight into the gene networks that maintain cellular diversity. Current avenues used to identify the genetic markers of subtypes possess drawbacks, such as the classification of cell types following sequencing. This presents a challenge for data analysis and requires rigorous validation methods to ensure that clusters contain cells of the same function. We propose a technique for identifying the morphology and functionality of a cell prior to isolation and sequencing, which will allow for the easier identification of subtype-specific markers. This technique may be extended to non-neuronal cell types, as well as to rare populations of cells with minor variations. This protocol yields excellent-quality data, as many of the libraries have provided read depths greater than 20 million reads for single cells. This methodology overcomes many of the hurdles presented by Single-cell RNA-Seq and may be suitable for researchers aiming to profile cell types in a straightforward and highly efficient manner.

Single cell transcriptome profiling of developing chick retinal cells.

The vertebrate retina is a specialized photosensitive tissue comprised of six neuronal and one glial cell types, each of which develops in prescribed proportions at overlapping timepoints from a common progenitor pool. While each of these cells has a specific function contributing to proper vision in the mature animal, their differential representation in the retina as well as the presence of distinctive cellular subtypes makes identifying the transcriptomic signatures that lead to each retinal cell's fate determination and development challenging. We have analyzed transcriptomes from individual cells isolated from the chick retina throughout retinogenesis. While we focused our efforts on the retinal ganglion cells, our transcriptomes of developing chick cells also contained representation from multiple retinal cell types, including photoreceptors and interneurons at different stages of development. Most interesting was the identification of transcriptomes from individual mixed lineage progenitor cells in the chick as these cells offer a window into the cell fate decision-making process. Taken together, these data sets will enable us to uncover the most critical genes acting in the steps of cell fate determination and early differentiation of various retinal cell types.

Polo-Like Kinase 3 Appears Dispensable for Normal Retinal Development Despite Robust Embryonic Expression.

During retinogenesis seven different cell types are generated in distinct yet overlapping timepoints from a population of retinal progenitor cells. Previously, we performed single cell transcriptome analyses of retinal progenitor cells to identify candidate genes that may play roles in the generation of early-born retinal neurons. Based on its expression pattern in subsets of early retinal cells, polo-like kinase 3 (Plk3) was identified as one such candidate gene. Further characterization of Plk3 expression by in situ hybridization revealed that this gene is expressed as cells exit the cell cycle. We obtained a Plk3 deficient mouse and investigated changes in the retina's morphology and transcriptome through immunohistochemistry, in situ hybridization and gene expression profiling. These experiments have been performed initially on adult mice and subsequently extended throughout retinal development. Although morphological studies revealed no consistent changes in retinogenesis upon Plk3 loss, microarray profiling revealed potential candidate genes altered in Plk3-KO mice. Further studies will be necessary to understand the connection between these changes in gene expression and the loss of a protein kinase such as Plk3.

Disentangling neural cell diversity using single-cell transcriptomics.

Cellular specialization is particularly prominent in mammalian nervous systems, which are composed of millions to billions of neurons that appear in thousands of different 'flavors' and contribute to a variety of functions. Even in a single brain region, individual neurons differ greatly in their morphology, connectivity and electrophysiological properties. Systematic classification of all mammalian neurons is a key goal towards deconstructing the nervous system into its basic components. With the recent advances in single-cell gene expression profiling technologies, it is now possible to undertake the enormous task of disentangling neuronal heterogeneity. High-throughput single-cell RNA sequencing and multiplexed quantitative RT-PCR have become more accessible, and these technologies enable systematic categorization of individual neurons into groups with similar molecular properties. Here we provide a conceptual and practical guide to classification of neural cell types using single-cell gene expression profiling technologies.

Expression Profiling of Developing Zebrafish Retinal Cells.

During retinal development, a variety of different types of neurons are produced. Understanding how each of these types of retinal nerve cells is generated is important from a developmental biology perspective. It is equally important if one is interested in how to regenerate cells after an injury or a disease. To gain more insight into how retinal neurons develop in the zebrafish, we performed single-cell mRNA profiling and in situ hybridizations (ISHs) on retinal sections and whole-mount zebrafish. Through the series of ISHs, designed and performed solely by undergraduate students in the laboratory, we were able to retrospectively identify our single-cell mRNA profiles as most likely coming from developing amacrine cells. Further analysis of these profiles will reveal genes that can be mutated using genome editing techniques. Together these studies increase our knowledge of the genes driving development of different cell types in the zebrafish retina.

Transcriptomic analyses of Onecut1 and Onecut2 deficient retinas.

In this article, we further explore the data generated for the research article "Onecut1 and Onecut2 play critical roles in the development of the mouse retina". To better understand the functionality of the Onecut family of transcription factors in retinogenesis, we investigated the retinal transcriptomes of developing and mature mice to identify genes with differential expression. This data article reports the full transcriptomes resulting from these experiments and provides tables detailing the differentially expressed genes between wildtype and Onecut1 or 2 deficient retinas. The raw array data of our transcriptomes as generated using Affymetrix microarrays are available on the NCBI Gene Expression Omnibus (GEO) browser (Reference number GSE57917 and GSE57918GSE57917GSE57918).

Onecut1 and Onecut2 play critical roles in the development of the mouse retina.

The entire repertoire of intrinsic factors that control the cell fate determination process of specific retinal neurons has yet to be fully identified. Single cell transcriptome profiling experiments of retinal progenitor cells revealed considerable gene expression heterogeneity between individual cells, especially among different classes of transcription factors. In this study, we show that two of those factors, Onecut1 and Onecut2, are expressed during mouse retinal development. Using mice that are deficient for each of these transcription factors, we further demonstrate a significant loss (∼70-80%) of horizontal cells in the absence of either of these proteins, while the other retinal cells appear at normal numbers. Microarray profiling experiments performed on knockout retinas revealed defects in horizontal cell genes as early as E14.5. Additional profiling assays showed an upregulation of several stress response genes in the adult Onecut2 knockout, suggesting that the integrity of the retina is compromised in the absence of normal numbers of horizontal cells. Interestingly, melanopsin, the gene coding for the photopigment found in photosensitive ganglion cells, was observed to be upregulated in Onecut1 deficient retinas, pointing to a possible regulatory role for Onecut1. Taken together, our data show that similar to Onecut1, Onecut2 is also necessary for the formation of normal numbers of horizontal cells in the developing retina.