Blood transcriptomic comparison of individuals with and without autism spectrum disorder: A combined-samples mega-analysis.

Blood-based microarray studies comparing individuals affected with autism spectrum disorder (ASD) and typically developing individuals help characterize differences in circulating immune cell functions and offer potential biomarker signal. We sought to combine the subject-level data from previously published studies by mega-analysis to increase the statistical power. We identified studies that compared ex vivo blood or lymphocytes from ASD-affected individuals and unrelated comparison subjects using Affymetrix or Illumina array platforms. Raw microarray data and clinical meta-data were obtained from seven studies, totaling 626 affected and 447 comparison subjects. Microarray data were processed using uniform methods. Covariate-controlled mixed-effect linear models were used to identify gene transcripts and co-expression network modules that were significantly associated with diagnostic status. Permutation-based gene-set analysis was used to identify functionally related sets of genes that were over- and under-expressed among ASD samples. Our results were consistent with diminished interferon-, EGF-, PDGF-, PI3K-AKT-mTOR-, and RAS-MAPK-signaling cascades, and increased ribosomal translation and NK-cell related activity in ASD. We explored evidence for sex-differences in the ASD-related transcriptomic signature. We also demonstrated that machine-learning classifiers using blood transcriptome data perform with moderate accuracy when data are combined across studies. Comparing our results with those from blood-based studies of protein biomarkers (e.g., cytokines and trophic factors), we propose that ASD may feature decoupling between certain circulating signaling proteins (higher in ASD samples) and the transcriptional cascades which they typically elicit within circulating immune cells (lower in ASD samples). These findings provide insight into ASD-related transcriptional differences in circulating immune cells. © 2016 Wiley Periodicals, Inc.

Altered Expression of Long Noncoding RNAs in Blood After Ischemic Stroke and Proximity to Putative Stroke Risk Loci.

BACKGROUND AND PURPOSE: Although peripheral blood mRNA and micro-RNA change after ischemic stroke, any role for long noncoding RNA (lncRNA), which comprise most of the genome and have been implicated in various diseases, is unknown. Thus, we hypothesized that lncRNA expression also changes after stroke. METHODS: lncRNA expression was assessed in 266 whole-blood RNA samples drawn once per individual from patients with ischemic stroke and matched with vascular risk factor controls. Differential lncRNA expression was assessed by ANCOVA (P<0.005; fold change>|1.2|), principal components analysis, and hierarchical clustering on a derivation set (n=176) and confirmed on a validation set (n=90). Poststroke temporal lncRNA expression changes were assessed using ANCOVA with confounding factor correction (P<0.005; partial correlation with time since event >|0.4|). Because sexual dimorphism exists in stroke, analyses were performed for each sex separately. RESULTS: A total of 299 lncRNAs were differentially expressed between stroke and control males, whereas 97 lncRNAs were differentially expressed between stroke and control females. Significant changes of lncRNA expression with time after stroke were detected for 49 lncRNAs in men and 31 lncRNAs in women. Some differentially expressed lncRNAs mapped close to genomic locations of previously identified putative stroke-risk genes, including lipoprotein, lipoprotein(a)-like 2, ABO (transferase A, α1-3-N-acetylgalactosaminyltransferase; transferase B, α1-3-galactosyltransferase) blood group, prostaglandin 12 synthase, and α-adducins. CONCLUSIONS: This study provides evidence of altered and sexually dimorphic lncRNA expression in peripheral blood of patients with stroke compared with that of controls and suggests that lncRNAs have potential for stroke biomarker development. Some regulated lncRNA could regulate some previously identified putative stroke-risk genes.

Alternative Splicing of Putative Stroke/Vascular Risk Factor Genes Expressed in Blood Following Ischemic Stroke Is Sexually Dimorphic and Cause-Specific.

