Cigarette Smoke Specifically Affects Small Airway Epithelial Cell Populations and Triggers the Expansion of Inflammatory and Squamous Differentiation Associated Basal Cells.

Smoking is a major risk factor for chronic obstructive pulmonary disease (COPD) and causes remodeling of the small airways. However, the exact smoke-induced effects on the different types of small airway epithelial cells (SAECs) are poorly understood. Here, using air-liquid interface (ALI) cultures, single-cell RNA-sequencing reveals previously unrecognized transcriptional heterogeneity within the small airway epithelium and cell type-specific effects upon acute and chronic cigarette smoke exposure. Smoke triggers detoxification and inflammatory responses and aberrantly activates and alters basal cell differentiation. This results in an increase of inflammatory basal-to-secretory cell intermediates and, particularly after chronic smoke exposure, a massive expansion of a rare inflammatory and squamous metaplasia associated KRT6A+ basal cell state and an altered secretory cell landscape. ALI cultures originating from healthy non-smokers and COPD smokers show similar responses to cigarette smoke exposure, although an increased pro-inflammatory profile is conserved in the latter. Taken together, the in vitro models provide high-resolution insights into the smoke-induced remodeling of the small airways resembling the pathological processes in COPD airways. The data may also help to better understand other lung diseases including COVID-19, as the data reflect the smoke-dependent variable induction of SARS-CoV-2 entry factors across SAEC populations.

Monocyte-derived alveolar macrophages are key drivers of smoke-induced lung inflammation and tissue remodeling.

Smoking is a leading risk factor of chronic obstructive pulmonary disease (COPD), that is characterized by chronic lung inflammation, tissue remodeling and emphysema. Although inflammation is critical to COPD pathogenesis, the cellular and molecular basis underlying smoking-induced lung inflammation and pathology remains unclear. Using murine smoke models and single-cell RNA-sequencing, we show that smoking establishes a self-amplifying inflammatory loop characterized by an influx of molecularly heterogeneous neutrophil subsets and excessive recruitment of monocyte-derived alveolar macrophages (MoAM). In contrast to tissue-resident AM, MoAM are absent in homeostasis and characterized by a pro-inflammatory gene signature. Moreover, MoAM represent 46% of AM in emphysematous mice and express markers causally linked to emphysema. We also demonstrate the presence of pro-inflammatory and tissue remodeling associated MoAM orthologs in humans that are significantly increased in emphysematous COPD patients. Inhibition of the IRAK4 kinase depletes a rare inflammatory neutrophil subset, diminishes MoAM recruitment, and alleviates inflammation in the lung of cigarette smoke-exposed mice. This study extends our understanding of the molecular signaling circuits and cellular dynamics in smoking-induced lung inflammation and pathology, highlights the functional consequence of monocyte and neutrophil recruitment, identifies MoAM as key drivers of the inflammatory process, and supports their contribution to pathological tissue remodeling.

Characterization and improvement of RNA-Seq precision in quantitative transcript expression profiling.

MOTIVATION: Measurement precision determines the power of any analysis to reliably identify significant signals, such as in screens for differential expression, independent of whether the experimental design incorporates replicates or not. With the compilation of large-scale RNA-Seq datasets with technical replicate samples, however, we can now, for the first time, perform a systematic analysis of the precision of expression level estimates from massively parallel sequencing technology. This then allows considerations for its improvement by computational or experimental means. RESULTS: We report on a comprehensive study of target identification and measurement precision, including their dependence on transcript expression levels, read depth and other parameters. In particular, an impressive recall of 84% of the estimated true transcript population could be achieved with 331 million 50 bp reads, with diminishing returns from longer read lengths and even less gains from increased sequencing depths. Most of the measurement power (75%) is spent on only 7% of the known transcriptome, however, making less strongly expressed transcripts harder to measure. Consequently, <30% of all transcripts could be quantified reliably with a relative error<20%. Based on established tools, we then introduce a new approach for mapping and analysing sequencing reads that yields substantially improved performance in gene expression profiling, increasing the number of transcripts that can reliably be quantified to over 40%. Extrapolations to higher sequencing depths highlight the need for efficient complementary steps. In discussion we outline possible experimental and computational strategies for further improvements in quantification precision. CONTACT: rnaseq10@boku.ac.at

Hybridization thermodynamics of NimbleGen microarrays.

BACKGROUND: While microarrays are the predominant method for gene expression profiling, probe signal variation is still an area of active research. Probe signal is sequence dependent and affected by probe-target binding strength and the competing formation of probe-probe dimers and secondary structures in probes and targets. RESULTS: We demonstrate the benefits of an improved model for microarray hybridization and assess the relative contributions of the probe-target binding strength and the different competing structures. Remarkably, specific and unspecific hybridization were apparently driven by different energetic contributions: For unspecific hybridization, the melting temperature Tm was the best predictor of signal variation. For specific hybridization, however, the effective interaction energy that fully considered competing structures was twice as powerful a predictor of probe signal variation. We show that this was largely due to the effects of secondary structures in the probe and target molecules. The predictive power of the strength of these intramolecular structures was already comparable to that of the melting temperature or the free energy of the probe-target duplex. CONCLUSIONS: This analysis illustrates the importance of considering both the effects of probe-target binding strength and the different competing structures. For specific hybridization, the secondary structures of probe and target molecules turn out to be at least as important as the probe-target binding strength for an understanding of the observed microarray signal intensities. Besides their relevance for the design of new arrays, our results demonstrate the value of improving thermodynamic models for the read-out and interpretation of microarray signals.

