Single-cell RNA sequencing of liver fine-needle aspirates captures immune diversity in the blood and liver in chronic hepatitis B patients.

BACKGROUND AND AIMS: HBV infection is restricted to the liver, where it drives exhaustion of virus-specific T and B cells and pathogenesis through dysregulation of intrahepatic immunity. Our understanding of liver-specific events related to viral control and liver damage has relied almost solely on animal models, and we lack useable peripheral biomarkers to quantify intrahepatic immune activation beyond cytokine measurement. Our objective was to overcome the practical obstacles of liver sampling using fine-needle aspiration and develop an optimized workflow to comprehensively compare the blood and liver compartments within patients with chronic hepatitis B using single-cell RNA sequencing. APPROACH AND RESULTS: We developed a workflow that enabled multi-site international studies and centralized single-cell RNA sequencing. Blood and liver fine-needle aspirations were collected, and cellular and molecular captures were compared between the Seq-Well S 3 picowell-based and the 10× Chromium reverse-emulsion droplet-based single-cell RNA sequencing technologies. Both technologies captured the cellular diversity of the liver, but Seq-Well S 3 effectively captured neutrophils, which were absent in the 10× dataset. CD8 T cells and neutrophils displayed distinct transcriptional profiles between blood and liver. In addition, liver fine-needle aspirations captured a heterogeneous liver macrophage population. Comparison between untreated patients with chronic hepatitis B and patients treated with nucleoside analogs showed that myeloid cells were highly sensitive to environmental changes while lymphocytes displayed minimal differences. CONCLUSIONS: The ability to electively sample and intensively profile the immune landscape of the liver, and generate high-resolution data, will enable multi-site clinical studies to identify biomarkers for intrahepatic immune activity in HBV and beyond.

Timing and localization of myasthenia gravis-related gene expression.

Myasthenia gravis (MG) is an acquired autoimmune disorder caused by autoantibodies binding acetylcholine receptors (AChR), muscle-specific kinase (MuSK), agrin or low-density lipoprotein receptor-related protein 4 (Lrp4). These autoantibodies inhibit neuromuscular transmission by blocking the function of these proteins and thereby cause fluctuating skeletal muscle weakness. Several reports suggest that these autoantibodies might also affect the central nervous system (CNS) in MG patients. A comprehensive overview of the timing and localization of the expression of MG-related antigens in other organs is currently lacking. To investigate the spatio-temporal expression of MG-related genes outside skeletal muscle, we used in silico tools to assess public expression databases. Acetylcholine esterase, nicotinic AChR α1 subunit, agrin, collagen Q, downstream of kinase-7, Lrp4, MuSK and rapsyn were included as MG-related genes because of their well-known involvement in either congenital or autoimmune MG. We investigated expression of MG-related genes in (1) all human tissues using GTEx data, (2) specific brain regions, (3) neurodevelopmental stages, and (4) cell types using datasets from the Allen Institute for Brain Sciences. MG-related genes show heterogenous spatio-temporal expression patterns in the human body as well as in the CNS. For each of these genes, several (new) tissues, brain areas and cortical cell types with (relatively) high expression were identified suggesting a potential role for these genes outside skeletal muscle. The possible presence of MG-related antigens outside skeletal muscle suggests that autoimmune MG, congenital MG or treatments targeting the same proteins may affect MG-related protein function in other organs.

Cingulate networks associated with gray matter loss in Parkinson's disease show high expression of cholinergic genes in the healthy brain.

Structural covariance networks are able to identify functionally organized brain regions by gray matter volume covariance across a population. We examined the transcriptomic signature of such anatomical networks in the healthy brain using postmortem microarray data from the Allen Human Brain Atlas. A previous study revealed that a posterior cingulate network and anterior cingulate network showed decreased gray matter in brains of Parkinson's disease patients. Therefore, we examined these two anatomical networks to understand the underlying molecular processes that may be involved in Parkinson's disease. Whole brain transcriptomics from the healthy brain revealed upregulation of genes associated with serotonin, GPCR, GABA, glutamate, and RAS-signaling pathways. Our results also suggest involvement of the cholinergic circuit, in which genes NPPA, SOSTDC1, and TYRP1 may play a functional role. Finally, both networks were enriched for genes associated with neuropsychiatric disorders that overlap with Parkinson's disease symptoms. The identified genes and pathways contribute to healthy functions of the posterior and anterior cingulate networks and disruptions to these functions may in turn contribute to the pathological and clinical events observed in Parkinson's disease.

