cerebroViz: an R package for anatomical visualization of spatiotemporal brain data

Multifaceted impact of specialized neuropeptide-intensive neurons on the selective vulnerability in Alzheimer's disease

INTRODUCTION

Widespread disruption of neuropeptide (NP) networks in Alzheimer's disease (AD) and disproportionate absence of neurons expressing high NP-producing, coined as HNP neurons, have been reported for the entorhinal cortex (EC) of AD brains. Hypothesizing that functional features of HNP neurons are involved in the early pathogenesis of AD, we aim to understand the molecular mechanisms underlying these observations.

METHODS

Multiscale and spatiotemporal transcriptomic analysis was used to investigate AD-afflicted and healthy brains. Our focus encompassed NP expression dynamics in AD, AD-associated NPs (ADNPs) trajectories with aging, and the neuroanatomical distribution of HNP neuron.

RESULTS

Findings include that 1) HNP neurons exhibited heightened metabolic needs and an upregulation of gene expressions linked to protein misfolding; 2) dysfunctions of ADNP production occurred in aging and mild cognitive decline; 3) HNP neurons co-expressing ADNPs were preferentially distributed in brain regions susceptible to AD.

DISCUSSION

We identified potential mechanisms that contribute to the selective vulnerability of HNP neurons to AD. Our results indicate that the functions of HNP neurons predispose them to oxidative stress and protein misfolding, potentially serving as inception sites for misfolded proteins in AD.

Genetically Regulated Gene Expression in the Brain Associated With Chronic Pain: Relationships With Clinical Traits and Potential for Drug Repurposing

BACKGROUND

Chronic pain is a common, poorly understood condition. Genetic studies including genome-wide association studies have identified many relevant variants, which have yet to be translated into full understanding of chronic pain. Transcriptome-wide association studies using transcriptomic imputation methods such as S-PrediXcan can help bridge this genotype-phenotype gap.

METHODS

We carried out transcriptomic imputation using S-PrediXcan to identify genetically regulated gene expression associated with multisite chronic pain in 13 brain tissues and whole blood. Then, we imputed genetically regulated gene expression for over 31,000 Mount Sinai BioMe participants and performed a phenome-wide association study to investigate clinical relationships in chronic pain-associated gene expression changes.

RESULTS

We identified 95 experiment-wide significant gene-tissue associations (p < 7.97 × 10-7), including 36 unique genes and an additional 134 gene-tissue associations reaching within-tissue significance, including 53 additional unique genes. Of the 89 unique genes in total, 59 were novel for multisite chronic pain and 18 are established drug targets. Chronic pain genetically regulated gene expression for 10 unique genes was significantly associated with cardiac dysrhythmia, metabolic syndrome, disc disorders/dorsopathies, joint/ligament sprain, anemias, and neurologic disorder phecodes. Phenome-wide association study analyses adjusting for mean pain score showed that associations were not driven by mean pain score.

CONCLUSIONS

We carried out the largest transcriptomic imputation study of any chronic pain trait to date. Results highlight potential causal genes in chronic pain development and tissue and direction of effect. Several gene results were also drug targets. Phenome-wide association study results showed significant associations for phecodes including cardiac dysrhythmia and metabolic syndrome, thereby indicating potential shared mechanisms.

Pleiotropy with sex-specific traits reveals genetic aspects of sex differences in Parkinson's disease

Parkinson's disease is an age-related neurodegenerative disorder with a higher incidence in males than females. The causes for this sex difference are unknown. Genome-wide association studies (GWAS) have identified 90 Parkinson's disease risk loci, but the genetic studies have not found sex-specific differences in allele frequency on autosomal chromosomes or sex chromosomes. Genetic variants, however, could exert sex-specific effects on gene function and regulation of gene expression. To identify genetic loci that might have sex-specific effects, we studied pleiotropy between Parkinson's disease and sex-specific traits. Summary statistics from GWASs were acquired from large-scale consortia for Parkinson's disease (n cases = 13 708; n controls = 95 282), age at menarche (n = 368 888 females) and age at menopause (n = 69 360 females). We applied the conditional/conjunctional false discovery rate (FDR) method to identify shared loci between Parkinson's disease and these sex-specific traits. Next, we investigated sex-specific gene expression differences in the superior frontal cortex of both neuropathologically healthy individuals and Parkinson's disease patients (n cases = 61; n controls = 23). To provide biological insights to the genetic pleiotropy, we performed sex-specific expression quantitative trait locus (eQTL) analysis and sex-specific age-related differential expression analysis for genes mapped to Parkinson's disease risk loci. Through conditional/conjunctional FDR analysis we found 11 loci shared between Parkinson's disease and the sex-specific traits age at menarche and age at menopause. Gene-set and pathway analysis of the genes mapped to these loci highlighted the importance of the immune response in determining an increased disease incidence in the male population. Moreover, we highlighted a total of nine genes whose expression or age-related expression in the human brain is influenced by genetic variants in a sex-specific manner. With these analyses we demonstrated that the lack of clear sex-specific differences in allele frequencies for Parkinson's disease loci does not exclude a genetic contribution to differences in disease incidence. Moreover, further studies are needed to elucidate the role that the candidate genes identified here could have in determining a higher incidence of Parkinson's disease in the male population.

