Gene expression variation in the adult human retina

Translatability barriers between preclinical and clinical trials of AAV gene therapy in inherited retinal diseases

Inherited retinal diseases (IRDs) are progressive degenerative diseases which cause gradual vision loss or complete blindness. As over 270 gene mutations have been identified in the underlying pathology of IRDs, gene therapy as a treatment modality has been an increasingly active realm of investigation. Currently, the most common vehicle of ocular gene delivery is the adeno-associated virus (AAV) vector. This is injected into the immune-privileged subretinal space to mediate transgene expression in retinal cells. Although numerous animal models of IRDs have demonstrated successful outcomes following AAV-mediated gene delivery, many of these studies fail to translate into successful outcomes in clinical trials. The purpose of this review is to A) comparatively assess preclinical and clinical IRD trials in which the success of AAV-mediated therapy failed to translate between animal and human participants B) discuss factors which may complicate the translatability of gene therapy in animals to results in humans.

Mural Serum Response Factor (SRF) Deficiency Provides Insights into Retinal Vascular Functionality and Development

Serum response factor (SRF) controls the expression of muscle contraction and motility genes in mural cells (MCs) of the vasculature. In the retina, MC-SRF is important for correct angiogenesis during development and the continuing maintenance of the vascular tone. The purpose of this study was to provide further insights into the effects of MC SRF deficiency on the vasculature and function of the mature retina in SrfiMCKO mice that carry a MC-specific deletion of Srf. Retinal morphology and vascular integrity were analyzed in vivo via scanning laser ophthalmoscopy (SLO), angiography, and optical coherence tomography (OCT). Retinal function was evaluated with full-field electroretinography (ERG). We found that retinal blood vessels of these mutants exhibited different degrees of morphological and functional alterations. With increasing severity, we found vascular bulging, the formation of arteriovenous (AV) anastomoses, and ultimately, a retinal detachment (RD). The associated irregular retinal blood pressure and flow distribution eventually induced hypoxia, indicated by a negative ERG waveform shape. Further, the high frequency of interocular differences in the phenotype of individual SrfiMCKO mice points to a secondary nature of these developments far downstream of the genetic defect and rather dependent on the local retinal context.

Single-cell RNA sequencing in vision research: Insights into human retinal health and disease

Gene expression provides valuable insight into cell function. As such, vision researchers have frequently employed gene expression studies to better understand retinal physiology and disease. With the advent of single-cell RNA sequencing, expression experiments provide an unparalleled resolution of information. Instead of studying aggregated gene expression across all cells in a heterogenous tissue, single-cell technology maps RNA to an individual cell, which facilitates grouping of retinal and choroidal cell types for further study. Single-cell RNA sequencing has been quickly adopted by both basic and translational vision researchers, and single-cell level gene expression has been studied in the visual systems of animal models, retinal organoids, and primary human retina, RPE, and choroid. These experiments have generated detailed atlases of gene expression and identified new retinal cell types. Likewise, single-cell RNA sequencing investigations have characterized how gene expression changes in the setting of many retinal diseases, including how choroidal endothelial cells are altered in age-related macular degeneration. In addition, this technology has allowed vision researchers to discover drivers of retinal development and model rare retinal diseases with induced pluripotent stem cells. In this review, we will overview the growing number of single-cell RNA sequencing studies in the field of vision research. We will summarize experimental considerations for designing single-cell RNA sequencing experiments and highlight important advancements in retinal, RPE, choroidal, and retinal organoid biology driven by this technology. Finally, we generalize these findings to genes involved in retinal degeneration and outline the future of single-cell expression experiments in studying retinal disease.

Inherited Eye Diseases with Retinal Manifestations through the Eyes of Homeobox Genes

Retinal development is under the coordinated control of overlapping networks of signaling pathways and transcription factors. The paper was conceived as a review of the data and ideas that have been formed to date on homeobox genes mutations that lead to the disruption of eye organogenesis and result in inherited eye/retinal diseases. Many of these diseases are part of the same clinical spectrum and have high genetic heterogeneity with already identified associated genes. We summarize the known key regulators of eye development, with a focus on the homeobox genes associated with monogenic eye diseases showing retinal manifestations. Recent advances in the field of genetics and high-throughput next-generation sequencing technologies, including single-cell transcriptome analysis have allowed for deepening of knowledge of the genetic basis of inherited retinal diseases (IRDs), as well as improve their diagnostics. We highlight some promising avenues of research involving molecular-genetic and cell-technology approaches that can be effective for IRDs therapy. The most promising neuroprotective strategies are aimed at mobilizing the endogenous cellular reserve of the retina.

The Common ABCA4 Variant p.Asn1868Ile Shows Nonpenetrance and Variable Expression of Stargardt Disease When Present in trans With Severe Variants

Purpose

To assess the occurrence and the disease expression of the common p.Asn1868Ile variant in patients with Stargardt disease (STGD1) harboring known, monoallelic causal ABCA4 variants.

Methods