Genome-wide association studies have identified putative ischemic stroke risk genes, yet, their expression after stroke is unexplored in spite of growing interest in elucidating their specific role and identifying candidate genes for stroke treatment. Thus, we took an exploratory approach to investigate sexual dimorphism, alternative splicing, and etiology in putative risk gene expression in blood following cardioembolic, atherosclerotic large vessel disease and small vessel disease/lacunar causes of ischemic stroke in each sex compared to controls. Whole transcriptome arrays assessed 71 putative stroke/vascular risk factor genes for blood RNA expression at gene-, exon-, and alternative splicing-levels. Male (n = 122) and female (n = 123) stroke and control volunteers from three university medical centers were matched for race, age, vascular risk factors, and blood draw time since stroke onset. Exclusion criteria included: previous stroke, drug abuse, subarachnoid or intracerebral hemorrhage, hemorrhagic transformation, infection, dialysis, cancer, hematological abnormalities, thrombolytics, anticoagulants or immunosuppressants. Significant differential gene expression (fold change > |1.2|, p < 0.05, partial correlation > |0.4|) and alternative splicing (false discovery rate p < 0.3) were assessed. At gene level, few were differentially expressed: ALDH2, ALOX5AP, F13A1, and IMPA2 (males, all stroke); ITGB3 (females, cardioembolic); ADD1 (males, atherosclerotic); F13A1, IMPA2 (males, lacunar); and WNK1 (females, lacunar). GP1BA and ITGA2B were alternatively spliced in both sexes (all patients vs. controls). Six genes in males, five in females, were alternatively spliced in all stroke compared to controls. Alternative splicing and exon-level analyses associated many genes with specific etiology in either sex. Of 71 genes, 70 had differential exon-level expression in stroke patients compared to control subjects. Among stroke patients, 24 genes represented by differentially expressed exons were male-specific, six were common between sexes, and two were female-specific. In lacunar stroke, expression of 19 differentially expressed exons representing six genes (ADD1, NINJ2, PCSK9, PEMT, SMARCA4, WNK1) decreased in males and increased in females. Results demonstrate alternative splicing and sexually dimorphic expression of most putative risk genes in stroke patients' blood. Since expression was also often cause-specific, sex, and etiology are factors to consider in stroke treatment trials and genetic association studies as society trends toward more personalized medicine.

Analysis of the wheat endosperm transcriptome.

Among the cereals, wheat is the most widely grown geographically and is part of the staple diet in much of the world. Understanding how the cereal endosperm develops and functions will help generate better tools to manipulate grain qualities important to end-users. We used a genomics approach to identify and characterize genes that are expressed in the wheat endosperm. We analyzed the 17,949 publicly available wheat endosperm EST sequences to identify genes involved in the biological processes that occur within this tissue. Clustering and assembly of the ESTs resulted in the identification of 6,187 tentative unique genes, 2,358 of which formed contigs and 3,829 remained as singletons. A BLAST similarity search against the NCBI non-redundant sequence database revealed abundant messages for storage proteins, putative defense proteins, and proteins involved in starch and sucrose metabolism. The level of abundance of the putatively identified genes reflects the physiology of the developing endosperm. Half of the identified genes have unknown functions. Approximately 61% of the endosperm ESTs has been tentatively mapped in the hexaploid wheat genome. Using microarrays for global RNA profiling, we identified endosperm genes that are specifically up regulated in the developing grain.

Evidence for differential alternative splicing in blood of young boys with autism spectrum disorders.