Genome-wide analysis of anorexia nervosa and major psychiatric disorders and related traits reveals genetic overlap and identifies novel risk loci for anorexia nervosa

Anorexia nervosa (AN) is a heritable eating disorder (50-60%) with an array of commonly comorbid psychiatric disorders and related traits. Although significant genetic correlations between AN and psychiatric disorders and related traits have been reported, their shared genetic architecture is largely understudied. We investigated the shared genetic architecture of AN and schizophrenia (SCZ), bipolar disorder (BIP), major depression (MD), mood instability (Mood), neuroticism (NEUR), and intelligence (INT). We applied the conditional false discovery rate (FDR) method to identify novel risk loci for AN, and conjunctional FDR to identify loci shared between AN and related phenotypes, to summarize statistics from relevant genome-wide association studies (GWAS). Individual GWAS samples varied from 72,517 to 420,879 participants. Using conditional FDR we identified 58 novel AN loci. Furthermore, we identified 38 unique loci shared between AN and major psychiatric disorders (SCZ, BIP, and MD) and 45 between AN and psychological traits (Mood, NEUR, and INT). In line with genetic correlations, the majority of shared loci showed concordant effect directions. Functional analyses revealed that the shared loci are involved in 65 unique pathways, several of which overlapped across analyses, including the "signal by MST1" pathway involved in Hippo signaling. In conclusion, we demonstrated genetic overlap between AN and major psychiatric disorders and related traits, and identified novel risk loci for AN by leveraging this overlap. Our results indicate that some shared characteristics between AN and related disorders and traits may have genetic underpinnings.

The contributions of mitochondrial and nuclear mitochondrial genetic variation to neuroticism

Neuroticism is a heritable trait composed of separate facets, each conferring different levels of protection or risk, to health. By examining mitochondrial DNA in 269,506 individuals, we show mitochondrial haplogroups explain 0.07-0.01% of variance in neuroticism and identify five haplogroup and 15 mitochondria-marker associations across a general factor of neuroticism, and two special factors of anxiety/tension, and worry/vulnerability with effect sizes of the same magnitude as autosomal variants. Within-haplogroup genome-wide association studies identified H-haplogroup-specific autosomal effects explaining 1.4% variance of worry/vulnerability. These H-haplogroup-specific autosomal effects show a pleiotropic relationship with cognitive, physical and mental health that differs from that found when assessing autosomal effects across haplogroups. We identify interactions between chromosome 9 regions and mitochondrial haplogroups at P < 5 × 10^-8, revealing associations between general neuroticism and anxiety/tension with brain-specific gene co-expression networks. These results indicate that the mitochondrial genome contributes toward neuroticism and the autosomal links between neuroticism and health.

Multiomic neuropathology improves diagnostic accuracy in pediatric neuro-oncology

The large diversity of central nervous system (CNS) tumor types in children and adolescents results in disparate patient outcomes and renders accurate diagnosis challenging. In this study, we prospectively integrated DNA methylation profiling and targeted gene panel sequencing with blinded neuropathological reference diagnostics for a population-based cohort of more than 1,200 newly diagnosed pediatric patients with CNS tumors, to assess their utility in routine neuropathology. We show that the multi-omic integration increased diagnostic accuracy in a substantial proportion of patients through annotation to a refining DNA methylation class (50%), detection of diagnostic or therapeutically relevant genetic alterations (47%) or identification of cancer predisposition syndromes (10%). Discrepant results by neuropathological WHO-based and DNA methylation-based classification (30%) were enriched in histological high-grade gliomas, implicating relevance for current clinical patient management in 5% of all patients. Follow-up (median 2.5 years) suggests improved survival for patients with histological high-grade gliomas displaying lower-grade molecular profiles. These results provide preliminary evidence of the utility of integrating multi-omics in neuropathology for pediatric neuro-oncology.