The coding and noncoding regions of ABCA4 were sequenced in 67 and 63 STGD1 probands respectively, harboring monoallelic ABCA4 variants. In case p.Asn1868Ile was detected, segregation analysis was performed whenever possible. Probands and affected siblings harboring p.Asn1868Ile without additional variants in cis were clinically evaluated retrospectively. Two asymptomatic siblings carrying the same ABCA4 variants as their probands were clinically examined. The penetrance of p.Asn1868Ile was calculated using allele frequency data of ABCA4 variants in non-Finnish European individuals.

Results

The p.Asn1868Ile variant was found in cis with known variants in 14/67 probands. In 27/67 probands, we identified p.Asn1868Ile without additional variants in cis, in combination with known, mainly severe ABCA4 variants. In 23/27 probands, the trans configuration was established. Among 27 probands and 6/7 STGD1 siblings carrying p.Asn1868Ile, 42% manifested late-onset disease (>44 years). We additionally identified four asymptomatic relatives carrying a combination of a severe variant and p.Asn1868Ile; ophthalmologic examination in two persons did not reveal STGD1. Based on ABCA4 allele frequency data, we conservatively estimated the penetrance of p.Asn1868Ile, when present in trans with a severe variant, to be below 5%.

Conclusions

A significant fraction of genetically unexplained STGD1 cases carries p.Asn1868Ile as a second variant. Our findings suggest exceptional differences in disease expression or even nonpenetrance of this ABCA4 variant, pointing toward an important role for genetic or environmental modifiers in STGD1.

Transcriptome profiling of aging Drosophila photoreceptors reveals gene expression trends that correlate with visual senescence

BACKGROUND

Aging is associated with functional decline of neurons and increased incidence of both neurodegenerative and ocular disease. Photoreceptor neurons in Drosophila melanogaster provide a powerful model for studying the molecular changes involved in functional senescence of neurons since decreased visual behavior precedes retinal degeneration. Here, we sought to identify gene expression changes and the genomic features of differentially regulated genes in photoreceptors that contribute to visual senescence.

RESULTS

To identify gene expression changes that could lead to visual senescence, we characterized the aging transcriptome of Drosophila sensory neurons highly enriched for photoreceptors. We profiled the nuclear transcriptome of genetically-labeled photoreceptors over a 40 day time course and identified increased expression of genes involved in stress and DNA damage response, and decreased expression of genes required for neuronal function. We further show that combinations of promoter motifs robustly identify age-regulated genes, suggesting that transcription factors are important in driving expression changes in aging photoreceptors. However, long, highly expressed and heavily spliced genes are also more likely to be downregulated with age, indicating that other mechanisms could contribute to expression changes at these genes. Lastly, we identify that circular RNAs (circRNAs) strongly increase during aging in photoreceptors.

CONCLUSIONS

Overall, we identified changes in gene expression in aging Drosophila photoreceptors that could account for visual senescence. Further, we show that genomic features predict these age-related changes, suggesting potential mechanisms that could be targeted to slow the rate of age-associated visual decline.

Prospective purification and characterization of Müller glia in the mouse retina regeneration assay

Reactive gliosis is an umbrella term for various glia functions in neurodegenerative diseases and upon injury. Specifically, Müller glia (MG) in some species readily regenerate retinal neurons to restore vision loss after insult, whereas mammalian MG respond by reactive gliosis-a heterogeneous response which frequently includes cell hypertrophy and proliferation. Limited regeneration has been stimulated in mammals, with a higher propensity in young MG, and in vitro compared to in vivo, but the underlying processes are unknown. To facilitate studies on the mechanisms regulating and limiting glia functions, we developed a strategy to purify glia and their progeny by fluorescence-activated cell sorting. Dual-transgenic nuclear reporter mice, which label neurons and glia with red and green fluorescent proteins, respectively, have enabled MG enrichment up to 93% purity. We applied this approach to MG in a mouse retina regeneration ex vivo assay. Combined cell size and cell cycle analysis indicates that most MG hypertrophy and a subpopulation proliferates which, over time, become even larger in cell size than the ones that do not proliferate. MG undergo timed differential genomic changes in genes controlling stemness and neurogenic competence; and glial markers are downregulated. Genes that are potentially required for, or associated with, regeneration and reactive gliosis are differentially regulated by retina explant culture time, epidermal growth factor stimulation, and animal age. Thus, MG enrichment facilitates cellular and molecular studies which, in combination with the mouse retina regeneration assay, provide an experimental approach for deciphering mechanisms that possibly regulate reactive gliosis and limit regeneration in mammals.

Design and Sampling Plan Optimization for RT-qPCR Experiments in Plants: A Case Study in Blueberry