QuantTB - a method to classify mixed Mycobacterium tuberculosis infections within whole genome sequencing data.

BACKGROUND: Mixed infections of Mycobacterium tuberculosis and antibiotic heteroresistance continue to complicate tuberculosis (TB) diagnosis and treatment. Detection of mixed infections has been limited to molecular genotyping techniques, which lack the sensitivity and resolution to accurately estimate the multiplicity of TB infections. In contrast, whole genome sequencing offers sensitive views of the genetic differences between strains of M. tuberculosis within a sample. Although metagenomic tools exist to classify strains in a metagenomic sample, most tools have been developed for more divergent species, and therefore cannot provide the sensitivity required to disentangle strains within closely related bacterial species such as M. tuberculosis. Here we present QuantTB, a method to identify and quantify individual M. tuberculosis strains in whole genome sequencing data. QuantTB uses SNP markers to determine the combination of strains that best explain the allelic variation observed in a sample. QuantTB outputs a list of identified strains, their corresponding relative abundances, and a list of drugs for which resistance-conferring mutations (or heteroresistance) have been predicted within the sample. RESULTS: We show that QuantTB has a high degree of resolution and is capable of differentiating communities differing by less than 25 SNPs and identifying strains down to 1× coverage. Using simulated data, we found QuantTB outperformed other metagenomic strain identification tools at detecting strains and quantifying strain multiplicity. In a real-world scenario, using a dataset of 50 paired clinical isolates from a study of patients with either reinfections or relapses, we found that QuantTB could detect mixed infections and reinfections at rates concordant with a manually curated approach. CONCLUSION: QuantTB can determine infection multiplicity, identify hetero-resistance patterns, enable differentiation between relapse and re-infection, and clarify transmission events across seemingly unrelated patients - even in low-coverage (1×) samples. QuantTB outperforms existing tools and promises to serve as a valuable resource for both clinicians and researchers working with clinical TB samples.

Transcriptional and cell type profiles of cortical brain regions showing ultradian cortisol rhythm dependent responses to emotional face stimulation.

The characteristic endogenous circadian rhythm of plasma glucocorticoid concentrations is made up from an underlying ultradian pulsatile secretory pattern. Recent evidence has indicated that this ultradian cortisol pulsatility is crucial for normal emotional response in man. In this study, we investigate the anatomical transcriptional and cell type signature of brain regions sensitive to a loss of ultradian rhythmicity in the context of emotional processing. We combine human cell type and transcriptomic atlas data of high spatial resolution with functional magnetic resonance imaging (fMRI) data. We show that the loss of cortisol ultradian rhythm alters emotional processing response in cortical brain areas that are characterized by transcriptional and cellular profiles of GABAergic function. We find that two previously identified key components of rapid non-genomic GC signaling - the ANXA1 gene and retrograde endocannabinoid signaling - show most significant differential expression (q = 3.99e-10) and enrichment (fold enrichment = 5.56, q = 9.09e-4). Our results further indicate that specific cell types, including a specific NPY-expressing GABAergic neuronal cell type, and specific G protein signaling cascades underly the cerebral effects of a loss of ultradian cortisol rhythm. Our results provide a biological mechanistic underpinning of our fMRI findings, indicating specific cell types and cascades as a target for manipulation in future experimental studies.

Transcriptomic Signatures Associated With Regional Cortical Thickness Changes in Parkinson's Disease.

Cortical atrophy is a common manifestation in Parkinson's disease (PD), particularly in advanced stages of the disease. To elucidate the molecular underpinnings of cortical thickness changes in PD, we performed an integrated analysis of brain-wide healthy transcriptomic data from the Allen Human Brain Atlas and patterns of cortical thickness based on T1-weighted anatomical MRI data of 149 PD patients and 369 controls. For this purpose, we used partial least squares regression to identify gene expression patterns correlated with cortical thickness changes. In addition, we identified gene expression patterns underlying the relationship between cortical thickness and clinical domains of PD. Our results show that genes whose expression in the healthy brain is associated with cortical thickness changes in PD are enriched in biological pathways related to sumoylation, regulation of mitotic cell cycle, mitochondrial translation, DNA damage responses, and ER-Golgi traffic. The associated pathways were highly related to each other and all belong to cellular maintenance mechanisms. The expression of genes within most pathways was negatively correlated with cortical thickness changes, showing higher expression in regions associated with decreased cortical thickness (atrophy). On the other hand, sumoylation pathways were positively correlated with cortical thickness changes, showing higher expression in regions with increased cortical thickness (hypertrophy). Our findings suggest that alterations in the balanced interplay of these mechanisms play a role in changes of cortical thickness in PD and possibly influence motor and cognitive functions.

