Therapy-Induced Evolution of Human Lung Cancer Revealed by Single-Cell RNA Sequencing.

Lung cancer, the leading cause of cancer mortality, exhibits heterogeneity that enables adaptability, limits therapeutic success, and remains incompletely understood. Single-cell RNA sequencing (scRNA-seq) of metastatic lung cancer was performed using 49 clinical biopsies obtained from 30 patients before and during targeted therapy. Over 20,000 cancer and tumor microenvironment (TME) single-cell profiles exposed a rich and dynamic tumor ecosystem. scRNA-seq of cancer cells illuminated targetable oncogenes beyond those detected clinically. Cancer cells surviving therapy as residual disease (RD) expressed an alveolar-regenerative cell signature suggesting a therapy-induced primitive cell-state transition, whereas those present at on-therapy progressive disease (PD) upregulated kynurenine, plasminogen, and gap-junction pathways. Active T-lymphocytes and decreased macrophages were present at RD and immunosuppressive cell states characterized PD. Biological features revealed by scRNA-seq were biomarkers of clinical outcomes in independent cohorts. This study highlights how therapy-induced adaptation of the multi-cellular ecosystem of metastatic cancer shapes clinical outcomes.

Early prediction of preeclampsia in pregnancy with cell-free RNA.

Liquid biopsies that measure circulating cell-free RNA (cfRNA) offer an opportunity to study the development of pregnancy-related complications in a non-invasive manner and to bridge gaps in clinical care1-4. Here we used 404 blood samples from 199 pregnant mothers to identify and validate cfRNA transcriptomic changes that are associated with preeclampsia, a multi-organ syndrome that is the second largest cause of maternal death globally5. We find that changes in cfRNA gene expression between normotensive and preeclamptic mothers are marked and stable early in gestation, well before the onset of symptoms. These changes are enriched for genes specific to neuromuscular, endothelial and immune cell types and tissues that reflect key aspects of preeclampsia physiology6-9, suggest new hypotheses for disease progression and correlate with maternal organ health. This enabled the identification and independent validation of a panel of 18 genes that when measured between 5 and 16 weeks of gestation can form the basis of a liquid biopsy test that would identify mothers at risk of preeclampsia long before clinical symptoms manifest themselves. Tests based on these observations could help predict and manage who is at risk for preeclampsia-an important objective for obstetric care10,11.

A molecular cell atlas of the human lung from single-cell RNA sequencing.

Although single-cell RNA sequencing studies have begun to provide compendia of cell expression profiles1-9, it has been difficult to systematically identify and localize all molecular cell types in individual organs to create a full molecular cell atlas. Here, using droplet- and plate-based single-cell RNA sequencing of approximately 75,000 human cells across all lung tissue compartments and circulating blood, combined with a multi-pronged cell annotation approach, we create an extensive cell atlas of the human lung. We define the gene expression profiles and anatomical locations of 58 cell populations in the human lung, including 41 out of 45 previously known cell types and 14 previously unknown ones. This comprehensive molecular atlas identifies the biochemical functions of lung cells and the transcription factors and markers for making and monitoring them; defines the cell targets of circulating hormones and predicts local signalling interactions and immune cell homing; and identifies cell types that are directly affected by lung disease genes and respiratory viruses. By comparing human and mouse data, we identified 17 molecular cell types that have been gained or lost during lung evolution and others with substantially altered expression profiles, revealing extensive plasticity of cell types and cell-type-specific gene expression during organ evolution including expression switches between cell types. This atlas provides the molecular foundation for investigating how lung cell identities, functions and interactions are achieved in development and tissue engineering and altered in disease and evolution.

Ageing hallmarks exhibit organ-specific temporal signatures.

Ageing is the single greatest cause of disease and death worldwide, and understanding the associated processes could vastly improve quality of life. Although major categories of ageing damage have been identified-such as altered intercellular communication, loss of proteostasis and eroded mitochondrial function1-these deleterious processes interact with extraordinary complexity within and between organs, and a comprehensive, whole-organism analysis of ageing dynamics has been lacking. Here we performed bulk RNA sequencing of 17 organs and plasma proteomics at 10 ages across the lifespan of Mus musculus, and integrated these findings with data from the accompanying Tabula Muris Senis2-or 'Mouse Ageing Cell Atlas'-which follows on from the original Tabula Muris3. We reveal linear and nonlinear shifts in gene expression during ageing, with the associated genes clustered in consistent trajectory groups with coherent biological functions-including extracellular matrix regulation, unfolded protein binding, mitochondrial function, and inflammatory and immune response. Notably, these gene sets show similar expression across tissues, differing only in the amplitude and the age of onset of expression. Widespread activation of immune cells is especially pronounced, and is first detectable in white adipose depots during middle age. Single-cell RNA sequencing confirms the accumulation of T cells and B cells in adipose tissue-including plasma cells that express immunoglobulin J-which also accrue concurrently across diverse organs. Finally, we show how gene expression shifts in distinct tissues are highly correlated with corresponding protein levels in plasma, thus potentially contributing to the ageing of the systemic circulation. Together, these data demonstrate a similar yet asynchronous inter- and intra-organ progression of ageing, providing a foundation from which to track systemic sources of declining health at old age.

