Spatial Transcriptomics Resolve an Emphysema-Specific Lymphoid Follicle B Cell Signature in Chronic Obstructive Pulmonary Disease.

Rationale: Within chronic obstructive pulmonary disease (COPD), emphysema is characterized by a significant yet partially understood B cell immune component. Objectives: To characterize the transcriptomic signatures from lymphoid follicles (LFs) in ever-smokers without COPD and patients with COPD with varying degrees of emphysema. Methods: Lung sections from 40 patients with COPD and ever-smokers were used for LF proteomic and transcriptomic spatial profiling. Formalin- and O.C.T.-fixed lung samples obtained from biopsies or lung explants were assessed for LF presence. Emphysema measurements were obtained from clinical chest computed tomographic scans. High-confidence transcriptional target intersection analyses were conducted to resolve emphysema-induced transcriptional networks. Measurements and Main Results: Overall, 115 LFs from ever-smokers and Global Initiative for Chronic Obstructive Lung Disease (GOLD) 1-2 and GOLD 3-4 patients were analyzed. No LFs were found in never-smokers. Differential gene expression analysis revealed significantly increased expression of LF assembly and B cell marker genes in subjects with severe emphysema. High-confidence transcriptional analysis revealed activation of an abnormal B cell activity signature in LFs (q-value = 2.56E-111). LFs from patients with GOLD 1-2 COPD with emphysema showed significantly increased expression of genes associated with antigen presentation, inflammation, and B cell activation and proliferation. LFs from patients with GOLD 1-2 COPD without emphysema showed an antiinflammatory profile. The extent of centrilobular emphysema was significantly associated with genes involved in B cell maturation and antibody production. Protein-RNA network analysis showed that LFs in emphysema have a unique signature skewed toward chronic B cell activation. Conclusions: An off-targeted B cell activation within LFs is associated with autoimmune-mediated emphysema pathogenesis.

Spatial Transcriptomics Resolve an Emphysema-specific Lymphoid Follicle B Cell Signature in COPD.

RATIONALE: Within chronic obstructive pulmonary disease (COPD), emphysema is characterized by a significant yet partially understood B cell immune component. OBJECTIVE: To characterize the transcriptomic signatures from lymphoid follicles (LFs) in ever-smokers without COPD and COPD patients with varying degrees of emphysema. METHODS: Lung sections from 40 COPD patients and ever-smokers were used for LF proteomic and transcriptomic spatial profiling. Formalin and OCT-fixed lung samples obtained from biopsies or lung explants, were assessed for LF presence. Emphysema measurements were obtained from clinical chest CT scans. High confidence transcriptional (HCT) target intersection analyses were conducted to resolve emphysema-induced transcriptional networks. MEASUREMENTS AND MAIN RESULTS: Overall, 115 LFs from ever-smokers and GOLD 1-2 and GOLD 3-4 patients were analyzed. No LFs were found in never-smokers. Differential gene expression analysis revealed significantly increased expression of LF assembly and B cell markers genes in subjects with severe emphysema. HCT analysis revealed activation of abnormal B cell activity signature in LFs (q-value: 2.56E-111). LFs from GOLD 1-2 COPD patients with emphysema showed significantly increased expression of genes associated with antigen presentation, inflammation, and B cell activation and proliferation. LFs from GOLD 1-2 COPD patients without emphysema showed an anti-inflammatory profile. The extent of centrilobular emphysema was significantly associated with genes involved in B cell maturation and antibody production. Protein-RNA network analysis showed that LFs in emphysema have a unique signature skewed towards chronic B cell activation. CONCLUSIONS: An off-targeted B cell activation within LFs is associated with autoimmune-mediated emphysema pathogenesis.

Ultra High-plex Spatial Proteogenomic Investigation of Giant Cell Glioblastoma Multiforme Immune Infiltrates Reveals Distinct Protein and RNA Expression Profiles.

A deeper understanding of complex biological processes, including tumor development and immune response, requires ultra high-plex, spatial interrogation of multiple "omes". Here we present the development and implementation of a novel spatial proteogenomic (SPG) assay on the GeoMx Digital Spatial Profiler platform with next-generation sequencing readout that enables ultra high-plex digital quantitation of proteins (>100-plex) and RNA (whole transcriptome, >18,000-plex) from a single formalin-fixed paraffin-embedded (FFPE) sample. This study highlighted the high concordance, R > 0.85 and <15% change in sensitivity between the SPG assay and the single-analyte assays on various cell lines and tissues from human and mouse. Furthermore, we demonstrate that the SPG assay was reproducible across multiple users. When used in conjunction with advanced cellular neighborhood segmentation, distinct immune or tumor RNA and protein targets were spatially resolved within individual cell subpopulations in human colorectal cancer and non-small cell lung cancer. We used the SPG assay to interrogate 23 different glioblastoma multiforme (GBM) samples across four pathologies. The study revealed distinct clustering of both RNA and protein based on pathology and anatomic location. The in-depth investigation of giant cell glioblastoma multiforme (gcGBM) revealed distinct protein and RNA expression profiles compared with that of the more common GBM. More importantly, the use of spatial proteogenomics allowed simultaneous interrogation of critical protein posttranslational modifications alongside whole transcriptomic profiles within the same distinct cellular neighborhoods. Significance: We describe ultra high-plex spatial proteogenomics; profiling whole transcriptome and high-plex proteomics on a single FFPE tissue section with spatial resolution. Investigation of gcGBM versus GBM revealed distinct protein and RNA expression profiles.