Molecular characterization of the stress network in individuals at risk for schizophrenia.

The biological mechanisms underlying inter-individual differences in human stress reactivity remain poorly understood. We aimed to identify the molecular underpinning of aberrant neural stress sensitivity in individuals at risk for schizophrenia. Linking mRNA expression data from the Allen Human Brain Atlas to task-based fMRI revealed 201 differentially expressed genes in cortex-specific brain regions differentially activated by stress in individuals with low (healthy siblings of schizophrenia patients) or high (healthy controls) stress sensitivity. These genes are associated with stress-related psychiatric disorders (e.g. schizophrenia and anxiety) and include markers for specific neuronal populations (e.g. ADCYAP1, GABRB1, SSTR1, and TNFRSF12A), neurotransmitter receptors (e.g. GRIN3A, SSTR1, GABRB1, and HTR1E), and signaling factors that interact with the corticosteroid receptor and hypothalamic-pituitary-adrenal axis (e.g. ADCYAP1, IGSF11, and PKIA). Overall, the identified genes potentially underlie altered stress reactivity in individuals at risk for schizophrenia and other psychiatric disorders and play a role in mounting an adaptive stress response in at-risk individuals, making them potentially druggable targets for stress-related diseases.

Transcriptomic signatures of brain regional vulnerability to Parkinson's disease.

The molecular mechanisms underlying caudal-to-rostral progression of Lewy body pathology in Parkinson's disease remain poorly understood. Here, we identified transcriptomic signatures across brain regions involved in Braak Lewy body stages in non-neurological adults from the Allen Human Brain Atlas. Among the genes that are indicative of regional vulnerability, we found known genetic risk factors for Parkinson's disease: SCARB2, ELOVL7, SH3GL2, SNCA, BAP1, and ZNF184. Results were confirmed in two datasets of non-neurological subjects, while in two datasets of Parkinson's disease patients we found altered expression patterns. Co-expression analysis across vulnerable regions identified a module enriched for genes associated with dopamine synthesis and microglia, and another module related to the immune system, blood-oxygen transport, and endothelial cells. Both were highly expressed in regions involved in the preclinical stages of the disease. Finally, alterations in genes underlying these region-specific functions may contribute to the selective regional vulnerability in Parkinson's disease brains.

Co-expression Patterns between ATN1 and ATXN2 Coincide with Brain Regions Affected in Huntington's Disease.

Cytosine-adenine-guanine (CAG) repeat expansions in the coding regions of nine polyglutamine (polyQ) genes (HTT, ATXN1, ATXN2, ATXN3, CACNA1A, ATXN7, ATN1, AR, and TBP) are the cause of several neurodegenerative diseases including Huntington's disease (HD), six different spinocerebellar ataxias (SCAs), dentatorubral-pallidoluysian atrophy, and spinobulbar muscular atrophy. The expanded CAG repeat length in the causative gene is negatively related to the age-at-onset (AAO) of clinical symptoms. In addition to the expanded CAG repeat length in the causative gene, the normal CAG repeats in the other polyQ genes can affect the AAO, suggesting functional interactions between the polyQ genes. However, there is no detailed assessment of the relationships among polyQ genes in pathologically relevant brain regions. We used gene co-expression analysis to study the functional relationships among polyQ genes in different brain regions using the Allen Human Brain Atlas (AHBA), a spatial map of gene expression in the healthy brain. We constructed co-expression networks for seven anatomical brain structures, as well as a region showing a specific pattern of atrophy in HD patients detected by magnetic resonance imaging (MRI) of the brain. In this HD-associated region, we found that ATN1 and ATXN2 were co-expressed and shared co-expression partners which were enriched for DNA repair genes. We observed a similar co-expression pattern in the frontal lobe, parietal lobe, and striatum in which this relation was most pronounced. Given that the co-expression patterns for these anatomical structures were similar to those for the HD-associated region, our results suggest that their disruption is likely involved in HD pathology. Moreover, ATN1 and ATXN2 also shared many co-expressed genes with HTT, the causative gene of HD, across the brain. Although this triangular relationship among these three polyQ genes may also be dysregulated in other polyQ diseases, stronger co-expression patterns between ATN1 and ATXN2 observed in the HD-associated region, especially in the striatum, may be more specific to HD.