The qPCR assay has become a routine technology in plant biotechnology and agricultural research. It is unlikely to be technically improved, but there are still challenges which center around minimizing the variability in results and transparency when reporting technical data in support of the conclusions of a study. There are a number of aspects of the pre- and post-assay workflow that contribute to variability of results. Here, through the study of the introduction of error in qPCR measurements at different stages of the workflow, we describe the most important causes of technical variability in a case study using blueberry. In this study, we found that the stage for which increasing the number of replicates would be the most beneficial depends on the tissue used. For example, we would recommend the use of more RT replicates when working with leaf tissue, while the use of more sampling (RNA extraction) replicates would be recommended when working with stems or fruits to obtain the most optimal results. The use of more qPCR replicates provides the least benefit as it is the most reproducible step. By knowing the distribution of error over an entire experiment and the costs at each step, we have developed a script to identify the optimal sampling plan within the limits of a given budget. These findings should help plant scientists improve the design of qPCR experiments and refine their laboratory practices in order to conduct qPCR assays in a more reliable-manner to produce more consistent and reproducible data.

Vision from next generation sequencing: multi-dimensional genome-wide analysis for producing gene regulatory networks underlying retinal development, aging and disease

Genomics and genetics have invaded all aspects of biology and medicine, opening uncharted territory for scientific exploration. The definition of "gene" itself has become ambiguous, and the central dogma is continuously being revised and expanded. Computational biology and computational medicine are no longer intellectual domains of the chosen few. Next generation sequencing (NGS) technology, together with novel methods of pattern recognition and network analyses, has revolutionized the way we think about fundamental biological mechanisms and cellular pathways. In this review, we discuss NGS-based genome-wide approaches that can provide deeper insights into retinal development, aging and disease pathogenesis. We first focus on gene regulatory networks (GRNs) that govern the differentiation of retinal photoreceptors and modulate adaptive response during aging. Then, we discuss NGS technology in the context of retinal disease and develop a vision for therapies based on network biology. We should emphasize that basic strategies for network construction and analyses can be transported to any tissue or cell type. We believe that specific and uniform guidelines are required for generation of genome, transcriptome and epigenome data to facilitate comparative analysis and integration of multi-dimensional data sets, and for constructing networks underlying complex biological processes. As cellular homeostasis and organismal survival are dependent on gene-gene and gene-environment interactions, we believe that network-based biology will provide the foundation for deciphering disease mechanisms and discovering novel drug targets for retinal neurodegenerative diseases.

Signaling pathway cloud regulation for in silico screening and ranking of the potential geroprotective drugs

The major challenges of aging research include absence of the comprehensive set of aging biomarkers, the time it takes to evaluate the effects of various interventions on longevity in humans and the difficulty extrapolating the results from model organisms to humans. To address these challenges we propose the in silico method for screening and ranking the possible geroprotectors followed by the high-throughput in vivo and in vitro validation. The proposed method evaluates the changes in the collection of activated or suppressed signaling pathways involved in aging and longevity, termed signaling pathway cloud, constructed using the gene expression data and epigenetic profiles of young and old patients' tissues. The possible interventions are selected and rated according to their ability to regulate age-related changes and minimize differences in the signaling pathway cloud. While many algorithmic solutions to simulating the induction of the old into young metabolic profiles in silico are possible, this flexible and scalable approach may potentially be used to predict the efficacy of the many drugs that may extend human longevity before conducting pre-clinical work and expensive clinical trials.

Intra- and inter-individual variance of gene expression in clinical studies

Background

Variance in microarray studies has been widely discussed as a critical topic on the identification of differentially expressed genes; however, few studies have addressed the influence of estimating variance.

Methodology/Principal Findings

To break intra- and inter-individual variance in clinical studies down to three levels–technical, anatomic, and individual–we designed experiments and algorithms to investigate three forms of variances. As a case study, a group of “inter-individual variable genes” were identified to exemplify the influence of underestimated variance on the statistical and biological aspects in identification of differentially expressed genes. Our results showed that inadequate estimation of variance inevitably led to the inclusion of non-statistically significant genes into those listed as significant, thereby interfering with the correct prediction of biological functions. Applying a higher cutoff value of fold changes in the selection of significant genes reduces/eliminates the effects of underestimated variance.

Conclusions/Significance

Our data demonstrated that correct variance evaluation is critical in selecting significant genes. If the degree of variance is underestimated, “noisy” genes are falsely identified as differentially expressed genes. These genes are the noise associated with biological interpretation, reducing the biological significance of the gene set. Our results also indicate that applying a higher number of fold change as the selection criteria reduces/eliminates the differences between distinct estimations of variance.

Microarray analysis of sexually dimorphic gene expression in human minor salivary glands

OBJECTIVE

We hypothesized that differential mRNA transcription between the sexes may be linked to the 9:1 female-to-male gender-related relative risk for the development of Sjögren's syndrome (SS), an autoimmune disease that leads to inflammation and dysfunction in the lachrymal and salivary glands.

SUBJECTS AND METHODS

RNA from minor salivary glands was collected from nine healthy volunteers (four men and five women) and analyzed using the Agilent 4 × 44K human microarray platform. Differential expression was confirmed by qRT-PCR.

RESULTS

Comparison of the transcriptome of minor salivary glands from normal male and female volunteers with that of salivary glands and secretory epithelia identified a number of gender, species, and tissue-specific gene expression patterns. These differences include, but are not limited to, a diverse set of genes involved in immune modulation, chemotactic control, inhibition of complement, metabolism, and neurogenesis.

