The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

Despite the critical role of the cardiovascular system, our understanding of its cellular and transcriptional diversity remains limited. We therefore sought to characterize the cellular composition, phenotypes, molecular pathways, and communication networks between cell types at the tissue and sub-tissue level across the cardiovascular system of the healthy Wistar rat, an important model in preclinical cardiovascular research. We obtained high quality tissue samples under controlled conditions that reveal a level of cellular detail so far inaccessible in human studies.

CTLA-4 immune checkpoint blockade is clinically effective in a subset of patients with metastatic melanoma. We identify a subcluster of MAGE-A cancer-germline antigens, located within a narrow 75 kb region of chromosome Xq28, that predicts resistance uniquely to blockade of CTLA-4, but not PD-1. We validate this gene expression signature in an independent anti-CTLA-4-treated cohort and show its specificity to the CTLA-4 pathway with two independent anti-PD-1-treated cohorts. Autophagy, a process critical for optimal anti-cancer immunity, has previously been shown to be suppressed by the MAGE-TRIM28 ubiquitin ligase in vitro. We now show that the expression of the key autophagosome component LC3B and other activators of autophagy are negatively associated with MAGE-A protein levels in human melanomas, including samples from patients with resistance to CTLA-4 blockade. Our findings implicate autophagy suppression in resistance to CTLA-4 blockade in melanoma, suggesting exploitation of autophagy induction for potential therapeutic synergy with CTLA-4 inhibitors.

Mycobacterium tuberculosis (Mtb) can adopt a non-growing dormant state during infection that may be critical to both active and latent tuberculosis. During dormancy, Mtb is widely tolerant toward antibiotics, a significant obstacle in current anti-tubercular drug regimens, and retains the ability to persist in its environment. We aimed to identify novel mechanisms that permit Mtb to survive dormancy in an in vitro carbon starvation model using transposon insertion sequencing and gene expression analysis. We identified a previously uncharacterized component of the lipid transport machinery, omamC, which was upregulated and required for survival during carbon starvation. We show that OmamC plays a role both in increasing fatty acid stores during growth in rich media and enhancing fatty acid utilization during starvation. Besides its involvement in lipid metabolism, OmamC levels affected the expression of the anti-anti-sigma factor rv0516c and other genes to improve Mtb survival during carbon starvation and increase its tolerance toward rifampicin, a first-line drug effective against non-growing Mtb. Importantly, we show that Mtb can be eradicated during carbon starvation, in an OmamC-dependent manner, by inhibiting lipid metabolism with the lipase inhibitor tetrahydrolipstatin. This work casts new light into the survival processes of non-replicating, drug-tolerant Mtb by identifying new proteins involved in lipid metabolism required for the survival of dormant bacteria and exposing a potential vulnerability that could be exploited for antibiotic discovery.IMPORTANCETuberculosis is a global threat, with ~10 million yearly active cases. Many more people, however, live with "latent" infection, where Mycobacterium tuberculosis survives in a non-replicative form. When latent bacteria activate and regrow, they elicit immune responses and result in significant host damage. Replicating and non-growing bacilli can co-exist; however, non-growing bacteria are considerably less sensitive to antibiotics, thus complicating treatment by necessitating long treatment durations. Here, we sought to identify genes important for bacterial survival in this non-growing state using a carbon starvation model. We found that a previously uncharacterized gene, omamC, is involved in storing and utilizing fatty acids as bacteria transition between these two states. Importantly, inhibiting lipid metabolism using a lipase inhibitor eradicates non-growing bacteria. Thus, targeting lipid metabolism may be a viable strategy for treating the non-growing population in strategies to shorten treatment durations of tuberculosis.

