Identifying RNA Biomarkers and Molecular Pathways Involved in Multiple Subtypes of Uveitis.

PURPOSE: Uveitis is a heterogeneous collection of diseases. We tested the hypothesis that despite the diversity of uveitides, there could be common mechanisms shared by multiple subtypes, and that evidence of these common mechanisms may be detected as gene expression profiles in whole blood. DESIGN: Cohort study. METHODS: Ninety subjects with uveitis including axial spondyloarthritis (n = 17), sarcoidosis (n = 13), inflammatory bowel disease (n = 12), tubulointerstitial nephritis with uveitis (n = 10), or idiopathic uveitis (n = 38) as well as 18 healthy controls were enrolled, predominantly at Oregon Health & Science University. RNA-Seq data generated from peripheral, whole blood identified 19,859 unique transcripts. We analyzed gene expression pathways via Kyoto Encyclopedia of Genes and Genomes and Gene Ontology (GO). We validated our list of upregulated genes by comparison to a previously published study on peripheral blood gene expression among 50 subjects with diverse forms of uveitis. RESULTS: Both the Kyoto Encyclopedia of Genes and Genomes and GO analysis identified multiple shared pathways or GO terms with a P value of <.0001. Almost all pathways related to the immune response and/or response to an infection. A total of 119 individual transcripts were upregulated by at least 1.5-fold and false discovery rate <.05, and 61 were downregulated by similar criteria. Comparing mRNA from our study with a false discovery rate <.05 and the prior report, we identified 10 common gene transcripts: ICAM1, IL15RA, IL15, IRF1, IL10RB, GSK3A, TYK2, MEF2A, MEF2B, and MEF2D. CONCLUSIONS: Many forms of uveitis share overlapping mechanisms. These data support the concept that a single therapeutic approach could benefit diverse forms of this disease.

Revising the Diagnosis of Idiopathic Uveitis by Peripheral Blood Transcriptomics.

PURPOSE: To test the hypothesis that idiopathic uveitis can be categorized into subtypes based on gene expression from blood. DESIGN: Case control study. METHODS: We applied RNA-Seq to peripheral blood from patients with uveitis associated with 1 of 4 systemic diseases, including axial spondyloarthritis (n = 17), sarcoidosis (n = 13), inflammatory bowel disease (n = 12), tubulo-interstitial nephritis with uveitis (n = 10), or idiopathic uveitis (n = 38) as well as 18 healthy control subjects evaluated predominantly at Oregon Health and Science University. A high-dimensional negative binomial regression model implemented in the edgeR R package compared each disease group with the control subjects. The 20 most distinctive genes for each diagnosis were extracted. Of 80 genes, there were 75 unique genes. A classification algorithm was developed by fitting a gradient boosting tree with 5-fold cross-validation. Messenger RNA from subjects with idiopathic uveitis were analyzed to see if any fit clinically and by gene expression pattern with one of the diagnosable entities. RESULTS: For uveitis associated with a diagnosable systemic disease, gene expression profiling achieved an overall accuracy of 85% (balanced average of sensitivity plus specificity, P < .001). Although most patients with idiopathic uveitis presumably have none of these 4 associated systemic diseases, gene expression profiles helped to reclassify 11 of 38 subjects. CONCLUSIONS: Peripheral blood gene expression profiling is a potential adjunct in accurate differential diagnosis of the cause of uveitis. Validation of these results and characterization of the gene expression profile from additional discrete diagnoses could enhance the value of these observations.

Phenotypic and gene expression features associated with variation in chronic ethanol consumption in heterogeneous stock collaborative cross mice.

Of the more than 100 studies that have examined relationships between excessive ethanol consumption and the brain transcriptome, few rodent studies have examined chronic consumption. Heterogeneous stock collaborative cross mice freely consumed ethanol vs. water for 3 months. Transcriptional differences were examined for the central nucleus of the amygdala, a brain region known to impact ethanol preference. Early preference was modestly predictive of final preference and there was significant escalation of preference in females only. Genes significantly correlated with female preference were enriched in annotations for the primary cilium and extracellular matrix. A single module in the gene co-expression network was enriched in genes with an astrocyte annotation. The key hub node was the master regulator, orthodenticle homeobox 2 (Otx2). These data support an important role for the extracellular matrix, primary cilium and astrocytes in ethanol preference and consumption differences among individual female mice of a genetically diverse population.

