miRglmm: a generalized linear mixed model of isomiR-level counts improves estimation of miRNA-level differential expression and uncovers variable differential expression between isomiRs.

MicroRNA-seq data is produced by aligning small RNA sequencing reads of different miRNA transcript isoforms, called isomiRs, to known microRNAs. Aggregation to microRNA-level counts discards information and violates core assumptions of differential expression (DE) methods developed for mRNA-seq data. We establish miRglmm, a DE method for microRNA-seq data, that uses a generalized linear mixed model of isomiR-level counts, facilitating detection of miRNA with differential expression or differential isomiR usage. We demonstrate that miRglmm outperforms current DE methods in estimating DE for miRNA, whether or not there is significant isomiR variability, and simultaneously provides estimates of isomiR-level DE.

Association of disease severity and genetic variation during primary Respiratory Syncytial Virus infections.

BACKGROUND: Respiratory Syncytial Virus (RSV) disease in young children ranges from mild cold symptoms to severe symptoms that require hospitalization and sometimes result in death. Studies have shown a statistical association between RSV subtype or phylogenic lineage and RSV disease severity, although these results have been inconsistent. Associations between variation within RSV gene coding regions or residues and RSV disease severity has been largely unexplored. METHODS: Nasal swabs from children (< 8 months-old) infected with RSV in Rochester, NY between 1977-1998 clinically presenting with either mild or severe disease during their first cold-season were used. Whole-genome RSV sequences were obtained using overlapping PCR and next-generation sequencing. Both whole-genome phylogenetic and non-phylogenetic statistical approaches were performed to associate RSV genotype with disease severity. RESULTS: The RSVB subtype was statistically associated with disease severity. A significant association between phylogenetic clustering of mild/severe traits and disease severity was also found. GA1 clade sequences were associated with severe disease while GB1 was significantly associated with mild disease. Both G and M2-2 gene variation was significantly associated with disease severity. We identified 16 residues in the G gene and 3 in the M2-2 RSV gene associated with disease severity. CONCLUSION: These results suggest that phylogenetic lineage and the genetic variability in G or M2-2 genes of RSV may contribute to disease severity in young children undergoing their first infection.

Patterns of Unwanted Biological and Technical Expression Variation Among 49 Human Tissues.

Tissue gene expression studies are impacted by biological and technical sources of variation, which can be broadly classified into wanted and unwanted variation. The latter, if not addressed, results in misleading biological conclusions. Methods have been proposed to reduce unwanted variation, such as normalization and batch correction. A more accurate understanding of all causes of variation could significantly improve the ability of these methods to remove unwanted variation while retaining variation corresponding to the biological question of interest. We used 17,282 samples from 49 human tissues in the Genotype-Tissue Expression data set (v8) to investigate patterns and causes of expression variation. Transcript expression was transformed to z-scores, and only the most variable 2% of transcripts were evaluated and clustered based on coexpression patterns. Clustered gene sets were assigned to different biological or technical causes based on histologic appearances and metadata elements. We identified 522 variable transcript clusters (median: 11 per tissue) among the samples. Of these, 63% were confidently explained, 16% were likely explained, 7% were low confidence explanations, and 14% had no clear cause. Histologic analysis annotated 46 clusters. Other common causes of variability included sex, sequencing contamination, immunoglobulin diversity, and compositional tissue differences. Less common biological causes included death interval (Hardy score), disease status, and age. Technical causes included blood draw timing and harvesting differences. Many of the causes of variation in bulk tissue expression were identifiable in the Tabula Sapiens data set of single-cell expression. This is among the largest explorations of the underlying sources of tissue expression variation. It uncovered expected and unexpected causes of variable gene expression and demonstrated the utility of matched histologic specimens. It further demonstrated the value of acquiring meaningful tissue harvesting metadata elements to use for improved normalization, batch correction, and analysis of both bulk and single-cell RNA-seq data.

A model-based hierarchical Bayesian approach to Sholl analysis.

MOTIVATION: Due to the link between microglial morphology and function, morphological changes in microglia are frequently used to identify pathological immune responses in the central nervous system. In the absence of pathology, microglia are responsible for maintaining homeostasis, and their morphology can be indicative of how the healthy brain behaves in the presence of external stimuli and genetic differences. Despite recent interest in high throughput methods for morphological analysis, Sholl analysis is still widely used for quantifying microglia morphology via imaging data. Often, the raw data are naturally hierarchical, minimally including many cells per image and many images per animal. However, existing methods for performing downstream inference on Sholl data rely on truncating this hierarchy so rudimentary statistical testing procedures can be used. RESULTS: To fill this longstanding gap, we introduce a parametric hierarchical Bayesian model-based approach for analyzing Sholl data, so that inference can be performed without aggressive reduction of otherwise very rich data. We apply our model to real data and perform simulation studies comparing the proposed method with a popular alternative. AVAILABILITY AND IMPLEMENTATION: Software to reproduce the results presented in this article is available at: https://github.com/vonkaenelerik/hierarchical_sholl. An R package implementing the proposed models is available at: https://github.com/vonkaenelerik/ShollBayes.