Population-level variation in enhancer expression identifies disease mechanisms in the human brain

Identification of risk variants for neuropsychiatric diseases within enhancers underscores the importance of understanding population-level variation in enhancer function in the human brain. Besides regulating tissue-specific and cell-type-specific transcription of target genes, enhancers themselves can be transcribed. By jointly analyzing large-scale cell-type-specific transcriptome and regulome data, we cataloged 30,795 neuronal and 23,265 non-neuronal candidate transcribed enhancers. Examination of the transcriptome in 1,382 brain samples identified robust expression of transcribed enhancers. We explored gene-enhancer coordination and found that enhancer-linked genes are strongly implicated in neuropsychiatric disease. We identified expression quantitative trait loci (eQTLs) for both genes and enhancers and found that enhancer eQTLs mediate a substantial fraction of neuropsychiatric trait heritability. Inclusion of enhancer eQTLs in transcriptome-wide association studies enhanced functional interpretation of disease loci. Overall, our study characterizes the gene-enhancer regulome and genetic mechanisms in the human cortex in both healthy and diseased states.

Distinct tau neuropathology and cellular profiles of an APOE3 Christchurch homozygote protected against autosomal dominant Alzheimer's dementia

We describe in vivo follow-up PET imaging and postmortem findings from an autosomal dominant Alzheimer's disease (ADAD) PSEN1 E280A carrier who was also homozygous for the APOE3 Christchurch (APOE3ch) variant and was protected against Alzheimer's symptoms for almost three decades beyond the expected age of onset. We identified a distinct anatomical pattern of tau pathology with atypical accumulation in vivo and unusual postmortem regional distribution characterized by sparing in the frontal cortex and severe pathology in the occipital cortex. The frontal cortex and the hippocampus, less affected than the occipital cortex by tau pathology, contained Related Orphan Receptor B (RORB) positive neurons, homeostatic astrocytes and higher APOE expression. The occipital cortex, the only cortical region showing cerebral amyloid angiopathy (CAA), exhibited a distinctive chronic inflammatory microglial profile and lower APOE expression. Thus, the Christchurch variant may impact the distribution of tau pathology, modulate age at onset, severity, progression, and clinical presentation of ADAD, suggesting possible therapeutic strategies.

Molecular Modelling Hurdle in the Next-Generation Sequencing Era

There are challenges in the genetic diagnosis of rare diseases, and pursuing an optimal strategy to identify the cause of the disease is one of the main objectives of any clinical genomics unit. A range of techniques are currently used to characterize the genomic variability within the human genome to detect causative variants of specific disorders. With the introduction of next-generation sequencing (NGS) in the clinical setting, geneticists can study single-nucleotide variants (SNVs) throughout the entire exome/genome. In turn, the number of variants to be evaluated per patient has increased significantly, and more information has to be processed and analyzed to determine a proper diagnosis. Roughly 50% of patients with a Mendelian genetic disorder are diagnosed using NGS, but a fair number of patients still suffer a diagnostic odyssey. Due to the inherent diversity of the human population, as more exomes or genomes are sequenced, variants of uncertain significance (VUSs) will increase exponentially. Thus, assigning relevance to a VUS (non-synonymous as well as synonymous) in an undiagnosed patient becomes crucial to assess the proper diagnosis. Multiple algorithms have been used to predict how a specific mutation might affect the protein’s function, but they are far from accurate enough to be conclusive. In this work, we highlight the difficulties of genomic variability determined by NGS that have arisen in diagnosing rare genetic diseases, and how molecular modelling has to be a key component to elucidate the relevance of a specific mutation in the protein’s loss of function or malfunction. We suggest that the creation of a multi-omics data model should improve the classification of pathogenicity for a significant amount of the detected genomic variability. Moreover, we argue how it should be incorporated systematically in the process of variant evaluation to be useful in the clinical setting and the diagnostic pipeline.

Association of APOE Genotypes and Chronic Traumatic Encephalopathy

Importance