CONCLUSION

Analysis of these changes provides insight into the protective and predisposing molecular factors that may be involved in the development of Sjögren's syndrome. Some of the gene changes observed in this study correlate with previously observed sexual dimorphisms in salivary gland function and also illustrate several new targets for further investigation.

Validation and characterization of DNA microarray gene expression data distribution and associated moments

Background

The data from DNA microarrays are increasingly being used in order to understand effects of different conditions, exposures or diseases on the modulation of the expression of various genes in a biological system. This knowledge is then further used in order to generate molecular mechanistic hypotheses for an organism when it is exposed to different conditions. Several different methods have been proposed to analyze these data under different distributional assumptions on gene expression. However, the empirical validation of these assumptions is lacking.

Results

Best fit hypotheses tests, moment-ratio diagrams and relationships between the different moments of the distribution of the gene expression was used to characterize the observed distributions. The data are obtained from the publicly available gene expression database, Gene Expression Omnibus (GEO) to characterize the empirical distributions of gene expressions obtained under varying experimental situations each of which providing relatively large number of samples for hypothesis testing. All data were obtained from either of two microarray platforms - the commercial Affymetrix mouse 430.2 platform and a non-commercial Rosetta/Merck one. The data from each platform were preprocessed in the same manner.

Conclusions

The null hypotheses for goodness of fit for all considered univariate theoretical probability distributions (including the Normal distribution) are rejected for more than 50% of probe sets on the Affymetrix microarray platform at a 95% confidence level, suggesting that under the tested conditions a priori assumption of any of these distributions across all probe sets is not valid. The pattern of null hypotheses rejection was different for the data from Rosetta/Merck platform with only around 20% of the probe sets failing the logistic distribution goodness-of-fit test. We find that there are statistically significant (at 95% confidence level based on the F-test for the fitted linear model) relationships between the mean and the logarithm of the coefficient of variation of the distributions of the logarithm of gene expressions. An additional novel statistically significant quadratic relationship between the skewness and kurtosis is identified. Data from both microarray platforms fail to identify with any one of the chosen theoretical probability distributions from an analysis of the l-moment ratio diagram.

Distinct signature of altered homeostasis in aging rod photoreceptors: implications for retinal diseases

Background

Advanced age contributes to clinical manifestations of many retinopathies and represents a major risk factor for age-related macular degeneration, a leading cause of visual impairment and blindness in the elderly. Rod photoreceptors are especially vulnerable to genetic defects and changes in microenvironment, and are among the first neurons to die in normal aging and in many retinal degenerative diseases. The molecular mechanisms underlying rod photoreceptor vulnerability and potential biomarkers of the aging process in this highly specialized cell type are unknown.

Methodology/Principal Findings

To discover aging-associated adaptations that may influence rod function, we have generated gene expression profiles of purified rod photoreceptors from mouse retina at young adult to early stages of aging (1.5, 5, and 12 month old mice). We identified 375 genes that showed differential expression in rods from 5 and 12 month old mouse retina compared to that of 1.5 month old retina. Quantitative RT-PCR experiments validated expression change for a majority of the 25 genes that were examined. Macroanalysis of differentially expressed genes using gene class testing and protein interaction networks revealed overrepresentation of cellular pathways that are potentially photoreceptor-specific (angiogenesis and lipid/retinoid metabolism), in addition to age-related pathways previously described in several tissue types (oxidative phosphorylation, stress and immune response).

Conclusions/Significance

Our study suggests a progressive shift in cellular homeostasis that may underlie aging-associated functional decline in rod photoreceptors and contribute to a more permissive state for pathological processes involved in retinal diseases.

Transcriptional profile analysis of RPGRORF15 frameshift mutation identifies novel genes associated with retinal degeneration

PURPOSE

To identify genes and molecular mechanisms associated with photoreceptor degeneration in a canine model of XLRP caused by an RPGR exon ORF15 microdeletion. Methods. Expression profiles of mutant and normal retinas were compared by using canine retinal custom cDNA microarrays. qRT-PCR, Western blot analysis, and immunohistochemistry (IHC) were applied to selected genes, to confirm and expand the microarray results.

RESULTS

At 7 and 16 weeks, respectively, 56 and 18 transcripts were downregulated in the mutant retinas, but none were differentially expressed (DE) at both ages, suggesting the involvement of temporally distinct pathways. Downregulated genes included the known retina-relevant genes PAX6, CHML, and RDH11 at 7 weeks and CRX and SAG at 16 weeks. Genes directly or indirectly active in apoptotic processes were altered at 7 weeks (CAMK2G, NTRK2, PRKCB, RALA, RBBP6, RNF41, SMYD3, SPP1, and TUBB2C) and 16 weeks (SLC25A5 and NKAP). Furthermore, the DE genes at 7 weeks (ELOVL6, GLOD4, NDUFS4, and REEP1) and 16 weeks (SLC25A5 and TARS2) are related to mitochondrial functions. qRT-PCR of 18 genes confirmed the microarray results and showed DE of additional genes not on the array. Only GFAP was DE at 3 weeks of age. Western blot and IHC analyses also confirmed the high reliability of the transcriptomic data.