Developmental Ethanol Exposure Impacts Purkinje Cells but Not Microglia in the Young Adult Cerebellum.

Fetal alcohol spectrum disorders (FASD) caused by developmental ethanol exposure lead to cerebellar impairments, including motor problems, decreased cerebellar weight, and cell death. Alterations in the sole output of the cerebellar cortex, Purkinje cells, and central nervous system immune cells, microglia, have been reported in animal models of FASD. To determine how developmental ethanol exposure affects adult cerebellar microglia and Purkinje cells, we used a human third-trimester binge exposure model in which mice received ethanol or saline from postnatal (P) days 4-9. In adolescence, cerebellar cranial windows were implanted and mice were aged to young adulthood for examination of microglia and Purkinje cells in vivo with two-photon imaging or in fixed tissue. Ethanol had no effect on microglia density, morphology, dynamics, or injury response. However, Purkinje cell linear frequency was reduced by ethanol. Microglia-Purkinje cell interactions in the Purkinje Cell Layer were altered in females compared to males. Overall, developmental ethanol exposure had few effects on cerebellar microglia in young adulthood and Purkinje cells appeared to be more susceptible to its effects.

Smad4 restricts injury-provoked biliary proliferation and carcinogenesis.

Cholangiocarcinoma (CCA) is a deadly and heterogeneous type of cancer characterized by a spectrum of epidemiologic associations as well as genetic and epigenetic alterations. We seek to understand how these features inter-relate in the earliest phase of cancer development and through the course of disease progression. For this, we studied murine models of liver injury integrating the most commonly occurring gene mutations of CCA - including Kras, Tp53, Arid1a and Smad4 - as well as murine hepatobiliary cancer models and derived primary cell lines based on these mutations. Among commonly mutated genes in CCA, we found that Smad4 functions uniquely to restrict reactive cholangiocyte expansion to liver injury through restraint of the proliferative response. Inactivation of Smad4 accelerates carcinogenesis, provoking pre-neoplastic biliary lesions and CCA development in an injury setting. Expression analyses of Smad4-perturbed reactive cholangiocytes and CCA lines demonstrated shared enriched pathways, including cell-cycle regulation, MYC signaling and oxidative phosphorylation, suggesting that Smad4 may act via these mechanisms to regulate cholangiocyte proliferation and progression to CCA. Overall, we showed that TGFß/SMAD4 signaling serves as a critical barrier restraining cholangiocyte expansion and malignant transformation in states of biliary injury.

Brooklyn plots to identify co-expression dysregulation in single cell sequencing.

Altered open chromatin regions, impacting gene expression, is a feature of some human disorders. We discovered it is possible to detect global changes in genomically-related adjacent gene co-expression within single cell RNA sequencing (scRNA-seq) data. We built a software package to generate and test non-randomness using 'Brooklyn plots' to identify the percent of genes significantly co-expressed from the same chromosome in ∼10 MB intervals across the genome. These plots establish an expected low baseline of co-expression in scRNA-seq from most cell types, but, as seen in dilated cardiomyopathy cardiomyocytes, altered patterns of open chromatin appear. These may relate to larger regions of transcriptional bursting, observable in single cell, but not bulk datasets.

Developmental ethanol exposure has minimal impact on cerebellar microglial dynamics, morphology, and interactions with Purkinje cells during adolescence.