Optimization of whole-genome sequencing of Plasmodium falciparum from low-density dried blood spot samples.

BACKGROUND: Whole-genome sequencing (WGS) is becoming increasingly useful to study the biology, epidemiology, and ecology of malaria parasites. Despite ease of sampling, DNA extracted from dried blood spots (DBS) has a high ratio of human DNA compared to parasite DNA, which poses a challenge for downstream genetic analyses. The effects of multiple methods for DNA extraction, digestion of methylated DNA, and amplification were evaluated on the quality and fidelity of WGS data recovered from DBS. METHODS: Low parasite density mock DBS samples were created, extracted either with Tween-Chelex or QIAamp, treated with or without McrBC, and amplified with one of three different amplification techniques (two sWGA primer sets and one rWGA). Extraction conditions were evaluated on performance of sequencing depth, percentiles of coverage, and expected SNP concordance. RESULTS: At 100 parasites/µL, Chelex-Tween-McrBC samples had higher coverage (5 × depth = 93% genome) than QIAamp extracted samples (5 × depth = 76% genome). The two evaluated sWGA primer sets showed minor differences in overall genome coverage and SNP concordance, with a newly proposed combination of 20 primers showing a modest improvement in coverage over those previously published. CONCLUSIONS: Overall, Tween-Chelex extracted samples that were treated with McrBC digestion and are amplified using 6A10AD sWGA conditions had minimal dropout rate, higher percentages of coverage at higher depth, and more accurate SNP concordance than QiaAMP extracted samples. These findings extend the results of previously reported methods, making whole genome sequencing accessible to a larger number of low density samples that are commonly encountered in cross-sectional surveys.

FLASH: a next-generation CRISPR diagnostic for multiplexed detection of antimicrobial resistance sequences.

The growing prevalence of deadly microbes with resistance to previously life-saving drug therapies is a dire threat to human health. Detection of low abundance pathogen sequences remains a challenge for metagenomic Next Generation Sequencing (NGS). We introduce FLASH (Finding Low Abundance Sequences by Hybridization), a next-generation CRISPR/Cas9 diagnostic method that takes advantage of the efficiency, specificity and flexibility of Cas9 to enrich for a programmed set of sequences. FLASH-NGS achieves up to 5 orders of magnitude of enrichment and sub-attomolar gene detection with minimal background. We provide an open-source software tool (FLASHit) for guide RNA design. Here we applied it to detection of antimicrobial resistance genes in respiratory fluid and dried blood spots, but FLASH-NGS is applicable to all areas that rely on multiplex PCR.

Investigating Transfusion-related Sepsis Using Culture-Independent Metagenomic Sequencing.

BACKGROUND: Transfusion-related sepsis remains an important hospital infection control challenge. Investigation of septic transfusion events is often restricted by the limitations of bacterial culture in terms of time requirements and low yield in the setting of prior antibiotic administration. METHODS: In 3 gram-negative septic transfusion cases, we performed metagenomic next-generation sequencing (mNGS) of direct clinical blood specimens in addition to standard culture-based approaches utilized for infection control investigations. Pathogen detection leveraged IDSeq, a new open-access microbial bioinformatics portal. Phylogenetic analysis was performed to assess microbial genetic relatedness and understand transmission events. RESULTS: mNGS of direct clinical blood specimens afforded precision detection of pathogens responsible for each case of transfusion-related sepsis and enabled discovery of a novel Acinetobacter species in a platelet product that had become contaminated despite photochemical pathogen reduction. In each case, longitudinal assessment of pathogen burden elucidated the temporal sequence of events associated with each transfusion-transmitted infection. We found that informative data could be obtained from culture-independent mNGS of residual platelet products and leftover blood specimens that were either unsuitable or unavailable for culture or that failed to grow due to prior antibiotic administration. We additionally developed methods to enhance accuracy for detecting transfusion-associated pathogens that share taxonomic similarity to contaminants commonly found in mNGS library preparations. CONCLUSIONS: Culture-independent mNGS of blood products afforded rapid and precise assessment of pathogen identity, abundance, and genetic relatedness. Together, these challenging cases demonstrated the potential for metagenomics to advance existing methods for investigating transfusion-transmitted infections.