CONCLUSIONS

Several DE genes were identified in mutant retinas. At 7 weeks, a combination of nonclassic anti- and proapoptosis genes appear to be involved in photoreceptor degeneration, whereas at both 7 and 16 weeks, the expression of mitochondria-related genes indicates that they may play a relevant role in the disease process.

Unraveling a multifactorial late-onset disease: from genetic susceptibility to disease mechanisms for age-related macular degeneration

Aging-associated neurodegenerative diseases significantly influence the quality of life of affected individuals. Genetic approaches, combined with genomic technology, have provided powerful insights into common late-onset diseases, such as age-related macular degeneration (AMD). Here, we discuss current findings on the genetics of AMD to highlight areas of rapid progress and new challenges. We also attempt to integrate available genetic and biochemical data with cellular pathways involved in aging to formulate an integrated model of AMD pathogenesis.

Age-dependent variability in gene expression in male Fischer 344 rat retina

Recent evidence suggests that older adults may be a sensitive population with regard to environmental exposure to toxic compounds. One source of this sensitivity could be an enhanced variability in response. Studies on phenotypic differences have suggested that variation in response does increase with age. However, few reports address the question of variation in gene expression as an underlying cause for increased variability of phenotypic response in the aged. In this study, we utilized global analysis to compare variation in constitutive gene expression in the retinae of young (4 months), middle-aged (11 months), and aged (23 months) Fischer 344 rats. Three hundred and forty transcripts were identified in which variance in expression increased from 4 to 23 months of age, while only 12 transcripts were found for which it decreased. Functional roles for identified genes were clustered in basic biological categories including cell communication, function, metabolism, and response to stimuli. Our data suggest that population stochastically induced variability should be considered in assessing sensitivity due to old age.

Design and analysis issues in genome-wide somatic mutation studies of cancer

The availability of the human genome sequence and progress in sequencing and bioinformatic technologies have enabled genome-wide investigation of somatic mutations in human cancers. This article briefly reviews challenges arising in the statistical analysis of mutational data of this kind. A first challenge is that of designing studies that efficiently allocate sequencing resources. We show that this can be addressed by two-stage designs and demonstrate via simulations that even relatively small studies can produce lists of candidate cancer genes that are highly informative for future research efforts. A second challenge is to distinguish mutated genes that are selected for by cancer (drivers) from mutated genes that have no role in the development of cancer and simply happened to mutate (passengers). We suggest that this question is best approached as a classification problem and discuss some of the difficulties of more traditional testing-based approaches. A third challenge is to identify biologic processes affected by the driver genes. This can be pursued by gene set analyses. These can reliably identify functional groups and pathways that are enriched for mutated genes even when the individual genes involved in those pathways or sets are not mutated at sufficient frequencies to provide conclusive evidence as drivers.

Identification of intra-group, inter-individual, and gene-specific variances in mRNA expression profiles in the rheumatoid arthritis synovial membrane

Introduction

Rheumatoid arthritis (RA) is a chronic inflammatory and destructive joint disease characterized by overexpression of pro-inflammatory/pro-destructive genes and other activating genes (for example, proto-oncogenes) in the synovial membrane (SM). The gene expression in disease is often characterized by significant inter-individual variances via specific synchronization/desynchronization of gene expression. To elucidate the contribution of the variance to the pathogenesis of disease, expression variances were tested in SM samples of RA patients, osteoarthritis (OA) patients, and normal controls (NCs).

Method

Analysis of gene expression in RA, OA, and NC samples was carried out using Affymetrix U133A/B oligonucleotide arrays, and the results were validated by real-time reverse transcription-polymerase chain reaction. For the comparison between RA and NC, 568 genes with significantly different variances in the two groups (P ≤ 0.05; Bonferroni/Holm corrected Brown-Forsythe version of the Levene test) were selected. For the comparison between RA and OA, 333 genes were selected. By means of the Kyoto Encyclopedia of Genes and Genomes, the pathways/complexes significantly affected by higher gene expression variances were identified in each group.

Results

Ten pathways/complexes significantly affected by higher gene expression variances were identified in RA compared with NC, including cytokine–cytokine receptor interactions, the transforming growth factor-beta pathway, and anti-apoptosis. Compared with OA, three pathways with significantly higher variances were identified in RA (for example, B-cell receptor signaling and vascular endothelial growth factor signaling). Functionally, the majority of the identified pathways are involved in the regulation of inflammation, proliferation, cell survival, and angiogenesis.

Conclusion

In RA, a number of disease-relevant or even disease-specific pathways/complexes are characterized by broad intra-group inter-individual expression variances. Thus, RA pathogenesis in different individuals may depend to a lesser extent on common alterations of the expression of specific key genes, and rather on individual-specific alterations of different genes resulting in common disturbances of key pathways.

Microarray-based gene expression analysis during retinal maturation of albino rats

BACKGROUND