Chronic Voluntary Ethanol Drinking in Cynomolgus Macaques Elicits Gene Expression Changes in Prefrontal Cortical Area 46.

BACKGROUND: Genome-wide profiling to examine brain transcriptional features associated with excessive ethanol (EtOH) consumption has been applied to a variety of species including rodents, nonhuman primates (NHPs), and humans. However, these data were obtained from cross-sectional samples which are particularly vulnerable to individual variation when obtained from small outbred populations typical of human and NHP studies. In the current study, a novel within-subject design was used to examine the effects of voluntary EtOH consumption on prefrontal cortex (PFC) gene expression in a NHP model. METHODS: Two cohorts of cynomolgus macaques (n = 23) underwent a schedule-induced polydipsia procedure to establish EtOH self-administration followed by 6 months of daily open access to EtOH (4% w/v) and water. Individual daily EtOH intakes ranged from an average of 0.7 to 3.7 g/kg/d. Dorsal lateral PFC area 46 (A46) brain biopsies were collected in EtOH-naïve and control monkeys; contralateral A46 biopsies were collected from the same monkeys following the 6 months of fluid consumption. Gene expression changes were assessed using RNA-Seq paired analysis, which allowed for correction of individual baseline differences in gene expression. RESULTS: A total of 675 genes were significantly down-regulated following EtOH consumption; these were functionally enriched for immune response, cell adhesion, plasma membrane, and extracellular matrix. A total of 567 genes that were up-regulated following EtOH consumption were enriched in microRNA target sites and included target sites associated with Toll-like receptor pathways. The differentially expressed genes were also significantly enriched in transcription factor binding sites. CONCLUSIONS: The data presented here are the first to use a longitudinal biopsy strategy to examine how chronic EtOH consumption affects gene expression in the primate PFC. Prominent effects were seen in both cell adhesion and neuroimmune pathways; the latter contained both pro- and antiinflammatory genes. The data also indicate that changes in miRNAs and transcription factors may be important epigenetic regulators of EtOH consumption.

Cross-species molecular dissection across alcohol behavioral domains.

This review summarizes the proceedings of a symposium presented at the "Alcoholism and Stress: A Framework for Future Treatment Strategies" conference held in Volterra, Italy on May 9-12, 2017. Psychiatric diseases, including alcohol-use disorders (AUDs), are influenced through complex interactions of genes, neurobiological pathways, and environmental influences. A better understanding of the common neurobiological mechanisms underlying an AUD necessitates an integrative approach, involving a systematic assessment of diverse species and phenotype measures. As part of the World Congress on Stress and Alcoholism, this symposium provided a detailed account of current strategies to identify mechanisms underlying the development and progression of AUDs. Dr. Sean Farris discussed the integration and organization of transcriptome and postmortem human brain data to identify brain regional- and cell type-specific differences related to excessive alcohol consumption that are conserved across species. Dr. Brien Riley presented the results of a genome-wide association study of DSM-IV alcohol dependence; although replication of genetic associations with alcohol phenotypes in humans remains challenging, model organism studies show that COL6A3, KLF12, and RYR3 affect behavioral responses to ethanol, and provide substantial evidence for their role in human alcohol-related traits. Dr. Rob Williams expanded upon the systematic characterization of extensive genetic-genomic resources for quantifying and clarifying phenotypes across species that are relevant to precision medicine in human disease. The symposium concluded with Dr. Robert Hitzemann's description of transcriptome studies in a mouse model selectively bred for high alcohol ("binge-like") consumption and a non-human primate model of long-term alcohol consumption. Together, the different components of this session provided an overview of systems-based approaches that are pioneering the experimental prioritization and validation of novel genes and gene networks linked with a range of behavioral phenotypes associated with stress and AUDs.

On the relationships in rhesus macaques between chronic ethanol consumption and the brain transcriptome.