BACKGROUND: Since RNA expression differences have been reported in autism spectrum disorder (ASD) for blood and brain, and differential alternative splicing (DAS) has been reported in ASD brains, we determined if there was DAS in blood mRNA of ASD subjects compared to typically developing (TD) controls, as well as in ASD subgroups related to cerebral volume. METHODS: RNA from blood was processed on whole genome exon arrays for 2-4-year-old ASD and TD boys. An ANCOVA with age and batch as covariates was used to predict DAS for ALL ASD (n=30), ASD with normal total cerebral volumes (NTCV), and ASD with large total cerebral volumes (LTCV) compared to TD controls (n=20). RESULTS: A total of 53 genes were predicted to have DAS for ALL ASD versus TD, 169 genes for ASD_NTCV versus TD, 1 gene for ASD_LTCV versus TD, and 27 genes for ASD_LTCV versus ASD_NTCV. These differences were significant at P <0.05 after false discovery rate corrections for multiple comparisons (FDR <5% false positives). A number of the genes predicted to have DAS in ASD are known to regulate DAS (SFPQ, SRPK1, SRSF11, SRSF2IP, FUS, LSM14A). In addition, a number of genes with predicted DAS are involved in pathways implicated in previous ASD studies, such as ROS monocyte/macrophage, Natural Killer Cell, mTOR, and NGF signaling. The only pathways significant after multiple comparison corrections (FDR <0.05) were the Nrf2-mediated reactive oxygen species (ROS) oxidative response (superoxide dismutase 2, catalase, peroxiredoxin 1, PIK3C3, DNAJC17, microsomal glutathione S-transferase 3) and superoxide radical degradation (SOD2, CAT). CONCLUSIONS: These data support differences in alternative splicing of mRNA in blood of ASD subjects compared to TD controls that differ related to head size. The findings are preliminary, need to be replicated in independent cohorts, and predicted alternative splicing differences need to be confirmed using direct analytical methods.

GABA- and acetylcholine-related gene expression in blood correlate with tic severity and microarray evidence for alternative splicing in Tourette syndrome: a pilot study.

Tourette syndrome (TS) is a complex childhood neurodevelopmental disorder characterized by motor and vocal tics. Recently, altered numbers of GABAergic-parvalbumin (PV) and cholinergic interneurons were observed in the basal ganglia of individuals with TS. Thus, we postulated that gamma-amino butyric acid (GABA)- and acetylcholine (ACh)-related genes might be associated with the pathophysiology of TS. Total RNA isolated from whole blood of 26 un-medicated TS subjects and 23 healthy controls (HC) was processed on Affymetrix Human Exon 1.0 ST arrays. Data were analyzed to identify genes whose expression correlated with tic severity in TS, and to identify genes differentially spliced in TS compared to HC subjects. Many genes (3627) correlated with tic severity in TS (p < 0.05) among which GABA- (p = 2.1 × 10⁻³) and ACh- (p = 4.25 × 10⁻8) related genes were significantly over-represented. Moreover, several GABA and ACh-related genes were predicted to be alternatively spliced in TS compared to HC including GABA receptors GABRA4 and GABRG1, the nicotinic ACh receptor CHRNA4 and cholinergic differentiation factor (CDF). This pilot study suggests that at least some of these GABA- and ACh-related genes observed in blood that correlate with tics or are alternatively spliced are involved in the pathophysiology of TS and tics.

Brain and blood microRNA expression profiling of ischemic stroke, intracerebral hemorrhage, and kainate seizures.

MicroRNAs (miRNAs) regulate gene expression and have a critical role in many biologic and pathologic processes. We hypothesized that miRNA expression profiles in injured brain (hippocampus) would show common as well as unique profiles when compared with those of blood. Adult, untouched, control rats were compared with rats with sham surgeries, ischemic strokes, brain hemorrhage (lysed blood, fresh blood, or thrombin), and kainate-induced seizures. Brain and whole-blood miRNA expression profiles were assessed 24 h later using TaqMan rodent miRNA arrays. MicroRNA response profiles were different for each condition. Many miRNAs changed more than 1.5-fold in brain and blood after each experimental manipulation, and several miRNAs were upregulated or downregulated in both brain and blood after a given injury. A few miRNAs (e.g., miR-298, miR-155, and miR-362-3p) were upregulated or downregulated more than twofold in both brain and blood after several different injuries. The results show the possible use of blood miRNAs as biomarkers for brain injury; that selected blood miRNAs may correlate with miRNA changes in the brain; and that many of the mRNAs, previously shown to be regulated in brain and blood after brain injury, are likely accounted for by changes in miRNA expression.