Introduction: Fetal alcohol spectrum disorders (FASD) are the most common cause of non-heritable, preventable mental disability, occurring in almost 5% of births in the United States. FASD lead to physical, behavioral, and cognitive impairments, including deficits related to the cerebellum. There is no known cure for FASD and their mechanisms remain poorly understood. To better understand these mechanisms, we examined the cerebellum on a cellular level by studying microglia, the principal immune cells of the central nervous system, and Purkinje cells, the sole output of the cerebellum. Both cell types have been shown to be affected in models of FASD, with increased cell death, immune activation of microglia, and altered firing in Purkinje cells. While ethanol administered in adulthood can acutely depress the dynamics of the microglial process arbor, it is unknown how developmental ethanol exposure impacts microglia dynamics and their interactions with Purkinje cells in the long term. Methods: To address this question, we used a mouse model of human 3rd trimester exposure, whereby L7cre/Ai9+/-/Cx3cr1G/+ mice (with fluorescently labeled microglia and Purkinje cells) of both sexes were subcutaneously treated with a binge-level dose of ethanol (5.0 g/kg/day) or saline from postnatal days 4-9. Cranial windows were implanted in adolescent mice above the cerebellum to examine the long-term effects of developmental ethanol exposure on cerebellar microglia and Purkinje cell interactions using in vivo two-photon imaging. Results: We found that cerebellar microglia dynamics and morphology were not affected after developmental ethanol exposure. Microglia dynamics were also largely unaltered with respect to how they interact with Purkinje cells, although subtle changes in these interactions were observed in females in the molecular layer of the cerebellum. Discussion: This work suggests that there are limited in vivo long-term effects of ethanol exposure on microglia morphology, dynamics, and neuronal interactions, so other avenues of research may be important in elucidating the mechanisms of FASD.

Patterns of unwanted biological and technical expression variation across 49 human tissues.

All tissue-based gene expression studies are impacted by biological and technical sources of variation. Numerous methods are used to normalize and batch correct these datasets. A more accurate understanding of all causes of variation could further optimize these approaches. We used 17,282 samples from 49 tissues in the Genotype Tissue Expression (GTEx) dataset (v8) to investigate patterns and causes of expression variation. Transcript expression was normalized to Z-scores and only the most variable 2% of transcripts were evaluated and clustered based on co-expression patterns. Clustered gene sets were solved to different biological or technical causes related to metadata elements and histologic images. We identified 522 variable transcript clusters (median 11 per tissue) across the samples. Of these, 64% were confidently explained, 15% were likely explained, 7% were low confidence explanations and 14% had no clear cause. Common causes included sex, sequencing contamination, immunoglobulin diversity, and compositional tissue differences. Less common biological causes included death interval (Hardy score), muscle atrophy, diabetes status, and menopause. Technical causes included brain pH and harvesting differences. Many of the causes of variation in bulk tissue expression were identifiable in the Tabula Sapiens dataset of single cell expression. This is the largest exploration of the underlying sources of tissue expression variation. It uncovered expected and unexpected causes of variable gene expression. These identified sources of variation will inform which metadata to acquire with tissue harvesting and can be used to improve normalization, batch correction, and analysis of both bulk and single cell RNA-seq data.

A Model-Based Hierarchical Bayesian Approach to Sholl Analysis.

Due to the link between microglial morphology and function, morphological changes in microglia are frequently used to identify pathological immune responses in the central nervous system. In the absence of pathology, microglia are responsible for maintaining homeostasis, and their morphology can be indicative of how the healthy brain behaves in the presence of external stimuli and genetic differences. Despite recent interest in high throughput methods for morphological analysis, Sholl analysis is still the gold standard for quantifying microglia morphology via imaging data. Often, the raw data are naturally hierarchical, minimally including many cells per image and many images per animal. However, existing methods for performing downstream inference on Sholl data rely on truncating this hierarchy so rudimentary statistical testing procedures can be used. To fill this longstanding gap, we introduce a fully parametric model-based approach for analyzing Sholl data. We generalize our model to a hierarchical Bayesian framework so that inference can be performed without aggressive reduction of otherwise very rich data. We apply our model to three real data examples and perform simulation studies comparing the proposed method with a popular alternative.

Protocol to use TopNet for gene regulatory network modeling using gene expression data from perturbation experiments.

Inference of gene regulatory networks from gene perturbation experiments is the most reliable approach for investigating interdependence between genes. Here, we describe the initial gene perturbations, expression measurements, and preparation steps, followed by network modeling using TopNet. Summarization and visualization of the estimated networks and optional genetic testing of dependencies revealed by the network model are demonstrated. While developed for gene perturbation experiments, TopNet models data in which nodes are both perturbed and measured. For complete details on the use and execution of this protocol, please refer to McMurray et al. (2021).

A curated human cellular microRNAome based on 196 primary cell types.

