Explore the vast range of titles available.

Spatial transcriptomics (ST) technology through in situ capturing has enabled topographical gene expression profiling of tumor tissues. However, each capturing spot may contain diverse immune and malignant cells, with different cell densities across tissue regions. Cell type deconvolution in tumor ST data remains challenging for existing methods designed to decompose general ST or bulk tumor data. We develop the Spatial Cellular Estimator for Tumors (SpaCET) to infer cell identities from tumor ST data. SpaCET first estimates cancer cell abundance by integrating a gene pattern dictionary of copy number alterations and expression changes in common malignancies. A constrained regression model then calibrates local cell densities and determines immune and stromal cell lineage fractions. SpaCET provides higher accuracy than existing methods based on simulation and real ST data with matched double-blind histopathology annotations as ground truth. Further, coupling cell fractions with ligand-receptor coexpression analysis, SpaCET reveals how intercellular interactions at the tumor-immune interface promote cancer progression. Cell type deconvolution in tumor spatial transcriptomics (ST) data remains challenging. Here, the authors develop Spatial Cellular Estimator for Tumors (SpaCET) to infer cell types and intercellular interactions from ST data in cancer across different platforms, with improved performance over similar methods.

Angiotensin-converting enzyme 2 (ACE2) and accessory proteases (TMPRSS2 and CTSL) are needed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cellular entry, and their expression may shed light on viral tropism and impact across the body. We assessed the cell-type-specific expression of ACE2 , TMPRSS2 and CTSL across 107 single-cell RNA-sequencing studies from different tissues. ACE2 , TMPRSS2 and CTSL are coexpressed in specific subsets of respiratory epithelial cells in the nasal passages, airways and alveoli, and in cells from other organs associated with coronavirus disease 2019 (COVID-19) transmission or pathology. We performed a meta-analysis of 31 lung single-cell RNA-sequencing studies with 1,320,896 cells from 377 nasal, airway and lung parenchyma samples from 228 individuals. This revealed cell-type-specific associations of age, sex and smoking with expression levels of ACE2 , TMPRSS2 and CTSL . Expression of entry factors increased with age and in males, including in airway secretory cells and alveolar type 2 cells. Expression programs shared by ACE2 + TMPRSS2 + cells in nasal, lung and gut tissues included genes that may mediate viral entry, key immune functions and epithelial–macrophage cross-talk, such as genes involved in the interleukin-6, interleukin-1, tumor necrosis factor and complement pathways. Cell-type-specific expression patterns may contribute to the pathogenesis of COVID-19, and our work highlights putative molecular pathways for therapeutic intervention. An integrated analysis of over 100 single-cell and single-nucleus transcriptomics studies illustrates severe acute respiratory syndrome coronavirus 2 viral entry gene coexpression patterns across different human tissues, and shows association of age, smoking status and sex with viral entry gene expression in respiratory cell populations.

Salmon is a computational tool that uses sample-specific models and a dual-phase inference procedure to correct biases in RNA-seq data and rapidly quantify transcript abundances. We introduce Salmon, a lightweight method for quantifying transcript abundance from RNA–seq reads. Salmon combines a new dual-phase parallel inference algorithm and feature-rich bias models with an ultra-fast read mapping procedure. It is the first transcriptome-wide quantifier to correct for fragment GC-content bias, which, as we demonstrate here, substantially improves the accuracy of abundance estimates and the sensitivity of subsequent differential expression analysis.

Gene expression fundamentally shapes the structural and functional architecture of the human brain. Open-access transcriptomic datasets like the Allen Human Brain Atlas provide an unprecedented ability to examine these mechanisms in vivo; however, a lack of standardization across research groups has given rise to myriad processing pipelines for using these data. Here, we develop the abagen toolbox, an open-access software package for working with transcriptomic data, and use it to examine how methodological variability influences the outcomes of research using the Allen Human Brain Atlas. Applying three prototypical analyses to the outputs of 750,000 unique processing pipelines, we find that choice of pipeline has a large impact on research findings, with parameters commonly varied in the literature influencing correlations between derived gene expression and other imaging phenotypes by as much as ρ ≥ 1.0. Our results further reveal an ordering of parameter importance, with processing steps that influence gene normalization yielding the greatest impact on downstream statistical inferences and conclusions. The presented work and the development of the abagen toolbox lay the foundation for more standardized and systematic research in imaging transcriptomics, and will help to advance future understanding of the influence of gene expression in the human brain.

Dissecting the relationship between genotype and phenotype is one of the central goals in developmental biology and medicine. Transcriptome analysis is a powerful strategy to connect genotype to phenotype of a cell. Here we review the history, progress, potential applications and future developments of single-cell transcriptome analysis. In combination with live cell imaging and lineage tracing, it will be possible to decipher the full gene expression network underlying physiological functions of individual cells in embryos and adults, and to study diseases.