In recent years, the rat has become a commonly-used animal model for the study of retinal diseases. Similar to other tissues, the retina undergoes significant functional changes during maturation. Aiming to gain knowledge on additional aspects of retinal maturation, we performed gene expression and histological analyses of the rat retina during maturation.

METHODS

Rat retinas were dissected at three time points. Histological examination of the samples was performed, and the expression levels of retinal genes were evaluated using the rat whole-genome microarray system. Quantitative real-time PCR analysis was used to validate selected expression patterns. Various statistical and bioinformatic tools were used to identify differentially expressed genes.

RESULTS

The microarray analysis revealed a relatively high number of highly expressed non-annotated genes. We identified 603 differentially expressed genes, which were grouped into six clusters based on changes in expression levels during the first 20 weeks of life. A bioinformatic analysis of these clusters revealed sets of genes encoding proteins with functions that are likely to be relevant to retinal maturation (potassium, sodium, calcium, and chloride channels, synaptic vesicle transport, and axonogenesis). The histological analysis revealed a significant reduction of outer nuclear layer thickness and retinal ganglion cell number during maturation.

CONCLUSIONS

These data, taken together with our previously reported electrophysiological data, contribute to our understanding of the retinal maturation processes of this widely-used animal model.

Limited dynamic range of immune response gene expression observed in healthy blood donors using RT-PCR

The use of quantitative gene expression analysis for the diagnosis, prognosis, and monitoring of disease requires the ability to distinguish pathophysiological changes from natural variations. To characterize these variations in apparently healthy subjects, quantitative real-time PCR was used to measure various immune response genes in whole blood collected from blood bank donors. In a single-time-point study of 131 donors, of 48 target genes, 43 were consistently expressed and 34 followed approximately log-normal distribution. Most transcripts showed a limited dynamic range of expression across subjects. Specifically, 36 genes had standard deviations (SDs) of 0.44 to 0.79 cycle threshold (C(T)) units, corresponding to less than a 3-fold variation in expression. Separately, a longitudinal study of 8 healthy individuals demonstrated a total dynamic range (> 2 standard error units) of 2- to 4-fold in most genes. In contrast, a study of whole blood gene expression in 6 volunteers injected with LPS showed 15 genes changing in expression 10- to 90-fold within 2 to 5 h and returning to within normal range within 21 hours. This work demonstrates that (1) the dynamic range of expression of many immune response genes is limited among healthy subjects; (2) expression levels for most genes analyzed are approximately log-normally distributed; and (3) individuals exposed to an infusion of bacterial endotoxin (lipopolysaccharide), show gene expression profiles that can be readily distinguished from those of a healthy population. These results suggest that normal reference ranges can be established for gene expression assays, providing critical standards for the diagnosis and management of disease.

Variability of DNA Microarray Gene Expression Profiles in Cultured Rat Primary Hepatocytes

DNA microarray is a powerful tool in biomedical research. However, transcriptomic profiling using DNA microarray is subject to many variations including biological variability. To evaluate the different sources of variation in mRNA gene expression profiles, gene expression profiles were monitored using the Affymetrix RatTox U34 arrays in cultured primary hepatocytes derived from six rats over a 26 hour period at 6 time points (0h, 2h, 5h, 8h, 14h and 26h) with two replicate arrays at each time point for each animal. In addition, the impact of sample size on the variability of differentially expressed gene lists and the consistency of biological responses were also investigated. Excellent intra-animal reproducibility was obtained at all time points with 0 out of 370 present probe sets across all time points showing significant difference between the 2 replicate arrays (3-way ANOVA, p ≤ 0.0001). However, large inter-animal biological variation in mRNA expression profiles was observed with 337 out of 370 present probe sets showing significant differences among 6 animals (3-way ANOVA, p ≤ 0.05). Principal Component Analysis (PCA) revealed that time effect (PC1) in this data set accounted for 47.4% of total variance indicating the dynamics of transcriptomics. The second and third largest effects came from animal difference, which accounted for 16.9% (PC2 and PC3) of the total variance. The reproducibility of gene lists and their functional classification was declined considerably when the sample size was decreased. Overall, our results strongly support that there is significant inter-animal variability in the time-course gene expression profiles, which is a confounding factor that must be carefully evaluated to correctly interpret microarray gene expression studies. The consistency of the gene lists and their biological functional classification are also sensitive to sample size with the reproducibility decreasing considerably under small sample size.

Defining the human macula transcriptome and candidate retinal disease genes using EyeSAGE

PURPOSE

To develop large-scale, high-throughput annotation of the human macula transcriptome and to identify and prioritize candidate genes for inherited retinal dystrophies, based on ocular-expression profiles using serial analysis of gene expression (SAGE).

METHODS

Two human retina and two retinal pigment epithelium (RPE)/choroid SAGE libraries made from matched macula or midperipheral retina and adjacent RPE/choroid of morphologically normal 28- to 66-year-old donors and a human central retina longSAGE library made from 41- to 66-year-old donors were generated. Their transcription profiles were entered into a relational database, EyeSAGE, including microarray expression profiles of retina and publicly available normal human tissue SAGE libraries. EyeSAGE was used to identify retina- and RPE-specific and -associated genes, and candidate genes for retina and RPE disease loci. Differential and/or cell-type specific expression was validated by quantitative and single-cell RT-PCR.