This is the first description of the relationship between chronic ethanol self-administration and the brain transcriptome in a non-human primate (rhesus macaque). Thirty-one male animals self-administered ethanol on a daily basis for over 12 months. Gene transcription was quantified with RNA-Seq in the central nucleus of the amygdala (CeA) and cortical Area 32. We constructed coexpression and cosplicing networks, and we identified areas of preservation and areas of differentiation between regions and network types. Correlations between intake and transcription included largely distinct gene sets and annotation categories across brain regions and between expression and splicing; positive and negative correlations were also associated with distinct annotation groups. Membrane, synaptic and splicing annotation categories were over-represented in the modules (gene clusters) enriched in positive correlations (CeA); our cosplicing analysis further identified the genes affected only at the exon inclusion level. In the CeA coexpression network, we identified Rab6b, Cdk18 and Igsf21 among the intake-correlated hubs, while in the Area 32, we identified a distinct hub set that included Ppp3r1 and Myeov2. Overall, the data illustrate that excessive ethanol self-administration is associated with broad expression and splicing mechanisms that involve membrane and synapse genes.

Anti-Müllerian hormone is produced heterogeneously in primate preantral follicles and is a potential biomarker for follicle growth and oocyte maturation in vitro.

PURPOSE: The main goals of this study were to investigate the expression of anti-Müllerian hormone (AMH) and its receptor (AMHR2) during follicular development in primates, and to evaluate the potential of AMH as a biomarker for follicle growth and oocyte maturation in vitro. METHODS: The mRNA and protein expression of AMH and AMHR2 were determined using isolated follicles and ovarian sections from rhesus macaques (n = 4) by real-time PCR and immunohistochemistry, respectively. Isolated secondary follicles were cultured individually. Follicle growth and media AMH concentrations were assessed by ELISA. The mRNA expression profiles, obtained from RNA sequencing, of in vitro- and in vivo-developed antral follicles were compared. Secondary follicles from additional animals (n = 35) were cultured. Follicle growth, oocyte maturation, and media AMH concentrations were evaluated for forecasting follicular development in vitro by AMH levels. RESULTS: AMH immunostaining was heterogeneous in the population of preantral follicles that were also stained for AMHR2. The mRNA expression profiles were comparable between in vivo- and in vitro-developed follicles. AMH levels produced by growing follicles were higher than those of nongrowing follicles in culture. With a cutoff value of 1.40 ng/ml, 85 % of nongrowing follicles could be identified while eliminating only 5 % of growing follicles. Growing follicles that generated metaphase II-stage oocytes secreted greater amounts of AMH than did those yielding immature germinal vesicle-stage oocytes. CONCLUSIONS: AMH, co-expressed with AMHR2, was produced heterogeneously by preantral follicles in macaques with levels correlated positively with follicle growth and oocyte maturation. AMH may serve as a biomarker for primate follicular development in vitro.

Splicing landscape of the eight collaborative cross founder strains.

BACKGROUND: The Collaborative Cross (CC) is a large panel of genetically diverse recombinant inbred mouse strains specifically designed to provide a systems genetics resource for the study of complex traits. In part, the utility of the CC stems from the extensive genome-wide annotations of founder strain sequence and structural variation. Still missing, however, are transcriptome-specific annotations of the CC founder strains that could further enhance the utility of this resource. RESULTS: We provide a comprehensive survey of the splicing landscape of the 8 CC founder strains by leveraging the high level of alternative splicing within the brain. Using deep transcriptome sequencing, we found that a majority of the splicing landscape is conserved among the 8 strains, with ~65% of junctions being shared by at least 2 strains. We, however, found a large number of potential strain-specific splicing events as well, with an average of ~3000 and ~500 with ≥3 and ≥10 sequence read coverage, respectively, within each strain. To better understand strain-specific splicing within the CC founder strains, we defined criteria for and identified high-confidence strain-specific splicing events. These splicing events were defined as exon-exon junctions 1) found within only one strain, 2) with a read coverage ≥10, and 3) defined by a canonical splice site. With these criteria, a total of 1509 high-confidence strain-specific splicing events were identified, with the majority found within two of the wild-derived strains, CAST and PWK. Strikingly, the overwhelming majority, 94%, of these strain-specific splicing events are not yet annotated. Strain-specific splicing was also located within genomic regions recently reported to be over- and under-represented within CC populations. CONCLUSIONS: Phenotypic characterization of CC populations is increasing; thus these results will not only aid in further elucidating the transcriptomic architecture of the individual CC founder strains, but they will also help in guiding the utilization of the CC populations in the study of complex traits. This report is also the first to establish guidelines in defining and identifying strain-specific splicing across different mouse strains.