Genomic Profiling of Evolving Daptomycin Resistance in a Patient with Recurrent Staphylococcus argenteus Sepsis.

Staphylococcus argenteus is a novel staphylococcal species associated with invasive disease. We report the first case of daptomycin/vancomycin-resistant S. argenteus, initially speciated as Staphylococcus aureus, that developed from repeated treatment with daptomycin for a complex vascular graft infection. Whole-genome sequencing of longitudinally collected isolates identified acquisition of MprF S337L, a mutation predicted to increase surface charge and repel cationic molecules.

B cell repertoires in HLA-sensitized kidney transplant candidates undergoing desensitization therapy.

BACKGROUND: Kidney transplantation is the most effective treatment for end-stage renal disease. Sensitization refers to pre-existing antibodies against human leukocyte antigen (HLA) protein and remains a major barrier to successful transplantation. Despite implementation of desensitization strategies, many candidates fail to respond. Our objective was to determine whether measuring B cell repertoires could differentiate candidates that respond to desensitization therapy. METHODS: We developed an assay based on high-throughput DNA sequencing of the variable domain of the heavy chain of immunoglobulin genes to measure changes in B cell repertoires in 19 highly HLA-sensitized kidney transplant candidates undergoing desensitization and 7 controls with low to moderate HLA sensitization levels. Responders to desensitization had a decrease of 5% points or greater in cumulated calculated panel reactive antibody (cPRA) levels, and non-responders had no decrease in cPRA. RESULTS: Dominant B cell clones were not observed in highly sensitized candidates, suggesting that the B cells responsible for sensitization are either not present in peripheral blood or present at comparable levels to other circulating B cells. Candidates that responded to desensitization therapy had pre-treatment repertoires composed of a larger fraction of class-switched (IgG and IgA) isotypes compared to non-responding candidates. After B cell depleting therapy, the proportion of switched isotypes increased and the mutation frequencies of the remaining non-switched isotypes (IgM and IgD) increased in both responders and non-responders, perhaps representing a shift in the repertoire towards memory B cells or plasmablasts. Conversely, after transplantation, non-switched isotypes with fewer mutations increased, suggesting a shift in the repertoire towards naïve B cells. CONCLUSIONS: Relative abundance of different B cell isotypes is strongly perturbed by desensitization therapy and transplantation, potentially reflecting changes in the relative abundance of memory and naïve B cell compartments. Candidates that responded to therapy experienced similar changes to those that did not respond. Further studies are required to understand differences between these two groups of highly sensitized kidney transplant candidates.

Proteome-wide protein interaction measurements of bacterial proteins of unknown function.

Despite the enormous proliferation of bacterial genome data, surprisingly persistent collections of bacterial proteins have resisted functional annotation. In a typical genome, roughly 30% of genes have no assigned function. Many of these proteins are conserved across a large number of bacterial genomes. To assign a putative function to these conserved proteins of unknown function, we created a physical interaction map by measuring biophysical interaction of these proteins. Binary protein--protein interactions in the model organism Streptococcus pneumoniae (TIGR4) are measured with a microfluidic high-throughput assay technology. In some cases, informatic analysis was used to restrict the space of potential binding partners. In other cases, we performed in vitro proteome-wide interaction screens. We were able to assign putative functions to 50 conserved proteins of unknown function that we studied with this approach.

Genetic measurement of memory B-cell recall using antibody repertoire sequencing.

Annual influenza vaccinations aim to protect against seasonal infections, and vaccine strain compositions are updated every year. This protection is based on antibodies that are produced by either newly activated or memory B cells recalled from previous encounters with influenza vaccination or infection. The extent to which the B-cell repertoire responds to vaccination and recalls antibodies has so far not been analyzed at a genetic level--which is to say, at the level of antibody sequences. Here, we developed a consensus read sequencing approach that incorporates unique barcode labels on each starting RNA molecule. These labels allow one to combine multiple sequencing reads covering the same RNA molecule to reduce the error rate to a desired level, and they also enable accurate quantification of RNA and isotype levels. We validated this approach and analyzed the differential response of the antibody repertoire to live-attenuated or trivalent-inactivated influenza vaccination. Additionally, we analyzed the antibody repertoire in response to repeated yearly vaccinations with trivalent-inactivated influenza vaccination. We found antibody sequences that were present in both years, providing a direct genetic measurement of B-cell recall.