BACKGROUND: An incomplete picture of the expression distribution of microRNAs (miRNAs) across human cell types has long hindered our understanding of this important regulatory class of RNA. With the continued increase in available public small RNA sequencing datasets, there is an opportunity to more fully understand the general distribution of miRNAs at the cell level. RESULTS: From the NCBI Sequence Read Archive, we obtained 6,054 human primary cell datasets and processed 4,184 of them through the miRge3.0 small RNA sequencing alignment software. This dataset was curated down, through shared miRNA expression patterns, to 2,077 samples from 196 unique cell types derived from 175 separate studies. Of 2,731 putative miRNAs listed in miRBase (v22.1), 2,452 (89.8%) were detected. Among reasonably expressed miRNAs, 108 were designated as cell specific/near specific, 59 as infrequent, 52 as frequent, 54 as near ubiquitous, and 50 as ubiquitous. The complexity of cellular microRNA expression estimates recapitulates tissue expression patterns and informs on the miRNA composition of plasma. CONCLUSIONS: This study represents the most complete reference, to date, of miRNA expression patterns by primary cell type. The data are available through the human cellular microRNAome track at the UCSC Genome Browser (https://genome.ucsc.edu/cgi-bin/hgHubConnect) and an R/Bioconductor package (https://bioconductor.org/packages/microRNAome/).

A systems genomics approach uncovers molecular associates of RSV severity.

Respiratory syncytial virus (RSV) infection results in millions of hospitalizations and thousands of deaths each year. Variations in the adaptive and innate immune response appear to be associated with RSV severity. To investigate the host response to RSV infection in infants, we performed a systems-level study of RSV pathophysiology, incorporating high-throughput measurements of the peripheral innate and adaptive immune systems and the airway epithelium and microbiota. We implemented a novel multi-omic data integration method based on multilayered principal component analysis, penalized regression, and feature weight back-propagation, which enabled us to identify cellular pathways associated with RSV severity. In both airway and immune cells, we found an association between RSV severity and activation of pathways controlling Th17 and acute phase response signaling, as well as inhibition of B cell receptor signaling. Dysregulation of both the humoral and mucosal response to RSV may play a critical role in determining illness severity.

Gene network modeling via TopNet reveals functional dependencies between diverse tumor-critical mediator genes.

Malignant cell transformation and the underlying reprogramming of gene expression require the cooperation of multiple oncogenic mutations. This cooperation is reflected in the synergistic regulation of non-mutant downstream genes, so-called cooperation response genes (CRGs). CRGs affect diverse hallmark features of cancer cells and are not known to be functionally connected. However, they act as critical mediators of the cancer phenotype at an unexpectedly high frequency >50%, as indicated by genetic perturbations. Here, we demonstrate that CRGs function within a network of strong genetic interdependencies that are critical to the malignant state. Our network modeling methodology, TopNet, takes the approach of incorporating uncertainty in the underlying gene perturbation data and can identify non-linear gene interactions. In the dense space of gene connectivity, TopNet reveals a sparse topological gene network architecture, effectively pinpointing functionally relevant gene interactions. Thus, among diverse potential applications, TopNet has utility for identification of non-mutant targets for cancer intervention.

Tissue, age, sex, and disease patterns of matrisome expression in GTEx transcriptome data.

The extracellular matrix (ECM) has historically been explored through proteomic methods. Whether or not global transcriptomics can yield meaningful information on the human matrisome is unknown. Gene expression data from 17,382 samples across 52 tissues, were obtained from the Genotype-Tissue Expression (GTEx) project. Additional datasets were obtained from The Cancer Genome Atlas (TCGA) program and the Gene Expression Omnibus for comparisons. Gene expression levels generally matched proteome-derived matrisome expression patterns. Further, matrisome gene expression properly clustered tissue types, with some matrisome genes including SERPIN family members having tissue-restricted expression patterns. Deeper analyses revealed 382 gene transcripts varied by age and 315 varied by sex in at least one tissue, with expression correlating with digitally imaged histologic tissue features. A comparison of TCGA tumor, TCGA adjacent normal and GTEx normal tissues demonstrated robustness of the GTEx samples as a generalized matrix control, while also determining a common primary tumor matrisome. Additionally, GTEx tissues served as a useful non-diseased control in a separate study of idiopathic pulmonary fibrosis (IPF) matrix changes, while identifying 22 matrix genes upregulated in IPF. Altogether, these findings indicate that the transcriptome, in general, and GTEx in particular, has value in understanding the state of organ ECM.

The effect of air pollution on the transcriptomics of the immune response to respiratory infection.

Combustion related particulate matter air pollution (PM) is associated with an increased risk of respiratory infections in adults. The exact mechanism underlying this association has not been determined. We hypothesized that increased concentrations of combustion related PM would result in dysregulation of the innate immune system. This epidemiological study includes 111 adult patients hospitalized with respiratory infections who underwent transcriptional analysis of their peripheral blood. We examined the association between gene expression at the time of hospitalization and ambient measurements of particulate air pollutants in the 28 days prior to hospitalization. For each pollutant and time lag, gene-specific linear models adjusting for infection type were fit using LIMMA (Linear Models For Microarray Data), and pathway/gene set analyses were performed using the CAMERA (Correlation Adjusted Mean Rank) program. Comparing patients with viral and/or bacterial infection, the expression patterns associated with air pollution exposure differed. Adjusting for the type of infection, increased concentrations of Delta-C (a marker of biomass smoke) and other PM were associated with upregulation of iron homeostasis and protein folding. Increased concentrations of black carbon (BC) were associated with upregulation of viral related gene pathways and downregulation of pathways related to antigen presentation. The pollutant/pathway associations differed by lag time and by type of infection. This study suggests that the effect of air pollution on the pathogenesis of respiratory infection may be pollutant, timing, and infection specific.

The role of P2Y12 in the kinetics of microglial self-renewal and maturation in the adult visual cortex in vivo.

Microglia are the brain's resident immune cells with a tremendous capacity to autonomously self-renew. Because microglial self-renewal has largely been studied using static tools, its mechanisms and kinetics are not well understood. Using chronic in vivo two-photon imaging in awake mice, we confirm that cortical microglia show limited turnover and migration under basal conditions. Following depletion, however, microglial repopulation is remarkably rapid and is sustained by the dynamic division of remaining microglia, in a manner that is largely independent of signaling through the P2Y12 receptor. Mathematical modeling of microglial division demonstrates that the observed division rates can account for the rapid repopulation observed in vivo. Additionally, newly born microglia resemble mature microglia within days of repopulation, although morphological maturation is different in newly born microglia in P2Y12 knock out mice. Our work suggests that microglia rapidly locally and that newly born microglia do not recapitulate the slow maturation seen in development but instead take on mature roles in the CNS.

Single cell RNA-seq analysis of the flexor digitorum brevis mouse myofibers.

BACKGROUND: Skeletal muscle myofibers can be separated into functionally distinct cell types that differ in gene and protein expression. Current single cell expression data is generally based upon single nucleus RNA, rather than whole myofiber material. We examined if a whole-cell flow sorting approach could be applied to perform single cell RNA-seq (scRNA-seq) in a single muscle type. METHODS: We performed deep, whole cell, scRNA-seq on intact and fragmented skeletal myofibers from the mouse fast-twitch flexor digitorum brevis muscle utilizing a flow-gated method of large cell isolation. We performed deep sequencing of 763 intact and fragmented myofibers. RESULTS: Quality control metrics across the different gates indicated only 171 of these cells were optimal, with a median read count of 239,252 and an average of 12,098 transcripts per cell. scRNA-seq identified three clusters of myofibers (a slow/fast 2A cluster and two fast 2X clusters). Comparison to a public skeletal nuclear RNA-seq dataset demonstrated a diversity in transcript abundance by method. RISH validated multiple genes across fast and slow twitch skeletal muscle types. CONCLUSION: This study introduces and validates a method to isolate intact skeletal muscle myofibers to generate deep expression patterns and expands the known repertoire of fiber-type-specific genes.

Airway gene-expression classifiers for respiratory syncytial virus (RSV) disease severity in infants.

BACKGROUND: A substantial number of infants infected with RSV develop severe symptoms requiring hospitalization. We currently lack accurate biomarkers that are associated with severe illness. METHOD: We defined airway gene expression profiles based on RNA sequencing from nasal brush samples from 106 full-tem previously healthy RSV infected subjects during acute infection (day 1-10 of illness) and convalescence stage (day 28 of illness). All subjects were assigned a clinical illness severity score (GRSS). Using AIC-based model selection, we built a sparse linear correlate of GRSS based on 41 genes (NGSS1). We also built an alternate model based upon 13 genes associated with severe infection acutely but displaying stable expression over time (NGSS2). RESULTS: NGSS1 is strongly correlated with the disease severity, demonstrating a naïve correlation (ρ) of ρ = 0.935 and cross-validated correlation of 0.813. As a binary classifier (mild versus severe), NGSS1 correctly classifies disease severity in 89.6% of the subjects following cross-validation. NGSS2 has slightly less, but comparable, accuracy with a cross-validated correlation of 0.741 and classification accuracy of 84.0%. CONCLUSION: Airway gene expression patterns, obtained following a minimally-invasive procedure, have potential utility for development of clinically useful biomarkers that correlate with disease severity in primary RSV infection.