Collective effects of human genomic variation on microbiome function.

Studies of the impact of host genetics on gut microbiome composition have mainly focused on the impact of individual single nucleotide polymorphisms (SNPs) on gut microbiome composition, without considering their collective impact or the specific functions of the microbiome. To assess the aggregate role of human genetics on the gut microbiome composition and function, we apply sparse canonical correlation analysis (sCCA), a flexible, multivariate data integration method. A critical attribute of metagenome data is its sparsity, and here we propose application of a Tweedie distribution to accommodate this. We use the TwinsUK cohort to analyze the gut microbiomes and human variants of 250 individuals. Sparse CCA, or sCCA, identified SNPs in microbiome-associated metabolic traits (BMI, blood pressure) and microbiome-associated disorders (type 2 diabetes, some neurological disorders) and certain cancers. Both common and rare microbial functions such as secretion system proteins or antibiotic resistance were found to be associated with host genetics. sCCA applied to microbial species abundances found known associations such as Bifidobacteria species, as well as novel associations. Despite our small sample size, our method can identify not only previously known associations, but novel ones as well. Overall, we present a new and flexible framework for examining host-microbiome genetic interactions, and we provide a new dimension to the current debate around the role of human genetics on the gut microbiome.

Linking plasmid-based beta-lactamases to their bacterial hosts using single-cell fusion PCR.

The horizonal transfer of plasmid-encoded genes allows bacteria to adapt to constantly shifting environmental pressures, bestowing functional advantages to their bacterial hosts such as antibiotic resistance, metal resistance, virulence factors, and polysaccharide utilization. However, common molecular methods such as short- and long-read sequencing of microbiomes cannot associate extrachromosomal plasmids with the genome of the host bacterium. Alternative methods to link plasmids to host bacteria are either laborious, expensive, or prone to contamination. Here we present the One-step Isolation and Lysis PCR (OIL-PCR) method, which molecularly links plasmid-encoded genes with the bacterial 16S rRNA gene via fusion PCR performed within an emulsion. After validating this method, we apply it to identify the bacterial hosts of three clinically relevant beta-lactamases within the gut microbiomes of neutropenic patients, as they are particularly vulnerable multidrug-resistant infections. We successfully detect the known association of a multi-drug resistant plasmid with Klebsiella pneumoniae, as well as the novel associations of two low-abundance genera, Romboutsia and Agathobacter. Further investigation with OIL-PCR confirmed that our detection of Romboutsia is due to its physical association with Klebsiella as opposed to directly harboring the beta-lactamase genes. Here we put forth a robust, accessible, and high-throughput platform for sensitively surveying the bacterial hosts of mobile genes, as well as detecting physical bacterial associations such as those occurring within biofilms and complex microbial communities.

What Is Metagenomics Teaching Us, and What Is Missed?

Shotgun metagenomic sequencing has revolutionized our ability to detect and characterize the diversity and function of complex microbial communities. In this review, we highlight the benefits of using metagenomics as well as the breadth of conclusions that can be made using currently available analytical tools, such as greater resolution of species and strains across phyla and functional content, while highlighting challenges of metagenomic data analysis. Major challenges remain in annotating function, given the dearth of functional databases for environmental bacteria compared to model organisms, and the technical difficulties of metagenome assembly and phasing in heterogeneous environmental samples. In the future, improvements and innovation in technology and methodology will lead to lowered costs. Data integration using multiple technological platforms will lead to a better understanding of how to harness metagenomes. Subsequently, we will be able not only to characterize complex microbiomes but also to manipulate communities to achieve prosperous outcomes for health, agriculture, and environmental sustainability.

A Large Panel of Drosophila simulans Reveals an Abundance of Common Variants.