Repetitive head impact (RHI) exposure is the chief risk factor for chronic traumatic encephalopathy (CTE). However, the occurrence and severity of CTE varies widely among those with similar RHI exposure. Limited evidence suggests that the APOEε4 allele may confer risk for CTE, but previous studies were small with limited scope.

Objective

To test the association between APOE genotype and CTE neuropathology and related endophenotypes.

Design, Setting, and Participants

This cross-sectional genetic association study analyzed brain donors from February 2008 to August 2019 from the Veterans Affairs-Boston University-Concussion Legacy Foundation Brain Bank. All donors had exposure to RHI from contact sports or military service. All eligible donors were included. Analysis took place between June 2020 and April 2022.

Exposures

One or more APOEε4 or APOEε2 alleles.

Main Outcomes and Measures

CTE neuropathological status, CTE stage (0-IV), semiquantitative phosphorylated tau (p-tau) burden in 11 brain regions (0-3), quantitative p-tau burden in the dorsolateral frontal lobe (log-transformed AT8+ pixel count per mm2), and dementia.

Results

Of 364 consecutive brain donors (100% male; 53 [14.6%] self-identified as Black and 311 [85.4%] as White; median [IQR] age, 65 [47-77] years) 20 years or older, there were 294 individuals with CTE and 70 controls. Among donors older than 65 years, APOEε4 status was significantly associated with CTE stage (odds ratio [OR], 2.34 [95% CI, 1.30-4.20]; false discovery rate [FDR]-corrected P = .01) and quantitative p-tau burden in the dorsolateral frontal lobe (β, 1.39 [95% CI, 0.83-1.94]; FDR-corrected P = 2.37 × 10-5). There was a nonsignificant association between APOEε4 status and dementia (OR, 2.64 [95% CI, 1.06-6.61]; FDR-corrected P = .08). Across 11 brain regions, significant associations were observed for semiquantitative p-tau burden in the frontal and parietal cortices, amygdala, and entorhinal cortex (OR range, 2.45-3.26). Among football players, the APOEε4 association size for CTE stage was similar to playing more than 7 years of football. Associations were significantly larger in the older half of the sample. There was no significant association for CTE status. Association sizes were similar when donors with an Alzheimer disease neuropathological diagnosis were excluded and were reduced but remained significant after adjusting for neuritic and diffuse amyloid plaques. No associations were observed for APOEε2 status. Models were adjusted for age at death and race.

Conclusions and Relevance

APOEε4 may confer increased risk for CTE-related neuropathological and clinical outcomes among older individuals with RHI exposure. Further work is required to validate these findings in an independent sample.

Quantitative prediction of variant effects on alternative splicing in MAPT using endogenous pre-messenger RNA structure probing

Splicing is highly regulated and is modulated by numerous factors. Quantitative predictions for how a mutation will affect precursor mRNA (pre-mRNA) structure and downstream function are particularly challenging. Here, we use a novel chemical probing strategy to visualize endogenous precursor and mature MAPT mRNA structures in cells. We used these data to estimate Boltzmann suboptimal structural ensembles, which were then analyzed to predict consequences of mutations on pre-mRNA structure. Further analysis of recent cryo-EM structures of the spliceosome at different stages of the splicing cycle revealed that the footprint of the B^act complex with pre-mRNA best predicted alternative splicing outcomes for exon 10 inclusion of the alternatively spliced MAPT gene, achieving 74% accuracy. We further developed a β-regression weighting framework that incorporates splice site strength, RNA structure, and exonic/intronic splicing regulatory elements capable of predicting, with 90% accuracy, the effects of 47 known and 6 newly discovered mutations on inclusion of exon 10 of MAPT. This combined experimental and computational framework represents a path forward for accurate prediction of splicing-related disease-causing variants.

A human iPSC-astroglia neurodevelopmental model reveals divergent transcriptomic patterns in schizophrenia

While neurodevelopmental abnormalities have been associated with schizophrenia (SCZ), the role of astroglia in disease pathophysiology remains poorly understood. In the present study, we used a human induced pluripotent stem cell (iPSC)-derived astrocyte model to investigate the temporal patterns of astroglia differentiation during developmental stages critical for SCZ using RNA sequencing. The model generated astrocyte-specific gene expression patterns during differentiation that corresponded well to astroglia-specific expression signatures of in vivo cortical fetal development. Using this model we identified SCZ-specific expression dynamics, and found that SCZ-associated differentially expressed genes were significantly enriched in the medial prefrontal cortex, striatum, and temporal lobe, targeting VWA5A and ADAMTS19. In addition, SCZ astrocytes displayed alterations in calcium signaling, and significantly decreased glutamate uptake and metalloproteinase activity relative to controls. These results implicate novel transcriptional dynamics in astrocyte differentiation in SCZ together with functional changes that are potentially important biological components of SCZ pathology.