Transcriptomic analysis of starch biosynthesis in the developing grain of hexaploid wheat.

The expression of genes involved in starch synthesis in wheat was analyzed together with the accumulation profiles of soluble sugars, starch, protein, and starch granule distribution in developing caryopses obtained from the same biological materials used for profiling of gene expression using DNA microarrays. Multiple expression patterns were detected for the different starch biosynthetic gene isoforms, suggesting their relative importance through caryopsis development. Members of the ADP-glucose pyrophosphorylase, starch synthase, starch branching enzyme, and sucrose synthase gene families showed different expression profiles; expression of some members of these gene families coincided with a period of high accumulation of starch while others did not. A biphasic pattern was observed in the rates of starch and protein accumulation which paralleled changes in global gene expression. Metabolic and regulatory genes that show a pattern of expression similar to starch accumulation and granule size distribution were identified, suggesting their coinvolvement in these biological processes.

Identification and validation of suitable endogenous reference genes for gene expression studies in human peripheral blood.

BACKGROUND: Gene expression studies require appropriate normalization methods. One such method uses stably expressed reference genes. Since suitable reference genes appear to be unique for each tissue, we have identified an optimal set of the most stably expressed genes in human blood that can be used for normalization. METHODS: Whole-genome Affymetrix Human 2.0 Plus arrays were examined from 526 samples of males and females ages 2 to 78, including control subjects and patients with Tourette syndrome, stroke, migraine, muscular dystrophy, and autism. The top 100 most stably expressed genes with a broad range of expression levels were identified. To validate the best candidate genes, we performed quantitative RT-PCR on a subset of 10 genes (TRAP1, DECR1, FPGS, FARP1, MAPRE2, PEX16, GINS2, CRY2, CSNK1G2 and A4GALT), 4 commonly employed reference genes (GAPDH, ACTB, B2M and HMBS) and PPIB, previously reported to be stably expressed in blood. Expression stability and ranking analysis were performed using GeNorm and NormFinder algorithms. RESULTS: Reference genes were ranked based on their expression stability and the minimum number of genes needed for nomalization as calculated using GeNorm showed that the fewest, most stably expressed genes needed for acurate normalization in RNA expression studies of human whole blood is a combination of TRAP1, FPGS, DECR1 and PPIB. We confirmed the ranking of the best candidate control genes by using an alternative algorithm (NormFinder). CONCLUSION: The reference genes identified in this study are stably expressed in whole blood of humans of both genders with multiple disease conditions and ages 2 to 78. Importantly, they also have different functions within cells and thus should be expressed independently of each other. These genes should be useful as normalization genes for microarray and RT-PCR whole blood studies of human physiology, metabolism and disease.

Transcriptional profiling of wheat caryopsis development using cDNA microarrays.

The expression of 7,835 genes in developing wheat caryopses was analyzed using cDNA arrays. Using a mixed model analysis of variance (ANOVA) method, 29% (2,237) of the genes on the array were identified to be differentially expressed at the 6 different time-points examined, which covers the developmental stages from coenocytic endosperm to physiological maturity. Comparison of genes differentially expressed between two time-points revealed a dynamic transcript accumulation profile with major re-programming events that occur at 3-7, 7-14 and 21-28 DPA. A k-means clustering algorithm grouped the differentially expressed genes into 10 clusters, revealing co-expression of genes involved in the same pathway such as carbohydrate and protein synthesis or preparation for desiccation. Functional annotation of genes that show peak expression at specific time-points correlated with the developmental events associated with the respective stages. Results provide information on the temporal expression during caryopsis development for a significant number of differentially expressed genes with unknown function.