The rapidly expanding availability of large NGS data sets provides an opportunity to investigate population genetics at an unprecedented scale. Drosophila simulans is the sister species of the model organism Drosophila melanogaster, and is often presumed to share similar demographic history. However, previous population genetic and ecological work suggests very different signatures of selection and demography. Here, we sequence a new panel of 170 inbred genotypes of a North American population of D. simulans, a valuable complement to the DGRP and other D. melanogaster panels. We find some unexpected signatures of demography, in the form of excess intermediate frequency polymorphisms. Simulations suggest that this is possibly due to a recent population contraction and selection. We examine the outliers in the D. simulans genome determined by a haplotype test to attempt to parse the contribution of demography and selection to the patterns observed in this population. Untangling the relative contribution of demography and selection to genomic patterns of variation is challenging, however, it is clear that although D. melanogaster was thought to share demographic history with D. simulans different forces are at work in shaping genomic variation in this population of D. simulans.

Disease-specific biases in alternative splicing and tissue-specific dysregulation revealed by multitissue profiling of lymphocyte gene expression in type 1 diabetes.

Genome-wide association studies (GWAS) have identified multiple, shared allelic associations with many autoimmune diseases. However, the pathogenic contributions of variants residing in risk loci remain unresolved. The location of the majority of shared disease-associated variants in noncoding regions suggests they contribute to risk of autoimmunity through effects on gene expression in the immune system. In the current study, we test this hypothesis by applying RNA sequencing to CD4+, CD8+, and CD19+ lymphocyte populations isolated from 81 subjects with type 1 diabetes (T1D). We characterize and compare the expression patterns across these cell types for three gene sets: all genes, the set of genes implicated in autoimmune disease risk by GWAS, and the subset of these genes specifically implicated in T1D. We performed RNA sequencing and aligned the reads to both the human reference genome and a catalog of all possible splicing events developed from the genome, thereby providing a comprehensive evaluation of the roles of gene expression and alternative splicing (AS) in autoimmunity. Autoimmune candidate genes displayed greater expression specificity in the three lymphocyte populations relative to other genes, with significantly increased levels of splicing events, particularly those predicted to have substantial effects on protein isoform structure and function (e.g., intron retention, exon skipping). The majority of single-nucleotide polymorphisms within T1D-associated loci were also associated with one or more cis-expression quantitative trait loci (cis-eQTLs) and/or splicing eQTLs. Our findings highlight a substantial, and previously underrecognized, role for AS in the pathogenesis of autoimmune disorders and particularly for T1D.

Neurons That Underlie Drosophila melanogaster Reproductive Behaviors: Detection of a Large Male-Bias in Gene Expression in fruitless-Expressing Neurons.

Male and female reproductive behaviors in Drosophila melanogaster are vastly different, but neurons that express sex-specifically spliced fruitless transcripts (fru P1) underlie these behaviors in both sexes. How this set of neurons can generate such different behaviors between the two sexes is an unresolved question. A particular challenge is that fru P1-expressing neurons comprise only 2-5% of the adult nervous system, and so studies of adult head tissue or whole brain may not reveal crucial differences. Translating Ribosome Affinity Purification (TRAP) identifies the actively translated pool of mRNAs from fru P1-expressing neurons, allowing a sensitive, cell-type-specific assay. We find four times more male-biased than female-biased genes in TRAP mRNAs from fru P1-expressing neurons. This suggests a potential mechanism to generate dimorphism in behavior. The male-biased genes may direct male behaviors by establishing cell fate in a similar context of gene expression observed in females. These results suggest a possible global mechanism for how distinct behaviors can arise from a shared set of neurons.

Sex Differences in Drosophila Somatic Gene Expression: Variation and Regulation by doublesex.

Sex differences in gene expression have been widely studied in Drosophila melanogaster Sex differences vary across strains, but many molecular studies focus on only a single strain, or on genes that show sexually dimorphic expression in many strains. How extensive variability is and whether this variability occurs among genes regulated by sex determination hierarchy terminal transcription factors is unknown. To address these questions, we examine differences in sexually dimorphic gene expression between two strains in Drosophila adult head tissues. We also examine gene expression in doublesex (dsx) mutant strains to determine which sex-differentially expressed genes are regulated by DSX, and the mode by which DSX regulates expression. We find substantial variation in sex-differential expression. The sets of genes with sexually dimorphic expression in each strain show little overlap. The prevalence of different DSX regulatory modes also varies between the two strains. Neither the patterns of DSX DNA occupancy, nor mode of DSX regulation explain why some genes show consistent sex-differential expression across strains. We find that the genes identified as regulated by DSX in this study are enriched with known sites of DSX DNA occupancy. Finally, we find that sex-differentially expressed genes and genes regulated by DSX are highly enriched on the fourth chromosome. These results provide insights into a more complete pool of potential DSX targets, as well as revealing the molecular flexibility of DSX regulation.