Virulence factors generally enhance a pathogen's fitness and thereby foster transmission. However, most studies of pathogen fitness have been performed by averaging the phenotypes over large populations. Here, we have analyzed the fitness costs of virulence factor expression by Salmonella enterica subspecies I serovar Typhimurium in simple culture experiments. The type III secretion system ttss-1, a cardinal virulence factor for eliciting Salmonella diarrhea, is expressed by just a fraction of the S. Typhimurium population, yielding a mixture of cells that either express ttss-1 (TTSS-1(+) phenotype) or not (TTSS-1(-) phenotype). Here, we studied in vitro the TTSS-1(+) phenotype at the single cell level using fluorescent protein reporters. The regulator hilA controlled the fraction of TTSS-1+ individuals and their ttss-1 expression level. Strikingly, cells of the TTSS-1(+) phenotype grew slower than cells of the TTSS-1(-) phenotype. The growth retardation was at least partially attributable to the expression of TTSS-1 effector and/or translocon proteins. In spite of this growth penalty, the TTSS-1(+) subpopulation increased from <10% to approx. 60% during the late logarithmic growth phase of an LB batch culture. This was attributable to an increasing initiation rate of ttss-1 expression, in response to environmental cues accumulating during this growth phase, as shown by experimental data and mathematical modeling. Finally, hilA and hilD mutants, which form only fast-growing TTSS-1(-) cells, outcompeted wild type S. Typhimurium in mixed cultures. Our data demonstrated that virulence factor expression imposes a growth penalty in a non-host environment. This raises important questions about compensating mechanisms during host infection which ensure successful propagation of the genotype.

The intestinal microbiota features intricate metabolic interactions involving the breakdown and reuse of host- and diet-derived nutrients. The competition for these resources can limit pathogen growth. Nevertheless, some enteropathogenic bacteria can invade this niche through mechanisms that remain largely unclear. Using a mouse model for Salmonella diarrhea and a transposon mutant screen, we discovered that initial growth of Salmonella Typhimurium (S. Tm) in the unperturbed gut is powered by S. Tm hyb hydrogenase, which facilitates consumption of hydrogen (H2), a central intermediate of microbiota metabolism. In competitive infection experiments, a hyb mutant exhibited reduced growth early in infection compared to wild-type S. Tm, but these differences were lost upon antibiotic-mediated disruption of the host microbiota. Additionally, introducing H2-consuming bacteria into the microbiota interfered with hyb-dependent S. Tm growth. Thus, H2 is an Achilles' heel of microbiota metabolism that can be subverted by pathogens and might offer opportunities to prevent infection.

Repeated emergence of SARS-CoV-2 variants with increased fitness necessitates rapid detection and characterization of new lineages. To address this need, we developed PyR 0 , a hierarchical Bayesian multinomial logistic regression model that infers relative prevalence of all viral lineages across geographic regions, detects lineages increasing in prevalence, and identifies mutations relevant to fitness. Applying PyR 0 to all publicly available SARS-CoV-2 genomes, we identify numerous substitutions that increase fitness, including previously identified spike mutations and many non-spike mutations within the nucleocapsid and nonstructural proteins. PyR 0 forecasts growth of new lineages from their mutational profile, identifies viral lineages of concern as they emerge, and prioritizes mutations of biological and public health concern for functional characterization.

Copy number variants (CNVs) are major contributors to genetic diversity and disease. While standardized methods, such as the genome analysis toolkit (GATK), exist for detecting short variants, technical challenges have confounded uniform large-scale CNV analyses from whole-exome sequencing (WES) data. Given the profound impact of rare and de novo coding CNVs on genome organization and human disease, we developed GATK-gCNV, a flexible algorithm to discover rare CNVs from sequencing read-depth information, complete with open-source distribution via GATK. We benchmarked GATK-gCNV in 7,962 exomes from individuals in quartet families with matched genome sequencing and microarray data, finding up to 95% recall of rare coding CNVs at a resolution of more than two exons. We used GATK-gCNV to generate a reference catalog of rare coding CNVs in WES data from 197,306 individuals in the UK Biobank, and observed strong correlations between per-gene CNV rates and measures of mutational constraint, as well as rare CNV associations with multiple traits. In summary, GATK-gCNV is a tunable approach for sensitive and specific CNV discovery in WES data, with broad applications.