Systematic reconstruction of RNA functional motifs with high-throughput microfluidics.

We present RNA-mechanically induced trapping of molecular interactions (RNA-MITOMI), a microfluidic platform that allows integrated synthesis and functional assays for programmable RNA libraries. The interaction of a comprehensive library of RNA mutants with stem-loop-binding protein precisely defined the RNA structural and sequence features that govern affinity. The functional motif reconstructed in a single experiment on our platform uncovers new binding specificities and enriches interpretation of phylogenetic data.

The RootChip: an integrated microfluidic chip for plant science.

Studying development and physiology of growing roots is challenging due to limitations regarding cellular and subcellular analysis under controlled environmental conditions. We describe a microfluidic chip platform, called RootChip, that integrates live-cell imaging of growth and metabolism of Arabidopsis thaliana roots with rapid modulation of environmental conditions. The RootChip has separate chambers for individual regulation of the microenvironment of multiple roots from multiple seedlings in parallel. We demonstrate the utility of The RootChip by monitoring time-resolved growth and cytosolic sugar levels at subcellular resolution in plants by a genetically encoded fluorescence sensor for glucose and galactose. The RootChip can be modified for use with roots from other plant species by adapting the chamber geometry and facilitates the systematic analysis of root growth and metabolism from multiple seedlings, paving the way for large-scale phenotyping of root metabolism and signaling.

Characterizing expression changes in noncoding RNAs during aging and heterochronic parabiosis across mouse tissues.

Molecular mechanisms of organismal and cell aging remain incompletely understood. We, therefore, generated a body-wide map of noncoding RNA (ncRNA) expression in aging (16 organs at ten timepoints from 1 to 27 months) and rejuvenated mice. We found molecular aging trajectories are largely tissue-specific except for eight broadly deregulated microRNAs (miRNAs). Their individual abundance mirrors their presence in circulating plasma and extracellular vesicles (EVs) whereas tissue-specific ncRNAs were less present. For miR-29c-3p, we observe the largest correlation with aging in solid organs, plasma and EVs. In mice rejuvenated by heterochronic parabiosis, miR-29c-3p was the most prominent miRNA restored to similar levels found in young liver. miR-29c-3p targets the extracellular matrix and secretion pathways, known to be implicated in aging. We provide a map of organism-wide expression of ncRNAs with aging and rejuvenation and identify a set of broadly deregulated miRNAs, which may function as systemic regulators of aging via plasma and EVs.

High-performance binary protein interaction screening in a microfluidic format.

The standard procedure to increase microfluidic chip performance is to grow the number of parallel test systems on the chip. This process is accompanied by miniaturizing biochemical workflows and micromechanical elements, which is often a major challenge for both engineering fields. In this work, we show that it is possible to substantially increase the runtime performance of a microfluidic affinity assay for protein interactions by simultaneously engineering fluid logics and assay chemistry. For this, synergistic effects between the micro- and chemical architecture of the chip are exploited. The presented strategy of reducing the runtime rather than size and volume of the mechanical elements and biological reagent compartments will, in general, be of importance for future analytical test systems on microfluidic chips to overcome performance barriers.

Metagenomic characterization of swine slurry in a North American swine farm operation.

Modern day large-scale, high-density farming environments are inherently susceptible to viral outbreaks, inadvertently creating conditions that favor increased pathogen transmission and potential zoonotic spread. Metagenomic sequencing has proven to be a useful tool for characterizing the microbial burden in both people, livestock, and environmental samples. International efforts have been successful at characterizing pathogens in commercial farming environments, especially swine farms, however it is unclear whether the full extent of microbial agents have been adequately captured or is representative of farms elsewhere. To augment international efforts we performed metagenomic next-generation sequencing on nine swine slurry and three environmental samples from a United States of America (U.S.A.) farm operation, characterized the microbial composition of slurry, and identified novel viruses. We assembled a remarkable total of 1792 viral genomes, of which 554 were novel/divergent. We assembled 1637 Picobirnavirus genome segments, of which 538 are novel. In addition, we discovered 10 new viruses belonging to a novel taxon: porcine Statoviruses; which have only been previously reported in human, macaques, mouse, and cows. We assembled 3 divergent Posaviruses and 3 swine Picornaviruses. In addition to viruses described, we found other eukaryotic genera such as Entamoeba and Blastocystis, and bacterial genera such as Listeria, Treponema, Peptoclostridium and Bordetella in the slurry. Of these, two species Entamoeba histolytica and Listeria monocytogenes known to cause human disease were detected. Further, antimicrobial resistance genes such as tetracycline and MLS (macrolide, lincosamide, streptogramin) were also identified. Metagenomic surveillance in swine fecal slurry has great potential for novel and antimicrobial resistant pathogen detection.