Characterising the shared genetic determinants of bipolar disorder, schizophrenia and risk-taking

Increased risk-taking is a central component of bipolar disorder (BIP) and is implicated in schizophrenia (SCZ). Risky behaviours, including smoking and alcohol use, are overrepresented in both disorders and associated with poor health outcomes. Positive genetic correlations are reported but an improved understanding of the shared genetic architecture between risk phenotypes and psychiatric disorders may provide insights into underlying neurobiological mechanisms. We aimed to characterise the genetic overlap between risk phenotypes and SCZ, and BIP by estimating the total number of shared variants using the bivariate causal mixture model and identifying shared genomic loci using the conjunctional false discovery rate method. Summary statistics from genome wide association studies of SCZ, BIP, risk-taking and risky behaviours were acquired (n = 82,315-466,751). Genomic loci were functionally annotated using FUMA. Of 8.6-8.7 K variants predicted to influence BIP, 6.6 K and 7.4 K were predicted to influence risk-taking and risky behaviours, respectively. Similarly, of 10.2-10.3 K variants influencing SCZ, 9.6 and 8.8 K were predicted to influence risk-taking and risky behaviours, respectively. We identified 192 loci jointly associated with SCZ and risk phenotypes and 206 associated with BIP and risk phenotypes, of which 68 were common to both risk-taking and risky behaviours and 124 were novel to SCZ or BIP. Functional annotation implicated differential expression in multiple cortical and sub-cortical regions. In conclusion, we report extensive polygenic overlap between risk phenotypes and BIP and SCZ, identify specific loci contributing to this shared risk and highlight biologically plausible mechanisms that may underlie risk-taking in severe psychiatric disorders.

Lithium increases mitochondrial respiration in iPSC-derived neural precursor cells from lithium responders

Lithium (Li), valproate (VPA) and lamotrigine (LTG) are commonly used to treat bipolar disorder (BD). While their clinical efficacy is well established, the mechanisms of action at the molecular level are still incompletely understood. Here we investigated the molecular effects of Li, LTG and VPA treatment in induced pluripotent stem cell (iPSC)-derived neural precursor cells (NPCs) generated from 3 healthy controls (CTRL), 3 affective disorder Li responsive patients (Li-R) and 3 Li non-treated patients (Li-N) after 6 h and 1 week of exposure. Differential expression (DE) analysis after 6 h of treatment revealed a transcriptional signature that was associated with all three drugs and most significantly enriched for ribosome and oxidative phosphorylation (OXPHOS) pathways. In addition to the shared DE genes, we found that Li exposure was associated with 554 genes uniquely regulated in Li-R NPCs and enriched for spliceosome, OXPHOS and thermogenesis pathways. In-depth analysis of the treatment-associated transcripts uncovered a significant decrease in intron retention rate, suggesting that the beneficial influence of these drugs might partly be related to splicing. We examined the mitochondrial respiratory function of the NPCs by exploring the drugs' effects on oxygen consumption rate (OCR) and glycolytic rate (ECAR). Li improved OCR levels only in Li-R NPCs by enhancing maximal respiration and reserve capacity, while VPA enhanced maximal respiration and reserve capacity in Li-N NPCs. Overall, our findings further support the involvement of mitochondrial functions in the molecular mechanisms of mood stabilizers and suggest novel mechanisms related to the spliceosome, which warrant further investigation.

Gene expression imputation and cell-type deconvolution in human brain with spatiotemporal precision and its implications for brain-related disorders

As the most complex organ of the human body, the brain is composed of diverse regions, each consisting of distinct cell types and their respective cellular interactions. Human brain development involves a finely tuned cascade of interactive events. These include spatiotemporal gene expression changes and dynamic alterations in cell-type composition. However, our understanding of this process is still largely incomplete owing to the difficulty of brain spatiotemporal transcriptome collection. In this study, we developed a tensor-based approach to impute gene expression on a transcriptome-wide level. After rigorous computational benchmarking, we applied our approach to infer missing data points in the widely used BrainSpan resource and completed the entire grid of spatiotemporal transcriptomics. Next, we conducted deconvolutional analyses to comprehensively characterize major cell-type dynamics across the entire BrainSpan resource to estimate the cellular temporal changes and distinct neocortical areas across development. Moreover, integration of these results with GWAS summary statistics for 13 brain-associated traits revealed multiple novel trait–cell-type associations and trait-spatiotemporal relationships. In summary, our imputed BrainSpan transcriptomic data provide a valuable resource for the research community and our findings help further studies of the transcriptional and cellular dynamics of the human brain and related diseases.