RESULTS

Cone photoreceptor-associated gene expression was elevated in the macula transcription profiles. Analysis of the longSAGE retina tags enhanced tag-to-gene mapping and revealed alternatively spliced genes. Analysis of candidate gene expression tables for the identified Bardet-Biedl syndrome disease gene (BBS5) in the BBS5 disease region table yielded BBS5 as the top candidate. Compelling candidates for inherited retina diseases were identified.

CONCLUSIONS

The EyeSAGE database, combining three different gene-profiling platforms including the authors' multidonor-derived retina/RPE SAGE libraries and existing single-donor retina/RPE libraries, is a powerful resource for definition of the retina and RPE transcriptomes. It can be used to identify retina-specific genes, including alternatively spliced transcripts and to prioritize candidate genes within mapped retinal disease regions.

Human lung project: evaluating variance of gene expression in the human lung

Nondiseased tissue is an important reference for microarray studies of pulmonary disease. We obtained 23 single lungs from multiorgan donors at time of procurement. Donors varied in age, sex, smoking history, and ethnicity. Lungs were dissected into upper and lower lobe peripheral sections for RNA extraction. Microarray analysis was performed using Affymetrix Hu-133 Plus 2.0 arrays. We observed that the relative variability of gene expression increased rapidly from technical (lowest), to regional, to population (highest). In addition, age and sex have measurable effects on gene expression. Gene expression variability is heterogeneously distributed among biologic categories. We conclude that gene expression variability is greater between individuals than within individuals and that population variability is the most important factor in the study design of microarray experiments of the human lung. Classes of genes with high population variability are biologically important and provide a novel perspective into lung physiology and pathobiology. Our study represents the first comprehensive analysis of nondiseased lung tissue. The generation of this robust dataset has important implications for the design and implementation of future comparative expression analysis with pulmonary disease states.

Thrombospondin 1, thrombospondin 2 and the eye

Thrombospondin 1 and thrombospondin 2 (TSP1 and TSP2), which comprise the subgroup A thrombospondins, are matricellular proteins. As matricellular proteins, they modulate interactions between cells and the cellular environment, regulate cell adhesion and typically are expressed during tissue formative processes. In general, TSP1 and TSP2 counter angiogenesis (including tumour angiogenesis) and play important but contrasting roles during cutaneous repair. The two proteins are involved in development, including that of the eye, although evidence suggests that they have their greatest impact during tissue production in the adult. In the normal adult eye, they tend to be found at sites of ongoing matrix synthesis or cell-matrix interactions. At these sites, the two proteins possibly influence cellular differentiation and/or basement membrane deposition. TSP1 is also present in the intraocular fluids and drainage pathway, where it may function in maintaining the anti-angiogenic environment and in intraocular pressure control, respectively. TSP1 could also be involved in ocular immune privilege. Unlike in skin wounds, where TSP1 is derived from the blood and is present only in the early phases of repair, ocular tissue damage appears to lead to protacted TSP1 synthesis by local cells. This response might help suppress angiogenesis in the transparent tissues of the eye and so lessen visual axis opacification following injury. However, TSP2, which is also produced by damaged ophthalmic tissue and may be especially important in matrix organisation, seems to augment contraction in anomalous intraocular fibrosis. Elucidating the roles of TSP1 and TSP2 in ocular physiology and pathobiology may lead to improved therapies for neovascular, neoplastic, reparative and other ophthalmic diseases.

Proteinase and growth factor alterations revealed by gene microarray analysis of human diabetic corneas

PURPOSE

To identify proteinases and growth factors abnormally expressed in human corneas of donors with diabetic retinopathy (DR), additional to previously described matrix metalloproteinase (MMP)-10 and -3 and insulin-like growth factor (IGF)-I.

METHODS

RNA was isolated from 35 normal, diabetic, and DR autopsy human corneas ex vivo or after organ culture. Amplified cRNA was analyzed using 22,000-gene microarrays (Agilent Technologies, Palo Alto, CA). Gene expression in each diabetic corneal cRNA was assessed against pooled cRNA from 7 to 9 normal corneas. Select differentially expressed genes were validated by quantitative real-time RT-PCR (QPCR) and immunohistochemistry. Organ cultures were treated with a cathepsin inhibitor, cystatin C, or MMP-10.

RESULTS

More than 100 genes were upregulated and 2200 were downregulated in DR corneas. Expression of cathepsin F and hepatocyte growth factor (HGF) genes was increased in ex vivo and organ-cultured DR corneas compared with normal corneas. HGF receptor c-met, fibroblast growth factor (FGF)-3, its receptor FGFR3, tissue inhibitor of metalloproteinase (TIMP)-4, laminin alpha4 chain, and thymosin beta(4) genes were downregulated. The data were corroborated by QPCR and immunohistochemistry analyses; main changes of these components occurred in corneal epithelium. In organ-cultured DR corneas, cystatin C increased laminin-10 and integrin alpha(3)beta(1), whereas in normal corneas MMP-10 decreased laminin-10 and integrin alpha(3)beta(1) expression.