Type III secretion systems (TTSS) are employed by numerous pathogenic and symbiotic bacteria to inject a cocktail of different "effector proteins" into host cells. These effectors subvert host cell signaling to establish symbiosis or disease.

The SARS-CoV-2 pandemic has caused over 1 million deaths globally, mostly due to acute lung injury and acute respiratory distress syndrome, or direct complications resulting in multiple-organ failures. Little is known about the host tissue immune and cellular responses associated with COVID-19 infection, symptoms, and lethality. To address this, we collected tissues from 11 organs during the clinical autopsy of 17 individuals who succumbed to COVID-19, resulting in a tissue bank of approximately 420 specimens. We generated comprehensive cellular maps capturing COVID-19 biology related to patients' demise through single-cell and single-nucleus RNA-Seq of lung, kidney, liver and heart tissues, and further contextualized our findings through spatial RNA profiling of distinct lung regions. We developed a computational framework that incorporates removal of ambient RNA and automated cell type annotation to facilitate comparison with other healthy and diseased tissue atlases. In the lung, we uncovered significantly altered transcriptional programs within the epithelial, immune, and stromal compartments and cell intrinsic changes in multiple cell types relative to lung tissue from healthy controls. We observed evidence of: alveolar type 2 (AT2) differentiation replacing depleted alveolar type 1 (AT1) lung epithelial cells, as previously seen in fibrosis; a concomitant increase in myofibroblasts reflective of defective tissue repair; and, putative TP63 + intrapulmonary basal-like progenitor (IPBLP) cells, similar to cells identified in H1N1 influenza, that may serve as an emergency cellular reserve for severely damaged alveoli. Together, these findings suggest the activation and failure of multiple avenues for regeneration of the epithelium in these terminal lungs. SARS-CoV-2 RNA reads were enriched in lung mononuclear phagocytic cells and endothelial cells, and these cells expressed distinct host response transcriptional programs. We corroborated the compositional and transcriptional changes in lung tissue through spatial analysis of RNA profiles in situ and distinguished unique tissue host responses between regions with and without viral RNA, and in COVID-19 donor tissues relative to healthy lung. Finally, we analyzed genetic regions implicated in COVID-19 GWAS with transcriptomic data to implicate specific cell types and genes associated with disease severity. Overall, our COVID-19 cell atlas is a foundational dataset to better understand the biological impact of SARS-CoV-2 infection across the human body and empowers the identification of new therapeutic interventions and prevention strategies.

Effective treatment of bloodstream infections (BSIs) is relying on rapid identification of the causing pathogen and its antibiotic susceptibility. Still, most commercially available antibiotic susceptibility testing (AST) methods are based on monitoring bacterial growth, thus impacting the time to results. The Resistell AST is based on a new technology measuring the nanomotion caused by physiologically active bacterial cells and detecting the changes in nanomotion caused by the exposure to a drug.

SARS-CoV-2 distribution and circulation dynamics are not well understood due to challenges in assessing genomic data from tissue samples. We develop experimental and computational workflows for high-depth viral sequencing and high-resolution genomic analyses from formalin-fixed, paraffin-embedded tissues and apply them to 120 specimens from six subjects with fatal COVID-19. To varying degrees, viral RNA is present in extrapulmonary tissues from all subjects. The majority of the 180 viral variants identified within subjects are unique to individual tissue samples. We find more high-frequency (>10%) minor variants in subjects with a longer disease course, with one subject harboring ten such variants, exclusively in extrapulmonary tissues. One tissue-specific high-frequency variant was a nonsynonymous mutation in the furin-cleavage site of the spike protein. Our findings suggest adaptation and/or compartmentalized infection, illuminating the basis of extrapulmonary COVID-19 symptoms and potential for viral reservoirs, and have broad utility for investigating human pathogens.