Integrated analyses of single-cell atlases reveal age, gender, and smoking status associations with cell type-specific expression of mediators of SARS-CoV-2 viral entry and highlights inflammatory programs in putative target cells

Christoph Muus; Malte D Luecken; Gokcen Eraslan; Avinash Waghray; Graham Heimberg; Lisa Sikkema; Yoshihiko Kobayashi; Eeshit Dhaval Vaishnav; Ayshwarya Subramanian; Christopher Smillie; Karthik Jagadeesh; Elizabeth Thu Duong; Evgenij Fiskin; Elena Torlai Triglia; Christophe Becavin; Meshal Ansari; Peiwen Cai; Brian Lin; Justin Buchanan; Sijia Chen; Jian Shu; Adam L Haber; Hattie Chung; Daniel T Montoro; Taylor Adams; Hananeh Aliee; Samuel J Allon; Zaneta Andrusivova; Ilias Angelidis; Orr Ashenberg; Kevin Bassler; Christophe Becavin; Inbal Benhar; Joseph Bergenstrahle; Ludvig Bergenstrahle; Liam Bolt; Emelie Braun; Linh T Bui; Mark Chaffin; Evgeny Chichelnitskiy; Joshua Chiou; Thomas M Conlon; Michael S Cuoco; Marie Deprez; David S Fischer; Astrid Gillich; Joshua Gould; Minzhe Guo; Austin J Gutierrez; Arun C Habermann; Tyler Harvey; Peng He; Xiaomeng Hou; Lijuan Hu; Alok Jaiswal; Peiyong Jiang; Theodoros Kapellos; Christin S Kuo; Ludvig Larsson; Michael A Leney-Greene; Kyungtae Lim; Monika Litvinukova; Ji Lu; Leif S Ludwig; Wendy Luo; Henrike Maatz; Elo Maddissoon; Lira Mamanova; Kasidet Manakongtreecheep; Charles-Hugo Marquette; Ian Mbano; Alexi M McAdams; Ross J Metzger; Ahmad N Nabhan; Sarah K Nyquist; Jose Ordovas-Montanes; Lolita Penland; Olivier B Poirion; Segio Poli; CanCan Qi; Daniel Reichart; Ivan Rosas; Jonas Schupp; Rahul Sinha; Rene V Sit; Kamil Slowikowski; Michal Slyper; Neal Smith; Alex Sountoulidis; Maximilian Strunz; Dawei Sun; Carlos Talavera-Lopez; Peng Tan; Jessica Tantivit; Kyle J Travaglini; Nathan R Tucker; Katherine Vernon; Marc H Wadsworth; Julia Waldman; Xiuting Wang; Wenjun Yan; Ali Onder Yildirim; William Zhao; Carly G K Ziegler; Aviv Regev

Here, "5 denotes the raw count expression of gene i in cell j and age:sex, sex:smoking, and age:smoking are interaction terms of the three modeled covariates. These terms model whether there is a difference in the smoking effect in men and women, and likewise whether the age effect is different for smokers and non-smokers. While we model these interaction terms, we only tested age, sex, and smoking effects individually to reduce the multiple testing burden. We included the dataset term to model the batch effects between the diverse datasets we obtained, and the log total counts per cell was used as an offset. Here, the total counts were scaled to have a mean of 1 across all cells before the log was taken. In order to fit this model we pruned the data to contain only datasets that have at least 2 donors and for which smoking status metadata was provided. This resulted in a dataset of 761,693 cells and 165 samples from 77 donors for adult lung data. Only a single dataset remained for adult nasal data after this filtering on which the model could not be fit. To obtain cell-type specific associations the above model was fit within each cell type for all cell types with at least 1,000 cells. We performed Wald tests over the age and sex covariates independently and corrected for multiple testing via Benjamini-Hochberg over all tests within a cell type.