CONCLUSIONS

Elevated cathepsin F and the ability of its inhibitor to produce a more normal phenotype in diabetic corneas suggest increased proteolysis in these corneas. Proteinase changes may result from abnormalities of growth factors, such as HGF and FGF-3, in DR corneas. Specific modulation of proteinases and growth factors could reduce diabetic corneal epitheliopathy.

Distinct gene expression profiles and reduced JNK signaling in retinitis pigmentosa caused by RP1 mutations

To understand the mechanisms underlying autosomal dominant progressive retinitis pigmentosa (RP) caused by the mutations of the RP1 gene and to identify molecules that play roles in the early disease process, we used Affymetrix U74Av2 microarrays to compare the gene expression profiles of retinas from Rp1-/- and Rp1+/+ mice at postnatal days (P) 7, 10, 14, 18 and 21. These profiles were independently verified by comparison with results of retinal serial analysis of gene expression, U74Av2 array studies of mouse retinas, real-time PCR and in situ hybridization. We found that the disruption of Rp1 significantly affected the expression of multiple clusters of genes whose products were involved in diverse biological pathways. The molecular responses to the disruption of Rp1 changed dramatically during development and were distinct from responses to the disruption of photoreceptor transcription factors (Crx-/- or Nrl-/-) and a phototransduction molecule (Pde6brd1). We found specific alterations of gene expression in the c-Jun N-terminal kinase (JNK) signaling cascades. Western analysis confirmed that the phosphorylation of key members in the JNK signaling cascades (i.e. JNK1, JNK2, MAP2, MKK4 and c-Jun) is reduced, whereas phospho-ERK and phospho-p38 are unchanged, in Rp1-/- retinas at P18-21. Immunostaining demonstrated that, like Rp1, phospho-JNKs and phospho-MAP2 are present in outer segments of photoreceptors. Our studies reveal unique molecular phenotypes in multiple biological pathways and the specific reduction of JNK signaling cascades in RP1 diseases, and suggest that RP1, a doublecortin-containing microtubule associated protein, and JNK signaling cascades play integral roles in photoreceptor development and maintenance. Our studies further suggest JNK-related therapeutic strategies for RP1 diseases.

Differential gene expression in anatomical compartments of the human eye

DNA microarrays (representing approximately 30,000 human genes) were used to analyze gene expression in six different human eye compartments, revealing candidate genes for diseases affecting the cornea, lens and retina.

Background

The human eye is composed of multiple compartments, diverse in form, function, and embryologic origin, that work in concert to provide us with our sense of sight. We set out to systematically characterize the global gene expression patterns that specify the distinctive characteristics of the various eye compartments.

Results

We used DNA microarrays representing approximately 30,000 human genes to analyze gene expression in the cornea, lens, iris, ciliary body, retina, and optic nerve. The distinctive patterns of expression in each compartment could be interpreted in relation to the physiology and cellular composition of each tissue. Notably, the sets of genes selectively expressed in the retina and in the lens were particularly large and diverse. Genes with roles in immune defense, particularly complement components, were expressed at especially high levels in the anterior segment tissues. We also found consistent differences between the gene expression patterns of the macula and peripheral retina, paralleling the differences in cell layer densities between these regions. Based on the hypothesis that genes responsible for diseases that affect a particular eye compartment are likely to be selectively expressed in that compartment, we compared our gene expression signatures with genetic mapping studies to identify candidate genes for diseases affecting the cornea, lens, and retina.

Conclusion

Through genome-scale gene expression profiling, we were able to discover distinct gene expression 'signatures' for each eye compartment and identified candidate disease genes that can serve as a reference database for investigating the physiology and pathophysiology of the eye.

The Retinome - defining a reference transcriptome of the adult mammalian retina/retinal pigment epithelium

BACKGROUND

The mammalian retina is a valuable model system to study neuronal biology in health and disease. To obtain insight into intrinsic processes of the retina, great efforts are directed towards the identification and characterization of transcripts with functional relevance to this tissue.

RESULTS

With the goal to assemble a first genome-wide reference transcriptome of the adult mammalian retina, referred to as the retinome, we have extracted 13,037 non-redundant annotated genes from nearly 500,000 published datasets on redundant retina/retinal pigment epithelium (RPE) transcripts. The data were generated from 27 independent studies employing a wide range of molecular and biocomputational approaches. Comparison to known retina-/RPE-specific pathways and established retinal gene networks suggest that the reference retinome may represent up to 90% of the retinal transcripts. We show that the distribution of retinal genes along the chromosomes is not random but exhibits a higher order organization closely following the previously observed clustering of genes with increased expression.

CONCLUSION

The genome wide retinome map offers a rational basis for selecting suggestive candidate genes for hereditary as well as complex retinal diseases facilitating elaborate studies into normal and pathological pathways. To make this unique resource freely available we have built a database providing a query interface to the reference retinome 1.