Molecular characterization of selectively vulnerable neurons in Alzheimer's disease.

Alzheimer's disease (AD) is characterized by the selective vulnerability of specific neuronal populations, the molecular signatures of which are largely unknown. To identify and characterize selectively vulnerable neuronal populations, we used single-nucleus RNA sequencing to profile the caudal entorhinal cortex and the superior frontal gyrus-brain regions where neurofibrillary inclusions and neuronal loss occur early and late in AD, respectively-from postmortem brains spanning the progression of AD-type tau neurofibrillary pathology. We identified RORB as a marker of selectively vulnerable excitatory neurons in the entorhinal cortex and subsequently validated their depletion and selective susceptibility to neurofibrillary inclusions during disease progression using quantitative neuropathological methods. We also discovered an astrocyte subpopulation, likely representing reactive astrocytes, characterized by decreased expression of genes involved in homeostatic functions. Our characterization of selectively vulnerable neurons in AD paves the way for future mechanistic studies of selective vulnerability and potential therapeutic strategies for enhancing neuronal resilience.

Chloride channels regulate differentiation and barrier functions of the mammalian airway.

The conducting airway forms a protective mucosal barrier and is the primary target of airway disorders. The molecular events required for the formation and function of the airway mucosal barrier, as well as the mechanisms by which barrier dysfunction leads to early onset airway diseases, remain unclear. In this study, we systematically characterized the developmental landscape of the mouse airway using single-cell RNA sequencing and identified remarkably conserved cellular programs operating during human fetal development. We demonstrated that in mouse, genetic inactivation of chloride channel Ano1/Tmem16a compromises airway barrier function, results in early signs of inflammation, and alters the airway cellular landscape by depleting epithelial progenitors. Mouse Ano1-/-mutants exhibited mucus obstruction and abnormal mucociliary clearance that resemble the airway defects associated with cystic fibrosis. The data reveal critical and non-redundant roles for Ano1 in organogenesis, and show that chloride channels are essential for mammalian airway formation and function.

Ageing compromises mouse thymus function and remodels epithelial cell differentiation.

Ageing is characterised by cellular senescence, leading to imbalanced tissue maintenance, cell death and compromised organ function. This is first observed in the thymus, the primary lymphoid organ that generates and selects T cells. However, the molecular and cellular mechanisms underpinning these ageing processes remain unclear. Here, we show that mouse ageing leads to less efficient T cell selection, decreased self-antigen representation and increased T cell receptor repertoire diversity. Using a combination of single-cell RNA-seq and lineage-tracing, we find that progenitor cells are the principal targets of ageing, whereas the function of individual mature thymic epithelial cells is compromised only modestly. Specifically, an early-life precursor cell population, retained in the mouse cortex postnatally, is virtually extinguished at puberty. Concomitantly, a medullary precursor cell quiesces, thereby impairing maintenance of the medullary epithelium. Thus, ageing disrupts thymic progenitor differentiation and impairs the core immunological functions of the thymus.

Pan-viral serology implicates enteroviruses in acute flaccid myelitis.

Since 2012, the United States of America has experienced a biennial spike in pediatric acute flaccid myelitis (AFM)1-6. Epidemiologic evidence suggests non-polio enteroviruses (EVs) are a potential etiology, yet EV RNA is rarely detected in cerebrospinal fluid (CSF)2. CSF from children with AFM (n = 42) and other pediatric neurologic disease controls (n = 58) were investigated for intrathecal antiviral antibodies, using a phage display library expressing 481,966 overlapping peptides derived from all known vertebrate and arboviruses (VirScan). Metagenomic next-generation sequencing (mNGS) of AFM CSF RNA (n = 20 cases) was also performed, both unbiased sequencing and with targeted enrichment for EVs. Using VirScan, the viral family significantly enriched by the CSF of AFM cases relative to controls was Picornaviridae, with the most enriched Picornaviridae peptides belonging to the genus Enterovirus (n = 29/42 cases versus 4/58 controls). EV VP1 ELISA confirmed this finding (n = 22/26 cases versus 7/50 controls). mNGS did not detect additional EV RNA. Despite rare detection of EV RNA, pan-viral serology frequently identified high levels of CSF EV-specific antibodies in AFM compared with controls, providing further evidence for a causal role of non-polio EVs in AFM.