Stratification of asthma at the molecular level, especially using accessible biospecimens, could greatly enable patient selection for targeted therapy. To determine the value of blood analysis to identify transcriptional differences between clinically defined asthma and nonasthma groups, identify potential patient subgroups based on gene expression, and explore biological pathways associated with identified differences. Transcriptomic profiles were generated by microarray analysis of blood from 610 patients with asthma and control participants in the U-BIOPRED (Unbiased Biomarkers in Prediction of Respiratory Disease Outcomes) study. Differentially expressed genes (DEGs) were identified by analysis of variance, including covariates for RNA quality, sex, and clinical site, and Ingenuity Pathway Analysis was applied. Patient subgroups based on DEGs were created by hierarchical clustering and topological data analysis. A total of 1,693 genes were differentially expressed between patients with severe asthma and participants without asthma. The differences from participants without asthma in the nonsmoking severe asthma and mild/moderate asthma subgroups were significantly related (r = 0.76), with a larger effect size in the severe asthma group. The majority of, but not all, differences were explained by differences in circulating immune cell populations. Pathway analysis showed an increase in chemotaxis, migration, and myeloid cell trafficking in patients with severe asthma, decreased B-lymphocyte development and hematopoietic progenitor cells, and lymphoid organ hypoplasia. Cluster analysis of DEGs led to the creation of subgroups among the patients with severe asthma who differed in molecular responses to oral corticosteroids. Blood gene expression differences between clinically defined subgroups of patients with asthma and individuals without asthma, as well as subgroups of patients with severe asthma defined by transcript profiles, show the value of blood analysis in stratifying patients with asthma and identifying molecular pathways for further study. Clinical trial registered with www.clinicaltrials.gov (NCT01982162).

Many computational methods have been developed to infer cell type proportions from bulk transcriptomics data. However, an evaluation of the impact of data transformation, pre-processing, marker selection, cell type composition and choice of methodology on the deconvolution results is still lacking. Using five single-cell RNA-sequencing (scRNA-seq) datasets, we generate pseudo-bulk mixtures to evaluate the combined impact of these factors. Both bulk deconvolution methodologies and those that use scRNA-seq data as reference perform best when applied to data in linear scale and the choice of normalization has a dramatic impact on some, but not all methods. Overall, methods that use scRNA-seq data have comparable performance to the best performing bulk methods whereas semi-supervised approaches show higher error values. Moreover, failure to include cell types in the reference that are present in a mixture leads to substantially worse results, regardless of the previous choices. Altogether, we evaluate the combined impact of factors affecting the deconvolution task across different datasets and propose general guidelines to maximize its performance. Inferring cell type proportions from transcriptomics data is affected by data transformation, normalization, choice of method and the markers used. Here, the authors use single-cell RNAseq datasets to evaluate the impact of these factors and propose guidelines to maximise deconvolution performance.

Quantitative PCR (qPCR) has become the gold standard technique to measure cDNA and gDNA levels but the resulting data can be highly variable, artifactual and non-reproducible without appropriate verification and validation of both samples and primers. The root cause of poor quality data is typically associated with inadequate dilution of residual protein and chemical contaminants that variably inhibit Taq polymerase and primer annealing. The most susceptible, frustrating and often most interesting samples are those containing low abundant targets with small expression differences of 2-fold or lower. Here, Droplet Digital PCR (ddPCR) and qPCR platforms were directly compared for gene expression analysis using low amounts of purified, synthetic DNA in well characterized samples under identical reaction conditions. We conclude that for sample/target combinations with low levels of nucleic acids (Cq ≥ 29) and/or variable amounts of chemical and protein contaminants, ddPCR technology will produce more precise, reproducible and statistically significant results required for publication quality data. A stepwise methodology is also described to choose between these complimentary technologies to obtain the best results for any experiment.

Normalization of single-cell RNA sequencing data is necessary to eliminate cell-specific biases prior to downstream analyses. However, this is not straightforward for noisy single-cell data where many counts are zero. We present a novel approach where expression values are summed across pools of cells, and the summed values are used for normalization. Pool-based size factors are then deconvolved to yield cell-based factors. Our deconvolution approach outperforms existing methods for accurate normalization of cell-specific biases in simulated data. Similar behavior is observed in real data, where deconvolution improves the relevance of results of downstream analyses.

Purpose Multi-gene signatures provide biological insight and risk stratification in breast cancer. Intrinsic molecular subtypes defined by mRNA expression of 50 genes (PAM50) are prognostic in hormone-receptor positive postmenopausal breast cancer. Yet, for 25–40% in the PAM50 intermediate risk group, long-term risk remains uncertain. Our study aimed to (i) test the long-term prognostic value of the PAM50 signature in pre- and post-menopausal breast cancer; (ii) investigate if the PAM50 model could be improved by addition of other mRNAs implicated in oncogenesis. Methods We used archived FFPE samples from 1723 breast cancer survivors; high quality reads were obtained on 1253 samples. Transcript expression was quantified using a custom codeset with probes for > 100 targets. Cox models assessed gene signatures for breast cancer relapse and survival. Results Over 15 + years of follow-up, PAM50 subtypes were ( P  < 0.01) associated with breast cancer outcomes after accounting for tumor stage, grade and age at diagnosis. Results did not differ by menopausal status at diagnosis. Women with Luminal B (versus Luminal A) subtype had a > 60% higher hazard. Addition of a 13-gene hypoxia signature improved prognostication with > 40% higher hazard in the highest vs lowest hypoxia tertiles. Conclusions PAM50 intrinsic subtypes were independently prognostic for long-term breast cancer survival, irrespective of menopausal status. Addition of hypoxia signatures improved risk prediction. If replicated, incorporating the 13-gene hypoxia signature into the existing PAM50 risk assessment tool, may refine risk stratification and further clarify treatment for breast cancer.