Common schizophrenia risk variants are enriched in open chromatin regions of human glutamatergic neurons

The chromatin landscape of human brain cells encompasses key information to understanding brain function. Here we use ATAC-seq to profile the chromatin structure in four distinct populations of cells (glutamatergic neurons, GABAergic neurons, oligodendrocytes, and microglia/astrocytes) from three different brain regions (anterior cingulate cortex, dorsolateral prefrontal cortex, and primary visual cortex) in human postmortem brain samples. We find that chromatin accessibility varies greatly by cell type and, more moderately, by brain region, with glutamatergic neurons showing the largest regional variability. Transcription factor footprinting implicates cell-specific transcriptional regulators and infers cell-specific regulation of protein-coding genes, long intergenic noncoding RNAs and microRNAs. In vivo transgenic mouse experiments validate the cell type specificity of several of these human-derived regulatory sequences. We find that open chromatin regions in glutamatergic neurons are enriched for neuropsychiatric risk variants, particularly those associated with schizophrenia. Integration of cell-specific chromatin data with a bulk tissue study of schizophrenia brains increases statistical power and confirms that glutamatergic neurons are most affected. These findings illustrate the utility of studying the cell-type-specific epigenome in complex tissues like the human brain, and the potential of such approaches to better understand the genetic basis of human brain function.

Primate-restricted KRAB zinc finger proteins and target retrotransposons control gene expression in human neurons

In the first days of embryogenesis, transposable element-embedded regulatory sequences (TEeRS) are silenced by Kruppel-associated box (KRAB) zinc finger proteins (KZFPs). Many TEeRS are subsequently co-opted in transcription networks, but how KZFPs influence this process is largely unknown. We identify ZNF417 and ZNF587 as primate-specific KZFPs repressing HERVK (human endogenous retrovirus K) and SVA (SINE-VNTR-Alu) integrants in human embryonic stem cells (ESCs). Expressed in specific regions of the human developing and adult brain, ZNF417/587 keep controlling TEeRS in ESC-derived neurons and brain organoids, secondarily influencing the differentiation and neurotransmission profile of neurons and preventing the induction of neurotoxic retroviral proteins and an interferon-like response. Thus, evolutionarily recent KZFPs and their TE targets partner up to influence human neuronal differentiation and physiology.

mTADA is a framework for identifying risk genes from de novo mutations in multiple traits

Joint analysis of multiple traits can result in the identification of associations not found through the analysis of each trait in isolation. Studies of neuropsychiatric disorders and congenital heart disease (CHD) which use de novo mutations (DNMs) from parent-offspring trios have reported multiple putatively causal genes. However, a joint analysis method designed to integrate DNMs from multiple studies has yet to be implemented. We here introduce multiple-trait TADA (mTADA) which jointly analyzes two traits using DNMs from non-overlapping family samples. We first demonstrate that mTADA is able to leverage genetic overlaps to increase the statistical power of risk-gene identification. We then apply mTADA to large datasets of >13,000 trios for five neuropsychiatric disorders and CHD. We report additional risk genes for schizophrenia, epileptic encephalopathies and CHD. We outline some shared and specific biological information of intellectual disability and CHD by conducting systems biology analyses of genes prioritized by mTADA.

A tissue level atlas of the healthy human virome

BACKGROUND

Human-resident microbes can influence both health and disease. Investigating the microbiome using next-generation sequencing technology has revealed examples of mutualism and conflict between microbes and humans. Comparing to bacteria, the viral component of the microbiome (i.e., the "virome") is understudied. Somatic tissues of healthy individuals are usually inaccessible for the virome sampling; therefore, there is limited understanding of the presence and distribution of viruses in tissues in healthy individuals and how virus infection associates with human gene expression and perturbs immunological homeostasis.

RESULTS