Evaluating gene expression in C57BL/6J and DBA/2J mouse striatum using RNA-Seq and microarrays.

C57BL/6J (B6) and DBA/2J (D2) are two of the most commonly used inbred mouse strains in neuroscience research. However, the only currently available mouse genome is based entirely on the B6 strain sequence. Subsequently, oligonucleotide microarray probes are based solely on this B6 reference sequence, making their application for gene expression profiling comparisons across mouse strains dubious due to their allelic sequence differences, including single nucleotide polymorphisms (SNPs). The emergence of next-generation sequencing (NGS) and the RNA-Seq application provides a clear alternative to oligonucleotide arrays for detecting differential gene expression without the problems inherent to hybridization-based technologies. Using RNA-Seq, an average of 22 million short sequencing reads were generated per sample for 21 samples (10 B6 and 11 D2), and these reads were aligned to the mouse reference genome, allowing 16,183 Ensembl genes to be queried in striatum for both strains. To determine differential expression, 'digital mRNA counting' is applied based on reads that map to exons. The current study compares RNA-Seq (Illumina GA IIx) with two microarray platforms (Illumina MouseRef-8 v2.0 and Affymetrix MOE 430 2.0) to detect differential striatal gene expression between the B6 and D2 inbred mouse strains. We show that by using stringent data processing requirements differential expression as determined by RNA-Seq is concordant with both the Affymetrix and Illumina platforms in more instances than it is concordant with only a single platform, and that instances of discordance with respect to direction of fold change were rare. Finally, we show that additional information is gained from RNA-Seq compared to hybridization-based techniques as RNA-Seq detects more genes than either microarray platform. The majority of genes differentially expressed in RNA-Seq were only detected as present in RNA-Seq, which is important for studies with smaller effect sizes where the sensitivity of hybridization-based techniques could bias interpretation.

High throughput sequencing in mice: a platform comparison identifies a preponderance of cryptic SNPs.

BACKGROUND: Allelic variation is the cornerstone of genetically determined differences in gene expression, gene product structure, physiology, and behavior. However, allelic variation, particularly cryptic (unknown or not annotated) variation, is problematic for follow up analyses. Polymorphisms result in a high incidence of false positive and false negative results in hybridization based analyses and hinder the identification of the true variation underlying genetically determined differences in physiology and behavior. Given the proliferation of mouse genetic models (e.g., knockout models, selectively bred lines, heterogeneous stocks derived from standard inbred strains and wild mice) and the wealth of gene expression microarray and phenotypic studies using genetic models, the impact of naturally-occurring polymorphisms on these data is critical. With the advent of next-generation, high-throughput sequencing, we are now in a position to determine to what extent polymorphisms are currently cryptic in such models and their impact on downstream analyses. RESULTS: We sequenced the two most commonly used inbred mouse strains, DBA/2J and C57BL/6J, across a region of chromosome 1 (171.6 - 174.6 megabases) using two next generation high-throughput sequencing platforms: Applied Biosystems (SOLiD) and Illumina (Genome Analyzer). Using the same templates on both platforms, we compared realignments and single nucleotide polymorphism (SNP) detection with an 80 fold average read depth across platforms and samples. While public datasets currently annotate 4,527 SNPs between the two strains in this interval, thorough high-throughput sequencing identified a total of 11,824 SNPs in the interval, including 7,663 new SNPs. Furthermore, we confirmed 40 missense SNPs and discovered 36 new missense SNPs. CONCLUSION: Comparisons utilizing even two of the best characterized mouse genetic models, DBA/2J and C57BL/6J, indicate that more than half of naturally-occurring SNPs remain cryptic. The magnitude of this problem is compounded when using more divergent or poorly annotated genetic models. This warrants full genomic sequencing of the mouse strains used as genetic models.

Rhesus macaque rhadinovirus-associated non-Hodgkin lymphoma: animal model for KSHV-associated malignancies.