DMSO cryopreservation is the method of choice to preserve cells for droplet-based single-cell RNA sequencing.

Combining single-cell RNA sequencing (scRNA-seq) with upstream cell preservation procedures such as cryopreservation or methanol fixation has recently become more common. By separating cell handling and preparation, from downstream library generation, scRNA-seq workflows are more flexible and manageable. However, the inherent transcriptomic changes associated with cell preservation and how they may bias further downstream analysis remain unknown. Here, we present a side-by-side droplet-based scRNA-seq analysis, comparing the gold standard - fresh cells - to three different cell preservation workflows: dimethyl sulfoxide based cryopreservation, methanol fixation and CellCover reagent. Cryopreservation proved to be the most robust protocol, maximizing both cell integrity and low background ambient RNA. Importantly, gene expression profiles from fresh cells correlated most with those of cryopreserved cells. Such similarities were consistently observed across the tested cell lines (R ≥ 0.97), monocyte-derived macrophages (R = 0.97) and immune cells (R = 0.99). In contrast, both methanol fixation and CellCover preservation showed an increased ambient RNA background and an overall lower gene expression correlation to fresh cells. Thus, our results demonstrate the superiority of cryopreservation over other cell preservation methods. We expect our comparative study to provide single-cell omics researchers invaluable support when integrating cell preservation into their scRNA-seq studies.

Novel insights into the unfolded protein response using Pichia pastoris specific DNA microarrays.

BACKGROUND: DNA Microarrays are regarded as a valuable tool for basic and applied research in microbiology. However, for many industrially important microorganisms the lack of commercially available microarrays still hampers physiological research. Exemplarily, our understanding of protein folding and secretion in the yeast Pichia pastoris is presently widely dependent on conclusions drawn from analogies to Saccharomyces cerevisiae. To close this gap for a yeast species employed for its high capacity to produce heterologous proteins, we developed full genome DNA microarrays for P. pastoris and analyzed the unfolded protein response (UPR) in this yeast species, as compared to S. cerevisiae. RESULTS: By combining the partially annotated gene list of P. pastoris with de novo gene finding a list of putative open reading frames was generated for which an oligonucleotide probe set was designed using the probe design tool TherMODO (a thermodynamic model-based oligoset design optimizer). To evaluate the performance of the novel array design, microarrays carrying the oligo set were hybridized with samples from treatments with dithiothreitol (DTT) or a strain overexpressing the UPR transcription factor HAC1, both compared with a wild type strain in normal medium as untreated control. DTT treatment was compared with literature data for S. cerevisiae, and revealed similarities, but also important differences between the two yeast species. Overexpression of HAC1, the most direct control for UPR genes, resulted in significant new understanding of this important regulatory pathway in P. pastoris, and generally in yeasts. CONCLUSION: The differences observed between P. pastoris and S. cerevisiae underline the importance of DNA microarrays for industrial production strains. P. pastoris reacts to DTT treatment mainly by the regulation of genes related to chemical stimulus, electron transport and respiration, while the overexpression of HAC1 induced many genes involved in translation, ribosome biogenesis, and organelle biosynthesis, indicating that the regulatory events triggered by DTT treatment only partially overlap with the reactions to overexpression of HAC1. The high reproducibility of the results achieved with two different oligo sets is a good indication for their robustness, and underlines the importance of less stringent selection of regulated features, in order to avoid a large number of false negative results.

Single amino acid repeats in signal peptides.

There has been an increasing interest in single amino acid repeats ever since it was shown that these are the cause of a variety of diseases. Although a systematic study of single amino acid repeats is challenging, they have subsequently been implicated in a number of functional roles. In general surveys, leucine runs were among the most frequent. In the present study, we present a detailed investigation of repeats in signal peptides of secreted and type I membrane proteins in comparison with their mature parts. We focus on eukaryotic species because single amino acid repeats are generally rather rare in archaea and bacteria. Our analysis of over 100 species shows that repeats of leucine (but not of other hydrophobic amino acids) are over-represented in signal peptides. This trend is most pronounced in higher eukaryotes, particularly in mammals. In the human proteome, although less than one-fifth of all proteins have a signal peptide, approximately two-thirds of all leucine repeats are located in these transient regions. Signal peptides are cleaved early from the growing polypeptide chain and then degraded rapidly. This may explain why leucine repeats, which can be toxic, are tolerated at such high frequencies. The substantial fraction of proteins affected by the strong enrichment of repeats in these transient segments highlights the bias that they can introduce for systematic analyses of protein sequences. In contrast to a general lack of conservation of single amino acid repeats, leucine repeats were found to be more conserved than the remaining signal peptide regions, indicating that they may have an as yet unknown functional role.