To characterize the human virome in a tissue-specific manner, here we performed meta-transcriptomic analysis using the RNA-sequencing dataset from the Genotype-Tissue Expression (GTEx) Project. We analyzed the 8991 RNA-sequencing data obtained from 51 somatic tissues from 547 individuals and successfully detected 39 viral species in at least one tissue. We then investigated associations between virus infection and human gene expression and human disease onset. We detected some expected relationships; for instance, hepatitis C virus infection in the liver was strongly associated with interferon-stimulated gene upregulation and pathological findings of chronic hepatitis. The presence of herpes simplex virus type 1 in one subject's brain strongly associated with immune gene expression. While torque teno virus was detected in a broad range of human tissues, it was not associated with interferon responses. Being notable in light of its association with lymphoproliferative disorders, Epstein-Barr virus infection in the spleen and blood was associated with an increase in plasma cells in healthy subjects. Human herpesvirus 7 was often detected in the stomach; intriguingly, it associated with the proportion of human leukocytes in the stomach as well as digestive gene expression. Moreover, virus infections in the local tissues associated with systemic immune responses in circulating blood.

CONCLUSIONS

To our knowledge, this study is the first comprehensive investigation of the human virome in a variety of tissues in healthy individuals through meta-transcriptomic analysis. Further investigation of the associations described here, and application of this analytical pipeline to additional datasets, will be useful to reveal the impact of viral infections on human health.

Modeling Brain Disorders Using Induced Pluripotent Stem Cells

Brain disorders, from neurodegenerative to psychiatric disorders, are among the most challenging conditions to study because of the intricate nature of the human brain and the limitations of existing model systems in recapitulating all these intricacies. However, innovations in stem cell technologies now allow us to reprogram patient somatic cells to induced pluripotent stem cells (iPSCs), which can then be differentiated to disease-relevant neural and glial cells. iPSCs are a valuable tool to model brain disorders, as they can be derived from patients with known symptom histories, genetics, and drug-response profiles. Here, we discuss the premise and validity of the iPSC-based in vitro model system and highlight key findings from the most commonly studied neurodegenerative and psychiatric disorders.

BrainImageR: spatiotemporal gene set analysis referencing the human brain

Motivation

Neuronal analyses such as transcriptomics, epigenetics and genome-wide association studies must be assessed in the context of the human brain to generate biologically meaningful inferences. It is often difficult to access primary human brain tissue; therefore, approximations are made using alternative sources such as peripheral tissues or in vitro-derived neurons. Gene sets from these studies are then assessed for their association with the post-mortem human brain. However, most analyses of post-mortem datasets are achieved by building new computational tools each time in-house, which can cause discrepancies from study to study. The field is in need of a user-friendly tool to examine spatiotemporal expression with respect to the postmortem brain. Such a tool will be of use in the molecular interrogation of neurological and psychiatric disorders, with direct advantages for the disease-modeling and human genetics communities.

Results

We have developed brainImageR, an R package that calculates both the spatial and temporal association of a dataset with post-mortem human brain. BrainImageR identifies anatomical regions enriched for candidate gene set expression. It further predicts the developmental time point of the sample, a task that has become increasingly important in the field of in vitro neuronal modeling. These functionalities of brainImageR enable a quick and efficient characterization of a given dataset across normal human brain development.

Availability and implementation

BrainImageR is released under the Creative Commons CC BY-SA 4.0 license and can be accessed directly at brainimager.salk.edu or the R code can be downloaded through github at https://github.com/saralinker/brainImageR.

Gene expression imputation across multiple brain regions provides insights into schizophrenia risk

Transcriptomic imputation approaches combine eQTL reference panels with large-scale genotype data in order to test associations between disease and gene expression. These genic associations could elucidate signals in complex genome-wide association study (GWAS) loci and may disentangle the role of different tissues in disease development. We used the largest eQTL reference panel for the dorso-lateral prefrontal cortex (DLPFC) to create a set of gene expression predictors and demonstrate their utility. We applied DLPFC and 12 GTEx-brain predictors to 40,299 schizophrenia cases and 65,264 matched controls for a large transcriptomic imputation study of schizophrenia. We identified 413 genic associations across 13 brain regions. Stepwise conditioning identified 67 non-MHC genes, of which 14 did not fall within previous GWAS loci. We identified 36 significantly enriched pathways, including hexosaminidase-A deficiency, and multiple porphyric disorder pathways. We investigated developmental expression patterns among the 67 non-MHC genes and identified specific groups of pre- and postnatal expression.