Rhesus macaque rhadinovirus (RRV) is closely related to Kaposi sarcoma-associated herpesvirus (KSHV) and is associated with the development of B-cell hyperplasia and persistent lymphadenopathy resembling multicentric Castleman disease in rhesus macaques (RMs) coinfected with simian immunodeficiency virus (SIV). Here we investigated whether RMs experimentally infected with SIV and RRV can develop other disease manifestations observed in HIV- and KSHV-infected patients. As reported earlier, inoculation of SIV-infected RMs with RRV results in persistent RRV infection, whereas immunocompetent animals infected with RRV exhibit viremia 2 weeks after infection, followed by a period of no virus detection until they are subsequently made immunodeficient by SIV infection. A subset of animals developed abnormal cellular proliferations characterized as extranodal lymphoma and a proliferative mesenchymal lesion. In situ hybridization and immunohistochemistry analysis indicate RRV is present in both malignancies, and DNA microarray analysis detected viral interleukin-6 (vIL-6) and viral FLICE-like inhibitory protein (vFLIP) transcripts. Reverse-transcriptase polymerase chain reaction analysis confirmed vIL-6 and vFLIP expression, and that of RRV open reading frames 72 and 73, homologs of KSHV open reading frames shown to be expressed in primary effusion lymphoma. These data support the utility of the RRV-/SIV-infected RM as an excellent animal model to investigate KSHV-like pathogenesis.

DNA multiplex hybridization on microarrays and thermodynamic stability in solution: a direct comparison.

Hybridization intensities of 30 distinct short duplex DNAs measured on spotted microarrays, were directly compared with thermodynamic stabilities measured in solution. DNA sequences were designed to promote formation of perfect match, or hybrid duplexes containing tandem mismatches. Thermodynamic parameters DeltaH degrees , DeltaS degrees and DeltaG degrees of melting transitions in solution were evaluated directly using differential scanning calorimetry. Quantitative comparison with results from 63 multiplex microarray hybridization experiments provided a linear relationship for perfect match and most mismatch duplexes. Examination of outliers suggests that both duplex length and relative position of tandem mismatches could be important factors contributing to observed deviations from linearity. A detailed comparison of measured thermodynamic parameters with those calculated using the nearest-neighbor model was performed. Analysis revealed the nearest-neighbor model generally predicts mismatch duplexes to be less stable than experimentally observed. Results also show the relative stability of a tandem mismatch is highly dependent on the identity of the flanking Watson-Crick (w/c) base pairs. Thus, specifying the stability contribution of a tandem mismatch requires consideration of the sequence identity of at least four base pair units (tandem mismatch and flanking w/c base pairs). These observations underscore the need for rigorous evaluation of thermodynamic parameters describing tandem mismatch stability.

Multiplex SNP discrimination.

Multiplex hybridization reactions of perfectly matched duplexes and duplexes containing a single basepair mismatch (SNPs) were investigated on DNA microarrays. Effects of duplex length, G-C percentage, and relative position of the SNP on duplex hybridization and SNP resolution were determined. Our theoretical model of multiplex hybridization accurately predicts observed results and implicates target concentration as a critical variable in multiplex SNP detection.

A time course microarray study of gene expression in the mouse lacrimal gland after acute corneal trauma.

PURPOSE: To investigate the effect of corneal trauma on gene expression in the lacrimal gland and to assess how many genes and what specific genes are regulated in response to corneal trauma. METHODS: A mouse model with acute corneal trauma was created with a chemical burn to the cornea with silver nitrate. Sixty-four female BALB/c mice at 12 weeks of age were randomly divided into eight groups, eight mice per group. The corneas of four mice in each group were bilaterally cauterized with silver nitrate, and another four time-matched mice were used as the control. The total RNA of the lacrimal gland was then extracted, at eight time points--0.5, 1, 3, 8, 24, 72, 120, and 360 hours--after the corneal burn, and gene expression was examined with using cDNA microarray technology. RESULTS: Evaluation of 15,065 genes with multiple array replications showed significantly altered expression in 3,799 genes at one or more of the eight time points. Of those, 1,528 were known genes and 2,271 were unknown. The analysis of known genes showed broad and long-lasting gene suppression in most functional gene groups, including housekeeping, energy metabolism, protein degradation, DNA and protein synthesis, and apoptosis-associated genes. Heat shock genes were upregulated beginning at the 8-hour time point, indicating a stress response. CONCLUSIONS: This study demonstrates that corneal trauma has profound effects on the regulation of gene expression in the lacrimal gland and may provide genetic evidence for a cornea-to-lacrimal gland feedback mechanism in dry eye.